Solvent reorganization triggers photo-induced solvated electron generation in phenol

The analysis of the absorption spectrum and density of states of a cluster of phenol solvated with 15 water molecules indicates that the reorganization of the water molecules, facilitating the formation of solvated electrons, is a plausible mechanism in the photodissociation of phenol. Using quantitative wavefunction analysis, we demonstrate that while charge-transfer states involving electron transfer from phenol to water are mainly dark, a considerable number of them exists below the maximum of the pp* absorption band and could be populated by internal conversion. These low-lying charge-transfer states do not show extended O–H distances, but are found for large electron–hole separations at which several water molecules can solvate and stabilize the transferred electron. Thus, charge-transfer states in solvated phenol can be stabilized by two factors: (i) elongation of the O–H bond, as was extensively discussed in the past, and (ii) reorganization of solvent molecules, as it is shown here.

Based on theoretical calculations, [6][7][8][9] it has been proposed that after excitation to the bright first excited state of pp* character, a conical intersection can be reached upon O-H elongation to finally populate a charge-transfer (CT) state of ps* character with low oscillator strength.This state is dissociative with respect to the O-H bond, leading to the hydrogen atom and phenoxyl radical PhO products.For gas-phase phenol, the second excited state in the Franck-Condon region corresponds to an electronic excitation to a Rydberg-type orbital that, upon elongation of the O-H bond, takes on O-H-antibonding s* character.However, in aqueous solution this orbital is located in the solvent and, therefore, upon excitation of solvated phenol a concerted proton and electron transfer from phenol to the solvent has been postulated 6 to take place.
Experimental studies in the gas phase, [26][27][28][29][30][31][32][33] as well as in the apolar aprotic solvent cyclohexane, 46 support the mechanism of photodissociation by an O-H elongation, without electron transfer to the solvent (as the latter do not exist or no polar protic solvent is available).In aqueous solution, evidence for a photodissociation mechanism different from gas-phase phenol was recently found by Oliver et al. 47 using transient absorption spectroscopy.After excitation with 200 nm light, solvated electrons and phenoxyl radicals PhO were observed within 200 fs.Additionally, a signal attributed to phenol radical cations PhOH + , decaying on a sub-picosecond time scale, was measured.The occurrence of this intermediate species suggests that the photodissociation of phenol in water does not occur in a concerted way, but rather that the electron is first transferred to the solvent, generating PhOH + , and then the proton follows sequentially in less than a picosecond, generating PhO .Based on previous theoretical work, 6 one could speculate that these observations agree with a mechanism of direct population of the ps* state with transfer of the electron to the solvent, followed by O-H bond dissociation in the dissociative ps* state.In contrast, upon excitation of the bright pp* state at 267 nm the signals of solvated electrons and phenoxyl radicals begin to rise after about 2 ns. 47In this case, no signal corresponding to the intermediate PhOH + was observed, which the authors attributed to the short lifetime of this species. 47However, the absence of the signal may also indicate that the intermediate is not formed and that the proton and the electron are transferred to the solvent in a concerted way.Additional experiments on deuterated phenol (phenol-d 1 ) 47 at 267 nm revealed only a small kinetic isotope effect of 1.0 AE 0.4, thus indicating that proton transfer is not the rate-limiting step of the reaction in solution.However, this conclusion stands only if proton transfer is dominated by tunneling effects.Since for excitation at 267 nm no intermediate phenol radical cations PhOH + were observed, it is not clear whether the transfer of proton and electron to the solvent is concerted, as suggested by previous theoretical calculations, [6][7][8][9] or sequential, as observed after illumination at 200 nm.
Further evidence for a decoupled rather than concerted transfer of proton and electron to the solvent was found in the phenol-(NH 3 ) 5 cluster. 36While the near-infrared signal corresponding to Rydberg-type transitions of the transferred electron rose within 3 ps after excitation at 281 nm, the IR signal of the N-H vibrations indicating completed transfer of the hydrogen atom had a lifetime of 20 ps.These results show that first the electron is transferred to the solvent molecules and the proton follows shortly afterwards.
Recent photoelectron spectroscopy 35 applied to phenol in water, using excitation at 235.5 nm, which corresponds to an energy just above the conical intersection between the pp* and ps* states, indicated the formation of solvated electrons after 150 fs.They were attributed to a direct excitation to the ps* state leading to O-H dissociation by concerted proton-coupled electron transfer.However, these experiments could not rule out the occurrence of a mechanism involving a sequential transfer of proton and electron.
Actually, theoretical calculations by some of us predicted that sequential instead of concerted proton-coupled electron transfer operates in the phenol derivative 5,6-dihydroxyindole. 48fter excitation to the bright pp* state, the molecule relaxes to the minimum of this state.Then, the solvent molecules reorganize, facilitating the transfer of an electron to the solvent, generating a CT state.The proton transfer from one O-H moiety of the chromophore to the solvent occurs sequentially after the formation of the solvated electron.Thus, the reaction coordinate that drives electron transfer from the chromophore to the solvent is not the elongation of the O-H bond but solvent reorganization.
Motivated by these findings, the aim of the present paper is to investigate whether the mechanism of photodissociation by solvent reorganization, facilitating the formation of a solvated electron, as predicted for 5,6-dihydroxyindole, is also operative in phenol.If the reorganization of solvent molecules can lower the energy of the CT state without elongation of the O-H bond, one would expect that for some solvent configurations the CT state has a lower excitation energy than the bright pp* state already in the Franck-Condon region.Therefore, after excitation to the pp* state, the CT state could be easily populated by a conical intersection without O-H elongation.In order to investigate the existence of such low-lying CT states in the Franck-Condon region we report here the density of states (DOS) and the absorption spectrum of solvated phenol.In addition, we have characterized the involved electronic states by quantitative wavefunction analysis.We clearly show that some particular solvent configurations induce an energetic lowering of the CT state, supporting a sequential proton-coupled electron transfer mechanism in the excited-state dynamics of phenol.

Computational details
In order to describe a chromophore in solution, the introduction of vibrational sampling to sample the configurational space of both the chromophore and solvent is indispensable. 492][53] Phenol was solvated by a truncated octahedral box extended to a distance of 15 Å from any solute atom using the leap module of AmberTools14.The box was filled with 1223 TIP3P 54 water molecules while phenol was described by the general AMBER force field (GAFF). 55A steepest descent minimization and a conjugate gradient minimization have been carried out for 500 steps each.After minimization, the system was heated to 300 K at constant volume and temperature (NVT), using a Langevin thermostat 56 with a collision frequency of 5 ps À1 , for 200 ps with a timestep of 1 fs.Once the system was at 300 K, an MD simulation at constant pressure and temperature (NPT) was performed for 5 ns with a timestep of 2 fs.The Berendsen barostat 57 with a pressure relaxation time of 5 ps was used to keep the pressure at a value of 1 bar.The nonbonded interactions were truncated at 10 Å and the particle mesh Ewald method [58][59][60] was employed for the calculation of the Coulomb interactions using a grid spacing of 0.911 Å in each direction for the charge grid, in which the reciprocal sums are computed by a fourth-order interpolation, and a direct sum tolerance of 10 À5 .During the heating and the NPT MD simulations the bonds involving hydrogen atoms were constrained employing the SHAKE 61,62 algorithm.In order to allow for a more efficient sampling of relevant conformations where a water molecule forms a hydrogen bond with the acidic hydrogen atom of phenol, the water molecule closest to the acidic hydrogen atom of phenol was restrained during the heating and the 5 ns NPT MD trajectories employing a harmonic potential at an equilibrium distance between the acidic H atom of phenol and the O atom of water of 1.58 Å with a force constant of 50 kcal mol À1 Å À2 .Since an analysis of an unconstrained MD simulation shows that there is a water molecule hydrogen-bonded to phenol more than 80% of the time, the introduction of the constraint does not significantly adulterate the simulation.We consider that a hydrogen bond is formed when the separation between the oxygen atoms of phenol and water is smaller than 3.0 Å and the angle formed by these two oxygen atoms and the acidic hydrogen of phenol is larger than 1351.Note that in the unbiased simulation different water molecules form hydrogen bonds with phenol, while in the biased simulation the same water molecule interacts with phenol along the whole simulation time.
In order to calculate the absorption spectrum and DOS, 101 snapshots from the last 2 ns of the MD simulation were used where phenol is solvated by clusters containing a different number of water molecules.Specifically, three different clusters with 1, 5 and 15 water molecules were built.For the smallest cluster, only the water molecule forming a hydrogen bond with the acidic proton of phenol was selected along the 101 chosen snapshots, while the rest of solvent molecules were removed.For the larger cluster sizes, the four and fourteen water molecules closest to the hydrogen-bonded water molecule were selected, respectively.][70] The absorption spectra were obtained by convolution of the excitation energies of the 101 snapshots with Gaussian functions of width of 0.2 eV and height proportional to the oscillator strength.Since the energetic region of the pp* band is of interest to this study, only transitions with energies below 5.2 eV have been included in the absorption spectra and DOS.For comparability, the height of the spectra has been scaled to one.The same procedure was used to obtain the DOS but using Gaussian functions of height equal to unity.For comparison, also the absorption spectrum of phenol in the gas phase was computed following the same recipe.To this aim, the same snapshots were used and all the water molecules were removed.
Furthermore, in order to investigate whether bulk solvent effects, which are not properly described by the above cluster models, affect the excitation energies of phenol, two sets of QM/classical calculations including bulk-solvation effects were performed: one where the solvent is treated implicitly and another with explicit water molecules.In the implicit solvation calculation, we use the one-water-molecule cluster model described at the ADC(2)/aug-cc-pVDZ level of theory together with the conductor-like screening model 71 (COSMO) as classical model for aqueous solvation.In the second set of calculations, bulk-solvent effects were included by a electrostatic-embedding quantum mechanics/molecular mechanics [72][73][74][75] (QM/MM) scheme, where explicit TIP3P water molecules located in the MM region polarize the one-water-molecule cluster, which is defined as the QM region.The absorption spectra and DOS for both sets of QM/classical calculations were generated by using the same 101 snapshots and Gaussian convolution of excitation energies described above.
The excited states composing any of the so-obtained absorption spectra and DOS were characterized by analyzing the one-electron transition density with the TheoDORE package, 76-81 using the Lo ¨wdin population analysis. 82,835][66][67] The results are shown in Fig. S1 and S2 of the ESI.† As both the gasphase spectrum as well as the spectrum computed for the 5 water molecule cluster with CC2 are blue-shifted by 0.02 eV with respect to the ADC(2) spectra, and the convoluted spectra agree very well, all the results discussed below refer to ADC(2).

Results and discussion
We start by discussing the ability of the different cluster models in reproducing environmental effects.Fig. 1a shows the absorption spectra computed with phenol solvated by a different number of water molecules, as well as the experimental spectrum 35 in gas phase and in aqueous solution.In order to evaluate the solvatochromic effects, also the computed gas-phase spectrum is included.To ease the discussion, the absorption peaks of the calculated bands are collected in Table 1, together with the experimental values (see also Tables S1 and S2, ESI †).
The experimental spectrum in aqueous solution has a maximum at 4.59 eV 35,86 while the one in the gas phase possesses few vibronic transitions 35,85    spectrum in aqueous solution.Thus, already the inclusion of a single water molecule reproduces reasonably well the experimental spectrum.Note, however, that the computed spectrum is not able to describe the shoulder seen experimentally at 4.49 eV, which is due to vibronic transitions that are not included in our calculations.An increase from one to five water molecules causes only a slight red-shift of 0.01 eV.Interestingly, further increasing the cluster size to 15 water molecules leads to a noticeable blue-shift of 0.05 eV, i.e., bringing the spectrum further from the experimental solvated one by 0.2 eV.This energy behaviour could be an indication that the energy of the bright pp* state is not converged with respect to the number of water molecules treated explicitly.However, as the inclusion of more explicit water molecules at the quantum mechanical level of theory is intractable, we turn to hybrid QM/COSMO and QM/MM calculations in order to further simulate the effects of bulk solvation.The results obtained are collected in Fig. 1b.
In the QM/MM calculations, the QM region comprises phenol and the water molecule that is hydrogen-bonded to phenol -as we learnt from Fig. 1a that this water molecule induces an important red-shift in the spectrum and it is known 23,87 that hydrogen bonding plays a relevant role in the excited-state dynamics of phenol.The QM/COSMO spectrum, peaking at 4.78 eV, agrees very well with the spectrum calculated with 15 water molecules, while the QM/MM spectrum with a maximum at 4.85 eV shows worse agreement with the experimental band.In order to understand the origin of these differences, it is useful to remember that the QM/COSMO scheme treats the electrostatic interactions between the QM and classical regions in a selfconsistent manner, i.e., the chromophore and solvent are mutually polarized.In contrast, in the electrostatic-embedding QM/MM calculations the solvent is not polarized by the chromophore.Therefore, it is very likely that the erroneous blue-shift observed in the QM/MM spectrum is caused by the lack of polarizability in the electrostatic interactions.
In summary, from all the calculations, we observe that in terms of energetics the bulk solvent effects are better described by the QM/COSMO calculations than by the electrostatic-embedding QM/MM ones.Moreover, the agreement between the QM/COSMO spectrum and the spectrum calculated for 15 water molecules suggests that the energy of the pp* state is already converged with this cluster size.Within this solvation model, there is a small bathochromic effect of 0.07 eV, which is reasonable for the employed level of theory, taking into account that experimentally there is basically no solvent effect on the absorption energy.
Accordingly, all further analyses of the spectra were done using the model of phenol microsolvated by a cluster of 15 water molecules.One reason why we do not employ the QM/COSMO model is that in order to investigate later whether solvated electrons can be formed at the Franck-Condon region or not, an explicit QM description of several water molecules able to solvate the electron is necessary.
In order to find the most suitable model for the description of solvent effects, we also attempted calculations using different cluster models up to 10 water molecules combined with a COSMO solvation model.The inspection of the s* orbital involved in the CT state revealed that it is located at the surface of the water cluster, and the corresponding electron density is likely to extend beyond the cavity.As continuum models may fail [88][89][90] when the electronic density goes outside, we refrained from using COSMO for the 15 water molecules cluster model.
As a next step, the absorption spectrum and DOS of phenol solvated by 15 water molecules is decomposed into local and CT excitations.Local excitations in the energy range investigated here are pp* transitions where both p orbitals are located in the phenol ring, see one example in Fig. 2a with the two pairs of natural transition orbitals contributing to it.Instead, CT excitations involve an electron transfer from a p orbital of phenol to a s* orbital located in the solvent, as exemplified in Fig. 2b.][78][79][80][81] For this analysis, the system is divided into fragments. 91rbitrarily, we define phenol as one fragment and each water molecule as an additional fragment, i.e. we have 16 fragments.Two descriptors are used to classify the states.One is the CT number, 78 where g 0I ( r hole , r electron ) is the one-electron transition density matrix between the ground state and state I.The CT number gives the fraction of electron which is transferred from fragment A to fragment B. A pure CT transition corresponds to a CT number of 1, whereas a local excitation has a CT number of 0, as schematically illustrated in Fig. 2c.However, the CT number cannot unambiguously characterize different types of local excitations since any type of local excitation (located in phenol or in any water molecule) has a CT number close to zero.Thus, an additional descriptor is needed to distinguish between these two types of local states.
The second descriptor, denominated POS value, 77 provides the average position of the excitation in terms of the fragment number, i.e. the arithmetic mean of the position of the hole, POS hole , and the position of the electron, POS electron .The positions of the hole and of the electron can be calculated from the CT numbers and the norm of the exciton wavefunction O exc , 81 which is the sum of the CT numbers over all fragments, as Then, the mean position of the excitation is We label phenol (arbitrarily) as the fragment number 1 in our partition system while all the waters are numbered from 2 to 16, see Fig. 2c.Thus, a POS value close to 1 indicates that the excitation is located in phenol, i.e. both hole and electron are located within phenol orbitals.If, for example, for a given CT state the hole is located in phenol (fragment 1) and the electron is located in the water labeled as 3, then the POS value would be close to (1 + 3)/2 = 2.
Using the two descriptors, it is possible to define local pp* states as those having small CT numbers -we arbitrarily fix the value as r0.5 -and a POS value around 1, say r1.2.The CT states are then those having a CT value 40.5, regardless of the POS value.Note that with this definition, contributions from CT among the waters themselves are not excluded.However, by comparing the CT numbers with those from an additional analysis where only two fragments, phenol and the whole collective of water molecules, were considered it is possible to discern the origin of the CT transitions (see below).
Following the definition of the 16 fragments, the spectrum and DOS could be decomposed in terms of local pp* and CT states as shown in Fig. 3a and b, respectively.Panel a shows that the CT states are mainly dark, contributing only 5% to the spectrum.In addition, they are found at the blue end of the spectrum, mostly, at higher excitation energies than the maximum of the bright pp* band.Therefore, population of CT states by direct excitation is unlikely.However, despite being dark, there is a significant number of CT states in the energy range of the absorption spectrum, as it can be best appreciated in the DOS (Fig. 3b): 36% are CT states.The number of CT states if only two fragments (i.e.phenol and all water molecules) are considered is 34%, indicating that the majority of CT states come from excitations from phenol to the water.Interestingly, a considerable number of these CT states exists at energies below the peak of the pp* band.Indeed, the red tail observed in the DOS is virtually completely formed by CT states.These low-lying states can, in principle, be populated by internal conversion after excitation to the bright pp* state.This is an important result for the photodetachment mechanism of phenol, as the existence of these CT states indicates that the elongation of the O-H bond of phenol is not required in order to have low-lying CT states.Indeed, this can be corroborated by the computation of the average O-H distance of the geometries corresponding to these CT states that are energetically below the absorption band maximum.The computed distance is 0.9739 AE 0.0004 Å, which is very similar to the average O-H distance of all the geometries of the MD ensemble of 0.9740 AE 0.0004 Å, confirming that geometries with low-lying CT states do not possess elongated O-H bonds.
The present results suggest that the formation of solvated electrons by electron transfer from phenol to solvent without being accompanied by proton transfer in a concerted mechanism seems plausible.Accordingly, the sequential mechanism observed for dihydroxyindole, 48 where proton transfer occurs only once a solvated electron is formed, is also possible for solvated phenol.
Since the O-H elongation of phenol is not needed to lower the energy of CT states, it is very likely that such a lowering is caused by the arrangement of solvent molecules around the solvated electron.In order to test this hypothesis, we have This journal is © the Owner Societies 2019 attempted to optimize the crossing points between the pp* and CT states at the SA3-CASSCF(6,7) level 92,93 using the aug-cc-pvdz basis set.Geometry optimizations for the cluster of 15 water molecules are quite computationally demanding.Thus, the smaller cluster of 3 water molecules was employed.Despite the small number of water molecules present in the model, we were not able to identify any of the crossing points due to the bad convergence of the calculations as a consequence of the large number of low-frequency intermolecular normal modes.
Since the localization of crossing points was unsuccessful, we came back to the more realistic cluster model composed by phenol and 15 water molecules, and we performed a geometric analysis to identify the features of the solvent orientation that lead to low-energy CT states.We have calculated the root-meansquare (RMS) electron-hole separation of all CT states of the DOS, and plotted it versus the excitation energy of the states, see Fig. 4a.The data points were sorted according to the excitation energy and for each consecutive 20 data points the average excitation energy and RMS electron-hole separation were calculated.The trend clearly shows that larger electron-hole separations are associated with lower excitation energies.In particular, excitation energies below the energy at the maximum of the absorption spectrum (4.8 eV) possess an electron-hole separation of 5.5-6.0Å.This is a curious result since in order to form an exciton with a large electron-hole separation the attractive Coulomb electronhole interactions have to be broken.Therefore, there should be another factor that compensates the Coulomb penalty to lower the excitation energy of the CT states.This factor could be the solvation energy of the solvated electron since a large electron-hole separation indicates that the electron is well separated from phenol and, therefore, a large number of water molecules can surround and stabilize the electron.To investigate this possibility, the number of water molecules solvating the electron was examined.We consider that a water molecule is solvating the electron when it contains more than 10% of the excited electron population.Fig. 4b shows that a larger electron-hole separation allows for more water molecules surrounding the electron corroborating, thus, our hypothesis.Specifically, for an electron-hole separation of 5.5-6.0Å, which is the separation needed to have CT excitation energies lower than the energy at the maximum of the absorption band, 3 to 4 water molecules solvate the electron.This amount of water molecules provides a favorable solvation energy, which is able to compensate the unfavorable Coulomb penalty of having large exciton sizes.
Finally, it is worth to compare the experimental time scale of formation of solvated electrons after excitation to the pp* state with our stationary computations.Solvated electrons and phenoxyl radicals were experimentally detected after 2 ns. 47Therefore, contrary to what happens upon excitation to the ps* state, the photodetachment from phenol to water after excitation to the bright pp* state is a slow process.As seen in Fig. 3b, a considerable number of CT states lie at the low-energy region of the DOS band.At first glance, one could conclude that the population of these low-energy CT states from the initially excited pp* state might occur in an ultrafast manner (assuming a relatively strong coupling between both electronic states).However, the existence of low-energy CT states does not imply that for a particular phenol geometry and solvent arrangement the CT state is lower in energy than the bright pp* state, which would be the ideal energetic scenario leading to ultrafast internal conversion from the pp* state to the CT state, accompanied by formation of solvated electrons.Fig. 5 displays the probability distribution of the energy difference between the CT and pp* states for the selected classical MD ensemble of geometries.As can be seen, for the vast majority of geometries the energy of the pp* state is smaller than that of the CT state.In fact, the position of the maximum at 0.6 eV indicates that the two states  are often energetically well separated.This means that after excitation to the pp* state the CT state is likely not energetically accessible by internal conversion in an ultrafast manner.Population of the CT state will only happen after a strong red-shift of this state, which is induced by the diffusion of water molecules away from phenol such that the electron-hole separation increases, as shown in Fig. 4. The diffusion of solvent molecules required to induce solvated electron formation is a slow processa fact, which agrees with the nanosecond time scale experimentally determined by Oliver et al. 47 This speculative conclusion could be corroborated by running nonadiabatic dynamics simulations.However, the relatively large size of the system and the slowness of the process precludes performing such simulations with our current methods.

Conclusions
We show that the sequential proton-coupled electron transfer mechanism observed in dihydroxyindole 48 is also possible in phenol.This conclusion is drawn from the characterization of the absorption spectrum and DOS of phenol solvated with 15 water molecules.Our analysis shows that the CT states, which correspond to the formation of the precursor of a solvated electron by electron transfer from phenol to the solvent without accompanying proton transfer, are considerably abundant at energies below the energy at the maximum of the absorbing pp* band.Despite mostly dark, these low-lying CT states can be populated by internal conversion after excitation to the bright pp* state.As the O-H distances of the geometries with low-lying CT states are not elongated compared to the ensemble, the main factor contributing to the lowering of the energy of CT states in phenol is not elongation of the O-H bond.Instead, the lower excitation energy of the CT state is associated with a large electron-hole separation of about 5.5-6.0Å.While the larger exciton size results in a loss of attractive Coulomb interactions between hole and electron, it allows for a larger number of water molecules to surround and stabilize the electron, thereby reducing the energy of the CT state.We thus conclude that a sequential proton-coupled electron transfer mechanism is operative in phenol, complementing the well established photodissociation mechanism involving a stabilization of CT states by O-H elongation.Our results demonstrate that the reorganization of solvent molecules may drive photodissociation of phenol and suggest that the investigation of this mechanism in other phenol derivatives is worthwhile.
around 4.6 eV (270 mm), indicating that the solvent barely affects the energy-range of the spectra.The calculated spectra show small shifts, depending on the model used, which are within the error of the method.Inclusion of just one water molecule leads to a red-shift of 0.11 eV compared to the calculated gas-phase spectrum, and a blue-shift of 0.16 eV compared to the experimental

Fig. 1
Fig. 1 Absorption spectra calculated for phenol (a) solvated by a different number of water molecules, or (b) with different environment models, as indicated.Experimental data from ref. 35.

Fig. 2
Fig. 2 (a) Natural transition orbitals corresponding to the pp* excitation and (b) the charge-transfer (CT) state for one solvent geometry.(c) Illustration of the CT number and the POS descriptors, with an excitation within phenol (left) or from phenol to a water molecule (right).

Fig. 3
Fig. 3 Decomposition of (a) absorption spectrum and (b) density of states (DOS) (b) into pp* and charge-transfer (CT) states.

Fig. 4
Fig. 4 (a) Root-mean-square (RMS) electron-hole separation vs. average excitation energy of the charge-transfer (CT) states.Standard deviation indicated by the grey area.(b) Average electron-hole separation vs. number of water molecules with electron population 410%.

Fig. 5
Fig. 5 Probability distribution of the energy difference between the CT and pp* states for the 101 snapshots from the classical MD simulation for which electronic excited states were computed.

Table 1
Excitation energies at the maxima of the absorption peaks of experimental and calculated spectra of phenol in different environments; sh = shoulder Environment Cluster size Excitation energy [eV] Exp. a Gas phase 4.51, 4.63, 4.73 Exp.b Gas phase 4.51, 4.61, 4.63, 4.67, 4.73 Exp.b,c Aqueous solution 4.59, 4.49sh Calc.