Quantum dynamics of the ππ*/nπ* decay of the epigenetic nucleobase 1,5-dimethyl-cytosine in the gas phase.

A partial ultrafast ππ* → nπ* transfer is predicted. Many vibrational modes are activated, but oscillations of bonds and angles are quickly damped

Quantum Dynamics of Decay Processes in Photoexcited Nucleobases and Small Oligomers

DNA strongly absorbs UV light and this can trigger harmful processes causing alteration of the genetic code and cellular death. However, DNA and its components, the nucleobases, have an intrinsic capability to dissipate the excess electronic energy into heat thanks to a combination of intra-molecular and inter-molecular non-radiative processes, which strongly reduce the risk of potentially dangerous photo reactions. Although deeply studied in literature, many aspects of these mechanisms remain to be clarified. DNA can also be seen as a prototype of many multichromophoric systems so the interest in its photophysics goes beyond its biological relevance. This thesis is devoted to a computational study of DNA photophysics with dynamical approaches. The ultra-fast decays in DNA occur through non-radiative transitions triggered by the coupling between electrons and nuclei motions. These are intrinsically quantum phenomena and therefore we chose to tackle them with fully quantum dynamical approaches. More specifically, in this thesis, we have proposed and applied a computational protocol based on the parameterization with time-dependent density functional theory of Linear Vibronic Coupling models able to describe the competition between intra-base and inter-base decay mechanisms, combined with non-radiative propagations of vibronic wavepackets with advanced multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) approaches. With this strategy we have tackled the study of the ultra-fast light-activated dynamics (∼ 100 fs) of DNA components of growing complexity, starting from isolated nucleobases and their derivatives up to double-stranded tetrads

A Fragment Diabatization Linear Vibronic Coupling Model for Quantum Dynamics of Multichromophoric Systems: Population of the Charge Transfer State in the Photoexcited Guanine Cytosine Pair

We introduce a method (FrD-LVC) based on a fragment diabatization (FrD) for the parametrization of a Linear Vibronic Coupling (LVC) model suitable for studying the photophysics of multichromophore systems. In combination with effective quantum dynamics (QD) propagations with multilayer multiconfigurational time-dependent Hartree (ML-MCTDH), the FrD-LVC approach gives access to the study of the competition between intra-chromophore decays, like those at conical intersections, and inter-chromophore processes, like exciton localization/delocalization and the involvement of charge transfer (CT) states. We used FrD-LVC parametrized with TD-DFT calculations, adopting either CAM-B3LYP or ωB97X-D functionals, to study the ultrafast photoexcited QD of a Guanine-Cytosine (GC) hydrogen bonded pair, within a Watson-Crick arrangement, considering up to 12 coupled diabatic electronic states and the effect of all the 99 vibrational coordinates. The bright excited states localized on C and, especially, on G are predicted to be strongly coupled to the G->C CT state which is efficiently and quickly populated after an excitation to any of the four lowest energy bright local excited states. Our QD simulations show that more than 80% of the excited population on G and ~50% of that on C decays to this CT state in less than 50 fs. We investigate the role of vibronic effects in the population of the CT state and show it depends mainly on its large reorganization energy so that it can occur even when it is significantly less stable than the bright states in the Franck-Condon region. At the same time, we document that the formation of the GC pair almost suppresses the involvement of dark nπ* excited states in the photoactivated dynamics.</div

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

Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine

[Image: see text] We present a viable protocol to compute vibrational resonance Raman (vRR) spectra for systems with several close-lying and potentially coupled electronic states. It is based on the parametrization of linear vibronic coupling (LVC) models from time-dependent density functional theory (TD-DFT) calculations and quantum dynamics propagations of vibronic wavepackets with the multilayer version of the multiconfiguration time-dependent Hartree (ML-MCTDH) method. Our approach is applied to thymine considering seven coupled electronic states, comprising the three lowest bright states, and all vibrational coordinates. Computed vRR at different excitation wavelengths are in good agreement with the available experimental data. Up to 250 nm the signal is dominated by the lowest HOMO → LUMO transition, whereas at 233 nm, in the valley between the two lowest energy absorption bands, the contributions of all the three bright states, and their interferences and couplings, are important. Inclusion of solvent (water) effects improves the agreement with experiment, reproducing the coalescence of vibrational bands due to CC and C=O stretchings. With our approach we disentangle and assess the effect of interferences between the contribution of different quasi-resonant states to the transition polarizability and the effect of interstate couplings. Our findings strongly suggest that in cases of close-lying and potentially coupled states a simple inclusion of interference effects is not sufficient, and a fully nonadiabatic computation should instead be performed. We also document that for systems with strong couplings and quasi-degenerate states, the use of HT perturbative approach, not designed for these cases, may lead to large artifacts