The deactivation of π-stacked cytosine molecules following excitation by ultrashort laser pulses was studied using semiclassical dynamics simulations. Another deactivation channel was found to compete with a previously reported path that led to dimerization. For both pathways, the initial excited state was found to form a charge-separated neutral exciton state, which forms an excimer state by charge transfer. When the interbase distance becomes less than 3 Å, charge recombination occurs due to strong intermolecular interaction, ultimately leading to an avoided crossing. Results indicate that the C2–N1–C6–C5 and C2′–N1′–C6′–C5′ dihedral angles play a significant role in the vibronic coupling between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). Vibrational energy distribution determines the fate of the excimer at the avoided crossing. Higher-amplitude vibration of C5 or C6 atoms leads to a nonadiabatic transition to the electronic ground state (a photophysical pathway); otherwise, a chemical reaction leading to the formation of cyclobutane type dimer occurs as found in earlier studies. The S1 and S0 potential energy surfaces calculated at TD-DFT level and the simulated trajectories were found to be consistent with CASPT2 results.
1. Introduction
Electronic excitation of DNA nucleobases induced by solar ultraviolet light has been extensively studied since the 1960s because it can lead to the formation of harmful photoproducts such as cyclobutane pyrimidine dimers (CPDs) between adjacent pyrimidine bases within the same DNA strand [1–3]. The formation of CPDs is considered as the main cause of cell death, mutagenesis, and development of skin cancers [4]. However, the quantum yield of photolesion formation is very low, suggesting a highly efficient nonradiative deactivation of excited DNA to the ground state consistent with the photostability of nucleobases [3].
Electronically excited single nucleobases decay to the ground state primarily by ultrafast internal conversion in less than 1 picosecond. However, in oligo- and polynucleotides the decay of excited nucleobases is more complicated, (showing multiexponential processes with different time constants [5–14]) because of interbase interactions, including base stacking and base pairing, and base-solvent interactions [15–18]. Longer-time decay observed in single-stranded oligonucleotides, where base paring is absent, is attributed to an intrastrand excimer/exciplex state resulting from π-stacking of two adjacent nucleobases [13, 19–24].
The formation of a delocalized excimer state slows down ultrafast nonradiative decay [19–24]. The ab initio calculations at CASPT2/CASSCF level [25] for two stacked adenines of different arrangements favor a mechanism involving two possible decay pathways. In this mechanism, unstacked or poorly stacked pairs of bases relax to the ground state by an ultrafast internal conversion through the conical intersection (CI) between the lowest excited state and ground state of the monomer. On the other hand, an excimer state formed between intrastrand stacked bases exhibits a longer decay time to the ground state because of an energy barrier along the path toward the CI of the monomer. It has been suggested [23] that the excimer state has a charge-transfer character and the state lifetime is correlated with the energy required to transfer an electron from one base to its stacked neighbor.
The formation of excimers is important not only to understand the distinct photophysics of oligonucleotides and DNA but also to account for the intrinsic and distinct photoinduced reaction of cytosine and thymine, which form CPDs. It has been proposed [14] that the singlet excimer state is a precursor to the photodimerization of DNA bases and the excited state dimerization reaction occurs in competition with an internal conversion to the electronic ground state of nucleobases. A recent theoretical calculation [26, 27] for cytosine photodimer formation at the CASSCF/CASPT2 level suggests that the nonadiabatic photodimerization of two stacked cytosine molecules proceeds through a CI involving the lowest singlet excited state and the ground state. This conical intersection, (S1/S0)CI, occurs below the energy of excimer 1(C*C). In order to reach (S1/S0)CI from 1(C*C), the system has to surmount an energy barrier of about 0.2 eV. The formation of the stable excimer along the lowest singlet excited-state for stacked cytosine molecules may make the formation of cyclobutane cytosine dimer (C〈〉C) less efficient compared to the formation of cyclobutane thymine dimer (T〈〉T) from a stacked pair of thymine molecules, which do not form a stable excimer. For the latter, the favored mechanism involves a barrierless nonadiabatic mechanism for the photodimerization along a singlet excited-state through a (S1/S0)CI, which is the funnel for ultrafast nonadiabatic decay leading the formation of cyclobutane thymine dimer (T〈〉T). The (S1/S0)CI leads to two possible processes of ultrafast internal conversion to the ground state. One pathway involves essentially a chemical change (dimerization), the other (photophysical) leads to two separated monomer bases.
In this paper we present a semiclassical dynamics simulation study for the laser-induced excitation and deactivation of two π-stacked cytosine molecules. The simulation results provide detailed information on the dynamics from photon excitation to deactivation and are expected to be helpful in understanding this process.
2. Methodology
The SERID method is used to carry out dynamics simulations. This technique is described in detail elsewhere [28, 29]. In this methodology, the valence electrons are calculated by the time-dependent Schrödinger equation while both the radiation field and the motion of the nuclei are treated classically. The one-electron states are obtained for each time step by solving the time-dependent Schrödinger equation in a nonorthogonal basis:
(1)iℏ∂Ψj∂t=S-1·H·Ψj,
where S is the overlap matrix of the atomic orbitals. The laser pulse is characterized by the vector potential A, which is coupled to the Hamiltonian via the time-dependent Peierls substitution [30]. Consider
(2)Ηab(X-X′)=Hab0(X-X′)exp{iqℏcA·(X-X′)},
where Hab(X-X′) is the Hamiltonian matrix element for basis functions a and b on atoms at X and X′, respectively, and q=-e is the charge of the electron.
In SERID, forces acting on nuclei or ions are computed by the Ehrenfest equation:
(3)Mld2Xlαdt2=-12∑jΨj+·{∂H∂Xlα-iℏ∂S∂Xlα·∂∂t}·Ψj-∂Urep∂Xlα,
where Urep is effective nuclear-nuclear repulsive potential and Xlα=〈X^lα〉 is the expectation value of the time-dependent Heisenberg operator for the α coordinate of the nucleus labeled by l (with α=x,y,z). Equation (3) is obtained by neglecting the second and higher order terms of the quantum fluctuations X^-〈X^lα〉 in the exact Ehrenfest theorem.
The time-dependent Schrödinger equation (1) is solved by using a unitary algorithm obtained from the equation for the time evolution operator [31]. Equation (3) is numerically integrated with the velocity Verlet algorithm. A time step of 0.05 femtoseconds is used for this study and energy conservation was then found to hold better than 1 part in 106 in a one picosecond simulation at 298 K.
The present “Ehrenfest” principle is complementary to other methods based on different approximations, such as the full multiple spawning model developed by the Martínez group [32]. The limitation of this method is that the simulation trajectory moves along a path dominated by averaging over all the terms in the Born-Oppenheimer expansion:
(4)Ψtotal(Xn,xe,t)=∑iΨin(Xn,t)Ψie(xe,t),
rather than following the time evolution of a single potential energy surface, which is approximately decoupled from all the others. (Here, Xn and xe represent the sets of nuclear and electronic coordinates, resp., and the Ψie are eigenstates of the electronic Hamiltonian at fixed Xn.) The strengths of the present approach include the retention of all of the 3N nuclear degrees of freedom and incorporation of both the excitation due to a laser pulse and the subsequent deexcitation at an avoided crossing near a CI.
3. Results and Discussions
In a recent paper [33], we described simulations that suggest the formation of cyclobutane cytosine dimer (C〈〉C) from a π-stacked pair of cytosine molecules irradiated by an ultrashort laser pulse (4.1 eV, fluence = 83.62 J·m−2, fwhm = 25 fs). Formation of an exciton intermediate upon excitation, followed by excimer formation, was involved in this dimerization channel. This paper describes another deactivation pathway found from additional simulations involving the same initial configuration. A typical trajectory (frequency of 4.1 eV, laser duration of 50 fs and fluence of 77.19 J·m2) is detailed below.
The variations with time of the HOMO-1, HOMO, LUMO, and LUMO+1 energies and the time-dependent population of these four orbitals are shown in Figures 1(a) and 1(b), respectively. Figure 1(a) shows that there is an abrupt fall in the LUMO energy soon after the application of the laser pulse. An intersection, which named avoided crossing, between the HOMO and LUMO levels is induced by coupling of orbitals, with an energy gap of 0.02 eV, is found at 976 fs. Figure 1(b) shows that by the end of the laser pulse (at 50 fs), about 1.4 electrons are excited from the HOMO to the LUMO, implying electronic excitation of one cytosine molecule. The coupling between the HOMO and LUMO, observed at 976 fs in Figure 1(a), suggests electronic transition from the LUMO to the HOMO. This deexcitation ultimately brings the molecules to the electronic ground state. It can be seen from Figure 1(a) that, shortly after the coupling, both the LUMO and HOMO levels move toward their initial values. After the electronic transition, these two energy levels slightly fluctuate about constants values which are essentially the same as their initial values. The excited state lifetime of 920 fs is about 3~4 times that which is for nonadiabatic deactivation of cytosine monomer [34] and is consistent with dimerization [33].
(a) The variations with time (b) the time-dependent populations of the HOMO-1, HOMO LUMO, and LUMO+1 energies of two stacked cytosine molecules.
The variations with time of the lengths of the C5–C6, C5′–C6′, C5–C5′ and C6–C5′ bonds are shown in Figure 2. Both C5–C6 and C5′–C6′ are double bonds at the beginning of simulation. Starting at 1.38 Å, which is the length of a typical C–C double bond in a conjugated aromatic system, the C5–C6 bond length rapidly elongates to about 1.5 Å after the laser pulse is applied while the C5′–C6′ bond length remains essentially constant. This suggests excitation of the C molecule, but not the C' molecule. After 1000 fs, both C5–C6 and C5′–C6′ bond lengths shorten to their initial values and remain at this length until the end of the simulation. In Figure 2(b), both C5–C5′ and C6–C6′ distances are more or less constant up to 600 fs. The C5–C5′ and C6–C6′ distances sharply fall at around 700 fs, both reaching a minimum value of 2.0 Å at about 1000 fs. The minima suggest that the C〈〉C is not formed in this trajectory. The structures inserted in Figure 2(b) represent the average geometries at different stages of the trajectory. The timescale from 100 to 700 fs is defined as stage I. In this stage, weaker intermolecular interaction is observed because of the relatively large distance (>3 Å) between the C and C' molecules. Excitation distorts the C molecule (bottom); note that the unexcited molecule (C') remains planar. With the shortening intermolecular distance, the C′ molecule deforms as the excimer evolves to stage II (900–1000 fs). After 1000 fs (stage III), the two cytosine molecules move far apart due to deactivation.
(a) The variations with time of the lengths between the C5 and C6 atoms and the C5′ and C6′ atoms. (b) The variations with time of the lengths between the C5 and C5′ atoms and the C6 and C6′ atoms in two stacked cytosines molecules.
Figure 2(a), stage I, shows previously unobserved variations in the C5–C6 and C5′–C6′ bond lengths. The C5–C6 bond length shortens from 1.53 Å to 1.48 Å; the C5′–C6′ bond length stretches from 1.38 Å to 1.42 Å. These variations correlate with intermolecular charge transfer, as seen in Figure 3. Before 100 fs (Figure 3(b)), the relatively flat curve for the net charge on C molecule suggests no intermolecular charge transfer. Thereafter, about 1.4 electrons transfer from C′ to C (Figure 3(a)), after about 900 fs (when interbase distance is 3.1 Å), charge recombination drives electrons back C′ and the two molecules return to electrically neutral states at the avoided crossing (976 fs). When electrons migrate to C, the increasing electronic density between C5 and C6 results in a reduced C5–C6 bond length because of the increasing C5–C6 bond order. Conversely, decreasing electronic density in C′ due to the weakening of C5′–C6′π bond is manifested by an increase in the C5′–C6′ bond length.
Variations with time of the net charge on C and the valence electron numbers of C5 and C6 atoms.
The variations in the valence electrons of the C5 and C6 atoms, which are most actively involved in the excitation-deactivation process of pyridines, are also shown in Figure 3; Figure 3(d) zooms in on the 0 to 100 fs region of Figure 3(c). The valence electrons on the atoms are calculated by projection of single electronic wave functions to the orbitals of these atoms. Because of bond polarities due heteroatoms in the aromatic nucleus, valence electrons number of carbon atoms usually fluctuate about the 4 instead of being exactly equal to 4. For example, the valence electrons number of C5 and C6 are 4.31 and 3.86, respectively. Excitation by laser results in redistribution of π electrons. Consequently, the covalence number for the C5 atom drops from 4.31 to 3.9, while increasing to an average of about 4.5 for the C6 atom during the first 100 fs. Since there is minimal electron transfer between the two cytosine molecules during this time period (Figure 3(b)), this indicates a charge separation within the excited C molecule. The charge separated neutral excited state with a lifetime of 100 fs is very similar to the “Frenkel exciton states” [14] formed by the coupling of 1ππ* states (which are localized on single bases) of proximal nucleobases proposed by Gustavsson and Markovitsi and coworkers [15–18]. According to some experimental and theoretical research [23, 35–37], exciton formation tends to precede the formation of an excimer/exciplex. The shortening of C5–C6 and elongation of C5′–C6′, as shown in Figure 2(a), suggest an exciton state evolving to a charge transfer state (i.e., excimer) after 100 fs.
The variation of dihedral angles, C2–N1–C6–C5 and C2′–N1′–C6′–C5′, which describe the twisting of the C5–C6 and C5′–C6′ bonds, respectively, is compared to the HOMO-LUMO energy gap in Figure 4. The dihedral C2–N1–C6–C5 decreases from 0° to −27° within 120 fs, then fluctuates around an average value of −20°. It indicates that the C molecule has been deformed due to excitation. Greater torsion results in a smaller HOMO-LUMO energy gap and a higher electron affinity for C, leading to an electron transfer from C′ to C. The C2–N1–C6–C5 angle reaches a maximum of −35° at 986 fs and returns to its initial value after that time. On the other hand, the dihedral C2′–N1′–C6′–C5′ varies little until 900 fs, then sharply rises to 15° at 986 fs. The variation of the C2′–N1′–C6′–C5′ dihedral angle suggests that interbase interaction leads to deformation of C′ after 900 fs. The strong torsion of C′ gives rise to nonadiabatic transition, which leads to the electronic ground state.
Variations with time of the dihedral angles C2–N1–C6–C5 and C2′–N1′–C6′–C5′ and energy gap between HOMO and LUMO.
The geometry at the avoided crossing is shown in Figure 5. The C2–N1–C6–C5 and C2′–N1′–C6′–C5′ dihedral angles reach their maximum values at 976 fs. This strong torsion causes a nonadiabatic transition to electronic ground state. For the molecular geometry taken at 976 fs, the C5–C5′ and C6–C6′ distances are 2.13 and 2.19 Å, respectively. These compare favorably with 2.27 and 2.17 Å obtained by CASSCF/CASPT2 calculation [26] at the CI between the lowest excited singlet state and the ground state. This indicates that the avoided crossing found at 976 fs is, indeed, in close proximity to the CI. The distances are also very similar to the dimerization channel that we previously reported. We believe that the avoided crossing may be identical for both pathways.
Geometries taken from trajectory at time of 976 fs, when avoided crossing occur.
It is worth noting that the average interbase distance is 3 Å at 900 fs as the C2′–N1′–C6′–C5′ dihedral angle reaches −30° and electrons transfer back. We obtained the same result in our previous work [33]. This suggests that 3 Å is a minimum value for both deactivation pathways. At larger distances, the excited system evolves to an excimer characterized by charge transfer; at shorter distances, the strong intermolecular interaction may change the redox potential of bases and result in charge recombination, which leads to two neutral molecules. We have previously studied the effect of interbase distance (from 3.4 to 4.4 Å) using SERID [38]. The results suggest that 4.2 Å is the upper limit for the deactivation pathways associated with π stacking.
Due to the substantial defects of the mean field theory used in SERID, the electronic states are not explicit. The potential energy surfaces for the S1 and S0 states were calculated by using long-range corrected time-dependent density functional (LRC-TD-DFT) theory [39–42] for comparison and the results are summarized in Figure 6. The Cam-B3LYP function, which was reported to describe charge transfer state well [43–45], was used in this research. In this calculation, energies of serial single points taken along the simulation trajectory were calculated. A CI with energy value of 2.5 eV was observed at 990 fs in the TDDFT potential energy surface, which is very close to avoided crossing in our simulation. The non-optimized initial ground state geometry results in a “quasi vertical excited energy” of 3.2 eV, which is much lower than CASPT2 result of 4.41 eV [26]. However, the TD-DFT and CASPT2 calculations still yield comparable values for the energy gap between the Frank-Condon point and CI (0.7 eV and 0.9 eV, resp.). Based on [26], we believe that an optimized ground state will be located at −1.21 eV level in the current coordinate system. Therefore, the energy of the CI will be 3.71 eV, which agrees well with 3.5 eV from CASPT2. In [26], the CASPT2 calculation also yields a minimum in S1 that is lower than the CI by 0.2 eV, as same as our TD-DFT results. Our simulation indicates that from the minimum, the (C*C) excimer has to overcome a barrier of 0.9 eV to reach the funnel. The stability of the excimer may decrease the effectiveness of internal conversion in the singlet manifold compared to the thymine excimer. This explains why stacked thymines have lifetime of (about 600 fs [38]) which is shorter than stacked cytosine; internal conversion of the thymine excimer goes through a barrierless pathway. The geometry of the CI is shown in Figure 5. The C5–C5′ and C6–C6′ are 2.13 and 2.19 Å respectively, compare to 2.27 and 2.17 Å respectively in results of CASPT2. The similar geometries suggest that this semiclassical dynamics simulation is credible.
PES of stacked cytosines system at Cam-B3LYP/6-311+g* level. Geometry of conical intersection is inserted.
In our previous study [33], π-stacked cytosines could form CPD by laser pulse irradiation. In the dimerization mechanism, an evolution of “ultrafast exciton” to “charge transfer excimer” and ultimately leading to formation of CPD was reported. It is hard to understand that the almost identical approach in this work follows a photophysical deactivation mechanism. In order to distinguish dimerization and photophysical deactivation mechanism of π-stacked cytosines, we counted multiple trajectories simulations results. Different initial stacking orientations, with average intermolecular distance of 3.4 Å and dihedral C5–C6–C6′–C5′ of 36°, were used in multiple trajectories simulations. Simulations gave rise to both dimerization and photophysical deactivation pathway. To compare the vibrational degree of freedom, seven dimerization pathways and seven excimer pathways were chosen and several geometrical parameters with respect to central C5–C6 bond, namely, the distance of C6–N1 and C6–H6 and dihedral angles of C2–N1–C6–C5 and C4–C5–C6–H6, are presented in Figure 7. The red solid lines represent dimerization channels and blue dashed lines represent the photophysical deactivations channels., respectively. Excited state lifetimes of the all fourteen trajectories are approximately 350 fs because of the relatively short intermolecular distance [38]. All of the four aforementioned degrees of freedom have larger vibrational amplitudes in the photophysical channel before nonadiabatic transition. It suggests that vibration decreases possibility of bonding between C5 or C6 atoms.
The Variations of some degrees of freedom around C5 and C6 atoms of dimerization channels and photophysical deactivation channels. The two channels are simplified as “dimer” and “excimer,” respectively.
Obviously, through the CI in PES, both evolution forward to photoproducts and back to reactants are possible. Therefore, two competitive channels exist in π-stacked cytosines system. Figure 8 summarizes the changes in molecular geometry for representative trajectories from the competing pathways. Shortly after application of the laser pulse (<100 fs), a neutral dipolar coupling excited state with ultrashort lifetime, an exciton, is formed. After 120 fs, the twisting of dihedral angle C2–N1–C6–C5 increases the electron affinity of C and leads to a transfer of 1.4 electrons from C′ to C. As a result, an charge transfer excimer is formed. When two bases approach to within a distance of about 3 Å, charge recombination occurs; the stacked cytosines revert to neutral. The system evolves to an avoided crossing due to vibrational coupling of the HOMO and LUMO induced by maximal deformation of C and C′ molecules. At avoided crossing, vibrational energy distribution determines the deactivation fate.
Mechanism of photodimerization reaction of π-stacked cytosines system. The abbreviations are interpreted as follows: CT, charge transfer; CR, charge recombination; AC, avoided crossing; gs, ground state.
4. Conclusion
This paper presents simulation results for the photophysical deactivation of the π-stacked cytosine system following excitation subjected by a laser pulse with 4.1 eV photon energy. Only one cytosine is excited at the beginning. Excitation initially leads to the formation of a charge separated neutral exciton state with a lifetime of 100 fs. An “excimer state” results from charge transfer that occurs when the interbase distance is 3.8 Å in this simulation. Charge recombination occurs when the interbase distance shortens to about 3 Å, leading to neutral stacked bases that evolve into an avoided crossing prior to deactivation. Geometries taken from avoided crossing are very similar to the S1/S0 CI obtained by CASSCT/CASPT2 [40, 41]. Torsional vibrations (C2–N1–C6–C5 and C2′–N1′–C6′–C5′) play a significant role in vibronic coupling between the HOMO and LUMO and lead to nonadiabatic transition of the molecules to the ground state.
The geometries of avoided crossing are almost the same as those for dimerization channel reported previously. It indicates there are two competitive de-exciton channels at avoided crossing. Higher amplitude vibrations of C5 or C6 atoms favor the photophysical deactivation pathway.
To compare with other theoretical results, S1 and S0 PESs were calculated at TD-DFT level along with simulation trajectory. The CI was found to be very near the avoided crossing obtained by SERID method. The TD-DFT calculation yields a 0.2 eV barrier between the minimum of S1 and CI, indicating a relatively stable excimer in π-stacked cytosines system, which lessens the effectiveness of internal conversion compared to thymine excimer, which has a barrierless pathway from S1 to CI.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the Joint Funds of the National Natural Science Foundation of China (Grant No. U1330138) and the Natural Science Foundation Project of CQ CSTC (cstc2011jjA00009) and Project of the Science Technology Foundation of Chongqing Education Committee, China (no. KJ120516 and KJ120521).
CadetJ.SageE.DoukiT.Ultraviolet radiation-mediated damage to cellular DNA20055711-231710.1016/j.mrfmmm.2004.09.0122-s2.0-14644433816BeukersR.EkerA. P. M.LohmanP. H. M.50 years thymine dimer20087353054310.1016/j.dnarep.2007.11.0102-s2.0-39149121838MouretS.BaudouinC.CharveronM.FavierA.CadetJ.DoukiT.Cyclobutane pyrimidine dimers are predominant DNA lesions in whole human skin exposed to UVA radiation200610337137651377010.1073/pnas.06042131032-s2.0-33748784343MelnikovaV. O.AnanthaswamyH. N.Cellular and molecular events leading to the development of skin cancer20055711-29110610.1016/j.mrfmmm.2004.11.0152-s2.0-14644405691McCullaghM.HariharanM.LewisF. D.MarkovitsiD.DoukiT.SchatzG. C.Conformational control of TT dimerization in DNA conjugates: a molecular dynamics study2010114155215522110.1021/jp100983t2-s2.0-77951136861Serrano-AndrésL.MerchánM.BorinA. C.A three-state model for the photophysics of adenine200612256559657110.1002/chem.2005015152-s2.0-33748596860PecourtJ. M.PeonJ.KohlerB.DNA excited-state dynamics: ultrafast internal conversion and vibrational cooling in a series of nucleosides200112342103701037810.1021/ja01614532-s2.0-0035944466IsmailN.BlancafortL.OlivucciM.KohlerB.RobbM. A.Ultrafast decay of electronically excited singlet cytosine via a π,π* to n0,π* state switch2002124246818681910.1021/ja02582732-s2.0-0037134912LangerH.DoltsinisN. L.MarxD.Excited-state dynamics and coupled proton-electron transfer of guanine2005691734173710.1002/cphc.2004005782-s2.0-24944446297BarbattiM.LischkaH.Can the nonadiabatic photodynamics of aminopyrimidine be a model for the ultrafast deactivation of adenine?2007111152852285810.1021/jp070089w2-s2.0-34247582809HudockH. R.LevineB. G.ThompsonA. L.SatzgerH.TownsendD.GadorN.UllrichS.StolowA.MartínezT. J.Ab Initio molecular dynamics and time-resolved photoelectron spectroscopy of electronically excited uracil and thymine2007111348500850810.1021/jp07236652-s2.0-34548572225KwokW.MaC.PhillipsD. L.Femtosecond time- and wavelength-resolved fluorescence and absorption spectroscopic study of the excited states of adenosine and an adenine oligomer200612836118941190510.1021/ja06220022-s2.0-33748493389Crespo-HernándezC. E.KohlerB.Influence of secondary structure on electronic energy relaxation in adenine homopolymers200410830111821118810.1021/jp04960462-s2.0-4043151600MiddletonC. T.de La HarpeK.SuC.LawY. K.Crespo-HernándezC. E.KohlerB.DNA excited-state dynamics: from single bases to the double helix20096021723910.1146/annurev.physchem.59.032607.0937192-s2.0-64549150256BouvierB.DognonJ. P.LaveryR.MarkovitsiD.MilliéP.OnidasD.ZakrzewskaK.Influence of conformational dynamics on the exciton states of DNA oligomers200310748135121352210.1021/jp036164u2-s2.0-0348207686EmanueleE.MarkovitsiD.MilliéP.ZakrzewskaK.UV spectra and excitation delocalization in DNA: influence of the spectral width2005671387139210.1002/cphc.2005000142-s2.0-22444433057OnidasD.GustavssonT.LazzarottoE.MarkovitsiD.Fluorescence of the DNA double helices (dAdT)n·(dAdT) n studied by femtosecond spectroscopy20079375143514810.1039/b707882j2-s2.0-34548810221BuchvarovI.WangQ.RaytchevM.TrifonovA.FiebigT.Electronic energy delocalization and dissipation in single- and double-stranded DNA2007104124794479710.1073/pnas.06067571042-s2.0-34247642158PanZ. Z.McCullaghM.SchatzG. C.LewisF. D.Conformational control of thymine photodimerization in purine-containing trinucleotides20112121432143810.1021/jz200532w2-s2.0-79959387410ZgierskiM. Z.FujiwaraT.LimE. C.Conical intersections and ultrafast intramolecular excited-state dynamics in nucleic acid bases and electron donor-acceptor molecules20084634-628929910.1016/j.cplett.2008.06.0242-s2.0-52049092514Crespo-HernándezC. E.CohenB.KohlerB.Molecular spectroscopy: complexity of excited-state dynamics in DNA (Reply)2006441artcile E810.1038/nature049042-s2.0-33744991779CohenB.HareP. M.KohlerB.Ultrafast excited-state dynamics of adenine and monomethylated adenines in solution: implications for the nonradiative decay mechanism200312544135941360110.1021/ja035628z2-s2.0-0242362228TakayaT.SuC.de la HarpeK.Crespo-HernándezC. E.KohlerB.UV excitation of single DNA and RNA strands produces high yields of exciplex states between two stacked bases200810530102851029010.1073/pnas.08020791052-s2.0-48749126017Crespo-HernándezC. E.CohenB.KohlerB.Base stacking controls excited-state dynamics in A·T DNA200543670541141114410.1038/nature039332-s2.0-24144459117Olaso-GonzálezG.MerchánM.Serrano-AndrésL.The role of adenine excimers in the photophysics of oligonucleotides2009131124368437710.1021/ja808280j2-s2.0-67949103470Roca-SanjuánD.Olaso-GonzálezG.González-RamírezI.Serrano-AndrésL.MerchánM.Molecular basis of DNA photodimerization: intrinsic production of cyclobutane cytosine dimers200813032107681077910.1021/ja803068n2-s2.0-49549088826Serrano-PérezJ. J.González-RamírezI.CotoP. B.MerchánM.Serrano-AndrésL.Theoretical insight into the intrinsic ultrafast formation of cyclobutane pyrimidine dimers in UV-irradiated DNA: thymine versus cytosine200811245140961409810.1021/jp806794x2-s2.0-57149108311DouY.TorralvaB. R.AllenR. E.Semiclassical electron-radiation-ion dynamics (SERID) and cis-trans photoisomerization of butadiene20035015–172615264310.1080/09500340320001207592-s2.0-0242456295DouY.TorralvaB. R.AllenR. E.Interplay of electronic and nuclear degrees of freedom in a femtosecond-scale photochemical reaction20043924–635235710.1016/j.cplett.2004.05.0872-s2.0-3042679776GrafM.VoglP.Electromagnetic fields and dielectric response in empirical tight-binding theory19955184940494910.1103/PhysRevB.51.49402-s2.0-0000067185AllenR. E.DumitricaT.TorralvaB. R.TsenK. T.Electronic and structural response of materials to fast, intense laser pulses2001chapter 7New York, NY, USAAcademic PressBen-NunM.MartínezT. J.AB initio quantum molecular dynamics200212143951210.1002/0471264318.ch72-s2.0-33947245222YuanS.ZhangW.LiuL.DouY.FangW.LoG. V.Detailed mechanism for photoinduced cytosine dimerization: a semiclassical dynamics simulation201111546132911329710.1021/jp207550a2-s2.0-81555209065SharonovA.GustavssonT.CarréV.RenaultE.MarkovitsiD.Cytosine excited state dynamics studied by femtosecond fluorescence upconversion and transient absorption spectroscopy20033801-217318010.1016/j.cplett.2003.09.0212-s2.0-0141768256KohlerB.Nonradiative decay mechanisms in DNA model systems20101132047205310.1021/jz100491x2-s2.0-77954255381BittnerE. R.Lattice theory of ultrafast excitonic and charge-transfer dynamics in DNA200612591209490910.1063/1.23354522-s2.0-33748582019ImprotaR.SantoroF.BaroneV.LamiA.Vibronic model for the quantum dynamical study of the competition between bright and charge-transfer excited states in single-strand polynucleotides: the adenine dimer case200911352153461535410.1021/jp906278t2-s2.0-73949102640ZhangW.2011Xi'an, ChinaNorthwest UniversityIikuraH.TsunedaT.YanaiT.HiraoK.A long-range correction scheme for generalized-gradient-approximation exchange functionals200111583540354410.1063/1.13835872-s2.0-0035934184SongJ. W.HirosawaT.TsunedaT.HiraoK.Long-range corrected density functional calculations of chemical reactions: redetermination of parameter20071261515410510.1063/1.27215322-s2.0-34247355539VydrovO. A.ScuseriaG. E.Intermittent search process and teleportation20071262323410910.1063/1.2741516HendersonT. M.JaneskoB. G.ScuseriaG. E.Generalized gradient approximation model exchange holes for range-separated hybrids20081281919410510.1063/1.29217972-s2.0-44349094803SantoroF.BaroneV.ImprotaR.Can TD-DFT calculations accurately describe the excited states behavior of stacked nucleobases? The cytosine dimer as a test case200829695796410.1002/jcc.208532-s2.0-41549147727SantoroF.BaroneV.ImprotaR.Excited states decay of the A-T DNA: a PCM/TD-DFT study in aqueous solution of the (9-methyl-adenine)2·(1-methyl-thymine)2 stacked tetramer20091314215232152452-s2.0-7035030720410.1021/ja904777hLangeA. W.HerbertJ. M.Both intra- And interstrand charge-transfer excited states in aqueous B-DNA are present at energies comparable to, or just above, the 1ππ* excitonic bright states2009131113913392210.1021/ja808998q2-s2.0-67849114552