Investigating Tunneling-Controlled Chemical Reactions through Ab Initio Ring Polymer Molecular Dynamics

Quantum rates in dissipative systems with spatially varying friction

We investigate whether making the friction spatially dependent on the reaction coordinate introduces quantum effects into the thermal reaction rates for dissipative reactions. Quantum rates are calculated using the numerically exact multi-configuration time-dependent Hartree method, as well as the approximate ring-polymer molecular dynamics (RPMD), ring-polymer instanton methods, and classical molecular dynamics. By conducting simulations across a wide range of temperatures and friction strengths, we can identify the various regimes that govern the reactive dynamics. At high temperatures, in addition to the spatial-diffusion and energy-diffusion regimes predicted by Kramer's rate theory, a (coherent) tunneling-dominated regime is identified at low friction. At low temperatures, incoherent tunneling dominates most of Kramer's curve, except at very low friction, when coherent tunneling becomes dominant. Unlike in classical mechanics, the bath's influence changes the equilibrium time-independent properties of the system, leading to a complex interplay between spatially dependent friction and nuclear quantum effects even at high temperatures. More specifically, a realistic friction profile can lead to an increase (or decrease) of the quantum (classical) rates with friction within the spatial-diffusion regime, showing that classical and quantum rates display qualitatively different behaviors. Except at very low frictions, we find that RPMD captures most of the quantum effects in the thermal reaction rates.

Beyond Cavity Born-Oppenheimer: On Nonadiabatic Coupling and Effective Ground State Hamiltonians in Vibro-Polaritonic Chemistry

The emerging field of vibro-polaritonic chemistry studies the impact of light-matter hybrid states known as vibrational polaritons on chemical reactivity and molecular properties. Here, we discuss vibro-polaritonic chemistry from a quantum chemical perspective beyond the cavity Born-Oppenheimer (CBO) approximation and examine the role of electron-photon correlation in effective ground state Hamiltonians. We first quantitatively review ab initio vibro-polaritonic chemistry based on the molecular Pauli-Fierz Hamiltonian in dipole approximation and a vibrational strong coupling (VSC) Born-Huang expansion. We then derive nonadiabatic coupling elements arising from both "slow" nuclei and cavity modes compared to "fast" electrons via the generalized Hellmann-Feynman theorem, discuss their properties, and reevaluate the CBO approximation. In the second part, we introduce a crude VSC Born-Huang expansion based on adiabatic electronic states, which provides a foundation for widely employed effective Pauli-Fierz Hamiltonians in ground state vibro-polaritonic chemistry. Those do not strictly respect the CBO approximation but an alternative scheme, which we name crude CBO approximation. We argue that the crude CBO ground state misses electron-photon correlation relative to the CBO ground state due to neglected cavity-induced nonadiabatic transition dipole couplings to excited states. A perturbative connection between both ground state approximations is proposed, which identifies the crude CBO ground state as a first-order approximation to its CBO counterpart. We provide an illustrative numerical analysis of the cavity Shin-Metiu model with a focus on nonadiabatic coupling under VSC and electron-photon correlation effects on classical activation barriers. We finally discuss the potential shortcomings of the electron-polariton Hamiltonian when employed in the VSC regime.

A highly accurate full-dimensional ab initio potential surface for the rearrangement of methylhydroxycarbene (H3C-C-OH)

We report here a full-dimensional machine learning global potential surface (PES) for the rearrangement of methylhydroxycarbene (H₃C-C-OH, 1t). The PES is trained with the fundamental invariant neural network (FI-NN) method on 91 564 ab initio energies calculated at the UCCSD(T)-F12a/cc-pVTZ level of theory, covering three possible product channels. FI-NN PES has the correct symmetry properties with respect to permutation of four identical hydrogen atoms and is suitable for dynamics studies of the 1t rearrangement. The averaged root mean square error (RMSE) is 11.4 meV. Six important reaction pathways, as well as the energies and vibrational frequencies at the stationary geometries on these pathways are accurately preproduced by our FI-NN PES. To demonstrate the capacity of the PES, we calculated the rate coefficient of hydrogen migration in -CH₃ (path A) and hydrogen migration of -OH (path B) with instanton theory on this PES. Our calculations predicted the half-life of 1t to be 95 min, which is excellent in agreement with experimental observations.