Excited-state structural characterization of a series of nanosecond-lived [Fe(terpy)2]2+ derivatives using x-ray solution scattering

[ F e ( t e r p y ) 2 ] 2 + (terpy = 2,2':6',2″-terpyridine) is a transition metal complex where the spin state is photoswitchable and where the properties of the metal-centered quintet excited state (5MC) can be tuned by substituting different electron withdrawing or electron donating groups on the 4' position of the terpyridine. To better understand the physics determining the photoswitching performance, a deeper insight into the positions of the relevant potential energy surfaces and the molecular structure of the 5MC state is needed. We present a structural investigation based on Time Resolved x-ray Solution Scattering (TR-XSS) by which we determine the average dFe-N bond-length elongation following population of the 5MC state as well as the lifetime of this state in a series of seven modified [Fe(terpy)2]2+ systems in aqueous solution following photo-excitation. The analysis of the TR-XSS data is supported by Density Functional Theory (DFT) and Molecular Dynamics calculations. The quintet state lifetime is determined to vary by more than a factor of 10 (from 1.5 to 16 ns) based on the electron withdrawing/donating properties of the substituting group. Both the DFT calculations and the structural analysis of the experimental data show that the main photo-induced change in metal-ligand bond lengths ΔdFe-N is ∼0.2 Å for all systems.

Capturing ultrafast molecular motions and lattice dynamics in spin crossover film using femtosecond diffraction methods

A comprehensive insight into ultrafast dynamics of photo-switchable materials is desired for efficient control of material properties through light excitation. Here, we study a polycrystalline spin crossover thin film as a prototypical example and reveal the sequential photo-switching dynamics, from local molecular rearrangement to global lattice deformation. On the earliest femtosecond timescale, the local molecular structural rearrangement occurs within a constant unit-cell volume through a two-step process, involving initial Fe-ligand bond elongation followed by ligand rotation. The highly-oriented structure of the nanocrystalline films and the experimental geometry enables resolving the full anisotropic lattice structural dynamics in and out of the sample plane separately. While both molecular switching and lattice heating influence lattice volume, they exert varying degrees of impact at disparate time scales following photoexcitation. This study highlights the opportunities provided by Mega-electron-volt electron and X-ray free electron laser to advance the understanding of ultrafast dynamics of photo-switchable materials.

Structure and spin of the low- and high-spin states of Fe2+(phen)3 studied by x-ray scattering and emission spectroscopy

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.

Water Solvent Reorganization upon Ultrafast Resonant Stimulated X-ray Raman Excitation of a Metalloporphyrin Dimer

We propose an X-ray Raman pump-X-ray diffraction probe scheme to follow solvation dynamics upon charge migration in a solute molecule. The X-ray Raman pump selectively prepares a valence electronic wavepacket in the solute, while the probe provides information about the entire molecular ensemble. A combination of molecular dynamics and ab initio quantum chemistry simulations is applied to a Zn-Ni porphyrin dimer in water. Using time-resolved X-ray diffraction and pair distribution functions, we extracted solvation shell dynamics.

Toward a quantitative description of solvation structure: a framework for differential solution scattering measurements

Appreciating that the role of the solute-solvent and other outer-sphere interactions is essential for understanding chemistry and chemical dynamics in solution, experimental approaches are needed to address the structural consequences of these interactions, complementing condensed-matter simulations and coarse-grained theories. High-energy X-ray scattering (HEXS) combined with pair distribution function analysis presents the opportunity to probe these structures directly and to develop quantitative, atomistic models of molecular systems in situ in the solution phase. However, at concentrations relevant to solution-phase chemistry, the total scattering signal is dominated by the bulk solvent, prompting researchers to adopt a differential approach to eliminate this unwanted background. Though similar approaches are well established in quantitative structural studies of macromolecules in solution by small- and wide-angle X-ray scattering (SAXS/WAXS), analogous studies in the HEXS regime-where sub-ångström spatial resolution is achieved-remain underdeveloped, in part due to the lack of a rigorous theoretical description of the experiment. To address this, herein we develop a framework for differential solution scattering experiments conducted at high energies, which includes concepts of the solvent-excluded volume introduced to describe SAXS/WAXS data, as well as concepts from the time-resolved X-ray scattering community. Our theory is supported by numerical simulations and experiment and paves the way for establishing quantitative methods to determine the atomic structures of small molecules in solution with resolution approaching that of crystallography.

Eliminating finite-size effects on the calculation of x-ray scattering from molecular dynamics simulations

Structural studies using x-ray scattering methods for investigating molecules in solution are shifting focus toward describing the role and effects of the surrounding solvent. However, forward models based on molecular dynamics (MD) simulations to simulate structure factors and x-ray scattering from interatomic distributions such as radial distribution functions (RDFs) face limitations imposed by simulations, particularly at low values of the scattering vector q. In this work, we show how the value of the structure factor at q = 0 calculated from RDFs sampled from finite MD simulations is effectively dependent on the size of the simulation cell. To eliminate this error, we derive a new scheme to renormalize the sampled RDFs based on a model of the excluded volume of the particle-pairs they were sampled from, to emulate sampling from an infinite system. We compare this new correction method to two previous RDF-correction methods, developed for Kirkwood-Buff theory applications. We present a quantitative test to assess the reliability of the simulated low-q scattering signal and show that our RDF-correction successfully recovers the correct q = 0 limit for neat water. We investigate the effect of MD-sampling time on the RDF-corrections, before advancing to a molecular example system, comprised of a transition metal complex solvated in a series of water cells with varying densities. We show that our correction recovers the correct q = 0 behavior for all densities. Furthermore, we employ a simple continuum scattering model to dissect the total scattering signal from the solvent-solvent structural correlations in a solute-solvent model system to find two distinct contributions: a non-local density-contribution from the finite, fixed cell size in NVT simulations, and a local contribution from the solvent shell. We show how the second contribution can be approximated without also including the finite-size contribution. Finally, we provide a "best-practices"-checklist for experimentalists planning to incorporate explicit solvation MD simulations in future work, offering guidance for improving the accuracy and reliability of structural studies using x-ray scattering methods in solution.

Exploring fingerprints of ultrafast structural dynamics in molecular solutions with an X-ray laser

We apply ultrashort X-ray laser pulses to track optically excited structural dynamics of [Ir2(dimen)4]2+ molecules in solution. In our exploratory study we determine angular correlations in the scattered X-rays, which comprise a complex fingerprint of the ultrafast dynamics. Model-assisted analysis of the experimental correlation data allows us to elucidate various aspects of the photoinduced changes in the excited molecular ensembles. We unambiguously identify that in our experiment the photoinduced transition dipole moments in [Ir2(dimen)4]2+ molecules are oriented perpendicular to the Ir-Ir bond. The analysis also shows that the ground state conformer of [Ir2(dimen)4]2+ with a larger Ir-Ir distance is mostly responsible for the formation of the excited state. We also reveal that the ensemble of solute molecules can be characterized with a substantial structural heterogeneity due to solvent influence. The proposed X-ray correlation approach offers an alternative path for studies of ultrafast structural dynamics of molecular ensembles in the liquid and gas phases.

Projection to extract the perpendicular component (PEPC) method for extracting kinetics from time-resolved data

Time-resolved x-ray liquidography (TRXL) is a potent method for investigating the structural dynamics of chemical and biological reactions in the liquid phase. It has enabled the extraction of detailed structural aspects of various dynamic processes, the molecular structures of intermediates, and kinetics of reactions across a wide range of systems, from small molecules to proteins and nanoparticles. Proper data analysis is key to extracting the information of the kinetics and structural dynamics of the studied system encrypted in the TRXL data. In typical TRXL data, the signals from solute scattering, solvent scattering, and solute-solvent cross scattering are mixed in the q-space, and the solute kinetics and solvent hydrodynamics are mixed in the time domain, thus complicating the data analysis. Various methods developed so far generally require prior knowledge of the molecular structures of candidate species involved in the reaction. Because such information is often unavailable, a typical data analysis often involves tedious trial and error. To remedy this situation, we have developed a method named projection to extract the perpendicular component (PEPC), capable of removing the contribution of solvent kinetics from TRXL data. The resulting data then contain only the solute kinetics, and, thus, the solute kinetics can be easily determined. Once the solute kinetics is determined, the subsequent data analysis to extract the structural information can be performed with drastically improved convenience. The application of the PEPC method is demonstrated with TRXL data from the photochemistry of two molecular systems: [Au(CN)₂^-]₃ in water and CHI₃ in cyclohexane.

Simulating the solvation structure of low- and high-spin [Fe(bpy)3]2+: long-range dispersion and many-body effects

When characterizing transition metal complexes and their functionalities, the importance of including the solvent as an active participant is becoming more and more apparent. Whereas many studies have evaluated long-range dispersion effects inside organic molecules and organometallics, less is known about their role in solvation. Here, we have analysed the components within solute-solvent and solvent-solvent interactions of one of the most studied iron-based photoswitch model systems, in two spin states. We find that long-range dispersion effects modulate the coordination significantly, and that this is accurately captured by density functional theory models including dispersion corrections. We furthermore correlate gas-phase relaxed complex-water clusters to thermally averaged molecular densities. This shows how the gas-phase interactions translate to solution structure, quantified through 3D molecular densities, angular distributions, and radial distribution functions. We show that finite-size simulation cells can cause the radial distribution functions to have artificially enlarged amplitudes. Finally, we quantify the effects of many-body interactions within the solvent shells, and find that almost a fifth of the total interaction energy of the solute-shell system in the high-spin state comes from many-body contributions, which cannot be captured by by pair-wise additive force field methods.

Capturing photochemical and photophysical transformations in iron complexes with ultrafast X-ray spectroscopy and scattering

Light-driven chemical transformations provide a compelling approach to understanding chemical reactivity with the potential to use this understanding to advance solar energy and catalysis applications. Capturing the non-equilibrium trajectories of electronic excited states with precision, particularly for transition metal complexes, would provide a foundation for advancing both of these objectives. Of particular importance for 3d metal compounds is characterizing the population dynamics of charge-transfer (CT) and metal-centered (MC) electronic excited states and understanding how the inner coordination sphere structural dynamics mediate the interaction between these states. Recent advances in ultrafast X-ray laser science has enabled the electronic excited state dynamics in 3d metal complexes to be followed with unprecedented detail. This review will focus on simultaneous X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS) studies of iron coordination and organometallic complexes. These simultaneous XES-XSS studies have provided detailed insight into the mechanism of light-induced spin crossover in iron coordination compounds, the interaction of CT and MC excited states in iron carbene photosensitizers, and the mechanism of Fe-S bond dissociation in cytochrome c.

Spin-state dependence of the structural and vibrational properties of solvated iron(ii) polypyridyl complexes from AIMD simulations: III. [Fe(tpen)]Cl2 in acetonitrile

In order to achieve an in-depth understanding of the role played by the solvent in the photoinduced low-spin (LS) → high-spin (HS) transition in solvated Fe(ii) complexes, an accurate description of the solvated complexes in the two spin states is required. To this end, we are applying state-of-the-art ab initio molecular dynamics (AIMD) simulations to the study of the structural and vibrational properties of iron(ii) polypyridyl complexes. Two aqueous LS complexes were investigated in this framework, namely, [Fe(bpy)3]2+ (bpy = 2,2'-bipyridine) [Lawson Daku and Hauser, J. Phys. Chem. Lett., 2010, 1, 1830; Lawson Daku, Phys. Chem. Chem. Phys., 2018, 20, 6236] and [Fe(tpy)2]2+ (tpy = 2,2':6',2''-ter-pyridine) [Lawson Daku, Phys. Chem. Chem. Phys., 2019, 21, 650]. For aqueous [Fe(bpy)3]2+, combining the results of forefront wide-angle X-ray scattering experiments with those of the AIMD simulations allowed the visualization of the interlaced coordination and solvation spheres of the photoinduced HS state [Khakhulin et al., Phys. Chem. Chem. Phys., 2019, 21, 9277]. In this paper, we report the extension of our AIMD studies to the spin-crossover complex [Fe(tpen)]2+ (tpen = N,N,N',N'-tetrakis(2-pyridylmethyl)ethylenediamine) in acetonitrile (ACN). The determined LS and HS solution structures of the complex are in excellent agreement with the experimental results obtained by high-resolution transient X-ray absorption spectroscopy [Zhang et al., ACS Omega, 2019, 4, 6375]. The first solvation shell of [Fe(tpen)]2+ consists of ACN molecules located in the grooves defined by the chelating coordination motif of the tpen ligand. Upon the LS → HS change of states, the solvation number of the complex is found to increase from ≈9.2 to ≈11.9 and an inner solvation shell is formed. This inner solvation shell originates from the occupancy by about one ACN molecule of the internal cavity which results from the arrangement of the 4 pyridine rings of the tpen ligand, and which becomes accessible to the solvent molecules in the HS state only thanks to the structural changes undergone by the complex. The presence of this inner solvation shell for the solvated HS complex probably plays a key role in the spin-state dependent reactivity of [Fe(tpen)]2+ in liquid solutions.

Characterizing electronic and atomic structures for amorphous and molecular metal oxide catalysts at functional interfaces by combining soft X-ray spectroscopy and high-energy X-ray scattering

Amorphous thin film materials and heterogenized molecular catalysts supported on electrode and other functional interfaces are widely investigated as promising catalyst formats for applications in solar and electrochemical fuels catalysis. However the amorphous character of these catalysts and the complexity of the interfacial architectures that merge charge transport properties of electrode and semiconductor supports with discrete sites for multi-step catalysis poses challenges for probing mechanisms that activate and tune sites for catalysis. This minireview discusses advances in soft X-ray spectroscopy and high-energy X-ray scattering that provide opportunities to resolve interfacial electronic and atomic structures, respectively, that are linked to catalysis. This review discusses how these techniques can be partnered with advances in nanostructured interface synthesis for combined soft X-ray spectroscopy and high-energy X-ray scattering analyses of thin film and heterogenized molecular catalysts. These combined approaches enable opportunities for the characterization of both electronic and atomic structures underlying fundamental catalytic function, and that can be applied under conditions relevant to device applications.

Exploring the light-induced dynamics in solvated metallogrid complexes with femtosecond pulses across the electromagnetic spectrum

Oligonuclear complexes of d⁴-d⁷ transition metal ion centers that undergo spin-switching have long been developed for their practical role in molecular electronics. Recently, they also have appeared as promising photochemical reactants demonstrating improved stability. However, the lack of knowledge about their photophysical properties in the solution phase compared to mononuclear complexes is currently hampering their inclusion into advanced light-driven reactions. In the present study, the ultrafast photoinduced dynamics in a solvated [2 × 2] iron(II) metallogrid complex are characterized by combining measurements with transient optical-infrared absorption and x-ray emission spectroscopy on the femtosecond time scale. The analysis is supported by density functional theory calculations. The photocycle can be described in terms of intra-site transitions, where the Fe^II centers in the low-spin state are independently photoexcited. The Franck-Condon state decays via the formation of a vibrationally hot high-spin (HS) state that displays coherent behavior within a few picoseconds and thermalizes within tens of picoseconds to yield a metastable HS state living for several hundreds of nanoseconds. Systematic comparison with the closely related mononuclear complex [Fe(terpy)₂]²⁺ reveals that nuclearity has a profound impact on the photoinduced dynamics. More generally, this work provides guidelines for expanding the integration of oligonuclear complexes into new photoconversion schemes that may be triggered by ultrafast spin-switching.

Direct observation of nuclear reorganization driven by ultrafast spin transitions

One of the most basic molecular photophysical processes is that of spin transitions and intersystem crossing between excited states surfaces. The change in spin states affects the spatial distribution of electron density through the spin orbit coupling interaction. The subsequent nuclear reorganization reports on the full extent of the spin induced change in electron distribution, which can be treated similarly to intramolecular charge transfer with effective reaction coordinates depicting the spin transition. Here, single-crystal [Fe^II(bpy)₃](PF₆)₂, a prototypical system for spin crossover (SCO) dynamics, is studied using ultrafast electron diffraction in the single-photon excitation regime. The photoinduced SCO dynamics are resolved, revealing two distinct processes with a (450 ± 20)-fs fast component and a (2.4 ± 0.4)-ps slow component. Using principal component analysis, we uncover the key structural modes, ultrafast Fe–N bond elongations coupled with ligand motions, that define the effective reaction coordinate to fully capture the relevant molecular reorganization.

Electron spin is a fundamental property of molecules, and changes in spin state affect both molecular structure and dynamics. Here, the authors resolve, by ultrafast electron diffraction, the nuclear reorganization stabilizing spin transitions in a [Fe^II(bpy)₃](PF₆)₂ crystal.

Ultrafast X-ray Photochemistry at European XFEL: Capabilities of the Femtosecond X-ray Experiments (FXE) Instrument

Time-resolved X-ray methods are widely used for monitoring transient intermediates over the course of photochemical reactions. Ultrafast X-ray absorption and emission spectroscopies as well as elastic X-ray scattering deliver detailed electronic and structural information on chemical dynamics in the solution phase. In this work, we describe the opportunities at the Femtosecond X-ray Experiments (FXE) instrument of European XFEL. Guided by the idea of combining spectroscopic and scattering techniques in one experiment, the FXE instrument has completed the initial commissioning phase for most of its components and performed first successful experiments within the baseline capabilities. This is demonstrated by its currently 115 fs (FWHM) temporal resolution to acquire ultrafast X-ray emission spectra by simultaneously recording iron Kα and Kβ lines, next to wide angle X-ray scattering patterns on a photoexcited aqueous solution of [Fe(bpy)3]2+, a transition metal model compound.

Vibrational wavepacket dynamics in Fe carbene photosensitizer determined with femtosecond X-ray emission and scattering

The non-equilibrium dynamics of electrons and nuclei govern the function of photoactive materials. Disentangling these dynamics remains a critical goal for understanding photoactive materials. Here we investigate the photoinduced dynamics of the [Fe(bmip)2]2+ photosensitizer, where bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine, with simultaneous femtosecond-resolution Fe Kα and Kβ X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS). This measurement shows temporal oscillations in the XES and XSS difference signals with the same 278 fs period oscillation. These oscillations originate from an Fe-ligand stretching vibrational wavepacket on a triplet metal-centered (3MC) excited state surface. This 3MC state is populated with a 110 fs time constant by 40% of the excited molecules while the rest relax to a 3MLCT excited state. The sensitivity of the Kα XES to molecular structure results from a 0.7% average Fe-ligand bond length shift between the 1 s and 2p core-ionized states surfaces.