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

Internal Conversion Cascade in a Carbon Nanobelt: A Multiconfigurational Quantum Dynamical Study

Carbon nanobelts feature intriguing photophysical properties, due to their high symmetry and structural rigidity. Here, we consider a (6,6) armchair carbon nanobelt, i.e., the very first carbon nanobelt to be synthesized [Povie et al., Science 2017, 356, 172] and characterize the internal conversion dynamics using multiconfigurational quantum dynamics via the multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) method. A symmetry-adapted linear vibronic coupling Hamiltonian for 26 electronic states and 210 vibrational modes is employed. Electronic excitations are found to decay through a dense manifold of excited states, which interact via multiple conical intersections, while inducing minimal geometry change. It is shown that a rapid coherent decay, exhibiting a nonvanishing quantum flux on a time scale of less than 50 fs, transitions toward a slower, decoherent decay at longer times. As previously suggested in the literature, electronic relaxation is hindered by phonon bottlenecks such that a stepwise internal conversion cascade is observed. The computed vibronic absorption spectrum is shown to be in good agreement with the experimental spectrum.

The Excited State Dynamics of a Mutagenic Guanosine Etheno Adduct Investigated by Femtosecond Fluorescence Spectroscopy and Quantum Mechanical Calculations

Femtosecond fluorescence upconversion experiments were combined with CASPT2 and time dependent DFT calculations to characterize the excited state dynamics of the mutagenic etheno adduct 1,N2-etheno-2'-deoxyguanosine (ϵdG). This endogenously formed lesion is attracting great interest because of its ubiquity in human tissues and its highly mutagenic properties. The ϵdG fluorescence is strongly modified with respect to that of the canonical nucleoside dG, notably by an about 6-fold increase in fluorescence lifetime and quantum yield at neutral pH. In addition, femtosecond fluorescence upconversion experiments reveal the presence of two emission bands with maxima at 335 nm for the shorter-lived and 425 nm for the longer-lived. Quantum mechanical calculations rationalize these findings and provide absorption and fluorescence spectral shapes similar to the experimental ones. Two different bright minima are located on the potential energy surface of the lowest energy singlet excited state. One planar minimum, slightly more stable, is associated with the emission at 335 nm, whereas the other one, with a bent etheno ring, is associated with the red-shifted emission.

Vibrational Electronic-Thermofield Coupled Cluster (VE-TFCC) Theory for Quantum Simulations of Vibronic Coupling Systems at Thermal Equilibrium

A vibrational electronic-thermofield coupled cluster (VE-TFCC) approach is developed to calculate thermal properties of non-adiabatic vibronic coupling systems. The thermofield (TF) theory and a mixed linear exponential ansatz based on second-quantized Bosonic construction operators are introduced to propagate the thermal density operator as a "pure state" in the Bogoliubov representation. Through this compact representation of the thermal density operator, the approach is basis-set-free and scales classically (polynomial) as the number of degrees of freedoms (DoF) in the system increases. The VE-TFCC approach is benchmarked with small test models and a real molecular compound (CoF4- anion) against the conventional sum over states (SOS) method and applied to calculate thermochemistry properties of a gas-phase reaction: CoF3 + F- → CoF4-. Results shows that the VE-TFCC approach, in conjunction with vibronic models, provides an effective protocol for calculating thermodynamic properties of vibronic coupling systems.

The photoactivated dynamics of dGpdC and dCpdG sequences in DNA: a comprehensive quantum mechanical study

Study of alternating DNA GC sequences by different time-resolved spectroscopies has provided fundamental information on the interaction between UV light and DNA, a process of great biological importance. Multiple decay paths have been identified, but their interplay is still poorly understood. Here, we characterize the photophysics of GC-DNA by integrating different computational approaches, to study molecular models including up to 6 bases described at a full quantum mechanical level. Quantum dynamical simulations, exploiting a nonadiabatic linear vibronic coupling (LVC) model, coupled with molecular dynamics sampling of the initial structures of a (GC)5 DNA duplex, provide new insights into the photophysics in the sub-picosecond time-regime. They indicate a substantial population transfer, within 50 fs, from the spectroscopic states towards G → C charge transfer states involving two stacked bases (CTintra), thus explaining the ultrafast disappearance of fluorescence. This picture is consistent with that provided by quantum mechanical geometry optimizations, using time dependent-density functional theory and a polarizable continuum model, which we use to parametrize the LVC model and to map the main excited state deactivation pathways. For the first time, the infrared and excited state absorption signatures of the various states along these pathways are comprehensively mapped. The computational models suggest that the main deactivation pathways, which, according to experiment, lead to ground state recovery on the 10-50 ps time scale, involve CTintra followed by interstrand proton transfer from the neutral G to C-. Our calculations indicate that CTintra is populated to a larger extent and more rapidly in GC than in CG steps and suggest the likely involvement of monomer-like and interstrand charge transfer decay routes for isolated and less stacked CG steps. These findings underscore the importance of the DNA sequence and thermal fluctuations for the dynamics. They will also aid the interpretation of experimental results on other sequences.

Linear Vibronic Coupling Approach for Surface-Enhanced Raman Scattering: Quantifying the Charge-Transfer Enhancement Mechanism

The outstanding amplification observed in surface-enhanced Raman scattering (SERS) is due to several enhancement mechanisms, and standing out among them are the plasmonic (PL) and charge-transfer (CT) mechanisms. The theoretical estimation of the enhancement factors of the CT mechanism is challenging because the excited-state coupling between bright plasmons and dark CT states must be properly introduced into the model to obtain reliable intensities. In this work, we aim at simulating electrochemical SERS spectra, considering models of pyridine on silver clusters subjected to an external electric field E⃗ that represents the effect of an electrode potential Vel. The method adopts quantum dynamical propagations of nuclear wavepackets on the coupled PL and CT states described with linear vibronic coupling models parametrized for each E⃗ through a fragment-based maximum-overlap diabatization. By presenting results at different values of E⃗, we show that indeed there is a relation between the population transferred to the CT states and the total scattered intensity. The tuning and detuning processes of the CT states with the bright PLs as a function of the electric field are in good agreement with those observed in experiments. Finally, our estimations for the CT enhancement factors predict values in the order of 105 to 106, meaning that when the CT and PL states are both in resonance with the excitation wavelength, the CT and PL enhancements are comparable, and vibrational bands whose intensity is amplified by different mechanisms can be observed together, in agreement with what was measured by typical experiments on silver electrodes.

Time dependent vibrational electronic coupled cluster (VECC) theory for non-adiabatic nuclear dynamics

A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born-Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV-VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)].

Quantum-Classical Protocol for Efficient Characterization of Absorption Lineshape and Fluorescence Quenching upon Aggregation: The Case of Zinc Phthalocyanine Dyes

A quantum-classical protocol that incorporates Jahn-Teller vibronic coupling effects and cluster analysis of molecular dynamics simulations is reported, providing a tool for simulations of absorption spectra and ultrafast nonadiabatic dynamics in large molecular photosystems undergoing aggregation in solution. Employing zinc phthalocyanine dyes as target systems, we demonstrated that the proposed protocol provided fundamental information on vibronic, electronic couplings and thermal dynamical effects that mostly contribute to the absorption spectra lineshape and the fluorescence quenching processes upon dye aggregation. Decomposing the various effects arising upon dimer formation, the structure-property relations associated with their optical responses have been deciphered at atomistic resolution.

Mind the GAP: quantifying the breakdown of the linear vibronic coupling Hamiltonian

Excited state dynamics play a critical role across a broad range of scientific fields. Importantly, the highly non-equilibrium nature of the states generated by photoexcitation means that excited state simulations should usually include an accurate description of the coupled electronic-nuclear motion, which often requires solving the time-dependent Schrödinger equation (TDSE). One of the biggest challenges for these simulations is the requirement to calculate the PES over which the nuclei evolve. An effective approach for addressing this challenge is to use the approximate linear vibronic coupling (LVC) Hamiltonian, which enables a model potential to be parameterised using relatively few quantum chemistry calculations. However, this approach is only valid provided there are no large amplitude motions in the excited state dynamics. In this paper we introduce and deploy a metric, the global anharmonicity parameter (GAP), which can be used to assess the accuracy of an LVC potential. Following its derivation, we illustrate its utility by applying it to three molecules exhibiting different rigidity in their excited states.

The Resonance Raman Spectrum of Cytosine in Water: Analysis of the Effect of Specific Solute-Solvent Interactions and Non-Adiabatic Couplings

In this contribution, we report a computational study of the vibrational Resonance Raman (vRR) spectra of cytosine in water, on the grounds of potential energy surfaces (PES) computed by time-dependent density functional theory (TD-DFT) and CAM-B3LYP and PBE0 functionals. Cytosine is interesting because it is characterized by several close-lying and coupled electronic states, challenging the approach commonly used to compute the vRR for systems where the excitation frequency is in quasi-resonance with a single state. We adopt two recently developed time-dependent approaches, based either on quantum dynamical numerical propagations of vibronic wavepackets on coupled PES or on analytical correlation functions for cases in which inter-state couplings were neglected. In this way, we compute the vRR spectra, considering the quasi-resonance with the eight lowest-energy excited states, disentangling the role of their inter-state couplings from the mere interference of their different contributions to the transition polarizability. We show that these effects are only moderate in the excitation energy range explored by experiments, where the spectral patterns can be rationalized from the simple analysis of displacements of the equilibrium positions along the different states. Conversely, at higher energies, interference and inter-state couplings play a major role, and the adoption of a fully non-adiabatic approach is strongly recommended. We also investigate the effect of specific solute-solvent interactions on the vRR spectra, by considering a cluster of cytosine, hydrogen-bonded by six water molecules, and embedded in a polarizable continuum. We show that their inclusion remarkably improves the agreement with the experiments, mainly altering the composition of the normal modes, in terms of internal valence coordinates. We also document cases, mostly for low-frequency modes, in which a cluster model is not sufficient, and more elaborate mixed quantum classical approaches, in explicit solvent models, need to be applied.

Computing the electronic circular dichroism spectrum of DNA quadruple helices of different topology: A critical test for a generalized excitonic model based on a fragment diabatization

In this study, we exploit a recently developed fragment diabatization-based excitonic model, FrDEx, to simulate the electronic circular dichroism (ECD) spectra of three guanine-rich DNA sequences arranged in guanine quadruple helices with different topologies: thrombin binding aptamer (antiparallel), c-Myc promoter (parallel), and human telomeric sequence (3+1 hybrid). Starting from time-dependent density functional theory (TD-DFT) calculations with the M052X functional, we apply our protocol to parameterize the FrDEX Hamiltonian, which accounts for electron density overlap and includes both the coupling with charge transfer transitions and the effect of the surrounding bases on the local excitation of each chromophore. The TD-DFT/M052X spectral shapes are in good agreement with the experimental ones, the main source of discrepancy being related to the intrinsic error on the computed transition energies of guanine monomer. FrDEx spectra are fairly close to the reference TD-DFT ones, allowing a significant advance with respect to a more standard excitonic Hamiltonian. We also show that the ECD spectra are sensitive to the inclusion of the inner K + $$ {}^{+} $$ cation in the calculation.

Effect of the Thermal Fluctuations of the Photophysics of GC and CG DNA Steps: A Computational Dynamical Study

Here we refine and assess two computational procedures aimed to include the effect of thermal fluctuations on the electronic spectra and the ultrafast excited state dynamics of multichromophore systems, focusing on DNA duplexes. Our approach is based on a fragment diabatization procedure that, from a given Quantum Mechanical (QM) reference method, can provide the parameters (energy and coupling) of the reference diabatic states on the basis of the isolated fragments, either for a purely electronic excitonic Hamiltonian (FrDEx) or a linear vibronic coupling Hamiltonian (FrD-LVC). After having defined the most cost-effective procedure for DNA duplexes on two smaller fragments, FrDEx is used to simulate the absorption and Electronic Circular Dichroism (ECD) spectra of (GC)₅ sequences, including the coupling with the Charge Transfer (CT) states, on a number of structures extracted from classical Molecular Dynamics (MD) simulations. The computed spectra are close to the reference TD-DFT calculations and fully consistent with the experimental ones. We then couple MD simulations and FrD-LVC to simulate the interplay between local excitations and CT transitions, both intrastrand and interstrand, in GC and CG steps when included in a oligoGC or in oligoAT DNA sequence. We predict that for both sequences a substantial part of the photoexcited population on G and C decays, within 50-100 fs, to the corresponding intrastrand CT states. This transfer is more effective for GC steps that, on average, are more closely stacked than CG ones.

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

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.

Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform

We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.

Computational Model for Electrochemical Surface-Enhanced Raman Scattering: Key Role of the Surface Charges and Synergy between Electromagnetic and Charge-Transfer Enhancement Mechanisms

We present a computational model for electrochemical surface-enhanced Raman scattering (EC-SERS). The surface excess of charge induced by the electrode potential (V_el) was introduced by applying an external electric field to a set of clusters [Ag_n]^q with (n, q) of (19, ±1) or (20, 0) on which a molecule adsorbs. Using DFT/TD-DFT calculations, these metal-molecule complexes were classified by the adsorbate partial charge, and the main V_el-dependent properties were simultaneously studied with the aid of vibronic resonance Raman computations, namely, changes on the vibrational wavenumbers, relative intensities, and enhancement factors (EFs) for all SERS mechanisms: chemical or nonresonant, and resonance Raman with bright states of the adsorbate, charge-transfer (CT) states, and plasmon-like excitations on the metal cluster. We selected two molecules to test our model, pyridine, for which V_el has a remarkable effect, and 9,10-bis((E)-2-(pyridin-4-yl)vinyl)anthracene, which is almost insensitive to the applied bias. The results nicely reproduced most of the experimental observations, while the limitations of our approach were critically evaluated. We detected that accounting explicitly for the surface charges is key for EC-SERS models and that the highest calculated EFs, up to 10⁷ to 10⁸, are obtained by interstate coupling of bright local excitations of the metal cluster and CT states. These results highlight the importance of nonadiabatic effects in SERS and the capabilities of EC-SERS as a technique with potential to study excited-state coupling by tuning the CT and plasmon-like states by manipulating V_el.

Nucleic Acids as a Playground for the Computational Study of the Photophysics and Photochemistry of Multichromophore Assemblies

ConspectusThe interaction between light and multichromophoric assemblies (MCAs) is the primary event of many fundamental processes, from photosynthesis to organic photovoltaics, and it triggers dynamical processes that share remarkable similarities at the molecular scale: light absorption, energy and charge transfer, internal conversions, emission, and so on. Those events often involve many chromophores and different excited electronic states that are coupled on an ultrafast time scale. This Account aims to discuss some of the chemical physical effects ruling these processes, a fundamental step toward their control, based on our experience on nucleic acids.In the last 15 years, we have, indeed, studied the photophysics and photochemistry of DNA and its components. By combining different quantum mechanical methods, we investigated the molecular processes responsible for the damage of the genetic code or, on the contrary, those preventing it by dissipating the excess energy deposited in the system by UV absorption. Independently of its fundamental biological role, DNA, with its fluctuating closely stacked bases stabilized by weak nonbonding interactions, can be considered a prototypical MCA. Therefore, it allows one to tackle within a single system many of the conceptual and methodological challenges involved in the study of photoinduced processes in MCA.In this Account, by using the outcome of our studies on oligonucleotides as a guideline, we thus highlight the most critical modellistic issues to be faced when studying, either experimentally or computationally, the interaction between UV light and DNA and, at the same time, bring out their general relevance for the study of MCAs.We first discuss the rich photoactivated dynamics of nucleobases (the chromophores), highlighting the main effects modulating the interplay between their excited states and how the latter can affect the photoactivated dynamics of the polynucleotides, either providing effective monomer-like nonradiative decay routes or triggering reactive processes (e.g., triplet generation).We then tackle the reaction paths involving multiple bases, showing that in the DNA duplex the most important ones involve two stacked bases, forming a neutral excimer or a charge transfer (CT) state, which exhibit different spectral signatures and photochemical reactivity. In particular, we analyze the factors affecting the dynamic equilibrium between the excimer and CT, such as the fluctuations of the backbone or the rearrangement of the solvent.Next, we highlight the importance of the effects not directly connected to the chromophores, such as the flexibility of the backbone or the solvent effect. The former, affecting the stacking geometry of the bases, can determine the preference between different deactivation paths. The latter is particularly influential for CT states, making very important an accurate treatment of dynamical solvation effects, involving both the solvent bulk and specific solute-solvent interactions.In the last section, we describe the main methodological challenges related to the study of polynucleotide excited states and stress the benefits derived by the integration of complementary approaches, both computational and experimental. Only exploiting different point of views, in our opinion, it is possible to shed light on the complex phenomena triggered by light absorption in DNA, as in every MCA.

Quantum dynamics of the photoinduced charge separation in a symmetric donor–acceptor–donor triad: The role of vibronic couplings, symmetry and temperature

The photoinduced charge separation in a symmetric donor–acceptor–donor (D–A–D) triad is studied quantum mechanically using a realistic diabatic vibronic coupling model. The model includes a locally excited DA*D state and two charge-transfer states D+A−D and DA−D+ and is constructed according to a procedure generally applicable to semirigid D–A–D structures and based on energies, forces, and force constants obtained by quantum chemical calculations. In this case, the electronic structure is described by time-dependent density functional theory, and the corrected linear response is used in conjunction with the polarizable continuum model to account for state-specific solvent effects. The multimode dynamics following the photoexcitation to the locally excited state are simulated by the hybrid Gaussian-multiconfigurational time-dependent Hartree method, and temperature effects are included using thermo field theory. The dynamics are connected to the transient absorption spectrum obtained in recent experiments, which is simulated and fully assigned from first principles. It is found that the charge separation is mediated by symmetry-breaking vibrations of relatively low frequency, which implies that temperature should be accounted for to obtain reliable estimates of the charge transfer rate.

The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory

The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.

Surface Hopping Dynamics on Vibronic Coupling Models

The simulation of photoinduced non-adiabatic dynamics is of great relevance in many scientific disciplines, ranging from physics and materials science to chemistry and biology. Upon light irradiation, different relaxation processes take place in which electronic and nuclear motion are intimately coupled. These are best described by the time-dependent molecular Schrödinger equation, but its solution poses fundamental practical challenges to contemporary theoretical chemistry. Two widely used and complementary approaches to this problem are multiconfigurational time-dependent Hartree (MCTDH) and trajectory surface hopping (SH). MCTDH is an accurate fully quantum-mechanical technique but often is feasible only in reduced dimensionality, in combination with approximate vibronic coupling (VC) Hamiltonians, or both (i.e., reduced-dimensional VC potentials). In contrast, SH is a quantum–classical technique that neglects most nuclear quantum effects but allows nuclear dynamics in full dimensionality by calculating potential energy surfaces on the fly. If nuclear quantum effects do not play a central role and a linear VC (LVC) Hamiltonian is appropriate—e.g., for stiff molecules that generally keep their conformation in the excited state—then it seems advantageous to combine the efficient LVC and SH techniques. In this Account, we describe how surface hopping based on an LVC Hamiltonian (SH/LVC)—as recently implemented in the SHARC surface hopping package—can provide an economical and automated approach to simulate non-adiabatic dynamics. First, we illustrate the potential of SH/LVC in a number of showcases, including intersystem crossing in SO₂, intra-Rydberg dynamics in acetone, and several photophysical studies on large transition-metal complexes, which would be much more demanding or impossible to perform with other methods. While all of the applications provide very useful insights into light-induced phenomena, they also hint at difficulties faced by the SH/LVC methodology that need to be addressed in the future. Second, we contend that the SH/LVC approach can be useful to benchmark SH itself. By the use of the same (LVC) potentials as MCTDH calculations have employed for decades and by relying on the efficiency of SH/LVC, it is possible to directly compare multiple SH test calculations with a MCTDH reference and ponder the accuracy of various correction algorithms behind the SH methodology, such as decoherence corrections or momentum rescaling schemes. Third, we demonstrate how the efficiency of SH/LVC can also be exploited to identify essential nuclear and electronic degrees of freedom to be employed in more accurate MCTDH calculations. Lastly, we show that SH/LVC is able to advance the development of SH protocols that can describe nuclear dynamics including explicit laser fields—a very challenging endeavor for trajectory-based schemes. To end, this Account compiles the typical costs of contemporary SH simulations, evidencing the great advantages of using parametrized potentials. The LVC model is a sleeping beauty that, kissed by SH, is fueling the field of excited-state molecular dynamics. We hope that this Account will stimulate future research in this direction, leveraging the advantages of the SH/VC schemes to larger extents and extending their applicability to uncharted territories.

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 (La and Lb) 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.

Electronic Circular Dichroism Spectra of DNA Quadruple Helices Studied by Molecular Dynamics Simulations and Excitonic Calculations including Charge Transfer States

We here investigate the Electronic Circular Dichroism (ECD) Spectra of two representative Guanine-rich sequences folded in a Quadruple helix (GQ), by using a recently developed fragment diabatisation based excitonic model (FrDEx). FrDEx can include charge transfer (CT) excited states and consider the effect of the surrounding monomers on the local excitations (LEs). When applied to different structures generated by molecular dynamics simulations on a fragment of the human telomeric sequence (Tel21/22), FrDEx provides spectra fully consistent with the experimental one and in good agreement with that provided by quantum mechanical (QM) method used for its parametrization, i.e., TD-M05-2X. We show that the ECD spectrum is moderately sensitive to the conformation adopted by the bases of the loops and more significantly to the thermal fluctuations of the Guanine tetrads. In particular, we show how changes in the overlap of the tetrads modulate the intensity of the ECD signal. We illustrate how this correlates with changes in the character of the excitonic states at the bottom of the La and Lb bands, with larger LE and CT involvement of bases that are more closely stacked. As an additional test, we utilised FrDEx to compute the ECD spectrum of the monomeric and dimeric forms of a GQ forming sequence T30695 (5′TGGGTGGGTGGGTGGG3′), i.e., a system containing up to 24 Guanine bases, and demonstrated the satisfactory reproduction of the experimental and QM reference results. This study provides new insights on the effects modulating the ECD spectra of GQs and, more generally, further validates FrDEx as an effective tool to predict and assign the spectra of closely stacked multichromophore systems.