Extended Lagrangian free energy molecular dynamics
Abstract
Extended free energy Lagrangians are proposed for first principles molecular dynamics simulations at finite electronic temperatures for plane-wave pseudopotential and local orbital density matrix-based calculations. Thanks to the extended Lagrangian description, the electronic degrees of freedom can be integrated by stable geometric schemes that conserve the free energy. For the local orbital representations both the nuclear and electronic forces have simple and numerically efficient expressions that are well suited for reduced complexity calculations. A rapidly converging recursive Fermi operator expansion method that does not require the calculation of eigenvalues and eigenfunctions for the construction of the fractionally occupied density matrix is discussed. An efficient expression for the Pulay force that is valid also for density matrices with fractional occupation occurring at finite electronic temperatures is also demonstrated.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- October 2011
- DOI:
- 10.1063/1.3656977
- arXiv:
- arXiv:1111.0261
- Bibcode:
- 2011JChPh.135p4111N
- Keywords:
-
- ab initio calculations;
- free energy;
- molecular dynamics method;
- pseudopotential methods;
- 31.15.xv;
- 02.70.Ns;
- 31.15.at;
- Molecular dynamics and other numerical methods;
- Molecular dynamics and particle methods;
- Molecule transport characteristics;
- molecular dynamics;
- electronic structure of polymers;
- Physics - Chemical Physics;
- Physics - Computational Physics
- E-Print:
- J. Chem. Phys. 135, 164111 (2011)