-
Atomic-scale identification of the active sites of nanocatalysts
Authors:
Yao Yang,
Jihan Zhou,
Zipeng Zhao,
Geng Sun,
Saman Moniri,
Colin Ophus,
Yongsoo Yang,
Ziyang Wei,
Yakun Yuan,
Cheng Zhu,
Qiang Sun,
Qingying Jia,
Hendrik Heinz,
Jim Ciston,
Peter Ercius,
Philippe Sautet,
Yu Huang,
Jianwei Miao
Abstract:
Alloy nanocatalysts have found broad applications ranging from fuel cells to catalytic converters and hydrogenation reactions. Despite extensive studies, identifying the active sites of nanocatalysts remains a major challenge due to the heterogeneity of the local atomic environment. Here, we advance atomic electron tomography to determine the 3D local atomic structure, surface morphology and chemi…
▽ More
Alloy nanocatalysts have found broad applications ranging from fuel cells to catalytic converters and hydrogenation reactions. Despite extensive studies, identifying the active sites of nanocatalysts remains a major challenge due to the heterogeneity of the local atomic environment. Here, we advance atomic electron tomography to determine the 3D local atomic structure, surface morphology and chemical composition of PtNi and Mo-doped PtNi nanocatalysts. Using machine learning trained by density functional theory calculations, we identify the catalytic active sites for the oxygen reduction reaction from experimental 3D atomic coordinates, which are corroborated by electrochemical measurements. By quantifying the structure-activity relationship, we discover a local environment descriptor to explain and predict the catalytic active sites at the atomic level. The ability to determine the 3D atomic structure and chemical species coupled with machine learning is expected to expand our fundamental understanding of a wide range of nanocatalysts.
△ Less
Submitted 18 February, 2022;
originally announced February 2022.
-
Prediction of Carbon Nanostructure Mechanical Properties and Role of Defects Using Machine Learning
Authors:
Qi Zhao,
Jordan J. Winetrout,
Yanxun Xu,
Yusu Wang,
Hendrik Heinz
Abstract:
Carbon fiber and graphene-based nanostructures such as carbon nanotubes (CNTs) and defective structures have extraordinary potential as strong and lightweight materials. A longstanding bottleneck has been lack of understanding and implementation of atomic-scale engineering to harness the theoretical limits of modulus and tensile strength, of which only a fraction is routinely reached today. Here w…
▽ More
Carbon fiber and graphene-based nanostructures such as carbon nanotubes (CNTs) and defective structures have extraordinary potential as strong and lightweight materials. A longstanding bottleneck has been lack of understanding and implementation of atomic-scale engineering to harness the theoretical limits of modulus and tensile strength, of which only a fraction is routinely reached today. Here we demonstrate accurate and fast predictions of mechanical properties for CNTs and arbitrary 3D graphitic assemblies based on a training set of over 1000 stress-strain curves from cutting-edge reactive MD simulation and machine learning (ML). Several ML methods are compared and show that our newly proposed hierarchically structured graph neural networks with spatial information (HS-GNNs) achieve predictions in modulus and strength for any 3D nanostructure with only 5-10% error across a wide range of possible values. The reliability is sufficient for practical applications and a great improvement over off-the shelf ML methods with up to 50% deviation, as well as over earlier models for specific chemistry with 20% deviation. The algorithms allow more than 10 times faster mechanical property predictions than traditional molecular dynamics simulations, the identification of the role of defects and random 3D morphology, and high-throughput screening of 3D structures for enhanced mechanical properties. The algorithms can be scaled to morphologies up to 100 nm in size, expanded for chemically similar compounds, and trained to predict a broader range of properties.
△ Less
Submitted 23 October, 2021; v1 submitted 1 October, 2021;
originally announced October 2021.
-
Implementing Reactivity in Molecular Dynamics Simulations with the Interface Force Field (IFF-R) and Other Harmonic Force Fields
Authors:
Jordan J. Winetrout,
Krishan Kanhaiya,
Geeta Sachdeva,
Ravindra Pandey,
Behzad Damirchi,
Adri van Duin,
Gregory Odegard,
Hendrik Heinz
Abstract:
The Interface force field (IFF) enables accurate simulations of bulk and interfacial properties of compounds and multiphase materials. However, the simulation of reactions and mechanical properties up to failure remains challenging and expensive. Here we introduce the Reactive Interface Force Field (IFF-R) to analyze bond breaking and failure of complex materials using molecular dynamics simulatio…
▽ More
The Interface force field (IFF) enables accurate simulations of bulk and interfacial properties of compounds and multiphase materials. However, the simulation of reactions and mechanical properties up to failure remains challenging and expensive. Here we introduce the Reactive Interface Force Field (IFF-R) to analyze bond breaking and failure of complex materials using molecular dynamics simulations. IFF-R uses a Morse potential instead of a harmonic potential as typically employed in molecular dynamics force fields to describe the bond energy, which can render any desired bond reactive by specification of the curve shape of the potential energy and the bond dissociation energy. This facile extension of IFF and other force fields that utilize a harmonic bond energy term allows the description of bond breaking without loss in functionality, accuracy, and speed. The method enables quantitative, on-the-fly computations of bond breaking and stress-strain curves up to failure in any material. We illustrate accurate predictions of mechanical behavior for a variety of material systems, including metals (iron), ceramics (carbon nanotubes), polymers (polyacrylonitrile and cellulose I\b{eta}), and include sample parameters for common bonds based on using experimental and high-level (MP2) quantum mechanical reference data. Computed structures, surface energies, elastic moduli, and tensile strengths are in excellent agreement with available experimental data. Non-reactive properties are shown to be essentially identical to IFF values. Computations are approximately 50 times faster than using ReaxFF and require only a single set of parameters. Compatibility of IFF and IFF-R with biomolecular force fields allows the quantitative analysis of the mechanics of proteins, DNA, and other biological molecules.
△ Less
Submitted 29 July, 2021;
originally announced July 2021.
-
Molecular Dynamics Modeling of Epoxy Resins using the Reactive Interface Force Field
Authors:
G. M. Odegard,
S. U. Patil,
P. Deshpande,
K. Kanhaiya,
J. Winetrout,
H. Heinz,
S. Shah,
M. Maiaru
Abstract:
Predictive computational modeling of polymer materials is necessary for the efficient design of composite materials and the corresponding processing methods. Molecular dynamics (MD) modeling is especially important for establishing accurate processing-structure-property relationships for neat resins. For MD modeling of amorphous polymer materials, an accurate force field is fundamental to reliable…
▽ More
Predictive computational modeling of polymer materials is necessary for the efficient design of composite materials and the corresponding processing methods. Molecular dynamics (MD) modeling is especially important for establishing accurate processing-structure-property relationships for neat resins. For MD modeling of amorphous polymer materials, an accurate force field is fundamental to reliable prediction of properties. Reactive force fields, in which chemical bonds can be formed or broken, offer further capability in predicting the mechanical behavior of amorphous polymers subjected to relatively large deformations. To this end, the Reactive Interface Force Field (IFF-R) has been recently developed to provide an efficient means to predict the behavior of materials under these conditions. Although IFF-R has been proven to be accurate for some crystalline organic and inorganic systems, it has not yet been proven to be accurate for amorphous polymer systems. The objective of this study is to use IFF-R to predict the thermo-mechanical properties of three different epoxy systems and validate with experimental measurements. The results indicate that IFF-R predicts thermo-mechanical properties that agree closely with experiment. Therefore, IFF-R can be used to reliably establish mechanical properties of polymers on the molecular level for future design of new composite materials and processing methods.
△ Less
Submitted 21 November, 2021; v1 submitted 29 July, 2021;
originally announced July 2021.
-
Capturing Nucleation at 4D Atomic Resolution
Authors:
Jihan Zhou,
Yongsoo Yang,
Yao Yang,
Dennis S. Kim,
Andrew Yuan,
Xuezeng Tian,
Colin Ophus,
Fan Sun,
Andreas K. Schmid,
Michael Nathanson,
Hendrik Heinz,
Qi An,
Hao Zeng,
Peter Ercius,
Jianwei Miao
Abstract:
Nucleation plays a critical role in many physical and biological phenomena ranging from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases. However, nucleation is a challenging process to study in experiments especially in the early stage when several atoms/molecules start to form a new phase from its parent phase. Here, we advance…
▽ More
Nucleation plays a critical role in many physical and biological phenomena ranging from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases. However, nucleation is a challenging process to study in experiments especially in the early stage when several atoms/molecules start to form a new phase from its parent phase. Here, we advance atomic electron tomography to study early stage nucleation at 4D atomic resolution. Using FePt nanoparticles as a model system, we reveal that early stage nuclei are irregularly shaped, each has a core of one to few atoms with the maximum order parameter, and the order parameter gradient points from the core to the boundary of the nucleus. We capture the structure and dynamics of the same nuclei undergoing growth, fluctuation, dissolution, merging and/or division, which are regulated by the order parameter distribution and its gradient. These experimental observations differ from classical nucleation theory (CNT) and to explain them we propose the order parameter gradient (OPG) model. We show the OPG model generalizes CNT and energetically favours diffuse interfaces for small nuclei and sharp interfaces for large nuclei. We further corroborate this model using molecular dynamics simulations of heterogeneous and homogeneous nucleation in liquid-solid phase transitions of Pt. We anticipate that the OPG model is applicable to different nucleation processes and our experimental method opens the door to study the structure and dynamics of materials with 4D atomic resolution.
△ Less
Submitted 28 February, 2019; v1 submitted 27 July, 2018;
originally announced July 2018.
-
Accurate Force Field Parameters and pH Resolved Surface Models for Hydroxyapatite to Understand Structure, Mechanics, Hydration, and Biological Interfaces
Authors:
Tzu-Jen,
Hendrik Heinz
Abstract:
Mineralization of bone and teeth involves interactions between biomolecules and hydroxyapatite. Associated complex interfaces and processes remain difficult to analyze at the 1 to 100 nm scale using current laboratory techniques, and prior models for atomistic simulations are limited in the representation of chemical bonding, surface chemistry, and interfacial interactions. This work introduces an…
▽ More
Mineralization of bone and teeth involves interactions between biomolecules and hydroxyapatite. Associated complex interfaces and processes remain difficult to analyze at the 1 to 100 nm scale using current laboratory techniques, and prior models for atomistic simulations are limited in the representation of chemical bonding, surface chemistry, and interfacial interactions. This work introduces an accurate force field along with pH-resolved surface models for hydroxyapatite to represent chemical bonding, structural, surface, interfacial, and mechanical properties in quantitative agreement with experiment. The accuracy is orders of magnitude higher in comparison to earlier models to facilitate quantitative monitoring of inorganic-biological assembly. The force field is integrated into the CHARMM, AMBER, OPLS-AA, PCFF, and INTERFACE force fields to enable realistic simulations of apatite-biological systems of any composition and ionic strength. Specifically, the parameters reproduce lattice constants (<0.5% deviation), IR spectrum, cleavage energies, immersion energies in water (<5% deviation), and elastic constants (<10% deviation) of hydroxyapatite in comparison to experiment. Interactions between mineral, water, and organic compounds are represented by standard combination rules in the force field without additional adjustable parameters and shown to achieve quantitative accuracy. Surface models for common (001), (010), (020), (101) facets and nanocrystals are introduced as a function of pH on the basis of extensive experimental data. New insight into surface and immersion energies, the structure of aqueous interfaces, density profiles, and superficial dissolution is described. Mechanisms of specific binding of peptides, drugs, and mineralization can be analyzed and the force field is extensible to substituted and defective apatites as well as to other calcium phosphate phases.
△ Less
Submitted 30 November, 2015;
originally announced December 2015.
-
Three-dimensional coordinates of individual atoms in materials revealed by electron tomography
Authors:
Rui Xu,
Chien-Chun Chen,
Li Wu,
M. C. Scott,
W. Theis,
Colin Ophus,
Matthias Bartels,
Yongsoo Yang,
Hadi Ramezani-Dakhel,
Michael R. Sawaya,
Hendrik Heinz,
Laurence D. Marks,
Peter Ercius,
Jianwei Miao
Abstract:
Crystallography, the primary method for determining the three-dimensional (3D) atomic positions in crystals, has been fundamental to the development of many fields of science. However, the atomic positions obtained from crystallography represent a global average of many unit cells in a crystal. Here, we report, for the first time, the determination of the 3D coordinates of thousands of individual…
▽ More
Crystallography, the primary method for determining the three-dimensional (3D) atomic positions in crystals, has been fundamental to the development of many fields of science. However, the atomic positions obtained from crystallography represent a global average of many unit cells in a crystal. Here, we report, for the first time, the determination of the 3D coordinates of thousands of individual atoms and a point defect in a material by electron tomography with a precision of ~19 picometers, where the crystallinity of the material is not assumed. From the coordinates of these individual atoms, we measure the atomic displacement field and the full strain tensor with a 3D resolution of ~1nm^3 and a precision of ~10^-3, which are further verified by density functional theory calculations and molecular dynamics simulations. The ability to precisely localize the 3D coordinates of individual atoms in materials without assuming crystallinity is expected to find important applications in materials science, nanoscience, physics and chemistry.
△ Less
Submitted 21 May, 2015;
originally announced May 2015.
-
Structure and Phase Transitions of Alkyl Chains on Mica
Authors:
Hendrik Heinz,
Hein J. Castelijns,
Ulrich W. Suter
Abstract:
We use molecular dynamics as a tool to understand the structure and phase transitions [Osman et. al. J. Phys. Chem. B 2000, 104, 4433; 2002, 106, 653] in alkylammonium micas. The consistent force field 91 is extended for accurate simulation of mica and related minerals. We investigate mica sheets with 12 octadecyltrimethylammonium (C18) ions or 12 dioctadecyldimethylammonium (2C18) ions, respect…
▽ More
We use molecular dynamics as a tool to understand the structure and phase transitions [Osman et. al. J. Phys. Chem. B 2000, 104, 4433; 2002, 106, 653] in alkylammonium micas. The consistent force field 91 is extended for accurate simulation of mica and related minerals. We investigate mica sheets with 12 octadecyltrimethylammonium (C18) ions or 12 dioctadecyldimethylammonium (2C18) ions, respectively, as single and layered structures at different temperatures with periodicity in the xy plane by NVT dynamics. The alkylammonium ions reside preferably above the cavities in the mica surface with an aluminum-rich boundary. The nitrogen atoms are 380 to 390 pm distant to the superficial silicon-aluminum plane. With increasing temperature, rearrangements of C18 ions on the mica surface are found, while 2C18 ions remain tethered due to geometric restraints. We present basal-plane spacings in the duplicate structures, tilt angles of the alkyl chains, and gauche-trans ratios to analyze the chain conformation. Also, the individual phase transitions of the two systems on heating are explained. Where experimental data are available, the agreement is very good. We propose a geometric parameter lamba for the saturation of the surface with alkyl chains, which determines the preferred self-assembly pattern, i.e., islands, intermediate, or continuous. Lambda also determines the tilt angles in continuous layers on mica or other surfaces. The thermal decomposition appears to be a Hofmann elimination with mica as a base-template.
△ Less
Submitted 24 November, 2003;
originally announced November 2003.
-
Surface Structure of Organoclays
Authors:
Hendrik Heinz,
Ulrich W. Suter
Abstract:
The properties of organically modified clay minerals determine essentially the quality of polymer-clay nanocomposites. We investigate alkylammonium-micas with different alkyl chain length by molecular dynamics simulation as homogeneous mixtures and separated island structures on the mica surface (80 percent alkali exchange, 20 percent of the alkali ions remain), including a detailed model of the…
▽ More
The properties of organically modified clay minerals determine essentially the quality of polymer-clay nanocomposites. We investigate alkylammonium-micas with different alkyl chain length by molecular dynamics simulation as homogeneous mixtures and separated island structures on the mica surface (80 percent alkali exchange, 20 percent of the alkali ions remain), including a detailed model of the Si --> Al...K charge defects on the surface of phyllosilicates. By comparison with experiment, we find that long chains (C18) lead to a mixed phase of alkali ions and organic ions on the surface, while chains of a medium length (~C12) give rise to phase separation. Very short chains, by thermodynamical arguments, prefer again a homogeneously mixed surface. We explain this interesting observation with a basic free energy model and present a diagram of the phase behaviour as a function of surface saturation with alkyl chains and the chain length. Besides, an order-disorder transition of the tethered alkyl chains is found on heating when ~20 percent of the torsional angles along the hydrocarbon backbones are gauche. At higher temperature, also rearrangements of the ammonium headgroups on the surface are possible, which lead to metastable structures after quenching.
△ Less
Submitted 24 November, 2003;
originally announced November 2003.
-
Analysis of the Phase Transitions in Alkyl- Mica by Density and Pressure Profiles
Authors:
Hendrik Heinz,
Wolfgang Paul,
Ulrich W. Suter,
Kurt Binder
Abstract:
In a previous work [H. Heinz, H. J. Castelijns, U. W. Suter, J. Am. Chem. Soc. 115, 9500 (2003)], we developed an accurate force field and simulated the phase transitions in C18-mica (octadecyltrimethylammonium-mica) as well as the absence of such transitions in 2C18-mica (dioctadecyldimethylammonium-mica) between room temperature and 100 deg C. Here we analyze (i) average z coordinates of the c…
▽ More
In a previous work [H. Heinz, H. J. Castelijns, U. W. Suter, J. Am. Chem. Soc. 115, 9500 (2003)], we developed an accurate force field and simulated the phase transitions in C18-mica (octadecyltrimethylammonium-mica) as well as the absence of such transitions in 2C18-mica (dioctadecyldimethylammonium-mica) between room temperature and 100 deg C. Here we analyze (i) average z coordinates of the carbon atoms and interdigitation of the hydrocarbon bilayers, (ii) density profiles and (iii) pressure profiles of the structures along all Cartesian axes. In C18-mica, the standard deviation in the z coordinate for the chain atoms is high and more than doubles in the disordered phase. The order-disorder transition is accompanied by a change in the orientation of the ammonium headgroup, as well as decreasing tensile and shear stress in the disordered phase. In 2C18-mica, the standard deviation in the z coordinate for the chain atoms is low and does not increase remarkably on heating. The backbones display a highly regular structure, which is slightly obscured by rotations in the C18 backbones and minor headgroup displacements at 100 deg C. Close contacts between the bulky headgroups with sidearms cause significant local pressure which is in part not relieved at 100 deg C. An increase of the basal-plane spacing at higher temperature is found in both systems due to larger separation between the two hydrocarbon layers and an increased z spacing between adjacent chain atoms (=decreased tilt of the chains relative to the surface normal), and, in C18-mica only, a stronger upward orientation of the C18 chain at the ammonium headgroup. The likelihood for chain interdigitation between the two hydrocarbon layers is 0.24 to 0.30 for C18-mica,and 0.65 to 0.26 for 2C18-mica (for 20 deg C to 100 deg C).
△ Less
Submitted 23 November, 2003;
originally announced November 2003.
-
Calculation of local pressure tensors in systems with many-body interactions
Authors:
Hendrik Heinz,
Wolfgang Paul,
Kurt Binder
Abstract:
Local pressures are important in the calculation of interface tensions and in analyzing micromechanical behavior. The calculation of local pressures in computer simulations has been limited to systems with pairwise interactions between the particles, which is not sufficient for chemically detailed systems with many-body potentials such as angles and torsions. We introduce a method to calculate l…
▽ More
Local pressures are important in the calculation of interface tensions and in analyzing micromechanical behavior. The calculation of local pressures in computer simulations has been limited to systems with pairwise interactions between the particles, which is not sufficient for chemically detailed systems with many-body potentials such as angles and torsions. We introduce a method to calculate local pressures in systems with n-body interactions (n=2,3,4, . . .) based on a micromechanical definition of the pressure tensor. The local pressure consists of a kinetic contribution from the linear momentum of the particles and an internal contribution from dissected many-body interactions by infinitesimal areas. To define dissection by a small area, respective n-body interactions are divided into two geometric centers, effectively reducing them to two-body interactions. Consistency with hydrodynamics-derived formulas for systems with two-body interactions (J. H. Irving and J. G. Kirkwood, J. Chem. Phys. 18, 817 (1950)), for average cross-sectional pressures (B. D. Todd, D. J. Evans, and P. J. Daivis, Phys. Rev. E 52, 1627 (1995)), and for volume averaged pressures (virial formula) is shown. As a simple numerical example, we discuss liquid propane in a cubic box. Local, crosssectional,and volume-averaged pressures as well as relative contributions from two-body and three-body forces are analyzed with the proposed method, showing full numerical equivalence with the existing approaches. The method allows computing local pressures in the presence of many-body interactions in atomistic simulations of complex materials and biological systems.
△ Less
Submitted 24 February, 2006; v1 submitted 31 August, 2003;
originally announced September 2003.