The structure and spin of photoexcited Fe2+(phen)3 in water are examined by x-ray scattering and x-ray emission spectroscopy with 100 ps time resolution. Excitation of the low-spin (LS) ground state (GS) to the charge transfer state 1MLCT* leads to the formation of a high-spin (HS) state that returns to the GS in 725 ps. Density functional theory (DFT) predicts a Fe–N bond elongation in HS by 0.19 Å in agreement with the scattering data. The angle between the ligands increases by 5.4° in HS, which allows the solvent to get 0.33 Å closer to Fe in spite of the expansion of the molecule. The rise in solvent temperature from the return of photoproducts to the GS is dominated by the formation dynamics of HS, 1MLCT* → HS, which is followed by a smaller rise from the HS → GS transition. The latter agrees with the 0.61 eV energy gap E(HS)−E(LS) calculated by DFT. However, the temperature rise from the 1MLCT → HS transition is greater than expected, by a factor of 2.1, which is explained by the re-excitation of nascent HS* by the 1.2 ps pump pulse. This hypothesis is supported by optical spectroscopy measurements showing that the 1.2 ps long pump pulse activates the HS* → 5MLCT* channel, which is followed by the ultrafast return to HS* via intersystem crossing. Finally, the spins of the photoproducts are monitored by the Kβ emission and the spectra confirm that the spins of LS and HS states are 0 and 2, respectively.

Electron transfer in transition metal complexes can convert photon energy into electric or chemical energy which has found applications in solar cells and photocatalysis.1 In photosynthesis, the transfer of an electron from a metal to a ligand is responsible for the creation of a chemical potential that activates ATP, the energy carrier in biological systems.2 In photocatalysis, the light driven splitting of water into dioxygen and dihydrogen, catalyzed by a photoactivated Ru(II) complex, is an example of the strive toward cleaner and renewable energy.3 

In transition metal complexes with open 3d, 4d, or 5d shells, a d electron is transferred upon photon absorption to a ligand orbital, and this state is called the metal-to-ligand charge transfer (MLCT) state. The process is reversible and the electron returns to the d orbital of the metal. In 4d and 5d complexes, the LMCT lifetime is long, from nanoseconds to microseconds. In 3d complexes it is of the order of sub-picoseconds in most cases. In 4d and 5d complexes, the outer electron density is predominantly 4d or 5d in character and the ligand bonding is stronger than for 3d complexes. The ligand field splitting of 4d or 5d orbitals is therefore greater than for 3d complexes. In fact, the lowest excited state for 4d and 5d complexes is the LMCT state whereas for 3d complexes it is a metal centered (MC) ligand field state. In the decay of the MLCT state in 3d complexes, the electron returns to the ground state (GS) via the MC state which greatly shortens the LMCT lifetime compared with the equivalent 4d and 5d systems.4 In the 4d and 5d complexes Ru2+(bpy)3 and Os2+(bpy)3, the MLCT lifetimes are 900 (Ref. 5) and 25 ns (Ref. 6), respectively, whereas it is 0.96 ps for Fe2+(bpy)3 .7 The short MLCT lifetime of 3d complexes makes them unsuitable for light conversion since the charge separated state does not have time to interact with the environment. Since 4d and 5d metals are rare and expensive, it is important to synthetize ligands for 3d systems with the strongest possible ligand field. The ligand field strength depends on the electronegativity of the ligand atoms in the metal–ligand bonds: O and halogens produce weak fields, N intermediate fields, and C and P strong fields. In the carbene complex Fe3+(btz)3(PF6)2, the strong Fe–C σ bonds generate a high electron density on Fe3+ and the MC states are inaccessible in the decay process. The MLCT state decays directly to the GS by photon emission with a long 528 ps lifetime.8 Another example is the Fe3+(phtmeimb)2 complex, which has a record MLCT lifetime of 2 ns.9 

In this work we present a time resolved study of the change in structure and spin of the HS state of Fe2+(phen)3 (phen = 1,10-phenanthroline) in water using wide-angle x-ray scattering (TRWAXS) and x-ray emission spectroscopy (TRXES). The aim is to measure the change in the average Fe–N bond length to the ligands, the relative ligand angles, the spin states of Fe as well as the cage structure and hydrodynamics of the solvent. The work was performed by laser/x-ray pump probe using 100 ps x-ray pulses from the European Synchrotron Radiation Facility (ESRF).10,11 The experimental details are given in Sec. II and in the supplementary material.35 

Fe2+(phen)3 has the stoichiometry Fe2+(N2C12H8)3 with 67 atoms per molecule. Phenanthroline, a medium strong ligand, has two Fe–N bonds per ligand in Fe2+(phen)3. The LS and HS structures were calculated by density functional theory (DFT) using the B3LYP function and the 6-311++G(3df,3pd) basis set. The effect of the solvent was included using the polarizable continuum model (CPCM) with Fe2+(phen)3 placed in a cavity in the solvent field. The LS structure was assumed to have D3 symmetry in the calculation. The LS calculation predicts six equal Fe–N bond lengths of 1.99 Å and relative ligand angles of 88.8°. The HS state was assumed to have C2 symmetry. The Fe–N bonds and the relative ligand angles were calculated to be 2.18 Å and 94.2° on average. The LS structure is shown in Fig. 1(a). The atomic positions and structural parameters are listed in Appendix 1 in the supplementary material.35 

FIG. 1.

(a) Structure of Fe2+(phen)3 with Fe2+ in red, N in blue, and C in gray. H atoms are omitted. (b) In an octahedral ligand field, the five basic 3d wavefunctions are split into t2g and eg energy levels. (c) LS and HS configurations for the six 3d electrons in Fe2+.

FIG. 1.

(a) Structure of Fe2+(phen)3 with Fe2+ in red, N in blue, and C in gray. H atoms are omitted. (b) In an octahedral ligand field, the five basic 3d wavefunctions are split into t2g and eg energy levels. (c) LS and HS configurations for the six 3d electrons in Fe2+.

Close modal

The 3d orbitals in the ligand field states t2g and eg are shown in Fig. 1(b) and the LS and HS configurations are shown in Fig. 1(c).

In the ground state t2g, the lopes of the 3dxy, 3dxz, and 3dyz wavefunctions are pointing along the diagonals in the (x, y), (y, z), and (z, x) planes towards the midpoints of the N-N lines of the ligands where the repulsive energy is lowest. In eg the dx2-y2 and dz2 lopes are pointing along the x, y, and z axes in the direction of the N atoms. The repulsive energy is higher and the Fe–N bonds are longer.

The electron configuration of Fe2+ is [Ar]3d6. The 3d orbitals are tenfold degenerate and can accommodate five up-spins and five down-spins. The six 3d electrons in Fe2+are distributed in t2g and eg according to Hund's rule: maximize the total spin at the lowest possible total energy. If the ligand field splitting is greater than the spin-pairing energy in an orbital, the case for Fe2+(phen)3, the three orbitals in t2g are filled with three up- and three down-spins. The total spin is S = 0. For the next 7–10 electrons in the following 3d metals, eg is filled by maximizing the spin. For ten 3d electrons, the total spin S is  0. If the ligand field splitting is smaller than the spin-pairing energy, the first five electrons have parallel spins in t2g and eg. The next five electrons are then put in t2g and eg with antiparallel spins.

As mentioned, phenanthroline is a medium strong ligand and the GS of Fe2+ (phen)3 has six 3d electrons in t2g with antiparallel spins and S = 0. When the eg orbitals are populated in HS at the end of the photocycle, eg receives the recombining electron which triggers the transfer of one t2g electron to eg with a change in spin. The two t2g electrons have parallel spins and the total spin S is  2.

The photocycle of Fe2+(phen)3 in acetonitrile was studied by ultrafast optical spectroscopy by Tribollet et al.12 and Gallé et al.,13 by time resolved EXAFS in aqueous solution by Nozawa et al.,14 and by XANES by Wang et al.15 The Fe–N bond lengths derived from the EXAFS work are 1.98(1) and 2.15(1) Å for the LS and HS states, respectively. The HS lifetime was reported to be 690 and 685 ps in the EXAFS and XANES works, respectively. For comparison, the HS lifetimes of the sister complexes Fe2+(terpy)2 and Fe2+(bpy)2 in aqueous solutions are 2.61 ns (Ref. 34) and 0.665 ns (Ref. 7), respectively, and the Fe2+(btbip)2 complex has a HS lifetime in acetonenitrile of 0.23 ns (Ref. 16).

The Fe2+(phen)3 wavefunctions, transition rates, and energy curves were calculated in acetonitrile by Sousa et al.17 and the energy curves and spin states are summarized in Fig. 2. The photocycle proceeds as follows. When Fe2+(phen)3 is excited from the GS 1A1, a 3d electron in t2g moves to one of the ligands to form the 1MLCT state. The 1MLCT transits in 36 fs to the triplet 3MLCT that in turn decays to the MC triplet 3T1 in 150 fs. 3T1 decays in turn to 5T2 in 102 fs by intersystem crossing (ISC). The 5T2 state is reached in ∼300 fs. 5T2 is initially in an excited vibrational state that relaxes to the bottom of the potential in a few picoseconds via vibrational cooling by the solvent. The final step 5T2 → 1A1 proceeds via ISC in 1.1 ns in acetonitrile. Note that 3T1 and 5T2 have one and two antibonding electrons in eg, respectively, which explains the progressive elongation of the Fe–N bonds.

FIG. 2.

(a) Calculated energy curves and spin states for Fe2+(phen)3 vs the Fe–N bond distance. The data are adopted from Sousa et al.17 

FIG. 2.

(a) Calculated energy curves and spin states for Fe2+(phen)3 vs the Fe–N bond distance. The data are adopted from Sousa et al.17 

Close modal

The HS(S = 2) state was first observed by Tribollet et al.12 Their Fe2+(phen)3/acetonitrile solution was excited by 30 ps pulses at 540 nm and probed by 30 ps pulses of white light. They observed a transient absorption feature around 340 nm with a 1.1 ns lifetime that they assigned to the HS state 5T2. The absorbing feature was assigned to the 5T2 → 5MLCT excitation.

TRWAXS is used in this work to monitor the structural changes in the LS to HS transition, i.e. the change in Fe–N bond lengths, the relative ligand angles, the effective radius of the solvent cage, and the change in the solvent temperature. As 1D scattering data cannot be inverted to electron density directly, the measured scattering curves were compared with the predictions from combined DFT and MD simulations. The change in the spin state of Fe2+ was probed by TRXES from the Kβ emission line. The experiments were performed on beamline ID09 at the European Synchrotron Radiation Facility (ESRF) and the time resolution was obtained by pump and probe. A 10 mM aqueous solution of Fe2+(phen)3 was circulated in a liquid jet, excited by 1.2 ps pulses at 400 nm and probed by 100 ps x-ray pulses from an undulator. In the TRWAXS study, 18 keV x-ray pulses were produced by a 2.1% bandwidth (BW) Ru multilayer (ML) to maximize the intensity while keeping the beam quasi monochromatic. The intensity of the ML beam was 4 × 108 ph/pulse and the pump & probe pulses were repeated at 986.3 Hz. In the TRXES experiment, the 5% BW pink beam from the fundamental of the U17 undulator at 15 keV was used to generate 1s core holes on the Fe K-edge. The intensity of the pink beam was 2 × 109 ph/pulse. With the 100 ps time resolution from the x-ray pulse length, the HS state is the only structure resolved in the present experiment. The setup and the experimental parameters are described in more detail in the supplementary material.35 

The structural information from TRWAXS data is obtained by fitting candidate structures from DFT/MD simulations against the laser induced change in the scattering ΔS(q, t). The LS and HS structures were calculated by density functional theory (DFT) as mentioned earlier. From the DFT atomic positions, the scattering from randomly oriented Fe2+(phen)3 molecules was calculated by the Debye function as follows:18 
(1)
where rmn is the distance between atom pairs and fn(q) is the atomic form factor of atom n. The forward intensity S(q = 0) is the coherent intensity from all the electrons in the molecule (Σ Zi)2, with the sum extending over all atoms in the molecule and Zi being the number of electrons in atom i. Fe2+(phen)3 has 306 electrons per molecule, S(q = 0) = (306 e.u.)2 = 93 636 (e.u.)2 where e.u. is the electron unit defined by the atomic form factor.18 As the stoichiometry of the LS and HS structures does not change, the change in the forward scattering ΔS(q = 0) is zero. In contrast, if a molecule dissociates or associates, ΔS(q = 0) is negative or positive, respectively.

The scattering functions S(q) for the LS and HS structures are shown in Fig. 3(a) and qΔS(q) are shown in Fig. 3(b). q multiplication is used to amplify the high-q spectrum, which is particularly sensitive to intra molecular structural changes. Note that qΔS(q) has 12 nodal points in the range 0 < q < 10.7 Å−1, which is a very stringent constraint for the comparison of calculated and measured scattering curves.

FIG. 3.

Scattering functions for Fe2+(phen)3. (a) The total scattering for the DFT structures of the LS and HS states. (b) Difference scattering qΔS(q) (HS-LS). The negative part between 0 and 1 Å−1 is the signature of the Fe–N bond extension in the HS state. (c) Fourier transforms rS[r] and rΔS[r]. The assignments of the LS and HS contributions are indicated on the HS-LS curve.

FIG. 3.

Scattering functions for Fe2+(phen)3. (a) The total scattering for the DFT structures of the LS and HS states. (b) Difference scattering qΔS(q) (HS-LS). The negative part between 0 and 1 Å−1 is the signature of the Fe–N bond extension in the HS state. (c) Fourier transforms rS[r] and rΔS[r]. The assignments of the LS and HS contributions are indicated on the HS-LS curve.

Close modal
The change in the radial electron density for an average excited atom is obtained from the Fourier transform,19 
(2)
ΔS[r,t] is the form factor biased change in the radial electron density which is biased by the high-Z atom in the molecule. The square brackets in S[r] is used to show that it is the Fourier transform of S(q). The rS[r] and rΔS[r] functions are shown in Fig. 3(c). The green curve rΔS[r] has a negative peak at 1.87 Å and a positive at 2.29 Å which is the signature of the Fe–N bond elongation in the HS state. There is a small difference between the r values of the peak and valley positions compared with the DFT values of 1.99 and 2.18 Å which is due to the interference of positive and negative peaks in the difference curve.

The measured curves ΔS(q, t) have two principal contributions: one from the change in the solute/cage structure and one from the change in the solvent temperature from the relaxation of the photoproducts. The measured qΔS(q, t) curves for 100 ps and 2 ns are shown in Fig. 4. The blue curve in Fig. 4(a) is the total change in the solute/cage/solvent scattering (HS-LS). The solvent part is produced by the adiabatic temperature rise in the solvent after 100 ps. The orange curve is the solvent contribution that was measured in the separate dye/water experiment described in the supplementary material.35 The difference between the curves is the intrinsic change in the solute and solute cage structure. The negative part between 0.5 and 1.0 Å−1 in the blue curve is the signature of the Fe-N bond extension in HS. The rise in the solvent temperature after 100 ps is produced by the ultrafast cascade 1MLCT → 3MLCT → 3T2 → 5T2 shown in Fig. 2. For the 2 ns curves in Fig. 4(b), the total solute/cage/solvent curve in blue coincides with the dye/water curve: the HS state has returned to the GS and the residual signal is from the adiabatic temperature rise of the solvent. Note that at fast timescales t < 10 ns, the solvent temperature rises monotonically and that the rise accumulates all relaxation processes in the photocycle, including the ultrafast ones that are not resolved structurally with 100 ps time resolution.

FIG. 4.

qΔS(q,t) scattering curves for time delays of 100 ps and 2 ns. (a) The blue curve probes the total difference scattering from the solute, cage, and solvent. The orange curve is the solvent part alone. (b) After 2 ns, the blue and orange curves coincide as the HS population has returned to the GS; the residual orange signal is purely thermal.

FIG. 4.

qΔS(q,t) scattering curves for time delays of 100 ps and 2 ns. (a) The blue curve probes the total difference scattering from the solute, cage, and solvent. The orange curve is the solvent part alone. (b) After 2 ns, the blue and orange curves coincide as the HS population has returned to the GS; the residual orange signal is purely thermal.

Close modal

The normalization of the intensity is done by scaling the non-excited total scattering to 1.0 at q0 = 6.7 Å−1 and using that calibration for the difference scattering. The total scattering from a 10 mM “unit-cell” solution of Fe2+(phen)3 in water, elastic plus Compton scattering, is 13.3 (e.u.)2/unit.20,21

The experimental data were analyzed with the following terms:
(3)
where γ(t) is the HS population vs time. The solute term is the Debye scattering from the DFT structures and the cage term is calculated from the atom–atom distribution functions gij(r) from the MD simulation,27 
(4)
where (n,m) represents all solute–solvent pairs. The solvent term has two contributions,22–24 
(5)
The second term in Eq. (5) is not needed here since the timescale in the experiment is limited to 0–2 ns where the solvent density is constant.

The MD simulations of the solvent structure around LS and HS states were performed by MOLDY employing the Lennard–Jones potentials from the TIP4P/2005 model of water.25 In the simulation, the LS and HS structures of Fe2+(phen)3 were placed in a (49.6 Å)3 box with 4096 water molecules. The atom–atom distribution function gFe,Ow(r), where OW is oxygen in water, is a measure of the effective cage radius. The gFe,Ow(r) and gN,Ow(r) functions are shown in Fig. 5.

FIG. 5.

(a) The MD calculated distribution function g(Fe, Ow) showing the position of oxygen in water in the LS and HS states. The oxygen in water moves about 1.0 Å closer to Fe in HS in spite of the 0.19 Å Fe–N bond elongation in HS. (b) The g(N, Ow) function shows as well that Ow moves closer to N in HS.

FIG. 5.

(a) The MD calculated distribution function g(Fe, Ow) showing the position of oxygen in water in the LS and HS states. The oxygen in water moves about 1.0 Å closer to Fe in HS in spite of the 0.19 Å Fe–N bond elongation in HS. (b) The g(N, Ow) function shows as well that Ow moves closer to N in HS.

Close modal

The first peak in the LS function gFe,Ow(r) at 5.5 Å shifts to ∼4.5 Å in HS. The latter has three sub-peaks separated by about 0.4 Å. The splitting is caused by the distortion of the HS structure, in particular the slightly different relative ligand angles in HS that vary from 90.4° to 96.0°. When the Fe2+(phen)3 structure is rendered with 3D van der Waals spheres, the increased solvent accessibility to Fe in HS is evident, see Fig. 6. This accounts for the counterintuitive fact that water moves closer to Fe in HS in spite of the Fe–N elongation.

FIG. 6.

(a) Position and angle of the ligand L1 in the LS and HS states from the DFT calculations. Fe is shown in orange, N in blue, and C in gray. For L1 the Fe–N bond expands from 1.99 to 2.18 Å in LS and HS respectively. The N–Fe–N angle ϕ gets smaller in the expansion, from 82.6° to 76.93°. (b) The angles between adjacent ligand planes increase by 5.4° on average in HS. The solvent moves closer to Fe in HS as illustrated by the orange van der Waals sphere.

FIG. 6.

(a) Position and angle of the ligand L1 in the LS and HS states from the DFT calculations. Fe is shown in orange, N in blue, and C in gray. For L1 the Fe–N bond expands from 1.99 to 2.18 Å in LS and HS respectively. The N–Fe–N angle ϕ gets smaller in the expansion, from 82.6° to 76.93°. (b) The angles between adjacent ligand planes increase by 5.4° on average in HS. The solvent moves closer to Fe in HS as illustrated by the orange van der Waals sphere.

Close modal

A splitting of the gFe,Ow(r) function in HS was also observed in MD simulations of Fe2+(bpy)3 in water by the Hauser group26 although the effect is less pronounced.

The solute, cage, and solvent contributions to the total scattering are shown in Fig. 7(a) together with the time dependence of the HS population γ(t) shown in Fig. 7(b). The red curve in Fig. 7(b) is a Gaussian–exponential fit to the HS population with a 725 ps time decay and a 100 ps probe pulse (FWHM). The onset and decay are fitted to the convolution of an exponential decay and a Gaussian x-ray pulse as,
(6)
where τ is the laser/x-ray delay, τR is the HS time constant, and σ the duration of the x-ray pulse. The extrapolated fraction of HS molecules at time zero is 24.3% of the 10 mM solution.
FIG. 7.

(a) The total scattering qΔS(q,100 ps) is shown with dots and the solute, cage, and solvent contributions are shown in red, green, and blue, respectively. The purple curve is the sum of components. The amplitude of the cage term is reduced to 33% of the MD value to fit the data optimally. The upturn in the cage term at low q is the signature of the cage contraction in HS. (b) The HS population inferred from the solute amplitude in Fig. 7(a) as a function of time. The red line is the fit to the data with an exponential 725 ps time constant and 100 ps x-ray pulse (FWHM).

FIG. 7.

(a) The total scattering qΔS(q,100 ps) is shown with dots and the solute, cage, and solvent contributions are shown in red, green, and blue, respectively. The purple curve is the sum of components. The amplitude of the cage term is reduced to 33% of the MD value to fit the data optimally. The upturn in the cage term at low q is the signature of the cage contraction in HS. (b) The HS population inferred from the solute amplitude in Fig. 7(a) as a function of time. The red line is the fit to the data with an exponential 725 ps time constant and 100 ps x-ray pulse (FWHM).

Close modal

To optimize the fit of the experimental data to the DFT/MD model, the amplitude of ΔScage(q) is reduced to 33% of the MD calculated value. MD tends to overestimate the cage amplitude which is due to the fact that MD simulations are using rigid solute structures that allows the solvent to get closer to the solute.

The amplitude of the total scattering curves ΔS(q,t) decreases exponentially with a 725 ps time constant. The shape of the spectra vs. time does not change which is consistent with the fact the HS → LS transition is ultrafast for single molecules. By contrast, the amplitude decreases vs time from the decay of the population of the HS state.

Note that the cage signal is greatest between 0-2.5 Å−1. and that the positive part signals the contraction of the cage. In the low q range 0.3 – 0.8 Å−1, a small shift dr in the atom–atom distribution leads to a proportional change in the amplitude of ΔS(q) as shown in the supplementary material.35 In that approximation, the cage radius contracts by 0.33 Å in response to the metal–ligand elongation in the HS state.

The cage dynamics for Fe2+(bpy)3 in water was studied by TRWAXS by Khakhulin et al.27 and Haldrup et al.28 In Khakhulin's work, the solute and cage were calculated by ab initio MD (AIMD), which is a more realistic account of the solute/cage structure in water. In this complex the Fe–N bond extends from 1.99 Å (LS) to 2.18 Å(HS) and the AIMD calculation showed that the first peak in g(Fe, Ow)moves from 5.29 Å in LS to 5.51 Å in HS corresponding to a 0.22 Å expansion of the cage. The difference in the cage dynamics in Fe2+(phen)3 and Fe2+(bpy)3 is ascribed to the difference in the geometry of the ligands and the different structural distortions in the HS states.

As the temperature term dominates the dS(q,t) curves it is essential to have accurate solvent differentials to isolate the solute/cage contribution from the thermal background. This was done in a separate water/dye experiment performed with the same experimental conditions as the main experiment. The dye experiment is described in the supplementary material.35 

The amplitudes of the solvent and solute/cage terms have to be consistent with the model of the reaction pathway, i.e. the rise in the solvent temperature has to match the change in energy of the photoproducts and their concentrations. The solvent temperature rise will now be examined. With the 100 ps solvent differential from the dye/water measurement, the temperature rise ΔT(t) is obtained from the equation
(7)
and the result is shown in Fig. 8.
FIG. 8.

(a) Solvent temperature rise ΔT(t) during the photocycle. In the time range −50 to 50 ps, the steep rise ΔT0 is assigned to the ultrafast relaxation 1MLCT → 3MLCT → 3T25T2. These transitions are broadened by the 100 ps x-ray pulse. The slow rise from 50 ps to 2 ns is from the 725 ps transition 5T2(HS) → GS(LS). ΔT0 is a factor 2.1 greater than expected from the energy steps in the reaction.

FIG. 8.

(a) Solvent temperature rise ΔT(t) during the photocycle. In the time range −50 to 50 ps, the steep rise ΔT0 is assigned to the ultrafast relaxation 1MLCT → 3MLCT → 3T25T2. These transitions are broadened by the 100 ps x-ray pulse. The slow rise from 50 ps to 2 ns is from the 725 ps transition 5T2(HS) → GS(LS). ΔT0 is a factor 2.1 greater than expected from the energy steps in the reaction.

Close modal
The temperature curve is fitted to the expression
(8)
where H(t′) is the Heaviside step function, ΔT0 the transient rise from the ultrafast cascade to 5T2, ΔT1 the temperature increase from the depopulation of HS, τR the HS lifetime time, and σ is the length of the x-ray pulse. The Gaussian pulse length is σ = 42.5 ps for the 100 ps x-ray pulse (FWHM). The fit gives ΔT0 = 330 mK, ΔT1 = 34 mK, and ΔTtot = 364 mK.

The slow rise ΔT1 = 34 mK is consistent with the expected rise from a 2.4 mM excited solution at time zero and that the 5T2 → GS decay deposits 0.62 eV of energy to the solvent. The excited state concentration derived from ΔT0, which was expected to deposit 2.48 eV (3.10–0.62) of energy per molecule, is 5.1 mM. If the photoexcitation was a single-excitation event in this experiment, the two concentrations should be equal. We will now examine the discrepancy in the temperature rise since this effect has been observed in other complexes as well.

The excess solvent temperature rise might be explained by the fact that the 400 nm 1.2 ps pump pulse re-excites nacent 5T2*(HS) molecules to the 5MLCT band during the vibrational relaxation toward the bottom of the 5T2 potential. Note that the formation of the 5T2 state was indeed monitored by the 5T2 → 5MLCT transition at 340 nm (3.65 eV) in acetonitrile by Tribollet et al.12 Our 400 nm (3.10 eV) pump pulse might therefore excite 5T2* into the 5MLCT band in the early stage of the vibrational relaxation process. Note that the excitation-during-relaxation technique was used by C. B. Harris et al. in their classical work that probed the vibrational relaxation of I2* in CCl4 and in other solvents.29 The 5T2* → 5MLCT transition is spin allowed and the transient absorption should increase for a 400 nm (3.10 eV) pump pulse once the molecule has relaxed to match the 3.10 eV excitation window. The 5MLCT state might then transit to 1MLCT by ISC whereby the photocycle resumes. The re-excitation pathway is illustrated in Fig. 9(b).

FIG. 9.

(a) Transient optical absorption of Fe2+(phen)3 in water after 400 nm excitation with a 150 fs pulse. (a) For the 500 nm probe pulse curve in blue, the decay around t = 0 is due to the depletion of the GS. For the 400 nm probe pulse curve in red, the peak at 0.19 ps is assigned to the 5T2*5MLCT transition. (b) Energy curves showing the proposed 5T2* excitation to the 5MLCT band followed by the transit to 1,3MLCT by ISC.

FIG. 9.

(a) Transient optical absorption of Fe2+(phen)3 in water after 400 nm excitation with a 150 fs pulse. (a) For the 500 nm probe pulse curve in blue, the decay around t = 0 is due to the depletion of the GS. For the 400 nm probe pulse curve in red, the peak at 0.19 ps is assigned to the 5T2*5MLCT transition. (b) Energy curves showing the proposed 5T2* excitation to the 5MLCT band followed by the transit to 1,3MLCT by ISC.

Close modal

To examine this hypothesis, an optical–optical experiment was performed in Mathias Bargheer's laboratory at the University of Potsdam. A 2 mM water solution of Fe2+(phen)3 was excited by 400 nm 150 fs pump pulses and probed with 150 fs pulses at 400 nm and 500 nm. The results are shown in Fig. 9(a).

When the solution is probed at 500 nm near the 512 nm absorption peak, the absorption decreases with time from the depletion of the GS. For the 400 nm probe however, a Gaussian-like peak with a FWHM of 0.36 ps appears 0.19 ps after the pump pulse. The peak is tentatively assigned to the transient absorption from the 5T2* → 5MLCT excitation with 5MLCT returning to 1MLCT by ISC. A second photocycle turn is thereby initiated. Double-excited molecules dissipate 2 × 3.1 eV to the solvent, triple-excited 3 × 3.1 eV and so on.

For the 400 nm probe, the amplitude of the depletion hole and the creation peak is −0.064 and +0.098, respectively, which means that the cross section σ(5T2* → 5MLCT) is 53% greater than σ(1A1 → 1MLCT). That shows the importance of the proposed re-excitation pathway.

While WAXS probes atom–atom distances, intra- and extra-molecular alike, WAXS is insensitive to the spin of the molecules. In contrast, XES is atom-specific and a powerful probe of the spin and valence states of transition metals. As shown below the XES signals from metal complexes in solution are feeble but essentially background free. XES is used here to show that the spin of Fe2+ is 0 and 2 in the LS and HS states respectively.

Upon the absorption of an x-ray photon at the K-edge of Fe2+, a 1s electron is ejected and leaving behind a 1s core hole. The core hole is filled by an outer shell electron accompanied by the emission of either an x-ray photon or an Auger electron. For the K-edge of Fe, the branching ratio is 34% x-ray emission and 66% electron emission.30 For the Kα and Kβ lines, the emission is produced by the transitions 2 p 1 s and 3 p 1 s , respectively. The spin–orbit coupling of the 2p state splits Kα into two sub-lines Kα2 and Kα1 from the transitions 2p1/2 (2e-) → 1s and 2p3/2 (4e-) → 1s with the parenthesis showing the electron occupancy in that state. The Kα2 and Kα1 lines peak at 6390.8 and 6403.8 eV, respectively. The energy splitting of the Kβ emission from 3p1/2 (2e-) → 1s (Kβ′) and 3p3/2 (4e-) → 1s (Kβ1,3) is dominated by the exchange interaction between the 3p core hole and the 3d electrons. The spin–orbit interaction for Kβ is much weaker than for Kα, about 0.7 eV (Ref. 31) due to the greater shielding of the nucleus. The central energy of Kβ is 7058 eV for Fe.

XES is insensitive to the solvent hydrodynamics and the technique is ideal for monitoring population kinetics. One nice feature of XES compared with other spectroscopies is that the core holes can be produced by the pink beam which is ∼200 times stronger than the beam from a Si monochromator. The high beam intensity of 2 × 109 ph/pulse from the undulator was essential in the present experiment due the low pump–probe frequency of 986.3 Hz.

The spin sensitivity of the Kβ emission is due to the strong exchange interaction between the 3p core holes and the 3d electrons. The different spin states of Fe, singlet, doublet, …, quintet, have a unique Kβ spectrum which is essentially independent of the ligand and the solvent.32 The valence-to-core (VtC) transitions 3d → 1s probe the ligand-modified 3d levels and VtC can therefore provide valuable information about the metal-to-ligand bonding.33 However VtC is two-orders of magnitude weaker than Kβ and too weak to be used here due to the low pump–probe frequency (986.3 Hz).

The Kβ spectra ON and OFFand their HS-LS difference are shown in Fig. 10.

FIG. 10.

(a) Kβ spectra for laser excited (ON, red) and non-excited (OFF, blue) Fe2+(phen)3 molecules. The ON data was recorded 100 ps after excitation. The difference spectrum ON–OFF is shown in green. (b) XES lines for Fe2+. The 3d spin configuration is shown in the exchange split HS state S = 2.

FIG. 10.

(a) Kβ spectra for laser excited (ON, red) and non-excited (OFF, blue) Fe2+(phen)3 molecules. The ON data was recorded 100 ps after excitation. The difference spectrum ON–OFF is shown in green. (b) XES lines for Fe2+. The 3d spin configuration is shown in the exchange split HS state S = 2.

Close modal

The shape of the difference spectrum HS(100 ps)−LS has the “pure” shape of the LS (S = 0) → HS (S = 2) transition since the non-excited molecules cancel out in the difference spectrum. The amplitude is proportional to the degree of excitation. The ON spectrum is shifted +0.2 eV from the OFF peak at 7058.3 eV. The appearance of the Kβ peak at 7044.6 eV is due to the increase in exchange splitting of the 3p level in the magnetic HS state. The VtC peak at 7110 eV is also detected; it is 100 times weaker than Kβ.

The K′β and Kβ1,3 peaks in the HS spectrum are at 7045 and 7061 eV, respectively, corresponding to an exchange/spin–orbit splitting of 16 eV. The spectra in Fig. 11 were fitted to theoretical LS and HS spectra calculated by the CRISPY software. To match the experimental data, the calculated spectra were shifted +53 eV and broadened by 1.6 eV Lorentzians. The XES spectra from Fe2+(phen)3 in water resemble closely those for Fe2+(terpy)2 in water as reported by G. Vanko et al.34 

FIG. 11.

Measured and calculated spectra of the Kβ emission from Fe2+(phen)3 in water. The curves in orange were calculated by CRISPY.

FIG. 11.

Measured and calculated spectra of the Kβ emission from Fe2+(phen)3 in water. The curves in orange were calculated by CRISPY.

Close modal

The photon-induced LS to HS transition of Fe2+(phen)3 in water was probed by TRWAXS and TRXES using synchrotron radiation with 100 ps resolution. The TRWAXS data, measured over the wide q range 0.5–10.7 Å−1, allowed to probe multiple length scales, from atom–atom distances inside the complex, the solvation structure around it, and the hydrodynamics of the solvent.

The difference scattering ΔS(q, t) is dominated by the change in the solvent temperature from the relaxation of the 1MLCT (< 1 ps) and HS (725 ps) states. After thermal correction of the scattering data, the HS and LS structures were compared with DFT/MD simulations that provided trial functions for the LS and HS states and their cage structures. The high-q range 4.0–10.7 Å−1 is in excellent agreement with the calculated Debye scattering from DFT with a 0.19 Å Fe–N bond elongation in the LS → HS transition. The elongation is caused by the antibonding character in the HS state (t2g)4(eg)2 that has two 3d electrons in the antibonding orbitals dx2-y2 and dz2.

The solvent cage around the LS and HS states were simulated by MD and the Scage(q) term was calculated from the gi,j(r) functions. MD predicts a 1.0 Å contraction of the first solvation shell and a positive cage amplitude in the low-q range 0–0.8 Å−1. The experiment shows that the cage amplitude is overestimated by MD and by reducing it to the 33% level, the fit to the data greatly improves. The amplitude reduction in q space corresponds to a 0.33 Å contraction of the cage radius in the HS state.

The time dependence of the solvent temperature was derived from the thermal solvent differentials from a separate dye/water experiment. The temperature rise from the slow 725 ps HS→LS relaxation is well explained by the energy difference of 0.61 eV between the HS and LS states and the initial HS concentration derived from the Debye fitting. In contrast, the temperature rise from the ultrafast (<1 ps) 1MLCT → HS transition is 2.1× greater than expected from the DFT calculated energy difference of 2.49 eV (3.10–0.61) and the HS concentration. The anomaly is caused by the re-excitation of nascent HS molecules by the 1.2 ps pump pulse at 400 nm. . Ultrafast optical/optical absorption measurements showed the appearance of a new absorption feature after 0.19 ps that is tentatively assigned to 5T2 → 5MLCT. ISC between 5MLCT and 1MLCT then provides a link that restarts the photocycle which explains the excess temperature rise of the solvent.

Finally, the spin states of LS and HS were probed by XES with the core-holes produced by a pink undulator beam. When the 100 ps difference spectrum is compared with reference spectra for putative spin changes, it is assigned to the S = 0 to S = 2 transition. In the HS state, the Kβ1,3 peak from the 3p3/2 → 1s transition has a weak side-peak Kβ′ from the transition 3p1/2 → 1s which is −16 eV from Kβ1,3 . The increase in the Kβ′ intensity is due to the greater exchange interaction in the magnetic HS state.

We would like to thank Professor Carmen Sousa from the University of Barcelona most sincerely for explaining her theoretical work for us in great detail. We are also indebted to Dmitry Khakhulin and Kristoffer Haldrup for sharing their insight and experience with us. We also wish to thank Denis Leshchev, Norman Kretzschmar and Martin Nors Pedersen for help during the early stages of the project. We would also like to thank Laurent Eybert and Wolfgang Reichenbach for their work on the experimental setups on beamline ID09 at ESRF. The use of the femtosecond transient absorption spectrometer in Mathias Bargheer's lab at the University of Potsdam is also gratefully acknowledged. Finally, we are grateful for the EU Grant X-PROBE (Grant 637295) that funded Victoria Kabanova's Ph.D. studentship.

The authors have no conflicts to disclose.

Victoria Kabanova: Investigation (lead). Mathias Sander: Investigation (supporting). Matteo Levantino: Investigation (supporting). Qingyu Kong: Formal analysis (equal). Sophie Canton: Investigation (supporting). Marius Retegan: Formal analysis (supporting). Marco Cammarata: Validation (equal). Philipp Lenzen: Validation (supporting). Latévi Max Daku Lawson: Methodology (supporting). Michael Wulff: Investigation (equal); Validation (equal); Writing – original draft (equal).

The data that support the findings of this study are available within the supplementary material.

1.
N. S.
Lewis
, “
Research opportunities to advance solar energy utilization
,”
Science
351
(
6271
),
aad1920
(
2016
).
2.
A. B.
Wöhri
et al, “
Light-induced structural changes in a photosynthetic reaction center caught by Laue diffraction
,”
Science
328
(
5978
),
630
633
(
2010
).
3.
S. W.
Kohl
et al, “
Consecutive thermal H2 and light-induced O2 evolution from water promoted by a metal complex
,”
Science
324
(
5923
),
74
77
(
2009
).
4.
J. K.
McCusker
, “
Electronic structure in the transition metal block and its implications for light harvesting
,”
Science
363
(
6426
),
484
488
(
2019
).
5.
A.
Vlček
, Jr.
,, “
The life and times of excited states of organometallic and coordination compounds
,”
Coord. Chem. Rev.
200–202
,
933
978
(
2000
).
6.
C.
Creutz
et al, “
Lifetimes, spectra, and quenching of the excited states of polypyridine complexes of iron (II), ruthenium (II), and osmium (II)
,”
J. Am. Chem. Soc.
102
(
4
),
1309
1319
(
1980
).
7.
W.
Gawelda
et al, “
Ultrafast nonadiabatic dynamics of [FeII (bpy)3]2+ in solution
,”
J. Am. Chem. Soc.
129
(
26
),
8199
8206
(
2007
).
8.
P.
Chábera
et al, “
FeII Hexa N-heterocyclic carbene complex with a 528 ps metal-to-ligand charge-transfer excited-state lifetime
,”
J. Phys. Chem. Lett.
2018
,
9
,
459
463
;
[PubMed]
P.
Chábera
et al, “
A low-spin Fe(iii) complex with 100-ps ligand-to-metal charge transfer photoluminescence
,”
Nature
543
,
695
699
(
2017
).
[PubMed]
9.
K. S.
Kjaer
et al, “
Luminescence and reactivity of a charge-transfer excited iron complex with nanosecond lifetime
,”
Science
363
(
6424
),
249
253
(
2019
).
10.
M.
Wulff
et al, “
The realization of sub-nanosecond pump and probe experiments at the ESRF
,”
Faraday Disc.
122
,
13
26
(
2003
).
11.
M.
Levantino
et al, “
Structural dynamics probed by X-ray pulses from synchrotrons and XFELs
,”
C. R. Phys.
22
(
S2
),
75
94
(
2021
).
12.
J.
Tribollet
et al, “
Transient absorption spectroscopy of the iron (II)[Fe (phen)3]2+ complex: Study of the non-radiative relaxation of an isolated iron (II) complex
,”
Chem. Phys. Lett.
513
(
1–3
),
42
47
(
2011
).
13.
G.
Gallé
et al, “
Transient absorption spectroscopy of the [Fe (2 CH3-phen)3]2+ complex: Study of the high spin↔ low spin relaxation of an isolated iron (II) complex
,”
Chem. Phys. Lett.
556
,
82
88
(
2013
).
14.
S.
Nozawa
et al, “
Direct probing of spin state dynamics coupled with electronic and structural modifications by picosecond time-resolved XAFS
,”
J. Am. Chem. Soc.
132
(
1
),
61
63
(
2010
).
15.
H.
Wang
et al, “
Development of picosecond time-resolved X-ray absorption spectroscopy by high-repetition-rate laser pump/X-ray probe at Beijing Synchrotron Radiation Facility
,”
J. Synchrotron Radiat.
24
(
3
),
667
673
(
2017
).
16.
D.
Leshchev
et al, “
Tracking the picosecond deactivation dynamics of a photoexcited iron carbene complex by time-resolved X-ray scattering
,”
Chem. Sci.
9
(
2
),
405
414
(
2018
).
17.
C.
Sousa
et al, “
Theoretical study of the light-induced spin crossover mechanism in [Fe (mtz)6]2+ and [Fe (phen)3]2+
,”
J. Phys. Chem. A
121
(
51
),
9720
9727
(
2017
).
18.
B. E.
Warren
,
X-Ray Diffraction
(
Dover Publication Inc
.,
New York
,
1990
), pp.
116
117
.
19.
S.
Bratos
et al, “
X-ray ‘filming’ of atomic motions in chemical reactions
,”
Chem. Phys.
304
(
3
),
245
251
(
2004
).
20.
G.
Hura
et al, “
A high-quality x-ray scattering experiment on liquid water at ambient conditions
,”
J. Chem. Phys.
113
(
20
),
9140
9148
(
2000
).
21.
F.
Hajdu
, “
Revised parameters of the analytic fits for coherent and incoherent scattered X-ray intensities of the first 36 atoms
,”
Acta Crystallogr., Sect. A
28
(
3
),
250
252
(
1972
).
22.
M.
Wulff
,
S.
Bratos
,
A.
Plech
,
R.
Vuilleumier
,
F.
Mirloup
,
M.
Lorenc
,
Q.
Kong
, and
H.
Ihee
,
J. Chem. Phys.
124
,
34501
34513
(
2006
).
23.
K. S.
Kjaer
et al, “
Introducing a standard method for experimental determination of the solvent response in laser pump, X-ray probe time-resolved wide-angle X-ray scattering experiments on systems in solution
,”
Phys. Chem. Chem. Phys.
15
(
36
),
15003
15016
(
2013
).
24.
M.
Cammarata
et al, “
Impulsive solvent heating probed by picosecond x-ray diffraction
,”
J. Chem. Phys.
124
(
12
),
124504
(
2006
).
25.
J. L. F.
Abascal
and
C.
Vega
, “
A general purpose model for the condensed phases of water: TIP4P/2005
,”
J. Chem. Phys.
123
(
23
),
234505
(
2005
).
26.
L.
Daku
,
L.
Max
, and
A.
Hauser
, “
Ab initio molecular dynamics study of an aqueous solution of [Fe (bpy)3](Cl)2 in the low-spin and in the high-spin states
,”
J. Phys. Chem. Lett.
1
(
12
),
1830
1835
(
2010
).
27.
D.
Khakhulin
et al, “
Visualizing the coordination-spheres of photoexcited transition metal complexes with ultrafast hard X-rays
,”
Phys. Chem. Chem. Phys.
21
(
18
),
9277
9284
(
2019
).
28.
K.
Haldrup
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)32+
,”
J. Phys. Chem. A
116
(
40
),
9878
9887
(
2012
).
29.
M.
Berg
,
A. L.
Harris
, and
C. B.
Harris
, “
Rapid solvent-induced recombination and slow energy relaxation in a simple chemical reaction: Picosecond studies of iodine photodissociation in CCl4
,”
Phys. Rev. Lett.
54
(
9
),
951
(
1985
).
30.
Lawrence Berkeley National Laboratory
,
X-Ray Data Booklet
(
University California
,
2001
).
31.
D.-J.
Huang
,
D. M.
Riffe
, and
J. L.
Erskine
, “
Simultaneous determination of Fe 3p spin-orbit and exchange splittings in photoemission
,”
Phys. Rev. B
51
(
21
),
15170
(
1995
).
32.
W.
Zhang
et al, “
Tracking excited-state charge and spin dynamics in iron coordination complexes
,”
Nature
509
(
7500
),
345
348
(
2014
).
33.
C. J.
Pollock
and
S.
DeBeer
, “
Valence-to-core X-ray emission spectroscopy: A sensitive probe of the nature of a bound ligand
,”
J. Am. Chem. Soc.
133
(
14
),
5594
5601
(
2011
).
34.
G.
Vanko
et al, “
Detailed characterization of a nanosecond-lived excited state: X-ray and theoretical investigation of the quintet state in photoexcited [Fe (terpy)2]2+
,”
J. Phys. Chem. C
119
(
11
),
5888
5902
(
2015
).
35.
See the supplementary material provides more information about the experimental parameters, the atom–atom gij(r) functions from MD and the scaling of the cage term. The change in the scattering from a uniform expansion of the Fe–N bonds is also presented. The calculation highlights the sensitivity of the WAXS data to changes in the Fe–N bond lengths. That is followed by a presentation of the dye based measurements of the thermal scattering differentials for water. Then follows a discussion of alternative explanations for the excess solvent heating, such as ultrafast combinations pathways and direct laser excitation of water. The XES intensity generated by a pink undulator beam is then presented followed by a description of the Von Hámos spectrometer and the fitting of the measured XES spectra to putative reference spins states.