Related Experiment Video
Updated: May 18, 2026
05:56
Exploring Caspase Mutations and Post-Translational Modification by Molecular Modeling Approaches
Published on: October 13, 2022
Energy conserving, linear scaling Born-Oppenheimer molecular dynamics
M J Cawkwell1, Anders M N Niklasson
1Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA. cawkwell@lanl.gov
The Journal of Chemical Physics
|October 9, 2012
Summary
Researchers developed a new method for Born-Oppenheimer molecular dynamics simulations. This approach achieves linear scaling computational cost and long-term energy conservation, enabling more efficient and accurate simulations of molecular systems.
Area of Science:
- Computational Chemistry
- Molecular Dynamics
- Quantum Chemistry
Background:
- Born-Oppenheimer molecular dynamics (MD) is crucial for simulating molecular behavior.
- Traditional methods face computational challenges with increasing system size.
- Accurate energy conservation is vital for reliable simulation results.
Purpose of the Study:
- To develop a computationally efficient MD method.
- To achieve long-term conservation of total energy in simulations.
- To enable accurate simulations of larger molecular systems.
Main Methods:
- Implemented linear scaling computational cost using density matrix purification.
- Utilized sparse matrix algebra and a numerical threshold for efficiency.
- Employed the extended Lagrangian Born-Oppenheimer MD formalism.
Main Results:
- Simulations demonstrated simultaneous linear scaling and energy conservation.
- The method achieved a low pre-factor for computational cost.
- Generated microcanonical trajectories indistinguishable from exact force calculations over hundreds of picoseconds.
Conclusions:
- The developed method offers a significant advancement in computational efficiency for MD.
- Long-term energy conservation is maintained with approximate forces.
- This approach paves the way for more extensive and accurate molecular dynamics studies.
Related Concept Videos
01:21Molecular Kinetic Energy
The word "gas" comes from the Flemish word meaning "chaos," first used to describe vapors by the chemist J. B. van Helmont. Consider a container filled with gas, with a continuous and random motion of molecules. During collisions, the velocity component parallel to the wall is unchanged, and the component perpendicular to the wall reverses direction but does not change in magnitude. If the molecule’s velocity changes in the x-direction, then its momentum is changed. During the short time of the...
01:12Fermi Level Dynamics
The vacuum level denotes the energy threshold required for an electron to escape from a material surface. It is usually positioned above the conduction band of a semiconductor and acts as a benchmark for comparing electron energies within various materials.
Electron affinity in semiconductors refers to the energy gap between the minimum of its conduction band and the vacuum level and it is a critical parameter in determining how easily a semiconductor can accept additional electrons.
The work...
Electron affinity in semiconductors refers to the energy gap between the minimum of its conduction band and the vacuum level and it is a critical parameter in determining how easily a semiconductor can accept additional electrons.
The work...
01:28Conservation of Linear Momentum for a System of Particles
In the dynamic realm of billiards, a fascinating interplay of forces governs the motion of cue balls and stationary balls. When the cue ball collides with a stationary ball, linear momentum is exchanged. The cue ball imparts a fraction of its linear momentum to the stationary ball, causing the cue ball to decelerate while initiating the motion of the stationary ball.
The impulsive force at play during this interaction is of extremely short duration, rendering its impulse negligible. When...
The impulsive force at play during this interaction is of extremely short duration, rendering its impulse negligible. When...
01:27The Molecular Nature of Internal Energy
The internal energy of a molecule is determined by its degrees of freedom, including translational, rotational, and vibrational motions. In addition to these kinetic activities, the energy of molecules is also shaped by electronic energy, intermolecular forces, and the rest-mass energy of electrons and nuclei. These factors collectively influence the energy state of the molecules. The equipartition theorem of classical mechanics provides insight into this energy distribution. It posits that the...
01:26Scaling
In designing and analyzing filters, resonant circuits, or circuit analysis at large, working with standard element values like 1 ohm, 1 henry, or 1 farad can be convenient before scaling these values to more realistic figures. This approach is widely utilized by not employing realistic element values in numerous examples and problems; it simplifies mastering circuit analysis through convenient component values. The complexity of calculations is thereby reduced, with the understanding that...
01:16Energy Conservation and Bernoulli's Equation
Applying the conservation of energy principle or the work-energy theorem to an incompressible, inviscid fluid in laminar, steady, irrotational flow leads to Bernoulli's equation. It states that the sum of the fluid pressure, potential, and kinetic energy per unit volume is constant along a streamline.
All the terms in the equation have the dimension of energy per unit volume. The kinetic energy per unit volume is called the kinetic energy density, and the potential energy per unit volume is...
All the terms in the equation have the dimension of energy per unit volume. The kinetic energy per unit volume is called the kinetic energy density, and the potential energy per unit volume is...