Acid-base equilibriums of lumichrome and its 1-methyl, 3-methyl, and 1,3-dimethyl derivatives

Lumichrome photophysical properties at different pH were characterized by UV-vis spectroscopy and steady-state and time-resolved fluorescence techniques, in four forms of protonation/deprotonation: neutral form, two monoanions, and dianion. The excited-state lifetimes of these forms of lumichrome were measured and discussed. The results were compared to those obtained for similar forms of alloxazine and/or isoalloxazine, and also to those of 1-methyl- and 3-methyllumichrome and 1,3-dimethyllumichrome. The absorption, emission, and synchronous spectra of lumichrome, 1-methyl- and 3-methyllumichrome, and 1,3-dimethyllumichrome at different pH were measured and used in discussion of fluorescence of neutral and deprotonated forms of lumichrome. The analysis of steady-state and time-resolved spectra and the DFT calculations both predict that the N(1) monoanion and the N(1,3) dianion of lumichrome have predominantly isoalloxazinic structures. Additionally, we confirmed that neutral lumichrome exists in its alloxazinic form only, in both the ground and the excited state. We also confirmed the existence and the alloxazinic structure of a second N(3) monoanion. The estimated values of pK(a) = 8.2 are for the equilibrium between neutral lumichrome and alloxazinic and isoalloxazinic monoanions, with proton dissociation from N(1)-H and N(3)-H groups proceeding at the almost the same pH, while the second value pK(a) = 11.4 refers to the formation of the isoalloxazinic dianion in the ground state.     

Photophysics,Excited‐state Double‐Proton Transfer and Hydrogen‐bonding Properties of 5‐Deazaalloxazines

The photophysical properties of 5‐deazaalloxazine and 1,3‐dimethyl‐5‐deazaalloxazine were studied in different solvents. These compounds have higher values of fluorescence quantum yields and longer fluorescence lifetimes, compared to those obtained for their alloxazine analogs. Electronic structure and S₀–S_i transitions were investigated using the ab initio methods [MP2, CIS(D), EOM‐CCSD] with the correlation‐consistent basis sets. Also the time‐dependent density functional theory (TD‐DFT) has been employed. The lowest singlet excited states of 5‐deazaalloxazine and 1,3‐dimethyl‐5‐deazaalloxazine are predicted to have the π, π* character, whereas similar alloxazines have two close‐lying π, π* and n, π* transitions. Experimental steady‐state and time‐resolved spectral studies indicate formation of an isoalloxazinic excited state via excited‐state double‐proton transfer (ESDPT) catalyzed by an acetic acid molecule that forms a hydrogen bond complex with the 5‐deazaalloxazine molecule. Solvatochromism of both 5‐deazaalloxazine and its 1,3‐dimethyl substituted derivative was analyzed using the Kamlet–Taft scale and four‐parameter Catalán solvent scale. The most significant result of our studies is that the both scales show a strong influence of solvent acidity (hydrogen bond donating ability) on the emission properties of these compounds, indicating the importance of intermolecular solute–solvent hydrogen‐bonding interactions in their excited state.     

Effect of Solvent Polarity on Excited-State Double Proton Transfer Process of 1, 5-Dihydroxyanthraquinone

The excited-state double proton transfer (ESDPT) properties of 1, 5-dihydroxyanthraquinone (1, 5-DHAQ) in various solvents were investigated using femtosecond transient absorption spectroscopy and the DFT/TDDFT method. The steady-state fluorescence spectra in toluene, tetrahydrofuran (THF) and acetonitrile (ACN) solvents presented that the solvent polarity has an effect on the position of the ESDPT fluorescence emission peak for the 1, 5-DHAQ system. Transient absorption spectra show that the increasing polarity of the solvent accelerates the rate of excited state dynamics. Calculated potential energy curves analysis further verified the experimental results. The ESDPT barrier decreases gradually with the increase of solvent polarity from toluene, THF to ACN solvent. It is convinced that the increase of solvent polarity can promote the occurrence of the ESDPT dynamic processes for the 1, 5-DHAQ system. This work clarifies the mechanism of the influence of solvent polarity on the ESDPT process of 1, 5-DHAQ, which provides novel ideas for design and synthesis of new hydroxyanthraquinone derivatives.     

Excited-state double-proton transfer in the 7-azaindole dimer in the gas phase. 3. Reaction mechanism studied by picosecond time-resolved REMPI spectroscopy

The excited-state double-proton-transfer (ESDPT) reaction in the jet-cooled 7-azaindole dimer (7AI2) has been investigated with picosecond time-resolved resonance-enhanced multiphoton ionization spectroscopy. The observed decay profiles of 7AI2 by exciting the origin and the intermolecular stretch fundamental in the S1 state are well reproduced by single-exponential functions with time constants of 1.9 +/- 0.9 ps and 860 +/- 300 fs, respectively. This result provides clear evidence of the concerted mechanism of ESDPT in 7AI2.     

In search of excited-state proton transfer in the lumichrome dimer in the solid state: theoretical and experimental approach

Quantum chemical density functional theory (DFT) calculations and spectral data were employed to investigate the possibility of the excited-state double proton transfer (ESDPT) in lumichrome crystals. The calculations in a lumichrome dimer predict a transfer of a proton in the first excited state, leading to a cation-anion pair. The presently reported X-ray structure of 1,3-dimethyllumichrome and its complex solid-state luminescence indicate that also in this molecule intermolecular hydrogen bonds might be involved in the photophysics. The long-wavelength emission in lumichrome crystals peaked at 530 nm is attributed to excited-state proton transfer, whereas a wider emission band in methylated lumichrome derivatives peaked at 560 nm is attributed to ions formed upon photoexcitation of the crystals.     

The invalidity of intermolecular proton transfer triggered twisted intramolecular charge transfer in excited state for 2-(4′-diethylamino-2′-hydroxyphenyl)-1H-imidazo-[4,5-b]pyridine

The excited-state dynamics of the excited-state proton transfer and intramolecular twisted charge transfer (TICT) reactions of a molecular photoswitch 2-(4′-diethylamino-2′-hydroxyphenyl)-1H-imidazo-[4,5-b]pyridine (DHP) in aprotic and alcoholic solvents have been theoretically investigated by using time-dependent density functional theory. The excited-state intramolecular proton transfer (ESIPT) reaction of DHP proceeding upon excitation in all the solvents has been confirmed, and the dual emission has been assigned to the enol and keto forms of DHP. However, for methanol and ethanol solvents within strong hydrogen-bonded capacity, the intermolecular hydrogen bonds between DHP and methanol/ethanol would promote an excited-state double proton transfer (ESDPT) along the hydrogen-bonded bridge. Importantly, the previous proposed ESDPT-triggered TICT mechanism of DHP in methanol and ethanol was not supported by our calculations. The twist motion would increase the total energy of the system for both the products of ESIPT and ESDPT. According to the calculations of the transition states, the ESDPT reaction occurs much easier in keto form generated by ESIPT. Therefore, a sequential ESIPT and ESDPT mechanism of DHP in methanol and ethanol has been reasonably proposed.     

Excited-state double-proton transfer in a model DNA base pair: Resolution for stepwise and concerted mechanism controversy in the 7-azaindole dimer revealed by frequency- and time-resolved spectroscopy

The excited-state double-proton transfer (ESDPT) reaction in the dual hydrogen-bonded 7-azaindole dimer (7AI₂) in a supersonic jet expansion has been extensively studied with various laser spectroscopic methods and quantum chemistry calculations by many groups. This article reviews the results and discussions associated with stepwise and concerted mechanism controversy on ESDPT of 7AI₂ together with the excited-state dynamics associated with ESDPT.     

Excited-State Double Proton Transfer in 1-Azacarbazole Hydrogen Bonded Complexes

Excited-state double proton transfer (ESDPT) has been studied in a variety of 1-azacarbazole (1AC) hosted hydrogen-bonded complexes. In 1 AC/carboxylic acids hydrogen bonded complexes, large association constants of > 1.0 × 10⁴ M^?1 are measured in the ground state and the rate of ESDPT is » 5.0 × 10⁹ s^?1, resulting in a dominant proton-transfer tautomer emission. In several 1 AC/lactam hydrogen bonded complexes, however, spectral and dynamic results show the existence of a fast excited-state equilibrium between normal and proton-transfer tautomer states. The result can be tentatively rationalized by a non-catalytic ESDPT mechanism incorporating tautomerization energy of the guest molecule.     

Theoretical study of the excited‐state double proton transfer in the (3‐methyl‐7‐azaindole)‐(7‐azaindole) heterodimer

Excited‐state double proton transfer (ESDPT) in the (3‐methyl‐7‐azaindole)‐(7‐azaindole) heterodimer is theoretically investigated by the long‐range corrected time‐dependent density functional theory method and the complete‐active‐space second‐order perturbation theory method. The calculated potential energy profiles exhibit a lower barrier for the concerted mechanism in the locally excited state than for the stepwise mechanism through the charge‐transfer state. This result suggests that the ESDPT in the isolated heterodimer is likely to follow the former mechanism, as has been exhibited for the ESDPT in the homodimer of 7‐azaindole. © 2012 Wiley Periodicals, Inc.     

Solvent effect on the excited-state proton transfer of 7-hydroxyquinoline along a hydrogen-bonded ethanol dimer

We have studied the solvent effect on structures and potential energy surfaces along proton transfer in the ground and the excited states of 7-hydroxyquinoline interacting with an ethanol dimer using ab initio calculations. The proton transfer is forbidden in the ground state not only in vacuum but also in solvents of n-heptane, ethanol, and dimethyl sulfoxide. In the excited state, although the proton transfer is forbidden in vacuum, it is possible in solvent due to its greatly reduced barrier (～10 kcal mol(-1)) and highly stabilized product. It has also been found from the calculations that the proton-transfer barrier in the excited state decreases as the dielectric constant of a solvent increases. Our calculations are consistent with experimental results that the proton transfer does not take place in the ground state and that the excited-state proton-transfer rate increases as the solvent polarity increases. Our calculated absorption and emission properties are in excellent agreement with experimental results. Projection factors (reflecting geometrical change from the ground state to the excited state) and reorganization energies for several low frequency vibrations in connection with the excited-state proton transfer are discussed as well.     

Cyano Analogues of 7‐Azaindole: Probing Excited‐State Charge‐Coupled Proton Transfer Reactions in Protic Solvents

The interplay between excited‐state charge and proton transfer reactions in protic solvents is investigated in a series of 7‐azaindole (7AI) derivatives: 3‐cyano‐7‐azaindole (3CNAI), 5‐cyano‐7‐azaindole (5CNAI), 3,5‐dicyano‐7‐azaindole (3,5CNAI) and dicyanoethenyl‐7‐azaindole (DiCNAI). Similar to 7AI, 3CNAI and 3,5CNAI undergo methanol catalyzed excited‐state double proton transfer (ESDPT), resulting in dual (normal and proton transfer) emission. Conversely, ESDPT is prohibited for 5CNAI and DiCNAI in methanol, as supported by a unique normal emission with high quantum efficiency. Instead, the normal emission undergoes prominent solvatochromism. Detailed relaxation dynamics and temperature dependent studies are carried out. The results conclude that significant excited‐state charge transfer (ESCT) takes place for both 5CNAI and DiCNAI. The charge‐transfer specie possesses a different dipole moment from that of the proton‐transfer tautomer species. Upon reaching the equilibrium polarization, there exists a solvent‐polarity induced barrier during the proton‐transfer tautomerization, and ESDPT is prohibited for 5CNAI and DiCNAI during the excited‐state lifespan. The result is remarkably different from 7AI, which is also unique among most excited‐state charge/proton transfer coupled systems studied to date.     

Theoretical studies for excited-state tautomerization in the 7-azaindole-(CH3OH)n (n = 1 and 2) complexes in the gas phase

The excited-state tautomerization of 7-azaindole (7AI) complexes bonded with either one or two methanol molecule(s) was studied by systematic quantum mechanical calculations in the gas phases. Electronic structures and energies for the reactant, transition state (TS), and product were computed at the complete active space self-consistent field (CASSCF) levels with the second-order multireference perturbation theory (MRPT2) to consider the dynamic electron correlation. The time-dependent density functional theory (TDDFT) was also used for comparison. The excited-state double proton transfer (ESDPT) in 7AI-CH(3)OH occurs in a concerted but asynchronous mechanism. Similarly, such paths are also found in the two transition states during the excited-state triple proton transfer (ESTPT) of the 7AI-(CH(3)OH)(2) complex. In the first TS, the pyrrole ring proton first migrated to methanol, while in the second the methanol proton moved first to the pyridine ring. The CASSCF level with the MRPT2 correction showed that the former path was much preferable to the latter, and the ESDPT is much slower than the ESTPT. Additionally, the vibrational-mode enhanced tautomerization in the 7AI-(CH(3)OH)(2) complex was also studied. We found that the excitation of the low-frequency mode shortens the reaction path to increase the tautomerization rate. Overall, most TDDFT methods used in this study predicted different TS structures and barriers from the CASSCF methods with MRPT2 corrections.     

Ultrafast optical pump-probe studies of the cytochrome b(6)f complex in solution and crystalline states

The cytochrome b6f complex of oxygenic photosynthesis contains a single chlorophyll a (Chl a) molecule whose function is presently unknown. The singlet excited state of the Chl a molecule is quenched by the surrounding protein matrix, and thus the Chl a molecule in the b6f complex may serve as an exceptionally sensitive probe of the protein structure. For the first time, singlet excited-state dynamics were measured in well-diffracting crystals using femtosecond time-resolved optical pump-probe methodology. Lifetimes of the Chl a molecule in crystals of the cytochrome b6f complex having different space groups were 3-6 times longer than those determined in detergent solutions of the b6f. The observed differences in excited state dynamics may arise from small (1-1.5 A) changes in the local protein structure caused by crystal packing. The Chl a excited state lifetimes measured in the dissolved cytochrome b6f complexes from several different species are essentially the same, in spite of differences in the local amino acid sequences around the Chl a. This supports an earlier hypothesis that the short excited state lifetime of Chl a is critical for the function of the b6f complex.     

PHOTOTAUTOMERISM OF LUMICHROME IN METHANOL-ACETIC ACID MIXTURES. A STEADY-STATE AND TIME-RESOLVED FLUORESCENCE STUDY

Abstract— Fluorescence lifetimes of 3-methyllumichrome dual fluorescence at 440 nm (alloxazinic) and 520 nm (isoalloxazinic) and of 1,3-dimethyllumichrome, which is unable to phototautomerize, at the same wavelengths have been measured in methanol-acetic acid mixtures. The fluorescence decays of both lumichromes studied are exponential and the phototautomeric fluorescence of 3-methyllumichrome is created within approximately 50 ps. From the initial values of about 0.9 ns (alloxazinic) and 6.4 ns (isoalloxazinic) in 5% acetic acid the lifetimes shorten considerably with increasing acid concentration reaching 0.2 ns for 1,3-dimethyllumichrome and 80 ps for the alloxazinic and 2.4 ns for the isoalloxazinic form of 3-methyllumichrome in pure acid. Static and dynamic fluorescence quenching constants were estimated. The differences in the quenching rate and equilibrium constants of both lumichromes are interpreted in terms of different equilibria of hydrogen bond formation in ground and excited state as influenced by steric effects of the methyl substituent at the N-l position in 1,3-dimethyl. The hydrogen bonding atN–10 of 3-methyllumichrome with acetic acid is a prerequisite for additional hydrogen bonding at N-l enabling excited state proton transfer.     

Excited-State Double Proton Transfer of 1, 8-Dihydroxy-2-Naphthaldehyde: a MS-CASPT2//CASSCF Study

Excited-state double proton transfer (ESDPT) is a controversial issue which has long been plagued with theoretical and experimental communities. Herein, we took 1, 8-dihydroxy-2-naphthaldehyde (DHNA) as a prototype and used combined complete active space self-consistent field (CASSCF) and multi-state complete active-space second-order perturbation (MS-CASPT2) methods to investigate ESDPT and excited-state deactivation pathways of DHNA. Three different tautomer minima of S1-ENOL, S1-KETO-1, and S1-KETO-2 and two crucial conical intersections of S1S0-KETO-1 and S1S0-KETO-2 in and between the S₀ and S₁ states were obtained. S1-KETO-1 and S1-KETO-2 should take responsibility for experimentally observing dual-emission bands. In addition, two-dimensional potential energy surfaces (2D-PESs) and linear interpolated internal coordinate paths connecting relevant structures were calculated at the MS-CASPT2//CASSCF level and confirmed a stepwise ESDPT mechanism. Specifically, the first proton transfer from S1-ENOL to S1-KETO-1 is barrierless, whereas the second one from S1-KETO-1 to S1-KETO-2 demands a barrier of ca. 6.0 kcal/mol. The linear interpolated internal coordinate path connecting S1-KETO-1 (S1-KETO-2) and S1S0-KETO-1 (S1S0-KETO-2) is uphill with a barrier of ca. 12.0 kcal/mol, which will trap DHNA in the S₁ state while therefore enabling dual-emission bands. On the other hand, the S₁/S₀ conical intersections would also prompt the S₁ system to decay to the S₀ state, which could be to certain extent suppressed by locking the rotation of the C5$-$C8$-$C9$-$O10 dihedral angle. These mechanistic insights are not only helpful for understanding ESDPT but also useful for designing novel molecular materials with excellent photoluminescent performances.     

The Host/Guest Type of Excited‐State Proton Transfer; a General Review

Contemporary progress regarding guest/host types of excited‐state double proton transfer has been reviewed, among which are the biprotonic transfer within doubly H‐bonded host/guest complexes, the transfer through a solvent bridge relay, the intramolecular double proton transfer and solvation dynamics coupled proton transfer. Of particular emphases are the photophysical and photochemical properties of excited‐state double proton transfer (ESDPT) in 7‐azaindole and its corresponding analogues. From the chemical aspect, two types of ESDPT reaction, namely the catalytic and non‐catalytic types of ESDPT, have been classified and reviewed separately. For the case of static host/guest hydrogen‐bonded complexes both hydrogen‐bonding strength and configuration (i.e. geometry) play key roles in accounting for the reaction dynamics. In addition to the dynamical concern, excited‐state thermodynamics are of importance to fine‐tune the proton transfer reaction in the non‐catalytic host/guest type of ESDPT. The mechanisms of protic solvent assisted ESDPT, depending on host molecules and proton‐transfer models, have been reviewed where the plausible resolution is deduced. Particular attention has been given to the excited‐state proton transfer dynamics in pure water, aiming at its future perspective in biological applications. Finally, the differentiation in mechanism between solvent diffusive reorganization and solvent relaxation to affect the host/guest ESPT dynamics is made and discussed in de tail.     

Coherent nuclear wavepacket motions in ultrafast excited-state intramolecular proton transfer: sub-30-fs resolved pump-probe absorption spectroscopy of 10-hydroxybenzo[h]quinoline in solution

The dynamics of the excited-state intramolecular proton transfer of 10-hydroxybenzo[h]quinoline (10-HBQ) and the associated coherent nuclear motion were investigated in solution by femtosecond absorption spectroscopy. Sub-picosecond transient absorption measurements revealed spectral features of the stimulated emission and absorption of the keto excited state (the product of the reaction). The stimulated emission band appeared in the 600-800-nm region, corresponding to the wavelength region of the steady-state keto fluorescence. It showed successive temporal changes with time constants of 350 fs and 8.3 ps and then disappeared with the lifetime of the keto excited state (260 ps). The spectral feature of the stimulated emission changed in the 350-fs dynamics, which was likely assignable to the intramolecular vibrational energy redistribution in the keto excited state. The 8.3-ps change caused a spectral blue shift and was attributed to the vibrational cooling process. The excited-state absorption was observed in the 400-600-nm region, and it also showed temporal changes characterized by the 350-fs and 8.3-ps components. To examine the coherent nuclear dynamics (nuclear wavepacket motion) in excited-state 10-HBQ, we carried out pump-probe measurements of the stimulated emission and absorption signals with time resolution as good as 27 fs. The obtained data showed substantially modulated signals due to the excited-state vibrational coherence up to a delay time of several picoseconds after photoexcitation. This means that the vibrational coherence created by photoexcitation in the enol excited state is transferred to the product. Fourier transform analysis indicated that four frequency components in the 200-700-cm(-1) region contribute to the oscillatory signal, corresponding to the coherent nuclear motions in excited-state 10-HBQ. Especially, the lowest-frequency mode at 242 cm(-1) is dephased significantly faster than the other three modes. This observation was regarded as a manifestation that the nuclear motion of the 242-cm(-1) mode is correlated with the structural change of the molecule associated with the reaction (the reaction coordinate). The 242-cm(-1) mode observed in excited-state 10-HBQ was assigned to a vibration corresponding to the ground-state vibration at 243 cm(-1) by referring to the results of resonance Raman measurements and density functional calculations. It was found that the nuclear motion of this lowest-frequency mode involves a large displacement of the OH group toward the nitrogen site as well as in-plane skeletal deformation that assists the oxygen and nitrogen atoms to come closer to each other. We discuss the importance of the nuclear wavepacket motion on a multidimensional potential-energy surface including the vibrational coordinate of the low-frequency modes.     

Excited-state proton transfer through water bridges and structure of hydrogen-bonded complexes in 1H-pyrrolo[3,2-h]quinoline: adiabatic time-dependent density functional theory study

Proton transfer reaction is studied for 1H-pyrrolo[3,2-h]quinoline-water complexes (PQ-(H(2)O)(n), n = 0-2) in the ground and the lowest excited singlet states at the density functional theory (DFT) level. Cyclic hydrogen-bonded complexes are considered, in which water molecules form a bridge connecting the proton donor (pyrrole NH group) and acceptor (quinoline nitrogen) atoms. To understand the effect of the structure and length of water bridges on the excited-state tautomerization in PQ, the potential energy profile of the lowest excited singlet state is calculated adiabatically by the time-dependent DFT (TDDFT) method. The S(0) --> S(1) excitation of PQ is accompanied by significant intramolecular transfer of electron density from the pyrrole ring to the quinoline fragment, so that the acidity of the N-H group and the basicity of the nitrogen atom of the quinoline moiety are increased. These excited-state acid-base changes introduce a driving force for the proton transfer reaction. The adiabatic TDDFT calculations demonstrate, however, that the phototautomerization requires a large activation energy in the isolated PQ molecule due to a high energy barrier separating the normal form and the tautomer. In the 1:1 cyclic PQ-H(2)O complex, the energy barrier is dramatically reduced, so that upon excitation of this complex the tautomerization can occur rapidly in one step as concerted asynchronous movements of the two protons assisted by the water molecule. In the PQ-(H(2)O)(2) solvate two water molecules form a cyclic bridge with sterically strained and unfavorable hydrogen bonds. As a result, some extra activation energy is needed for initiating the proton dislocation along the longer hydrogen-bond network. The full tautomerization in this complex is still possible; however, the cooperative proton transfer is found to be highly asynchronous. Large relaxation and reorganization of the hydrogen-bonded water bridge in PQ-(H(2)O)(2) are required during the proton translocation from the pyrrole NH group to the quinoline nitrogen; this may block the complete tautomerization in this type of solvate.     

Excited-state double-proton transfer in the 7-azaindole dimer in the gas phase. 2. Cooperative nature of double-proton transfer revealed by H/D kinetic isotopic effects

The dispersed fluorescence (DF) spectra of the 7-azaindole dimer (7AI2) and deuterated dimers 7AI2-hd and 7AI2-dd, where hd and dd indicate the deuteration of an imino proton and two imino protons, have been measured in a supersonic free jet expansion. The undeuterated 7AI2-hh dimer exhibits only the tautomer fluorescence, but both the normal and tautomer fluorescence have been detected by exciting the origins of 7AI2-h*d, 7AI2-hd* and 7AI2-dd in the S1-S0 region, where h* and d* indicate the localization of the excitation on 7AI-h or 7AI-d moiety. The DF spectra indicate that 7AI2-h*d and 7AI2-hd* undergo excited-state proton/deuteron transfer (ESPDT), while excited-state double-deuteron transfer (ESDDT) occurs in 7AI2-dd. The H/D kinetic isotopic effects on ESDPT have been investigated by measuring the intensity ratios of the normal fluorescence to the tautomer fluorescence. The ESPDT rate is about 1/60th of the ESDPT rate, and the ESDDT rate is about 1/12th of the ESPDT rate, where ESPDT rate is an average of the rates for 7AI2-h*d and 7AI2-hd*. The observed H/D kinetic isotope effects imply that the ESDPT reaction of 7AI2 has a "cooperative" nature; i.e., the motion of the two moving protons strongly couples each other through the electron motions. The difference in the estimated ESPDT reaction rates, 9.8 x 10(9) and 6.9 x 109 s(-1) for 7AI2-h*d and 7AI2-hd*, respectively, is consistent with the concerted mechanism rather than the stepwise mechanism.     

Nonradiative deactivation mechanisms of uracil, thymine, and 5-fluorouracil: a comparative ab initio study

The mechanisms of the ultrafast nonradiative deactivation of uracil and its substituted derivatives thymine (5-methyluracil) and 5-fluorouracil after absorption of UV light are explored and compared by means of ab initio multistate (MS) CASPT2 calculations. The MS-CASPT2 method is applied for the calculation of potential energy profiles, especially for the geometry optimization in the electronically excited state, with the aim of an accurate prediction of deactivation pathways. The resulting energy curves of each molecule exhibit that the conical intersection between the (1)ππ* and ground states is accessible via small energy barriers from the minimum in the (1)ππ* state as well as from that in the (1)nπ* state. The barrier of 5-fluorouracil in the (1)ππ* state is calculated to be definitely higher than those of uracil and thymine, which is consistent with experiments and suggests that the elongation of the excited-state lifetime of uracil by fluorine substitution is significantly contributed from intrinsic electronic effect of the molecule. However, no evidence of the experimentally observed longer excited-state lifetime of thymine than uracil is found in the presently calculated MS-CASPT2 potential energy curves in the (1)ππ* and (1)nπ* states, implying nonnegligible contribution of other factors such as solvation effect and substituent mass to the photoinduced dynamics of uracil derivatives.