Search papers, labs, and topics across Lattice.
The paper introduces Path-Integral Hybrid Monte Carlo with Enveloping Bridging Potentials (PIHMC-EBP), a novel method for computing numerically exact tunneling splittings in molecular systems by constructing barrierless bridging potentials that connect symmetry-related regions. Tailored nonlocal updates enhance sampling of slow collective motions, eliminating the need for quadrature and time-step convergence checks required by path-integral molecular dynamics. Applied to malonaldehyde, HCl dimer, and water dimer, PIHMC-EBP achieves the most precise tunneling splittings reported to date while significantly reducing computational cost.
Tunneling splittings, a notoriously difficult quantum property to compute, can now be calculated with orders-of-magnitude less compute thanks to a clever "bridging potential" trick.
A path-integral hybrid Monte Carlo approach with enveloping bridging potentials (PIHMC-EBP) is proposed for calculating numerically exact tunneling splittings in molecular systems. The central idea is to construct an approximately barrierless bridging potential that smoothly connects symmetry-related regions of ring-polymer phase space, enabling direct sampling of the free-energy profile from which the relevant splittings are obtained. Two tailored nonlocal updates are designed to enhance the sampling of slow collective motions. Compared with path-integral molecular dynamics using thermodynamic integration, PIHMC-EBP requires neither quadrature nor time-step convergence checks, thereby substantially reducing the manual effort required to analyze the results. Applications to malonaldehyde (and its deuterated isotopologue) and the HCl dimer using state-of-the-art potential energy surfaces provide the most precise tunneling splittings reported to date for both systems, while simultaneously reducing the overall computational cost by several times and three orders of magnitude, respectively. Finally, application to the water dimer yields the first numerically exact path-integral calculations of the ground-state tunneling splittings on three different potential energy surfaces, all obtained simultaneously by reweighting a single set of trajectories.
