Abstract:Rigid bodies, made of smaller composite beads, are commonly used to simulate anisotropic particles with molecular dynamics or Monte Carlo methods. To accurately represent the particle shape and to obtain smooth and realistic effective pair interactions between two rigid bodies, each body may need to contain hundreds of spherical beads. Given an interacting pair of particles, traditional MD methods calculate the inter-body distances between the beads of the rigid bodies within a certain distance. For a system containing many anisotropic particles, distance calculations are computationally costly and limit the attainable system size and simulation time. However, the effective interaction between two rigid particles only depends on the distance between their center of masses and their relative orientation. Therefore, a function directly mapping the center of mass distance and orientation to the interaction energy between the two rigid bodies, would completely bypass inter-bead distance calculations. It is challenging to derive such a general function analytically for most non-spherical rigid bodies. We have trained neural nets, powerful tools to fit nonlinear functions to complex datasets, to achieve this task. The pair configuration is taken as input and the energy, forces and torques between two rigid particles are predicted directly. We show that MD simulations of cubes and cylinders performed with forces and torques obtained from the gradients of the energy neural-nets quantitatively match traditional simulations that uses composite rigid bodies. Both structural quantities and dynamic measures are in agreement, while achieving up to 23 times speed up over traditional molecular dynamics, depending on hardware and system size. The method can, in principle, be applied to any irregular shape with any pair interaction, provided that sufficient training data can be obtained.

Interaction of complex particles: A framework for the rapid and accurate approximation of pair potentials using neural networks

Molecular Dynamics Simulations of Anisotropic Particles Accelerated by Neural-Net Predicted Interactions

Bayesian Analysis Reveals the Key to Extracting Pair Potentials from Neutron Scattering Data

PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations.

Machine-learned coarse-grained potentials for particles with anisotropic shapes and interactions

Building nonparametric $n$-body force fields using Gaussian process regression

Inferring interaction potentials from stochastic particle trajectories

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A neural network potential based on pairwise resolved atomic forces and energies

Fast and Sample-Efficient Interatomic Neural Network Potentials for Molecules and Materials Based on Gaussian Moments

Long-range-corrected hybrid density functionals including random phase approximation correlation: Application to noncovalent interactions

Highly Accurate Description of Long-Range Interactions through the Combination of Neural Networks and Physical Models

Pair Approximating the Action For Molecular Rotations in Path Integral Monte Carlo

A RANDOM BATCH EWALD METHOD FOR PARTICLE SYSTEMS WITH COULOMB INTERACTIONS

Generation of Pairwise Potentials Using Multi-Dimensional Data Mining.

Effects of non-pairwise repulsion on nanoparticle assembly

Fast proper orthogonal descriptors for many-body interatomic potentials

Enhanced Deep Potential Model for Fast and Accurate Molecular Dynamics; Application to the Hydrated Electron

Modeling of Metal Nanoparticles: Development of Neural-Network Interatomic Potential Inspired by Features of the Modified Embedded-Atom Method

Effective Many-body Interactions in Reduced-Dimensionality Spaces Through Neural Network Models

Pair-distribution-function guided optimization of fingerprints for atom-centered neural network potentials.