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.

Polarizable Embedding without Artificial Boundary Polarization.

We present a fully self-consistent polarizable embedding (PE) model that does not suffer from unphysical boundary polarization. This is achieved through the use of the minimum-image convention (MIC) in the induced electrostatics. It is a simple yet effective approach that includes a more physically accurate description of the polarization throughout the molecular system. Using PE with MIC (PE-MIC), we shed new light on the limitations of commonly employed cutoff models, such as the droplet model, when used in PE calculations. Specifically, we investigate the effects of the unphysical polarization at the outer boundary by comparing induced dipoles and the associated electrostatic potentials, as well as some optical properties of solute-solvent and biomolecular systems. We show that the magnitude of the inaccuracies caused by the unphysical polarization depends on multiple parameters: the nature of the quantum subsystem and of the environment, the cutoff model and distance, and the calculated property.

Ultrafast Jahn-Teller Photoswitching in Cobalt Single-Ion Magnets.

Single-ion magnets (SIMs) constitute the ultimate size limit in the quest for miniaturizing magnetic materials. Several bottlenecks currently hindering breakthroughs in quantum information and communication technologies could be alleviated by new generations of SIMs displaying multifunctionality. Here, ultrafast optical absorption spectroscopy and X-ray emission spectroscopy are employed to track the photoinduced spin-state switching of the prototypical complex [Co(terpy)2 ]2+ (terpy = 2,2':6',2″-terpyridine) in solution phase. The combined measurements and their analysis supported by density functional theory (DFT), time-dependent-DFT (TD-DFT) and multireference quantum chemistry calculations reveal that the complex undergoes a spin-state transition from a tetragonally elongated doublet state to a tetragonally compressed quartet state on the femtosecond timescale, i.e., it sustains ultrafast Jahn-Teller (JT) photoswitching between two different spin multiplicities. Adding new Co-based complexes as possible contenders in the search for JT photoswitching SIMs will greatly widen the possibilities for implementing magnetic multifunctionality and eventually controlling ultrafast magnetization with optical photons.

Tracking structural solvent reorganization and recombination dynamics following e- photoabstraction from aqueous I- with femtosecond x-ray spectroscopy and scattering.

We present a sub-picosecond resolved investigation of the structural solvent reorganization and geminate recombination dynamics following 400 nm two-photon excitation and photodetachment of a valence p electron from the aqueous atomic solute, I-(aq). The measurements utilized time-resolved X-ray Absorption Near Edge Structure (TR-XANES) spectroscopy and X-ray Solution Scattering (TR-XSS) at the Linac Coherent Light Source x-ray free electron laser in a laser pump/x-ray probe experiment. The XANES measurements around the L1-edge of the generated nascent iodine atoms (I0) yield an average electron ejection distance from the iodine parent of 7.4 ± 1.5 Å with an excitation yield of about 1/3 of the 0.1M NaI aqueous solution. The kinetic traces of the XANES measurement are in agreement with a purely diffusion-driven geminate iodine-electron recombination model without the need for a long-lived (I0:e-) contact pair. Nonequilibrium classical molecular dynamics simulations indicate a delayed response of the caging H2O solvent shell and this is supported by the structural analysis of the XSS data: We identify a two-step process exhibiting a 0.1 ps delayed solvent shell reorganization time within the tight H-bond network and a 0.3 ps time constant for the mean iodine-oxygen distance changes. The results indicate that most of the reorganization can be explained classically by a transition from a hydrophilic cavity with a well-ordered first solvation shell (hydrogens pointing toward I-) to an expanded cavity around I0 with a more random orientation of the H2O molecules in a broadened first solvation shell.

Disentangling the resonant Auger spectra of ozone: overlapping core-hole states and core-excited state dynamics.

We investigate the resonant and non-resonant Auger spectra of ozone with a newly implemented multi-reference protocol based on the one-center approximation [Tenorio et al., J. Chem. Theory Comput. 2022, 18, 4387-4407]. The results of our calculations are compared to existing experimental data, where we elucidate the resonant Auger spectrum measured at 530.8 and 536.7 eV, that correspond to the 1sOT → π*(2b1) and 1sOT → σ*(7a1) resonances, and at 542.3 eV, which lies near the 1sOC → σ*(7a1) excited state and above the 1sOT-1 ionization threshold. Using molecular dynamics simulations, we demonstrate the relevance of few-femtoseconds nuclear dynamics in the resonant Auger spectrum of ozone following the 1sOT → π*(2b1) core-excitation.

Resolving Femtosecond Solvent Reorganization Dynamics in an Iron Complex by Nonadiabatic Dynamics Simulations.

The ultrafast dynamical response of solute-solvent interactions plays a key role in transition metal complexes, where charge transfer states are ubiquitous. Nonetheless, there exist very few excited-state simulations of transition metal complexes in solution. Here, we carry out a nonadiabatic dynamics study of the iron complex [Fe(CN)4(bpy)]2- (bpy = 2,2'-bipyridine) in explicit aqueous solution. Implicit solvation models were found inadequate for reproducing the strong solvatochromism in the absorption spectra. Instead, direct solute-solvent interactions, in the form of hydrogen bonds, are responsible for the large observed solvatochromic shift and the general dynamical behavior of the complex in water. The simulations reveal an overall intersystem crossing time scale of 0.21 ± 0.01 ps and a strong reliance of this process on nuclear motion. A charge transfer character analysis shows a branched decay mechanism from the initially excited singlet metal-to-ligand charge transfer (1MLCT) states to triplet states of 3MLCT and metal-centered (3MC) character. We also find that solvent reorganization after excitation is ultrafast, on the order of 50 fs around the cyanides and slower around the bpy ligand. In contrast, the nuclear vibrational dynamics, in the form of Fe-ligand bond changes, takes place on slightly longer time scales. We demonstrate that the surprisingly fast solvent reorganizing should be observable in time-resolved X-ray solution scattering experiments, as simulated signals show strong contributions from the solute-solvent scattering cross term. Altogether, the simulations paint a comprehensive picture of the coupled and concurrent electronic, nuclear, and solvent dynamics and interactions in the first hundreds of femtoseconds after excitation.

Femtochemistry of bimolecular reactions from weakly bound complexes: computational study of the H + H'OD → H'OH + D or HOD + H' exchange reactions.

A full-dimensional wavepacket propagation describing the bimolecular exchange reactions H + H'OD → H'OH + D or HOD + H' initiated by photolysis of HCl in the hydrogen-bound complex (HCl)⋯(HOD) is reported. The dynamics of this reaction is carried out with the MCTDH method on an ab initio potential energy surface (PES) of H3O and the initial state is derived from the ground state wavefunction of the complex obtained by relaxation on its own electronic ground state ab initio PES. The description of the system makes use of polyspherical coordinates parametrizing a set of Radau and Jacobi vectors. The calculated energy- and time-resolved reaction probabilities show, owing to the large collision energies at play stemming from the (almost full) photolysis of HCl, that the repulsion between oxygen in the H'OD molecule and the incoming hydrogen atom is the main feature of the collision and leads to non-reactive scattering. No abstraction reaction products are observed. However, both exchange processes are still observable, with a preference in O-H' bond dissociation over that of O-D. The selectivity is reversed upon vibrational pre-excitation of the O-D stretching mode in the H'OD molecule. It is shown that, after the collision, the hydrogen atom of HCl does most likely not encounter the almost stationary chlorine atom again but we also consider the limit case where the H atom is forced to collide multiple times against H'OD as a result of being pushed back by the Cl atom.

Ultrafast Rotational and Translational Energy Relaxation in Neat Liquids.

The excess energy flow pathways during rotational and translational relaxation induced by rotational or translational excitation of a single molecule of and within each of four different neat liquids (H2O, MeOH, CCl4, and CH4) are studied using classical molecular dynamics simulations and energy flux analysis. For all four liquids, the relaxation processes for both types of excitation are ultrafast, but the energy flow is significantly faster for the polar, hydrogen-bonded (H-bonded) liquids H2O and MeOH. Whereas the majority of the initial excess energy is transferred into hindered rotations (librations) for rotational excitation in the H-bonded liquids, an almost equal efficiency for transfer to translational and rotational motions is observed in the nonpolar, non-H-bonded liquids CCl4 and CH4. For translational excitation, transfer to translational motions dominates for all liquids. In general, the energy flows are quite local; i.e., more than 70% of the energy flows directly to the first solvent shell molecules, reaching almost 100% for CCl4 and CH4. Finally, the determined validity of linear response theory for these nonequilibrium relaxation processes is quite solvent-dependent, with the deviation from linear response most marked for rotational excitation and for the nonpolar liquids.

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation.

Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized 1B2u (ππ*) (S2) and 1B3u (nπ*) (S1) states, the participation of the optically dark 1Au (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite 1Au (nπ*) and 1B3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The 1Au (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

Trajectory surface-hopping photoinduced dynamics from Rydberg states of trimethylamine.

We present a computational study on nonadiabatic excited-state dynamics initiated from the 3p Rydberg states of trimethylamine (TMA). We utilise a methodology based on full-dimensional (39 D) trajectory surface-hopping (TSH) simulations, in which propagation is carried out on on-the-fly density functional theory (DFT)/time-dependent DFT (TD-DFT) potentials. Both our electronic structure benchmarks to high-level ab initio methods (EOM-CCSD, CASPT2) and TSH simulations demonstrate high-accuracy of the applied CAM-B3LYP functional for the description of Rydberg excited states. Based on our excited-state simulations, we construct the following mechanistic picture: when pumped resonantly to the 3p Rydberg manifold, TMA coherently vibrates along the planarisation mode with a period of 104 fs and an exponential coherence decay time constant of 240 fs. Nonadiabatic dynamics occur on a faster (∼1 ps) and a slower (∼3 ps) timescale, along the N-C stretching mode by mixing with a dissociative σN-C* state. As a minor relaxation channel, 3p → 3s internal conversion occurs via branching at the σN-C*/3s intersection. We find that photodissociaton is hardly observable within 3 ps (1%), which is a failure of the range-separated hybrid CAM-B3LYP functional, as a consequence of its static electron correlation deficiency at long range. In contrast, pure DFT (GGA-BLYP) provides an accurate long-range description (19% dissociation yield), also supported by comparison to recent ultrafast experiments, even if the Rydberg state energies are significantly underestimated (>1 eV). Finally, we reveal the crucial role of vibrational coherence and energy transfer from the planarisation mode for N-C bond activation and resulting nonadiabatic dynamics. The present work illustrates the importance of nuclear-electronic coupling for excited-state dynamics and branching at conical intersections.

An assessment of different electronic structure approaches for modeling time-resolved x-ray absorption spectroscopy.

We assess the performance of different protocols for simulating excited-state x-ray absorption spectra. We consider three different protocols based on equation-of-motion coupled-cluster singles and doubles, two of them combined with the maximum overlap method. The three protocols differ in the choice of a reference configuration used to compute target states. Maximum-overlap-method time-dependent density functional theory is also considered. The performance of the different approaches is illustrated using uracil, thymine, and acetylacetone as benchmark systems. The results provide guidance for selecting an electronic structure method for modeling time-resolved x-ray absorption spectroscopy.

Publisher's Note: Femtosecond X-Ray Scattering Study of Ultrafast Photoinduced Structural Dynamics in Solvated [Co(terpy)_{2}]

This corrects the article DOI: 10.1103/PhysRevLett.117.013002.

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.

Hot Branching Dynamics in a Light-Harvesting Iron Carbene Complex Revealed by Ultrafast X-ray Emission Spectroscopy.

Iron N-heterocyclic carbene (NHC) complexes have received a great deal of attention recently because of their growing potential as light sensitizers or photocatalysts. We present a sub-ps X-ray spectroscopy study of an FeII NHC complex that identifies and quantifies the states involved in the deactivation cascade after light absorption. Excited molecules relax back to the ground state along two pathways: After population of a hot 3 MLCT state, from the initially excited 1 MLCT state, 30 % of the molecules undergo ultrafast (150 fs) relaxation to the 3 MC state, in competition with vibrational relaxation and cooling to the relaxed 3 MLCT state. The relaxed 3 MLCT state then decays much more slowly (7.6 ps) to the 3 MC state. The 3 MC state is rapidly (2.2 ps) deactivated to the ground state. The 5 MC state is not involved in the deactivation pathway. The ultrafast partial deactivation of the 3 MLCT state constitutes a loss channel from the point of view of photochemical efficiency and highlights the necessity to screen transition-metal complexes for similar ultrafast decays to optimize photochemical performance.

Time-resolved near-edge X-ray absorption fine structure of pyrazine from electronic structure and nuclear wave packet dynamics simulations.

As a demonstration of the analysis of the electronic structure and the nuclear dynamics from time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS), we present the TR-NEXAFS spectra of pyrazine following the excitation to the 1B2u(ππ*) state. The spectra are calculated combining the frozen-core/core-valence separated equation-of-motion coupled cluster singles and doubles approach for the spectral signatures and the multiconfiguration time-dependent Hartree method for the wave packet propagation. The population decay from the 1B2u(ππ*) state to the 1B3u(nπ*) and 1Au(nπ*) states, followed by oscillatory flow of population between the 1B3u(nπ*) and 1Au(nπ*) states, is interpreted by means of visualization of the potential energy curves and the reduced nuclear densities. By examining the electronic structure of the three valence-excited states and the final core-excited states, we observe that the population dynamics is explicitly reflected in the TR-NEXAFS spectra, especially when the heteroatoms are selected as the X-ray absorption sites. This work illustrates the feasibility of extracting fine details of molecular photophysical processes from TR-NEXAFS spectra by using currently available theoretical methods.

Simulation of ultrafast excited-state dynamics and elastic x-ray scattering by quantum wavepacket dynamics.

Simulation of the ultrafast excited-state dynamics and elastic X-ray scattering of the [Fe(bmip)2]2+ [bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-4-pyridine] complex is presented and analyzed. We employ quantum wavepacket dynamics simulations on a 5-dimensional potential energy surface (PES) calculated by time-dependent density functional theory with 26 coupled diabatic states. The simulations are initiated by explicit inclusion of a time-dependent electromagnetic field. In the case of resonant excitation into singlet metal-to-ligand charge transfer (1MLCT) states, kinetic (exponential) population dynamics are observed with small nuclear motion. In agreement with transient optical absorption spectroscopy experiments, we observe a subpicosecond 1MLCT → 3MLCT intersystem crossing and a subsequent decay into triplet metal-centered (3MC) states on a picosecond time scale. The simulated time-resolved difference scattering signal is dominated by the 3MC component, for which the structural distortions are significant. On the other hand, excitation into 1MC states leads to ballistic (nonexponential) population dynamics with strong nuclear motion. The reason for these ballistic dynamics is that in this case, the excitation occurs into a nonequilibrium region, i.e., far from the minimum of the 1MC PES. This results in wavepacket dynamics along the principal breathing mode, which is clearly visible in both the population dynamics and difference scattering. Finally, the importance of decomposing the difference scattering into components by electronic states is highlighted, information which is not accessible from elastic X-ray scattering experiments.

Finding intersections between electronic excited state potential energy surfaces with simultaneous ultrafast X-ray scattering and spectroscopy.

Light-driven molecular reactions are dictated by the excited state potential energy landscape, depending critically on the location of conical intersections and intersystem crossing points between potential surfaces where non-adiabatic effects govern transition probabilities between distinct electronic states. While ultrafast studies have provided significant insight into electronic excited state reaction dynamics, experimental approaches for identifying and characterizing intersections and seams between electronic states remain highly system dependent. Here we show that for 3d transition metal systems simultaneously recorded X-ray diffuse scattering and X-ray emission spectroscopy at sub-70 femtosecond time-resolution provide a solid experimental foundation for determining the mechanistic details of excited state reactions. In modeling the mechanistic information retrieved from such experiments, it becomes possible to identify the dominant trajectory followed during the excited state cascade and to determine the relevant loci of intersections between states. We illustrate our approach by explicitly mapping parts of the potential energy landscape dictating the light driven low-to-high spin-state transition (spin crossover) of [Fe(2,2'-bipyridine)3]2+, where the strongly coupled nuclear and electronic dynamics have been a source of interest and controversy. We anticipate that simultaneous X-ray diffuse scattering and X-ray emission spectroscopy will provide a valuable approach for mapping the reactive trajectories of light-triggered molecular systems involving 3d transition metals.

Mechanism of Photoinduced Dihydroazulene Ring-Opening Reaction.

The photoinduced ring-opening reaction is a key process in the functioning of dihydroazulene/vinylheptafulvene (DHA/VHF) photoswitches. Over the years, the mechanism of this reaction has been extensively debated. Herein, by means of nonadiabatic trajectory dynamics simulations and quantum chemistry calculations, we present the first detailed and comprehensive investigation on the mechanism of the photoinduced ring-opening reaction of DHA. The results show the crucial role of the excited-state ring planarization process for the bond breaking. Our dynamics simulations show that the DHA ring opening is an ultrafast reaction that does not follow exponential kinetics but exhibits ballistic dynamics. Upon photoexcitation, the planarization occurs within 300-500 fs. This leads to the ring-opening reaction and concurrent decay of the molecule to the ground state within 100 fs through an S1 → S0 internal conversion process toward forming the VHF isomer. These results are consistent with previous ultrafast time-resolved experiments and lead to a thorough understanding of the DHA/VHF photoconversion.

Ultrafast structural dynamics of photo-reactions observed by time-resolved x-ray cross-correlation analysis.

We applied angular X-ray Cross-Correlation analysis (XCCA) to scattering images from a femtosecond resolution X-ray free-electron laser pump-probe experiment with solvated PtPOP {[Pt2(P2O5H2)4]4-} metal complex molecules. The molecules were pumped with linear polarized laser pulses creating an excited state population with a preferred orientational (alignment) direction. Two time scales of 1.9 ± 1.5 ps and 46 ± 10 ps were revealed by angular XCCA associated with structural changes and rotational dephasing of the solvent molecules, respectively. These results illustrate the potential of XCCA to reveal hidden structural information in the analysis of time-resolved x-ray scattering data from molecules in solution.