Simulating Coherent Multidimensional Spectroscopy of Nonadiabatic Molecular Processes: From the Infrared to the X-ray Regime

Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation

The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.

Conical Intersections at Interfaces Revealed by Phase-Cycling Interface-Specific Two-Dimensional Electronic Spectroscopy (i2D-ES)

Conical intersections (CIs) hold significant stake in manipulating and controlling photochemical reaction pathways of molecules at interfaces and surfaces by affecting molecular dynamics therein. Currently, there is no tool for characterizing CIs at interfaces and surfaces. To this end, we have developed phase-cycling interface-specific two-dimensional electronic spectroscopy (i2D-ES) and combined it with advanced computational modeling to explore nonadiabatic CI dynamics of molecules at the air/water interface. Specifically, we integrated the phase locked pump pulse pair with an interface-specific electronic probe to obtain the two-dimensional interface-specific responses. We demonstrate that the nonadiabatic transitions of an interface-active azo dye molecule that occur through the CIs at the interface have different kinetic pathways from those in the bulk water. Upon photoexcitation, two CIs are present: one from an intersection of an optically active S2 state with a dark S1 state and the other from the intersection of the progressed S1 with the ground state S0. We find that the molecular conformations in the ground state are different for interfacial molecules. The interfacial molecules are intimately correlated with the locally populated excited state S2 being farther away from the CI region. This leads to slower nonadiabatic dynamics at the interface than in bulk water. Moreover, we show that the nonadiabatic transition from the S1 dark state to the ground state is significantly longer at the interface than that in the bulk, which is likely due to the orientationally restricted configuration of the excited state at the interface. Our findings suggest that orientational configurations of molecules manipulate reaction pathways at interfaces and surfaces.

The photodissociation dynamics and ultrafast electron diffraction image of cyclobutanone from the surface hopping dynamics simulation

The comprehension of nonadiabatic dynamics in polyatomic systems relies heavily on the simultaneous advancements in theoretical and experimental domains. The gas-phase ultrafast electron diffraction (UED) technique has attracted significant attention as a unique tool for monitoring photochemical and photophysical processes at the all-atomic level with high temporal and spatial resolutions. In this work, we simulate the UED spectra of cyclobutanone using the trajectory surface hopping method at the extended multi-state complete active space second order perturbation theory (XMS-CASPT2) level and thereby predict the results of the upcoming UED experiments in the Stanford Linear Accelerator Laboratory. The simulated results demonstrate that a few pathways, including the C2 and C3 dissociation channels, as well as the ring opening channel, play important roles in the nonadiabatic reactions of cyclobutanone. We demonstrate that the simulated UED signal can be directly interpreted in terms of atomic motions, which provides a unique way of monitoring the evolution of the molecular structure in real time. Our work not only provides numerical data that help to determine the accuracy of the well-known surface hopping dynamics at the high XMS-CASPT2 electronic-structure level but also facilitates the understanding of the microscopic mechanisms of the photoinduced reactions in cyclobutanone.

Stimulated X-ray emission spectroscopy

We describe an emerging hard X-ray spectroscopy technique, stimulated X-ray emission spectroscopy (S-XES). S-XES has the potential to characterize the electronic structure of 3d transition metal complexes with spectral information currently not reachable and might lead to the development of new ultrafast X-ray sources with properties beyond the state of the art. S-XES has become possible with the emergence of X-ray free-electron lasers (XFELs) that provide intense femtosecond X-ray pulses that can be employed to generate a population inversion of core-hole excited states resulting in stimulated X-ray emission. We describe the instrumentation, the various types of S-XES, the potential applications, the experimental challenges, and the feasibility of applying S-XES to characterize dilute systems, including the Mn4Ca cluster in the oxygen evolving complex of photosystem II.

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?

Evidence of Excited-State Vibrational Mode Governing the Photorelaxation of a Charge-Transfer Complex

Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical-computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ "rocking" and "bending" modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.

Coherent x-ray spontaneous emission spectroscopy of conical intersections

Conical intersections are known to play a vital role in many photochemical processes. The breakdown of the Born-Oppenheimer approximation in the vicinity of a conical intersection causes exciting phenomena, such as the ultrafast radiationless decay of excited states. The passage of a molecule through a conical intersection creates a coherent superposition of electronic states via nonadiabatic couplings. Detecting this coherent superposition may serve as a direct probe of the conical intersection. In this paper, we theoretically demonstrate the use of coherent spontaneous emission in samples with long-range order for probing the occurrence of a conical intersection in a molecule. Our simulations show that the spectrum contains clear signatures of the created coherent superposition of electronic states. We investigate the bandwidth requirements for the x-ray probes, which influence the observation of coherent superposition generated by the conical intersection.

Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra

Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.

A Quantum Algorithm from Response Theory: Digital Quantum Simulation of Two-Dimensional Electronic Spectroscopy

Multidimensional optical spectroscopies are powerful techniques to investigate energy transfer pathways in natural and artificial systems. Because of the high information content of the spectra, numerical simulations of the optical response are of primary importance to assist the interpretation of spectral features. However, the increasing complexity of the investigated systems and their quantum dynamics call for the development of novel simulation strategies. In this work, we consider using digital quantum computers. By combining quantum dynamical simulation and nonlinear response theory, we present a quantum algorithm for computing the optical response of molecular systems. The quantum advantage stems from the efficient quantum simulation of the dynamics governed by the molecular Hamiltonian, and it is demonstrated by explicitly considering exciton-vibrational coupling. The protocol is tested on a near-term quantum device, providing the digital quantum simulation of the linear and nonlinear response of simple molecular models.

Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study

This study leverages two-pulse femtosecond stimulated Raman spectroscopy (2FSRS) to characterize molecular systems with avoided crossings (ACs) and conical intersections (CIs) in their low-lying excited electronic states. By simulating 2FSRS spectra of microscopically inspired ACs and CIs models, we demonstrate that 2FSRS not only delivers valuable information on the molecular parameters characterizing ACs and CIs but also helps distinguish between these two systems.

Tuning ultrafast time-evolution of photo-induced charge-transfer states: A real-time electronic dynamics study in substituted indenotetracene derivatives

Photo-induced charge transfer (CT) states are pivotal in many technological and biological processes. A deeper knowledge of such states is mandatory for modeling the charge migration dynamics. Real-time time-dependent density functional theory (RT-TD-DFT) electronic dynamics simulations are employed to explicitly observe the electronic density time-evolution upon photo-excitation. Asymmetrically substituted indenotetracene molecules, given their potential application as n-type semiconductors in organic photovoltaic materials, are here investigated. Effects of substituents with different electron-donating characters are analyzed in terms of the overall electronic energy spacing and resulting ultrafast CT dynamics through linear response (LR-)TD-DFT and RT-TD-DFT based approaches. The combination of the computational techniques here employed provided direct access to the electronic density reorganization in time and to its spatial and rational representation in terms of molecular orbital occupation time evolution. Such results can be exploited to design peculiar directional charge dynamics, crucial when photoactive materials are used for light-harvesting applications.

Watching the Interplay between Photoinduced Ultrafast Charge Dynamics and Nuclear Vibrations

Here is presented the ultrafast hole-electron dynamics of photoinduced metal to ligand charge-transfer (MLCT) states in a Ru(II) complex, [Ru(dcbpy)2(NCS)2]4- (dcbpy = 4,4'-dicarboxy-2,2'-bipyridine), a photoactive molecule employed in dye sensitized solar cells. Via cutting-edge computational techniques, a tailored computational protocol is here presented and developed to provide a detailed analysis of the electronic manifold coupled with nuclear vibrations to better understand the nonradiative pathways and the resulting overall dye performances in light-harvesting processes (electron injection). Thus, the effects of different vibrational modes were investigated on both the electronic levels and charge transfer dynamics through a theoretical-computational approach. First, the linear response time-dependent density functional (LR-TDDFT) formalism was employed to characterize excitation energies and spacing among electronic levels (the electronic layouts). Then, to understand the ultrafast (femtosecond) charge dynamics on the molecular scale, we relied on the nonperturbative mean-field quantum electronic dynamics via real-time (RT-) TDDFT. Three vibrational modes were selected, representative for collective nuclear movements that can have a significant influence on the electronic structure: two involving NCS- ligands and one involving dcbpy ligands. As main results, we observed that such MLCT states, under vibrational distortions, are strongly affected and a faster interligand electron transfer mechanism is observed along with an increasing MLCT character of the adiabatic electronic states approaching closer in energy due to the vibrations. Such findings can help both in providing a molecular picture of multidimensional vibro-electronic spectroscopic techniques, used to characterize ultrafast coherent and noncoherent dynamics of complex systems, and to improve dye performances with particular attention to the study of energy or charge transport processes and vibronic couplings.

Active Q-Switched X-Ray Regenerative Amplifier Free-Electron Laser

Despite tremendous progress in x-ray free-electron laser (FEL) science over the last decade, future applications still demand fully coherent, stable x rays that have not been demonstrated in existing x-ray FEL facilities. In this Letter, we describe an active Q-switched x-ray regenerative amplifier FEL scheme to produce fully coherent, high-brightness, hard x rays at a high-repetition rate. By using simple electron-beam phase space manipulation, we show this scheme is flexible in controlling the x-ray cavity quality factor Q and hence the output radiation. We report both theoretical and numerical studies on this scheme with a wide range of accelerator, x-ray cavity, and undulator parameters.

Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore

By employing the doorway-window (DW) on-the-fly simulation protocol, we performed ab initio simulations of peak evolutions and beating maps of electronic two-dimensional (2D) spectra of a polyatomic molecule in the gas phase. As the system under study, we chose pyrazine, which is a paradigmatic example of photodynamics dominated by conical intersections (CIs). From the technical perspective, we demonstrate that the DW protocol is a numerically efficient methodology suitable for simulations of 2D spectra for a wide range of excitation/detection frequencies and population times. From the information content perspective, we show that peak evolutions and beating maps not only reveal timescales of transitions through CIs but also pinpoint the most relevant coupling and tuning modes active at these CIs.

X-ray/Extreme Ultraviolet Floquet State Multidimensional Spectroscopy, an Analogue of Multiple Quantum Nuclear Magnetic Resonance

There is great interest in developing fully coherent multidimensional X-ray/extreme ultraviolet (XUV) spectroscopic techniques because of their capability for achieving atomic spectral selectivity. Current proposals rest on using sequentially and coherently driven core excitations with multiple X-ray/XUV excitation pulses and measuring the output using time domain Fourier transform methods. In this paper, we propose an alternative method that creates an entanglement of core and optical transitions to form a Floquet state that creates directional and coherent output beams. Multidimensional spectra are obtained by measuring the intensity of output beams while tuning the optical frequencies across resonances. This approach expands on previous optical pump-XUV probe spectroscopy of MoTe₂ by theoretically demonstrating its multidimensional capabilities. Both parametric and non-parametric pathways are proposed to optimize the resolution of inhomogeneous broadening and k-selective features.

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.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Quantum-chemical calculation of two-dimensional infrared spectra using localized-mode VSCF/VCI

Computational protocols for the simulation of two-dimensional infrared (2D IR) spectroscopy usually rely on vibrational exciton models which require an empirical parameterization. Here, we present an efficient quantum-chemical protocol for predicting static 2D IR spectra that does not require any empirical parameters. For the calculation of anharmonic vibrational energy levels and transition dipole moments, we employ the localized-mode vibrational self-consistent field (L-VSCF)/vibrational configuration interaction (L-VCI) approach previously established for (linear) anharmonic theoretical vibrational spectroscopy [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365-3377 (2014)]. We demonstrate that with an efficient expansion of the potential energy surface using anharmonic one-mode potentials and harmonic two-mode potentials, 2D IR spectra of metal carbonyl complexes and dipeptides can be predicted reliably. We further show how the close connection between L-VCI and vibrational exciton models can be exploited to extract the parameters of such models from those calculations. This provides a novel route to the fully quantum-chemical parameterization of vibrational exciton models for predicting 2D IR spectra.

Generalized Energy-Based Fragmentation Approach for the Electronic Emission Spectra of Large Systems

The excited-state (ES) geometry optimization and electronic emission (fluorescence and phosphorescence) spectra and the ES vibrational spectra of large systems are great challenges in quantum chemistry. In this work, we develop a generalized energy-based fragmentation (GEBF) approach to compute the localized ES structures and vibrational frequencies of large systems. In this approach, the ES energy derivatives (gradients or Hessians) for a localized ES of a large system can be obtained by combining the ES energy derivatives of the corresponding active subsystems (including local excitation center) and the ground-state energy derivatives of inactive subsystems. Two strategies are adopted to overcome two difficulties from state-classification and state-tracking for treating specific ESs. First, for state-classification, we develop an improved density-based spatial clustering applied with noise algorithm with a modified transition orbital projection (TOP) algorithm, which allow a certain ES energy and energy derivatives of the whole system to be calculated with different ES energies and energy derivatives of active subsystems. Furthermore, we also employ the TOP algorithm for tracking the ESs in their geometry optimizations at the time-dependent density functional theory (TDDFT) level. Then, the GEBF approach is applied to investigate the optimized ES geometries or ES vibrational frequencies for two typical systems. Our results show that the cost-effective GEBF approach can accurately reproduce the TDDFT fluorescence spectra of the cytosine derivative and the experimental phosphorescence spectra of the β-cyclodextrin derivative. The GEBF approach is expected to be routinely applied to investigate the electronic emission spectra of very large systems with local chromophores.

Time-resolved X-ray and XUV based spectroscopic methods for nonadiabatic processes in photochemistry

The photochemistry of numerous molecular systems is influenced by conical intersections (CIs). These omnipresent nonadiabatic phenomena provide ultra-fast radiationless relaxation channels by creating degeneracies between electronic states and decide over the final photoproducts. In their presence, the Born-Oppenheimer approximation breaks down, and the timescales of the electron and nuclear dynamics become comparable. Due to the ultra-fast dynamics and the complex interplay between nuclear and electronic degrees of freedom, the direct experimental observation of nonadiabatic processes close to CIs remains challenging. In this article, we give a theoretical perspective on novel spectroscopic techniques capable of observing clear signatures of CIs. We discuss methods that are based on ultra-short laser pulses in the extreme ultraviolet and X-ray regime, as their spectral and temporal resolution allow for resolving the ultra-fast dynamics near CIs.

Accurate Simulation of Spectroscopic Signatures of Cavity-Assisted, Conical-Intersection-Controlled Singlet Fission Processes

A numerically accurate, fully quantum methodology has been developed for the simulation of the dynamics and nonlinear spectroscopic signals of cavity-assisted, conical-intersection-controlled singlet fission systems. The methodology is capable of handling several molecular systems strongly coupled to the photonic mode of the cavity and treats the intrinsic conical intersection and cavity-induced polaritonic conical intersections in a numerically exact manner. Contributions of higher-lying molecular electronic states are accounted for comprehensively. The intriguing process of cavity-modified fission dynamics, including all of its electronic, vibrational, and photonic degrees of freedom, together with its two-dimensional spectroscopic manifestation, is simulated for two rubrene dimers strongly coupled to the cavity mode.

Probing nonadiabatic dynamics with attosecond pulse trains and soft x-ray Raman spectroscopy

Linear off-resonant x-ray Raman techniques are capable of detecting the ultrafast electronic coherences generated when a photoexcited wave packet passes through a conical intersection. A hybrid femtosecond or attosecond probe pulse is employed to excite the system and stimulate the emission of the signal photon, where both fields are components of a hybrid pulse scheme. In this paper, we investigate how attosecond pulse trains, as provided by high-harmonic generation processes, perform as probe pulses in the framework of this spectroscopic technique, instead of single Gaussian pulses. We explore different combination schemes for the probe pulse as well as the impact of parameters of the pulse trains on the signals. Furthermore, we show how Raman selection rules and symmetry consideration affect the spectroscopic signal, and we discuss the importance of vibrational contributions to the overall signal. We use two different model systems, representing molecules of different symmetries, and quantum dynamics simulations to study the difference in the spectra. The results suggest that such pulse trains are well suited to capture the key features associated with the electronic coherence.

Generation of intense phase-stable femtosecond hard X-ray pulse pairs

Coherent nonlinear spectroscopies and imaging in the X-ray domain provide direct insight into the coupled motions of electrons and nuclei with resolution on the electronic length scale and timescale. The experimental realization of such techniques will strongly benefit from access to intense, coherent pairs of femtosecond X-ray pulses. We have observed phase-stable X-ray pulse pairs containing more than 3 × 107 photons at 5.9 keV (2.1 Å) with ∼1 fs duration and 2 to 5 fs separation. The highly directional pulse pairs are manifested by interference fringes in the superfluorescent and seeded stimulated manganese Kα emission induced by an X-ray free-electron laser. The fringes constitute the time-frequency X-ray analog of Young’s double-slit interference, allowing for frequency domain X-ray measurements with attosecond time resolution.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Towards Understanding Excited-State Properties of Organic Molecules Using Time-Resolved Soft X-ray Absorption Spectroscopy

The extension of the pump-probe approach known from UV/VIS spectroscopy to very short wavelengths together with advanced simulation techniques allows a detailed analysis of excited-state dynamics in organic molecules or biomolecular structures on a nanosecond to femtosecond time level. Optical pump soft X-ray probe spectroscopy is a relatively new approach to detect and characterize optically dark states in organic molecules, exciton dynamics or transient ligand-to-metal charge transfer states. In this paper, we describe two experimental setups for transient soft X-ray absorption spectroscopy based on an LPP emitting picosecond and sub-nanosecond soft X-ray pulses in the photon energy range between 50 and 1500 eV. We apply these setups for near-edge X-ray absorption fine structure (NEXAFS) investigations of thin films of a metal-free porphyrin, an aggregate forming carbocyanine and a nickel oxide molecule. NEXAFS investigations have been carried out at the carbon, nitrogen and oxygen K-edge as well as on the Ni L-edge. From time-resolved NEXAFS carbon, K-edge measurements of the metal-free porphyrin first insights into a long-lived trap state are gained. Our findings are discussed and compared with density functional theory calculations.

Simulation of Condensed-Phase Spectroscopy with Near-Term Digital Quantum Computers

Spectroscopy is an indispensable tool for understanding the structures and dynamics of molecular systems. However, computational modeling of spectroscopy is challenging due to the exponential scaling of computational complexity with system sizes unless drastic approximations are made. Quantum computers could potentially overcome these classically intractable computational tasks, but the existing approaches using quantum computers to simulate spectroscopy can only handle isolated and static molecules. In this work, we develop a workflow that combines multi-scale modeling and a time-dependent variational quantum algorithm to compute the linear spectroscopy of systems interacting with their condensed-phase environment via the relevant time correlation function. We demonstrate the feasibility of our approach by numerically simulating the UV-vis absorption spectra of organic semiconductors. We show that our dynamical approach captures several spectral features that are otherwise overlooked by static methods. Our method can be directly used for other linear condensed-phase spectroscopy and could potentially be extended to nonlinear multi-dimensional spectroscopy.

Modelling Photoionisation in Isocytosine: Potential Formation of Longer-Lived Excited State Cations in its Keto Form

Studying the effects of UV and VUV radiation on non-canonical DNA/RNA nucleobases allows us to compare how they release excess energy following absorption with respect to their canonical counterparts. This has attracted much research attention in recent years because of its likely influence on the origin of our genetic lexicon in prebiotic times. Here we present a CASSCF and XMS-CASPT2 theoretical study of the photoionisation of non-canonical pyrimidine nucleobase isocytosine in both its keto and enol tautomeric forms. We analyse their lowest energy cationic excited states including 2 π + , 2 n O + and 2 n N + and compare these to the corresponding electronic states in cytosine. Investigating lower-energy decay pathways we find - unexpectedly - that keto-isocytosine+ presents a sizeable energy barrier potentially inhibiting decay to its cationic ground state, whereas enol-isocytosine+ features a barrierless and consequently ultrafast pathway analogous to the one previously found for the canonical (keto) form of cytosine+ . Dynamic electron correlation reduces the energy barrier in the keto form substantially (by ∼1 eV) but it is nevertheless still present. We additionally compute the UV/Vis absorption signals of the structures encountered along these decay channels to provide spectroscopic fingerprints to assist future experiments in monitoring these intricate photo-processes.

Simulation and Analysis of the Transient Absorption Spectrum of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Acetonitrile

4-(N,N-Dimethylamino)benzonitrile (DMABN) is a well-known model compound for dual fluorescence-in sufficiently polar solvents, it exhibits two distinct fluorescence emission bands. The interpretation of its transient absorption (TA) spectrum in the visible range is the subject of a long-standing controversy. In the present study, we resolve this issue by calculating the TA spectrum on the basis of nonadiabatic molecular dynamics simulations. An unambiguous assignment of spectral signals to specific excited-state structures is achieved by breaking down the calculated spectrum into contributions from twisted and nontwisted molecular geometries. In particular, the much-discussed excited-state absorption band near 1.7 eV (ca. 700 nm) is attributed to the near-planar locally excited (LE) minimum on the S1 state. On the technical side, our study demonstrates that the second-order approximate coupled cluster singles and doubles (CC2) method can be used successfully to calculate the TA spectra of moderately large organic molecules, provided that the system in question does not approach a crossing between the lowest excited state and the singlet ground state within the time frame of the simulation.

Timing of charge migration in betaine by impact of fast atomic ions

The way molecules break after ion bombardment is intimately related to the early electron dynamics generated in the system, in particular, charge (or electron) migration. We exploit the natural positive-negative charge splitting in the zwitterionic molecule betaine to selectively induce double electron removal from its negatively charged side by impact of fast O6+ ions. The loss of two electrons in this localized region of the molecular skeleton triggers a competition between direct Coulomb explosion and charge migration that is examined to obtain temporal information from ion-ion coincident measurements and nonadiabatic molecular dynamics calculations. We find a charge migration time, from one end of the molecule to the other, of approximately 20 to 40 femtoseconds. This migration time is longer than that observed in molecules irradiated by ultrashort light pulses and is the consequence of charge migration being driven by adiabatic nuclear dynamics in the ground state of the molecular dication.

Time-Resolved Photoelectron Spectroscopy of Conical Intersections with Attosecond Pulse Trains

Conical Intersections (CIs), which are believed to be ubiquitous in molecular and biological systems, open up ultrafast nonradiative decay channels. A superposition of electronic states is created when a molecule passes through a CI and the nuclear wave packet branches. The resulting electronic coherence can be considered a unique signature of the CI. The involved electronic states can be resolved in the energy domain with photoelectron spectroscopy using a femtosecond pulse as a probe. However, the observation of the created electronic coherence in the time domain requires probe pulses with several electron volts of bandwidth. Attosecond pulses can probe the electronic coherence but are unable to resolve the involved electronic states. In this Letter, we propose to address this restriction by using time-resolved photoelectron spectroscopy with an attosecond pulse train as a probe. We theoretically demonstrate that the resulting photoelectron spectrum may yield energy resolution as well as the information on the created coherences in the time domain.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Polariton induced conical intersection and berry phase

We investigate the Polariton induced conical intersection (PICI) created from coupling a diatomic molecule with the quantized photon mode inside an optical cavity, and the corresponding Berry Phase effects. We use the rigorous Pauli-Fierz Hamiltonian to describe the quantum light-matter interactions between a LiF molecule and the cavity, and use the exact quantum propagation to investigate the polariton quantum dynamics. The molecular rotations relative to the cavity polarization direction play a role as the tuning mode of the PICI, resulting in an effective CI even within a diatomic molecule. To clearly demonstrate the dynamical effects of the Berry phase, we construct two additional models that have the same Born-Oppenheimer surface, but the effects of the geometric phase are removed. We find that when the initial wavefunction is placed in the lower polaritonic surface, the Berry phase causes a π phase-shift in the wavefunction after the encirclement around the CI, indicated from the nuclear probability distribution. On the other hand, when the initial wavefunction is placed in the upper polaritonic surface, the geometric phase significantly influences the couplings between polaritonic states and therefore, the population dynamics between them. These BP effects are further demonstrated through the photo-fragment angular distribution. PICI created from the quantized radiation field has the promise to open up new possibilities to modulate photochemical reactivities.

Site-specific generation of excited state wavepackets with high-intensity attosecond x rays

High-intensity attosecond x rays can produce coherent superpositions of valence-excited states through two-photon Raman transitions. The broad-bandwidth, high-field nature of the pulses results in a multitude of accessible excited states. Multiconfigurational quantum chemistry with the time-dependent Schrödinger equation is used to examine population transfer dynamics in stimulated x-ray Raman scattering of the nitric oxide oxygen and nitrogen K-edges. Two pulse schemes initiate wavepackets of different characters and demonstrate how chemical differences between core-excitation pathways affect the dynamics. The population transfer to valence-excited states is found to be sensitive to the electronic structure and pulse conditions, highlighting complexities attributed to the Rabi frequency. The orthogonally polarized two-color-pulse setup has increased selectivity while facilitating longer, less intense pulses than the one-pulse setup. Population transfer in the 1s → Rydberg region is more effective but less selective at the nitrogen K-edge; the selectivity is reduced by double core-excited states. Result interpretation is aided by resonant inelastic x-ray scattering maps.

Capturing fingerprints of conical intersection: Complementary information of non-adiabatic dynamics from linear x-ray probes

Many recent experimental ultrafast spectroscopy studies have hinted at non-adiabatic dynamics indicating the existence of conical intersections, but their direct observation remains a challenge. The rapid change of the energy gap between the electronic states complicated their observation by requiring bandwidths of several electron volts. In this manuscript, we propose to use the combined information of different x-ray pump-probe techniques to identify the conical intersection. We theoretically study the conical intersection in pyrrole using transient x-ray absorption, time-resolved x-ray spontaneous emission, and linear off-resonant Raman spectroscopy to gather evidence of the curve crossing.

Non-linear self-driven spectral tuning of Extreme Ultraviolet Femtosecond Pulses in monoatomic materials

Self-action nonlinearity is a key aspect - either as a foundational element or a detrimental factor - of several optical spectroscopies and photonic devices. Supercontinuum generation, wavelength converters, and chirped pulse amplification are just a few examples. The recent advent of Free Electron Lasers (FEL) fostered building on nonlinearity to propose new concepts and extend optical wavelengths paradigms for extreme ultraviolet (EUV) and X-ray regimes. No evidence for intrapulse dynamics, however, has been reported at such short wavelengths, where the light-matter interactions are ruled by the sharp absorption edges of core electrons. Here, we provide experimental evidence for self-phase modulation of femtosecond FEL pulses, which we exploit for fine self-driven spectral tunability by interaction with sub-micrometric foils of selected monoatomic materials. Moving the pulse wavelength across the absorption edge, the spectral profile changes from a non-linear spectral blue-shift to a red-shifted broadening. These findings are rationalized accounting for ultrafast ionization and delayed thermal response of highly excited electrons above and below threshold, respectively.

Selective Enhancement of Spectroscopic Features by Quantum Optimal Control

Tailored light can be used to steer atomic motions into selected quantum pathways. In optimal control theory (OCT), the target is usually expressed in terms of the molecular wave function, a quantity that is not directly observable in experiment. We present simulations using OCT that optimize the spectroscopic signal itself. By shaping the optical pump, the x-ray stimulated Raman signal, which occurs solely during the passage through conical intersections, is temporally controlled and amplified by up to 2 orders of magnitude. This enhancement can be crucial in order to bring small coherence-based signatures above the detectable threshold. Our approach is applicable to any signal that depends on the expectation value of a positive definite operator.

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

Monitoring molecular vibronic coherences in a bichromophoric molecule by ultrafast X-ray spectroscopy

The role of quantum-mechanical coherences in the elementary photophysics of functional optoelectronic molecular materials is currently under active study. Designing and controlling stable coherences arising from concerted vibronic dynamics in organic chromophores is the key for numerous applications. Here, we present fundamental insight into the energy transfer properties of a rigid synthetic heterodimer that has been experimentally engineered to study coherences. Quantum non-adiabatic excited state simulations are used to compute X-ray Raman signals, which are able to sensitively monitor the coherence evolution. Our results verify their vibronic nature, that survives multiple conical intersection passages for several hundred femtoseconds at room temperature. Despite the contributions of highly heterogeneous evolution pathways, the coherences are unambiguously visualized by the experimentally accessible X-ray signals. They offer direct information on the dynamics of electronic and structural degrees of freedom, paving the way for detailed coherence measurements in functional organic materials.

Multi-wave mixing in the high harmonic regime: monitoring electronic dynamics

It has been demonstrated that electronic coherences across many eV can be detected in pump-probe experiments involving high harmonic sources. An additional degree of control over the phase matching can be employed by investigating a more general class of multi-wave mixing. Non-collinear multi-wave mixing of high harmonics with energy (q₁ω₁ + q₂ω₂) can be selectively detected along the direction of (q₁k₁ + q₂k₂). Simulations based on a recently developed semi-perturbative approach show that only the specific harmonic signals with q₁ω₁ close to the energy difference between ground state and excited states are observable when the two input pulses are well separated in time. The coherent dynamics between different states can be selectively tracked by detecting the time-delay dependent signals with different q₁k₁, which can overcome the potential spectral congestion in real experiments. Additionally, such non-collinear geometry can be used to separate the dephasing induced decay and collision induced recovery behaviors of pump-probe high harmonic signal typically observed in the time-resolved high harmonic pump-probe signals.

Manipulating Core Excitations in Molecules by X-Ray Cavities

Core excitations on different atoms are highly localized and therefore decoupled. By placing molecules in an x-ray cavity the core transitions become coupled via the exchange of cavity photons and form delocalized hybrid light-matter excitations known as core polaritons. We demonstrate these effects for the two inequivalent carbon atoms in 1,1-difluoroethylene. Polariton signatures in the x-ray absorption, two-photon absorption, and multidimensional four-wave mixing signals are predicted.

Intermolecular conical intersections in molecular aggregates

Conical intersections (CoIns) of multidimensional potential energy surfaces are ubiquitous in nature and control pathways and yields of many photo-initiated intramolecular processes. Such topologies can be potentially involved in the energy transport in aggregated molecules or polymers but are yet to be uncovered. Here, using ultrafast two-dimensional electronic spectroscopy (2DES), we reveal the existence of intermolecular CoIns in molecular aggregates relevant for photovoltaics. Ultrafast, sub-10-fs 2DES tracks the coherent motion of a vibrational wave packet on an optically bright state and its abrupt transition into a dark state via a CoIn after only 40 fs. Non-adiabatic dynamics simulations identify an intermolecular CoIn as the source of these unusual dynamics. Our results indicate that intermolecular CoIns may effectively steer energy pathways in functional nanostructures for optoelectronics.

On the role of non-diagonal system-environment interactions in bridge-mediated electron transfer

Bridge-mediated electron transfer (ET) between a donor and an acceptor is prototypical for the description of numerous most important ET scenarios. While multi-step ET and the interplay of sequential and direct superexchange transfer pathways in the donor-bridge-acceptor (D-B-A) model are increasingly understood, the influence of off-diagonal system-bath interactions on the transfer dynamics is less explored. Off-diagonal interactions account for the dependence of the ET coupling elements on nuclear coordinates (non-Condon effects) and are typically neglected. Here, we numerically investigate with quasi-adiabatic propagator path integral simulations the impact of off-diagonal system-environment interactions on the transfer dynamics for a wide range of scenarios in the D-B-A model. We demonstrate that off-diagonal system-environment interactions can have profound impact on the bridge-mediated ET dynamics. In the considered scenarios, the dynamics itself does not allow for a rigorous assignment of the underlying transfer mechanism. Furthermore, we demonstrate how off-diagonal system-environment interaction mediates anomalous localization by preventing long-time depopulation of the bridge B and how coherent transfer dynamics between donor D and acceptor A can be facilitated. The arising non-exponential short-time dynamics and coherent oscillations are interpreted within an equivalent Hamiltonian representation of a primary reaction coordinate model that reveals how the complex vibronic interplay of vibrational and electronic degrees of freedom underlying the non-Condon effects can impose donor-to-acceptor coherence transfer on short timescales.