The femtosecond duration of the intense hard X-ray pulses generated by the LCLS (Linac Coherent Light Source) X-ray free-electron laser^{16, 17} creates the opportunity to study spin dynamics with iron 3p–1s (Kβ) X-ray fluorescence spectroscopy^{18, 19}. Figure 1 shows diagrams of the measurement technique and relevant energy levels (Fig. 1–c), a ‘ball-and-stick’ representation of the [Fe(2,2′-bipyridine)₃]²⁺ complex (Fig. 1d), and the dependence of photoexcited spin crossover dynamics on the Fe–ligand distance (Fig. 1e). Given the roughly 100 femtosecond (fs) time resolution of the measurement¹⁷, the subfemtosecond lifetime of the iron 1s core hole makes X-ray fluorescence an effectively instantaneous probe²⁰. A variety of distinct electronic excited states, including singlet and triplet metal-to-ligand charge transfer states (^1,3MLCT), triplet ligand field excited states (³T) and quintet ligand field excited states (⁵T₂) have been proposed to participate in the spin crossover mechanism^{6, 8, 10, 11, 21, 22} (Fig. 1e). Distinguishing electronic excited states with different charge and spin density, such as the ^1,3MLCT, ³T and ⁵T₂ states listed above, represents a critical step in characterizing the spin crossover mechanism.

Figure 2a shows the sensitivity of the iron Kβ fluorescence spectrum to the 3d spin moment, a sensitivity that results from the exchange interaction between the 3p and 3d electrons^{18, 19, 23, 24, 25}. Equally important, the ground-state spectra of iron coordination complexes with different ligation, but the same iron spin moment, exhibit similar Kβ fluorescence spectra. This insensitivity of Kβ fluorescence spectroscopy to the details of the coordinating ligands and the local symmetry of the complex has previously been used to characterize the electronic ground-state spin moment of a variety of molecular systems^{19, 25}. We note that the insensitivity of the Kβ fluorescence spectrum to the electronic properties of the ligand means that the spectrum cannot be used to distinguish between singlet and triplet MLCT states. We utilize these spectra of distinct spin configurations to model transient difference spectra—that is, the time and energy dependence of the fluorescent amplitude difference between excited-state and ground-state spectra. Figure 2b shows the model complex difference spectra generated from the ground-state spectra of the relevant excited-state spin configurations and the singlet ground state. These model complex difference spectra confirm that each excited-state spin moment generates a distinct difference spectrum that cannot be reproduced by a linear combination of the other difference spectra (see Fig. 2, Extended Data Fig. 1 and Methods for details).

The time-resolved Kβ fluorescence spectra provide the sensitivity to spin dynamics needed to answer a critical question regarding the spin crossover mechanism: does the ⁵T₂ state form directly from the ^1,3MLCT state^{6, 13, 26}, or does spin crossover involve a ³T transient^{8, 10}? Ultraviolet–visible transient absorption^{13, 14}, time-resolved luminescence¹³, and time-resolved iron K-edge XANES⁶ have been used to characterize the spin crossover dynamics of [Fe(2,2′-bipyridine)₃]²⁺, and the similar rates measured for ³MLCT decay and ⁵T₂ formation were attributed to the ³MLCT excited state converting directly to the ⁵T₂ excited state, although a conversion including transient triplet states was also considered⁶. Potential energy surfaces calculated for this system allow either mechanism to proceed with minimal reaction barriers^{21, 22}, but cannot explain why the ^1,3MLCT and ⁵T₂ states should be strongly coupled: the leading order spin–orbit interaction cannot couple the ^1,3MLCT and ⁵T₂ states because a transition between these states requires the excitation of two electrons on two distinct centres, whereas spin–orbit coupling is predominantly a single-centre, one-electron operator²².

Figure 2c, d shows the transient difference spectra for [Fe(2,2′-bipyridine)₃]²⁺ measured for a 50-fs and a 1-ps (picosecond) time delay. The spectrum in Fig. 2d clearly demonstrates the ease of identifying the ⁵T₂ state with the Kβ fluorescence spectrum. Determining whether spin crossover from the ^1,3MLCT to the ⁵T₂ proceeds through a transient ³T state proves more challenging because the relaxation dynamics do not lead to a time regime where the majority of the excited molecules reside in the ³T excited state. The significant difference between the spectra in Fig. 2c and d, however, clearly demonstrates the presence of excited-state species other than the ⁵T₂ state. With statistically rigorous kinetic modelling, ^1,3MLCT, ³T and ⁵T₂ states can be clearly distinguished in the relaxation dynamics probed with Kβ fluorescence.

The ability to spectroscopically distinguish between ^1,3MLCT, ³T and ⁵T₂ electronic excited states allows the spin crossover mechanism to be determined from the time evolution of the iron Kβ fluorescence spectrum. The time-resolved difference spectra, model fits of the difference spectra, and the parameters extracted from the fit can be found in Fig. 3, Extended Data Figs 2, 3, 4 and Extended Data Table 1. We have fitted the difference spectra to two distinct models: one where the ^1,3MLCT decays directly to a ⁵T₂ excited state and one where the ^1,3MLCT relaxes to the ⁵T₂ state via a ³T transient. Figure 3b, c shows the time-dependent difference signal measured at two X-ray fluorescence energies: 7,061 eV, where the difference signal is largest, and 7,054 eV, where the triplet model complex has a spectral signature clearly distinct from the ^1,3MLCT and ⁵T₂ states as shown in Fig. 2b. The fits in Fig. 3b, c have been determined from a global analysis of the full time-dependent spectra. The statistical significance of the more complex kinetic model involving the triplet transient can be determined from an F-test comparison of the two models (described in Methods). The reduction in residuals achieved with the model containing the triplet transient is sufficient to reject the direct ^1,3MLCT⁵T₂ model with greater than 95% confidence. Note that the successful use of a kinetic model to describe subpicosecond dynamics implies that the Kβ spectra do not depend significantly on the time-evolving nuclear structure, consistent with the insensitivity of the ground-state Kβ spectra to the ligand details.

The successful analysis of the experimental data relies on two constraints presented by the model spectra shown in Fig. 2b and two constraints derived from the kinetic models. We force (1) the shape and (2) the relative amplitudes of the difference signals for the ^1,3MLCT, ³T and ⁵T₂ electronic excited states to match the shape and relative amplitudes of the model complex difference spectra. We also require (3) all X-ray fluorescence energies to be fitted with a single time zero and (4) all MLCT excited states to undergo spin crossover, consistent with previous measurements of the spin crossover quantum yield¹³. The ultrafast rise of the difference signal shown in Fig. 3b greatly constrains the value of time zero and the final ⁵T₂ state population. For the fit to the direct spin crossover mechanism shown in Fig. 3b, the fast rise in signal at 7,061 eV requires a fast rise in ⁵T₂ population. As shown in Fig. 3c, the fast rise in the direct mechanism fit at 7,061 eV also leads to a fast drop in signal at 7,054 eV, because the ⁵T₂ state has a negative difference signal at 7,054 eV. For the fit to the sequential spin crossover mechanism also shown in Fig. 3b, the fast rise in signal at 7,061 eV can be accommodated initially by a rise in ³T population. Because the ³T state does not have a negative difference signal at 7,054 eV, the fast rise in ³T population does not lead to a fast drop at 7,054 eV. The stepwise transition through the ³T leads to a delayed onset of the drop in fluorescence amplitude at 7,054 eV relative to the rise in signal at 7,061 eV, consistent with the experimental data. For the direct model, a shift in time zero to fit the data in Fig. 3c would lead to a poor fit of the data in Fig. 3b.

Relaxation to the ⁵T₂ excited state via a ³T transient provides a more satisfying explanation for the relaxation dynamics. We speculate that the sequential relaxation occurs more promptly than the direct crossover from the ^1,3MCLT to the ⁵T₂ excited state because the sequential transition involves single electronic transitions coupled by a spin–orbit operator, whereas the direct transition involves the simultaneous transition of two distinct electrons on two centres and cannot occur with the first-order spin–orbit operator. The sequential relaxation, like the direct transition, provides an energetically feasible pathway with minimal reaction barriers between states that can be coupled with standard spin–orbit interactions²². The spin–orbit matrix elements in conjunction with the calculated potential energies of a variety of electronic excited states of [Fe(2,2′-bipyridine)₃]²⁺ as a function of the metal–ligand bond distance provide an approximate explanation for the fast intersystem crossing and the extremely short lifetime of the ³T excited state. A diagram of these potential energy surfaces can be found in Fig. 1e. In principle, the triplet ligand field excited state could be either a ³T₁ or a ³T₂ state. Computations indicate a crossing of the ³T₂ state in the Franck–Condon region of the ^1,3MLCT excited state and that the ^1,3MLCT³T₂⁵T₂ pathway dominates²⁷; however, relaxation trajectories involving the ³T₁ ligand field excited state remain plausible, and more definite conclusions will require a more complete calculation of the multidimensional potential energy surfaces, including the potentially important role of metal–ligand torsional motion²⁸. The sequential model fit in Fig. 3 gives a 150 ± 50 fs time constant for ^1,3MLCT decay to the ³T state and a 70 ± 30 fs time constant for ³T decay to the ⁵T₂ state. Although the mechanistic conclusions we have drawn from our measurements differ from the earlier interpretation²⁶, our experimental findings do not contradict the earlier results, but rather expand on them. The extracted decay time for the ^1,3MLCT excited state and the effective rise time for the ⁵T₂ excited state agree with the time constants observed previously within experimental error²⁶. The similarity of the ^1,3MLCT decay time and the ⁵T₂ rise time results from the rate of ³T decay being greater than that of ³T formation. This inhibits the build-up of molecules in the ³T excited state and challenges the temporal differentiation of the distinct electronic states involved in spin crossover (see Extended Data Fig. 2d). Only with a technique highly sensitive to the iron spin multiplicity can the presence of the ³T transient excited state in the relaxation dynamics be robustly resolved.

The complex excited-state electronic structure of molecules containing transition metals has inhibited the unambiguous interpretation of experimental measurements and the development of excited-state quantum dynamics simulations. We have demonstrated here that ultrafast X-ray fluorescence spectroscopy enables robust measurements of the charge and spin dynamics integral to excited-state relaxation in 3d transition-metal coordination complexes, which represents an important step towards an incisive mechanistic understanding of excited-state dynamics in 3d transition metal complexes.

Experimental procedures

We performed femtosecond hard X-ray fluorescence measurements on a 50 mM solution of [Fe(2,2′-bipyridine)₃]²⁺ in water at the X-ray pump-probe (XPP) instrument at the Linac Coherent Light Source (LCLS). The experiment used a 0.1 mm thick planar liquid jet oriented at an angle of 45° with respect to the direction of the incident X-ray beam. We measured the ultraviolet–visible absorption spectrum of the solution before and after the measurement to ensure no appreciable sample damage had occurred. The sample solution was collinearly excited with a 70 fs FWHM 520 nm laser beam (120 mJ cm⁻²) generated by optical parametric amplification of the 800 nm output of a Ti:sapphire regenerative amplifier laser system (Coherent, Legend). With 520 nm light, we excited [Fe(2,2′-bipyridine)₃]²⁺ at the peak of the MLCT band (Extended Data Fig. 5a). We set the pump laser fluence to maximize excitation yield, while avoiding other deleterious photophysical phenomena. Previous time-resolved hard X-ray spectroscopy measurements of iron spin crossover compounds have used higher, often significantly higher, optical laser fluence^{29, 30, 31}. We used an excitation laser fluence where the transient optical signal changes linearly with pump fluence, as shown in Extended Data Fig. 5b. The 8 keV X-ray laser pulses, with an average bandwidth of 0.3%, were focused using Be compound refractive lenses to a 50 μm diameter spot size at the sample position. Shot-to-shot fluctuations in the X-ray incidence energy and band width do not influence the X-ray fluorescence spectrum when the X-ray energy is well above the core ionization threshold. For iron, with a 1s ionization threshold of 7.112 keV, the 8 keV X-ray energy used in the experiment achieves this goal.

The incoming X-ray pulse energy was measured using non-invasive diagnostics before the sample³². A high-resolution energy dispersive X-ray emission spectrometer³³, based on the von Hamos geometry, was used to capture the iron 3p–1s (Kβ) fluorescence. The spectrometer was equipped with 4 cylindrically bent (0.5 m radius) Ge(620) crystal analysers and set to cover a Bragg angle range from 78.0° to 80.4°. The CSPAD 2D pixel array detector (388 × 370 pixels)³⁴ intersected the X-rays diffracted from the crystal analysers in an energy range from 7,033 to 7,084 eV.

The detector response calibration involved a pixel dependent dark current (pedestal) subtraction, a common mode offset, and an experimentally determined gain map. The gain map was built from histograms of each pixel response extracted from multiple images (after dark current and common mode offset corrections) collected over many minutes. Gaussians were fitted to the zero and one photon peaks of the histograms, enabling fine-tuned dark and gain corrections to the histograms directly from the data. The zero photon peaks were centred at zero analogue-to-digital units and the separation between the zero and one photon peaks were scaled to unity for all pixels. The counts for each pixel in a given time-step were obtained by averaging the analogue-to-digital values above a threshold of 2.5σ of the zero-photon peak and scaling to the incident X-ray intensity. The final 1D spectrum for each time-step was obtained by integrating the signal in the non-dispersive direction³³.

The shot-to-shot X-ray–optical relative time of arrival fluctuations were measured for every X-ray–optical pulse pair with a timing diagnostic tool based on optical detection of X-ray generated carriers in a Si₃N₄ thin film. A description of the time diagnostic tool and the demonstrated performance of the tool can be found elsewhere^{17, 35}. This experimental measure of the relative timing can be used to sort each experimental shot by the relative time of arrival. Although the timing tool provides an accurate measure of the shot-to-shot variation in the relative time of arrival between the X-ray and optical laser pulses, it does not provide an accurate measure of the instrument response function. The timing tool uses changes in the Si₃N₄ dielectric function to modify the transmission of a chirped white light pulse through the Si₃N₄ thin film. These changes in the dielectric function result from the increase in free carriers generated by X-ray ionization, Auger relaxation and impact ionization. The temporal response is the convolution of these complex dynamics with the cross-correlation of the X-ray and optical laser pulses. Without a detailed model of the carrier generation, the cross-correlation cannot be extracted from the timing tool. At present, no experimental means of cross-correlating the hard X-ray and optical pulses has been demonstrated.

The final time resolution of the experiment results from the convolution of the optical and X-ray pulse durations, the group velocity walk-off of the X-ray and optical pulses in the sample and the error in the relative time of arrival measurement. These factors would predict a cross-correlation of roughly 150 fs FWHM. In the data analysis, the instrument response function FWHM and time zero (coincident arrival of the X-ray and optical pulses) are fit parameters.

Kβ fluorescence spectra for model complexes

The 3p–1s X-ray (Kβ) fluorescence spectra of model complexes play an important role in our analysis of the time-dependent data. The Kβ fluorescence spectra of 3d transition-metal ions reflect the 3p−3d exchange interaction, which makes the line shapes sensitive to the spin state of the transition metal atom^{19, 23, 24, 36, 37}. Kβ fluorescence provides a powerful technique for spin state studies, particularly when there are advantages of working with penetrating hard X-rays. When a sample contains multiple spin states, the spin state distribution can be readily and precisely calculated from the line shape variations¹⁹.

We measured the Kβ fluorescence spectra of a series of iron complexes with different spin states at beamline 6-2 of the Stanford Synchrotron Radiation Lightsource (SSRL). All the samples were cooled to 10 K to reduce the influence of X-ray damage. The static spectra, collected with a multi-crystal high-resolution X-ray emission spectrometer, are shown in Fig. 2a.

We use the model complex difference spectra generated from molecules that have different spin multiplicities in their electronic ground state to model the time-dependent populations of electron excited states with different spin multiplicities. We verify the validity of using the model complex difference spectra generated from the quintet [Fe(phenanthroline)₂(NCS)₂] and the singlet [Fe(2,2′-bipyridine)₃]²⁺ model compounds for the quintet excited state by comparing it with the transient difference spectra of [Fe(2,2′-bipyridine)₃]²⁺ after a 1-ps time delay (see Fig. 2d). The validity of model complex difference spectra for the ^1,3MLCT and ³T excited states proves more challenging to demonstrate because we do not isolate these excited states at any time delay in our pump probe measurements (the fit to the 50-fs time delay spectra shown in Fig. 2c indicates a population ratio of 5:1.3:1 for the ^1,3MLCT:³T:⁵T₂ excited states).

Despite this limitation, the model for the ^1,3MLCT excited generated from doublet [Fe(2,2′-bipyridine)₃]³⁺ and singlet [Fe(2,2′-bipyridine)₃]²⁺ compounds should be robust since the only distinction is the presence of the electron on the 2,2′-bipyridine ligand which should have minimal impact on the Kβ fluorescence spectrum. For the ³T transient, no long-lived triplet excited state can be used to extract an excited state Kβ fluorescence difference spectrum as an internal reference. Instead, we use the ground state model complex difference spectrum obtained from triplet Fe(ii) phthalocyanine (FePc) and singlet [Fe(2,2′-bipyridine)₃]²⁺ Kβ spectra as our reference difference spectra. We used the four-coordinate FePc, rather than an octahedral model complex, because octahedral Fe(ii) complexes cannot have a triplet ground state. While de Beer et al. have shown that tetrahedral, octahedral, and square planar molecules in the same quintet or sextet spin state have very similar Kβ spectra²⁵, this cannot be demonstrated experimentally for intermediate spin states. Instead, we use theoretical calculations to demonstrate this point. We theoretically calculated the Kβ fluorescence spectra of a four-coordinate square planar and a six-coordinate octahedral ferrous complex using atomic multiplet theory³⁸. This theory is the standard method for calculating and interpreting hard X-ray fluorescence spectra³⁸. For all calculations, the Slater–Condon parameters were reduced to 80% of their atomic value and the 3d orbital and spin angular momentum (LS) coupling was switched off for simplicity. The Kβ spectra were calculated as a 3p1s fluorescence following 1s ionization. For FePc, we use the previously published crystal field parameters (10Dq = 2.7 eV, Ds = 0.86 and Dt = 0.247) in our calculations³⁹. For the six-coordinate octahedral complex calculation, we used a 10Dq = 1.5 eV, consistent with the experimental 10Dq ≈ 1.5 eV measured for [Fe(2,2′-bipyridine)₃]²⁺ (ref. 9). This value also ensures a low spin (S = 0) ground state, a high spin (S = 2) first excited state and an intermediate spin (S = 1) second excited state.

Extended Data Fig. 1a shows the calculated Kβ fluorescence spectra for both the four- and six-coordinate complexes. The square planar and octahedral symmetries have similar triplet state Kβ fluorescence spectra, consistent with prior experimental and theoretical findings for high spin complexes^{19, 25}. The accuracy of the calculations can also be assessed by comparing calculated and experimental difference spectra. In Extended Data Fig. 1b we show a comparison between the calculated difference spectrum generated when subtracting an octahedral crystal field singlet state from the square planar triplet ground state and the experimental difference spectrum generated by subtracting singlet [Fe(2,2′-bipyridine)₃]²⁺ spectrum from the triplet FePc spectrum. The calculated difference spectrum reproduces the qualitative features of the experimental difference spectrum. The insensitivity of the calculated spectra to the coordination geometry and the ability of the calculations to reproduce the main features of the experimental difference spectrum validate the use of the FePc fluorescence spectrum as a model for the triplet excited state of [Fe(2,2′-bipyridine)₃]²⁺.

Using model complex difference spectra has proven more fruitful for the kinetic modelling than singular value decomposition (SVD). The model complex difference spectra demonstrate that differentiation of the ^1,3MLCT and the ³T excited states depends upon both the shape of the difference spectra and the relative amplitudes of the difference spectra. To first order, the integrated area of the Kβ fluorescence spectra do not change with spin state. The integral of the absolute value of the difference spectrum, however, depends linearly on the magnitude of the spin change³⁹. This robust and reproducible aspect of Kβ fluorescence spectroscopy makes the relative amplitudes of the difference spectra an important distinction. SVD, however, struggles to differentiate species when a difference in relative amplitude is a key distinguishing feature of the difference spectra. For this reason, we have used model complex difference spectra, rather than SVD to model the time resolved data.

Kinetic modelling of the [Fe(2,2′-bipyridine)₃]²⁺ experimental population dynamics

We have used two distinct kinetic models to analyse the time-dependent electron dynamics in [Fe(2,2′-bipyridine)₃]²⁺. For the direct transition between ^1,3MLCT and ⁵T₂, without a ³T transient state, the relaxation mechanism can be expressed as follows:

where ^1,3MLCT corresponds to the electronic excited state populated by optical excitation, ⁵T₂ corresponds to the quintet ligand field state, and ¹A₁ represents the electronic ground state. The differential rate equations for each species are given by the following mass balance simultaneous equations,

The integrated rate equations provide the following time-dependent populations for the three species,

From prior ultrafast measurements, we know that the lifetime of the ⁵T₂ excited state is roughly 660 ps (refs 6, 13, 15). The long lifetime of the ⁵T₂ excited state allows us to set k₃ ≈ 0 when we model the kinetics in the first couple of picoseconds. The integrated rate equations can be reduced to:

For the sequential kinetic model with a ³T transient state, the relaxation mechanism can be expressed as follows:

where ^1,3MLCT corresponds to the electronic excited state populated by optical excitation, ³T corresponds to the triplet ligand field excited state, and ⁵T₂ corresponds to the quintet ligand field excited state, and ¹A₁ represents the electronic ground state. The differential rate equations for each species are given by the following mass balance simultaneous equations:

The integrated rate equations provide the following time-dependent populations for the four species:

The long lifetime of the ⁵T₂ excited state allows us to set k₃ ≈ 0 when we model the kinetics in the first couple of picoseconds. The integrated rate equations can be reduced to:

To fit the experimental data to a kinetic model, we must convolve the kinetic model with the instrument response function which we describe with a Gaussian function. Taking the example of , which is an exponential decay starting at time zero (t₀), it will be expressed as

where H is the Heaviside step function and σ is the temporal width of the instrument response function.

Statistical determination of the correct kinetic model

Given two distinct kinetic models, we must determine which model best represents the experimental data. Choosing the model with smaller residual sum of squares (RSS) is not sufficient because the two models do not have the same number of fit parameters. We have used the statistical F-test to determine whether the model with or without a ³T transient provides the best fit of the experimental data⁴⁰.

The F-test provides a statistically robust method for comparing the quality of two models with a different number of fit parameters when the simpler model 1 can be ‘nested’ within the more complex model 2. Model 1 has p₁ parameters, and model 2 has p₂ parameters, where p₂ > p₁. For any choice of parameters in model 1, the model 2 should always be able to fit the data at least as well as the model 1. Thus, model 2 typically will have a lower RSS than model 1. The F-test allows us to determine the statistical significance of this variance in RSS. The F statistic can be calculated by

where n is the number of data points (time delays) fitted by the two models. For the null hypothesis that model 2 does not provide a fit statistically superior to that provided by model 1, the F will have an F distribution defined by the degrees of freedom, (p₂ − p₁) and (n − p₂). To reject the null hypothesis, F must exceed a critical value that depends upon the degrees of freedom and the level of confidence⁴⁰.

[Fe(2,2′-bipyridine)₃]²⁺ experimental data modelling

Using the reference difference spectra with the kinetic model, we fit the time-dependent difference Kβ fluorescence spectra for optically excited [Fe(2,2′-bipyridine)₃]²⁺ in water. The parameters extracted from the fit of the two kinetic models can be found in Extended Data Table 1. We compute the time constants and uncertainties reported in Extended Data Table 1 by fitting multiple runs of the same experiment and then calculating the mean and the standard deviation. The experimental two-dimensional transient difference spectra, fit spectra, residuals, and excited electronic state populations extracted from the best fit for each model can be found in Extended Data Figs 2 and 3. Given the very short lifetime of the ³T excited state, the deviations between the fits of the two models predominantly occur within the first 500 fs. The two-dimensional plots of the residuals in Extended Data Figs 2c and 3c highlight the regions where the ^1,3MLCT³T⁵T₂ model provides a fit superior to that of the ^1,3MLCT⁵T₂ model. Unsurprisingly, this corresponds to time delays with larger ³T populations and spectral regions with the largest difference between the ³T and ⁵T₂ spectra (7,053–7,056 eV).

The residual sum of squares quantifies the variable quality in the fits. The residual sum of squares for each model is: RSS₁ = 3.77 and RSS₂ = 3.21. In this situation, we have p₁ = 5, p₂ = 6 and n = 45. To be 95% confident that the complex model is better than the nested model, the calculated F value must be larger than the F distribution value that captures 95% of the distribution for F(p₂ − p₁, n − p₂) which is 4.09. The calculated F value is 6.71 which exceeded 4.09. So with 95% confidence we conclude that the model containing the ³T transient provides a better description of the experimental data.

Influence of instrument response function parameters on the data analysis

We utilize the instrument response function (IRF) as a variable since the technology does not yet exist to measure the instrument response time accurately. This leads to an increase in the number of parameters in the data analysis. This increase in fit parameters makes statistically differentiating the robustness of alternative kinetic models more difficult, rather than easier.

To ensure that the statistical superiority of the kinetic model possessing the ³T transient does not result from our uncertainty about the instrument response function parameters, we have investigated how variation of time zero and FWHM values differentially influence the RSS for the direct ^1,3MLCT⁵T₂ model and the ^1,3MLCT³T⁵T₂ model. For the range of time zero and FWHM values reported in Extended Data Table 1 that adequately fit the experimental data with either model, the model containing the ³T transient always provides a significantly superior fit to the experimental data. We have used the instrument response function values that minimize the RSS for the ^1,3MLCT⁵T₂ model to fit the data with the ^1,3MLCT³T⁵T₂ model. Using this sub-optimal instrument response function only increases the RSS₂ from 3.21 to 3.27, both significantly less than the direct model RSS₁ = 3.77. Using the definition for F given above and p₁ = 5, p₂ = 6 and n = 45, we calculate F = 5.98, in excess of the 4.09 value needed to conclude with 95% confidence that the complex model provides a better representation of the experimental data than the nested model.

Experimental time resolution can also influence the ability to identify a distinct excited state. For the case of the triplet transient, the temporal resolution of 150 fs has little impact on the characterization of the triplet excited state dynamics. To demonstrate that the roughly 150 fs FWHM IRF does not inhibit our ability to characterize the triplet population dynamics, we have simulated the ^1,3MLCT³T⁵T₂ population kinetics using the time constants extracted from the best fit to the experimental data listed in Extended Data Table 1 with an IRF possessing a 150 fs FWHM and a 5 fs FWHM. The initial time dependence of the ^1,3MLCT state signal depends significantly on the time resolution (though the decays for time delays longer than 200 fs look similar), but the shape and amplitude of the triplet population is similar. The convolution of the IRF and the lifetime of the ^1,3MLCT excited state determine the time dependence of the ³T transient state formation observed experimentally. The low transient population of the triplet state results primarily from the fact that the decay rate of the ³T state exceeds that of the ^1,3MLCT state by a factor of two.

1. Gust, D., Moore, T. A. & Moore, A. L. Mimicking photosynthetic solar energy transduction. Acc. Chem. Res. 34, 40–48 (2001)
   • CAS
   • OpenURL
   • ISI
   • PubMed
   • Article
2. Sato, O., Iyoda, T., Fujishima, A. & Hashimoto, K. Photoinduced magnetization of a cobalt-iron cyanide. Science 272, 704–705 (1996)
   • CAS
   • OpenURL
   • ADS
   • ISI
   • PubMed
   • Article
3. Ferrere, S. & Gregg, B. A. Photosensitization of TiO₂ by Fe-II(2,2′-bipyridine-4,4′-dicarboxylic acid)₂(CN)₂: band selective electron injection from ultra-short-lived excited states. J. Am. Chem. Soc. 120, 843–844 (1998)
   • CAS
   • OpenURL
   • ISI
   • Article
4. Heyduk, A. F. & Nocera, D. G. Hydrogen produced from hydrohalic acid solutions by a two-electron mixed-valence photocatalyst. Science 293, 1639–1641 (2001)
   • CAS
   • OpenURL
   • ADS
   • PubMed
   • Article
5. Goldbeck, R. A., Kim-Shapiro, D. B. & Kliger, D. S. Fast natural and magnetic circular dichroism spectroscopy. Annu. Rev. Phys. Chem. 48, 453–479 (1997)
   • CAS
   • OpenURL
   • ADS
   • PubMed
   • Article
6. Bressler, C. et al. Femtosecond XANES study of the light-induced spin crossover dynamics in an iron(II) complex. Science 323, 489–492 (2009)
   • CAS
   • OpenURL
   • ADS
   • ISI
   • PubMed
   • Article
7. Huse, N. et al. Femtosecond soft X-ray spectroscopy of solvated transition-metal complexes: deciphering the interplay of electronic and structural dynamics. J. Phys. Chem. Lett. 2, 880–884 (2011)
   • CAS
   • OpenURL
   • Article
8. Gutlich, P. & Goodwin, H. A. in Spin Crossover in Transition Metal Compounds I Vol. 233, Topics in Current Chemistry (eds Gutlich, P. & Goodwin, H. A.) 1–47 (Springer, 2004)
9. Creutz, C., Chou, M., Netzel, T. L., Okumura, M. & Sutin, N. Lifetimes, spectra, and quenching of the excited-states of polypyridine complexes of iron(II), ruthenium(II), and osmium(II). J. Am. Chem. Soc. 102, 1309–1319 (1980)
   • CAS
   • OpenURL
   • Article
10. Hauser, A. Intersystem crossing in the Fe(PTZ)₆ (BF₄)₂ spin crossover system (PTZ = 1-propyltetrazole). J. Chem. Phys. 94, 2741–2748 (1991)
    • CAS
    • OpenURL
    • Article
11. McCusker, J. K. et al. Subpicosecond ¹MLCT-⁵T₂ intersystem crossing of low-spin polypyridyl ferrous complexes. J. Am. Chem. Soc. 115, 298–307 (1993)
    • CAS
    • OpenURL
    • Article
12. Monat, J. E. & McCusker, J. K. Femtosecond excited-state dynamics of an iron(II) polypyridyl solar cell sensitizer model. J. Am. Chem. Soc. 122, 4092–4097 (2000)
    • CAS
    • OpenURL
    • ISI
    • Article
13. Gawelda, W. et al. Ultrafast nonadiabatic dynamics of [Fe(II)(bpy)₃]²⁺ in solution. J. Am. Chem. Soc. 129, 8199–8206 (2007)
    • CAS
    • OpenURL
    • PubMed
    • Article
14. Consani, C. et al. Vibrational coherences and relaxation in the high-spin state of aqueous [Fe-II(bpy)₃]²⁺. Angew. Chem. Int. Edn 48, 7184–7187 (2009)
    • CAS
    • OpenURL
    • Article
15. Lemke, H. T. et al. Femtosecond X-ray absorption spectroscopy at a hard X-ray free electron laser: application to spin crossover dynamics. J. Phys. Chem. A 117, 735–740 (2013)
    • CAS
    • OpenURL
    • PubMed
    • Article
16. Emma, P. et al. First lasing and operation of an angstrom-wavelength free-electron laser. Nature Photon. 4, 641–647 (2010)
    • CAS
    • OpenURL
    • ADS
    • ISI
    • Article
17. Harmand, M. et al. Achieving few-femtosecond time-sorting at hard X-ray free-electron lasers. Nature Photon. 7, 215–218 (2013)
    • CAS
    • OpenURL
    • ADS
    • Article
18. Haldrup, K. et al. Guest-host interactions investigated by time-resolved X-ray spectroscopies and scattering at MHz rates: solvation dynamics and photoinduced spin transition in aqueous [Fe(bipy)₃]²⁺. J. Phys. Chem. A 116, 9878–9887 (2012)
    • CAS
    • OpenURL
    • PubMed
    • Article
19. Vankó, G. et al. Probing the 3d spin momentum with X-ray emission spectroscopy: the case of molecular-spin transitions. J. Phys. Chem. B 110, 11647–11653 (2006)
    • CAS
    • OpenURL
    • PubMed
    • Article
20. Krause, M. O. & Oliver, J. H. Natural widths of atomic K and L levels, K-alpha X-ray lines and several KLL Auger lines. J. Phys. Chem. Ref. Data 8, 329–338 (1979)
    • CAS
    • OpenURL
    • ADS
    • ISI
    • Article
21. de Graaf, C. & Sousa, C. Study of the light-induced spin crossover process of the [Fe(II)(bpy)₃]²⁺ complex. Chemistry 16, 4550–4556 (2010)
    • CAS
    • OpenURL
    • PubMed
    • Article
22. de Graaf, C. & Sousa, C. On the role of the metal-to-ligand charge transfer states in the light-induced spin crossover in Fe-II(bpy)₃. Int. J. Quantum Chem. 111, 3385–3393 (2011)
    • CAS
    • OpenURL
    • Article
23. Glatzel, P. & Bergmann, U. High resolution 1s core hole X-ray spectroscopy in 3d transition metal complexes — electronic and structural information. Coord. Chem. Rev. 249, 65–95 (2005)
    • CAS
    • OpenURL
    • Article
24. de Groot, F. High resolution X-ray emission and X-ray absorption spectroscopy. Chem. Rev. 101, 1779–1808 (2001)
    • CAS
    • OpenURL
    • ISI
    • PubMed
    • Article
25. Lee, N., Petrenko, T., Bergmann, U., Neese, F. & DeBeer, S. Probing valence orbital composition with iron K β X-ray emission spectroscopy. J. Am. Chem. Soc. 132, 9715–9727 (2010)
    • CAS
    • OpenURL
    • ISI
    • PubMed
    • Article
26. Cannizzo, A. et al. Light-induced spin crossover in Fe(II)-based complexes: the full photocycle unraveled by ultrafast optical and X-ray spectroscopies. Coord. Chem. Rev. 254, 2677–2686 (2010)
    • CAS
    • OpenURL
    • Article
27. Sousa, C. et al. Ultrafast deactivation mechanism of the excited singlet in the light-induced spin crossover of [Fe(2,2′-bipyridine)₃]²⁺. Chemistry 19, 17541–17551 (2013)
    • CAS
    • OpenURL
    • PubMed
    • Article
28. Alvarez, S. Relationships between temperature, magnetic moment, and continuous symmetry measures in spin crossover complexes. J. Am. Chem. Soc. 125, 6795–6802 (2003)
    • CAS
    • OpenURL
    • PubMed
    • Article
29. Khalil, M. et al. Picosecond X-ray absorption spectroscopy of a photoinduced iron(II) spin crossover reaction in solution. J. Phys. Chem. A 110, 38–44 (2006)
    • CAS
    • OpenURL
    • PubMed
    • Article
30. Nozawa, S. et al. Direct probing of spin state dynamics coupled with electronic and structural modifications by picosecond time-resolved XAFS. J. Am. Chem. Soc. 132, 61–63 (2010)
    • CAS
    • OpenURL
    • PubMed
    • Article
31. Gawelda, W. Time-Resolved X-Ray Absorption Spectroscopy of Transition Metal Complexes. Ph.D. thesis, École Polytechnique Fédérale de Lausanne. (2006)
32. Feng, Y. P. et al. A single-shot intensity-position monitor for hard X-ray FEL sources. Proc. SPIE 8140, 81400Q (2011)
33. Alonso-Mori, R. et al. A multi-crystal wavelength dispersive X-ray spectrometer. Rev. Sci. Instrum. 83, 9 (2012)
    • CAS
    • OpenURL
    • Article
34. Koerner, L. J., Philipp, H. T., Hromalik, M. S., Tate, M. W. & Gruner, S. M. X-ray tests of a pixel array detector for coherent X-ray imaging at the Linac Coherent Light Source. J. Instrum. 4, P03001 (2009)
    • OpenURL
    • Article
35. Bionta, M. R. et al. Spectral encoding of X-ray/optical relative delay. Opt. Express 19, 21855–21865 (2011)
    • CAS
    • OpenURL
    • PubMed
    • Article
36. Vankó, G. et al. Picosecond time-resolved X-ray emission spectroscopy: ultrafast spin-state determination in an iron complex. Angew. Chem. Int. Edn 49, 5910–5912 (2010)
    • CAS
    • OpenURL
    • Article
37. Vankó, G. et al. Spin-state studies with XES and RIXS: From static to ultrafast. J. Elec. Spec. Relat. Phenom. 188, 166–171 (2013)
    • CAS
    • OpenURL
    • Article
38. de Groot, F. M. F. & Kotani, A. Core Level Spectroscopy of Solids (CRC Press, Boca Raton, 2008)
39. Stepanow, S. et al. Mixed-valence behavior and strong correlation effects of metal phthalocyanines adsorbed on metals. Phys. Rev. B 83, 220401 (2011)
    • CAS
    • OpenURL
    • Article
40. Kutner, M. H., Nachtsheim, C. J. & Neter, J. Applied Linear Regression Models (McGraw-Hill/Irwin, 2004)

Download references

We thank P. Frank, B. Lin and S. DeBeer for discussion, S. DeBeer for some model iron complex X-ray fluorescence spectra, and D. Stanbury for providing some iron complexes. Experiments were carried out at LCLS and SSRL, which are National User Facilities operated for DOE and OBES respectively by Stanford University. W.Z., R.W.H., H.W.L., D.A.M., Z.S. and K.J.G. acknowledge support from the AMOS programme within the Chemical Sciences, Geosciences and Biosciences Division of the Office of Basic Energy Sciences, Office of Science, US Department of Energy. E.I.S. acknowledges support from the NSF (CHE-0948211). R.G.H. acknowledges a Gerhard Casper Stanford Graduate Fellowship and the Achievements Rewards for College Scientists (ARCS) Foundation. T.K. acknowledges the German Research Foundation (DFG), grant KR3611/2-1. K.S.K., M.M.N. and T.B.v.D. acknowledge support from the Danish National Research Foundation and from DANSCATT. K.K. thanks the Volkswagen Foundation for support under the Peter Paul Ewald fellowship program (I/85832). G.V. acknowledges support from the European Research Council (ERC-StG-259709) and the Lendület Programme of the Hungarian Academy of Sciences. C.B., W.G. and A.G. thank the DFG (SFB925), as well as the European XFEL, for financial support.