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.