Comparison of ultrafast electron and X-ray diffraction – A computational study

Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe

Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.

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.

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

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.

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.

Ultrafast X-Ray Probes of Elementary Molecular Events

Elementary events that determine photochemical outcomes and molecular functionalities happen on the femtosecond and subfemtosecond timescales. Among the most ubiquitous events are the nonadiabatic dynamics taking place at conical intersections. These facilitate ultrafast, nonradiative transitions between electronic states in molecules that can outcompete slower relaxation mechanisms such as fluorescence. The rise of ultrafast X-ray sources, which provide intense light pulses with ever-shorter durations and larger observation bandwidths, has fundamentally revolutionized our spectroscopic capabilities to detect conical intersections. Recent theoretical studies have demonstrated an entirely new signature emerging once a molecule traverses a conical intersection, giving detailed insights into the coupled nuclear and electronic motions that underlie, facilitate, and ultimately determine the ultrafast molecular dynamics. Following a summary of current sources and experiments, we survey these techniques and provide a unified overview of their capabilities. We discuss their potential to dramatically increase our understanding of ultrafast photochemistry.

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.

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.

Ultrafast x-ray and electron scattering of free molecules: A comparative evaluation

Resolving gas phase molecular motions with simultaneous spatial and temporal resolution is rapidly coming within the reach of x-ray Free Electron Lasers (XFELs) and Mega-electron-Volt (MeV) electron beams. These two methods enable scattering experiments that have yielded fascinating new results, and while both are important methods for determining transient molecular structures in photochemical reactions, it is important to understand their relative merits. In the present study, we evaluate the respective scattering cross sections of the two methods and simulate their ability to determine excited state molecular structures in light of currently existing XFEL and MeV source parameters. Using the example of optically excited N-methyl morpholine and simulating the scattering patterns with shot noise, we find that the currently achievable signals are superior with x-ray scattering for equal samples and on a per-shot basis and that x-ray scattering requires fewer detected signal counts for an equal fidelity structure determination. Importantly, within the independent atom model, excellent structure determinations can be achieved for scattering vectors only to about 5 Å-1, leaving larger scattering vector ranges for investigating vibrational motions and wavepackets. Electron scattering has a comparatively higher sensitivity toward hydrogen atoms, which may point to applications where electron scattering is inherently the preferred choice, provided that excellent signals can be achieved at large scattering angles that are currently difficult to access.

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.

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.