Ultrafast Twisting Dynamics of Thioflavin‐T: Spectroscopy of the Twisted Intramolecular Charge‐Transfer State

Understanding the excited‐state properties of thioflavin‐T (ThT) has been of immense importance, because of its efficient amyloid‐sensing ability related to neurodegenerative disorders. The excited‐state dynamics of ThT is studied by using sub‐pico‐ and nanosecond time‐resolved transient absorption techniques as well as density functional theory (DFT)/time‐dependent DFT calculations. Barrierless twisting around the central C?C bond between two aromatic moieties is the dominant process that contributes to the ultrafast dynamics of the S₁ state. The spectroscopic properties of the intramolecular charge‐transfer state are characterized for the first time. The energetics of the S₀ and S₁ states has also been correlated with the experimentally observed spectroscopic parameters and structural dynamics. A longer‐lived transient state populated with a very low yield has been characterized as the triplet state.     

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.     

Time-dependent density functional theory (TDDFT) study of the excited charge-transfer state formation of a series of aromatic donor-acceptor systems

Singlet excitation energy calculations for a series of acceptor para-substituted N,N-dimethyl-anilines that are dual (4-(N,N-dimethylamino)benzonitrile, 4DMAB-CN, 4-(N,N-dimethylamino)benzaldhyde, 4DMAB-CHO, 1-methyl-7-cyano-2,3,4,5-tetrahydro-1H-1-benzazepine, NMC7) and nondual (4-aminobenzonitrile, 4AB-CN, 3-(N,N-dimethylamino)benzonitrile, 3DMAB-CN, and 4-nitro(N,N-dimethyl) aniline, 4DMAB-NO(2)) fluorescent have been performed using time-dependent density functional theory (TDDFT). The B3LYP and MPW1PW91 functionals with a 6-311+G(2d,p) (Bg) basis set have been used to compute excitation energies. Ground-state geometries were optimized using density functional theory (DFT) with both B3LYP and MPW1PW91 functionals combined with a 6-31G(d) basis set. For most of the molecules presented in this study, potential energy surfaces have been computed according to the coordinates related to the three following mechanisms proposed in the literature: twisting, wagging, and planar intramolecular charge transfer (ICT). Comparison of the three models for the different molecules leads to the conclusion that only the twisting ICT model is able to explain the low frequency, strongly solvent-dependent energy band present in the fluorescence spectra. According to this model, the 4AB-CN molecule is calculated to be nondual fluorescent in agreement with the experimental spectra. The single band observed in the fluorescence spectra of TMAB-CN (4-(N,N-dimethylamino)-3,5-(dimethyl)benzonitrile) is due to a large stabilization of the charge-transfer excited state along the twisting coordinate. The nondual fluorescence of the 4DMAB-NO(2) molecule is explained by the same mechanism. In the case of 3DMAB-CN, the single observed emission, which is solvent-dependent, has been assigned to the lowest charge-transfer excited state. The dual fluorescence of 4DMAB-CN and 4DMAB-CHO is explained within the twisting ICT model by a double mechanism (already proposed by Serrano et al.: Serrano-Andrés, L.; Merchán, M.; Roos, B. J.; Lindh, R. J. Am. Chem. Soc. 1995, 117, 3189) that involves the presence of two low-lying states close enough in energy. The observation of dual fluorescence in NMC7, that has been one of the origins of the planar ICT model put forward by Zachariasse et al. (Zachariasse, K.; van der Haar, T.; Hebecker, A.; Leinhos, U.; Kühnle, W. Pure Appl. Chem. 1993, 65, 1745), could be fully understood by a double mechanism within the twisting ICT model. Within the set of investigated molecules, our calculations confirm that the twisting ICT model is the only mechanism acceptable to explain the dual and nondual fluorescence phenomenon. Our calculations are in complete agreement with experimental data.     

The role of the beta-ionone ring in the photochemical reaction of rhodopsin

Time-dependent density functional theory (TDDFT) calculations on the photoabsorption process of the 11-cis retinal protonated Schiff base (PSB) chromophore show that the Franck-Condon relaxation of the first excited state of the chromophore involves a torsional twist motion of the beta-ionone ring relative to the conjugated retinyl chain. For the ground state, the beta-ionone ring and the retinyl chain of the free retinal PSB chromophore form a -40 degrees dihedral angle as compared to -94 degrees for the first excited state. The double bonds of the retinal are shorter for the fully optimized structure of the excited state than for the ground state suggesting a higher cis-trans isomerization barrier for the excited state than for the ground state. According to the present TDDFT calculations, the excitation of the retinal PSB chromophore does not primarily lead to a reaction along the cis-trans torsional coordinate at the C11-C12 bond. The activation of the isomerization center seems to occur at a later stage of the photo reaction. The results obtained at the TDDFT level are supported by second-order M?ller-Plesset (MP2) and approximate singles and doubles-coupled cluster (CC2) calculations on retinal chromophore models; the MP2 and CC2 calculations yield for them qualitatively the same ground state and excited-state structures as obtained in the density functional theory and TDDFT calculations.     

Mechanism of ligand photodissociation in photoactivable [Ru(bpy)2L2]2+ complexes: a density functional theory study

A series of four photodissociable Ru polypyridyl complexes of general formula [Ru(bpy)2L2](2+), where bpy = 2,2'-bipyridine and L = 4-aminopyridine (1), pyridine (2), butylamine (3), and gamma-aminobutyric acid (4), was studied by density functional theory (DFT) and time-dependent density functional theory (TDDFT). DFT calculations (B3LYP/LanL2DZ) were able to predict and elucidate singlet and triplet excited-state properties of 1-4 and describe the photodissociation mechanism of one monodentate ligand. All derivatives display a Ru --> bpy metal-to-ligand charge transfer (MLCT) absorption band in the visible spectrum and a corresponding emitting triplet (3)MLCT state (Ru --> bpy). 1-4 have three singlet metal-centered (MC) states 0.4 eV above the major (1)MLCT states. The energy gap between the MC states and lower-energy MLCT states is significantly diminished by intersystem crossing and consequent triplet formation. Relaxed potential energy surface scans along the Ru-L stretching coordinate were performed on singlet and triplet excited states for all derivatives employing DFT and TDDFT. Excited-state evolution along the reaction coordinate allowed identification and characterization of the triplet state responsible for the photodissociation process in 1-4; moreover, calculation showed that no singlet state is able to cause dissociation of monodentate ligands. Two antibonding MC orbitals contribute to the (3)MC state responsible for the release of one of the two monodentate ligands in each complex. Comparison of theoretical triplet excited-state energy diagrams from TDDFT and unrestricted Kohn-Sham data reveals the experimental photodissociation yields as well as other structural and spectroscopic features.     

Excited state geometry optimizations by analytical energy gradient of long-range corrected time-dependent density functional theory

An analytical excitation energy gradient of long-range corrected time-dependent density functional theory (LC-TDDFT) is presented. This is based on a previous analytical TDDFT gradient formalism, which avoids solving the coupled-perturbed Kohn-Sham equation for each nuclear degree of freedom. In LC-TDDFT, exchange interactions are evaluated by combining the short-range part of a DFT exchange functional with the long-range part of the Hartree-Fock exchange integral. This LC-TDDFT gradient was first examined by calculating the excited state geometries and adiabatic excitation energies of small typical molecules and a small protonated Schiff base. As a result, we found that long-range interactions play a significant role even in valence excited states of small systems. This analytical LC-TDDFT gradient was also applied to the investigations of small twisted intramolecular charge transfer (TICT) systems. By comparing with calculated ab initio multireference perturbation theory and experimental results, we found that LC-TDDFT gave much more accurate absorption and fluorescence energies of these systems than those of conventional TDDFTs using pure and hybrid functionals. For optimized excited state geometries, LC-TDDFT provided fairly different twisting and wagging angles of these small TICT systems in comparison with conventional TDDFT results.     

DFT and TDDFT study related to electron transfer in nonbonded porphine...C60 complexes

Spectroscopic properties of a ground state nonbonded porphine-buckminsterfullerene (H2P...C60) complex are studied in several different relative orientations of C60 with respect to the porphine plane by using the density functional (DFT) and time-dependent density functional (TDDFT) theories. The geometries and electronic structures of the ground states are optimized with the B3LYP and PBE functionals and a SVP basis set. Excitation energies and oscillator strengths are obtained from the TDDFT calculations. The relative orientation of C60 is found to affect the equilibrium distance between H2P and C60 especially in the case of the PBE functional. The excitation energies of different H2P...C60 complexes are found to be practically the same for the same excitations when the B3LYP functional is used but to differ notably when PBE is used in calculations. Existence of the states related to a photoinduced electron transfer within a porphyrin-fullerene dyad is also studied. All calculations predict a formation of an excited charge-transfer complex state, a locally excited donor (porphine) state, as well as a locally excited acceptor (fullerene) state in the investigated H2P...C60 complexes.     

Structural and spectral properties of 4-bromo-1-naphthyl chalcones: a quantum chemical study

The structural and optical properties of 4-bromo-1-naphthyl chalcones (BNC) have been studied by using quantum chemical methods. The density functional theory (DFT) and the singly excited configuration interaction (CIS) methods were employed to optimize the ground and excited state geometries of unsubstituted and substituted BNC with different electron withdrawing and donating groups in both gas and solvent phases. Based on the ground and excited state geometries, the absorption and emission spectra of BNC molecules were calculated using the time-dependent density functional theory (TDDFT) method. The solvent phase calculations were performed using the polarizable continuum model (PCM). The geometrical parameters, vibrational frequencies, and relative stability of cis- and trans-isomers of unsubstituted and substituted BNC molecules have been studied. The results from the TDDFT calculations reveal that the substitution of electron withdrawing and electron donating groups affects the absorption and emission spectra of BNC.     

Identifying the Bond Responsible for the Fluorescence Modulation in an Amyloid Fibril Sensor

An ultrafast intramolecular bond twisting process is known to be the responsible mechanism for the sensing activity of the extensively used amyloid fibril sensor thioflavin T (ThT). However, it is not yet known which one of the two possible single bonds in ThT is actually involved in the twisting process. To resolve this fundamental issue, two derivatives of ThT have been designed and synthesized and subsequently their photophysical properties have been studied in different solvents. It is understood from the present study that the rotation around the central C? C single bond, and not that around the C? N single bond, is primarily responsible for the sensor activity of ThT. Detailed viscosity‐dependent fluorescence studies revealed that the ThT derivative with restricted C? N bond rotation acts as a better sensor than the derivative with free C? N bond rotation. The better sensory activity is directly correlated with a shorter excited‐state lifetime. Results obtained from the photophysical studies of the ThT derivatives have also been supported by the results obtained from quantum chemical calculations.     

Theoretical Study on Resonance Raman Spectra of Tetraoxaporphyrin Dication by TDDFT Calculation

The B state excited resonance Raman scattering of tetraoxaporphyrin dication (TOP²⁺) was theoretically studied with DFT/TDDFT calculations and the sum-over-states approach of polarizability including both the A and B terms contributions. The resonance Raman spectra calculated with PBE1PBE, B3LYP, Cam-B3LYP, and B3LYP-D3 functionals are similar to each other in general, with PBE1PBE and B3LYP being better in reproducing resonance Raman intensities in comparison with the experiment. The calculated relative intensities of the totally symmetric modes are excellently consistent with the experiment. The TDDFT calculations manifested a considerable deformation of the B state along theυ2,υ6, υ7, and υ8 modes, which is responsible for the strong resonance Raman intensities of these modes. The resonance Raman intensities of non-totally symmetric modes were calculated to be weaker than the totally symmetric modes by one or two order of magnitude, whichqualitatively agrees with the experiment. However, the resonance Raman intensity of the υ10 mode (C_βC_β stretch, B_1g symmetry) predicted by TDDFT calculations is unexpectedly small whereas that of the υ11 mode (symmetric C_αC_m stretch, B_1g symmetry) is too large, which is assumed to be caused by the Jahn-Teller instability for the B state of TOP²⁺.     

Theoretical analysis of excited states and energy transfer mechanism in conjugated dendrimers

The excited states of the phenylene ethynylene dendrimer are investigated comprehensively by various electronic‐structure methods. Several computational methods, including SCS‐ADC(2), TDHF, TDDFT with different functionals (B3LYP, BH&HLYP, CAM‐B3LYP), and DFT/MRCI, are applied in systematic calculations. The theoretical approach based on the one‐electron transition density matrix is used to understand the electronic characters of excited states, particularly the contributions of local excitations and charge‐transfer excitations within all interacting conjugated branches. Furthermore, the potential energy curves of low‐lying electronic states as the functions of ethynylene bonds are constructed at different theoretical levels. This work provides us theoretical insights on the intramolecular excited‐state energy transfer mechanism of the dendrimers at the state‐of‐the‐art electronic‐structure theories. © 2014 Wiley Periodicals, Inc.     

Computational study of thioflavin T torsional relaxation in the excited state

Quantum-chemical calculations of the Thioflavin T (ThT) molecule in the ground S0 and first excited singlet S1 states were carried out. It has been established that ThT in the ground state has a noticeable nonplanar conformation: the torsion angle phi between the benzthiazole and the dimethylaminobenzene rings has been found to be approximately 37 degrees. The energy barriers of the intramolecular rotation appearing at phi = 0 and 90 degrees are quite low: semiempirical AM1 and PM3 methods predict values approximately 700 cm-1 and ab initio methods approximately 1000-2000 cm(-1). The INDO/S calculations of vertical transitions to the S1(abs) excited state have revealed that energy ES1(abs) is minimal for the twisted conformation with phi = 90 degrees and that the intramolecular charge-transfer takes place upon the ThT fragments' rotation from phi = 0 to 90 degrees. Ab initio CIS/RHF calculations were performed to find optimal geometries in the excited S1 state for a series of conformers having fixed phi values. The CIS calculations have predicted a minimum of the S1 state energy at phi approximately 21 degrees; however, the energy values are 1.5 times overestimated in comparison to experimental data. Excited state energy dependence on the torsion angle phi, obtained by the INDO/S method, reveals that ES1(fluor) is minimal at phi = approximately 80-100 degrees, and a plateau is clearly observed for torsion angles ranging from 20 to 50 degrees. On the basis of the calculation results, the following scheme of photophysical processes in the excited S1 state of the ThT is suggested. According to the model, a twisted internal charge-transfer (TICT) process takes place for the ThT molecule in the excited singlet state, resulting in a transition from the fluorescent locally excited (LE) state to the nonfluorescent TICT state, accompanied by torsion angle phi growth from 37 to 90 degrees. The TICT process effectively competes with radiative transition from the LE state and is responsible for significant quenching of the ThT fluorescence in low-viscosity solvents. For viscous solvents or when the ThT molecule is located in a rather rigid microenvironment, for example, when it is bound to amyloid fibrils, internal rotation in the dye molecule is blocked due to steric hindrance, which results in suppression of the LE --> TICT quenching process and in a high quantum yield of fluorescence.     

Computational Assessment of MLCT versus MC Stabilities in First-to-Third-Row d6 Pseudo-Octahedral Transition Metal Complexes

An accurate modeling of metal-to-ligand-charge-transfer (MLCT) and metal-centered (MC) excited state energies is key to predict the photoinduced response in transition metal complexes (TMCs). Herein, the importance of the ground state and excited state reference geometries is addressed for three-prototype d⁶ pseudo-octahedral TMCs, each displaying a different potential energy landscape of MLCT versus MC relative stabilities. Several functionals are used within the time-dependent density functional theory (TDDFT), as well as multireference wave-function theory (MS-CASPT2), applied to [Mn(im)(CO)₃(phen)]⁺, [Ru(im)₂(bpy)₂]²⁺, and [Re(im)(CO)₃(phen)]⁺, (im: imidazole, bpy: bypiridine, phen: phenantroline). The results revel that TDDFT is robust except when using B3LYP functional for first-row d⁶ TMCs. In contrast, MS-CASPT2 calculations are strongly biased in those cases with competitive MLCT/MC states. The results reinforce the reliability of B3LYP to describe the excited states in d⁶ TMCs, but question the validity of assessing the density functional theory (DFT)/TDDFT performance via direct comparison with MS-CASPT2 performed at the same DFT reference geometry as a standard strategy. © 2019 Wiley Periodicals, Inc.     

Low-lying excited states of C120 and C151: a multireference perturbation theory study

Excited states of two 7-aminocoumarin derivatives, coumarin 120 (7-amino-4-methylcoumarin) and coumarin 151 (7-amino-4-trifluoromethylcoumarin), were investigated using generalized multiconfigurational quasidegenerate perturbation theory (GMC-QDPT), multiconfigurational quasidegenerate perturbation theory (MC-QDPT) and time-dependent density functional theory (TDDFT) with the B3LYP and CAM-B3LYP functionals. The absorption and fluorescence spectra of C120 and C151 were calculated. We elucidated the characters of the low-lying states of C120 and C151. The absorption spectra calculated with GMC-QDPT and TDDFT B3LYP agreed well with the experimental data, while for the fluorescence spectra, the TDDFT calculations overestimated the fluorescence spectra compared to GMC-QDPT calculations. Utilizing active spaces with large numbers of electrons and orbitals for reference functions, GMC-QDPT showed a better performance than MC-QDPT with a complete active space self-consistent field (CASSCF) reference of active space with smaller number of electrons and orbitals. In our gas phase calculation, we found that the optimized structures for the first excited states have a planar amino group with a CN single bond, while the amino group is pyramidal in the ground state.     

Torsion-wagging tunneling and vibrational states in hydrazine determined from its ab initio potential energy surface

Geometries, anharmonic vibrations, and torsion-wagging (TW) multiplets of hydrazine and its deuterated species are studied using high-level ab initio methods employing the second-order Mo?ller-Plesset perturbation theory (MP2) as well as the coupled cluster singles and doubles model including connected triple corrections, CCSD(T), in conjunction with extended basis sets containing diffuse and core functions. To describe the splitting patterns caused by tunneling in TW states, the 3D potential energy surface (PES) for the large-amplitude TW modes is constructed. Stationary points in the 3D PES, including equivalent local minima and saddle points are characterized. Using this 3D PES, a flexible Hamiltonian is built numerically and then employed to solve the vibrational problem for TW coupled motion. The calculated ground state r(av) structure is expected to be more reliable than the experimental one that has been determined using a simplified structural model. The calculated fundamental frequencies allowed resolution of the assignment problems discussed earlier in the literature. The determined energy barriers, including the contributions from the small-amplitude vibrations, to the tunneling of the symmetric and antisymmetric wagging mode of 1997 cm(-1) and 3454 cm(-1), respectively, are in reasonable agreement with the empirical estimates of 2072 cm(-1) and 3312 cm(-1), respectively [W. ?odyga et al. J. Mol. Spectrosc. 183, 374 (1997)]. However, the empirical torsion barrier of 934 cm(-1) appears to be overestimated. The ab initio calculations yield two torsion barriers: cis and trans of 744 cm(-1) and 2706 cm(-1), respectively. The multiplets of the excited torsion states are predicted from the refined 3D PES.     

On the ordering of the first two excited electronic states in all-trans linear polyenes

Reported experimental evidence of the relative position of the first two excited electronic states in linear polyenes was carefully examined and compared with that derived from time dependent density functional theory (TDDFT) theoretical calculations performed at the B3LYP level on optimized geometries. The energy values for the first two triplet states 3Bu and 3Ag, obtained from TDDFT calculations, were found to be highly strongly correlated with the experimental values. Also, the theoretical calculations for the electronic transition 1 1Ag --> 1 1Bu were also extremely well correlated with their experimental counterparts; even more important, the three reported experimental data for 1 1Ag --> 2 1Ag transitions in these systems conformed to the correlation for the TDDFT 1 1Ag --> 1 1Bu transition. The first excited electronic state in the linear polyenes studied (from ethene to the compound consisting of 40 ethene units, P40) was found to be 1Bu. The energy gap between the excited states 2 1Ag and 1 1Bu decreased with increasing length of the polyene chain, but not to the extent required to cause inversion, at least up to P40. In the all-trans linear polyenes studied, the widely analyzed energy gap from the ground electronic state to the first excited singlet state for infinitely long chains may be meaningless as, even in P40, it is uncertain whether the ground electronic state continues to be a singlet.     

On the excited states involved in the luminescent probe [Ru(bpy)2dppz]2+

The nature of the excited states of [Ru(bpy)2dppz]2+ has been investigated using density functional theory with the hybrid functional B3LYP. The excitations were studied via linear response theory (TDDFT) and DeltaSCF calculations and the solvent effects were introduced by embedding the molecule in a continuum dielectric medium. It was found that the solvent effects are critical in understanding the nature of the excitations. For the molecule in ethanol, the lowest absorption predicted by TDDFT is a dark state 3pi --> pi with the electron and hole spread over the dppz ligand. Next come the excitations of 3MLCT between the ruthenium and the dppz and finally the 3MLCT excitations between the ruthenium and the bpy ligands not associated with the phenazine. Using deltaSCF calculations two low-lying excited states were identified and the geometry optimized in the presence of the continuum medium. At the optimal geometry the lowest excited state is 3MLCT (Ru --> dppz). The 3pi --> pi state is found only 0.026 eV higher.     

Time-dependent density functional theory study on intramolecular charge transfer and solvent effect of dimethylaminobenzophenone

The lower singlet excited states for dimethylaminobenzophenone have been investigated as a function of the twisting motion with inclusion of solvent effects. Theoretical calculations have been performed using time-dependent density functional theory. The B3LYP and MPW1PW91 functionals with a 6-311+G(2d,p) basis set have been used to compute transition energies. The solvent effects have been described within the polarizable continuum model. Ground-state geometries are optimized using density functional theory with both B3LYP and MPW1PW91 functionals combined with 6-31G(d) basis sets. Vertical absorption energy calculations characterize the lower singlet excited states both in vacuum and in different kinds of solvents. A large redshift of the absorption maximum in the polar solvent suggests an intramolecular charge transfer character of the excited state. We have constructed the potential energy curves of two possible twisting motions of the excited states both in vacuum and in the polar solvent of acetonitrile: the twisting of only the dimethylamino group and the twisting of the dimethylaminophenyl group with respect to the benzoyl group. Both twisting processes predict the formation of the twisted intramolecular charge transfer state associated with the crossing of a low barrier. The presence of the polar solvent significantly changes the shape of the energy curves. Calculated emission energies for both the isolated and the solvated systems show a large Stokes shift between the absorption and fluorescence maxima. Two possible twisting motions produce similar fluorescence spectroscopic consequences. Our results including solvent effects explain the weak "dual-fluorescence" feature of dimethylaminobenzophenone, and imply that the two possible twisting motions may occur in the excited-state relaxation dynamics, but the twisting of the dimethylamino group seems to take place easier.     

Excited-state proton transfer in 7-hydroxy-4-methylcoumarin along a hydrogen-bonded water wire

TDDFT, RI-CC2, and CIS calculations have been performed for the nondissociative excited-state proton transfer (ESPT) in the S1 state of 7-hydroxy-4-methylcoumarin (7H4MC) along a H-bonded water wire of three water molecules bridging the proton donor (OH) and the proton acceptor (C[double bond]O) groups (7H4MC.(H2O)3). The observed structural reorganization in the water-wire cluster is interpreted as a proton-transfer (PT) reaction along the H2O solvent wire. The shift of electron density within the organic chromophore 7H4MC due to the optical excitation appears to be the driving force for ESPT. All the methods used show that the reaction path occurs in the 1pipi* state, and no crossing with a Rydberg-type 1pisigma* state is found. TDDFT and RI-CC2 calculations predict an exoergic reaction of the excited-state enol-to-keto transformation. The S1 potential energy curve reveals well-defined Cs minima of enol- and keto-clusters, separated by a single barrier with a height of 17-20 kcal/mol. After surmounting this barrier, spontaneous PT along the water wire is observed, leading without any further barrier to the keto structure. The TDDFT and RI-CC2 methods appear to be reliable approaches to describe the energy surfaces of ESPT. The CIS method predicts an endoergic ESPT reaction and an energy barrier, which is too high.