Non-adiabatic Origin of Roaming OH Dynamics in the Formic Acid Dimer Dication

Abstract

Ionization of molecular clusters can trigger chemical reactions and drive chemical evolution even at very low temperatures, influencing chemistry in interstellar, atmospheric, and planetary environments exposed to ionizing radiation. To investigate such processes involving the dissociation of both intramolecular and intermolecular bonds under controlled conditions, we examined the dynamics of the formic acid (FA) dimer dication in an ultrafast extreme-ultraviolet (EUV) pump and near-infrared (NIR) probe experiment, combined with ab initio molecular dynamics simulations. The dissociation of the intermolecular bond and formation of the two-body FA⁺ + FA⁺ channel could be explained by ground-state dynamics, whereas the three-body breakup channels required a more detailed description.We developed a simplified dimer model for the breakup process that enabled carrying out non-adiabatic molecular dynamics simulations on excited-state CASPT2 potential energy surfaces that capture intermolecular as well as intramolecular dynamics. The simulations showed that soon after the dimer dissociation and non-adiabatic decay to the ground electronic state, a roaming-OH mechanism develops, accounting for the observed kinetic-energy-release distributions and momentum correlations in the FA⁺ + CHO⁺ + OH, and FA⁺ + H₂O⁺ + CO three-body breakup channels. The simplified modeling approach may serve as a practical framework for studying excited-state dynamics in molecular dimer and cluster breakup processes.