Ab Initio Calculation of Total 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.

A comparative review of time-resolved x-ray and electron scattering to probe structural dynamics

The structure of molecules, particularly the dynamic changes in structure, plays an essential role in understanding physical and chemical phenomena. Time-resolved (TR) scattering techniques serve as crucial experimental tools for studying structural dynamics, offering direct sensitivity to molecular structures through scattering signals. Over the past decade, the advent of x-ray free-electron lasers (XFELs) and mega-electron-volt ultrafast electron diffraction (MeV-UED) facilities has ushered TR scattering experiments into a new era, garnering significant attention. In this review, we delve into the basic principles of TR scattering experiments, especially focusing on those that employ x-rays and electrons. We highlight the variations in experimental conditions when employing x-rays vs electrons and discuss their complementarity. Additionally, cutting-edge XFELs and MeV-UED facilities for TR x-ray and electron scattering experiments and the experiments performed at those facilities are reviewed. As new facilities are constructed and existing ones undergo upgrades, the landscape for TR x-ray and electron scattering experiments is poised for further expansion. Through this review, we aim to facilitate the effective utilization of these emerging opportunities, assisting researchers in delving deeper into the intricate dynamics of molecular structures.

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.

Femtosecond Electronic and Hydrogen Structural Dynamics in Ammonia Imaged with Ultrafast Electron Diffraction

Directly imaging structural dynamics involving hydrogen atoms by ultrafast diffraction methods is complicated by their low scattering cross sections. Here we demonstrate that megaelectronvolt ultrafast electron diffraction is sufficiently sensitive to follow hydrogen dynamics in isolated molecules. In a study of the photodissociation of gas phase ammonia, we simultaneously observe signatures of the nuclear and corresponding electronic structure changes resulting from the dissociation dynamics in the time-dependent diffraction. Both assignments are confirmed by ab initio simulations of the photochemical dynamics and the resulting diffraction observable. While the temporal resolution of the experiment is insufficient to resolve the dissociation in time, our results represent an important step towards the observation of proton dynamics in real space and time.

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.

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.

Efficient Computation of Two-Electron Reduced Density Matrices via Selected Configuration Interaction

We create an approach to efficiently calculate two-electron reduced density matrices (2-RDMs) using selected configuration interaction wavefunctions. This is demonstrated using the specific example of Monte Carlo configuration interaction (MCCI). The computation of the 2-RDMs is accelerated by using ideas from fast implementations of full configuration interaction (FCI) and recent advances in implementing the Slater–Condon rules using hardware bitwise operations. This method enables a comparison of MCCI and truncated CI 2-RDMs with FCI values for a range of molecules, which includes stretched bonds and excited states. The accuracy in energies, wavefunctions, and 2-RDMs is seen to exhibit a similar behavior. We find that MCCI can reach sufficient accuracy of the 2-RDM using significantly fewer configurations than truncated CI, particularly for systems with strong multireference character.

Direct Monitoring of Conical Intersection Passage via Electronic Coherences in Twisted X-Ray Diffraction

Quantum coherences in electronic motions play a critical role in determining the pathways and outcomes of virtually all photophysical and photochemical molecular processes. However, the direct observation of electronic coherences in the vicinity of conical intersections remains a formidable challenge. We propose a novel time-resolved twisted x-ray diffraction technique that can directly monitor the electronic coherences created as the molecule passes through a conical intersection. We show that the contribution of electronic populations to this signal is canceled out when using twisted x-ray beams that carry a light orbital angular momentum, providing a direct measurement of transient electronic coherences in gas-phase molecules.

Optically Induced Anisotropy in Time-Resolved Scattering: Imaging Molecular-Scale Structure and Dynamics in Disordered Media with Experiment and Theory

Time-resolved scattering experiments enable imaging of materials at the molecular scale with femtosecond time resolution. However, in disordered media they provide access to just one radial dimension thus limiting the study of orientational structure and dynamics. Here we introduce a rigorous and practical theoretical framework for predicting and interpreting experiments combining optically induced anisotropy and time-resolved scattering. Using impulsive nuclear Raman and ultrafast x-ray scattering experiments of chloroform and simulations, we demonstrate that this framework can accurately predict and elucidate both the spatial and temporal features of these experiments.

Multichannel photodissociation dynamics in CS2 studied by ultrafast electron diffraction

The structural dynamics of photoexcited gas-phase carbon disulfide (CS₂) molecules are investigated using ultrafast electron diffraction. The dynamics were triggered by excitation of the optically bright ¹B₂(¹Σ_u⁺) state by an ultraviolet femtosecond laser pulse centred at 200 nm. In accordance with previous studies, rapid vibrational motion facilitates a combination of internal conversion and intersystem crossing to lower-lying electronic states. Photodissociation via these electronic manifolds results in the production of CS fragments in the electronic ground state and dissociated singlet and triplet sulphur atoms. The structural dynamics are extracted from the experiment using a trajectory-fitting filtering approach, revealing the main characteristics of the singlet and triplet dissociation pathways. Finally, the effect of the time-resolution on the experimental signal is considered and an outlook to future experiments provided.

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.

Determination of excited state molecular structures from time-resolved gas-phase X-ray scattering

We present a comprehensive investigation of a recently introduced method to determine transient structures of molecules in excited electronic states with sub-ångstrom resolution from time-resolved gas-phase scattering signals. The method, which is examined using time-resolved X-ray scattering data measured on the molecule N-methylmorpholine (NMM) at the Linac Coherent Light Source (LCLS), compares the experimentally measured scattering patterns against the simulated patterns corresponding to a large pool of molecular structures to determine the full set of structural parameters. In addition, we examine the influence of vibrational state distributions and find the effect negligible within the current experimental detection limits, despite that the molecules have a comparatively high internal vibrational energy. The excited state structures determined using three structure pools generated using three different computational methods are in good agreement, demonstrating that the procedure is largely independent of the computational chemistry method employed as long as the pool is sufficiently expansive in the vicinity of the sought structure and dense enough to yield good matches to the experimental patterns.

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.

Ultraintense, ultrashort pulse X-ray scattering in small molecules

We examine X-ray scattering from an isolated organic molecule from the linear to nonlinear absorptive regime. In the nonlinear regime, we explore the importance of both the coherent and incoherent channels and observe the onset of nonlinear behavior as a function of pulse duration and energy. In the linear regime, we test the sensitivity of the scattering signal to molecular bonding and electronic correlation via calculations using the independent atom model (IAM), Hartree-Fock (HF) and density functional theory (DFT). Finally, we describe how coherent X-ray scattering can be used to directly visualize femtosecond charge transfer and dissociation within a single molecule undergoing X-ray multiphoton absorption.

Ultrafast X-ray scattering offers a structural view of excited-state charge transfer

Intramolecular charge transfer and the associated changes in molecular structure in N,N'-dimethylpiperazine are tracked using femtosecond gas-phase X-ray scattering. The molecules are optically excited to the 3p state at 200 nm. Following rapid relaxation to the 3s state, distinct charge-localized and charge-delocalized species related by charge transfer are observed. The experiment determines the molecular structure of the two species, with the redistribution of electron density accounted for by a scattering correction factor. The initially dominant charge-localized state has a weakened carbon-carbon bond and reorients one methyl group compared with the ground state. Subsequent charge transfer to the charge-delocalized state elongates the carbon-carbon bond further, creating an extended 1.634 Å bond, and also reorients the second methyl group. At the same time, the bond lengths between the nitrogen and the ring-carbon atoms contract from an average of 1.505 to 1.465 Å. The experiment determines the overall charge transfer time constant for approaching the equilibrium between charge-localized and charge-delocalized species to 3.0 ps.

Simultaneous observation of nuclear and electronic dynamics by ultrafast electron diffraction

Simultaneous observation of nuclear and electronic motion is crucial for a complete understanding of molecular dynamics in excited electronic states. It is challenging for a single experiment to independently follow both electronic and nuclear dynamics at the same time. Here we show that ultrafast electron diffraction can be used to simultaneously record both electronic and nuclear dynamics in isolated pyridine molecules, naturally disentangling the two components. Electronic state changes (S₁→S₀ internal conversion) were reflected by a strong transient signal in small-angle inelastic scattering, and nuclear structural changes (ring puckering) were monitored by large-angle elastic diffraction. Supported by ab initio nonadiabatic molecular dynamics and diffraction simulations, our experiment provides a clear view of the interplay between electronic and nuclear dynamics of the photoexcited pyridine molecule.

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.

Observation of the molecular response to light upon photoexcitation

When a molecule interacts with light, its electrons can absorb energy from the electromagnetic field by rapidly rearranging their positions. This constitutes the first step of photochemical and photophysical processes that include primary events in human vision and photosynthesis. Here, we report the direct measurement of the initial redistribution of electron density when the molecule 1,3-cyclohexadiene (CHD) is optically excited. Our experiments exploit the intense, ultrashort hard x-ray pulses of the Linac Coherent Light Source (LCLS) to map the change in electron density using ultrafast x-ray scattering. The nature of the excited electronic state is identified with excellent spatial resolution and in good agreement with theoretical predictions. The excited state electron density distributions are thus amenable to direct experimental observation.

Photoabsorption is a fundamental process that leads to changes in the electron density in matter. Here, the authors show a direct measurement of the distribution of electron density when a cyclohexadine molecule is excited by pulsed UV radiation and probed by a time delayed X-ray pulse generated at LCLS.

On the use of the Obara–Saika recurrence relations for the calculation of structure factors in quantum crystallography

Modern methods of quantum crystallography are techniques firmly rooted in quantum chemistry and, as in many quantum chemical strategies, electron densities are expressed as two-centre expansions that involve basis functions centred on atomic nuclei. Therefore, the computation of the necessary structure factors requires the evaluation of Fourier transform integrals of basis function products. Since these functions are usually Cartesian Gaussians, in this communication it is shown that the Fourier integrals can be efficiently calculated by exploiting an extension of the Obara–Saika recurrence formulas, which are successfully used by quantum chemists in the computation of molecular integrals. Implementation and future perspectives of the technique are also discussed.

Scattering off molecules far from equilibrium

Pump-probe gas phase X-ray scattering experiments, enabled by the development of X-ray free electron lasers, have advanced to reveal scattering patterns of molecules far from their equilibrium geometry. While dynamic displacements reflecting the motion of wavepackets can probe deeply into the reaction dynamics, in many systems, the thermal excitation embedded in the molecules upon optical excitation and energy randomization can create systems that encompass structures far from the ground state geometry. For polyatomic molecular systems, large amplitude vibrational motions are associated with anharmonicity and shifts of interatomic distances, making analytical solutions using traditional harmonic approximations inapplicable. More generally, the interatomic distances in a polyatomic molecule are not independent and the traditional equations commonly used to interpret the data may give unphysical results. Here, we introduce a novel method based on molecular dynamic trajectories and illustrate it on two examples of hot, vibrating molecules at thermal equilibrium. When excited at 200 nm, 1,3-cyclohexadiene (CHD) relaxes on a subpicosecond time scale back to the reactant molecule, the dominant pathway, and to various forms of 1,3,5-hexatriene (HT). With internal energies of about 6 eV, the energy thermalizes quickly, leading to structure distributions that deviate significantly from their vibrationless equilibrium. The experimental and theoretical results are in excellent agreement and reveal that a significant contribution to the scattering signal arises from transition state structures near the inversion barrier of CHD. In HT, our analysis clarifies that previous inconsistent structural parameters determined by electron diffraction were artifacts that might have resulted from the use of inapplicable analytical equations.

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.