NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface

Abstract

An algorithm is derived for computer simulation of geodesics on the constant potential-energy hypersurface of a system of N classical particles. First, a basic time-reversible geodesic algorithm is derived by discretizing the geodesic stationarity condition and implementing the constant potential energy constraint via standard Lagrangian multipliers. The basic NVU algorithm is tested by single-precision computer simulations of the Lennard-Jones liquid. Excellent numerical stability is obtained if the force cutoff is smoothed and the two initial configurations have identical potential energy within machine precision. Nevertheless, just as for NVE algorithms, stabilizers are needed for very long runs in order to compensate for the accumulation of numerical errors that eventually lead to "entropic drift" of the potential energy towards higher values. A modification of the basic NVU algorithm is introduced that ensures potential-energy and step-length conservation; center-of-mass drift is also eliminated. Analytical arguments confirmed by simulations demonstrate that the modified NVU algorithm is absolutely stable. Finally, simulations show that the NVU algorithm and the standard leap-frog NVE algorithm have identical radial distribution functions for the Lennard-Jones liquid.

References in corpus (14)

• Pressure-energy correlations in liquids. IV. "Isomorphs" in liquid state diagrams
• "Isomorphs" in liquid state diagrams
• Pressure-energy correlations in liquids. I. Results from computer simulations
• Strong pressure-energy correlations in van der Waals liquids
• Pressure-energy correlations in liquids. II. Analysis and consequences
• Pressure-energy correlations in liquids. III. Statistical mechanics and thermodynamics of liquids with hidden scale invariance
• Pressure-energy correlations in liquids. V. Isomorphs in generalized Lennard-Jones systems
• A repulsive reference potential reproducing the dynamics of a liquid with attractions
• Feasibility of single-order parameter description of equilibrium viscous liquid dynamics
• Predicting The Effective Temperaure of a Glass
• Strongly correlating liquids and their isomorphs
• Global perspectives on the energy landscapes of liquids, supercooled liquids, and glassy systems: The potential energy landscape ensemble
• Global perspectives on the energy landscapes of liquids, supercooled liquids, and glassy systems: Geodesic pathways through the potential energy landscape
• NVU dynamics. II. Comparing to four other dynamics

Cited by in corpus (15)

• Simplicity of condensed matter at its core: Generic definition of a Roskilde-simple system
• RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles
• Isomorphs, hidden scale invariance, and quasiuniversality
• Impact of size polydispersity on the nature of Lennard-Jones liquids
• The EXP pair-potential system. II. Fluid phase isomorphs
• NVU perspective on simple liquids' quasiuniversality
• Effect of Energy Polydispersity on the Nature of Lennard-Jones Liquids
• The EXP pair-potential system. I. Fluid phase isotherms, isochores, and quasiuniversality
• Isomorph theory of physical aging
• NVU dynamics. II. Comparing to four other dynamics
• NVU dynamics. III. Simulating molecules at constant potential energy
• Discrete dynamics versus analytic dynamics
• $NVU$ view on energy polydisperse Lennard-Jones systems
• Interpolating between pair-potential systems
• Viscous liquid dynamics modeled as random walks within overlapping hyperspheres