Origin, nature, and fate of the fluorescent state of the green fluorescent protein chromophore at the CASPT2//CASSCF resolution

Ab initio CASPT2//CASSCF relaxation path computations are employed to determine the intrinsic (e.g., in vacuo) mechanism underlying the rise and decay of the luminescence of the anionic form of the green fluorescent protein (GFP) fluorophore. Production and decay of the fluorescent state occur via a two-mode reaction coordinate. Relaxation along the first (totally symmetric) mode leads to production of the fluorescent state that corresponds to a planar species. The second (out-of-plane) mode controls the fluorescent state decay and mainly corresponds to a barrierless twisting of the fluorophore phenyl moiety. While a "space-saving" hula-twist conical intersection decay channel is found to lie only 5 kcal mol(-1) above the fluorescent state, the direct involvement of a hula-twist deformation in the decay is not supported by our data. The above results indicate that the ultrafast fluorescence decay observed for the GFP chromophore in solution is likely to have an intrinsic origin. The possible effects of the GFP protein cavity on the fluorescence lifetime of the investigated chromophore model are discussed.     

Ultrafast nonradiative decay of electronically excited States of malachite green: ab initio calculations

We have investigated the nonradiative deactivation process of malachite green in the singlet excited states, S(1) and S(2), by high-level ab initio quantum chemical calculations using the CASPT2//CASCF approach. The deactivation pathways connecting the Franck-Condon region and conical intersection regions are identified. The initial population in the S(1) state is on a flat surface and the relaxation involves a rotation of phenyl rings, which leads the molecule to reach the conical intersection between the S(1) and S(0) states, where it efficiently decays back to the ground state. There exists a small barrier connecting the Franck-Condon and conical intersection regions on the S(1) potential energy surface. The decay mechanism from the S(2) state also involves the twisting motion of phenyl rings. In contrast to the excitation to the S(1) state, the initial population is on a downhill ramp potential and the barrierless relaxation through the rotation of substituted phenyl rings is expected. During the course of relaxation, the molecule switches to the S(1) state at the conical intersection between S(2) and S(1), and then it decays back to the ground state through the intersection between S(1) and S(0). In relaxation from both S(1) and S(2), large distortion of phenyl rings is required for the ultrafast nonradiative decay to the ground state.     

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.     

Theoretical study of the isomerization mechanism of azobenzene and disubstituted azobenzene derivatives

A series of azobenzenes was studied using ab initio methods to determine the substituent effects on the isomerization pathways. Energy barriers were determined from three-dimensional potential energy surfaces of the ground and electronically excited states. In the ground state (S(0)), the inversion pathway was found to be preferred. Our results show that electron donating substituents increase the isomerization barrier along the inversion pathway, whereas electron withdrawing substituents decrease it. The inversion pathway of the first excited state (S(1)) showed trans --> cis barriers with no curve crossing between S(0) and S(1). In contrast, a conical intersection was found between the ground and first excited states along the rotation pathway for each of the azobenzenes studied. No barriers were found in this pathway, and we therefore postulate that after n --> pi (S(1) <-- S(0)) excitation, the rotation mechanism dominates. Upon pi --> pi (S(2) <-- S(0)) excitation, there may be sufficient energy to open an additional pathway (concerted-inversion) as proposed by Diau. Our potential energy surface explains the experimentally observed difference in trans-to-cis quantum yields between S(1) and S(2) excitations. The concerted inversion channel is not available to the remaining azobenzenes, and so they must employ the rotation pathway for both n --> pi and pi --> pi excitations.     

Photoisomerization of disperse red 1 studied with transient absorption spectroscopy and quantum chemical calculations

The photoisomerization of the push-pull substituted azo dye Disperse Red 1 is studied using femtosecond time-resolved absorption spectroscopy and other spectroscopic and computational techniques. In comparison with azobenzene, the pipi* state is more stabilized by the effects of push-pull substitution than the npi* state, but the latter is probably still the lowest in energy. This conclusion is based on the kinetics, anisotropy of the excited state absorption spectrum, the spectra of the ground states, and quantum chemical calculations. The S(1)(npi*) state is formed from the initially excited pipi* state in <0.2 ps, and decays to the ground state with time constants of 0.9 ps in toluene, 0.5 ps in acetonitrile, and 1.4 ps in ethylene glycol. Thermal isomerization transforms the Z isomer produced to the more stable E isomer with time constants of 29 s (toluene), 28 ms (acetonitrile), and 2.7 ms (ethylene glycol). The pathway of photoisomerization is likely to be rotation about the N=N bond. Quantum chemical calculations indicate that along the inversion pathway ground and excited state energy surfaces remain well separated, whereas rotation leads to a region where conical intersections can occur. For the ground-state Z to E isomerization, conclusive evidence is lacking, but inversion is more probably the favored pathway in the push-pull substituted systems than in the parent azobenzene.     

How Does the Trans–Cis Photoisomerization of Azobenzene Take Place in Organic Solvents?

The trans–cis photoisomerization of azobenzene‐containing materials is key to a number of photomechanical applications, but the actual conversion mechanism in condensed phases is still largely unknown. Herein, we study the ${{\rm{n}}{\rm{,{\rm \pi} ^\ast }}}$ isomerization in a vacuum and in various solvents via a modified molecular dynamics simulation adopting an ab initio torsion–inversion force field in the ground and excited states, while allowing for electronic transitions and a stochastic decay to the fundamental state. We determine the trans–cis photoisomerization quantum yield and decay times in various solvents (n‐hexane, anisole, toluene, ethanol, and ethylene glycol), and obtain results comparable with experimental ones where available. A profound difference between the isomerization mechanism in vacuum and in solution is found, with the often neglected mixed torsional–inversion pathway being the most important in solvents.     

Intramolecular charge transfer in 4-aminobenzonitriles does not necessarily need the twist

In electron donor/acceptor species such as 4-(dimethylamino)benzonitrile (DMABN), the excitation to the S(2) state is followed by internal conversion to the locally excited (LE) state. Dual fluorescence then becomes possible from both the LE and the twisted intramolecular charge-transfer (TICT) states. A detailed mechanism for the ICT of DMABN and 4-aminobenzonitrile (ABN) is presented in this work. The two emitting S(1) species are adiabatically linked along the amino torsion reaction coordinate. However, the S(2)/S(1) CT-LE radiationless decay occurs via an extended conical intersection "seam" that runs almost parallel to this torsional coordinate. At the lowest energy point on this conical intersection seam, the amino group is untwisted; however, the seam is accessible for a large range of torsional angles. Thus, the S(1) LE-TICT equilibration and dual fluorescence will be controlled by (a) the S(1) torsional reaction path and (b) the position along the amino group twist coordinate where the S(2)/S(1) CT-LE radiationless decay occurs. For DMABN, population of LE and TICT can occur because the two species have similar stabilities. However, in ABN, the equilibrium lies in favor of LE, as a TICT state was found at much higher energy with a low reaction barrier toward LE. This explains why dual fluorescence cannot be observed in ABN. The S(1)-->S(0) deactivation channel accessible from the LE state was also studied.     

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.     

Solvent effect on the singlet excited-state dynamics of 5-fluorouracil in acetonitrile as compared with water

The excited-state dynamics of 5-fluorouracil in acetonitrile has been investigated by femtosecond fluorescence upconversion spectroscopy in combination with quantum chemistry TD-DFT calculations ((PCM/TD-PBE0). Experimentally, it was found that when going from water to acetonitrile solution the fluorescence decay of 5FU becomes much faster. The calculations show that this is related to the opening of an additional decay channel in acetonitrile solution since the dark n/pi* excited state becomes near degenerate with the bright pi/pi* state, forming a conical intersection close to the Franck-Condon region. In both solvents, a S1-S0 conical intersection, governed by the out-of-plane motion of the fluorine atom, is active, allowing an ultrafast internal conversion to the ground state.     

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.     

Ultrafast deactivation of an excited cytosine-guanine base pair in DNA

Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy-lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited-state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier. The origin of this effect lies in the temporal electrostatic stabilization of dipole in the charge-transfer state in DNA.     

Mechanism of the N-cyclopropylimine-1-pyrroline photorearrangement

We present here a combined experimental and computational investigation into the photorearrangement of N-cyclopropylimines to yield pyrrolines. We show that the photochemistry, regiochemistry, and stereochemistry of the reaction can be understood in terms of a mechanism involving barrierless evolution in three different (S(2), S(1), S(0)) singlet states and sequential decay through two different (S(2)/S(1), and S(1)/S(0)) conical intersection funnels. We provide evidence that the reaction mechanism involves the generation of a nonequilibrated (i.e., transient) excited state diradical, whose decay can lead not only to pyrrolines but also to cyclopropylimine isomers. It is concluded that the reaction outcome depends on the details of the structure of such transient diradical and on the nature of the dynamics of its decay through the S(1)/S(0) conical intersection.     

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.     

Nonadiabatic dynamics of uracil: population split among different decay mechanisms

Nonadiabatic dynamics simulations performed at the state-averaged CASSCF method are reported for uracil. Supporting calculations on stationary points and minima on the crossing seams have been performed at the MR-CISD and CASPT2 levels. The dominant mechanism is characterized by relaxation into the S(2) minimum of ππ* character followed by the relaxation to the S(1) minimum of nπ* character. This mechanism contributes to the slower relaxation with a decay constant larger than 1.5 ps, in good agreement with the long time constants experimentally observed. A minor fraction of trajectories decay to the ground state with a time constant of about 0.7 ps, which should be compared to the experimentally observed short constant. The major part of trajectories decaying with this time constant follows the ππ* channel and hops to the ground state via an ethylenic conical intersection. A contribution of the relaxation proceeding via a ring-opening conical intersection was also observed. The existence of these two latter channels together with a reduced long time constant is responsible for a significantly shorter lifetime of uracil compared to that of thymine.     

A unified perspective on the hydrogen atom transfer and proton-coupled electron transfer mechanisms in terms of topographic features of the ground and excited potential energy surfaces as exemplified by the reaction between phenol and radicals

The relation between the hydrogen atom transfer (HAT) and proton-coupled electron transfer (PCET) mechanisms is discussed and is illustrated by multiconfigurational electronic structure calculations on the ArOH + R(*) --> ArO(*) + RH reactions. The key topographic features of the Born-Oppenheimer potential energy surfaces that determine the predominant reaction mechanism are the conical intersection seam of the two lowest states and reaction saddle points located on the shoulders of this seam. The saddle point corresponds to a crossing of two interacting valence bond states corresponding to the reactant and product bonding patterns, and the conical intersection corresponds to the noninteracting intersection of the same two diabatic states. The locations of mechanistically relevant conical intersection structures and relevant saddle point structures are presented for the reactions between phenol and the N- and O-centered radicals, (*)NH2 and (*)OOCH3. Points on the conical intersection of the ground doublet D0 and first excited doublet D1 states are found to be in close geometric and energetic proximity to the reaction saddle points. In such systems, either the HAT mechanism or both the HAT mechanism and the proton-coupled electron transfer (PCET) mechanism can take place, depending on the relative energetic accessibility of the reaction saddle points and the D0/D1 conical intersection seams. The discussion shows how the two mechanisms are related and how they blend into each other along intermediate reaction paths. The recognition that the saddle point governing the HAT mechanism is on the shoulder of the conical intersection governing the PCET mechanism is used to provide a unified view of the competition between the two mechanisms (and the blending of the two mechanisms) in terms of the prominent and connected features of the potential energy surface, namely the saddle point and the conical intersection. The character of the dual mechanism may be understood in terms of the dominant valence bond configurations of the intersecting states, which are zero-order approximations to the diabatic states.     

Stairway to the conical intersection: a computational study of the retinal isomerization

The potential-energy surface of the first excited state of the 11-cis-retinal protonated Schiff base (PSB11) chromophore has been studied at the density functional theory (DFT) level using the time-dependent perturbation theory approach (TDDFT) in combination with Becke's three-parameter hybrid functional (B3LYP). The potential-energy curves for torsion motions around single and double bonds of the first excited state have also been studied at the coupled-cluster approximate singles and doubles (CC2) level. The corresponding potential-energy curves for the ground state have been calculated at the B3LYP DFT and second-order M?ller-Plesset (MP2) levels. The TDDFT study suggests that the electronic excitation initiates a turn of the beta-ionone ring around the C6-C7 bond. The torsion is propagating along the retinyl chain toward the cis to trans isomerization center at the C11=C12 double bond. The torsion twist of the C10-C11 single bond leads to a significant reduction in the deexcitation energy indicating that a conical intersection is being reached by an almost barrierless rotation around the C10-C11 single bond. The energy released when passing the conical intersection can assist the subsequent cis to trans isomerization of the C11=C12 double bond. The CC2 calculations also show that the torsion barrier for the twist of the retinyl C10-C11 single bond adjacent to the isomerization center almost vanishes for the excited state. Because of the reduced torsion barriers of the single bonds, the retinyl chain can easily deform in the excited state. Thus, the CC2 and TDDFT calculations suggest similar reaction pathways on the potential-energy surface of the excited state leading toward the conical intersection and resulting in a cis to trans isomerization of the retinal chromophore. According to the CC2 calculations the cis to trans isomerization mechanism does not involve any significant torsion motion of the beta-ionone ring.     

The photochemistry of formaldehyde: internal conversion and molecular dissociation in a single step?

A novel, nonadiabatic reaction path for H2 + CO molecular dissociation of formaldehyde via an extended S1/S0 conical intersection seam has been mapped out using the CAS-SCF method with a full valence active space (10 electrons, 9 orbitals). Two conical intersection geometries have been optimized, CsCoIn, a saddle point in the intersection space, and C1CoIn, which is the lowest-energy crossing point. A minimum-energy path connecting these points along a seam has also been characterized. In addition to the conventional and "roaming-atom" mechanisms--where internal conversion takes place before ground-state dissociation--we suggest that a strictly nonadiabatic mechanism can operate, where internal conversion and dissociation take place in concert.     

Photochemical E/Z isomerization of aryl-substituted methanimines. AM1-CI potential energy curves along the CN double bond twisting and the N atom in-plane inversion

Potential energy curves of aryl-substituted methanimines along the C = N bond twisting and the in-plane inversion of the N atom were obtained by AM1-SDCI calculations, and the photochemical E/Z isomerization paths were elucidated. An aromatic group introduced at the C atom of the C = N group has little effect on the S₀ surface for the inversion and rotation paths, while it has a significant effect on the shape of the T₁ curve along the rotation path. It is suggested that phenylmethanimine and 2-anthrylmethanimine undergo photoisomerization by the rotation mechanism. The methoxy group introduced at the N atom raises the inversion barrier on the S₀ state, while it lowers the 90°-twistedT₁ energy. TheT₁ potential energy curves of N-methoxy-1-phenylethanimine and N-methoxy-1-(2-anthryl)ethanimine indicate that the former undergoes a two-way isomerization and the latter a one-way isomerization by the rotation mechanism, which is in accordance with experiment.     

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.