Crucial to many light-driven processes in transition metal complexes is the absorption and dissipation of energy by 3d electrons1–4. But a detailed understanding of such non-equilibrium excited-state dynamics and their interplay with structural changes is challenging: a multitude of excited states and possible transitions result in phenomena too complex to unravel when faced with the indirect sensitivity of optical spectroscopy to spin dynamics5 and the flux limitations of ultrafast X-ray sources6,7. Such a situation exists for archetypal polypyridyl iron complexes, such as [Fe(2,2′-bipyridine)3]2+, where the excited-state charge and spin dynamics involved in the transition from a low- to a high-spin state (spin crossover) have long been a source of interest and controversy6–15. Here we demonstrate that femtosecond resolution X-ray fluorescence spectroscopy, with its sensitivity to spin state, can elucidate the spin crossover dynamics of [Fe(2,2′-bipyridine)3]2+ on photoinduced metal-to-ligand charge transfer excitation. We are able to track the charge and spin dynamics, and establish the critical role of intermediate spin states in the crossover mechanism. We anticipate that these capabilities will make our method a valuable tool for mapping in unprecedented detail the fundamental electronic excited-state dynamics that underpin many useful light-triggered molecular phenomena involving 3d transition metal complexes.

We develop a general procedure to analyse the pre-edges in 1s x-ray absorption near edge structure (XANES) of transition metal oxides and coordination complexes. Transition metal coordination complexes can be described from a local model with one metal ion. The 1s 3d quadrupole transitions are calculated with the charge-transfer multiplet program. Tetrahedral coordination complexes have more intense pre-edge structures due to the local mixing of 3d and 4p states, implying a combination of 1s 3d quadrupole and 1s 4p dipole transitions. Divalent transition metal oxides can be described similar to coordination complexes, but for trivalent and tetravalent oxides, additional structures are visible in the pre-edge region due to non-local dipole transitions. The 1s 4p dipole transitions have large cross section at the 3d-band region due to the strong metal-metal interactions, which are oxygen mediated. This yields large intensity in the 3d-band region but at a different energy than the local 1s 3d quadrupole transitions because of smaller core-hole effects due to the delocalization of the excited electron.

The performance of six frequently used density functional theory (DFT) methods (RPBE, OLYP, TPSS, B3LYP, B3LYP*, and TPSSh) in the prediction of Mössbauer isomer shifts(δ) and quadrupole splittings (ΔEQ) is studied for an extended and diverse set of Fe complexes. In addition to the influence of the applied density functional and the type of the basis set, the effect of the environment of the molecule, approximated with the conducting-like screening solvation model (COSMO) on the computed Mössbauer parameters, is also investigated. For the isomer shifts the COSMO-B3LYP method is found to provide accurate δ values for all 66 investigated complexes, with a mean absolute error (MAE) of 0.05 mm s–1 and a maximum deviation of 0.12 mm s–1. Obtaining accurate ΔEQ values presents a bigger challenge; however, with the selection of an appropriate DFT method, a reasonable agreement can be achieved between experiment and theory. Identifying the various chemical classes of compounds that need different treatment allowed us to construct a recipe for ΔEQ calculations; the application of this approach yields a MAE of 0.12 mm s–1 (7% error) and a maximum deviation of 0.55 mm s–1 (17% error). This accuracy should be sufficient for most chemical problems that concern Fe complexes. Furthermore, the reliability of the DFT approach is verified by extending the investigation to chemically relevant case studies which include geometric isomerism, phase transitions induced by variations of the electronic structure (e.g., spin crossover and inversion of the orbital ground state), and the description of electronically degenerate triplet and quintet states. Finally, the immense and often unexploited potential of utilizing the sign of the ΔEQ in characterizing distortions or in identifying the appropriate electronic state at the assignment of the spectral lines is also shown.