We determined the complete relaxation dynamics of pyrene in ethanol from the second bright state, employing experimental and theoretical broadband heterodyne detected transient grating and two-dimensional photon echo (2DPE) spectroscopy, using pulses with duration of 6 fs and covering a spectral range spanning from 250 to 300 nm. Multiple lifetimes are assigned to conical intersections through a cascade of electronic states, eventually leading to a rapid population of the lowest long-living excited state and subsequent slow vibrational cooling. The lineshapes in the 2DPE spectra indicate that the efficiency of the population transfer depends on the kinetic energy deposited into modes required to reach a sloped conical intersection, which mediates the decay to the lowest electronic state. The presented experimental-theoretical protocol paves the way for studies on deep-ultraviolet-absorbing biochromophores ubiquitous in genomic and proteic systems.

We outline a computational approach for nonlinear electronic spectra, which accounts for the electronic energy fluctuations due to nuclear degrees of freedom and explicitly incorporates the fluctuations of higher excited states, induced by the dynamics in the photoactive state(s). This approach is based on mixed quantum-classical dynamics simulations. Tedious averaging over multiple trajectories is avoided by employing the linearly displaced Brownian harmonic oscillator to model the correlation functions. The present strategy couples accurate computations of the high-lying excited state manifold with dynamics simulations. The application is made to the two-dimensional electronic spectra of pyrene, a polycyclic aromatic hydrocarbon characterized by an ultrafast (few tens of femtoseconds) decay from the bright S2 state to the dark S1 state. The spectra for waiting times t2 = 0 and t2 = 1 ps demonstrate the ability of this approach to model electronic state fluctuations and realistic lineshapes. Comparison with experimental spectra [Krebs et al., New Journal of Physics, 2013, 15, 085016] shows excellent agreement and allows us to unambiguously assign the excited state absorption features.

Pump-probe electronic spectroscopy using femtosecond laser pulses has evolved into a standard tool for tracking ultrafast excited state dynamics. Its two-dimensional (2D) counterpart is becoming an increasingly available and promising technique for resolving many of the limitations of pump-probe caused by spectral congestion. The ability to simulate pump-probe and 2D spectra from ab initio computations would allow one to link mechanistic observables like molecular motions and the making/breaking of chemical bonds to experimental observables like excited state lifetimes and quantum yields. From a theoretical standpoint, the characterization of the electronic transitions in the visible (Vis)/ultraviolet (UV), which are excited via the interaction of a molecular system with the incoming pump/probe pulses, translates into the determination of a computationally challenging number of excited states (going over 100) even for small/medium sized systems. A protocol is therefore required to evaluate the fluctuations of spectral properties like transition energies and dipole moments as a function of the computational parameters and to estimate the effect of these fluctuations on the transient spectral appearance. In the present contribution such a protocol is presented within the framework of complete and restricted active space self-consistent field theory and its second-order perturbation theory extensions. The electronic excited states of adenine have been carefully characterized through a previously presented computational recipe [Nenov et al., Comput. Theor. Chem. 1040-1041, 295-303 (2014)]. A wise reduction of the level of theory has then been performed in order to obtain a computationally less demanding approach that is still able to reproduce the characteristic features of the reference data. Foreseeing the potentiality of 2D electronic spectroscopy to track polynucleotide ground and excited state dynamics, and in particular its expected ability to provide conformational dependent fingerprints in dimeric systems, the performances of the selected reduced level of calculations have been tested in the construction of 2D electronic spectra for the in vacuo adenine monomer and the unstacked adenine homodimer, thereby exciting the Lb/La transitions with the pump pulse pair and probing in the Vis to near ultraviolet spectral window.

A computational strategy to simulate two-dimensional electronic spectra (2DES) is introduced, which allows us to analyse ground state dynamics and to sample and measure different conformations attained by flexible molecular systems in solution. An explicit mixed quantum mechanics/molecular mechanics (QM/MM) approach is employed for the evaluation of the necessary electronic excited state energies and transition dipole moments. The method is applied towards a study of the highly flexible water-solvated adenine-adenine monophosphate (ApA), a system featuring two interacting adenine moieties that display various intermolecular arrangements, known to deeply affect their photochemical outcome. Molecular dynamics simulations and cluster analysis have been used to select the molecular conformations, reducing the complexity of the flexible ApA conformational space. By using our sum-over-states (SOS) approach to obtain the 2DES spectra for each of these selected conformations, we can discern spectral changes and relate them to specific nuclear arrangements: close lying π-stacked bases exhibit a splitting of their respective 1La signal traces; T-stacked bases exhibit the appearance of charge transfer states in the low-energy Vis probing window while displaying no 1La splitting, being particularly favoured when promoting amino to 5-ring interactions; unstacked and distant adenine moieties exhibit signals similar to those of the adenine monomer, as is expected for non-interacting nucleobases. 2DES maps reveal the spectral fingerprints associated with specific molecular conformations, and are thus a promising option to enable their quantitative spectroscopic detection beyond standard 1D pump-probe techniques. This is expected to aid the understanding of how nucleobase aggregation controls and modulates the photostability and photo-damage of extended DNA/RNA systems.

The decay channels of singlet excited adenine uracil monophosphate (ApU) in water are studied with CASPT2//CASSCF:MM potential energy calculations and simulation of the 2D-UV spectroscopic fingerprints with the aim of elucidating the role of the different electronic states of the stacked conformer in the excited state dynamics. The adenine (1) La state can decay without a barrier to a conical intersection with the ground state. In contrast, the adenine (1) Lb and uracil S(U) states have minima that are separated from the intersections by sizeable barriers. Depending on the backbone conformation, the CT state can undergo inter-base hydrogen transfer and decay to the ground state through a conical intersection, or it can yield a long-lived minimum stabilized by a hydrogen bond between the two ribose rings. This suggests that the (1) Lb , S(U) and CT states of the stacked conformer may all contribute to the experimental lifetimes of 18 and 240 ps. We have also simulated the time evolution of the 2D-UV spectra and provide the specific fingerprint of each species in a recommended probe window between 25 000 and 38 000 cm(-1) in which decongested, clearly distinguishable spectra can be obtained. This is expected to allow the mechanistic scenarios to be discerned in the near future with the help of the corresponding experiments. Our results reveal the complexity of the photophysics of the relatively small ApU system, and the potential of 2D-UV spectroscopy to disentangle the photophysics of multichromophoric systems.

The SOS//QM/MM [Rivalta et al., Int. J. Quant. Chem., 2014, 114, 85] method consists of an arsenal of computational tools allowing accurate simulation of one-dimensional (1D) and bi-dimensional (2D) electronic spectra of monomeric and dimeric systems with unprecedented details and accuracy. Prominent features like doubly excited local and excimer states, accessible in multi-photon processes, as well as charge-transfer states arise naturally through the fully quantum-mechanical description of the aggregates. In this contribution the SOS//QM/MM approach is extended to simulate time-resolved 2D spectra that can be used to characterize ultrafast excited state relaxation dynamics with atomistic details. We demonstrate how critical structures on the excited state potential energy surface, obtained through state-of-the-art quantum chemical computations, can be used as snapshots of the excited state relaxation dynamics to generate spectral fingerprints for different de-excitation channels. The approach is based on high-level multi-configurational wavefunction methods combined with non-linear response theory and incorporates the effects of the solvent/environment through hybrid quantum mechanics/molecular mechanics (QM/MM) techniques. Specifically, the protocol makes use of the second-order Perturbation Theory (CASPT2) on top of Complete Active Space Self Consistent Field (CASSCF) strategy to compute the high-lying excited states that can be accessed in different 2D experimental setups. As an example, the photophysics of the stacked adenine-adenine dimer in a double-stranded DNA is modeled through 2D near-ultraviolet (NUV) spectroscopy.