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

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.

References in corpus (18)

• Towards Exact Molecular Dynamics Simulations with Machine-Learned Force Fields
• An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics
• Efficient stochastic thermostatting of path integral molecular dynamics
• Nuclear quantum effects in solids using a colored-noise thermostat
• Colored-noise thermostats à la carte
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• High Order Path Integrals Made Easy
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• Portable implementation of a quantum thermal bath for molecular dynamics simulations
• Zero-Point Energy Leakage in Quantum Thermal Bath Molecular Dynamics Simulations
• Nuclear quantum effects in ab initio dynamics: theory and experiments for lithium imide
• Quantum Ring-Polymer Contraction Method: Including nuclear quantum effects at no additional computational cost in comparison to ab-initio molecular dynamics
• Quantum effects in graphene monolayers: Path-integral simulations
• Simulation of quantum zero-point effects in water using a frequency-dependent thermostat
• High-dimensional neural network potentials for solvation: The case of protonated water clusters in helium

Cited by in corpus (11)

• Machine Learning Force Fields
• Committee neural network potentials control generalization errors and enable active learning
• Automated Fitting of Neural Network Potentials at Coupled Cluster Accuracy: Protonated Water Clusters as Testing Ground
• Dynamical Strengthening of Covalent and Non-Covalent Molecular Interactions by Nuclear Quantum Effects at Finite Temperature
• Reduced Rovibrational Coupling Cartesian Dynamics for Semiclassical Calculations: Application to the Spectrum of the Zundel Cation
• Transferability of machine learning potentials: Protonated water neural network potential applied to the protonated water hexamer
• Quantum Nature of the Hydrogen Bond from Ambient Conditions down to Ultra-low Temperatures
• Random sampling versus active learning algorithms for machine learning potentials of quantum liquid water
• Deciphering High-order Structural Correlations within Fluxional Molecules from Classical and Quantum Configurational Entropy
• Machine-learning interatomic potentials achieving CCSD(T) accuracy for systems with extended covalent networks and van der Waals interactions
• Free energy profiles for chemical reactions in solution from high-dimensional neural network potentials: The case of the Strecker synthesis