The decay from the dark npi* excited state in uracil: an integrated CASPT2/CASSCF and PCM/TD-DFT study in the gas phase and in water

By integrating the results of MS-CASPT2/CASSCF and TD-PBE0 calculations, we propose a mechanism for the decay of the excited dark state in pyrimidine, fully consistent with all the available experimental results. An effective conical intersection (CI-npi) exists between the spectroscopic pi/pi* excited state (Spi) and a dark n/pi* state (Sn), and a fraction of the population decays to the minimum of Sn (Sn-min). The conical intersection between Sn and the ground-state is not involved in the decay mechanism, because of its high energy gap with respect to Sn-min. On the other hand, especially in hydrogen bonding solvents, the energy gap between Sn-min and CI-npi is rather small. After thermalization in Sn-min, the system can thus recross CI-npi and then quickly proceed on the Spi barrierless path toward the conical intersection with the ground state.     

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.     

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.     

Solvent effect on the singlet excited-state lifetimes of nucleic acid bases: A computational study of 5-fluorouracil and uracil in acetonitrile and water

The first comprehensive quantum mechanical study of solvent effects on the behavior of the two lowest energy excited states of uracil derivatives is presented. The absorption and emission spectra of uracil and 5-fluorouracil in acetonitrile and aqueous solution have been computed at the time-dependent density-functional theory level, using the polarizable continuum model (PCM) to take into account bulk solvent effects. The computed spectra and the solvent shifts provided by our method are close to their experimental counterpart. The S0/S1 conical intersection, located in the presence of hydrogen-bonded solvent molecules by CASSCF (8/8) calculations, indicates that the mechanism of ground-state recovery, involving out-of-plane motion of the 5 substituent, does not depend on the nature of the solvent. Extensive explorations of the excited-state surfaces in the Franck-Condon (FC) region show that solvent can modulate the accessibility of an additional decay channel, involving a dark n/pi* excited state. This finding provides the first unifying explanation for the experimental trend of 5-fluorouracil excited-state lifetime in different solvents. The microscopic mechanisms underlying solvent effects on the excited-state behavior of nucleobases are discussed.     

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.     

Photoisomerization dynamics of azobenzene in solution with S1 excitation: a femtosecond fluorescence anisotropy study

Measurements of anisotropy of femtosecond fluorescence after direct excitation of the S1(n,pi*) state of azobenzene in hexane and ethylene glycol solutions have been carried out to address the controversy about inversion and rotation in the mechanism of photoisomerization. The observed anisotropies in hexane decay to a nonzero asymptotic level with a relaxation period the same as that for slow decay of the corresponding biexponential transient; this effect is attributed to involvement of the out-of-plane CNNC-torsional motion on approach to a twisted conical intersection along the "rotation channel" that depolarizes the original in-plane transition moment. In contrast, when the rotational channel becomes substantially hindered in ethylene glycol, the anisotropies show no discernible decay feature, but the corresponding transients show prominent decays attributed to involvement of in-plane symmetric motions; the latter approach a planar-sloped conical intersection along a "concerted inversion channel" for efficient internal conversion through vibronic coupling. The proposed mechanism is consistent with theoretical calculations and rationalizes both results on quantum yields and ultrafast observations.     

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.     

Quenching of tryptophan (1)(pi,pi*) fluorescence induced by intramolecular hydrogen abstraction via an aborted decarboxylation mechanism

CASSCF computations show that the hydrogen-transfer-induced fluorescence quenching of the (1)(pi,pi*) excited state of zwitterionic tryptophan occurs in three steps: (1) formation of an intramolecular excited-state complex, (2) hydrogen transfer from the amino acid side chain to the indole chromophore, and (3) radiationless decay through a conical intersection, where the reaction path bifurcates to a photodecarboxylation and a phototautomerization route. We present a general model for fluorescence quenching by hydrogen donors, where the radiationless decay occurs at a conical intersection (real state crossing). At the intersection, the reaction responsible for the quenching is aborted, because the reaction path bifurcates and can proceed forward to the products or backward to the reactants. The position of the intersection along the quenching coordinate depends on the nature of the states and, in turn, affects the formation of photoproducts during the quenching. For a (1)(n,pi*) model system reported earlier (Sinicropi, A.; Pogni, R.; Basosi, R.; Robb, M. A.; Gramlich, G.; Nau, W. M.; Olivucci, M. Angew. Chem., Int. Ed. 2001, 40, 4185-4189), the ground and the excited state of the chromophore are hydrogen acceptors, and the excited-state hydrogen transfer is nonadiabatic and leads directly to the intersection point. There, the hydrogen transfer is aborted, and the reaction can return to the reactant pair or proceed further to the hydrogen-transfer products. In the tryptophan case, the ground state is not a hydrogen acceptor, and the excited-state hydrogen transfer is an adiabatic, sequential proton and electron transfer. The decay to the ground state occurs along a second reaction coordinate associated with decarboxylation of the amino acid side chain and the corresponding aborted conical intersection. The results show that, for (1)(pi,pi*) states, the hydrogen transfer alone is not sufficient to induce the quenching, and explain why fluorescence quenching induced by hydrogen donors is less general for (1)(pi,pi*) than for (1)(n,pi*) states.     

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.     

Theoretical characterization of photoisomerization channels of dimethylpyridines on the singlet and triplet potential energy surfaces

Photoexcitations and photoisomerizations due to low-lying n pi* and pi pi* excited states of dimethylpyridines are investigated by density functional theory, CASSCF, CASPT2 and MRCI methodologies. Mechanistic details for the formation of Dewar dimethylpyridines and the interconversions of dimethylpyridines are rationalized through the characterization of minima and transition states on the singlet and triplet potential energy surfaces of relevant intermediates. Our present theoretical schemes suggest that Mobius dimethylpyridine intermediate 14 and azabenzvalene intermediate 10 can serve as possible precursors to Dewar dimethylpyridines and singlet phototransposition products, respectively. The calculations suggest that an S1(pi pi*)/S0 conical intersection in dimethylpyridines 2 is involved in the formation of 14. An azabenzvalene 10 might be formed through S2(pi pi*)/S1(n pi*) interaction followed by an S1/S0 decay in dimethylpyridine 6. Calculated barriers of isomerizations from 14 to Dewar dimethylpyridine 7 and from 10 to 4 are 8.4 and 28.5 kcal mol(-1) at the B3LYP/6-311 G** level, respectively. In the suggested triplet multistage transposition mechanism, an out-of-plane distorted geometry 19 due to vibrational relaxation of the T1(3B1) excited state of 3,5-dimethylpyridine 6 is a precursor of the interconversion of 6 to 2.4-dimethylpyridine 4. The formation of a triplet azaprefulvene 21 with a barrier of 20.7 kcal mol(-1) is a key step during the triplet migration process leading to another out-of-plane distorted structure 27. Subsequent rearomatization of 27 completes the interconversion of 6 with 4. Present calculations provide some insight into the photochemistry of dimethylpyridines at 254 nm.     

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.     

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.     

Ultrafast decay of electronically excited singlet cytosine via a pi,pi* to n(O),pi* state switch

Singlet fluorescence lifetimes of adenosine, cytidine, guanosine, and thymidine, determined by femtosecond pump-probe spectroscopy (Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2000, 122, 9348. Pecourt, J.-M. L.; Peon, J.; Kohler, B. J. Am. Chem. Soc. 2001, 123, 10370), show that the excited states produced by 263 nm light in these nucleosides decay in the subpicosecond range (290-720 fs). Ultrafast radiationless decay to the ground state greatly reduces the probability of photochemical damage. In this work we present a theoretical study of isolated cytosine, the chromophore of cytidine. The experimental lifetime of 720 fs indicates that there must be an ultrafast decay channel for this species. We have documented the possible decay channels and approximate energetics, using a valence-bond derived analysis to rationalize the structural details of the paths. The mechanism favored by our calculations and the experimental data involves (1) a two-mode decay coordinate composed of initial bond length inversion followed by internal vibrational energy redistribution (IVR) to populate a carbon pyramidalization mode, (2) a state switch between the pi,pi* and nO,pi* (excitation from oxygen lone pair) excited states, and (3) decay to the ground state through a conical intersention. A second decay path through the nN,pi* state (excitation from the nitrogen lone pair), with a higher barrier, involves out-of-plane bending of the amino substituent.     

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.     

Conical for stepwise, glancing for concerted: the role of the excited-state topology in the three-body dissociation of sym-triazine

The highly debated three-body dissociation of sym-triazine to three HCN products has been investigated by translational spectroscopy and high-level ab initio calculations. Dissociation was induced by charge exchange between the sym-triazine radical cation and cesium. Calculated state energies and electronic couplings suggest that sym-triazine is produced in the 3s Rydberg and pi* <-- n manifolds. Analysis of the topology of these manifolds along with momentum correlation in the dissociation products suggest that the 3s Rydberg manifold characterized by a conical intersection of two potential energy surfaces leads to stepwise dissociation, while the pi* <-- n manifold consisting of a four-fold glancing intersection leads to a symmetric concerted reaction.     

Water effect on the excited-state decay paths of singlet excited cytosine

Two low-energy deactivation paths for singlet excited cytosine, one through a S₁/S₀ conical intersection of the ethylene type, and one through a conical intersection that involves the (n_N, π^*) state, are calculated in the presence of water. Water is included explicitly for several cytosine monohydrates, and as a bulk solvent, and the calculations are carried out at the complete active space self-consistent field (CASSCF) and complete active space second order perturbation (CASPT2) levels of theory. The effect of water on the lowest-energy path through the ethylenic conical intersection is a lowering of the energy barrier. This is explained by stabilization of the excited state, which has zwitterionic character in the vicinity of the conical intersection due to its similarity with the conical intersection of ethylene. In contrast to this, the path that involves the (n_N, π^*) state is destabilized by hydrogen bonding, although the bulk solvent effect reduces the destabilization. Overall, this path should remain energetically accessible.     

Fluorescence on-off switching mechanism of benzofurazans

Many fluorescent reagents with a benzofurazan (2,1,3-benzoxadiazole) skeleton have been developed and widely used in bio-analyses. In this study, we try to elucidate the fluorescence on-off switching mechanism of three fluorogenic reagents and their derivatives. Ten 4,7-disubstituted benzofurazans were used for this purpose and the measurements of their fluorescence, phosphorescence, photolysis, and time-resolved thermal lensing signal in acetonitrile were obtained in order to understand the relaxation processes of these compounds. These results indicate that the competition of fluorescence with a fast intersystem crossing or fast photoreaction plays a key role in the fluorescence on-off switching. Semi-empirical molecular orbital calculations show that the existence of the triplet n pi* state is responsible for the fast intersystem crossing while the proximity of the reactive second single pi pi* state to the first singlet pi pi* state contributes to the fast photoreaction in the excited states.     

Computational studies of the photophysics of hydrogen-bonded molecular systems

The role of electron- and proton-transfer processes in the photophysics of hydrogen-bonded molecular systems has been investigated with ab initio electronic-structure calculations. Adopting indole, pyridine, and ammonia as molecular building blocks, we discuss generic mechanisms of the photophysics of isolated aromatic chromophores (indole), complexes of pi systems with solvent molecules (indole-ammonia, pyridine-ammonia), hydrogen-bonded aromatic pairs (indole-pyridine), and intramolecularly hydrogen-bonded pi systems (7-(2'-pyridyl)indole). The reaction mechanisms are discussed in terms of excited-state minimum-energy paths, conical intersections, and the properties of frontier orbitals. A common feature of the photochemistry of the various systems is the electron-driven proton-transfer (EDPT) mechanism: highly polar charge-transfer states of 1pipi*, 1npi*, or 1pisigma* character drive the proton transfer, which leads, in most cases, to a conical intersection of the S1 and S0 surfaces and thus ultrafast internal conversion. In intramolecularly hydrogen-bonded aromatic systems, out-of-plane torsion is additionally needed for barrierless access to the S1-S0 conical intersection. The EDPT process plays an essential role in diverse photophysical phenomena, such as fluorescence quenching in protic solvents, the function of organic photostabilizers, and the photostability of biological molecules.     

Direct versus indirect H atom elimination from photoexcited phenol molecules

The active role of the optically dark pi sigma* state, following UV absorption, has been implicated in the photochemistry of a number of biomolecules. This work focuses on the role of the pi sigma* state in the photochemistry of phenol upon excitation at 200 nm. By probing the neutral hydrogen following UV excitation, we show that hydrogen elimination along the dissociative pi sigma* potential energy surface occurs within 103 +/- 30 fs, indicating efficient coupling at the S1/S2 and S0/S2 conical intersections, with no identifiable role of statistical unimolecular decay of vibronically excited (S0) phenol in the timeframe of our measurements.     

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.