Abstract
The response of a solvation shell to molecular solute photoexcitation is an ubiquitous phenomenon of great relevance in chemistry. This response can occur within just few tens of femtoseconds, making it very challenging to resolve experimentally. Thus, the details of the (an)isotropy of the solvent response around a solute, the presence of coherent solvent fluctuations, hydrogen bond reorganization mechanisms, and the intricate interplay between electronic, spin, nuclear, and solvent dynamics remain elusive. Here, we report large-scale nonadiabatic molecular dynamics simulations of [Fe(CN)(4)(bipy)](2-) (bipy=2,2'-bipyridine) in water, where the electronic evolution from singlet metal-to-ligand charge transfer (MLCT) states to triplet MLCT and metal-centered (MC) states overlaps temporally with the molecule's nuclear motion and a strong solvent shell response. We leverage vibronic coupling model potentials combined with electrostatic embedding, within our so-called vibronic coupling/molecular mechanics (VC/MM) method, to be able to compute several thousand nonadiabatic excited-state trajectories, including all relevant singlet and triplet states as well as several thousand explicit water molecules. This superior statistics affords an unprecedented view on the three-dimensional solvent distribution dynamics at few-fs and sub-Å resolution. The results reveal a direct solvent migration mechanism, where excitation to the MLCT states leads to the breaking of hydrogen bonds to the cyanide ligands within less than 100 fs, followed by the formation of hydrogen bonds with the negatively charged bipyridyl ligand by the same water molecules. Furthermore, the MLCT and MC states show very distinct solvent responses, which are overlapping in time, as governed by the electronic dynamics. More broadly, this work demonstrates how VC/MM nonadiabatic dynamics simulations can resolve anisotropic solvent dynamics around a photoexcited solute with unprecedented detail, offering a new perspective that could stimulate the development of time-resolved experimental techniques capable of probing such solvent behaviour.