Besides the direct description of the temporal evolution of the solvent shell around [Fe(CN) 4 (bipy)] 2--as discernable from the 3D-SDFs-our simulations can shed light on a number of additional questions, as posed in the introduction. We begin with "How inhomogeneous is the solvent response around the molecule?", where, in particular, we refer to the time constants of solvent response around the different moieties. An analysis of the time-dependent hydrogen bond counts (Section 2.5 and Figure ) reveals that the breaking of the hydrogen bonds to the cyanide ligands occurs biexponentially with time constants of about 15 and 240 fs. The new hydrogen bonds to the bipyridyl are formed more slowly, with time constants of 30 and 500 fs. Reasons for this behaviour are the more localized change in charge on the cyanides (compared to the delocalized charge on the bipyridyl), as well as that the waters forming the bipyridyl hydrogen bonds must first be released from the cyanides, as we will discuss in more detail below. The different time scales of breaking and forming of hydrogen bonds can also be extracted from the time-dependent RDFs (Section S2.3, Figures and). Quantifying the solvent response near the cyanide and bipyridine ligands provides time constants of 50-75 fs and 280 fs, respectively. Moreover, the obtained high-fidelity RDFs afford to identify the differences in the solvent dynamics around the equatorial and axial cyanide ligands (Fig. ), an endeavor that was not feasible before. We find that the axial cyanides show a stronger response in their first solvation shell compared to the equatorial ones. Interestingly, both kinds of cyanide solvent shells show some oscillatory response in the RDFs, but with different oscillation periods of 80 fs (axial) and 160 fs (equatorial). All these results support the finding that the solvent dynamics is very inhomogeneous around the molecule.