Exact tunneling splittings from symmetrized path integrals
arXiv:2305.11002 · doi:10.1063/5.0158879
Abstract
We develop a new simulation technique based on path-integral molecular dynamics for calculating ground-state tunneling splitting patterns from ratios of symmetrized partition functions. In particular, molecular systems are rigorously projected onto their J = 0 rotational state by an "Eckart spring" that connects two adjacent beads in a ring polymer. Using this procedure, the tunneling splitting can be obtained from thermodynamic integration at just one (sufficiently low) temperature. Converged results are formally identical to the values that would have been obtained by solving the full rovibrational Schrödinger equation on a given Born--Oppenheimer potential energy surface. The new approach is showcased with simulations of hydronium and methanol, which are in good agreement with wavefunction-based calculations and experimental measurements. The method will be of particular use for the study of low-barrier methyl rotations and other floppy modes, where instanton theory is not valid.
arXiv comment15 pages, 4 figures
References in corpus (12)
- Competing quantum effects in the dynamics of a flexible water model
- Accelerating the convergence of path integral dynamics with a generalized Langevin equation
- Elucidating the NuclearQuantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene
- Path Integral Molecular Dynamics for Bosons
- Full- and reduced-dimensionality instanton calculations of the tunnelling splitting in the formic acid dimer
- High Order Path Integrals Made Easy
- Path-integral dynamics of water using curvilinear centroids
- Perturbatively corrected ring-polymer instanton theory for accurate tunneling splittings
- Instanton theory of tunnelling in molecules with asymmetric isotopic substitutions
- ExoMol line lists -- XL. Ro-vibrational molecular line list for the hydronium ion (HO)
- Quantum symmetry from enhanced sampling methods
- Tunneling splittings of vibrationally excited states using general instanton paths