2017/05/03 by Teruo Hirakawa, Teppei Suzuki, David R. Bowler +1
Physics and Astronomy · #cond-mat.mtrl-sci #physics.comp-ph
paper · pdf · doi:10.1088/1361-648x/aa810d
22 pages, 4 figures, submitted to JPCM
arxiv created 2017/05/03 · arxiv updated 2017/10/11
We discuss the development and implementation of a constant temperature (NVT) molecular dynamics scheme that combines the Nosé-Hoover chain thermostat with the extended Lagrangian Born-Oppenheimer molecular dynamics (BOMD) scheme, using a linear scaling density functional theory (DFT) approach. An integration scheme for this canonical-ensemble extended Lagrangian BOMD is developed and discussed in the context of the Liouville operator formulation. Linear scaling DFT canonical-ensemble extended Lagrangian BOMD simulations are tested on bulk silicon and silicon carbide systems to evaluate our integration scheme. The results show that the conserved quantity remains stable with no systematic drift even in the presence of the thermostat.