Abstract:For a long time, performing converged path integral simulations at ultra-low, but finite temperatures of a few Kelvin has been a nearly impossible task. However, recent developments in advanced colored noise thermostatting schemes for path integral simulations, namely the Path Integral Generalized Langevin Equation Thermostat (PIGLET) and the Path Integral Quantum Thermal Bath (PIQTB), have been able to greatly reduce the computational cost of these simulations, thus making the ultra-low temperature regime accessible in practice. In this work, we investigate the influence of these two thermostatting schemes on the description of hydrogen-bonded systems at temperatures down to a few Kelvin as encountered, for example, in helium nanodroplet isolation or tagging photodissociation spectroscopy experiments. For this purpose, we analyze the prototypical hydrogen bond in the Zundel cation (H$_5$O$_{2}^{+}$) as a function of both, oxygen-oxygen distance and temperature in order to elucidate how the anisotropic quantum delocalization and, thus, the shape of the shared proton adapts depending on the donor-acceptor distance. The underlying electronic structure of the Zundel cation is described in terms of Behler's Neural Network Potentials of essentially converged Coupled Cluster accuracy, CCSD(T*)-F12a/AVTZ. In addition, the performances of the PIQTB and PIGLET methods for energetic, structural, and quantum delocalization properties are assessed and directly compared. Overall, our results emphasize the validity and practical usefulness of these two modern thermostatting approaches for path integral simulations of hydrogen-bonded systems even at ultra-low temperatures.

Path Integral Calculations of the Hydrogen Hugoniot Using Augmented Nodes

An ab initio approach to the Hugoniot

Study of Equation of State of the Dense Fluid He+H_2 Mixtures

A Method for Calculating the Heats of Formation of Medium-Sized and Large-Sized Molecules

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

Path Integral Simulation of Hydrogen Adsorption in Single-Walled Carbon Nanotubes at Low Temperatures

Stable liquid Hydrogen at high pressure by a novel ab-initio molecular dynamics

Mie–FH: A quantum corrected pair potential in the LAMMPS simulation package for hydrogen mixtures

Approaching exact hyperpolarizabilities via sum-over-states Monte Carlo configuration interaction

Ab initio path integral Monte Carlo simulations of warm dense two-component systems without fixed nodes: Structural properties

Direct free energy calculation from ab initio path integral Monte Carlo simulations of warm dense matter

Converged Colored Noise Path Integral Molecular Dynamics Study of the Zundel Cation down to Ultra-low Temperatures at Coupled Cluster Accuracy

On the Room-Temperature Phase Diagram of High Pressure Hydrogen: an Ab Initio Molecular Dynamics Perspective and A Diffusion Monte Carlo Study

Density Functional Study of Hydrogen Adsorption at Low Temperatures

A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory

Stochastic path-integral approach for predicting the superconducting temperatures of anharmonic solids

Towards the same line of liquid-liquid phase transition of dense hydrogen from various theoretical predictions

High-precision solution of the Dirac Equation for the hydrogen molecular ion using a basis-set expansion

Analyses of hydrogen local environments in metals and intermetallic compounds using inelastic neutron scattering calculations based on first-principles hydrogen adiabatic potentials

Ab initio Density Response and Local Field Factor of Warm Dense Hydrogen