Modeling Excited-State Proton Transfer to Solvent: A Dynamics Study of a Super Photoacid with a Hybrid Implicit/Explicit Solvent Model

The rapid growth of time-resolved spectroscopies and the theoretical advances in ab initio molecular dynamics (AIMD) pave the way to look at the real-time molecular motion following the electronic excitation. Here, we exploited the capabilities of AIMD combined with a hybrid implicit/explicit model of solvation to investigate the ultrafast excited-state proton transfer (ESPT) reaction of a super photoacid, known as QCy9, in water solution. QCy9 transfers a proton to a water solvent molecule within 100 fs upon the electronic excitation in aqueous solution, and it is the strongest photoacid reported in the literature so far. Because of the ultrafast kinetics, it has been experimentally hypothesized that the ESPT escapes the solvent dynamics control (Huppert et al., J. Photochem. Photobiol. A 2014, 277, 90). The sampling of the solvent configuration space on the ground electronic state is the first key step toward the simulation of the ESPT event. Therefore, several configurations in the Franck-Condon region, describing an average solvation, were chosen as starting points for the excited-state dynamics. In all cases, the excited-state evolution spontaneously leads to the proton transfer event, whose rate is strongly dependent on the hydrogen bond network around the proton acceptor solvent molecule. Our study revealed that the explicit representation at least of three solvation shells around the proton acceptor molecule is necessary to stabilize the excess proton. Furthermore, the analysis of the solvent molecule motions in proximity of the reaction site suggested that even in the case of the strongest photoacid, the ESPT is actually assisted by the solvation dynamics of the first and second solvation shells of the water accepting molecule.