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.


Title: A neural network-assisted open boundary molecular dynamics simulation method
A neural network-assisted molecular dynamics method is developed to reduce the computational cost of open boundary simulations. Particle influxes and neural network-derived forces are applied at the boundaries of an open domain consisting of explicitly modeled Lennard-Jones atoms in order to represent the effects of the unmodeled surrounding fluid. Canonical ensemble simulations with periodic boundaries are used to train the neural network and to sample boundary fluxes. The method, as implemented in the LAMMPS, yields temperature, kinetic energy, potential energy, and pressure values within 2.5% of those calculated using periodic molecular dynamics and runs two orders of magnitude faster than a comparable grand canonical molecular dynamics system.  more » « less
Award ID(s):
1720530
PAR ID:
10414964
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
156
Issue:
18
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We present a new boundary condition for simulations of magnetic reconnection based on the topology of a Klein bottle. When applicable, the new condition is computationally cheaper than fully periodic boundary conditions, reconnects more flux than systems with conducting boundaries, and does not require assumptions about regions external to the simulation as is necessary for open boundaries. The new condition reproduces the expected features of reconnection but cannot be straightforwardly applied in systems with asymmetric upstream plasmas. 
    more » « less
  2. Many technologically useful materials are polycrystals composed of small monocrystalline grains that are separated by grain boundaries of crystallites with different lattice orientations. The energetics and connectivities of the grain boundaries play an essential role in defining the effective properties of materials across multiple scales. In this paper we derive a Fokker–Planck model for the evolution of the planar grain boundary network. The proposed model considers anisotropic grain boundary energy which depends on lattice misorientation and takes into account mobility of the triple junctions, as well as independent dynamics of the misorientations. We establish long time asymptotics of the Fokker–Planck solution, namely the joint probability density function of misorientations and triple junctions, and closely related the marginal probability density of misorientations. Moreover, for an equilibrium configuration of a boundary network, we derive explicit local algebraic relations, a generalized Herring Condition formula, as well as formula that connects grain boundary energy density with the geometry of the grain boundaries that share a triple junction. Although the stochastic model neglects the explicit interactions and correlations among triple junctions, the considered specific form of the noise, under the fluctuation–dissipation assumption, provides partial information about evolution of a grain boundary network, and is consistent with presented results of extensive grain growth simulations. 
    more » « less
  3. Motivated by biological transport phenomena that involve the motion of interacting molecular motors along linear filaments, we developed a theoretical framework to analyze the dynamics of interacting oligomers (extended size particles) on one-dimensional lattices. Our method extends the asymmetric simple exclusion processes for interacting monomers to particles of arbitrary size, and it utilizes cluster mean-field calculations supplemented by extensive Monte Carlo computer simulations. Interactions between particles are accounted for by a thermodynamically consistent method that views the formation and breaking bonds between particles as a chemical process. The dynamics of the system are analyzed for both periodic and open boundary conditions. It is found that the nature of the current-density relation depends on the strength of interactions, on the size of oligomers and on the way interactions influence particles transition rates. Stationary phase diagram is also fully evaluated, and it is shown how the dynamic properties depend on the interactions and on the sizes of the particles. To explain the dynamic behavior of the system particles density correlations are explicitly analyzed for different ranges of parameters. Theoretical calculations generally agree well with the results from the computer simulations, suggesting that our method correctly describes the main features of the molecular mechanisms of the transport of interacting oligomers. 
    more » « less
  4. null (Ed.)
    Machine learning (ML) is transforming all areas of science. The complex and time-consuming calculations in molecular simulations are particularly suitable for an ML revolution and have already been profoundly affected by the application of existing ML methods. Here we review recent ML methods for molecular simulation, with particular focus on (deep) neural networks for the prediction of quantum-mechanical energies and forces, on coarse-grained molecular dynamics, on the extraction of free energy surfaces and kinetics, and on generative network approaches to sample molecular equilibrium structures and compute thermodynamics. To explain these methods and illustrate open methodological problems, we review some important principles of molecular physics and describe how they can be incorporated into ML structures. Finally, we identify and describe a list of open challenges for the interface between ML and molecular simulation. 
    more » « less
  5. Machine learning (ML) is transforming all areas of science.The complex and time-consuming calculations in molecular simulations are particularly suitable for an ML revolution and have already been profoundly affected by the application of existing ML methods. Here we review recent ML methods for molecular simulation, with particular focus on (deep) neural networks for the prediction of quantum-mechanical energies and forces, on coarse-grained molecular dynamics, on the extraction of free energy surfaces and kinetics, and on generative network approaches to sample molecular equilibrium structures and compute thermodynamics. To explain these methods and illustrate open methodological problems,we review some important principles of molecular physics and describe how they can be incorporated into ML structures. Finally,we identify and describe a list of open challenges for the interface between ML and molecular simulation. 
    more » « less