Natalie G. K.
Wong
a,
Conor D.
Rankine
b,
Cate S.
Anstöter
a and
Caroline E. H.
Dessent
*a
aDepartment of Chemistry, University of York, Heslington, YO10 5DD, UK. E-mail: caroline.dessent@york.ac.uk
bChemistry – School of Natural and Environmental Sciences, Newcastle University, Newcastle, upon Tyne, UK
First published on 27th June 2022
While common molecular anions show a strong propensity to undergo electron detachment upon UV excitation, this process often occurs in competition with molecular ion dissociation. The factors that affect the balance between these two major possible decay pathways have not been well understood to date. Laser photodissociation spectroscopy of the deprotonated forms of the UV filter molecules, Homosalate (HS) and Octyl Salicylate (OS), i.e. [HS − H]− and [OS − H]−, was used to acquire gas-phase UV absorption spectra for [HS − H]− and [OS − H]−via photodepletion from 3.0–5.8 eV. No photofragmentation (i.e. dissociation of the ionic molecular framework) was observed for either [HS − H]− and [OS − H]− following photoexcitation, revealing that electron loss entirely dominates the electronic decay pathways for these systems. High-level quantum chemical calculations were used to map out the excited states associated with [HS − H]− and [OS − H]−, revealing that the minimum-energy crossing points (MECPs) between the S1 and S0 states are located in elevated regions of the potential energy surface, making internal conversion unlikely. These results are consistent with our experimental observation that electron detachment out-competes hot ground state molecular fragmentation. More generally, our results reveal that the competition between molecular dissociation and electron detachment following anion photoexcitation can be determined by the magnitude of the energy gap between the excitation energy and the MECPs, rather than being a simple function of whether the excitation energy lies above the anion's vertical detachment energy.
In this work, we extend our previous LIMS studies of deprotonated UV filters to characterize the intrinsic UV absorption profile and photodegradation pathways of deprotonated homosalate (HS; Scheme 1) and octyl salicylate (OS; Scheme 1), i.e. [HS − H]− and [OS − H]−, respectively. Salicylates such as HS and OS are relatively weak UV absorbers but are extremely commonly found in commercially available sunscreens and cosmetic skin care products. They are generally employed in these products to augment SPF (Sun Protection Factor) and enhance the water resistance of formulations, and function as excellent solubilizers for crystalline chemical UV filters like oxybenzone and avobenzone.5,6
The experiments we perform are supported by quantum chemical calculations to provide direct physical insight into how the photofragment identity maps the non-radiative relaxation channels. We are particularly interested here in investigating the extent to which photodissociation competes with electron loss following UV excitation. Although electron detachment is frequently the dominant decay pathway following UVA/UVB photoexcitation for isolated (i.e. gaseous) anions, branching into ionic molecular photodissociation channels can also be significant.3,4,7 To provide two examples from the UV filter molecules we have studied, deprotonated oxybenzone decays predominantly by electron detachment following UV photoexcitation whereas deprotonated 4-benzophenone photodecays with considerably higher yields of molecular photofragments.3,7 With the aid of the quantum chemical calculations performed as part of this work, we aim to obtain a better understanding of the factors that control electron detachment versus molecular dissociation branching. This is potentially an important consideration in organic sunscreen molecule design, since an enhanced propensity for electron detachment is associated with production of reactive free electrons which could be harmful in a biological context.
HS and OS are ortho-disubstituted molecules that possess a spatial arrangement allowing the formation of an intramolecular hydrogen bond which can facilitate excited-state intramolecular proton transfer (ESIPT) upon photoexcitation.8,9 In salicylates, ESIPT develops between the hydroxyl hydrogen and the carbonyl oxygen, and has been observed, in the case of methyl salicylate, to be maintained in the gas and solution phase.10 As seen previously in the case of the sunscreen oxybenzone,3,11 ESIPT is understood to facilitate the fast and efficient non-radiative decay to the electronic ground state, with excess UV energy being dissipated harmlessly as heat.12–17 This intramolecular hydrogen bond effectively lowers the excited-state energy, and accounts for its UV-B absorbance.8,18 Deprotonation of a salicylate molecule at the hydroxyl hydrogen therefore has the potential to strongly perturb its ability to undergo ESIPT, and in turn, lower their capacity to absorb in the UV-B, and perform effectively within a bulk sunscreen formulation.1 For that reason, it is important to examine the influence of protonation state on the intrinsic photochemical behaviors of HS and OS.
[HS − H]− and [OS − H]− were mass selected and isolated in the ion trap prior to UV laser irradiation. UV-Vis photons were produced by a 10 Hz Nd:YAG (Surelite™, Amplitude Laser Group, San Jose, CA, USA) pumped OPO laser (Horizon™, Amplitude Laser Group), giving ∼0.3 mJ across the range 400–216 nm (3.10–5.74 eV). Laser step sizes of 2 nm were used throughout. The laser beam was focused as described previously.19 Photofragmentation experiments were conducted with a set ion accumulation time of 10 ms and a corresponding fragmentation time of 100 ms, allowing for each mass-selected ion packet to interact with one laser pulse, and minimize the likelihood of multiphoton events. Note that photoproducts are not detected until the end of the 100 ms window, allowing us to identify photoproducts at the end of any decay processes. When fluorescence is negligible, the UV-excited gaseous ion will fragment upon excited-state relaxation, yielding an action absorption spectrum by photodepletion.19–21
(1a) |
(1b) |
As shown in eqn (1a) and (1b), IntOFF and IntON represent the parent ion intensities with laser off and on, respectively; IntFRAG is the photofragment intensity with the laser on; λ is the excitation wavelength (nm); and P is the tunable laser pulse energy (mJ). The photodepletion intensities were taken from an average of three repeat runs at each wavelength within the range studied. We note that photofragment ions with m/z < 50 are not detectable within our mass spectrometer since low masses fall outside of the mass window of the ion trap.
Minimum-energy crossing points (MECP) between the first electronically-excited singlet (S1) and ground electronic (S0) states of [HS − H]− and [OS − H]− were located via a home-built external objective function (Fij) optimizer following the approach of Martinez et al.;38 convergence was indicated when the change in the value of the objective function, ΔFij, was below 1.0 × 10−6 a.u., provided that the convergence criterion of 3.0 × 10−4 for the RMS gradient was simultaneously satisfied. Initial values of σ and α were set to 3.50 and 0.02, respectively. The quality of reference wavefunction was assessed via D1 diagnostic inspection and reporting throughout. Ref. 39 and 40 provide further background on the importance of MECPs in the mechanism of internal conversion.39,40
Fig. 1 Negative ion electrospray ionization mass spectra of (a) [HS − H]− (m/z 261) and (b) [OS − H]− (m/z 249), respectively. |
The photodepletion (gas-phase absorption) spectrum of mass-selected [HS − H]− across the 3.10–5.74 eV (400–216 nm) range is shown in Fig. 2a. Mass selection is a key advantage of the experimental approach we employ here as it allows us to directly probe the intrinsic properties of the deprotonated systems. The gaseous absorption spectrum of [HS − H]− displays strong absorption in the UV-A region (with an onset below the lowest measured photon energy of 3.1 eV) through a band I which then reduces in intensity through the higher UV-B region (3.6–4.2 eV). Absorption again increases between 4.2–4.5 eV (II), and subsequently plateaus through the high UV-C wavelength region (III). Band IV therefore marks the onset of an intense broad photodepletion region which extends into the lower UV-C.
Fig. 2 Gas-phase UV absorption (photodepletion) spectra of (a) [HS − H]− and (b) [OS − H]−, respectively. The solid line is a five-point adjacent average of the data points. |
The photodepletion spectrum of [OS − H]− (m/z 249), over the range 3.10–5.74 eV, is displayed in Fig. 2b. The gas-phase absorbance profile can be seen to be extremely similar to that of the spectra of [HS − H]− (Fig. 2a), reflecting the similarity of the chromophores for both molecules.
An important aspect of our gas-phase LIMS experiment is that it allows us to identify the ionic photofragments produced following UV excitation of an isolated precursor ion.3,4,7,19 However, in this case, no photofragments were observed for either [HS − H]− or [OS − H]− across the entire UV region scanned. (We note that only fragments with approximately m/z > 50 can be detected in our ion-trap mass spectrometer.) This situation is highly unusual, with [HS − H]− or [OS − H]− representing the first two molecules we have studied by LIMS which do not produce ionic photofragments following photoexcitation.1,3,4,41 To further explore the possible ionic breakdown pathways for these ions, we proceeded to study them using collision-induced dissociation (CID).
The fact that stable fragments can be observed from [HS − H]− and [OS − H]− following CID illustrates that these systems can undergo molecular dissociation with production of stable ionic fragments, i.e. electron detachment is not the only potential breakdown pathway available to these anions. We note that in other organic sunscreen molecules we have studied,3,4,7,41 the ions associated with ground electronic state fragmentation (i.e. the CID fragments) have been observed following photoexcitation. This behavior has been traced to rapid internal conversion occurring on the excited state surface, followed by molecular dissociation of the system on the hot electronic ground state.7 To better understand why this does not occur for [HS − H]− and [OS − H]−, we turn to quantum chemical calculations.
The C1-symmetric S0 minimum-energy geometries of [HS − H]− and [OS − H]− were located at the [RI-MP2/ADC(2)]/aug-cc-pVDZ level (Tables S1 and S2, ESI†), and key parameters relating to electron detachment [vertical detachment energies (VDEs) and vertical dipole moments (VDMs)] at these geometries are tabulated in Table 2. The C1-symmetric D0 minimum-energy geometries of the neutral counterparts, [HS − H]˙ and [OS − H]˙, were also located at the [RI-MP2/ADC(2)]/aug-cc-pVDZ level (Tables S3 and S4, ESI†), and these have been used to determine supplementary adiabatic electron detachment energies (ADEs) which are, additionally, tabulated in Table 2. The greater difference between the VDE and ADE in [OS − H]− (ca. 0.4 eV) compared to [HS − H]− (<0.1 eV) is reflective of the greater geometric reorganization in the former (a consequence of the greater flexibility of the alkyl ester group to respond to the charge state) compared to the latter; [OS − H]− and [OS − H]˙ are separated in space by 18.3 Å Da−1/2, while [HS − H]− and [HS − H]˙ are only separated in space by 2.6 Å Da−1/2.
System | VDEa/eV | ADEb/eV | VDMc/D |
---|---|---|---|
a VDE = Eneutral – Eanion (at the optimized anion geometry). b ADE = Eneutral – Eanion (at the optimized neutral and anion geometries, respectively). c VDM = μneutral (at the optimized anion geometry). Vertical dipole moments are important quantities in assessing whether a dipole-bound excited state is possible for these anions. See ref. 43 for more details.42 | |||
[HS − H]− | 3.90 | 3.88 | 2.33 |
[OS − H]− | 4.34 | 3.90 | 5.97 |
Bands II–IV in [HS − H]− and [OS − H]− are embedded in the electron detachment continuum (i.e., they are located above the VDEs), and it is reasonable to assume that electronic excitation at these bands produces excited states that decay primarily via electron detachment, i.e., without photofragmentation, while Band I is located below (yet in proximity to) the electron detachment threshold; consequently, other decay mechanisms, e.g., S0 ← S1 internal conversion (IC) or fluorescence/phosphorescence, cannot be ruled out in such a straightforward sense. However, the calculated gas-phase absorption spectra of [HS − H]− and [OS − H]− are in good agreement with the photodepletion spectra (Fig. 4), which demonstrates that luminescence processes are not competitive for either of these systems. (If they were, there would be a mismatch between the absorption spectrum and photodepletion since luminescence decay is not evident as photodepletion of the precursor ion.)20
A naïve search for S0 ← S1 MECP in [HS − H]− and [OS − H]−, starting from the Franck–Condon point, reveals a false S0 ← S1 IC channel along the C2O10 (carbonyl) stretching coordinate; recent work by Marsili, Prlj, and Curchod43 has demonstrated these artificial S0 ← S1 MECP to be the product of systematic weaknesses of RI-MP2/ADC(2) for this class of systems. The real S0 ← S1 MECP in [HS − H]− and [OS − H]− are quasi-ethylenic, i.e., they are associated with twisting about a bond and pyramidalisation of one of the bond termini. These S0 ← S1 MECP were located at the [RI-MP2/ADC(2)]/aug-cc-pVDZ level (Tables S5 and S6, ESI†); their geometries are illustrated in Fig. 5, and key geometric parameters are tabulated in Table 3.
Fig. 5 Illustrations of the geometries of the S0 ← S1 MECP of (a) [HS − H]− and (b) [OS − H]−. Atom labels for both molecules are given in Fig. S1 of the ESI.† |
System | ΔEa/eV | d /Å Da−1/2 | ϕC2–C1–C7–O8/° | ϕC2–C1–C7–O9/° | θ pyr.b | θ w.c |
---|---|---|---|---|---|---|
a Relative to the S1-state Franck–Condon point. b θ pyr. is a pyramidalisation index quantifying the extent of the pyramidalisation at the C7 site; it is defined here as 1.0 − ((α + β + γ)/360.0), where α, β, and γ are the angles around the C7 site. c θ w. is an index quantifying the out-of-plane displacement of O10; it is defined here as the angle between the C2–O10 vector and the plane containing the aromatic ring (C1–C2–C3–C4–C5–C6). | ||||||
[HS − H]− | +0.61 | 16.6 | 88.2 | 32.6 | 0.091 | 8.2 |
[OS − H]− | +0.94 | 21.9 | 82.8 | 26.0 | 0.098 | 8.6 |
The S0 ← S1 MECP in [HS − H]− and [OS − H]− are located 0.61 and 0.94 eV above the S1-state Franck–Condon point (Table 3), respectively, and 0.23 and 0.57 eV above the peak of Band I, i.e., they are not readily accessible and are, to an extent, barrier-activated. To map the S0 ← S1 IC channel, potential energy surfaces have been constructed between the Franck–Condon points and S0 ← S1 MECP of [HS − H]− and [OS − H]−via linear interpolation of internal coordinates (LIIC). Independent single-point energy calculations have been carried out at each one of 15 interpolated geometries, with the calculated potential energy surfaces for [HS − H]− and [OS − H]− presented in Fig. 6a and b, respectively.
Although barrier heights determined via LIIC represent an upper bound to the true barrier height at the applied level of theory [i.e., a reduction in the barrier height might be expected following a minimum-energy pathway (MEP) scan, for example], it is nonetheless apparent that these S0 ← S1 MECP are located in elevated regions of the potential energy surface and are not likely to be accessed readily without a surplus of excitation energy that would place the excitation into the electron detachment continuum. This is qualitatively consistent with the recent time-resolved laser spectroscopic study on HS carried out by Holt and Stavros et al.;44 the authors recorded decays associated with two time constants – one ultrafast process corresponding to excited-state intramolecular proton transfer (ESIPT; not observed here in the deprotonated system) and a much longer process corresponding to excited-state decay.45 The time constant for the latter process was found to have an associated energy dependence, consistent with a barrier being presented to S0 ← S1 IC.45,46 The authors did not locate the S0 ← S1 MECP for HS, but we expect it closely resembles the one located here for [HS − H]−.
The consequence of these calculations is that excited states formed via excitation at Band I are unlikely to decay via the identified S0 ← S1 IC channel; the generation of molecular fragments on the hot ground state (observed via CID; Fig. 3) is not observed post-photoexcitation in our LIMS experiments because the systems do not return non-radiatively to the hot ground state via IC. This means that the only channels available to these ions following photoexcitation at Band I are either fluorescence/phosphorescence, or electron detachment. As discussed earlier in this section, the excellent agreement between the calculation absorption spectrum (Fig. 4) and the experimental photodepletion spectra means that luminescence processes are not significant, allowing us to conclude that electron detachment is dominant.
Finally, we note that the S1 states of both [HS − H]− and [OS − H]− occur at energies below the respective calculated ADEs, so that electron detachment would not be predicted to be an energetically decay pathway. We have, however unambiguously observed such behavior in several other systems, with electron detachment being observed experimentally significantly below the expected electron detachment energies.47,48 This is likely due to the considerable residual internal energy of the anions in our experiment, which are produced at ambient temperatures through the electrospray ionization source.
Our results for [HS − H]− contrast to those observed by Holt et al. for neutral HS, where time-resolved ion-yield spectroscopy showed that UV photoexcited gaseous HS undergoes ultrafast decay with a time constant of >170 fs as the dominant excited-state decay pathway.44 Calculations performed to support these measurements suggest that photoexcitation is followed by excited state intramolecular proton transfer (ESIPT), where the enol form of HS is converted to the keto tautomer, followed by rapid non-radiative dissipation of the excess energy. This represents prototypical behavior for a well-behaved organic sunscreen molecule, involving the conversion of UV photons to thermal energy. Our results for [HS − H]− indicate that the ultrafast decay pathway available to neutral HS is lost upon deprotonation, due to removal of the most acidic proton occurring at the enol group. Hence, ESIPT can no longer occur. This evolution from a functioning sunscreen molecule to one which displays poor sunscreen action upon deprotonation mirrors behavior we have seen previously in oxybenzone,3 another common sunscreen molecule.
The pathways available following photoexcitation of an anionic molecule or cluster are complex.54,55 Below the electron detachment energy, photoexcitation can only access valence excited states, which can decay radiatively or through dissociation into an anionic and neutral fragment pair. Above the detachment energy, electronic excited states are coupled to the electron detachment continuum,56 and can thus be considered autodetaching resonances. Electron detachment can either occur directly from the excited state or following excited state decay back to the ground electronic state and ensuing thermionic emission.57–60 Similar pathways are possible for molecular photodissociation,4,7 with formation of two neutral photofragments following electron detachment being an additional possibility.61,62 The high-level calculations performed as part of this work reveal that the competition between molecular dissociation and electron detachment following anion photoexcitation of the deprotonated UV filters [HS − H]− and [OS − H]− is determined by the magnitude of the energy gap between the excitation energy and the MECPs, rather than being a simple function of whether the excitation energy lies above the anion's vertical detachment energy. Our current understanding is sparse of how the potential energy surfaces will change in the oily, condensed-phase environment of a sunscreen mixture.63 Nonetheless, the issue is potentially important given the harmful nature of free electrons in a biological context, and merits further investigation for the systems studied here.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2cp01612e |
This journal is © the Owner Societies 2022 |