An official website of the United States government
Abstract Machine learning interatomic potentials (MLIPs) are a promising technique for atomic modeling. While small errors are widely reported for MLIPs, an open concern is whether MLIPs can accurately reproduce atomistic dynamics and related physical properties in molecular dynamics (MD) simulations. In this study, we examine the state-of-the-art MLIPs and uncover several discrepancies related to atom dynamics, defects, and rare events (REs), compared to ab initio methods. We find that low averaged errors by current MLIP testing are insufficient, and develop quantitative metrics that better indicate the accurate prediction of atomic dynamics by MLIPs. The MLIPs optimized by the RE-based evaluation metrics are demonstrated to have improved prediction in multiple properties. The identified errors, the evaluation metrics, and the proposed process of developing such metrics are general to MLIPs, thus providing valuable guidance for future testing and improvements of accurate and reliable MLIPs for atomistic modeling.
more »
« less
Machine‐learning interatomic potentials for pyrolysis of polysiloxanes and properties of SiCO ceramics
Abstract We present machine‐learning interatomic potentials (MLIPs) for simulations of Si–C–O–H compounds. The MLIPs are constructed from moment tensor potentials (MTPs) and were trained to a library of configurations that included polysiloxane structures, hypothetical crystalline and amorphous SiCOH structures, and trajectories of Si–C–O–H systems obtained via ab initio molecular dynamic (aiMD) simulations at elevated temperatures. Passive, active, and hybrid learning strategies were implemented to develop the MLIPs. The MLIPs reproduce vibrational properties of polymers and SiCOH structures obtained from aiMD simulations, thus providing a tool to identify chemical units and distinct structural characteristics through their vibrational properties. Simulations of the polymer‐to‐ceramic transformation show the development of mixed tetrahedra in SiCO ceramics and align with experimental observations. Million‐atom simulations for several nanoseconds highlight the precipitation of graphitic nanosheets from a carbon‐rich SiCO precursor. Atomistic simulations with the MLIPs deliver details of chemical reaction mechanisms during the pyrolysis of polysiloxanes, including methane abstraction and Kumada‐like rearrangements that transform the siloxane backbone. While the MLIPs still leave room for systematic improvement, they deliver simulations with “density functional theory (DFT)‐like” quality at low and high temperatures.
more »
« less
- Award ID(s):
- 1743701
- PAR ID:
- 10549730
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Journal of the American Ceramic Society
- ISSN:
- 0002-7820
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Atomistic simulation has a broad range of applications from drug design to materials discovery. Machine learning interatomic potentials (MLIPs) have become an efficient alternative to computationally expensive ab initio simulations. For this reason, chemistry and materials science would greatly benefit from a general reactive MLIP, that is, an MLIP that is applicable to a broad range of reactive chemistry without the need for refitting. Here we develop a general reactive MLIP (ANI-1xnr) through automated sampling of condensed-phase reactions. ANI-1xnr is then applied to study five distinct systems: carbon solid-phase nucleation, graphene ring formation from acetylene, biofuel additives, combustion of methane and the spontaneous formation of glycine from early earth small molecules. In all studies, ANI-1xnr closely matches experiment (when available) and/or previous studies using traditional model chemistry methods. As such, ANI-1xnr proves to be a highly general reactive MLIP for C, H, N and O elements in the condensed phase, enabling high-throughput in silico reactive chemistry experimentation.more » « less
-
null (Ed.)The origin in deshielding of 29 Si NMR chemical shifts in R 3 Si–X, where X = H, OMe, Cl, OTf, [CH 6 B 11 X 6 ], toluene, and O X (O X = surface oxygen), as well as i Pr 3 Si + and Mes 3 Si + were studied using DFT methods. At the M06-L/6-31G(d,p) level of theory the geometry optimized structures agree well with those obtained experimentally. The trends in 29 Si NMR chemical shift also reproduce experimental trends; i Pr 3 Si–H has the most shielded 29 Si NMR chemical shift and free i Pr 3 Si + or isolable Mes 3 Si + have the most deshielded 29 Si NMR chemical shift. Natural localized molecular orbital (NLMO) analysis of the contributions to paramagnetic shielding ( σ p ) in these compounds shows that Si–R (R = alkyl, H) bonding orbitals are the major contributors to deshielding in this series. The Si–R bonding orbitals are coupled to the empty p-orbital in i Pr 3 Si + or Mes 3 Si + , or to the orbital in R 3 Si–X. This trend also applies to surface bound R 3 Si–O X . This model also explains chemical shift trends in recently isolated t Bu 2 SiH 2 + , t BuSiH 2 + , and SiH 3 + that show more shielded 29 Si NMR signals than R 3 Si + species. There is no correlation between isotropic 29 Si NMR chemical shift and charge at silicon.more » « less
-
Abstract The symmetry breaking in octahedral silsesquioxane and its germanium analogues (Si8O12H8and Ge8O12H8) has been investigated using the M06-2X/6-31++G(3df, 3pd) method and group theory. Both structures undergo$${O}_{h}\downarrow {T}_{h}$$symmetry breaking, characterized by pseudo-Jahn−Teller stabilization energies of 0.22 kcal/mol for Si-POSS and 9.82 kcal/mol for Ge-POSS. Under the influence of the pseudo-Jahn–Teller effect, the distortion vector involves the vibrational a2gmode with imaginary frequency. The distortion forces in Oh-POSS are predominantly localized on the oxygen atoms and driven by the coupling between the lowest unoccupied molecular orbital (a1g) and the highest occupied molecular orbital (a2g). The symmetry breaking is attributed to a pseudo-Jahn–Teller mechanism of type (a2gx a1g) = a2g. The symmetrical substitution of oxygen atoms by X (where X = C, N, P) results in viable Th-Si8X12H8and Th-Ge8X12H8structures. The observed pseudo-Jahn–Teller distortion and substitutional symmetry breaking caused by X indicates a consistent electronic relaxation mechanism, characterized by the formation of C=C, N=N and P=P bonds on the POSS cubic faces, which serves as hallmarks of stability. Additionally, we find that the volume of substituted Th-symmetrical POSS is sufficiently large to accommodate small ions.more » « less
-
Abstract This work systematically investigates the thermodynamic stability of SiaOb(M)cCdstructures derived from polymeric precursors incorporating metal fillers: Ta, Nb, and Hf, at 1200 and 1500°C. Structural characterization of the polymer derived ceramics (PDCs) employs X‐ray diffraction, Fourier transform infrared spectroscopy, and X‐ray photoelectron spectroscopy. Enthalpies of formation relative to crystalline components (metal oxide, silica, silicon carbide, and graphite) are obtained from thermodynamic measurements by high temperature oxide melt solution calorimetry. The enthalpies of formation (∆H°f, comp) of Ta‐1200, Hf‐1200, Nb‐1200, Ta‐1500, Hf‐1500, and Nb‐1500 specimens are −137.82 ± 9.72, −256.31 ± 8.97, −82.80 ± 9.82, −182.80 ± 7.85, −292.54 ± 9.38, −224.98 ± 9.60 kJ/mol, respectively. Overall incorporation of Hf results in most thermodynamically stable structures at all synthesis temperatures. SiaOb(M)cCdspecimens employing Nb fillers undergo the most stable structural evolution in this temperature range. The results indicate strong thermodynamic drive for carbothermal reduction of metal oxide domains. Incorporation of Ta provides the greatest stabilization of SiO3C mixed bonding environments. Ultimately, the choice of metal filler influences composition, structural evolution, and thermodynamic stability in PDCs.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1111/jace.20056
