How well do Car-Parrinello simulations reproduce the Born-Oppenheimer surface ? Theory and Examples

Abstract

We derive an analytic expression for the average difference between the forces on the ions in a Car-Parrinello simulation and the forces obtained at the same ionic positions when the electrons are at their ground state. We show that for common values of the fictitious electron mass, a systematic bias may affect the Car-Parrinello forces in systems where the electron-ion coupling is large. We show that in the limit where the electronic orbitals are rigidly dragged by the ions the difference between the two dynamics amounts to a rescaling of the ionic masses, thereby leaving the thermodynamics intact. We study the examples of crystalline magnesium oxide and crystalline and molten silicon. We find that for crystalline silicon the errors are very small. For crystalline MgO the errors are very large but the dynamics can be quite well corrected within the rigid-ion model. We conclude that it is important to control the effect of the electron mass parameter on the quantities extracted from Car-Parrinello simulations.

Submitted to the Journal of Chemical Physics

Cited by in corpus (18)

• Towards an assessment of the accuracy of density functional theory for first principles simulations of water II
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• Ab-Initio Molecular Dynamics
• Time-reversible Born-Oppenheimer molecular dynamics
• Free energy of liquid water on the basis of quasi-chemical theory and ab initio molecular dynamics
• A unified electrostatic and cavitation model for first-principles molecular dynamics in solution
• Coupled Electron Ion Monte Carlo Calculations of Dense Metallic Hydrogen
• Interatomic potentials: Achievements and challenges
• The hydration state of HO${}^{-}$(aq)
• On the theory underlying the Car-Parrinello method and the role of the fictitious mass parameter
• A many-body interatomic potential for ionic systems: application to MgO
• Next generation extended Lagrangian first principles molecular dynamics
• Liquid Water through Density-Functional Molecular Dynamics: Plane-Wave vs Atomic-Orbital Basis Sets
• Analysis of the Time Reversible Born-Oppenheimer Molecular Dynamics
• Second generation wave-function thermostat for ab-initio molecular dynamics
• Melting slope of MgO from molecular dynamics and density functional theory
• A Quantum-Mechanics Molecular-Mechanics scheme for extended systems
• Ultrafast polarization switching in BaTiO${}_3$ by photoactivation of its ferroelectric and central modes