Time-resolved x-ray spectroscopy of nucleobases and their thionated analogs

Direct Observation of the ππ* to nπ* Transition in 2-Thiouracil via Time-Resolved NEXAFS Spectroscopy

The photophysics of nucleobases has been the subject of both theoretical and experimental studies over the past decades due to the challenges posed by resolving the steps of their radiationless relaxation dynamics, which cannot be described in the framework of the Born-Oppenheimer approximation (BOA). In this context, the ultrafast dynamics of 2-thiouracil has been investigated with a time-resolved NEXAFS study at the Free Electron Laser FLASH. Near Edge X-ray Absorption Fine Structure spectroscopy (NEXAFS) can be used to observe electronic transitions in ultrafast molecular relaxation. We performed time-resolved UV-pump/X-ray probe absorption measurements at the sulfur 2s (L1) and 2p (L2/3) edges. We are able to identify absorption features corresponding to the S2 (ππ*) and S1 (nπ*) electronic states. We observe a delay of 102 ± 11 fs in the population of the nπ* state with respect to the initial optical excitation and interpret the delay as the time scale for the S2 → S1 internal conversion. We furthermore identify oscillations in the absorption signal that match a similar observation in our previous X-ray photoelectron spectroscopy study on the same molecule.

Method to Compute the Interaction Energy of a Molecule in Ground and Excited States with a Discrete Environment: The Case of Uracil in Water

In this work, we present a method that is able to compute the interaction energy of a system of interest, in the ground or excited state, with an arbitrary number of water molecules representing the environment. As a test case, we take uracil. We considered five clusters containing 1, 12, 24, 26, and 37 water molecules. The method is a first step toward a more general approach to determining the interaction energy between a molecule, treated at a high level of theory, and a complex molecular environment that can be described as an explicit solvent model. Ground and excited electronic states of uracil were optimized in free space at the variational quantum Monte Carlo (QMC) level. In this way, we sampled electronic configurations that are used to compute all the contributions to the interaction energy with the environment. Excitation energies from the ground state were computed at the diffusion Monte Carlo (DMC) level. Numerical results are in agreement with available literature data on the solvatochromic effect on the n → π* and π → π* vertical transitions of uracil in water. Our method provides specific contributions arising from Pauli repulsion, electrostatic, polarization, and dispersion interactions.

Photoinduced hydrogen dissociation in thymine predicted by coupled cluster theory

The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate. Today, it is widely believed that its ultrafast excited state gas phase decay stems from a radiationless transition from the bright ππ* state to a dark nπ* state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the early gas phase ultrafast dynamics in thymine at the highest level of theory to date. This is made possible by performing wavepacket dynamics with a recently developed coupled cluster method. Our simulation confirms an ultrafast ππ* to nπ* transition (τ = 41 ± 14 fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with experiment. We also predict an as-yet uncharacterized πσ* channel that leads to hydrogen dissociation at one of the two N-H bonds. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel.

Coupled Cluster Theory for Nonadiabatic Dynamics: Nuclear Gradients and Nonadiabatic Couplings in Similarity Constrained Coupled Cluster Theory

Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply the method in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8(19), 4801-4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method.