Kate
Robertson
a,
William G.
Fortune
a,
Julia A.
Davies
a,
Anton N.
Boichenko
b,
Michael S.
Scholz
a,
Omri
Tau
a,
Anastasia V.
Bochenkova
b and
Helen H.
Fielding
*a
aDepartment of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK. E-mail: h.h.fielding@ucl.ac.uk
bDepartment of Chemistry, Lomonosov Moscow State University, 119991 Moscow, Russia
First published on 20th February 2023
Phenolate photooxidation is integral to a range of biological processes, yet the mechanism of electron ejection has been disputed. Here, we combine femtosecond transient absorption spectroscopy, liquid-microjet photoelectron spectroscopy and high-level quantum chemistry calculations to investigate the photooxidation dynamics of aqueous phenolate following excitation at a range of wavelengths, from the onset of the S0–S1 absorption band to the peak of the S0–S2 band. We find that for λ ≥ 266 nm, electron ejection occurs from the S1 state into the continuum associated with the contact pair in which the PhO˙ radical is in its ground electronic state. In contrast, we find that for λ ≤ 257 nm, electron ejection also occurs into continua associated with contact pairs containing electronically excited PhO˙ radicals and that these contact pairs have faster recombination times than those containing PhO˙ radicals in their ground electronic state.
An early picosecond transient absorption spectroscopy (TAS) study found that electron ejection was essentially instantaneous following photoexcitation at 265 nm.11 The first detailed femtosecond TAS study, by Chen et al.,12 found that 266 nm photoexcitation, high in the first absorption band (Fig. 1), generated solvated electrons, e(aq)−, on a timescale that was longer than typical charge-transfer-to-solvent (CTTS) processes in inorganic anions.13 A detailed global target analysis identified two competing pathways for e(aq)− formation: fast (sub-picosecond) electron ejection from vibrationally hot S1(1ππ*), hereafter termed simply S1, and vibrational relaxation followed by slow (few picosecond) electron ejection from vibrationally cold S1. A more recent femtosecond TAS study by Tyson and Verlet14 found that 257 nm photoexcitation, between the first and second absorption bands (Fig. 1), generated e(aq)− on a fast, sub-picosecond, timescale that was comparable with typical charge-transfer-to-solvent (CTTS) processes in inorganic anions. They concluded that the fast electron ejection was from vibrationally hot S1 and, to explain the difference between their observations and those of Chen et al.,12 proposed a Marcus model in which e(aq)− were formed via a tunnelling mechanism. We were curious as to why electron detachment from higher lying S1n(1nπ*) and S2(2ππ*) states, hereafter termed simply S1n and S2, had not been considered. This led us to undertake a systematic femtosecond TAS study of the photooxidation dynamics of aqueous phenolate following photoexcitation at a range of wavelengths, from the onset of the S0–S1 absorption band to the peak of the S0–S2 band, complemented by liquid-microjet photoelectron spectroscopy (LJ-PES) measurements and high-level quantum chemistry calculations.
The TA spectra were recorded using a commercial transient absorption spectrometer (Ultrafast Systems Helios Fire) following photoexcitation at seven pump wavelengths over the range 300–235 nm (4.13–5.28 eV). The UV pump beam was generated from one OPA and was attenuated to achieve a pulse energy of around 100 nJ at the sample. The broadband probe beam was created via white-light generation by focussing the 800 nm fundamental beam into a sapphire (visible) or calcium fluoride (UV) plate, giving a combined probe range spanning 350–720 nm. The relative polarisations of the pump and probe beam were set at the magic angle 54.7° to suppress polarisation effects. Solutions of 20 mM aqueous phenolate were prepared from aqueous phenol with concentrated sodium hydroxide solution added to achieve pH 13. The phenolate sample was flowed continuously at 10 mL min−1 through a Harrick cell with a 100 μm path length using a liquid diaphragm pump (KNF, Simdos 02). The instrument response function was determined at each pump wavelength using a fit to solvent-only spectra (Fig. S1†).
Photoelectron spectra of aqueous phenolate were recorded using our liquid microjet photoelectron spectrometer that has been described in detail elsewhere.15 Solutions of aqueous phenolate were prepared by addition of 0.1 mM phenol to ∼2 mM NaOH solution; all solutions were made up with pure water (resistivity >15 MΩ cm−1). The NaOH concentration was adjusted slightly to flatten the potential in the interaction region of the spectrometer,16 with complete formation of the phenolate anion verified by UV-visible absorption spectroscopy. Aqueous phenolate solution was passed into the spectrometer by a high-performance liquid-chromatography pump, operating with backing pressures ∼70 bar. The resulting microjet had a nominal diameter of 20 μm at a flow rate of 0.7 mL min−1. The solution was intersected with femtosecond pulses of ultraviolet light around 2 mm downstream from the liquid microjet nozzle and the resulting photoelectrons were guided by a strong inhomogeneous magnetic field into the time-of-flight region. The interaction region was around 1 mm from the skimmer, and the polarisation of the laser light was parallel to the time-of-flight axis. Photoelectron spectra of NO and Xe were recorded to measure the instrument function, potential gradient in the interaction region, and the vacuum level offset between the interaction region and the analyser.1,17
The VEEs of aqueous phenolate were calculated using the extended multi-configuration quasi-degenerate perturbation theory, XMCQDPT2,26,27 combined with the EFP method. The zeroth-order wave functions were constructed using the state-averaged complete active space self-consistent field method, SA-CASSCF. The active space included all valence π and π* orbitals as well as the n-orbital located at the oxygen atom. The XMCQDPT2 calculations employed a one-particle DFT/PBE0-based Fock matrix to obtain energies of all CASSCF semi-canonical orbitals.
The optimized structure of PhO− + 2H2O inside the water sphere of 1043 EFPs was also used in large-scale hybrid calculations of the first vertical detachment energy (D0 VDE) of aqueous phenolate. A series of structures of increasing size with additional outer-shell water molecules was constructed to test convergence of the calculated VDE with respect to the system size. The fully optimized QM/EFP system was additionally solvated with a large box of water molecules with a length of 100 Å. MD simulations with periodic boundary conditions were performed with the QM/EFP core kept frozen. The same MD protocol was used as in the initial MD step. Following equilibration, the system was slowly cooled down and then optimized. The final structure was cut into smaller-sized systems, with distances from water molecules to the chromophore ranging from 10 Å to 40 Å. The VDEs were calculated at the PBE0/(aug)-cc-pVDZ//EFP level of theory as a difference in energy between the chromophore anion and its neutral radical. Both the inner and outer water shells were treated as EFPs in these calculations. The calculated VDE reached a converged value of 7.3 eV for a system with ∼11250 water molecules. This value was further validated by high-level XMCQDPT2/SA(10)-CASSCF(8,8)/(aug)-cc-pVDZ+//EFP calculations, with both detached and valence states of the anion being treated simultaneously using a state-averaging procedure. The computational details are described elsewhere.1 All XMCQDPT2/EFP calculations were performed omitting the fragment polarization term, thus only allowing the QM electron density to be polarized in the field of effective fragments. The solvent response was calculated separately in a fully self-consistent manner at the PBE0/EFP level of theory and used as a correction to the XMCQDPT2/EFP results.
The electronic structure theory calculations were performed within the Firefly computational package,27 which is partially based on the GAMESS (US) source code.28
Fig. 2 480–720 nm transient absorption spectra of 20 mM aqueous phenolate at displayed pump-probe delays following photoexcitation at 287 nm (top) and 235 nm (bottom). The broad absorption centered around 700 nm is attributed to solvated electrons and the absorption feature with a peak around 540 nm is attributed to the S1 state of phenolate.12 |
Fig. 3 shows normalised kinetic traces of 20 mM aqueous phenolate at 540 nm (the S1 ESA) and 700 nm (e(aq)− absorption) following photoexcitation at seven wavelengths, from the origin of the S0–S1 absorption band at 300 nm, to the peak of the S0–S2 band at 235 nm, including the wavelengths employed in the earlier femtosecond TAS studies (266 nm and 257 nm).12,14 Kinetic traces before normalisation are presented in Fig. S9.† Even though the spectral profiles associated with S1 and e(aq)− extend over the full probe range (350–720 nm), their relative contributions to the transient absorption spectra can be deconvoluted due to their substantially different spectral profiles.
The 700 nm kinetic traces show that e(aq)− formation is more rapid following photoexcitation of S2 (λ ≤ 248 nm) compared to photoexcitation of S1 (λ ≥ 266 nm). It is also clear that there are two distinct timescales for e(aq)− formation following photoexcitation of S1: an initial few picosecond rise followed by a slower approximately 20 ps rise matching the fluorescence lifetime of S1 recorded by Chen et al.12 Interestingly, the timescales for e(aq)− decay also appear to depend on the photoexcitation wavelength: following photoexcitation of S2, the e(aq)− absorbance begins to decay within a few picoseconds, whereas following photoexcitation of S1, the absorbance only begins to decay after around 10 ps.
The 540 nm kinetic traces are also dependent on photoexcitation wavelength. Following photoexcitation with λ ≤ 248 nm, the 540 nm kinetic traces are virtually identical to the 700 nm kinetic traces. In contrast, following photoexcitation of S1 (λ ≥ 266 nm), the 540 nm transient absorption appears effectively instantaneously, then after a few picoseconds it begins to resemble the 700 nm transient absorption.
For 257 nm photoexcitation, it is clear from the 700 nm kinetic trace that there is a slower decay component in the 2–20 ps region that is characteristic of the S1 relaxation dynamics, and that the long time dynamics (100–1000 ps) has a contribution from the slower decay characteristic of S1. The contribution of S1 is confirmed by comparing kinetic traces at probe wavelengths of 515 and 680 nm (Fig. S16†).
Time constants obtained from fits to the transient absorption spectra at the six pump wavelengths corresponding to photoexcitation of either S1 (λ ≥ 266 nm) or S2 (λ ≤ 248 nm) are listed in Table 1. Following photoexcitation of S1, the spectra were fit using a procedure based on that employed by Chen et al.12 Following photoexcitation of S2, the spectra were best fit with four exponentials (one rise and three decays). Both fitting procedures are described in the ESI (Fig. S10–S16†). Photoexcitation at 257 nm is more complex and will be discussed later.
The time constants show that following photoexcitation to S1 (λ ≥ 266 nm), solvated electrons are formed on both few-picosecond and ∼20 ps timescales ( and τS1→e−(aq), respectively). This is consistent with the earlier assignment of two timescales to fast electron ejection from vibrationally hot S1, and vibrational relaxation followed by slower electron ejection from vibrationally cold S1.12 Following photoexcitation at 300 nm, close to the bottom of the well, the slow electron emission occurs on a timescale of 21 ± 0.8 ps which is very similar to the fluorescence lifetime of S1 measured by Chen et al.12 (22 ± 2 ps). In contrast, following photoexcitation to S2 (λ ≤ 248 nm), solvated electrons are formed on a sub-picosecond timescale and there is no evidence of the slower electron formation associated with relaxation to vibrationally cold S1. A similar observation was reported by Tyson and Verlet following photoexcitation at 257 nm.14
Following photoexcitation to S1 with λ in the range 287–266 nm, the e(aq)− decays are best described by fitting with two exponential decays, τr1 and τr2. Although geminate recombination is a complex process, and fitting exponentials only estimates the timescales, the τr1 timescales obtained from our fits are of the same order of magnitude as the faster recombination timescale observed in the earlier 266 nm TAS measurement (50 ps).12 It is worth noting that the e(aq)− absorbance does not decay to zero within our experimental window of 1 ns. Interestingly, we found that following photoexcitation at 300 nm, τr2 > 1000 ps, which is longer than our experimental window and thus has a large associated error (ESI†). Following photoexcitation to S2 (248 and 235 nm), the e(aq)− decays are best described by fitting with three exponential decays, τr1, τr2 and τr3. The last of these has a timescale longer than that of our experimental window of 1 ns. τr1 and τr2 are noticeably shorter than the equivalent timescales observed following photoexcitation of S1. This suggests that the PhO˙ radical could be formed in a higher lying electronically excited state, where recombination would be expected to be faster since the corresponding electron would have less kinetic energy, following photoexcitation at 248 and 235 nm compared to photoexcitation at λ ≥ 266 nm.
It is interesting to note that the fast time constant, , decreases following S1 photoexcitation from 300–266 nm. The decrease is most marked between 275 and 266 nm, which could suggest that there is some contribution from the S2 state at 266 nm. However, the long time dynamics (100–1000 ps) in the 700 nm kinetic trace for 266 nm photoexcitation follows that of S1 and there is no noticeable contribution from S2. Therefore, we attribute this decrease in the fast time constant to increasing vibrational energy in S1. In the Marcus picture,14 the tunnelling rate from S1 to the contact pair is very sensitive to barrier height and width, so the marked decrease in fast time constant between 275 and 266 nm photoexcitation could be because 266 nm excites close to the top of the barrier.
At 257 nm, simultaneous photoexcitation of S1 and S2 makes it difficult to separate the contributions from the S1 and S2 photooxidation mechanisms. The presence of the S1 ESA (Fig. S16†) as well as the faster electron formation time (consistent with the 257 nm TAS measurements reported by Tyson and Verlet14) and faster electron recombination time, confirm that 257 nm photoexcitation populates both S1 and S2.
To determine which photodetachment continua are accessible from the different electronically excited states of phenolate accessed in our measurements, we recorded photoelectron spectra following 1 + 1 resonance-enhanced photodetachment of aqueous phenolate at 285 nm (via S1), 236 nm (via S2), and 266 and 257 nm (wavelengths employed in the earlier TAS studies12,14). The spectra are plotted as a function of electron kinetic energy (eKE) in Fig. 4. The measured spectra are slightly distorted as a result of inelastic scattering of the photoelectrons in liquid water,16,17,29,30 so the true photoelectron spectra were retrieved using the algorithm described in ref. 31 and are overlaid on the measured spectra, in red. The photoelectron spectra obtained following photoexcitation at 285 nm (close to the peak of the S0–S1 absorption band) and 266 nm, correspond to resonance-enhanced detachment via S1 and are best fit by single Gaussians. In contrast, the broader spectra obtained following photoexcitation at 236 nm (close to the peak of the S0–S2 absorption band) and 257 nm, are best fit by three Gaussians (Fig. S17 and Table S1†). The increased broadening and number of Gaussians required to fit the two shorter wavelength photoelectron spectra compared to the two longer wavelength photoelectron spectra, indicates that higher lying electronically excited states of phenolate are accessible at 257 nm as well as 236 nm, and that these may be coupled to solvated electron continua associated with higher lying electronic states of the PhO˙ radical, supporting our interpretation of the different timescales for geminate recombination observed in our femtosecond TAS measurements.
Fig. 4 Resonance-enhanced 1 + 1 photoelectron spectra (PES) of aqueous phenolate at 236, 257, 266 and 285 nm (black) together with retrieved true PES (solid red lines).31 The 285 nm spectrum is adapted according to CC-BY license from ref. 31. Copyright 2022 the authors. Published by American Chemical Society under a Creative Commons CC BY 4.0 license. |
To determine which resonance-enhanced photodetachment pathways are accessible energetically, we carried out high-level quantum chemistry calculations of the electronic structure of the singlet states of PhO− and doublet states of PhO˙ (Fig. 5). Using Koopmans' arguments, we are able to determine that S0 is most likely to detach to D0 and D1n, S1 to D0, S1n to D1n, and S2 to D0, D1 and D2. Thus, it is clear that the peaks which are modelled well with single Gaussian functions in the 285 and 266 nm photoelectron spectra correspond to S1–D0 detachment. The situation is more complicated for the 236 and 257 nm photoelectron spectra, as their lineshapes could only be satisfactorily reproduced using a linear combination of Gaussian functions. The high eKE tails of both spectra lie at similar two-photon binding energies close to the D0 binding energy, suggesting that all four spectra contain contributions from detachment to D0. The 236 nm photoelectron spectrum likely has contributions from resonance-enhanced photodetachment via S2 to D0, D1 and D2, although we cannot rule out ultrafast internal conversion to S1 followed by photodetachment to D1. The 257 nm photoelectron spectrum likely has contributions from excitation of both S1 and S2, as discussed above. A detailed analysis of the photoelectron spectra of aqueous phenolate at a wide range of photon energies will be reported in a future publication.
Having explored which photodetachment pathways are accessible from the different electronically excited states, and knowing that the binding energy of the solvated electron is between 3.7 and 3.8 eV,29,32 we can explore which photooxidation pathways might be accessible (Fig. 6). The S1 state is presumably coupled to the solvated electron continuum associated with the PhO˙ radical in the D0 state, and the S2 state is likely coupled with the continua associated with the PhO˙ radical in the D0 and D1 states. Although Tyson and Verlet suggested that it was not possible to access a higher lying detachment continuum following photoexcitation at 257 nm, based on our calculations and the second vertical detachment energy of phenolate reported in ref. 33, we find that it is possible to access the D1 continuum from S2 following 257 nm photoexcitation. Thus, having determined the possible photooxidation pathways, and observed that higher lying detachment continua may be accessible following 1 + 1 photodetachment at 257 and 236 nm, we conclude that the faster timescale observed for electron recombination following photoexcitation at λ ≤ 257 nm could be the result of the PhO˙ radical being formed in an excited electronic state.
Fig. 6 Schematic energy level diagram illustrating the different photooxidation pathways and timescales following photoexcitation of the S1, S1n and S2 electronically excited states of phenolate in aqueous solution. Numbers not in parentheses are λmax (in eV) from the UV-vis spectrum presented in Fig. 1 and VDEs obtained from UV (D0) and X-ray (D1n) liquid-jet photoelectron spectroscopy measurements.31,33 Numbers in parentheses are the calculated VEEs and VDEs presented in Fig. 5. Horizontal arrows represent coupling of S2, hot S1 and cold S1 to the different solvated electron continua (shaded blocks) with representative timescales (see text). |
To confirm this, as Chen et al. observed the formation of the PhO˙ radical in its D0 state between 330 and 400 nm,12 we recorded the 350–450 nm transient absorption spectrum of aqueous phenolate following photoexcitation at 235 nm to see whether we could find direct evidence of the formation of the PhO˙ radical in a higher lying state (Fig. 7). In addition to the structured band system between 350 and 400 nm reported by Chen et al., we observe an additional peak around 427 nm and a shoulder on the longer wavelength edge of the peak around 400 nm. These additional features decay on a subpicosecond timescale and are not observed in the 350–450 nm transient absorption spectrum of aqueous phenolate following photoexcitation at 287 nm (Fig. 7). Similar features were observed for aqueous phenol following excitation at 200 nm and the small peak at 427 nm was assigned to the phenol radical cation.34 Such an assignment would not be consistent with our measurements as they were carried out at pH 13. Instead, we attribute the weak transient absorption around 427 nm to absorption of an excited electronic state of the PhO˙ radical, slightly red-shifted from the absorption band attributed to the radical in its ground electronic state. This is supported by our XMCQDPT2/SA(15)-CASSCF(9,8)/(aug)-cc-pVDZ//EFP calculations which show that in this region, the most bright D1 → D3 vertical transition is at 411 nm with an oscillator strength of 0.06, red-shifted from the most bright D0 → D2 transition located at 386 nm with an oscillator strength of 0.1. This is consistent with the earlier observation of this feature, as photolysis of phenol at 200 nm could generate phenoxy radicals in the D1 state.34 The shoulder on the long wavelength side of the 400 nm feature may arise from vibrational structure associated with the 427 nm feature or to hot/combination bands of the radical ground-state absorption. These UV TAS spectra reinforce our conclusion that photoexcitation at shorter wavelengths generates PhO˙ radicals in higher lying states that undergo more rapid geminate recombination than the PhO˙ radical in its ground electronic state.
Understanding the mechanism of phenolate photooxidation following photoexcitation of higher lying excited electronic states is directly relevant to the photooxidation of protein chromophores since the higher lying electronically excited states of both GFP and PYP chromophores have been implicated in the formation of solvated electrons.1,35–37 To the best of our knowledge this work represents one of the first direct observations of photooxidation of an organic chromophore in aqueous solution yielding an oxidised species in an electronically-excited state. The subsequent relaxation dynamics of such electronically-excited oxidised species represent an enticing target for future studies, particularly since chromophores containing phenolate motifs are ubiquitous in biology.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3sc00016h |
This journal is © The Royal Society of Chemistry 2023 |