Parts | Days | Selection | Search | Updates | Downloads | Help
HL: Fachverband Halbleiterphysik
HL 33: Frontiers of Electronic Structure Theory: Nuclear Dynamics, Methods
HL 33.4: Talk
Tuesday, March 17, 2015, 11:30–11:45, MA 004
Converged Nuclear Quantum Statistics from Semi-Classical Path Integrals — •Igor Poltavskyi and Alexandre Tkatchenko — Fritz-Haber-Institut der Max-Planck-Gesellschaft, Berlin, Germany
The quantum nature of nuclear motions plays a vital role in the structure, stability, and thermodynamics of molecular systems. The standard approach to take nuclear quantum effects (NQE) into account is the Feynman-Kac imaginary-time path-integral molecular dynamics (PIMD). Conventional PIMD simulations require exceedingly large number of classical subsystems (beads) to accurately capture NQE, resulting in considerable computational cost even at room temperature due to the rather high internal vibrational frequencies of many molecules of interest.
We propose a novel parameter-free form for the PI partition function and estimators to calculate converged thermodynamic averages. Our approach requires the same ingredients as the conventional PIMD simulations, but decreases the number of required beads by roughly an order of magnitude. This greatly extends the applicability of ab initio PIMD for realistic molecular systems. The developed method has been applied to study the thermodynamics of N2, H2O, CO2, and C6H6 molecules. For all of the considered systems at room temperature, 4 to 8 beads are enough to recover the NQE contribution to the total energy within 2% of the fully converged quantum result.