Related Experiment Video
Updated: May 29, 2026

Visualizing Hyporheic Flow Through Bedforms Using Dye Experiments and Simulation
Published on: November 18, 2015
NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface
Trond S Ingebrigtsen1, Søren Toxvaerd, Ole J Heilmann
1DNRF Centre Glass and Time, IMFUFA, Department of Sciences, Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark.
A new algorithm simulates particle motion on constant potential energy surfaces. This modified algorithm ensures stability and conserves energy, matching standard methods for accurate simulations of liquids.
Area of Science:
- Computational physics
- Molecular dynamics simulations
Background:
- Simulating particle dynamics requires accurate algorithms.
- Constant potential energy simulations are crucial for understanding liquid properties.
Purpose of the Study:
- To develop a stable algorithm for simulating geodesics on constant potential energy hypersurfaces.
- To address numerical errors and entropic drift in existing algorithms.
Main Methods:
- Discretizing geodesic stationarity conditions with Lagrangian multipliers.
- Implementing a modified algorithm for potential-energy and step-length conservation.
- Testing with single-precision simulations of Lennard-Jones liquid.
Main Results:
- The basic algorithm shows good stability with smoothed force cutoffs and precise initial conditions.
- A modified algorithm eliminates entropic drift and ensures absolute stability.
- The modified algorithm yields identical radial distribution functions to standard NVE algorithms.
Conclusions:
- The developed NVU algorithm provides a stable and accurate method for constant potential energy simulations.
- This advancement is crucial for reliable molecular dynamics studies of liquids.
- The algorithm's stability and accuracy are validated through analytical and simulation-based evidence.
Related Concept Videos
Curvilinear Motion: Normal and Tangential Components
The positive direction of the t-axis aligns with the increasing position of the car along the curved path, denoted by the unit vector ut. Simultaneously, the n-axis, perpendicular to the t-axis, dissects the curved path into differential arc segments, each forming the arc of a circle with a radius of...
Velocity Potential
Dynamics Of Circular Motion: Applications
Dynamics of Circular Motion
Any acceleration must be produced by some force. Therefore, any force or combination of forces can cause centripetal acceleration. A few examples include the tension in the rope on a...
Equations of Motion: Normal and Tangetial Components
Newton's second law of motion is employed to articulate the equation...
Curvilinear Motion: Rectangular Components
As the car advances, its position evolves over time. Quantifying the car's velocity involves computing the time...

