Verlet-like algorithms for Car-Parrinello molecular dynamics with unequal electronic occupations

Area of Science:

• Computational Chemistry
• Materials Science
• Quantum Mechanics

Background:

• Ab initio molecular dynamics (AIMD) simulations are crucial for understanding metallic, charged, and electrochemical systems.
• These systems often involve electronic states with unequal occupations, posing challenges for standard AIMD methods.
• Maintaining electronic orthonormality is essential for accurate dynamical simulations.

Purpose of the Study:

• To revisit and generalize the approach for handling fixed but arbitrary electronic occupations within the Car-Parrinello molecular dynamics (CPMD) scheme.
• To develop and validate algorithms for maintaining electronic orthonormality constraints in AIMD simulations with unequal occupations.
• To ensure accurate and energy-conserving dynamics for complex electronic systems.

Main Methods:

• Revisiting the general approach for fixed, arbitrary electronic occupations in Car-Parrinello molecular dynamics.
• Deriving expressions for constraining orbital velocities to satisfy orthonormality within velocity-Verlet integrators.
• Developing and comparing generalized unequal-occupation SHAKE and RATTLE algorithms.
• Validating the proposed algorithms using microcanonical ensemble simulations.

Main Results:

• The study presents a generalized unequal-occupation SHAKE algorithm for damped dynamics (energy optimization).
• A novel unequal-occupation RATTLE algorithm is derived and validated for microcanonical ensemble simulations.
• Proper orthogonalization methods are shown to be critical for correct state ordering and energy conservation.

Conclusions:

• The developed algorithms successfully maintain orthonormality constraints in AIMD simulations with unequally occupied electronic states.
• Accurate simulations of metallic, charged, and electrochemical systems are achievable with these generalized methods.
• Proper orthogonalization is key to achieving reliable and energy-conserving dynamics in complex electronic systems.