Related papers: Shadow molecular dynamics and atomic cluster expan…
Born-Oppenheimer dynamics is shown to provide an accurate approximation of time-independent Schr\"odinger observables for a molecular system with an electron spectral gap, in the limit of large ratio of nuclei and electron masses, without…
Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular…
We discuss the development and implementation of a constant temperature (NVT) molecular dynamics scheme that combines the Nos\'e-Hoover chain thermostat with the extended Lagrangian Born-Oppenheimer molecular dynamics (BOMD) scheme, using a…
A Lagrangian generalization of time-reversible Born-Oppenheimer molecular dynamics [Niklasson et al., Phys. Rev. Lett. vol. 97, 123001 (2006)] is proposed. The Lagrangian includes extended electronic degrees of freedom as auxiliary…
Constant potential method molecular dynamics simulation (CPM MD) enables the accurate modelling of atomistic electrode charges when studying the electrode-electrolyte interface at the nanoscale. Here we extend the theoretical framework of…
Variable charge models (e.g., EEM, QEq, ES+) in reactive molecular dynamics simulations often inherently impose a global charge transfer between atoms (approximating each system as ideal metal). Consequently, most surface processes (e.g.,…
Quantum computing has the potential to reduce the computational cost required for quantum dynamics simulations. However, existing quantum algorithms for coupled electron-nuclear dynamics simulation either require fault-tolerant devices, or…
The Non-Markovian Stochastic Schrodinger Equation (NMSSE) offers a promising approach for open quantum simulations, especially in large systems, owing to its low scaling complexity and suitability for parallel computing. However, its…
We present an improved field-theoretic approach to the grand-canonical potential suitable for linear scaling molecular dynamics simulations using forces from self-consistent electronic structure calculations. It is based on an exact…
Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is…
Extended Lagrangian Born-Oppenheimer molecular dynamics based on Kohn-Sham density functional theory is generalized in the limit of vanishing self-consistent field optimization prior to the force evaluations. The equations of motion are…
We introduce a scheme for molecular simulations, the Deep Potential Molecular Dynamics (DeePMD) method, based on a many-body potential and interatomic forces generated by a carefully crafted deep neural network trained with ab initio data.…
Ab initio Born-Oppenheimer molecular dynamics (AIMD) is a valuable method for simulating physico-chemical processes of complex systems, including reactive systems, and for training machine learning models and force fields. Speed and…
We introduce a combined molecular dynamics (MD) and quantum trajectories (QT) code to simulate the effects of near-resonant optical fields on state-vector evolution and particle motion in a collisional system. In contrast to collisionless…
We present a first principles molecular dynamics approach that is based on time-reversible ex- tended Lagrangian Born-Oppenheimer molecular dynamics [Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field…
Nuclear quantum effects such as zero-point energy and hydrogen tunnelling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei…
Accurate modelling of electrostatic interactions and charge transfer is fundamental to computational chemistry, yet most machine learning interatomic potentials (MLIPs) rely on local atomic descriptors that cannot capture long-range…
The computational bottleneck of molecular dynamics is the pairwise additive long-range interactions between particles. The random batch Ewald (RBE) method provides a highly efficient and superscalable solver for long-range interactions, but…
The force field describing the calculated interaction between atoms or molecules is the key to the accuracy of many molecular dynamics (MD) simulation results. Compared with traditional or semi-empirical force fields, machine learning force…
We develop a random batch Ewald (RBE) method for molecular dynamics simulations of particle systems with long-range Coulomb interactions, which achieves an $O(N)$ complexity in each step of simulating the $N$-body systems. The RBE method is…