Ultrafast X-ray Scattering from Molecules

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Ultrafast electron diffraction of photoexcited gas-phase cyclobutanone predicted by ab initio multiple cloning simulations

We present the result of our calculations of ultrafast electron diffraction (UED) for cyclobutanone excited into the S2 electronic state, which is based on the non-adiabatic dynamics simulations with the Ab Initio Multiple Cloning (AIMC) method with the electronic structure calculated at the SA(3)-CASSCF(12,12)/aug-cc-pVDZ level of theory. The key features in the UED pattern were identified, which can be used to distinguish between the reaction pathways observed in the AIMC dynamics, although there is a significant overlap between representative signals due to the structural similarity of the products. The calculated UED pattern can be compared with the experiment.

Extracting the electronic structure signal from X-ray and electron scattering in the gas phase

X-ray and electron scattering from free gas-phase molecules is examined using the independent atom model (IAM) and ab initio electronic structure calculations. The IAM describes the effect of the molecular geometry on the scattering, but does not account for the redistribution of valence electrons due to, for instance, chemical bonding. By examining the total, i.e. energy-integrated, scattering from three molecules, fluoroform (CHF3), 1,3-cyclohexadiene (C6H8) and naphthalene (C10H8), the effect of electron redistribution is found to predominantly reside at small-to-medium values of the momentum transfer (q ≤ 8 Å-1) in the scattering signal, with a maximum percent difference contribution at 2 ≤ q ≤ 3 Å-1. A procedure to determine the molecular geometry from the large-q scattering is demonstrated, making it possible to more clearly identify the deviation of the scattering from the IAM approximation at small and intermediate q and to provide a measure of the effect of valence electronic structure on the scattering signal.

Automatic Clustering of Excited-State Trajectories: Application to Photoexcited Dynamics

We introduce automatic clustering as a computationally efficient tool for classifying and interpreting trajectories from simulations of photo-excited dynamics. Trajectories are treated as time-series data, with the features for clustering selected by variance mapping of normalized data. The L2-norm and dynamic time warping are proposed as suitable similarity measures for calculating the distance matrices, and these are clustered using the unsupervised density-based DBSCAN algorithm. The silhouette coefficient and the number of trajectories classified as noise are used as quality measures for the clustering. The ability of clustering to provide rapid overview of large and complex trajectory data sets, and its utility for extracting chemical and physical insight, is demonstrated on trajectories corresponding to the photochemical ring-opening reaction of 1,3-cyclohexadiene, noting that the clustering can be used to generate reduced dimensionality representations in an unbiased manner.

Deciphering the Influence of Ground-State Distributions on the Calculation of Photolysis Observables

Nonadiabatic molecular dynamics offers a powerful tool for studying the photochemistry of molecular systems. Key to any nonadiabatic molecular dynamics simulation is the definition of its initial conditions (ICs), ideally representing the initial molecular quantum state of the system of interest. In this work, we provide a detailed analysis of how ICs may influence the calculation of experimental observables by focusing on the photochemistry of methylhydroperoxide (MHP), the simplest and most abundant organic peroxide in our atmosphere. We investigate the outcome of trajectory surface hopping simulations for distinct sets of ICs sampled from different approximate quantum distributions, namely harmonic Wigner functions and ab initio molecular dynamics using a quantum thermostat (QT). Calculating photoabsorption cross-sections, quantum yields, and translational kinetic energy maps from the results of these simulations reveals the significant effect of the ICs, in particular when low-frequency (∼ a few hundred cm-1) normal modes are connected to the photophysics of the molecule. Overall, our results indicate that sampling ICs from ab initio molecular dynamics using a QT is preferable for flexible molecules with photoactive low-frequency modes. From a photochemical perspective, our nonadiabatic dynamics simulations offer an explanation for a low-energy tail observed at high excitation energy in the translational kinetic energy map of MHP.

Mapping the photochemistry of cyclopentadiene: from theory to ultrafast X-ray scattering

The photoinduced ring-conversion reaction when cyclopentadiene (CP) is excited at 5.10 eV is simulated using surface-hopping semiclassical trajectories with XMS(3)-CASPT2(4,4)/cc-pVDZ electronic structure theory. In addition, PBE0/def2-SV(P) is employed for ground state propagation of the trajectories. The dynamics is propagated for 10 ps, mapping both the nonadiabatic short-time dynamics (<300 fs) and the increasingly statistical dynamics on the electronic ground state. The short-time dynamics yields a mixture of hot CP and bicyclo[2.1.0]pentene (BP), with the two products reached via different regions of the same conical intersection seam. On the ground state, we observe slow conversion from BP to CP which is modelled by RRKM theory with a transition state determined using PBE0/def2-TZVP. The CP products are furthermore associated with ground state hydrogen shifts and some H-atom dissociation. Finally, the prospects for detailed experimental mapping using novel ultrafast X-ray scattering experiments are discussed and observables for such experiments are predicted. In particular, we assess the possibility of retrieving electronic states and their populations alongside the structural dynamics.

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.

Monitoring vibronic coherences and molecular aromaticity in photoexcited cyclooctatetraene with an X-ray probe: a simulation study

Understanding conical intersection (CI) dynamics and subsequent conformational changes is key for exploring and controlling photo-reactions in aromatic molecules. Monitoring of their time-resolved dynamics remains a formidable experimental challenge. In this study, we simulate the photoinduced S₃ to S₁ non-adiabatic dynamics of cyclooctatetraene (COT), involving multiple CIs with relaxation times in good agreement with experiment. We further investigate the possibility to directly probe the CI passages in COT by off-resonant X-ray Raman spectroscopy (TRUECARS) and time-resolved X-ray diffraction (TRXD). We find that these signals sensitively monitor key chemical features during the ultrafast dynamics. First, we distinguish two CIs by using TRUECARS signals with their appearances at different Raman shifts. Second, we demonstrate that TRXD, where X-ray photons scatter off electron densities, can resolve ultrafast changes in the aromaticity of COT. It can further distinguish between planar and non-planar geometries explored during the dynamics, as e.g. two different tetraradical-type CIs. The knowledge gained from these measurements can give unique insight into fundamental chemical properties that dynamically change during non-adiabatic passages.

Transient vibration and product formation of photoexcited CS2 measured by time-resolved x-ray scattering

We have observed details of the internal motion and dissociation channels in photoexcited carbon disulfide (CS2) using time-resolved x-ray scattering (TRXS). Photoexcitation of gas-phase CS2 with a 200 nm laser pulse launches oscillatory bending and stretching motion, leading to dissociation of atomic sulfur in under a picosecond. During the first 300 fs following excitation, we observe significant changes in the vibrational frequency as well as some dissociation of the C-S bond, leading to atomic sulfur in the both 1D and 3P states. Beyond 1400 fs, the dissociation is consistent with primarily 3P atomic sulfur dissociation. This channel-resolved measurement of the dissociation time is based on our analysis of the time-windowed dissociation radial velocity distribution, which is measured using the temporal Fourier transform of the TRXS data aided by a Hough transform that extracts the slopes of linear features in an image. The relative strength of the two dissociation channels reflects both their branching ratio and differences in the spread of their dissociation times. Measuring the time-resolved dissociation radial velocity distribution aids the resolution of discrepancies between models for dissociation proposed by prior photoelectron spectroscopy work.

A laboratory frame density matrix for ultrafast quantum molecular dynamics

In most cases, the ultrafast dynamics of resonantly excited molecules are considered and almost always computed in the molecular frame, while experiments are carried out in the laboratory frame. Here, we provide a formalism in terms of a lab frame density matrix, which connects quantum dynamics in the molecular frame to those in the laboratory frame, providing a transparent link between computation and measurement. The formalism reveals that in any such experiment, the molecular frame dynamics vary for molecules in different orientations and that certain coherences, which are potentially experimentally accessible, are rejected by the orientation-averaged reduced vibronic density matrix. Instead, molecular angular distribution moments are introduced as a more accurate representation of experimentally accessible information. Furthermore, the formalism provides a clear definition of a molecular frame quantum tomography and specifies the requirements to perform such a measurement enabling the experimental imaging of molecular frame vibronic dynamics. Successful completion of such a measurement fully characterizes the molecular frame quantum dynamics for a molecule at any orientation in the laboratory frame.

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.

The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes

This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.

Imaging conical intersection dynamics during azobenzene photoisomerization by ultrafast X-ray diffraction

X-ray diffraction is routinely used for structure determination of stationary molecular samples. Modern X-ray photon sources, e.g., from free-electron lasers, enable us to add temporal resolution to these scattering events, thereby providing a movie of atomic motions. We simulate and decipher the various contributions to the X-ray diffraction pattern for the femtosecond isomerization of azobenzene, a textbook photochemical process. A wealth of information is encoded besides real-time monitoring of the molecular charge density for the cis to trans isomerization. In particular, vibronic coherences emerge at the conical intersection, contributing to the total diffraction signal by mixed elastic and inelastic photon scattering. They cause distinct phase modulations in momentum space, which directly reflect the real-space phase modulation of the electronic transition density during the nonadiabatic passage. To overcome the masking by the intense elastic scattering contributions from the electronic populations in the total diffraction signal, we discuss how this information can be retrieved, e.g., by employing very hard X-rays to record large scattering momentum transfers.

Signatures of electronic and nuclear coherences in ultrafast molecular x-ray and electron diffraction

Femtosecond x-ray and electron diffraction hold promise to image the evolving structures of single molecules. We present a unified quantum-electrodynamical formulation of diffraction signals, based on the exact many-body nuclear + electronic wavefunction that can be extracted from quantum chemistry simulations. This gives a framework for analyzing various approximate molecular dynamics simulations. We show that the complete description of ultrafast diffraction signals contains interesting contributions involving mixed elastic and inelastic scattered photons that are usually masked by other larger contributions and are neglected. These terms include overlaps of nuclear wavepackets between different electronic states that provide an electronic decoherence mechanism and are important for the time-resolved imaging of conical intersections.

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.

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 Fragment Method for Calculating Molecular X-ray Diffraction

A fragment-based approach for the prediction of elastic X-ray scattering is presented. The total diffraction pattern is assembled from anisotropic form factors calculated for individual molecular fragments, optionally including corrections for pairwise interactions between fragments. The approach is evaluated against full ab initio scattering calculations in the peptide diphenylalanine, and the optimal selection of fragments is examined in the ethanol molecule. The approach is found to improve significantly on the independent atom model while remaining conceptually simple and computationally efficient. It is expected to be particularly useful for macromolecules with repeated subunits, such as peptides, proteins, DNA, or RNA and other polymers, where it is straightforward to define appropriate fragments.

Electronic Coherence in Ultrafast X-Ray Scattering from Molecular Wave Packets

Simulations of nonresonant ultrafast x-ray scattering from a molecular wave packet in H_{2} are used to examine and classify the components that contribute to the total scattering signal. The elastic component, which can be used to determine the structural dynamics of the molecule, is also found to carry a strong signature of an adiabatic electron transfer that occurs in the simulated molecule. The inelastic component, frequently assumed to be constant, is found to change with the geometry of the molecule. Finally, a coherent mixed component due to interferences between different inelastic transitions is identified and shown to provide a direct probe of transient electronic coherences.

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.

Correspondence between electronic structure calculations and simulations: nonadiabatic dynamics in CS2

The choice of ab initio electronic structure method is an important factor in determining the fidelity of nonadiabatic dynamics simulations.

Effects of probe energy and competing pathways on time-resolved photoelectron spectroscopy: the ring-opening of 1,3-cyclohexadiene

Photoelectron spectra for the ring-opening dynamics of 1,3-cyclohexadiene are studied using a model based on quantum molecular dynamics and the Dyson orbital approach.

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.

Time-resolved X-ray scattering by electronic wave packets: analytic solutions to the hydrogen atom

This paper demonstrates how the time-dependent scattering signal of electronic wave packets in the hydrogen atom can be expressed analytically.

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.

Beyond structure: ultrafast X-ray absorption spectroscopy as a probe of non-adiabatic wavepacket dynamics

The excited state non-adiabatic dynamics of polyatomic molecules, leading to the coupling of structural and electronic dynamics, is a fundamentally important yet challenging problem for both experiment and theory. Ongoing developments in ultrafast extreme vacuum ultraviolet (XUV) and soft X-ray sources present new probes of coupled electronic-structural dynamics because of their novel and desirable characteristics. As one example, inner-shell spectroscopy offers localized, atom-specific probes of evolving electronic structure and bonding (via chemical shifts). In this work, we present the first on-the-fly ultrafast X-ray time-resolved absorption spectrum simulations of excited state wavepacket dynamics: photo-excited ethylene. This was achieved by coupling the ab initio multiple spawning (AIMS) method, employing on-the-fly dynamics simulations, with high-level algebraic diagrammatic construction (ADC) X-ray absorption cross-section calculations. Using the excited state dynamics of ethylene as a test case, we assessed the ability of X-ray absorption spectroscopy to project out the electronic character of complex wavepacket dynamics, and evaluated the sensitivity of the calculated spectra to large amplitude nuclear motion. In particular, we demonstrate the pronounced sensitivity of the pre-edge region of the X-ray absorption spectrum to the electronic and structural evolution of the excited-state wavepacket. We conclude that ultrafast time-resolved X-ray absorption spectroscopy may become a powerful tool in the interrogation of excited state non-adiabatic molecular dynamics.

Femtosecond photodissociation dynamics of 1,4-diiodobenzene by gas-phase X-ray scattering and photoelectron spectroscopy

We present a multifaceted investigation into the initial photodissociation dynamics of 1,4-diiodobenzene (DIB) following absorption of 267 nm radiation. We combine ultrafast time-resolved photoelectron spectroscopy and X-ray scattering experiments performed at the Linac Coherent Light Source (LCLS) to study the initial electronic excitation and subsequent rotational alignment, and interpret the experiments in light of Complete Active Space Self-Consistent Field (CASSCF) calculations of the excited electronic landscape. The initially excited state is found to be a bound ¹B₁ surface, which undergoes ultrafast population transfer to a nearby state in 35 ± 10 fs. The internal conversion most likely leads to one or more singlet repulsive surfaces that initiate the dissociation. This initial study is an essential and prerequisite component of a comprehensive study of the complete photodissociation pathway(s) of DIB at 267 nm. Assignment of the initially excited electronic state as a bound state identifies the mechanism as predissociative, and measurement of its lifetime establishes the time between excitation and initiation of dissociation, which is crucial for direct comparison of photoelectron and scattering experiments.