Probing ultrafast ππ*/nπ* internal conversion in organic chromophores via K-edge resonant absorption

The Role of Momentum Partitioning in Covariance Ion Imaging Analysis

We present results from a covariance ion imaging study, which employs extensive filtering, on the relationship between fragment momenta to gain deeper insight into photofragmentation dynamics. A new data analysis approach is introduced that considers the momentum partitioning between the fragments of the breakup of a molecular polycation to disentangle concurrent fragmentation channels, which yield the same ion species. We exploit this approach to examine the momentum exchange relationship between the products, which provides direct insight into the dynamics of molecular fragmentation. We apply these techniques to extensively characterize the dissociation of 1-iodopropane and 2-iodopropane dications prepared by site-selective ionization of the iodine atom using extreme ultraviolet intense femtosecond laser pulses with a photon energy of 95 eV. Our assignments are supported by classical simulations, using parameters largely obtained directly from the experimental data.

A variational reformulation of molecular properties in electronic-structure theory

Conventional quantum-mechanical calculations of molecular properties, such as dipole moments and electronic excitation energies, give errors that depend linearly on the error in the wave function. An exception is the electronic energy, whose error depends quadratically on the error in wave function. We here describe how all properties may be calculated with a quadratic error, by setting up a variational Lagrangian for the property of interest. Because the construction of the Lagrangian is less expensive than the calculation of the wave function, this approach substantially improves the accuracy of quantum-chemical calculations without increasing cost. As illustrated for excitation energies, this approach enables the accurate calculation of molecular properties for larger systems, with a short time-to-solution and in a manner well suited for modern computer architectures.

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.

Probing Conical Intersection in the Multipathway Isomerization of CH3Cl Using Coulomb Explosion

Ubiquitous ultrafast isomerization is paramount in photoexcited molecules, in which non-adiabatic coupling among multiple electronic states can occur. We use the pump-probe Coulomb explosion imaging method to study the isomerization of CH3Cl molecules. We find that the isomerization under our strong field pump-probe scheme proceeds along multiple pathways, which are encoded in several distinct branches of the time-resolved kinetic energy release spectra for the CH2++HCl+ Coulomb explosion channel. Apart from the isomerized dissociative pathway in neutral and cationic excited states, the pump laser can also induce coherent vibrational dynamics in two coupled intermediate states and set up the initial conditions for the two concurrently proceeding isomerization pathways. The isomerization of CH3Cl provides an intriguing example of a chemical reaction consisting of multiple pathways and non-adiabatic dynamics.

Raman Spectroscopy of Conical Intersections Using Entangled Photons

Ultrafast Raman spectroscopy with attosecond pulses in the extreme ultraviolet and X-ray regime has been proposed theoretically for tracking the non-adiabatic dynamics of molecules in great detail. The large bandwidth of these pulses, which span several electronvolts within a couple of femtoseconds, provides a unique tool for tracking non-adiabatic phenomena. However, spectroscopy with classical light is limited by the time-bandwidth product of the probe laser pulse. In this work, we theoretically investigate an ultrafast Raman spectroscopy scheme that utilizes pairs of entangled photons. Our model simulations demonstrate that the dynamics in the vicinity of a conical intersection can be resolved with unprecedented resolution in the time and frequency domain.

Spin Coupling Effect on Geometry-Dependent X-Ray Absorption of Diradicals

We theoretically investigate the influence of diradical electron spin coupling on the time-resolved X-ray absorption spectra of the photochemical ring opening of furanone. We predict geometry-dependent carbon K-edge signals involving transitions from core orbitals to both singly and unoccupied molecular orbitals. The most obvious features of the ring opening come from the carbon atom directly involved in the bond breaking through its transition to both the newly formed singly occupied and the available lowest unoccupied molecular orbitals (SOMO and LUMO, respectively). In addition to this primary feature, the singlet spin coupling of four unpaired electrons that arises in the core-to-LUMO states creates additional geometry dependence in some spectral features with both oscillator strengths and relative excitation energies varying observably as a function of the ring opening. We attribute this behavior to a spin-occupancy-induced selection rule, which occurs when singlet spin coupling is enforced in the diradical state. Notably, one of these geometry-sensitive core-to-LUMO transitions excites core electrons from a backbone carbon not involved in the bond breaking, providing a novel nonlocal X-ray probe of chemical dynamics arising from electron spin coupling.

Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone

We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.

Toward ultrafast soft x-ray spectroscopy of organic photovoltaic devices

Novel ultrafast x-ray sources based on high harmonic generation and at x-ray free electron lasers are opening up new opportunities to resolve complex ultrafast processes in condensed phase systems with exceptional temporal resolution and atomic site specificity. In this perspective, we present techniques for resolving charge localization, transfer, and separation processes in organic semiconductors and organic photovoltaic devices with time-resolved soft x-ray spectroscopy. We review recent results in ultrafast soft x-ray spectroscopy of these systems and discuss routes to overcome the technical challenges in performing time-resolved x-ray experiments on photosensitive materials with poor thermal conductivity and low pump intensity thresholds for nonlinear effects.

Characterizing Conical Intersections in DNA/RNA Nucleobases with Multiconfigurational Wave Functions of Varying Active Space Size

We characterize the photochemically relevant conical intersections between the lowest-lying accessible electronic excited states of the different DNA/RNA nucleobases using Cholesky decomposition-based complete active space self-consistent field (CASSCF) algorithms. We benchmark two different basis set contractions and several active spaces for each nucleobase and conical intersection type, measuring for the first time how active space size affects conical intersection topographies in these systems and the potential implications these may have toward their description of photoinduced phenomena. Our results show that conical intersection topographies are highly sensitive to the electron correlation included in the model: by changing the amount (and type) of correlated orbitals, conical intersection topographies vastly change, and the changes observed do not follow any converging pattern toward the topographies obtained with the largest and most correlated active spaces. Comparison across systems shows analogous topographies for almost all intersections mediating population transfer to the dark 1nO/Nπ* states, while no similarities are observed for the "ethylene-like" conical intersection ascribed to mediate the ultrafast decay component to the ground state in all DNA/RNA nucleobases. Basis set size seems to have a minor effect, appearing to be relevant only for purine-based derivatives. We rule out structural changes as a key factor in classifying the different conical intersections, which display almost identical geometries across active space and basis set change, and we highlight instead the importance of correctly describing the electronic states involved at these crossing points. Our work shows that careful active space selection is essential to accurately describe conical intersection topographies and therefore to adequately account for their active role in molecular photochemistry.

Following the Nonadiabatic Ultrafast Dynamics of Uracil via Simulated X-ray Absorption Spectra

The nucleobase uracil exhibits high photostability due to ultrafast relaxation processes mediated by conical intersections (CoIns), where the interplay between nuclear and electron dynamics becomes crucial. In our previous study, we observed seemingly long-lived traces of electronic coherence for the relaxation process through the S2/S1 CoIn by applying our ansatz for coupled nuclear and electron dynamics in molecules (NEMol). In this work, we theoretically investigate how time-dependent transient X-ray absorption spectroscopy can be used to observe this ultrafast dynamics. Therefore, we calculated X-ray absorption spectra (XAS) for the oxygen K-edge, using a multireference protocol in combination with NEMol dynamics. Thus, we have access to both the transient XAS based on the nuclear wavepacket dynamics and the modulation of the signals caused by the electronic coherence induced by the excitation process and the presence of a CoIn seam. In both cases, we were able to qualitatively predict its influence on the resulting XAS.

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

High-Throughput Screening for Ultrafast Photochemical Reaction Discovery

High-repetition-rate lasers present an opportunity to extend ultrafast spectroscopy from a detailed probe of singular model photochemical systems to a routine analysis technique in training machine learning models to aid the design cycle of photochemical syntheses. We bring together innovations in line scan cameras and micro-electro-mechanical grating modulators with sample delivery via high-pressure liquid chromatography pumps to demonstrate a transient absorption spectrometer that can characterize photoreactions initiated with ultrashort ultraviolet pulses in a time scale of minutes. Furthermore, we demonstrate that the ability to rapidly screen an important class of photochemical system, pyrimidine nucleosides, can be used to explore the effect of conformational modification on the evolution of excited-state processes.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory

The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.

Tuning UV Pump X-ray Probe Spectroscopy on the Nitrogen K Edge Reveals the Radiationless Relaxation of Pyrazine: Ab Initio Simulations Using the Quasiclassical Doorway-Window Approximation

Transient absorption UV pump X-ray probe spectroscopy has been established as a versatile technique for the exploration of ultrafast photoinduced dynamics in valence-excited states. In this work, an ab initio theoretical framework for the simulation of time-resolved UV pump X-ray probe spectra is presented. The method is based on the description of the radiation-matter interaction in the classical doorway-window approximation and a surface-hopping algorithm for the nonadiabatic nuclear excited-state dynamics. Using the second-order algebraic-diagrammatic construction scheme for excited states, UV pump X-ray probe signals were simulated for the carbon and nitrogen K edges of pyrazine, assuming a duration of 5 fs of the UV pump and X-ray probe pulses. It is predicted that spectra measured at the nitrogen K edge carry much richer information about the ultrafast nonadiabatic dynamics in the valence-excited states of pyrazine than those measured at the carbon K edge.

A theoretical study of the time-resolved x-ray absorption spectrum of the photoionized BT-1T cation

The time-resolved x-ray absorption spectrum of the BT-1T cation (BT-1T+) is theoretically simulated in order to investigate the charge transfer reaction of the system. We employ both trajectory surface hopping and quantum dynamics to simulate the structural evolution over time and the changes in the state populations. To compute the static x-ray absorption spectra (XAS) of the ground and excited states, we apply both the time-dependent density functional theory and the coupled cluster singles and doubles method. The results obtained are in good agreement between the methods. It is, furthermore, found that the small structural changes that occur during the reaction have little effect on the static XAS. Hence, the tr-XAS can be computed based on the state populations determined from a nuclear dynamics simulation and one set of static XAS calculations, utilizing the ground state optimized geometry. This approach can save considerable computational resources, as the static spectra need not to be calculated for all geometries. As BT-1T is a relatively rigid molecule, the outlined approach should only be considered when investigating non-radiative decay processes in the vicinity of the Franck-Condon point.

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.

Attosecond Monitoring of Nonadiabatic Molecular Dynamics by Transient X-ray Transmission Spectroscopy

Tracing the evolution of molecular coherences can provide a direct, unambiguous probe of nonadiabatic molecular processes, such as the passage through conical intersections of electronic states. Two techniques, attosecond transient absorption spectroscopy (ATAS) and Transient Redistribution of Ultrafast Electronic Coherences in Attosecond Raman Signals (TRUECARS), have been used or proposed for monitoring nonadiabatic molecular dynamics. Both techniques employ the transmission of a weak attosecond extreme-ultraviolet or X-ray probe to interrogate the molecule at controllable time delays with respect to an optical pump, thereby extracting dynamical information from transient spectral features. The connection between these techniques has not been firmly established yet. In this theoretical study, we provide a unified description of both transient transmission techniques, establishing their relationship as limits of the same pump-probe spectroscopy technique for different pulse parameter regimes. We demonstrate this by quantum dynamical simulations of thiophenol photodissociation and show how complementary coherence information can be revealed by the two techniques.

Ultrafast Electronic Relaxation in 6-Methyluracil and 5-Fluorouracil in Isolated and Aqueous Conditions: Substituent and Solvent Effects

We report ultrafast extreme ultraviolet photoelectron spectroscopy of 6-methyluracil (6mUra) and 5-fluorouracil (5FUra) in the gas phase and 6mUra and 5-fluorouridine in an aqueous environment. In the gas phase, internal conversion (IC) occurs from ¹ππ* to ¹nπ* states in tens of femtoseconds, followed by intersystem crossing to the ³ππ* state in several picoseconds. In an aqueous solution, 6mUra undergoes IC almost exclusively to the ground state (S₀) in about 100 fs, which is essentially the same process as that for unsubstituted uracil, but much faster than that for thymine (5-methyluracil). The different dynamics for C5 and C6 methylation suggest that IC from ¹ππ* to S₀ is facilitated by out-of-plane (OOP) motion of the C5 substituent. The slow IC for C5-substituted molecules in an aqueous environment is ascribed to the solvent reorganization that is required for this OOP motion to occur. The slow rate for 5FUrd may arise in part from an increased barrier height due to C5 fluorination.

Photophysics of uracil: an explicit time-dependent generating function-based method combining both nonadiabatic and spin-orbit coupling effects

We present a composite framework for calculating the rates of non-radiative deactivation processes, namely internal conversion (IC) and intersystem crossing (ISC), on an equal footing by explicitly computing the non-adiabatic coupling (NAC) and spin-orbit coupling (SOC) constants, respectively. The stationary-state approach uses a time-dependent generating function based on Fermi's golden rule. We validate the applicability of the framework by computing the rate of IC for azulene, obtaining comparable rates to experimental and previous theoretical results. Next, we investigate the photophysics associated with the complex photodynamics of the uracil molecule. Interestingly, our simulated rates corroborate experimental observations. Detailed analyses using Duschinsky rotation matrices, displacement vectors and NAC matrix elements are presented to interpret the findings alongside testing the suitability of the approach for such molecular systems. The suitability of the Fermi's golden rule based method is explained qualitatively in terms of single-mode potential energy surfaces.

Formation of Long-Lived Dark States during Electronic Relaxation of Pyrimidine Nucleobases Studied Using Extreme Ultraviolet Time-Resolved Photoelectron Spectroscopy

Ultrafast electronic relaxation of nucleobases from ¹ππ* states to the ground state (S₀) is considered essential for the photostability of DNA. However, transient absorption spectroscopy (TAS) has indicated that some nucleobases in aqueous solutions create long-lived ¹nπ*/³ππ* dark states from the ¹ππ* states with a high quantum yield of 0.4-0.5. We investigated electronic relaxation in pyrimidine nucleobases in both aqueous solutions and the gas phase using extreme ultraviolet (EUV) time-resolved photoelectron spectroscopy. Femtosecond EUV probe pulses cause ionization from all electronic states involved in the relaxation process, providing a clear overview of the electronic dynamics. The ¹nπ* quantum yields for aqueous cytidine and uracil (Ura) derivatives were found to be considerably lower (<0.07) than previous estimates reported by TAS. On the other hand, aqueous thymine (Thy) and thymidine exhibited a longer ¹ππ* lifetime and a higher quantum yield (0.12-0.22) for the ¹nπ* state. A similar trend was found for isolated Thy and Ura in the gas phase: the ¹ππ* lifetimes are 39 and 17 fs and the quantum yield for ¹nπ* are 1.0 and 0.45 for Thy and Ura, respectively. The result indicates that single methylation to the C₅ position hinders the out-of-plane deformation that drives the system to the conical intersection region between ¹ππ* and S₀, providing a large impact on the photophysics/photochemistry of a pyrimidine nucleobase. The significant reduction of ¹nπ* yield in aqueous solution is ascribed to the destabilization of the ¹nπ* state induced by hydrogen bonding.

Electron and ion spectroscopy of azobenzene in the valence and core shells

Azobenzene is a prototype and a building block of a class of molecules of extreme technological interest as molecular photo-switches. We present a joint experimental and theoretical study of its response to irradiation with light across the UV to x-ray spectrum. The study of valence and inner shell photo-ionization and excitation processes combined with measurement of valence photoelectron-photoion coincidence and mass spectra across the core thresholds provides a detailed insight into the site- and state-selected photo-induced processes. Photo-ionization and excitation measurements are interpreted via the multi-configurational restricted active space self-consistent field method corrected by second order perturbation theory. Using static modeling, we demonstrate that the carbon and nitrogen K edges of azobenzene are suitable candidates for exploring its photoinduced dynamics thanks to the transient signals appearing in background-free regions of the NEXAFS and XPS.

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the ¹ππ* (S₂) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the ¹nπ* (S₁) and ³ππ* (T₁) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576-16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fittingfitted (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) nonadiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed accelerated deactivation decay pathways, which is consistent with the experimental decay time of ∼400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory for Accurate X-ray Absorption Spectroscopy

It is demonstrated that the challenging core-hole particle (CHP) orbital relaxation for core electron spectra can be readily achieved by the mixed-reference spin-flip (MRSF)-time-dependent density functional theory (TDDFT). With the additional scalar relativistic effects on K-edge excitation energies of 24 second- and 17 third-row molecules, the particular ΔCHP-MRSF(R) exhibited near perfect predictions with RMSE ∼0.5 eV, featuring a median value of 0.3 and an interquartile range of 0.4. Overall, the CHP effect is 2-4 times stronger than relativistic ones, contributing more than 20 eV in the cases of sulfur and chlorine third-row atoms. Such high precision allows to explain the splitting and spectral shapes of O, N, and C atom K-edges in the ground state of thymine with atom as well as orbital specific accuracy. The same protocol with a double hole particle relaxation also produced remarkably accurate K-edge spectra of core to valence hole excitation energies from the first (n_O8π*) and second (ππ*) excited states of thymine, confirming the assignment of 1s → n excitation for the experimentally observed 526.4 eV peak. Regarding both accuracy and practicality, therefore, MRSF-TDDFT provides a promising protocol for core electron spectra of both ground and excited electronic states alike.

Pair Natural Orbital Equation-of-Motion Coupled-Cluster Method for Core Binding Energies: Theory, Implementation, and Benchmark

We present the theory and implementation of a lower scaling core-valence separated equation-of-motion coupled-cluster approach based on domain-based local pair natural orbitals for core binding energies. The accuracy of the new method has been compared with that of the standard equation-of-motion coupled-cluster method and experimentally measured results. The use of pair natural orbitals significantly reduces the computation cost and can be applied to large molecules.

The time-resolved atomic, molecular and optical science instrument at the Linac Coherent Light Source

The newly constructed time-resolved atomic, molecular and optical science instrument (TMO) is configured to take full advantage of both linear accelerators at SLAC National Accelerator Laboratory, the copper accelerator operating at a repetition rate of 120 Hz providing high per-pulse energy as well as the superconducting accelerator operating at a repetition rate of about 1 MHz providing high average intensity. Both accelerators power a soft X-ray free-electron laser with the new variable-gap undulator section. With this flexible light source, TMO supports many experimental techniques not previously available at LCLS and will have two X-ray beam focus spots in line. Thereby, TMO supports atomic, molecular and optical, strong-field and nonlinear science and will also host a designated new dynamic reaction microscope with a sub-micrometer X-ray focus spot. The flexible instrument design is optimized for studying ultrafast electronic and molecular phenomena and can take full advantage of the sub-femtosecond soft X-ray pulse generation program.

Hole Dynamics in Photoexcited Hematite Studied with Femtosecond Oxygen K-edge X-ray Absorption Spectroscopy

Hematite (α-Fe₂O₃) is a photoelectrode for the water splitting process because of its relatively narrow bandgap and abundance in the earth's crust. In this study, the photoexcited state of a hematite thin film was investigated with femtosecond oxygen K-edge X-ray absorption spectroscopy (XAS) at the PAL-XFEL in order to follow the dynamics of its photoexcited states. The 200 fs decay time of the hole state in the valence band was observed via its corresponding XAS feature.

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.

Time-Resolved Optical Pump-Resonant X-ray Probe Spectroscopy of 4-Thiouracil: A Simulation Study

We theoretically monitor the photoinduced ππ* → nπ* internal conversion process in 4-thiouracil (4TU), triggered by an optical pump. The element-sensitive spectroscopic signatures are recorded by a resonant X-ray probe tuned to the sulfur, oxygen, or nitrogen K-edge. We employ high-level electronic structure methods optimized for core-excited electronic structure calculation combined with quantum nuclear wavepacket dynamics computed on two relevant nuclear modes, fully accounting for their quantum nature of nuclear motions. We critically discuss the capabilities and limitations of the resonant technique. For sulfur and nitrogen, we document a pre-edge spectral window free from ground-state background and rich with ππ* and nπ* absorption features. The lowest sulfur K-edge shows strong absorption for both ππ* and nπ*. In the lowest nitrogen K-edge window, we resolve a state-specific fingerprint of the ππ* and an approximate timing of the conical intersection via its depletion. A spectral signature of the nπ* transition, not accessible by UV-vis spectroscopy, is identified. The oxygen K-edge is not sensitive to molecular deformations and gives steady transient absorption features without spectral dynamics. The ππ*/nπ* coherence information is masked by more intense contributions from populations. Altogether, element-specific time-resolved resonant X-ray spectroscopy provides a detailed picture of the electronic excited-state dynamics and therefore a sensitive window into the photophysics of thiobases.

Exciton Dynamics of a Diketo-Pyrrolopyrrole Core for All Low-Lying Electronic Excited States Using Density Functional Theory-Based Methods

Ab initio treatments of interexcited state internal conversion (IC) are more often than not missing from exciton dynamic descriptions, because of their inherent complexity. Here, we define "interexcited state IC" as a same-spin nonradiative transition between states i and j, where i ≠ j ≠ 0. Competing directly with multiexciton processes such as singlet fission or triplet photoupconversion, inclusion of this mechanism in the narrative of molecular photophysics would allow for strategic synthesis of chromophores for more efficient photon-harvesting applications. Herein, we present a robust formalism which can model these rates using density functional theory (DFT)-based methods within the Franck-Condon and Herzberg-Teller regime. Using an unsubstituted diketo-pyrrolopyrrole (DPP) core as a case study, we illustrate the exciton dynamics along the first four excited states for both singlet and triplet manifolds, showing ultrafast same-spin transfer mechanisms due to all excited states, excluding the first triplet level, being in close energetic proximity (within 0.8 eV of each other). The resulting electron same-spin rates outcompete the electron spin-flipping intersystem crossing (ISC) rates, with excitons firmly obeying Kasha's rule as they cascade down from the high-lying excited states toward the lower states. Furthermore, we calculated that only the first singlet excited state displayed a reasonable probability of triplet exciton generation, of ∼40%, with a near-zero chance of the exciton reverting to the singlet manifold once the electron-hole pair are of parallel spin.

Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy

The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.

Imaging the charge flow in photoexcited molecules would provide key information on photophysical and photochemical processes. Here the authors demonstrate tracking in real time after photoexcitation the change in charge density at a specific site of 2-thiouracil using time-resolved X-ray photoelectron spectroscopy.

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.

Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory

Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S₂ (ππ*) and S₁ (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N²) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.

Ultrafast Imaging of Molecules with Electron Diffraction

Photoexcited molecules convert light into chemical and mechanical energy through changes in electronic and nuclear structure that take place on femtosecond timescales. Gas phase ultrafast electron diffraction (GUED) is an ideal tool to probe the nuclear geometry evolution of the molecules and complements spectroscopic methods that are mostly sensitive to the electronic state. GUED is a passive probing tool that does not alter the molecular properties during the probing process and is sensitive to the spatial distribution of charge in the molecule, including both electrons and nuclei. Improvements in temporal resolution have enabled GUED to capture coherent nuclear motions in molecules in the excited and ground electronic states with femtosecond and subangstrom resolution. Here we present the basic theory of GUED and explain what information is encoded in the diffraction signal, review how GUED has been used to observe coherent structural dynamics in recent experiments, and discuss the advantages and limitations of the method. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Core-Level Spectroscopy of 2-Thiouracil at the Sulfur L1- and L2,3-Edges Utilizing a SASE Free-Electron Laser

In this paper, we report X-ray absorption and core-level electron spectra of the nucleobase derivative 2-thiouracil at the sulfur L₁- and L_2,3-edges. We used soft X-rays from the free-electron laser FLASH2 for the excitation of isolated molecules and dispersed the outgoing electrons with a magnetic bottle spectrometer. We identified photoelectrons from the 2p core orbital, accompanied by an electron correlation satellite, as well as resonant and non-resonant Coster–Kronig and Auger–Meitner emission at the L₁- and L_2,3-edges, respectively. We used the electron yield to construct X-ray absorption spectra at the two edges. The experimental data obtained are put in the context of the literature currently available on sulfur core-level and 2-thiouracil spectroscopy.

The Ultrafast Quantum Dynamics of Photoexcited Adenine-Thymine Basepair Investigated with a Fragment-based Diabatization and a Linear Vibronic Coupling Model

In this contribution we present a quantum dynamical study of the photoexcited hydrogen bonded base pair adenine–thymine (AT) in a Watson–Crick arrangement. To that end, we parametrize Linear Vibronic Coupling (LVC) models with Time-Dependent Density Functional Theory (TD-DFT) calculations, exploiting a fragment diabatization scheme (FrD) we have developed to define diabatic states on the basis of individual chromophores in a multichromophoric system. Wavepacket propagations were run with the multilayer extension of the Multiconfiguration Time-Dependent Hartree method. We considered excitations to the three lowest bright states, a ππ* state of thymine and two ππ* states (L_a and L_b) of adenine, and we found that on the 100 fs time scale the main decay pathways involve intramonomer population transfers toward nπ* states of the same nucleobase. In AT this transfer is less effective than in the isolated nucleobases, because hydrogen bonding destabilizes the nπ* states. The population transfer to the A → T charge transfer state is negligible, making the ultrafast (femtosecond) decay through the proton coupled electron transfer mechanism unlikely, in line with experimental results in apolar solvents. The excitation energy transfer is also very small. We carefully compare the predictions of LVC Hamiltonians obtained with different sets of diabatic states, defined so to match either local states of the two separated monomers or the base pair adiabatic states in the Franck–Condon region. To that end we also extend the flexibility of the FrD-LVC approach, introducing a new strategy to define fragments diabatic states that account for the effect of the rest of the multichromohoric system through a Molecular Mechanics potential.

Controlled ultrafast ππ*-πσ* dynamics in tryptophan-based peptides with tailored micro-environment

Ultrafast charge, energy and structural dynamics in molecules are driven by the topology of the multidimensional potential energy surfaces that determines the coordinated electronic and nuclear motion. These processes are also strongly influenced by the interaction with the molecular environment, making very challenging a general understanding of these dynamics on a microscopic level. Here we use electrospray and mass spectrometry technologies to produce isolated molecular ions with a controlled micro-environment. We measure ultrafast photo-induced ππ*-πσ* dynamics in tryptophan species in the presence of a single, charged adduct. A striking increase of the timescale by more than one order of magnitude is observed when changing the added adduct atom. A model is proposed to rationalize the results, based on the localized and delocalized effects of the adduct on the electronic structure of the molecule. These results offer perspectives to control ultrafast molecular processes by designing the micro-environment on the Angström length scale.

Solvent Effects in the Ultraviolet and X-ray Absorption Spectra of Pyridazine in Aqueous Solution

Electrostatic interaction of the solvent with the solute and fluctuations of the solvent configurations may make excitation energies of the solute different from those in the gas phase. These effects may dominate photoinduced or chemical reaction dynamics in solution systems and can be observed as shifts or broadening of peaks in absorption spectra. In this work, the nitrogen K-edge X-ray absorption spectra were measured for pyridazine in the gas phase and in aqueous solution. The ultraviolet and X-ray absorption spectra of pyridazine in aqueous solution, as well as those in the gas phase, were then calculated with models based on the algebraic-diagrammatic construction through second order [ADC(2)] with the resolution-of-identity (RI) approximation and compared with the spectra obtained in experiments. For aqueous solution, explicit local solvation structures were extracted from an ab initio molecular dynamics (AIMD) trajectory of pyridazine in bulk water, and RI-ADC(2) was combined with the conductor-like screening model (COSMO). The experimental absorption spectra of pyridazine in aqueous solution were reproduced with good accuracy by theoretical treatment of an ensemble containing the explicit local solvation structures of pyridazine with relevant water molecules combined with the COSMO solvation model of water for long-range solvation.

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.

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation

Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized ¹B_2u (ππ*) (S₂) and ¹B_3u (nπ*) (S₁) states, the participation of the optically dark ¹A_u (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite ¹A_u (nπ*) and ¹B_3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The ¹A_u (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

Modeling radiative and non-radiative pathways at both the Franck-Condon and Herzberg-Teller approximation level

Here, we present a concise model that can predict the photoluminescent properties of a given compound from first principles, both within and beyond the Franck-Condon approximation. The formalism required to compute fluorescence, Internal Conversion (IC), and Inter-System Crossing (ISC) is discussed. The IC mechanism, in particular, is a difficult pathway to compute due to difficulties associated with the computation of required bosonic configurations and non-adiabatic coupling elements. Here, we offer a discussion and breakdown on how to model these pathways at the Density Functional Theory (DFT) level with respect to its computational implementation, strengths, and current limitations. The model is then used to compute the photoluminescent quantum yield (PLQY) of a number of small but important compounds: anthracene, tetracene, pentacene, diketo-pyrrolo-pyrrole (DPP), and Perylene Diimide (PDI) within a polarizable continuum model. Rate constants for fluorescence, IC, and ISC compare well for the most part with respect to experiment, despite triplet energies being overestimated to a degree. The resulting PLQYs are promising with respect to the level of theory being DFT. While we obtained a positive result for PDI within the Franck-Condon limit, the other systems require a second order correction. Recomputing quantum yields with Herzberg-Teller terms yields PLQYs of 0.19, 0.08, 0.04, 0.70, and 0.99 for anthracene, tetracene, pentacene, DPP, and PDI, respectively. Based on these results, we are confident that the presented methodology is sound with respect to the level of quantum chemistry and presents an important stepping stone in the search for a tool to predict the properties of larger coupled systems.

Near-Infrared Light-Regulated Drug-Food Homologous Bioactive Molecules and Photothermal Collaborative Precise Antibacterial Therapy Nanoplatform with Controlled Release Property

Herein, a collaborative precise antibacterial wound healing therapy nanoplatform integrating drug-food homologous bioactive molecule (cinnamaldehyde, CA) with photothermal therapy (PTT) is presented. Copper-gallic acids-cinnamaldehyde-polydopamine nanorods (Cu-GA-CA-PDA NRs) with near-infrared light (NIR)-controlled CA release property are fabricated, which also integrate CA and photothermal synergistic sterilization, as well as antioxidant, anti-inflammatory, and anti-infection capacities. The characteristics of NIR-mediated CA release and photothermal response of Cu-GA-CA-PDA NRs support their excellent sterilization performance in vitro/in vivo. In addition, under the guidance of NIR, Cu-GA-CA-PDA NRs can hinder the formation of inflammatory cells, reduce oxidative stress damage, accelerate the regeneration of skin tissues in S. aureus-infected wound sites, and achieve the goal of promoting wound healing. Therefore, NIR-mediated Cu-GA-CA-PDA NRs with multifunctional biological activities provide a highly competitive strategy for curing bacteria-infected wounds.

Transient resonant Auger-Meitner spectra of photoexcited thymine

We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a ππ* excited state, thymine is known to undergo internal conversion to an nπ* excited state with a strong resonance at the oxygen K-edge, red-shifted from the ground state π* resonances of thymine (see our previous study Wolf, et al., Nat. Commun., 2017, 8, 29). We resolve and compare the Auger-Meitner electron spectra associated both with the excited state and ground state resonances, and distinguish participator and spectator decay contributions. Furthermore, we observe simultaneously with the decay of the nπ* state signatures the appearance of additional resonant Auger-Meitner contributions at photon energies between the nπ* state and the ground state resonances. We assign these contributions to population transfer from the nπ* state to a ππ* triplet state via intersystem crossing on the picosecond timescale based on simulations of the X-ray absorption spectra in the vibrationally hot triplet state. Moreover, we identify signatures from the initially excited ππ* singlet state which we have not observed in our previous study.

Impact of the Dynamic Electron Correlation on the Unusually Long Excited-State Lifetime of Thymine

Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S₁-trapping, S₂-trapping, and S₁&S₂-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S₂-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S₁-trapping mechanism with two lifetimes, τ₁ = 30 ± 1 fs and τ₂ = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S₂ (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S₂→S₁ internal conversion and resides around the minimum on the S₁ (n_Oπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S₂→S₁ internal conversion is mediated by fast bond length alternation distortion, the subsequent S₁→S₀ occurs through several conical intersections, involving a slow puckering motion.

Ultrafast strong-field dissociation of vinyl bromide: An attosecond transient absorption spectroscopy and non-adiabatic molecular dynamics study

Attosecond extreme ultraviolet (XUV) and soft x-ray sources provide powerful new tools for studying ultrafast molecular dynamics with atomic, state, and charge specificity. In this report, we employ attosecond transient absorption spectroscopy (ATAS) to follow strong-field-initiated dynamics in vinyl bromide. Probing the Br M edge allows one to assess the competing processes in neutral and ionized molecular species. Using ab initio non-adiabatic molecular dynamics, we simulate the neutral and cationic dynamics resulting from the interaction of the molecule with the strong field. Based on the dynamics results, the corresponding time-dependent XUV transient absorption spectra are calculated by applying high-level multi-reference methods. The state-resolved analysis obtained through the simulated dynamics and related spectral contributions enables a detailed and quantitative comparison with the experimental data. The main outcome of the interaction with the strong field is unambiguously the population of the first three cationic states, D ₁, D ₂, and D ₃. The first two show exclusively vibrational dynamics while the D ₃ state is characterized by an ultrafast dissociation of the molecule via C-Br bond rupture within 100 fs in 50% of the analyzed trajectories. The combination of the three simulated ionic transient absorption spectra is in excellent agreement with the experimental results. This work establishes ATAS in combination with high-level multi-reference simulations as a spectroscopic technique capable of resolving coupled non-adiabatic electronic-nuclear dynamics in photoexcited molecules with sub-femtosecond resolution.

Using X-ray free-electron lasers for spectroscopy of molecular catalysts and metalloenzymes

The metal centres in metalloenzymes and molecular catalysts are responsible for the rearrangement of atoms and electrons during complex chemical reactions, and they enable selective pathways of charge and spin transfer, bond breaking/making and the formation of new molecules. Mapping the electronic structural changes at the metal sites during the reactions gives a unique mechanistic insight that has been difficult to obtain to date. The development of X-ray free-electron lasers (XFELs) enables powerful new probes of electronic structure dynamics to advance our understanding of metalloenzymes. The ultrashort, intense and tunable XFEL pulses enable X-ray spectroscopic studies of metalloenzymes, molecular catalysts and chemical reactions, under functional conditions and in real time. In this Technical Review, we describe the current state of the art of X-ray spectroscopy studies at XFELs and highlight some new techniques currently under development. With more XFEL facilities starting operation and more in the planning or construction phase, new capabilities are expected, including high repetition rate, better XFEL pulse control and advanced instrumentation. For the first time, it will be possible to make real-time molecular movies of metalloenzymes and catalysts in solution, while chemical reactions are taking place.

Nonadiabatic Absorption Spectra and Ultrafast Dynamics of DNA and RNA Photoexcited Nucleobases

We have recently proposed a protocol for Quantum Dynamics (QD) calculations, which is based on a parameterisation of Linear Vibronic Coupling (LVC) Hamiltonians with Time Dependent (TD) Density Functional Theory (TD-DFT), and exploits the latest developments in multiconfigurational TD-Hartree methods for an effective wave packet propagation. In this contribution we explore the potentialities of this approach to compute nonadiabatic vibronic spectra and ultrafast dynamics, by applying it to the five nucleobases present in DNA and RNA. For all of them we computed the absorption spectra and the dynamics of ultrafast internal conversion (100 fs timescale), fully coupling the first 2-3 bright states and all the close by dark states, for a total of 6-9 states, and including all the normal coordinates. We adopted two different functionals, CAM-B3LYP and PBE0, and tested the effect of the basis set. Computed spectra are in good agreement with the available experimental data, remarkably improving over pure electronic computations, but also with respect to vibronic spectra obtained neglecting inter-state couplings. Our QD simulations indicate an effective population transfer from the lowest energy bright excited states to the close-lying dark excited states for uracil, thymine and adenine. Dynamics from higher-energy states show an ultrafast depopulation toward the more stable ones. The proposed protocol is sufficiently general and automatic to promise to become useful for widespread applications.