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

Prediction of photodynamics of 200 nm excited cyclobutanone with linear response electronic structure and ab initio multiple spawning

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S2). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. We use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

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.

Predicting the photodynamics of cyclobutanone triggered by a laser pulse at 200 nm and its MeV-UED signals-A trajectory surface hopping and XMS-CASPT2 perspective

This work is part of a prediction challenge that invited theoretical/computational chemists to predict the photochemistry of cyclobutanone in the gas phase, excited at 200 nm by a laser pulse, and the expected signal that will be recorded during a time-resolved megaelectronvolt ultrafast electron diffraction (MeV-UED). We present here our theoretical predictions based on a combination of trajectory surface hopping with XMS-CASPT2 (for the nonadiabatic molecular dynamics) and Born-Oppenheimer molecular dynamics with MP2 (for the athermal ground-state dynamics following internal conversion), coined (NA+BO)MD. The initial conditions were sampled from Born-Oppenheimer molecular dynamics coupled to a quantum thermostat. Our simulations indicate that the main photoproducts after 2 ps of dynamics are CO + cyclopropane (50%), CO + propene (10%), and ethene and ketene (34%). The photoexcited cyclobutanone in its second excited electronic state S2 can follow two pathways for its nonradiative decay: (i) a ring-opening in S2 and a subsequent rapid decay to the ground electronic state, where the photoproducts are formed, or (ii) a transfer through a closed-ring conical intersection to S1, where cyclobutanone ring opens and then funnels to the ground state. Lifetimes for the photoproduct and electronic populations were determined. We calculated a stationary MeV-UED signal [difference pair distribution function-ΔPDF(r)] for each (interpolated) pathway as well as a time-resolved signal [ΔPDF(r,t) and ΔI/I(s,t)] for the full swarm of (NA+BO)MD trajectories. Furthermore, our analysis provides time-independent basis functions that can be used to fit the time-dependent experimental UED signals [both ΔPDF(r,t) and ΔI/I(s,t)] and potentially recover the population of photoproducts. We also offer a detailed analysis of the limitations of our model and their potential impact on the predicted experimental signals.

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?

Quantum tomography of molecules using ultrafast electron diffraction

We propose a quantum tomography (QT) approach to retrieve the temporally evolving reduced density matrix in electronic state basis, where the populations and coherence between the ground state and excited state are reconstructed from the ultrafast electron diffraction signal. In order to showcase the capability of the proposed QT approach, we simulate the nuclear wavepacket dynamics and ultrafast electron diffraction of photoexcited pyrrole molecules using the ab initio quantum chemical CASSCF method. From the simulated time-resolved diffraction data, we retrieve the evolving density matrix in a crude diabatic representation basis and reveal the symmetry of the excited pyrrole wavepacket. Our QT approach opens the route to make a quantum version of "molecular movie" that covers the electronic degree of freedom and equips ultrafast electron diffraction with the power to reveal the coherence between electronic states, relaxation, and dynamics of population transfer.

Monitoring the Evolution of Relative Product Populations at Early Times during a Photochemical Reaction

Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps toward understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species among the emerging reaction products. Here, we show that mega-electronvolt ultrafast electron diffraction in combination with ab initio molecular dynamics calculations offer a powerful route to determining time-resolved populations of the various isomeric products formed after UV (266 nm) excitation of the five-membered heterocyclic molecule 2(5H)-thiophenone. This strategy provides experimental validation of the predicted high (∼50%) yield of an episulfide isomer containing a strained three-membered ring within ∼1 ps of photoexcitation and highlights the rapidity of interconversion between the rival highly vibrationally excited photoproducts in their ground electronic state.

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.

Molecular Chirality and Its Monitoring by Ultrafast X-ray Pulses

Major advances in X-ray sources including the development of circularly polarized and orbital angular momentum pulses make it possible to probe matter chirality at unprecedented energy regimes and with Ångström and femtosecond spatiotemporal resolutions. We survey the theory of stationary and time-resolved nonlinear chiral measurements that can be carried out in the X-ray regime using tabletop X-ray sources or large scale (XFEL, synchrotron) facilities. A variety of possible signals and their information content are discussed.

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.

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 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.

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.

Probing Delocalized Current Densities in Selenophene by Resonant X-ray Sum-Frequency Generation

Time-resolved, resonant X-ray sum-frequency generation in aligned selenophene molecules is calculated. A wave packet of valence-excited states, prepared by an extreme-ultraviolet pump pulse, is probed by two 12-keV X-ray probe pulses resonant with the Se core-excited states for variable time delays. At these hard-X-ray frequencies, the angström wavelength of the X-ray probe is comparable to the molecular size. We thus employ a nonlocal description of the light-matter interaction based on the minimal-coupling Hamiltonian. The wavevector-resolved resonant stimulated sum-frequency-generation signal, obtained by varying the propagation direction of hard-X-ray pulses, can thus directly monitor the transition current densities between core and ground/valence states. This is in contrast to off-resonant diffraction, which detects the transition charge densities.

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.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

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.

Transient X-ray Absorption Spectral Fingerprints of the S1 Dark State in Uracil

Low-lying dark nπ* states play an important role in many photophysical and photochemical processes of organic chromophores. Transient X-ray absorption spectroscopy (TXAS) provides a powerful technique for probing the dynamics of valence states by exciting the electrons into high-lying core excited states. We employ multiconfigurational self-consistent field calculations to investigate the TXAS of uracil along its nonradiative photodecay pathways. An open issue is whether dark nπ* state S₁ (n is the lone pair localized on an oxygen atom) is accessible when bright ππ* state S₂ is selectively excited. Vertical core excitations were calculated along the potential energy surfaces of the three lowest states, S₀-S₂, interpolated between two minima and two minimum-energy conical intersections. Computed TXAS data from the C, N, and O K edges show distinct spectral fingerprints of the dark state in all spectral regimes. At the O 1s edge, the nπ* state has a very strong absorption at 526-527 eV, while at the C (N) 1s edge, by contrast, there is almost zero (very weak) absorption at 279-282 eV (397-398 eV). All K-edge spectra can be used to sensitively detect the dark states. Our proposed O 1s feature has already been observed in a recent TXAS experiment with thymine. Natural transition orbital analysis is used to interpret all dominant features of the three lowest-valence states along the reaction coordinate and reveal some important valence fine-structure information from the core excitation.

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 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.