Ab initio studies on the photophysics of guanine tautomers: out-of-plane deformation and NH dissociation pathways to conical intersections

The radiationless decay mechanisms of the S1 excited states of the 7H-keto-amino, 7H-enol-amino, and 7H-keto-imino tautomers of guanine have been investigated with the CASPT2//CASSCF method. Out-of-plane deformation of the six-membered ring or the imino group as well as dissociation of NH bonds have been considered as photochemical pathways leading to conical intersections with the electronic ground state. It has been found that all three tautomers can reach S0-S1 conical intersections by out-of-plane deformation. However, only in the 7H-keto-amino tautomer the reaction path leading to the conical intersection is barrierless. This tautomer also has the lowest energy barrier for hydrogen detachment via the (1)pi sigma* state, whose potential energy surface intersects that of the (1)pi pi* state as well as that of the ground state. The other tautomers of guanine exhibit substantial energy barriers on their S1 potential energy surfaces with respect to both reaction mechanisms. These findings suggest that the 7H-keto-amino tautomer exhibits the shortest excited-state lifetime of the three tautomers due to particularly fast nonradiative deactivation processes through S0-S1 conical intersections. The computational results explain the remarkable observation that the energetically most stable 7H-keto-amino tautomer is missing in the resonant two-photon ionization spectrum of guanine in a supersonic jet. The results also explain that the energetically less stable 7H-enol-amino and 7H-keto-imino tautomers have longer excited-state lifetimes and are thus detectable by resonant two-photon ionization.     

Theoretical study on the singlet excited state of pterin and its deactivation pathway

The excited-state properties and related photophysical processes of the acidic and basic forms of pterin have been investigated by the density functional theory and ab initio methodologies. The solvent effects on the low-lying states have been estimated by the polarized continuum model and combined QM/MM calculations. Calculations reveal that the observed two strong absorptions arise from the strong pi --> pi* transitions to 1(pipi*L(a)) and 1(pipi*L(b)) in the acidic and basic forms of pterin. The first 1(pipi*L(a)) excited state is exclusively responsible for the experimental emission band. The vertical 1(n(N)pi*) state with a small oscillator strength, slightly higher in energy than the 1(pipi*L(a)) state, is less accessible by the direct electronic transition. The 1(n(N)pi*) state may be involved in the photophysical process of the excited pterin via the 1(pipi*L(a)/n(N)pi*) conical intersection. The radiationless decay of the excited PT to the ground state experiences a barrier of 13.8 kcal/mol for the acidic form to reach the (S(1)/S(0)) conical intersection. Such internal conversion can be enhanced with the increase in excitation energy, which will reduce the fluorescence intensity as observed experimentally.     

The fluorescence mechanism of 5-methyl-2-pyrimidinone: an ab initio study of a fluorescent pyrimidine analog

The photophysically important potential energy surfaces of the fluorescent pyrimidine analog 5-methyl-2-pyrimidinone have been explored using multireference configuration-interaction ab initio methods at three levels of dynamical correlation, all of which support a fluorescence mechanism. At vertical excitation S1 (dark, n(N)pi*) and S2 (bright, pipi*) are almost degenerate at 4.4 eV, with S3 (dark, n(O)pi*) at 5.1 eV. The excited system can follow the S1-S2 seam of conical intersections, accessible from the Franck-Condon region, to its minimum and then evolve from this conical intersection on the S1 (pipi*) surface to a global minimum. At lower levels of correlation, the S1 surface shows two minima separated by a barrier of up to 0.18 eV. The secondary minimum found at the lower levels of correlation becomes the global minimum with higher correlation. The S1 population at this minimum can be trapped from accessing the lowest energy S0-S1 (pipi*/gs) conical intersection by an energy gap at least 0.3-0.4 eV higher than the S1 minimum. The calculated emission energy from this minimum is 2.80 eV. Gradient pathways connecting important S1 geometries are presented, as well as other excited state conical intersections.     

Radiationless decay mechanism of cytosine: an ab initio study with comparisons to the fluorescent analogue 5-methyl-2-pyrimidinone

The ultrafast radiationless decay mechanism of photoexcited cytosine has been theoretically supported by exploring the important potential energy surfaces using multireference configuration-interaction ab initio methods for the gas-phase keto-tautomer free base. At vertical excitation, the bright state is S1 (pipi*) at 5.14 eV, with S2 (nNpi*) and S3 (nOpi*) being dark states at 5.29 and 5.93 eV, respectively. Minimum energy paths connect the Franck-Condon region to a shallow minimum on the pipi* surface at 4.31 eV. Two different energetically accessible conical intersections with the ground state surface are shown to be connected to this minimum. One pathway involves N3 distorting out of plane in a sofa conformation, and the other pathway involves a dihedral twist about the C5-C6 bond. Each of these pathways from the minimum contains a low barrier of 0.14 eV, easily accessed by low vibronic levels. The path involving the N3 sofa distortion leads to a conical intersection with the ground state at 4.27 eV. The other pathway leads to an intersection with the ground state at 3.98 eV, lower than the minimum by about 0.3 eV. Comparisons with our previously reported study of the fluorescent cytosine analogue 5-methyl-2-pyrimidinone (5M2P) reveal remarkably similar conformational distortions throughout the decay pathways of both bases. The different photophysical behavior between the two molecules is attributed to energetic differences. Vertical excitation in cytosine occurs at a much higher energy initially, creating more vibrational energy than 5M2P in the Franck-Condon region, and the minimum S1 energy for 5M2P is too low to access an intersection with the ground state, causing population trapping and fluorescence. Calculations of vertical excitation energies of 5-amino-2-pyrimidinone and 2-pyrimidinone reveal that the higher excitation energy of cytosine is likely due to the presence of the amino group at the 4-position.     

Ab initio study on deactivation pathways of excited 9H-guanine

The complete active space with second-order perturbation theory/complete active space self-consistent-field method was used to explore the nonradiative decay mechanism for excited 9H-guanine. On the 1pipi* (1L(a)) surface we determined a conical intersection (CI), labeled (S0pipi*)(CI), between the 1pipi* (1L(a)) excited state and the ground state, and a minimum, labeled (pipi*)min. For the 1pipi* (1L(a)) state, its probable deactivation path is to undergo a spontaneous relaxation to (pipi*)min first and then decay to the ground state through (S0pipi*)(CI), during which a small activation energy is required. On the 1n(N)pi* surface a CI between the 1n(N)pi* and 1pipi* (1L(a)) states was located, which suggests that the 1n(N)pi* excited state could transform to the 1pipi* (1L(a)) excited state first and then follow the deactivation path of the 1pipi* (1L(a)) state. This CI was also possibly involved in the nonradiative decay path of the second lowest 1pipi* (1L(b)) state. On the 1n(O)pi* surface a minimum was determined. The deactivation of the 1n(O)pi* state to the ground state was estimated to be energetically unfavorable. On the 1pisigma* surface, the dissociation of the N-H bond of the six-membered ring is difficult to occur due to a significant barrier.     

A three-state model for the photophysics of adenine

An ab initio theoretical study at the CASPT2 level is reported on minimum energy reaction paths, state minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of two tautomers of adenine: 9H- and 7H-adenine. The obtained results led to a complete interpretation of the photophysics of adenine and derivatives, both under jet-cooled conditions and in solution, within a three-state model. The ultrafast subpicosecond fluorescence decay measured in adenine is attributed to the low-lying conical intersection (gs/pipi* La)(CI), reached from the initially populated 1(pipi* La) state along a path which is found to be barrierless only in 9H-adenine, while for the 7H tautomer the presence of an intermediate plateau corresponding to an NH2-twisted conformation may explain the absence of ultrafast decay in 7-substituted compounds. A secondary picosecond decay is assigned to a path involving switches towards two other states, 1(pipi* Lb) and 1(npi*), ultimately leading to another conical intersection with the ground state, (gs/npi*), with a perpendicular disposition of the amino group. The topology of the hypersurfaces and the state properties explain the absence of secondary decay in 9-substituted adenines in water in terms of the higher position of the 1(npi*) state and also that the 1(pipi* Lb) state of 7H-adenine is responsible for the observed fluorescence in water. A detailed discussion comparing recent experimental and theoretical findings is given. As for other nucleobases, the predominant role of a pipi*-type state in the ultrafast deactivation of adenine is confirmed.     

Singlet excited-state dynamics of 5-fluorocytosine and Cytosine: an experimental and computational study

The photophysics of singlet excited 5-fluorocytosine (5FC) was studied in steady-state and time-resolved experiments and theoretically by quantum chemical calculations. Femtosecond transient absorption measurements show that replacement of the C5 hydrogen of cytosine by fluorine increases the excited-state lifetime by 2 orders of magnitude from 720 fs to 73 +/- 4 ps. Experimental evidence indicates that emission in both compounds originates from a single tautomeric form. The lifetime of 5FC is the same within experimental uncertainty in the solvents ethanol and dimethyl sulfoxide. The insensitivity of the S(1) lifetime to the protic nature of the solvent suggests that proton transfer is not the principal quenching mechanism for the excited state. Excited-state calculations were carried out for the amino-keto tautomer of 5FC, the dominant species in polar environments, in order to understand its longer excited-state lifetime. CASSCF and CAS-PT2 calculations of the excited states show that the minimum energy path connecting the minimum of the (1)pi,pi state with the conical intersection responsible for internal conversion has essentially the same energetics for cytosine and 5FC, suggesting that both bases decay nonradiatively by the same mechanism. The dramatic difference in lifetimes may be due to subtle changes along the decay coordinate. A possible reason may be differences in the intramolecular vibrational redistribution rate from the Franck-Condon active, in-plane modes to the out-of-plane modes that must be activated to reach the conical intersection region.     

Energetics of cytosine singlet excited-state decay paths--a difficult case for CASSCF and CASPT2

Three deactivation paths for singlet excited cytosine are calculated at the CASPT2//CASSCF (complete active space second-order perturbation//complete active space self-consistent field) level of theory, using extended active spaces that allow for a reliable characterization of the paths and their energies. The lowest energy path, with a barrier of approximately 0.1 eV, corresponds to torsion of the C5-C6 bond, and the decay takes place at a conical intersection analogous to the one found for ethylene and its derivatives. There is a further path with a low energy barrier of approximately 0.2 eV associated with the (n(N),pi*) state which could also be populated with a low energy excitation. The path associated with a conical intersection between the ground and (n(O),pi*) states is significantly higher in energy (> 1 eV). The presence of minima on the potential energy surface for the (n,pi*) states that could contribute to the biexponential decay found in the gas phase was investigated, but could not be established unequivocally.     

Theoretical study of the ground and excited states of 7-methyl guanine and 9-methyl guanine: comparison with experiment

The keto-enol tautomerization of 7-methyl-guanine and 9-methyl-guanine in the excited state was investigated using the time-dependent DFT (TDDFT) method. For both species, the potential energy surfaces of the ground state and two lowest singlet excited states (due to pi-->pi* and n-->pi* transitions) have been investigated and their features discussed in terms of consequences on the excited state dynamics. The findings suggest that, for both species, the state due to the n-->pi* transition, suspected to be an intermediate in the excited state deactivation, exhibits two minima with the second minimum characterized by an elongated N1-H distance. This structure, intermediate between enol and keto tautomers, might play a role in the excited state relaxation. The existence of this second well, however, is observed in both 7- and 9-methyl-guanine, which suggests that it cannot account alone for the different photophysical behavior of these species.     

Excited states of thiophene: ring opening as deactivation mechanism

(Time-dependent) Kohn-Sham density functional theory and a combined density functional/multi-reference configuration interaction method (DFT/MRCI) were employed to explore the ground and low-lying electronically excited states of thiophene. Spin-orbit coupling was taken into account using an efficient, nonempirical mean-field Hamiltonian. Phosphorescence lifetimes were calculated by means of spock.ci, a selecting direct multi-reference spin-orbit configuration interaction program. Throughout this paper, we use the following nomenclature: S1, S2,..., T1, T2,..., denominate electronic structures in their energetic order at the ground state minimum geometry, whereas S1, S2,..., T1, T2,..., refers to the actual order of electronic states at a given nuclear geometry. Multiple minima were found on the first excited singlet (S1) potential energy hypersurface with electronic structures S1 (piHOMO-1-->pi+piHOMO-->pi), S2 (piHOMO-->pi), and S3 (piHOMO-->sigma*) corresponding to the 2 1A1 (S1), 1 1B2 (S2), and 1 1B1 (S3) states in the vertical absorption spectrum, respectively. The S1 and S2 minimum geometries show out-of-plane deformations of the ring. The S3 electronic structure yields the global minimum on the S1 surface with an adiabatic excitation energy of merely 3.81 eV. It exhibits an asymmetric planar nuclear arrangement with one significantly elongated C-S bond. A constrained minimum energy path calculation connecting the S1 and S3 minima suggests that even low-lying vibrational levels of the S1 potential well can access the global minimum of the S1 surface. Nonradiative decay of the electronically excited singlet population to the electronic ground state via a close-by conical intersection will be fast. According to our work, this ring opening mechanism is most likely responsible for the lack of fluorescence in thiophene and the ultrafast decay of the S1 vibrational levels, as observed in time-resolved pump-probe femtosecond multiphoton ionization experiments. An alternative relaxation pathway leads from the S1 minimum via vibronic coupling to the S2 potential well followed by fast inter-system crossing to the T2 state. For an estimate of individual rate constants a quantum dynamical treatment will be required. The global minimum of the T1 surface has a chair-like nuclear conformation and corresponds to the T1 (1 3B2, piHOMO-->pi) electronic structure. Phosphorescence is weak here with a calculated radiative lifetime of 0.59 s. For the second potential well on the T1 surface with T3 (1 3B1, piHOMO-->sigma*) electronic structure, nonradiative processes are predicted to dominate the triplet decay.     

Enol-keto tautomerism of aromatic photochromic Schiff base N,N'-bis(salicylidene)-p-phenylenediamine: ground state equilibrium and excited state deactivation studied by solvatochromic measurements on ultrafast time scale

A photochromic symmetric Schiff base, N,N'-bis(salicylidene)-p-phenylenediamine, is proposed as a probe for the study of solvent dependent enol-keto tautomerism in the ground and excited states. The ground state equilibrium between the enol-keto tautomers is found to depend mainly not on polarity but on the proton donating ability of the solvent. Upon selective excitation of each of these tautomers, the same excited state of a keto tautomer is created: in enol, after the ultrafast excited state intramolecular proton transfer (ESIPT), reaction, and in keto tautomer, directly. Then some part (<30%) of excited molecules are transferred to the photochromic form in its ground state. The evidence of another ultrafast deactivation channel in the excited enol tautomer competing with ESIPT has been found. The solvent does not influence the ESIPT dynamics nor the efficiency of the creation of the photochrome.     

A new pathway for the rapid decay of electronically excited adenine

Combined density functional and multireference configuration interaction methods have been used to calculate the electronic spectrum of 9H-adenine, the most stable tautomer of 6-aminopurine. In addition, constrained minimum energy paths on excited potential energy hypersurfaces have been determined along several relaxation coordinates. The minimum of the first (1)[n-->pi*] state has been located at an energy of 4.54 eV for a nuclear arrangement in which the amino group is pyramidal whereas the ring system remains planar. Close by, another minimum on the S(1) potential energy hypersurface has been detected in which the C(2) center is deflected out of the molecular plane and the electronic character of S(1) corresponds to a nearly equal mixture of (1)[pi-->pi*] and (1)[n-->pi*] configurations. The adiabatic excitation energy of this minimum amounts to 4.47 eV. Vertical and adiabatic excitation energies of the lowest n-->pi* and pi-->pi* transitions as well as transition moments and their directions are in very good agreement with experimental data and lend confidence to the present quantum chemical treatment. On the S(1) potential energy hypersurface, an energetically favorable path from the singlet n-->pi* minimum toward a conical intersection with the electronic ground state has been identified. Close to the conical intersection, the six-membered ring of adenine is strongly puckered and the electronic structure of the S(1) state corresponds to a pi-->pi* excitation. The energetic accessibility of this relaxation path at about 0.1 eV above the singlet n-->pi* minimum is presumably responsible for the ultrafast decay of 9H-adenine after photoexcitation and explains why sharp vibronic peaks can only be observed in a rather narrow wavelength range above the origin. The detected mechanism should be equally applicable to adenosine and 9-methyladenine because it involves primarily geometry changes in the six-membered ring whereas the nuclear arrangement of the five-membered ring (including the N(9) center) is largely preserved.     

The electronic spectrum of protonated adenine: theory and experiment

In this work we present the results of a combined experimental and theoretical study concerned with the question how a proton changes the electronic spectrum and dynamics of adenine. In the experimental part, isolated adenine ions have been formed by electro-spray ionisation, stored, mass-selected and cooled in a Paul trap and dissociated by resonant photoexcitation with ns UV laser pulses. The S(0)-S1 spectrum of protonated adenine recorded by fragment ion detection lies in a similar energy range as the first pipi* transition of neutral 9H-adenine. It shows a flat onset with a broad substructure, indicating a large S(0)-S1 geometry shift and an ultra-short lifetime. In the theoretical part, relative energies of the ground and the excited states of the most important tautomers have been calculated by means of a combined density functional theory and multi-reference configuration interaction approach. Protonation at the nitrogen in position 1 of the neutral 9H-adenine tautomer yields the most stable protonated adenine species, 1H-9H-A+. The 3H-7H-A+ and the 3H-9H-A+ tautomers, formed by protonation of 7H- and 9H-adenine in 3-position, are higher in energy by 162 cm(-1) and 688 cm(-1), respectively. Other tautomers lie at considerably higher energies. Calculated vertical absorption spectra are reported for all investigated tautomers whereas geometry optimisations of excited states have been carried out only for the most interesting ones. The S1 state energies and geometries are found to depend on the protonation site. The theoretical data match best with the experimental onset of the spectrum for the 1H-9H-A+ tautomer although we cannot definitely exclude contributions to the experimental spectrum from the 3H-7H-A+ tautomer at higher energies. The vertical S(0)--> S1 excitation energy is similar to the one in neutral 9H-adenine. As for the neutral adenine, we find a conical intersection of the S1 of protonated adenine with the ground state in an out-of-plane coordinate but at lower energies and accessible without barrier.     

Conical intersections in thymine

The mechanisms which are responsible for the radiationless deactivation of the npi* and pipi* excited singlet states of thymine have been investigated with multireference ab initio methods (the complete-active-space self-consistent-field (CASSCF) method and second-order perturbation theory with respect to the CASSCF reference (CASPT2)) as well as with the CC2 (approximated singles and doubles coupled-cluster) method. The vertical excitation energies, the equilibrium geometries of the 1npi*and 1pipi* states, as well as their adiabatic excitation energies have been determined. Three conical intersections of the S1 and S0 energy surfaces have been located. The energy profiles of the excited states and the ground state have been calculated with the CASSCF method along straight-line reaction paths leading from the ground-state equilibrium geometry to the conical intersections. All three conical intersections are characterized by strongly out-of-plane distorted geometries. The lowest-energy conical intersection (CI1) arises from a crossing of the lowest 1pipi* state with the electronic ground state. It is found to be accessible in a barrierless manner from the minimum of the 1pipi* state, providing a direct and fast pathway for the quenching of the population of the lowest optically allowed excited states of thymine. This result explains the complete diffuseness of the absorption spectrum of thymine in supersonic jets. The lowest vibronic levels of the optically nearly dark 1npi* state are predicted to lie below CI1, explaining the experimental observation of a long-lived population of dark excited states in gas-phase thymine.     

Modeling of HeN+ clusters. II. Calculation of He3+ vibrational spectrum

We have computed the vibrational spectrum of the helium ionized trimer He(3)(+) using three different potential energy surfaces [D. T. Chang and G. L. Gellene, J. Chem. Phys. 119, 4694 (2003); E. Scifoni et al., ibid. 125, 164304 (2006); I. Paidarova et al., Chem. Phys. 342, 64 (2007)]. Differences in the details of these potential energy surfaces induce discrepancies between bound state energies of the order of 0.01 eV. The effects of the geometric phase induced by the conical intersection between the ground electronic potential energy surface and the first excited one are studied by computing vibrational spectra with and without this phase. The six lowest vibrational bound states are negligibly affected by the geometric phase. Indeed, they correspond to wavefunctions localized in the vicinity of the linear symmetric configurations and can be assigned well defined vibrational quantum numbers. On the other hand, higher excited states are delocalized, cannot be assigned definite vibrational quantum numbers, and the geometric phase shifts their energies by approximately 0.005 eV.     

Disentangling intrinsic ultrafast excited-state dynamics of cytosine tautomers

Gas-phase ultrafast excited-state dynamics of cytosine, 1-methylcytosine, and 5-fluorocytosine were investigated in molecular beams using femtosecond pump-probe photoionization spectroscopy to identify the intrinsic dynamics of the major cytosine tautomers. The results indicate that, upon photoexcitation in the first absorption band, the cytosine enol tautomer exhibits a significantly longer excited-state lifetime than its keto and imino counterparts. The initially excited states of the cytosine keto and imino tautomers decay with sub-picosecond dynamics for excitation wavelengths shorter than 300 nm, whereas that of the cytosine enol tautomer decays with time constants ranging from 3 to 45 ps for excitation between 260 and 285 nm.     

Time-dependent quantum wave-packet description of the 1pi sigma* photochemistry of phenol

The photoinduced hydrogen elimination reaction in phenol via the conical intersections of the dissociative 1pi sigma* state with the 1pi pi* state and the electronic ground state has been investigated by time-dependent quantum wave-packet calculations. A model including three intersecting electronic potential-energy surfaces (S0, 1pi sigma*, and 1pi pi*) and two nuclear degrees of freedom (OH stretching and OH torsion) has been constructed on the basis of accurate ab initio multireference electronic-structure data. The electronic population transfer processes at the conical intersections, the branching ratio between the two dissociation channels, and their dependence on the initial vibrational levels have been investigated by photoexciting phenol from different vibrational levels of its ground electronic state. The nonadiabatic transitions between the excited states and the ground state occur on a time scale of a few tens of femtoseconds if the 1pi pi*-1pi sigma* conical intersection is directly accessible, which requires the excitation of at least one quantum of the OH stretching mode in the 1pi pi* state. It is shown that the node structure, which is imposed on the nuclear wave packet by the initial preparation as well as by the transition through the first conical intersection (1pi pi*-1pi sigma*), has a profound effect on the nonadiabatic dynamics at the second conical intersection (1pi sigma*-S0). These findings suggest that laser control of the photodissociation of phenol via IR mode-specific excitation of vibrational levels in the electronic ground state should be possible.     

A three-state model for the photophysics of guanine

The nonadiabatic photochemistry of the guanine molecule (2-amino-6-oxopurine) and some of its tautomers has been studied by means of the high-level theoretical ab initio quantum chemistry methods CASSCF and CASPT2. Accurate computations, based by the first time on minimum energy reaction paths, states minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of the molecules lead to interpret the photochemistry of guanine and derivatives within a three-state model. As in the other purine DNA nucleobase, adenine, the ultrafast subpicosecond fluorescence decay measured in guanine is attributed to the barrierless character of the path leading from the initially populated 1(pi pi* L(a)) spectroscopic state of the molecule toward the low-lying methanamine-like conical intersection (gs/pi pi* L(a))CI. On the contrary, other tautomers are shown to have a reaction energy barrier along the main relaxation profile. A second, slower decay is attributed to a path involving switches toward two other states, 1(pi pi* L(b)) and, in particular, 1(n(O) pi*), ultimately leading to conical intersections with the ground state. A common framework for the ultrafast relaxation of the natural nucleobases is obtained in which the predominant role of a pi pi*-type state is confirmed.     

Nonradiative decay mechanisms of the biologically relevant tautomer of guanine

We present the excited-state potential energy profiles of the biologically relevant 9H-keto-amino tautomer of guanine with respect to the radiationless decay via the out-of-plane deformation of the six-membered ring as well as the dissociation of NH bonds. The CASPT2//CASSCF method is employed for the reaction-path calculations. The reaction path for the out-of-plane deformation in the (1)pi pi* state leads in a barrierless way to a conical intersection with the electronic ground state. For the NH dissociation via the (1)pi sigma* state, the 9H-keto-amino tautomer is shown to have lower energy barriers than the 7H tautomers which we have studied recently. These two radiationless decay mechanisms explain the unexpected missing of the biologically relevant form in the resonant two-photon ionization spectrum of guanine in a supersonic jet. It is suggested that these ultrafast deactivation processes provide the biologically relevant tautomer of guanine with a high degree of photostability.