Competing deactivation channels for excited [pi]-stacked cytosines.
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 . The formation of CPDs is considered as the main cause of cell death, mutagenesis, and development of skin cancers . 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 .
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 re-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  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  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  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, ([S.sub.1]/[[S.sub.0]).sub.CI], occurs below the energy of excimer [sup.1](C*C). In order to reach ([S.sub.1]/[[S.sub.0]).sub.CI] from [sup.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 ([S.sub.1]/[[S.sub.0]).sub.CI], which is the funnel for ultrafast nonadiabatic decay leading the formation of cyclobutane thymine dimer (T<>T). The ([S.sub.1]/[[S.sub.0]).sub.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 [pi]-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.
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 Schrodinger 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 Schroodinger equation in a nonorthogonal basis:
1[??][[partial derivative][[psi].sub.j]]/[partial derivative]t = [S.sup.-1] x H x [[psi].sub.j], (1)
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 . Consider
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)
where [H.sub.ab](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:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)
where [U.sub.rep] is effective nuclear-nuclear repulsive potential and [X.sub.la] = ([[??].sub.la]) is the expectation value of the time-dependent Heisenberg operator for the a coordinate of the nucleus labeled by l (with [alpha] = x, y, z). Equation (3) is obtained by neglecting the second and higher order terms of the quantum fluctuations [??] - <[[??].sub.l[alpha]]> in the exact Ehrenfest theorem.
The time-dependent Schrodinger 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 [10.sup.6] 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 Martinez group . 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:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)
rather than following the time evolution of a single potential energy surface, which is approximately decoupled from all the others. (Here, [X.sub.n] and [x.sub.e] represent the sets of nuclear and electronic coordinates, resp., and the [[psi].sup.e.sub.i] are eigenstates of the electronic Hamiltonian at fixed [X.sub.n].) 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 , we described simulations that suggest the formation of cyclobutane cytosine dimer (C<>]C) from a [pi]-stacked pair of cytosine molecules irradiated by an ultrashort laser pulse (4.1 eV, fluence = 83.62 J x [m.sup.-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 x [m.sup.2]) 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), suggest selectronic 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  and is consistent with dimerization .
The variations with time of the lengths of the [C.sub.5]-[C.sub.6], [C'.sub.5]-[C'.sub.6], [C.sub.5]-[C'.sub.5] and [C.sub.6]-[C'.sub.5] bonds are shown in Figure 2. Both [C.sub.5]-[C.sub.6] and [C'.sub.5]-[C'.sub.6] are double bonds at the beginning of simulation. Starting at 1.38 [Angstrom], which is the length of a typical C-C double bond in a conjugated aromatic system, the [C.sub.5]-[C.sub.6] bond length rapidly elongates to about 1.5 [Angstrom] after the laser pulse is applied while the [C'.sub.5]-[C'.sub.6] bond length remains essentially constant. This suggests excitation of the C molecule, but not the C molecule. After 1000 fs, both [C.sub.5]-[C.sub.6] and [C'.sub.5]-[C'.sub.6] bond lengths shorten to their initial values and remain at this length until the end of the simulation. In Figure 2(b), both [C.sub.5]-[C'.sub.5] and [C.sub.6]-[C'.sub.6] distances are more or less constant up to 600 fs. The [C.sub.5]-[C'.sub.5] and [C.sub.6]-[C'.sub.6] distances sharply fall at around 700 fs, both reaching a minimum value of 2.0 [Angstrom] 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[Angstrom]) 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.
Figure 2(a), stage I, shows previously unobserved variations in the [C.sub.5]-[C.sub.6] and [C'.sub.5]-[C'.sub.6] bond lengths. The [C.sub.5]-[C.sub.6] bond length shortens from 1.53 [Angstrom] to 1.48 [Angstrom]; the [C.sub.5]-[C.sub.6] bond length stretches from 1.38 A to 1.42 A. 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[Angstrom]), 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 [C.sub.5] and [C.sub.6] results in a reduced [C.sub.5]-[C.sub.6] bond length because of the increasing [C.sub.5]-[C.sub.6] bond order. Conversely, decreasing electronic density in C' due to the weakening of [C'.sub.5]-[C'.sub.6][pi] bond is manifested by an increase in the [C.sub.5]-[C.sub.6] bond length.
The variations in the valence electrons of the [C.sub.5] and [C.sub.6] 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 [C.sub.5] and [C.sub.6] are 4.31 and 3.86, respectively. Excitation by laser results in redistribution of [pi] electrons. Consequently, the covalence number for the [C.sub.5] atom drops from 4.31 to 3.9, while increasing to an average of about 4.5 for the [C.sub.6] 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"  formed by the coupling of [sup.1][pi][[pi].sup.*] 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 [C.sub.5]-[C.sub.6] and elongation of [C'.sub.5]-[C'.sub.6], 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, [C.sub.2]-[N.sub.1]-[C.sub.6]-[C.sub.5] and [C'.sub.2]-[N'.sub.1]-[C'.sub.6]-[C'.sub.5], which describe the twisting of the [C.sub.5]-[C.sub.6] and [C'.sub.5]-[C'.sub.6] bonds, respectively, is compared to the HOMO-LUMO energy gap in Figure 4. The dihedral [C.sub.2]-[N.sub.1]-[C.sub.6]-[C.sub.5] decreases from 0[degrees] to -27[degrees] within 120 fs, then fluctuates around an average value of -20[degrees]. 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'.sub.2] to C. The [C.sub.2]-[N.sub.1]-[C.sub.6]-[C.sub.5] angle reaches a maximum of-35[degrees] at 986 fs and returns to its initial value after that time. On the other hand, the dihedral [C'.sub.2]-[N'.sub.1]-[C'.sub.6]-[C'.sub.5] varies little until 900 fs, then sharply rises to 15[degrees] at 986 fs. The variation of the [C'.sub.2]-[N'.sub.1]-[C'.sub.6]-[C'.sub.5] 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.
The geometry at the avoided crossing is shown in Figure 5. The [C.sub.2]-Nx-[C.sub.6]-[C.sub.5] and [C'.sub.2]-[N'.sub.1]-[C'.sub.6]-[C'.sub.5] 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 [C.sub.5]-[C'.sub.5] and [C.sub.6]-[C'.sub.6] distances are 2.13 and 2.19 A, respectively. These compare favorably with 2.27 and 2.17 A obtained by CASSCF/CASPT2 calculation  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 maybe identical for both pathways.
It is worth noting that the average interbase distance is 3 [Angstrom] at 900 fs as the [C'.sub.2]-[N'.sub.1]-[C'.sub.6]-[C.sub.5]-dihedral angle reaches -30[degrees] and electrons transfer back. We obtained the same result in our previous work . This suggests that 3 [Angstrom] 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 [Angstrom]) using SERID . The results suggest that 4.2 A is the upper limit for the deactivation pathways associated with [pi] 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 [S.sub.1] and [S.sub.0] 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 CamB3LYP 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 .
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 , 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 , the CASPT2 calculation also yields a minimum in [S.sub.1] 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 ) 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 [C.sub.5]-[C'.sub.5] and [C.sub.6]-[C'.sub.6] are 2.13 and 2.19 [Angstrom] respectively, compare to 2.27 and 2.17 [Angstrom] respectively in results of CASPT2. The similar geometries suggest that this semiclassical dynamics simulation is credible.
In our previous study , [pi][pi]-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 [pi]-stacked cytosines, we counted multiple trajectories simulations results. Different initial stacking orientations, with average intermolecular distance of 3.4 [Angstrom] and dihedral [C.sub.5]-[C.sub.6]-[C'.sub.6]-[C'.sub.5] of 36[degrees], 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 [C.sub.5]-[C.sub.6] bond, namely, the distance of [C.sub.6]-Nx and [C.sub.6]-H6 and dihedral angles of [C.sub.2]-[N.sub.1]-[C.sub.6]-[C.sub.5] and [C.sub.4]-[C.sub.5]-[C.sub.6]-[H.sub.6], 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 . 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 [C.sub.5] or [C.sub.6] atoms.
Obviously, through the CI in PES, both evolution forward to photoproducts and back to reactants are possible. Therefore, two competitive channels exist in [pi]-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 [C.sub.2]-[N.sub.1]-[C.sub.6]-[C.sub.5] 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 [Angstrom], 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.
This paper presents simulation results for the photophysical deactivation of the [pi]-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 A in this simulation. Charge recombination occurs when the interbase distance shortens to about 3 A, leading to neutral stacked bases that evolve into an avoided crossing prior to deactivation. Geometries taken from avoided crossing are very similar to the [S.sub.1] /[S.sub.0] CI obtained by CASSCT/CASPT2 [40, 41]. Torsional vibrations ([C.sub.2]-[N.sub.1] -[C.sub.6]-[C.sub.5] and [C'.sub.2]-[N'.sub.1]-[C'.sub.6]-[C'.sub.5]) 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 [C.sub.5] or [C.sub.6] atoms favor the photophysical deactivation pathway.
To compare with other theoretical results, [S.sub.1] and [S.sub.0] 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 [S.sub.1] and CI, indicating a relatively stable excimer in [pi]-stacked cytosines system, which lessens the effectiveness of internal conversion compared to thymine excimer, which has a barrierless pathway from [S.sub.1] to CI.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
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).
 J. Cadet, E. Sage, and T. Douki, "Ultraviolet radiation-mediated damage to cellular DNA," Mutation Research/Fundamental and Molecular Mechanisms of Mutagenesis, vol. 571, no. 1-2, pp. 3-17, 2005.
 R. Beukers, A. P. M. Eker, and P. H. M. Lohman, "50 years thymine dimer," DNA Repair, vol. 7, no. 3, pp. 530-543, 2008.
 S. Mouret, C. Baudouin, M. Charveron, A. Favier, J. Cadet, and T. Douki, "Cyclobutane pyrimidine dimers are predominant DNA lesions in whole human skin exposed to UVA radiation,"
Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 37, pp. 13765-13770, 2006.
 V. O. Melnikova and H. N. Ananthaswamy, "Cellular and molecular events leading to the development of skin cancer," Mutation Research/Fundamental and Molecular Mechanisms of Mutagenesis, vol. 571, no. 1-2, pp. 91-106, 2005.
 M. McCullagh, M. Hariharan, F. D. Lewis, D. Markovitsi, T. Douki, and G. C. Schatz, "Conformational control of TT dimerization in DNA conjugates: a molecular dynamics study," Journal of Physical Chemistry B, vol. 114, no. 15, pp. 5215-5221, 2010.
 L. Serrano-Andres, M. Merchan, and A. C. Borin, "A three-state model for the photophysics of adenine," Chemistry: A European Journal, vol. 12, no. 25, pp. 6559-6571, 2006.
 J. M. Pecourt, J. Peon, and B. Kohler, "DNA excited-state dynamics: ultrafast internal conversion and vibrational cooling in a series of nucleosides," Journal of the American Chemical Society, vol. 123, no. 42, pp. 10370-10378, 2001.
 N. Ismail, L. Blancafort, M. Olivucci, B. Kohler, and M. A. Robb, "Ultrafast decay of electronically excited singlet cytosine via a [pi], [[pi].sup.*] to [n.sub.0], [[pi].sup.*] state switch," Journal of the American Chemical Society, vol. 124, no. 24, pp. 6818-6819, 2002.
 H. Langer, N. L. Doltsinis, and D. Marx, "Excited-state dynamics and coupled proton-electron transfer of guanine," ChemPhysChem, vol. 6, no. 9, pp. 1734-1737, 2005.
 M. Barbatti and H. Lischka, "Can the nonadiabatic photodynamics of aminopyrimidine be a model for the ultrafast deactivation of adenine?" Journal of Physical Chemistry A, vol. 111, no. 15, pp. 2852-2858, 2007.
 H. R. Hudock, B. G. Levine, A. L. Thompson et al., "Ab Initio molecular dynamics and time-resolved photoelectron spectroscopy of electronically excited uracil and thymine," Journal of Physical Chemistry A, vol. 111, no. 34, pp. 8500-8508, 2007
 W. Kwok, C. Ma, and D. L. Phillips, "Femtosecond timeand wavelength-resolved fluorescence and absorption spectroscopic study of the excited states of adenosine and an adenine oligomer," Journal of the American Chemical Society, vol. 128, no. 36, pp. 11894-11905, 2006.
 C. E. Crespo-Hernandez and B. Kohler, "Influence of secondary structure on electronic energy relaxation in adenine homopolymers," Journal of Physical ChemistryB, vol. 108, no. 30, pp. 11182-11188, 2004.
 C. T. Middleton, K. de La Harpe, C. Su, Y. K. Law, C. E. Crespo-Hernandez, and B. Kohler, "DNA excited-state dynamics: from single bases to the double helix," Annual Review of Physical Chemistry, vol. 60, pp. 217-239, 2009.
 B. Bouvier, J. P. Dognon, R. Lavery et al., "Influence of conformational dynamics on the exciton states of DNA oligomers," Journal of Physical ChemistryB, vol. 107, no. 48, pp. 13512-13522, 2003.
 E. Emanuele, D. Markovitsi, P. Millie, and K. Zakrzewska, "UV spectra and excitation delocalization in DNA: influence of the spectral width," ChemPhysChem, vol. 6, no. 7, pp. 1387-1392, 2005.
 D. Onidas, T. Gustavsson, E. Lazzarotto, and D. Markovitsi, "Fluorescence of the DNA double helices (dAdT)n-(dAdT) n studied by femtosecond spectroscopy," Physical Chemistry Chemical Physics, vol. 9, no. 37, pp. 5143-5148, 2007
 I. Buchvarov, Q. Wang, M. Raytchev, A. Trifonov, and T. Fiebig, "Electronic energy delocalization and dissipation in single- and double-stranded DNA," Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 12, pp. 4794-4797, 2007.
 Z. Z. Pan, M. McCullagh, G. C. Schatz, and F. D. Lewis, "Conformational control of thymine photo dimerization in purine-containing trinucleotides," The Journal of Physical Chemistry Letters, vol. 2, no. 12, pp. 1432-1438, 2011.
 M. Z. Zgierski, T. Fujiwara, and E. C. Lim, "Conical intersections and ultrafast intramolecular excited-state dynamics in nucleic acid bases and electron donor-acceptor molecules," Chemical Physics Letters, vol. 463, no. 4-6, pp. 289-299, 2008.
 C. E. Crespo-Hernandez, B. Cohen, and B. Kohler, "Molecular spectroscopy: complexity of excited-state dynamics in DNA (Reply)," Nature, vol. 441, artcile E8, 2006.
 B. Cohen, P. M. Hare, and B. Kohler, "Ultrafast excited-state dynamics of adenine and monomethylated adenines in solution: implications for the nonradiative decay mechanism," Journal of the American Chemical Society, vol. 125, no. 44, pp. 13594-13601, 2003.
 T. Takaya, C. Su, K. de la Harpe, C. E. Crespo-Hernandez, and B. Kohler, "UV excitation of single DNA and RNA strands produces high yields of exciplex states between two stacked bases," Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 30, pp. 10285-10290, 2008.
 C. E. Crespo-Hernandez, B. Cohen, and B. Kohler, "Base stacking controls excited-state dynamics in A-T DNA," Nature, vol. 436, no. 7054, pp. 1141-1144, 2005.
 G. Olaso-Gonzalez, M. Merchan, and L. Serrano-Andres, "The role of adenine excimers in the photophysics of oligonucleotides," Journal of the American Chemical Society, vol. 131, no. 12, pp. 4368-4377, 2009.
 D. Roca-Sanjuan, G. Olaso-Gonzalez, I. Gonzalez-Ramirez, L. Serrano-Andreas, and M. Merchaan, "Molecular basis of DNA photodimerization: intrinsic production of cyclobutane cytosine dimers," Journal of the American Chemical Society, vol. 130, no. 32, pp. 10768-10779, 2008.
 J. J. Serrano-Perez, I. Gonzalez-Ramirez, P. B. Coto, M. Merchan, and L. Serrano-Andres, "Theoretical insight into the intrinsic ultrafast formation of cyclobutane pyrimidine dimers in UV-irradiated DNA: thymine versus cytosine," Journal of Physical Chemistry B, vol. 112, no. 45, pp. 14096-14098, 2008.
 Y. Dou, B. R. Torralva, and R. E. Allen, "Semiclassical electron-radiation-ion dynamics (SERID) and cis-trans photoisomerization of butadiene," Journal of Modern Optics, vol. 50, no. 15-17, pp. 2615-2643, 2003.
 Y. Dou, B. R. Torralva, and R. E. Allen, "Interplay of electronic and nuclear degrees of freedom in a femtosecond-scale photochemical reaction," Chemical Physics Letters, vol. 392, no. 4-6, pp. 352-357, 2004.
 M. Graf and P. Vogl, "Electromagnetic fields and dielectric response in empirical tight-binding theory," Physical Review B, vol. 51, no. 8, pp. 4940-4949, 1995.
 R. E. Allen, T. Dumitrica, and B. R. Torralva, "Electronic and structural response of materials to fast, intense laser pulses," in Ultrafast Physical Processes in Semiconductors, Tsen, K. T. Tsen, Ed., chapter 7, Academic Press, New York, NY, USA, 2001.
 M. Ben-Nun and T. J. Martinez, "AB initio quantum molecular dynamics," Advances in Chemical Physics, vol. 121, pp. 439-512, 2002.
 S. Yuan, W. Zhang, L. Liu, Y. Dou, W. Fang, and G. V. Lo, "Detailed mechanism for photoinduced cytosine dimerization: a semiclassical dynamics simulation," Journal of Physical Chemistry A, vol. 115, no. 46, pp. 13291-13297, 2011.
 A. Sharonov, T. Gustavsson, V. Carre, E. Renault, and D. Markovitsi, "Cytosine excited state dynamics studied by femtosecond fluorescence upconversion and transient absorption spectroscopy," Chemical Physics Letters, vol. 380, no. 1-2, pp. 173-180, 2003.
 B. Kohler, "Nonradiative decay mechanisms in DNA model systems," Journal of Physical Chemistry Letters, vol. 1, no. 13, pp. 2047-2053, 2010.
 E. R. Bittner, "Lattice theory of ultrafast excitonic and charge-transfer dynamics in DNA," Journal of Chemical Physics, vol. 125, no. 9, Article ID 094909, 12 pages, 2006.
 R. Improta, F. Santoro, V. Barone, and A. Lami, "Vibronic model for the quantum dynamical study of the competition between bright and charge-transfer excited states in single-strand polynucleotides: the adenine dimer case," The Journal of Physical Chemistry A, vol. 113, no. 52, pp. 15346-15354, 2009.
 W. Zhang, Semiclassical dynamics Simulation for deactivation of stacked bases of DNA [Ph.D. dissertation], Northwest University, Xi'an, China, 2011.
 H. Iikura, T. Tsuneda, T. Yanai, and K. Hirao, "A long-range correction scheme for generalized-gradient-approximation exchange functionals," Journal of Chemical Physics, vol. 115, no. 8, pp. 3540-3544, 2001.
 J. W. Song, T. Hirosawa, T. Tsuneda, and K. Hirao, "Long-range corrected density functional calculations of chemical reactions: redetermination of parameter," Journal of Chemical Physics, vol. 126, no. 15, Article ID 154105, 2007.
 O. A. Vydrov and G. E. Scuseria, "Intermittent search process and teleportation," Journal of Chemical Physics, vol. 126, no. 23, Article ID 234109, 2007.
 T. M. Henderson, B. G. Janesko, and G. E. Scuseria, "Generalized gradient approximation model exchange holes for range-separated hybrids," Journal of Chemical Physics, vol. 128, no. 19, Article ID 194105, 2008.
 F. Santoro, V. Barone, and R. Improta, "Can TD-DFT calculations accurately describe the excited states behavior of stacked nucleobases? The cytosine dimer as a test case," Journal of Computational Chemistry, vol. 29, no. 6, pp. 957-964, 2008.
 F. Santoro, V. Barone, and R. Improta, "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 tetramer," Journal of the American Chemical Society, vol. 131, no. 42, pp. 15232-15245, 2009.
 A. W. Lange and J. M. Herbert, "Both intra- And interstrand charge-transfer excited states in aqueous B-DNA are present at energies comparable to, or just above, the 1[pi][pi]* excitonic bright states," Journal of the American Chemical Society, vol. 131, no. 11, pp. 3913-3922, 2009.
Shuai Yuan, (1) Huiling Hong, (1) Gang Wang, (1) Wenying Zhang, (1) Yusheng Dou, (1,2) and Glenn V. Lo (2)
(1) College of Bio-information, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
(2) Department of Physical Sciences, Nicholls State University, P.O.Box 2022, Thibodaux, LA 70310, USA
Correspondence should be addressed to Yusheng Dou; email@example.com
Received 13 April 2014; Accepted 30 May 2014; Published 3 August 2014
Academic Editor: Fuli Li
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Yuan, Shuai; Hong, Huiling; Wang, Gang; Zhang, Wenying; Dou, Yusheng; Lo, Glenn V.|
|Publication:||International Journal of Photoenergy|
|Date:||Jan 1, 2014|
|Previous Article:||A simple and efficient MPPT method for low-power PV cells.|
|Next Article:||Optical, electrical and photocatalytic properties of the ternary semiconductors [Zn.sub.x][Cd.sub.1-x]S, [Cu.sub.x][Cd.sub.1-x]S, and...|