A Path Integral Molecular Dynamics Simulation of a Harpoon-Type Redox Reaction in a Helium Nanodroplet

We present path integral molecular dynamics (PIMD) calculations of an electron transfer from a heliophobic Cs2 dimer in its (3 sigma u) state, located on the surface of a He droplet, to a heliophilic, fully immersed C60 molecule. Supported by electron ionization mass spectroscopy measurements (Renzler et al., J. Chem. Phys. 2016, 145, 181101), this spatially quenched reaction was characterized as a harpoon-type or long-range electron transfer in a previous high-level ab initio study (de Lara-Castells et al., J. Phys. Chem. Lett. 2017, 8, 4284). To go beyond the static approach, classical and quantum PIMD simulations are performed at 2 K, slightly below the critical temperature for helium superfluidity (2.172 K). Calculations are executed in the NVT ensemble as well as the NVE ensemble to provide insights into real-time dynamics. A droplet size of 2090 atoms is assumed to study the impact of spatial hindrance on reactivity. By changing the number of beads in the PIMD simulations, the impact of quantization can be studied in greater detail and without an implicit assumption of superfluidity. We find that the reaction probability increases with higher levels of quantization. Our findings confirm earlier, static predictions of a rotational motion of the Cs2 dimer upon reacting with the fullerene, involving a substantial displacement of helium. However, it also raises the new question of whether the interacting species are driven out-of-equilibrium after impurity uptake, since reactivity is strongly quenched if a full thermal equilibration is assumed. More generally, our work points towards a novel mechanism for long-range electron transfer through an interplay between nuclear quantum delocalization within the confining medium and delocalized electronic dispersion forces acting on the two reactants.