1
|
Witek J, Heindel JP, Guan X, Leven I, Hao H, Naullage P, LaCour A, Sami S, Menger MFSJ, Cofer-Shabica DV, Berquist E, Faraji S, Epifanovsky E, Head-Gordon T. M-Chem: a Modular Software Package for Molecular Simulation that Spans Scientific Domains. Mol Phys 2022; 121:e2129500. [PMID: 37470065 PMCID: PMC10353727 DOI: 10.1080/00268976.2022.2129500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2022] [Accepted: 09/06/2022] [Indexed: 10/10/2022]
Abstract
We present a new software package called M-Chem that is designed from scratch in C++ and parallelized on shared-memory multi-core architectures to facilitate efficient molecular simulations. Currently, M-Chem is a fast molecular dynamics (MD) engine that supports the evaluation of energies and forces from two-body to many-body all-atom potentials, reactive force fields, coarse-grained models, combined quantum mechanics molecular mechanics (QM/MM) models, and external force drivers from machine learning, augmented by algorithms that are focused on gains in computational simulation times. M-Chem also includes a range of standard simulation capabilities including thermostats, barostats, multi-timestepping, and periodic cells, as well as newer methods such as fast extended Lagrangians and high quality electrostatic potential generation. At present M-Chem is a developer friendly environment in which we encourage new software contributors from diverse fields to build their algorithms, models, and methods in our modular framework. The long-term objective of M-Chem is to create an interdisciplinary platform for computational methods with applications ranging from biomolecular simulations, reactive chemistry, to materials research.
Collapse
Affiliation(s)
- Jagna Witek
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | - Joseph P Heindel
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
| | - Xingyi Guan
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
| | - Itai Leven
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | - Hongxia Hao
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | | | - Allen LaCour
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
| | - Selim Sami
- Kenneth S. Pitzer Theory Center and Department of Chemistry
| | - M F S J Menger
- Stratingh Institute for Chemistry, University of Groningen, 9747 AG Groningen, The Netherlands
| | - D Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA, 19128 USA
| | - Eric Berquist
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Shirin Faraji
- Stratingh Institute for Chemistry, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry
- Chemical Sciences Division, Lawrence Berkeley National Laboratory
- Department of Bioengineering and Chemical and Biomolecular Engineering University of California, Berkeley, CA, USA
| |
Collapse
|
2
|
Haghighatlari M, Li J, Guan X, Zhang O, Das A, Stein CJ, Heidar-Zadeh F, Liu M, Head-Gordon M, Bertels L, Hao H, Leven I, Head-Gordon T. NewtonNet: a Newtonian message passing network for deep learning of interatomic potentials and forces. Digit Discov 2022; 1:333-343. [PMID: 35769203 PMCID: PMC9189860 DOI: 10.1039/d2dd00008c] [Citation(s) in RCA: 36] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Accepted: 04/26/2022] [Indexed: 04/14/2023]
Abstract
We report a new deep learning message passing network that takes inspiration from Newton's equations of motion to learn interatomic potentials and forces. With the advantage of directional information from trainable force vectors, and physics-infused operators that are inspired by Newtonian physics, the entire model remains rotationally equivariant, and many-body interactions are inferred by more interpretable physical features. We test NewtonNet on the prediction of several reactive and non-reactive high quality ab initio data sets including single small molecules, a large set of chemically diverse molecules, and methane and hydrogen combustion reactions, achieving state-of-the-art test performance on energies and forces with far greater data and computational efficiency than other deep learning models.
Collapse
Affiliation(s)
- Mojtaba Haghighatlari
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
| | - Jie Li
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
| | - Xingyi Guan
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA USA
| | - Oufan Zhang
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
| | - Akshaya Das
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
| | - Christopher J Stein
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA USA
- Theoretical Physics and Center for Nanointegration Duisburg-Essen (CENIDE), University of Duisburg-Essen 47048 Duisburg Germany
| | - Farnaz Heidar-Zadeh
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Center for Molecular Modeling (CMM), Ghent University B-9052 Ghent Belgium
- Department of Chemistry, Queen's University Kingston Ontario K7L 3N6 Canada
| | - Meili Liu
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Department of Chemistry, Beijing Normal University Beijing 100875 China
| | - Martin Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA USA
| | - Luke Bertels
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
| | - Hongxia Hao
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA USA
| | - Itai Leven
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA USA
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California Berkeley CA USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA USA
- Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California Berkeley CA USA
| |
Collapse
|
3
|
Guan X, Das A, Stein CJ, Heidar-Zadeh F, Bertels L, Liu M, Haghighatlari M, Li J, Zhang O, Hao H, Leven I, Head-Gordon M, Head-Gordon T. A benchmark dataset for Hydrogen Combustion. Sci Data 2022; 9:215. [PMID: 35581204 PMCID: PMC9114378 DOI: 10.1038/s41597-022-01330-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Accepted: 04/20/2022] [Indexed: 11/21/2022] Open
Abstract
The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ∼290,000 potential energies and ∼1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction. Measurement(s) | ab initio energies and forces of hydrogen combustion | Technology Type(s) | density functional theory • ab initio molecular dynamics • normal modes | Factor Type(s) | cartesian coordinates |
Collapse
Affiliation(s)
- Xingyi Guan
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Akshaya Das
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA
| | - Christopher J Stein
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA.,Theoretical Physics and Center for Nanointegration Duisburg-Essen (CENIDE), University of Duisburg-Essen, 47048, Duisburg, Germany
| | - Farnaz Heidar-Zadeh
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Department of Chemistry, Queen's University, Kingston, Ontario, K7L 3N6, Canada
| | - Luke Bertels
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA
| | - Meili Liu
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Department of Chemistry, Beijing Normal University, Beijing, 100875, China
| | - Mojtaba Haghighatlari
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jie Li
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA
| | - Oufan Zhang
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA
| | - Hongxia Hao
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Itai Leven
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, CA, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA. .,Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, CA, USA.
| |
Collapse
|
4
|
Abstract
Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O-H bonds at the surface are ~16 MV/cm larger on average than that found for O-H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.
Collapse
Affiliation(s)
- Hongxia Hao
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA, 94720, USA
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, CA, 94720, USA
- Departments of Chemistry, University of California, Berkeley, CA, 94720, USA
| | - Itai Leven
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA, 94720, USA
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, CA, 94720, USA
- Departments of Chemistry, University of California, Berkeley, CA, 94720, USA
| | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA, 94720, USA.
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, CA, 94720, USA.
- Departments of Chemistry, University of California, Berkeley, CA, 94720, USA.
- Departments of Bioengineering, University of California, Berkeley, CA, 94720, USA.
- Departments of Chemical and Biomolecular Engineering, University of California, Berkeley, CA, 94720, USA.
| |
Collapse
|
5
|
Adams EM, Hao H, Leven I, Rüttermann M, Wirtz H, Havenith M, Head‐Gordon T. Rücktitelbild: Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism (Angew. Chem. 48/2021). Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202111736] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Ellen M. Adams
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Hongxia Hao
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
| | - Itai Leven
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
| | | | - Hanna Wirtz
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Martina Havenith
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Teresa Head‐Gordon
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley California 94720 USA
- Department of Bioengineering University of California Berkeley California 94720 USA
| |
Collapse
|
6
|
Adams EM, Hao H, Leven I, Rüttermann M, Wirtz H, Havenith M, Head‐Gordon T. Back Cover: Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism (Angew. Chem. Int. Ed. 48/2021). Angew Chem Int Ed Engl 2021. [DOI: 10.1002/anie.202111736] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Ellen M. Adams
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Hongxia Hao
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
| | - Itai Leven
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
| | | | - Hanna Wirtz
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Martina Havenith
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Teresa Head‐Gordon
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley California 94720 USA
- Department of Bioengineering University of California Berkeley California 94720 USA
| |
Collapse
|
7
|
Adams EM, Hao H, Leven I, Rüttermann M, Wirtz H, Havenith M, Head‐Gordon T. Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202108766] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Ellen M. Adams
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Hongxia Hao
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
| | - Itai Leven
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
| | | | - Hanna Wirtz
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Martina Havenith
- Lehrstuhl für Physkalische Chemie II Ruhr Universität Bochum 44801 Bochum Germany
| | - Teresa Head‐Gordon
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California 94720 USA
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley California 94720 USA
- Department of Chemistry University of California Berkeley California 94720 USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley California 94720 USA
- Department of Bioengineering University of California Berkeley California 94720 USA
| |
Collapse
|
8
|
Guan X, Leven I, Heidar-Zadeh F, Head-Gordon T. Protein C-GeM: A Coarse-Grained Electron Model for Fast and Accurate Protein Electrostatics Prediction. J Chem Inf Model 2021; 61:4357-4369. [PMID: 34490776 DOI: 10.1021/acs.jcim.1c00388] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The electrostatic potential (ESP) is a powerful property for understanding and predicting electrostatic charge distributions that drive interactions between molecules. In this study, we compare various charge partitioning schemes including fitted charges, density-based quantum mechanical (QM) partitioning schemes, charge equilibration methods, and our recently introduced coarse-grained electron model, C-GeM, to describe the ESP for protein systems. When benchmarked against high quality density functional theory calculations of the ESP for tripeptides and the crambin protein, we find that the C-GeM model is of comparable accuracy to ab initio charge partitioning methods, but with orders of magnitude improvement in computational efficiency since it does not require either the electron density or the electrostatic potential as input.
Collapse
Affiliation(s)
- Xingyi Guan
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Itai Leven
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Farnaz Heidar-Zadeh
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States.,Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Havenith-Newen M, Adams EM, Head-Gordon T, Hao H, Rüttermann M, Leven I, Wirtz H. Proton Traffic Jam: Effect of Nanoconfinement and Acid Concentration on Proton Hopping Mechanism. Angew Chem Int Ed Engl 2021; 60:25419-25427. [PMID: 34402145 PMCID: PMC9293324 DOI: 10.1002/anie.202108766] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Indexed: 11/06/2022]
Abstract
The properties of the water network in concentrated HCl acid pools in nanometer-sized reverse non-ionic micelles were probed with TeraHertz absorption, dielectric relaxation spectroscopy, and reactive force field simulations capable of describing proton hopping mechanisms. We identify that only at a critical micelle size of W0=9 do solvated proton complexes form in the water pool, accompanied by a change in mechanism from Grotthuss forward shuttling to one that favors local oscillatory hopping. This is due to a preference for H+ and Cl- ions to adsorb to the micelle interface, together with an acid concentration effect that causes a "traffic jam" in which the short-circuiting of the hydrogen-bonding motif of the hydronium ion decreases the forward hopping rate throughout the water interior even as the micelle size increases. These findings have implications for atmospheric chemistry, biochemical and biophysical environments, and energy materials, as transport of protons vital to these processes can be suppressed due to confinement, aggregation, and/or concentration.
Collapse
Affiliation(s)
- Martina Havenith-Newen
- Ruhr-Universit�t Bochum, Physical Chemistry, Universit�tsstr. 150, 44780, Bochum, GERMANY
| | - Ellen M Adams
- Ruhr-Universität Bochum: Ruhr-Universitat Bochum, Chemistry and Biochemistry, GERMANY
| | - Teresa Head-Gordon
- UC Berkeley: University of California Berkeley, Chemistry, UNITED STATES
| | - Hongxia Hao
- Berkeley Laboratory: E O Lawrence Berkeley National Laboratory, Chemistry, UNITED STATES
| | | | - Itai Leven
- Lawrence Livermore National Laboratory, chemistry, GERMANY
| | - Hanna Wirtz
- Ruhr-Universität Bochum: Ruhr-Universitat Bochum, Chemistry, GERMANY
| |
Collapse
|
10
|
Leven I, Hao H, Tan S, Guan X, Penrod KA, Akbarian D, Evangelisti B, Hossain MJ, Islam MM, Koski JP, Moore S, Aktulga HM, van Duin ACT, Head-Gordon T. Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields. J Chem Theory Comput 2021; 17:3237-3251. [PMID: 33970642 DOI: 10.1021/acs.jctc.1c00118] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. In this Perspective, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.
Collapse
Affiliation(s)
- Itai Leven
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National LaboratoryBerkeley, California 94720, United States
| | - Hongxia Hao
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National LaboratoryBerkeley, California 94720, United States
| | - Songchen Tan
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xingyi Guan
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National LaboratoryBerkeley, California 94720, United States
| | - Katheryn A Penrod
- Department of Mechanical Engineering, Chemical Engineering, Engineering Science and Mechanics, Chemistry, Materials Science and Engineering, Penn State University, 240 Research East, University Park, Pennsylvania 16802, United States
| | - Dooman Akbarian
- Department of Mechanical Engineering, Chemical Engineering, Engineering Science and Mechanics, Chemistry, Materials Science and Engineering, Penn State University, 240 Research East, University Park, Pennsylvania 16802, United States
| | - Benjamin Evangelisti
- Department of Mechanical Engineering, Chemical Engineering, Engineering Science and Mechanics, Chemistry, Materials Science and Engineering, Penn State University, 240 Research East, University Park, Pennsylvania 16802, United States
| | - Md Jamil Hossain
- Department of Mechanical Engineering, Chemical Engineering, Engineering Science and Mechanics, Chemistry, Materials Science and Engineering, Penn State University, 240 Research East, University Park, Pennsylvania 16802, United States
| | - Md Mahbubul Islam
- Department of Mechanical Engineering, Wayne State University, Detroit, Michigan 48202, United States
| | - Jason P Koski
- Sandia National Laboratories, Albuquerque, New Mexico 87185-1315, United States
| | - Stan Moore
- Sandia National Laboratories, Albuquerque, New Mexico 87185-1315, United States
| | - Hasan Metin Aktulga
- Department of Computer Science and Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| | - Adri C T van Duin
- Department of Mechanical Engineering, Chemical Engineering, Engineering Science and Mechanics, Chemistry, Materials Science and Engineering, Penn State University, 240 Research East, University Park, Pennsylvania 16802, United States
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National LaboratoryBerkeley, California 94720, United States.,Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
11
|
Abstract
Nonreactive force fields are defined by perturbations of electron density that are relatively small, whereas chemical reactivity involves wholesale electronic rearrangements that make and break bonds. Thus, reactive force fields are incredibly difficult to develop compared to nonreactive force fields, yet at the same time, they fill a critical need when ab initio molecular dynamics methods are not affordable. We introduce a new reactive force field model for water that combines modified nonbonded terms of the ReaxFF model and its embedding in the electrostatic interactions described by our recently introduced coarse-grained electron model (C-GeM). The ReaxFF/C-GeM force field is characterized for many energetic and dissociative water properties for water clusters, structure, and dynamical properties under ambient conditions in the condensed phase, as well as the temperature dependence of density and water diffusion, with very good agreement with experiment. The ReaxFF/C-GeM force field should be more transferable and more broadly applicable to a range of reactive systems involving both proton and electron transfer in the condensed phase.
Collapse
Affiliation(s)
- Itai Leven
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Hongxia Hao
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | | | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
12
|
Tan S, Leven I, An D, Lin L, Head-Gordon T. Stochastic Constrained Extended System Dynamics for Solving Charge Equilibration Models. J Chem Theory Comput 2020; 16:5991-5998. [PMID: 32956587 DOI: 10.1021/acs.jctc.0c00514] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a new stochastic extended Lagrangian molecular dynamics solution to charge equilibration that eliminates self-consistent field (SCF) calculations, thus eliminating the computational bottleneck in solving the charge distribution with standard SCF solvers. By formulating both charges and chemical potential as latent variables and introducing a holonomic constraint that satisfies charge conservation, the SC-XLMD method accurately reproduces thermodynamic, dynamic, and structural properties within the framework of ReaxFF for a bulk water system and highly reactive RDX molecules simulated at high temperature. The SC-XLMD method shows excellent computational performance and is available in the publicly available LAMMPS package.
Collapse
Affiliation(s)
- Songchen Tan
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Kenneth S. Pitzer Theory Center, University of California, Berkeley, California 94720, United States
| | - Itai Leven
- Kenneth S. Pitzer Theory Center, University of California, Berkeley, California 94720, United States.,Chemical Sciences, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Dong An
- Department of Mathematics, University of California, Berkeley, California 94720, United States.,Computational Research, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Lin Lin
- Department of Mathematics, University of California, Berkeley, California 94720, United States.,Computational Research, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center, University of California, Berkeley, California 94720, United States.,Chemical Sciences, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States.,Departments of Chemistry, Bioengineering, and Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
13
|
Abstract
We have developed a new coarse-grained electron model, C-GeM, in which atoms are represented by a positive core and an electron shell described by Gaussian charge distributions, with the interaction energy between the core and shell reflecting the electronegativity of a given atomic element. By minimizing the electronic shell positions in the field of atomic core positions, the model can provide accurate electrostatic properties of molecules and their interactions. We have tested the performance of the C-GeM model for a set of molecules containing H, C, O, and Cl atoms to show that it can predict the electrostatic potential with high accuracy, and correctly describe the dissociation of HCl into ionic fragments in solution and to neutral atoms in the gas phase. The resulting C-GeM approach offers many advantages over expensive ab initio methods and reactive force field charge equilibration methodologies: it can rapidly predict the electrostatic potential surfaces of molecules, molecules dissociate into integer charge fragments so that redox reactions are easily described, there is no unphysical long-range charge transfer, it can account for out-of-plane polarization, and charges are not required to be centered on atoms, thereby accounting for electrostatic features such as sigma holes.
Collapse
Affiliation(s)
- Itai Leven
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Teresa Head-Gordon
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
14
|
Das AK, Urban L, Leven I, Loipersberger M, Aldossary A, Head-Gordon M, Head-Gordon T. Development of an Advanced Force Field for Water Using Variational Energy Decomposition Analysis. J Chem Theory Comput 2019; 15:5001-5013. [DOI: 10.1021/acs.jctc.9b00478] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
15
|
Leven I, Levy Y. Quantifying the two-state facilitated diffusion model of protein-DNA interactions. Nucleic Acids Res 2019; 47:5530-5538. [PMID: 31045207 PMCID: PMC6582340 DOI: 10.1093/nar/gkz308] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2019] [Revised: 04/13/2019] [Accepted: 04/17/2019] [Indexed: 01/13/2023] Open
Abstract
The current report extends the facilitated diffusion model to account for conflict between the search and recognition binding modes adopted by DNA-binding proteins (DBPs) as they search DNA and subsequently recognize and bind to their specific binding site. The speed of the search dynamics is governed by the energetic ruggedness of the protein-DNA landscape, whereas the rate for the recognition process is mostly dictated by the free energy barrier for the transition between the DBP's search and recognition binding modes. We show that these two modes are negatively coupled, such that fast 1D sliding and rapid target site recognition probabilities are unlikely to coexist. Thus, a tradeoff occurs between optimizing the timescales for finding and binding the target site. We find that these two kinetic properties can be balanced to produce a fast timescale for the total target search and recognition process by optimizing frustration. Quantification of the facilitated diffusion model by including a frustration term enables it to explain several experimental observations concerning search and recognition speeds. The extended model captures experimental estimate of the energetic ruggedness of the protein-DNA landscape and predicts how various molecular properties of protein-DNA binding affect recognition kinetics. Particularly, point mutations may change the frustration and so affect protein association with DNA, thus providing a means to modulate protein-DNA affinity by manipulating the protein's association or dissociation reactions.
Collapse
Affiliation(s)
- Itai Leven
- Department of Structural Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Yaakov Levy
- Department of Structural Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| |
Collapse
|
16
|
Abstract
The inertial EL/SCF method is developed to solve charge equilibration models for molecular dynamics, reducing the number of SCFs by 50–80% at each time step when compared to a conjugate gradient SCF solver and tested on diverse reactive systems.
Collapse
Affiliation(s)
- Itai Leven
- Kenneth S. Pitzer Center for Theoretical Chemistry
- University of California Berkeley
- Berkeley
- USA
- Department of Chemistry
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry
- University of California Berkeley
- Berkeley
- USA
- Department of Chemistry
| |
Collapse
|
17
|
Guerra R, Leven I, Vanossi A, Hod O, Tosatti E. Smallest Archimedean Screw: Facet Dynamics and Friction in Multiwalled Nanotubes. Nano Lett 2017; 17:5321-5328. [PMID: 28795813 PMCID: PMC5600185 DOI: 10.1021/acs.nanolett.7b01718] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2017] [Revised: 08/10/2017] [Indexed: 05/28/2023]
Abstract
We identify a new material phenomenon, where minute mechanical manipulations induce pronounced global structural reconfigurations in faceted multiwalled nanotubes. This behavior has strong implications on the tribological properties of these systems and may be the key to understand the enhanced interwall friction recently measured for boron-nitride nanotubes with respect to their carbon counterparts. Notably, the fast rotation of helical facets in these systems upon coaxial sliding may serve as a nanoscale Archimedean screw for directional transport of physisorbed molecules.
Collapse
Affiliation(s)
- Roberto Guerra
- Center
for Complexity and Biosystems, Department of Physics, University of Milan, 20133 Milan, Italy
| | - Itai Leven
- Department
of Physical Chemistry, School of Chemistry, The Raymond and Beverly
Sackler Faculty of Exact Sciences and The Sackler Center for Computational
Molecular and Materials Science, Tel Aviv
University, Tel Aviv 6997801, Israel
| | - Andrea Vanossi
- International
School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- CNR-IOM
Democritos National Simulation Center, Via Bonomea 265, 34136 Trieste, Italy
| | - Oded Hod
- Department
of Physical Chemistry, School of Chemistry, The Raymond and Beverly
Sackler Faculty of Exact Sciences and The Sackler Center for Computational
Molecular and Materials Science, Tel Aviv
University, Tel Aviv 6997801, Israel
| | - Erio Tosatti
- International
School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- CNR-IOM
Democritos National Simulation Center, Via Bonomea 265, 34136 Trieste, Italy
- The
Abdus Salam International Centre for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
18
|
Mandelli D, Leven I, Hod O, Urbakh M. Sliding friction of graphene/hexagonal -boron nitride heterojunctions: a route to robust superlubricity. Sci Rep 2017; 7:10851. [PMID: 28883489 PMCID: PMC5589749 DOI: 10.1038/s41598-017-10522-8] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2017] [Accepted: 08/09/2017] [Indexed: 11/09/2022] Open
Abstract
The origin of ultra-low friction exhibited by heterogeneous junctions of graphene and hexagonal boron nitride (h-BN) is revealed. For aligned interfaces, we identify a characteristic contact size, below which the junction behaves like its homogeneous counterparts with friction forces that grow linearly with the contact area. Superlubricity sets in due to the progressive appearance of Moiré patterns resulting in a collective stick-slip motion of the elevated super-structure ridges that turns into smooth soliton-like gliding with increasing contact size. Incommensurability effects are enhanced in misaligned contacts, where the friction coefficients further drop by orders of magnitude. Our fully atomistic simulations show that the superlubric regime in graphene/h-BN heterostructures persists up to significantly higher loads compared to the well-studied twisted homogeneous graphene interface. This indicates the potential of achieving robust superlubricity in practical applications using two-dimensional layered materials heterojunctions.
Collapse
Affiliation(s)
- D Mandelli
- School of Chemistry and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801, Tel Aviv, Israel
| | - I Leven
- School of Chemistry and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801, Tel Aviv, Israel
| | - O Hod
- School of Chemistry and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801, Tel Aviv, Israel
| | - M Urbakh
- School of Chemistry and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, 6997801, Tel Aviv, Israel.
| |
Collapse
|
19
|
Leven I, Guerra R, Vanossi A, Tosatti E, Hod O. Multiwalled nanotube faceting unravelled. Nat Nanotechnol 2016; 11:1082-1086. [PMID: 27548359 DOI: 10.1038/nnano.2016.151] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Accepted: 07/15/2016] [Indexed: 05/28/2023]
Abstract
Nanotubes show great promise for miniaturizing advanced technologies. Their exceptional physical properties are intimately related to their morphological and crystal structure. Circumferential faceting of multiwalled nanotubes reinforces their mechanical strength and alters their tribological and electronic properties. Here, the nature of this important phenomenon is fully rationalized in terms of interlayer registry patterns. Regardless of the nanotube identity (that is, diameter, chirality, chemical composition), faceting requires the matching of the chiral angles of adjacent layers. Above a critical diameter that corresponds well with experimental results, achiral multiwalled nanotubes display evenly spaced extended axial facets whose number equals the interlayer difference in circumferential unit cells. Elongated helical facets, commonly observed in experiment, appear in nanotubes that exhibit small interlayer chiral angle mismatch. When the wall chiralities are uncorrelated, faceting is suppressed and outer layer corrugation, which is induced by the Moiré superlattice, is obtained in agreement with experiments. Finally, we offer an explanation for the higher incidence of faceting in multiwalled boron nitride nanotubes with respect to their carbon-based counterparts.
Collapse
Affiliation(s)
- Itai Leven
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Raymond and Beverly Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Roberto Guerra
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- CNR-IOM Democritos National Simulation Center, Via Bonomea 265, 34136 Trieste, Italy
| | - Andrea Vanossi
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- CNR-IOM Democritos National Simulation Center, Via Bonomea 265, 34136 Trieste, Italy
| | - Erio Tosatti
- International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy
- CNR-IOM Democritos National Simulation Center, Via Bonomea 265, 34136 Trieste, Italy
- The Abdus Salam International Centre for Theoretical Physics (ICTP), Strada Costiera 11, 34151 Trieste, Italy
| | - Oded Hod
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Raymond and Beverly Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
20
|
Koren E, Leven I, Lörtscher E, Knoll A, Hod O, Duerig U. Coherent commensurate electronic states at the interface between misoriented graphene layers. Nat Nanotechnol 2016; 11:752-7. [PMID: 27271963 DOI: 10.1038/nnano.2016.85] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2015] [Accepted: 04/26/2016] [Indexed: 05/13/2023]
Abstract
Graphene and layered materials in general exhibit rich physics and application potential owing to their exceptional electronic properties, which arise from the intricate π-orbital coupling and the symmetry breaking in twisted bilayer systems. Here, we report room-temperature experiments to study electrical transport across a bilayer graphene interface with a well-defined rotation angle between the layers that is controllable in situ. This twisted interface is artificially created in mesoscopic pillars made of highly oriented pyrolytic graphite by mechanical actuation. The overall measured angular dependence of the conductivity is consistent with a phonon-assisted transport mechanism that preserves the electron momentum of conduction electrons passing the interface. The most intriguing observations are sharp conductivity peaks at interlayer rotation angles of 21.8° and 38.2°. These angles correspond to a commensurate crystalline superstructure leading to a coherent two-dimensional (2D) electronic interface state. Such states, predicted by theory, form the basis for a new class of 2D weakly coupled bilayer systems with hitherto unexplored properties and applications.
Collapse
Affiliation(s)
- Elad Koren
- IBM Research - Zurich, Rueschlikon 8803, Switzerland
| | - Itai Leven
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | | | - Armin Knoll
- IBM Research - Zurich, Rueschlikon 8803, Switzerland
| | - Oded Hod
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Urs Duerig
- IBM Research - Zurich, Rueschlikon 8803, Switzerland
| |
Collapse
|
21
|
Affiliation(s)
- Itai Leven
- Department
of Physical Chemistry, School of Chemistry, The Raymond and Beverly
Sackler Faculty of Exact Sciences and The Sackler Center for Computational
Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Tal Maaravi
- Department
of Physical Chemistry, School of Chemistry, The Raymond and Beverly
Sackler Faculty of Exact Sciences and The Sackler Center for Computational
Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Ido Azuri
- Department
of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Leeor Kronik
- Department
of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Oded Hod
- Department
of Physical Chemistry, School of Chemistry, The Raymond and Beverly
Sackler Faculty of Exact Sciences and The Sackler Center for Computational
Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
22
|
|
23
|
Abstract
The sliding energy landscape of the heterogeneous graphene/h-BN interface is studied by means of the registry index. For a graphene flake sliding on top of h-BN, the anisotropy of the sliding energy corrugation with respect to the misfit angle between the two naturally mismatched lattices is found to reduce with the flake size. For sufficiently large flakes, the sliding energy corrugation is expected to be at least an order of magnitude lower than that obtained for matching lattices regardless of the relative interlayer orientation. Therefore, in contrast to the case of the homogeneous graphene interface where flake reorientations are known to eliminate superlubricty, here, a stable low-friction state is expected to occur. Our results mark heterogeneous layered interfaces as promising candidates for dry lubrication purposes.
Collapse
Affiliation(s)
- Itai Leven
- Department of Chemical Physics, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel
| | - Dana Krepel
- Department of Chemical Physics, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel
| | - Ortal Shemesh
- Department of Chemical Physics, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel
| | - Oded Hod
- Department of Chemical Physics, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel
| |
Collapse
|
24
|
Garel J, Leven I, Zhi C, Nagapriya KS, Popovitz-Biro R, Golberg D, Bando Y, Hod O, Joselevich E. Ultrahigh torsional stiffness and strength of boron nitride nanotubes. Nano Lett 2012; 12:6347-6352. [PMID: 23130892 DOI: 10.1021/nl303601d] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
We report the experimental and theoretical study of boron nitride nanotube (BNNT) torsional mechanics. We show that BNNTs exhibit a much stronger mechanical interlayer coupling than carbon nanotubes (CNTs). This feature makes BNNTs up to 1 order of magnitude stiffer and stronger than CNTs. We attribute this interlayer locking to the faceted nature of BNNTs, arising from the polarity of the B-N bond. This property makes BNNTs superior candidates to replace CNTs in nanoelectromechanical systems (NEMS), fibers, and nanocomposites.
Collapse
Affiliation(s)
- Jonathan Garel
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovot 76100, Israel
| | | | | | | | | | | | | | | | | |
Collapse
|