Accelerated path integral methods for atomistic simulations at ultralow temperatures
Abstract
Path integral methods provide a rigorous and systematically convergent framework to include the quantum mechanical nature of atomic nuclei in the evaluation of the equilibrium properties of molecules, liquids, or solids at finite temperature. Such nuclear quantum effects are often significant for light nuclei already at room temperature, but become crucial at cryogenic temperatures such as those provided by superfluid helium as a solvent. Unfortunately, the cost of converged path integral simulations increases significantly upon lowering the temperature so that the computational burden of simulating matter at the typical superfluid helium temperatures becomes prohibitive. Here we investigate how accelerated path integral techniques based on colored noise generalized Langevin equations, in particular the socalled path integral generalized Langevin equation thermostat (PIGLET) variant, perform in this extreme quantum regime using as an example the quasirigid methane molecule and its highly fluxional protonated cousin, CH{sub 5}{sup +}. We show that the PIGLET technique gives a speedup of two orders of magnitude in the evaluation of structural observables and quantum kinetic energy at ultralow temperatures. Moreover, we computed the spatial spread of the quantum nuclei in CH{sub 4} to illustrate the limits of using such colored noise thermostats close to the many body quantummore »
 Authors:
 Lehrstuhl für Theoretische Chemie, Ruhr–Universität Bochum, 44780 Bochum (Germany)
 Laboratory of Computational Science and Modelling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne (Switzerland)
 Publication Date:
 OSTI Identifier:
 22679015
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Chemical Physics; Journal Volume: 145; Journal Issue: 5; Other Information: (c) 2016 Author(s); Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; EXPERIMENTAL DATA; GROUND STATES; KINETIC ENERGY; LANGEVIN EQUATION; MANYBODY PROBLEM; PATH INTEGRALS; QUANTUM MECHANICS; SIMULATION
Citation Formats
Uhl, Felix, Email: felix.uhl@rub.de, Marx, Dominik, and Ceriotti, Michele. Accelerated path integral methods for atomistic simulations at ultralow temperatures. United States: N. p., 2016.
Web. doi:10.1063/1.4959602.
Uhl, Felix, Email: felix.uhl@rub.de, Marx, Dominik, & Ceriotti, Michele. Accelerated path integral methods for atomistic simulations at ultralow temperatures. United States. doi:10.1063/1.4959602.
Uhl, Felix, Email: felix.uhl@rub.de, Marx, Dominik, and Ceriotti, Michele. Sun .
"Accelerated path integral methods for atomistic simulations at ultralow temperatures". United States.
doi:10.1063/1.4959602.
@article{osti_22679015,
title = {Accelerated path integral methods for atomistic simulations at ultralow temperatures},
author = {Uhl, Felix, Email: felix.uhl@rub.de and Marx, Dominik and Ceriotti, Michele},
abstractNote = {Path integral methods provide a rigorous and systematically convergent framework to include the quantum mechanical nature of atomic nuclei in the evaluation of the equilibrium properties of molecules, liquids, or solids at finite temperature. Such nuclear quantum effects are often significant for light nuclei already at room temperature, but become crucial at cryogenic temperatures such as those provided by superfluid helium as a solvent. Unfortunately, the cost of converged path integral simulations increases significantly upon lowering the temperature so that the computational burden of simulating matter at the typical superfluid helium temperatures becomes prohibitive. Here we investigate how accelerated path integral techniques based on colored noise generalized Langevin equations, in particular the socalled path integral generalized Langevin equation thermostat (PIGLET) variant, perform in this extreme quantum regime using as an example the quasirigid methane molecule and its highly fluxional protonated cousin, CH{sub 5}{sup +}. We show that the PIGLET technique gives a speedup of two orders of magnitude in the evaluation of structural observables and quantum kinetic energy at ultralow temperatures. Moreover, we computed the spatial spread of the quantum nuclei in CH{sub 4} to illustrate the limits of using such colored noise thermostats close to the many body quantum ground state.},
doi = {10.1063/1.4959602},
journal = {Journal of Chemical Physics},
number = 5,
volume = 145,
place = {United States},
year = {Sun Aug 07 00:00:00 EDT 2016},
month = {Sun Aug 07 00:00:00 EDT 2016}
}

We derive pathintegral expressions for the second and third virial coefficients of monatomic quantum gases. Unlike previous work that considered only Boltzmann statistics, we include exchange effects (BoseEinstein or FermiDirac statistics). We use stateoftheart pair and threebody potentials to calculate the third virial coefficient of {sup 3}He and {sup 4}He in the temperature range 2.624.5561 K. We obtain uncertainties smaller than those of the limited experimental data. Inclusion of exchange effects is necessary to obtain accurate results below about 7 K.

Lattice thermal conductivity of ultra high temperature ceramics ZrB{sub 2} and HfB{sub 2} from atomistic simulations
Atomistic GreenKubo simulations are performed to evaluate the lattice thermal conductivity for single crystals of the ultra high temperature ceramics ZrB{sub 2} and HfB{sub 2}. Recently developed interatomic potentials are used for these simulations. Heat current correlation functions show rapid oscillations, which can be identified with mixed metalBoron optical phonon modes. Results for temperatures from 300K to 1000K are presented. 
Comparisons Between Integral Equation Theory and Molecular Dynamics Simulations for Atomistic Models of Polyethylene Liquids
Molecular dynamics (MD) simulations were performed on dense liquids of polyethylene chains of 24 and 66 united atom CH{sub 2} units. A series of models was studied ranging in atomistic detail from coarsegrained, freelyjointed, tangent site chains to realistic, overlapping site models subjected to bond angle restrictions and torsional potentials. These same models were also treated with the selfconsistent, polymer reference interaction site model (PRISM) theory. The intramolecular and total structure factors, as well as, the intermolecular radial distribution functions g(r) and direct correlation functions C(r) were obtained from theory and simulation. Angular correlation functions were also simulation obtained frommore » 
''Direct'' calculation of quantum mechanical rate constants via path integral methods: Application to the reaction path Hamiltonian, with numerical test for the H+H/sub 2/ reaction in 3D
The method recently proposed by Miller, Schwartz, and Tromp for determining Boltzmann rate constants ''directly'': by the path integral evaluation of a reactive flux correlation function: is developed within the framework of the reaction path Hamiltonian model for a general polyatomic reaction. The expression for the correlation function, the time integral of which is the rate constant, is reduced to a single path integral over only one degree of freedom (the reaction coordinate). Effects of tunneling, ''frictional'' effects on the reaction coordinate due to coupling to other degrees of freedom, and the effects of recrossing the transition state dividing surfacemore »