skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Machine learning quantum-chemical bond scission in thermosets under extreme deformation
Despite growing interest in polymers under extreme conditions, most atomistic molecular dynamics simulations cannot describe the bond scission events underlying failure modes in polymer networks undergoing large strains. In this work, we propose a physics-based machine learning approach that can detect and perform bond breaking with near quantum-chemical accuracy on-the-fly in atomistic simulations. Particularly, we demonstrate that by coarse-graining highly correlated neighboring bonds, the prediction accuracy can be dramatically improved. By comparing with existing quantum mechanics/molecular mechanics methods, our approach is approximately two orders of magnitude more efficient and exhibits improved sensitivity toward rare bond breaking events at low strain. The proposed bond breaking molecular dynamics scheme enables fast and accurate modeling of strain hardening and material failure in polymer networks and can accelerate the design of polymeric materials under extreme conditions.  more » « less
Award ID(s):
2154916
PAR ID:
10498487
Author(s) / Creator(s):
;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
Applied Physics Letters
Volume:
122
Issue:
21
ISSN:
0003-6951
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. (, ArXivorg)
    null (Ed.)
    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. 
    more » « less
  2. ; ; ; (, Journal of the American Ceramic Society)
    Abstract Grain boundaries, ubiquitous in real materials, play an important role in the mechanical properties of ceramics. Using boron carbide as a typical superhard but brittle material under hypervelocity impact, we report atomistic reactive molecular dynamics simulations using the ReaxFF reactive force field fitted to quantum mechanics to examine grain‐boundary engineering strategies aimed at improving the mechanical properties. In particular, we examine the dynamical mechanical response of two grain‐boundary models with or without doped Si as a function of finite shear deformation. Our simulations show that doping Si into the grain boundary significantly increases the shear strength and stress threshold for amorphization and failure for both grain‐boundary structures. These results provide validation of our suggestions that Si doping provides a promising approach to mitigate amorphous band formation and failure in superhard boron carbide. 
    more » « less
  3. ; ; ; (, Physical Chemistry Chemical Physics)
    null (Ed.)
    We developed a new coarse-grained (CG) molecular dynamics force field for polyacrylamide (PAM) polymer based on fitting to the quantum mechanics (QM) equation of state (EOS). In this method, all nonbond interactions between representative beads are parameterized using a series of QM-EOS, which significantly improves the accuracy in comparison to common CG methods derived from atomistic molecular dynamics. This CG force-field has both higher accuracy and improved computational efficiency with respect to the OPLS atomistic force field. The nonbond components of the EOS were obtained from cold-compression curves on PAM crystals with rigid chains, while the covalent terms that contribute to the EOS were obtained using relaxed chains. For describing PAM gels we developed water–PAM interaction parameters using the same method. We demonstrate that the new CG-PAM force field reproduces the EOS of PAM crystals, isolated PAM chains, and water–PAM systems, while successfully predicting such experimental quantities as density, specific heat capacity, thermal conductivity and melting point. 
    more » « less
  4. ; (, Macromolecules)
    Understanding bond rupture in polymer networks remains a fundamental challenge due to the interplay of network topology and condensed phase effects. In this work, we introduce a predictive approach for determining bond rupture in thermosets based on shortest paths (SPs) analysis of the network topology. This method enumerates SP sets in networks with periodic boundary conditions, with applications to both all-atom and coarse-grained simulations. We find that bond rupture is most likely to initiate on the first (shortest) SP in the thermoset network and the strain at which the first bond ruptures is linearly correlated with the topological path length. As a result, one can predict the first bond rupture by computing the first SP directly from the initial thermoset topology without resorting to MD simulations. Furthermore, the specific bond rupture location along the first SP follows a probability distribution associated with the SP-based betweenness centrality. Subsequent bond rupture events are dictated by the instantaneous SP of partially broken networks. Moreover, we quantify the length scale dependence of SP distributions, introducing a means of partially bridging the observed ductile rupture in molecular simulations and brittle rupture in experiments. 
    more » « less
  5. ; ; ; (, Communications Chemistry)
    Abstract The potential energy landscape (PEL) formalism is a powerful tool within statistical mechanics to study the thermodynamic properties of classical low-temperature liquids and glasses. Recently, the PEL formalism has been extended to liquids/glasses that obey quantum mechanics, but applications have been limited to atomistic model liquids. In this work, we extend the PEL formalism to liquid/glassy water using path-integral molecular dynamics (PIMD) simulations, where nuclear quantum effects (NQE) are included. Our PIMD simulations, based on the q-TIP4P/F water model, show that the PEL of quantum water is both Gaussian and anharmonic. Importantly, the ring-polymers associated to the O/H atoms in the PIMD simulations, collapse at the local minima of the PEL (inherent structures, IS) for both liquid and glassy states. This allows us to calculate, analytically, the IS vibrational density of states (IS-VDOS) of the ring-polymer system using the IS-VDOS of classical water (obtained from classical MD simulations). The role of NQE on the structural properties of liquid/glassy water at various pressures are discussed in detail. Overall, our results demonstrate that the PEL formalism can effectively describe the behavior of molecular liquids at low temperatures and in the glass states, regardless of whether the liquid/glass obeys classical or quantum mechanics. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email