Ab initiotreatment of time-resolved x-ray scattering: Application to the photoisomerization of stilbene

Robust Inversion of Time-Resolved Data via Forward-Optimization in a Trajectory Basis

An inversion method for time-resolved data from ultrafast experiments is introduced, based on forward-optimization in a trajectory basis. The method is applied to experimental data from X-ray scattering of the photochemical ring-opening reaction of 1,3-cyclohexadiene and electron diffraction of the photodissociation of CS₂. In each case, inversion yields a model that reproduces the experimental data, identifies the main dynamic motifs, and agrees with independent experimental observations. Notably, the method explicitly accounts for continuity constraints and is robust even for noisy data.

Towards high-resolution X-ray scattering as a probe of electron correlation

We demonstrate that X-ray scattering can be used as a probe of electron–electron correlation.

Mapping static core-holes and ring-currents with X-ray scattering

Measuring the attosecond movement of electrons in molecules is challenging due to the high temporal and spatial resolutions required. X-ray scattering-based methods are promising, but many questions remain concerning the sensitivity of the scattering signals to changes in density, as well as the means of reconstructing the dynamics from these signals. In this paper, we present simulations of stationary core-holes and electron dynamics following inner-shell ionization of the oxazole molecule. Using a combination of time-dependent density functional theory simulations along with X-ray scattering theory, we demonstrate that the sudden core-hole ionization produces a significant change in the X-ray scattering response and how the electron currents across the molecule should manifest as measurable modulations to the time dependent X-ray scattering signal. This suggests that X-ray scattering is a viable probe for measuring electronic processes at time scales faster than nuclear motion.

Femtosecond X-ray Liquidography Visualizes Wavepacket Trajectories in Multidimensional Nuclear Coordinates for a Bimolecular Reaction

ConspectusVibrational wavepacket motions on potential energy surfaces are one of the critical factors that determine the reaction dynamics of photoinduced reactions. The motions of vibrational wavepackets are often discussed in the interpretation of observables measured with various time-resolved vibrational or electronic spectroscopies but mostly in terms of the frequencies of wavepacket motions, which are approximated by normal modes, rather than the actual positions of the wavepacket. Although the time-dependent positions (that is, the trajectory) of wavepackets are hypothesized or drawn in imagined or calculated potential energy surfaces, it is not trivial to experimentally determine the trajectory of wavepackets, especially in multidimensional nuclear coordinates for a polyatomic molecule. Recently, we performed a femtosecond X-ray liquidography (solution scattering) experiment on a gold trimer complex (GTC), [Au(CN)₂^-]₃, in water at X-ray free-electron lasers (XFELs) and elucidated the time-dependent positions of vibrational wavepackets from the Franck-Condon region to equilibrium structures on both excited and ground states in the course of the formation of covalent bonds between gold atoms.Bond making is an essential process in chemical reactions, but it is challenging to keep track of detailed atomic movements associated with bond making because of its bimolecular nature that requires slow diffusion of two reaction parties to meet each other. Bond formation in the solution phase has been elusive because the diffusion of the reactants limits the reaction rate of a bimolecular process, making it difficult to initiate and track the bond-making processes with an ultrafast time resolution. In principle, if the bimolecular encounter can be controlled to overcome the limitation caused by diffusion, the bond-making processes can be tracked in a time-resolved manner, providing valuable insight into the bimolecular reaction mechanism. In this regard, GTC offers a good model system for studying the dynamics of bond formation in solution. Au(I) atoms in GTC exhibit a noncovalent aurophilic interaction, making GTC an aggregate complex without any covalent bond. Upon photoexcitation of GTC, an electron is excited from an antibonding orbital to a bonding orbital, leading to the formation of covalent bonds among Au atoms. Since Au atoms in the ground state of GTC are located in close proximity within the same solvent cage, the formation of Au-Au covalent bonds occurs without its reaction rate being limited by diffusion through the solvent.Femtosecond time-resolved X-ray liquidography (fs-TRXL) data revealed that the ground state has an asymmetric bent structure. From the wavepacket trajectory determined in three-dimensional nuclear coordinates (two internuclear distances and one bond angle), we found that two covalent bonds are formed between three Au atoms of GTC asynchronously. Specifically, one covalent bond is formed first for the shorter Au-Au pair (of the asymmetric and bent ground-state structure) in 35 fs, and subsequently, the other covalent bond is formed for the longer Au-Au pair within 360 fs. The resultant trimer complex has a symmetric and linear geometry, implying the occurrence of bent-to-linear transformation concomitant with the formation of two equivalent covalent bonds, and exhibits vibrations that can be unambiguously assigned to specific normal modes based on the wavepacket trajectory, even without the vibrational frequencies provided by quantum calculation.

Mapping the emergence of molecular vibrations mediating bond formation

Fundamental studies of chemical reactions often consider the molecular dynamics along a reaction coordinate using a calculated or suggested potential energy surface^1-5. But fully mapping such dynamics experimentally, by following all nuclear motions in a time-resolved manner-that is, the motions of wavepackets-is challenging and has not yet been realized even for the simple stereotypical bimolecular reaction^6-8: A-B + C → A + B-C. Here we track the trajectories of these vibrational wavepackets during photoinduced bond formation of the gold trimer complex [Au(CN)₂^-]₃ in an aqueous monomer solution, using femtosecond X-ray liquidography^9-12 with X-ray free-electron lasers^13,14. In the complex, which forms when three monomers A, B and C cluster together through non-covalent interactions^15,16, the distance between A and B is shorter than that between B and C. Tracking the wavepacket in three-dimensional nuclear coordinates reveals that within the first 60 femtoseconds after photoexcitation, a covalent bond forms between A and B to give A-B + C. The second covalent bond, between B and C, subsequently forms within 360 femtoseconds to give a linear and covalently bonded trimer complex A-B-C. The trimer exhibits harmonic vibrations that we map and unambiguously assign to specific normal modes using only the experimental data. In principle, more intense X-rays could visualize the motion not only of highly scattering atoms such as gold but also of lighter atoms such as carbon and nitrogen, which will open the door to the direct tracking of the atomic motions involved in many chemical reactions.

Ab Initio Calculation of Total X-ray Scattering from Molecules

We present a method to calculate total X-ray scattering cross sections directly from ab initio electronic wave functions in atoms and molecules. The approach can be used in conjunction with multiconfigurational wave functions and exploits analytical integrals of Gaussian-type functions over the scattering operator, which leads to accurate and efficient calculations. The results are validated by comparison to experimental results and previous theory for the molecules H₂ and CO₂. Importantly, we find that the inelastic component of the total scattering varies strongly with molecular geometry. The method is appropriate for use in conjunction with quantum molecular dynamics simulations for the analysis of new ultrafast X-ray scattering experiments and to interpret accurate gas-phase scattering experiments.

Ab Initio Computation of Rotationally-Averaged Pump-Probe X-ray and Electron Diffraction Signals

We develop a new algorithm for the computation of the rotationally averaged elastic molecular diffraction signal for the cases of perpendicular or parallel pump-probe geometries. The algorithm first collocates the charge density from an arbitrary ab initio wave function onto a Becke quadrature grid [A. Becke, J. Chem. Phys. 1988 , 88 , 2457 ], providing a high-fidelity multiresolution representation of the charge density. A double sum is then performed over the Becke grid points, and the interaction between points computed using the scattering kernels of Williamson and Zewail [J. C. Williamson and A. H. Zewail, J. Phys. Chem. 1994 , 98 , 2766 ]. These kernels analytically average over the molecular orientations with the cos² γ selection factor appropriate for one-photon dipole absorption in a perpendicular pump-probe geometry. We show that the method is converged with small grids containing <500 points/atom. We implement the algorithm on a GPU for increased efficiency and emonstrate the algorithm for molecules with up to a few dozen atoms. We explore the accuracy of the independent atom model (IAM) by comparison with our new and more accurate method. We also investigate the possibility of detecting signatures of electronic transitions in polyatomic pump-probe diffraction experiments.

Ab initio calculation of inelastic scattering

We calculate nonresonant inelastic electron and X-ray scattering cross sections for bound-to-bound transitions in atoms and molecules from ab initio electronic wavefunctions.

Imaging rotations and vibrations in polyatomic molecules with X-ray scattering

An approach for calculating elastic X-ray scattering from polyatomic molecules in specific electronic, vibrational, and rotational states is presented, and is used to consider the characterization of specific states in polyatomic molecules using elastic X-ray scattering.

Elastic X-ray scattering from state-selected molecules

The characterization of electronic, vibrational, and rotational states using elastic (coherent) X-ray scattering is considered. The scattering is calculated directly from complete active space self-consistent field level ab initio wavefunctions for H₂ molecules in the ground-state X¹Σ_g⁺ and first-excited EF¹Σ_g⁺ electronic states. The calculated scattering is compared to recent experimental measurements [Y.-W. Liu et al., Phys. Rev. A 89, 014502 (2014)], and the influence of vibrational and rotational states on the observed signal is examined. The scaling of the scattering calculations with basis set is quantified, and it is found that energy convergence of the ab initio calculations is a good indicator of the quality of the scattering calculations.

Ab Initio Calculation of Molecular Diffraction

We discuss the application of ab initio X-ray diffraction (AIXRD) to the interpretation of time-resolved and static X-ray diffraction. In our approach, elastic X-ray scattering is calculated directly from the ab initio multiconfigurational wave function via a Fourier transform of the electron density, using the first Born approximation for elastic scattering. Significant gains in efficiency can be obtained by performing the required Fourier transforms analytically, making it possible to combine the calculation of ab initio X-ray diffraction with expensive quantum dynamics simulations. We show that time-resolved X-ray diffraction can detect not only changes in molecular geometry but also changes in the electronic state of a molecule. Calculations for cis-, trans-, and cyclo-butadiene, as well as benzene and 1,3-cyclohexadiene are included.

Tracking of azobenzene isomerization by X-ray emission spectroscopy

Cis-trans isomerizations are among the fundamental processes in photochemistry. In azobenzene or its derivatives this dynamics is, due to its reversibility, one of the reactions widely used in photostimulation of molecular motors or in molecular electronics. Though intensively investigated in the optical regime, no detailed study exists in the X-ray regime so far. Because the X-ray emission spectroscopy echoes the electronic structure sensitive to the geometry, this theoretical report based on the density functional theory and its time-dependent version presents different nitrogen K-edge X-ray emission spectra for cis and trans isomers with close interrelation to their electron configuration. Considering the spectrum along the isomerization path, these structural signatures can be utilized to probe the isomerization dynamics in the excited molecule. The scheme can further be generalized to the element specific photoreactions.

Photoisomerization of Stilbene: The Detailed XMCQDPT2 Treatment

We report the detailed XMCQDPT2/cc-pVTZ study of trans-cis photoisomerization in one of the core systems of both experimental and computational photochemistry-the stilbene molecule. For the first time, the potential energy surface (PES) of the S1 state has been directly optimized and scanned using a multistate multiconfiguration second-order perturbation theory. We characterize the trans-stilbene, pyramidalized (phantom), and DHP-cis-stilbene geometric domains of the S1 state and describe their stationary points including the transition states between them, as well as S1/S0 intersections. Also reported are the minima and the activation barriers in the ground state. Our calculations correctly predict the kinetic isotope effect due to H/D exchange at ethylenic hydrogens, the dynamic behavior of excited cis-stilbene, and trans-cis branching ratio after relaxation to S0 through a rather unsymmetric conical intersection. In general, the XMCQDPT2 results confirm the qualitative adequacy of the TDDFT (especially SF-TDDFT) picture of the excited stilbene but also reveal quantitative discrepancies that deserve further exploration.

Photoisomerization of stilbene: a spin-flip density functional theory approach

The photoisomerization process of 1,2-diphenylethylene (stilbene) is investigated using the spin-flip density functional theory (SFDFT), which has recently been shown to be a promising approach for locating conical intersection (CI) points (Minezawa, N.; Gordon, M. S. J. Phys. Chem. A2009, 113, 12749). The SFDFT method gives valuable insight into twisted stilbene to which the linear response time-dependent DFT approach cannot be applied. In contrast to the previous SFDFT study of ethylene, a distinct twisted minimum is found for stilbene. The optimized structure has a sizable pyramidalization angle and strong ionic character, indicating that a purely twisted geometry is not a true minimum. In addition, the SFDFT approach can successfully locate two CI points: the twisted-pyramidalized CI that is similar to the ethylene counterpart and another CI that possibly lies on the cyclization pathway of cis-stilbene. The mechanisms of the cis--trans isomerization reaction are discussed on the basis of the two-dimensional potential energy surface along the twisting and pyramidalization angles.

Ab initio studies of ultrafast x-ray scattering of the photodissociation of iodine

We computationally examine various aspects of the reaction dynamics of the photodissociation and recombination of molecular iodine. We use our recently proposed formalism to calculate time-dependent x-ray scattering signal changes from first principles. Different aspects of the dynamics of this prototypical reaction are studied, such as coherent and noncoherent processes, features of structural relaxation that are periodic in time versus nonperiodic dissociative processes, as well as small electron density changes caused by electronic excitation, all with respect to x-ray scattering. We can demonstrate that wide-angle x-ray scattering offers a possibility to study the changes in electron densities in nonperiodic systems, which render it a suitable technique for the investigation of chemical reactions from a structural dynamics point of view.

Computational studies of the x-ray scattering properties of laser aligned stilbene

The enhancement of the x-ray scattering signal from partially aligned molecular samples is investigated. The alignment properties of the studied molecular system are modeled based on the method of laser alignment. With the advances in the area of laser alignment of molecules, the application of this sample manipulation technique promises a great potential for x-ray scattering measurements. Preferential alignment of molecules in an otherwise amorphous sample leads to constructive interference and thus increases the scattering intensity. This enhances the structural information encoded in the scattering images and enables improved resolution in studies of reaction dynamics, as in this work is shown for the example of the photo-isomerization of stilbene. We demonstrate that the scattering signal is strongly influenced by the alignment axis. Even the most basic one-dimensional alignment offers significant improvement compared to the structural information provided by a randomly oriented sample. Although the signal is sensitive to the uncertainty in the alignment angle, it offers encouraging results even at realistic alignment uncertainties.

Photocrystallography

This review describes the development and application of a new crystallographic technique that is starting to enable the three-dimensional structural determination of molecules in their photo-activated states. So called `photocrystallography' has wide applicability, particularly in the currently exciting area of photonics, and a discussion of this applied potential is put into context in this article. Studies are classified into four groups: photo-structural changes that are (i) irreversible; (ii) long-lived but reversible under certain conditions; (iii) transient with photo-active lifetimes of the order of microseconds; (iv) very short lived, existing at the nanosecond or even picosecond level. As photo-structural changes relative to the `ground state' can be subtle, this article necessarily concentrates on small-molecule single-crystal X-ray diffraction given that high atomic resolution is possible. That said, where it is pertinent, references are also made to related major advances in photo-induced macromolecular crystallography. The review concludes with an outlook on this new research area, including the future possibility of studying even more ephemeral, femtosecond-lived, photo-active species.