Mateo
Forjan
a,
Goran
Zgrablić
a,
Silvije
Vdović
*a,
Marina
Šekutor
b,
Nikola
Basarić
b,
Piotr
Kabacinski
c,
Maryam
Nazari Haghighi Pashaki
d,
Hans-Martin
Frey
d,
Andrea
Cannizzo
*d and
Giulio
Cerullo
c
aInstitute of Physics, Bijenička cesta 46, 10 000 Zagreb, Croatia. E-mail: silvije@ifs.hr
bDepartment of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, Bijenička cesta 54, 10 000 Zagreb, Croatia
cIFN-CNR, Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, I-20133 Milano, Italy
dInstitute of Applied Physics, University of Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland
First published on 25th January 2022
The ultrafast photochemical reaction of quinone methide (QM) formation from adamantylphenol was monitored in real time using femtosecond transient absorption spectroscopy and fluorescence upconversion in solution at room temperature. Experiments were complemented by theoretical studies simulating the reaction pathway and elucidating its mechanism. Excitation with sub-20 fs UV pulses and broadband probing revealed ultrafast formation of the long-lived QM intermediate directly in the ground state, occurring with a time constant of around 100 fs. UV-vis transient absorption data covering temporal dynamics from femtoseconds to hundreds of milliseconds revealed persistence of the absorption band assigned to QM and partially overlapped with other contributions tentatively assigned to triplet excited states of the adamantyl derivative and the phenoxyl radical that are clearly distinguished by their evolution on different time scales. Our data, together with the computations, provide evidence of a non-adiabatic photodehydration reaction, which leads to the formation of QM in the ground state via a conical intersection, circumventing the generation of a transient QM excited state.
The most commonly used photochemical method to generate QMs in biological systems is photodeamination of suitably substituted phenols.16 The mechanism of photodeamination has therefore been investigated by a variety of approaches, including nanosecond laser flash photolysis (LFP),17 computations,18 and femtosecond transient absorption (fs-TA) spectroscopy.19 The major discovery was that photodeamination takes place in an adiabatic ultrafast photochemical reaction on the S1 potential energy surface, delivering QMs in the excited state and without any detectable intermediate en route.19 However, one serious drawback of photodeamination reactions is the fact that QMs react efficiently with water molecules,20,21 resulting in hydration products and loss of the starting material.
Photodehydration reactions22,23 are therefore a better option, since the reaction of QMs with water regenerates the starting molecules and therefore does not act as a material loss channel. For this reason, photodehydration reactions have been studied on several substrates, extended by investigations of antiproliferative activity of photogenerated QMs24–27 that highlighted some anthracene derivatives as potent agents against cancer stem cells.28 The mechanism of the photodehydration reaction has been fully unraveled on the simple o-hydroxymethylphenol (1) by computations and fs-TA measurements. It was demonstrated that the corresponding QM 2 (Scheme 1) is generated in an ultrafast adiabatic reaction on the S1 potential energy surface, taking place in less than 1 ps and generating the QM intermediate in its excited S1 state.29 However, it is not a general mechanism for the photodehydration, since the formation of QMs may involve benzoxete derivatives as intermediates,23 or it may proceed via phenoxyl radicals that undergo a homolytic cleavage of the alcohol OH group.27
Here we present a joint experimental and computational investigation of the photodehydration mechanism of the adamantyl derivative 3 (Scheme 2). A previous photochemical investigation demonstrated that 3 undergoes dehydration more efficiently (ΦR = 0.51)30 when compared to 1 (ΦR = 0.23),22 which raises questions about the mechanism and the timescale for the formation of the corresponding QM 4. Investigation of the photoreactivity of adamantylphenol derivatives is particularly important since the corresponding QMs generally show better antiproliferative activity compared to non-substituted QMs.24–27 To this end we combine fs-TA and fs fluorescence upconversion (FlUC) measurements performed on 3 in CH3CN:
H2O (1
:
1) with time-dependent density functional theory (TD-DFT) computations which support the observations and fully elucidate the photochemical reaction mechanism. Our main finding is that 1 and 3 do not follow the same photochemical pathway. Unlike the parent molecule 1, photodehydration of 3 appears to be a non-adiabatic reaction that proceeds via a conical intersection (CI), delivering directly QM in its ground electronic state, 4(S0). Our results demonstrate that changes in the molecular structure not altering the actual chromophore may nevertheless result in a profound change of the photochemical reactivity.
The second, high temporal resolution broadband fs-TA spectrometer used in these experiments has been described in detail elsewhere.31 Briefly, the setup is based on an amplified Ti:sapphire laser (Coherent Libra). The pump pulses are generated by SHG of a visible non-collinear optical parametric amplifier (NOPA) in a 20 μm-thick BBO crystal, resulting in sub-20 fs pulses in the UV tunable between 260 and 290 nm. A 35 μJ fraction of the laser output is sent to a computer-controlled delay line, followed by the WLC generation stage. The energy for WLC generation is adjusted to ∼1 μJ. For WLC a small fraction of the fundamental laser output is focused onto a 3 mm-thick CaF2 plate. Optionally, the WLC can be driven by the SH of Ti:sapphire, generated in a 1 mm-thick type 1 BBO crystal (θ = 29°). To avoid optical damage, a two-axis translation stage (Zaber Technologies, Vancouver, Canada, model T-LS28M) is used to slowly scan the plate in a plane perpendicular to the propagation direction. The pump and the probe beams are separately and non-collinearly focused on the sample at approximately 5° by two UV enhanced aluminum spherical mirrors with focal lengths fpump = 250 mm and fprobe = 200 mm, corresponding to focal spot diameters of 120 and 80 μm, respectively. The fluence of the pump beam with energy of 113 nJ was therefore 1 mJ cm−2.
Parallel polarizations of pump and probe beams were used for early time fs-TA measurements. The transmitted probe, selected by an iris, is focused on the entrance slit of a spectrometer (SP2300 Princeton Instruments, Acton, MA, USA) equipped with a CCD detector, with the spectral resolution of 0.5 nm, and a read-out electronics (Entwicklungsbuero EB Stresing, Berlin, Germany) capable of single-shot measurement of the probe spectrum at the full 1 kHz laser repetition rate.
In both experimental setups the pump is modulated by a mechanical chopper at 500 Hz and the absorbance difference spectrum as a function of probe wavelength λpr and pump–probe delay t is calculated as:
Synthesis of 3 has been described previously.30 For the first measurements over longer timescales, adamantylphenol 3 was dissolved in CH3CN–H2O (1:
1) in the concentration c = 5.0 × 10−4 M and the absorbance at the 267 nm excitation wavelength measured in a cuvette with 1 cm path was 0.7, while for the second set of measurements where sub-20 fs NOPA pump pulse resulted in lower fluence in order to have a reasonable level of the fs-TA signal the concentration was increased to c = 5.0 × 10−3 M. In these and subsequent time resolved fluorescence and laser flash photolysis measurements acetonitrile used for the measurements was of HPLC purity. Additionally, mQ H2O (Millipore) was used. All measurements were performed at ambient temperature (22 °C).
Based on previously published spectra for QM 4, the PA band at 425 nm most probably results from a transition from the ground state (S0) of the intermediate 4.30 This assignment is also in agreement with the results of our TD-DFT simulations that reveal a corresponding peak in the photoproduct absorption spectrum at 368 nm (see ESI,† Fig. S25c and Table S1). Although it is common that the computed spectra often do not precisely match the experimental spectra (here 368 nm vs. 425 nm) due to the nature of the used TD-DFT method, the qualitative match is in our case sufficient to reliably assign the photoactive species.
At longer wavelengths, a second PA band peaking at wavelengths longer than 600 nm is observed. Based on previous results,27 we assign it to solvated electrons. Transients due to the solvated electrons decay faster and with two characteristic time constants. From the bi-exponential fitting of data for this wavelength range, we obtain an initial fast decay of around 35 ps, followed by a slower decay of around 1.5 ns (see ESI† Fig. S4 to S6). Comparison with results of other studies focused on photochemistry of solvated electrons in different solvents is not straightforward since in our case we have a mixture of two different solvents and a solute. However, it was reported previously that in pure acetonitrile excess electrons take on two distinct forms – solvated electrons and solvated molecular anions – that are in equilibrium with each other.40 It is known that the timescale for ground-state solvated electrons to convert into solvated dimer anions is around 80 ps.41 These two forms have distinct absorption spectra and relaxation dynamics which in our case could be further modified by the presence of water and solute. To clarify the origin of solvated electrons, we measured excitation intensity dependence on the transient absorption spectra for both the sample and pure solvent at fixed delay of 10 ps that show quadratic dependence in the case of pure solvent, as expected for two-photon ionization (see ESI† Fig. S12). Intensity dependence for sample is almost linear, suggesting that it stems from free electrons resulting from one photon radical cation formation with only minor contribution of two-photon ionization of solvent molecules (see ESI† Fig. S11). Absence of signature of radical cation formation in our fs-TA experiments could be explained by the overlap with stronger solvated electron PA signal.
The observation of exclusively positive TA signal in the whole spectral range is due to the fact that the ground state absorption of 3 lies below 300 nm (see ESI† Fig. S13) and suggests that stimulated emission (SE) from the S1 state of 3 is either overlapped with a stronger PA signal, or decays on an ultrafast, sub-picosecond timescale.
The absence of short-lived spectral features that would correspond to the initially populated Franck–Condon state of the reactant 3 or to the excited state S1 of the photoproduct QM prompted us to perform another set of measurements with a higher, sub-50 fs temporal resolution. Fig. 2 shows kinetic traces at selected wavelengths following photoexcitation with sub-20 fs pulses at 280 nm. Due to the short duration of the pump pulses, strong artifacts (two-photon absorption, cross phase modulation) could not be avoided that complicated the analysis of the early reaction dynamics. The analysis of kinetics at relevant wavelengths showed that, apart from a faster rise of spectral features due to the increased temporal resolution, no evidence of intermediates was found. Measurements with WLC probe driven by SH of Ti:sapphire at 400 nm, extending the probed spectral range deeper into UV, revealed no spectral signatures either of additional photoexcited states or SE, and showed only the ground state bleach signal from adamantylphenol 3 in the 265–285 nm spectral range (see ESI† Fig. S14). The lack of any clear evidence of a SE, despite a detected steady-state emission, agrees with a signal dominated by the PA contributions, likely due to non-emissive species, and urges on carrying out time resolved fluorescence measurements.
![]() | ||
Fig. 2 Kinetic traces for the selected wavelengths using sub-20 fs UV NOPA pump pulse revealing early photochemistry of 3 in CH3CN–H2O (1![]() ![]() |
Careful analysis of the TA dynamics at selected wavelengths convinced us that it is not possible to reveal any intermediate state with sub-picosecond dynamics. We therefore decided to perform detailed fitting of single wavelengths using the sum of two exponential functions convoluted with a Gaussian instrument response function (IRF), where the first exponential function describes a fast rise, while the second one is connected with a slow decay (fixed at 1 ns). This procedure, which was employed to analyse the results of both experiments, allowed the use of a minimum number of free parameters. Details on the fitting procedure are given in the ESI.†
Results of fits of kinetic traces at several relevant wavelengths for the sub-50 fs TA measurements are summarized in Table 1. We further verified rise time of the TA signal obtained with sub-20 fs pumping by fitting the signal integrated in the 400–500 nm spectral window, which corresponds to the PA band ascribed to the QM (see ESI† Fig. S10). This procedure is justified by the fact that our analysis indicates absence of intermediate species with rapid dynamics but it also reduces the noise, resulting in better fitting of the rise time. In addition, in this spectral range the contribution from solvated electrons is less pronounced than at longer wavelengths. The 100 fs rise time obtained for the integrated signal is in good agreement with the rise times from single wavelength kinetic profiles (Table 1).
Wavelength [nm] | Rise time [fs] |
---|---|
425 | 78 ± 4 |
500 | 82 ± 6 |
600 | 125 ± 10 |
To provide additional evidence for assignment of the PA band at 425 nm to the ground state of the QM 4 and to confirm the absence of any prior intermediates, we collected nanosecond transient absorption data covering broad time-window up to hundreds of milliseconds.
Transient absorption spectra of 3 in Ar-purged and O2-purged CH3CN–H2O (1:
1), shown in Fig. 3, revealed additional transient species absorbing in the 300–450 nm spectral range. On the nanosecond time-scale, in the Ar-purged solution a short-lived transient was detected absorbing at 300–420 nm that was fitted to unimolecular kinetics at 400 nm with the decay time τ = 750 ± 30 ns (to be compared to the previously reported value τ = 3 μs).26 The transient was quenched by O2; in the O2-purged solution the decay time is 30 ± 10 ns. Based on the quenching by O2, this transient species may tentatively be assigned to the triplet excited state of 3.
![]() | ||
Fig. 3 Nanosecond transient absorption spectra of 3 (c = 2.0 × 10−4 M) in Ar-purged (left) and O2-purged (right) CH3CN–H2O (1![]() ![]() |
On a long time-scale, two transients were detected, a shorter-lived one absorbs at 350–400 nm and decays at 400 nm with the lifetime τ = 260 ± 10 μs in Ar-purged solution and τ = 420 ± 20 μs in O2-purged solution (Fig. 4). Since it is not quenched by O2, and according to the comparison with the published spectrum,39 it is assigned to phenoxyl radical (previously not detected). A radical cation could be formed in an ionization process, which would also explain the presence of solvated electrons detected as a red positive absorption band in fs-TA. Deprotonation of the radical cation then leads to the phenoxyl radical formation. However, since we found no evidence of radical formation in fs-TA measurements, it seems that more plausible interpretation is that the radical cation decays via back electron transfer while the phenoxyl radical could be formed directly in a H-transfer reaction that is known to evolve on nanosecond time scales.42 In this case, formation of phenoxyl radical cannot be resolved with fs-TA that covers only nanosecond time-window or in LFP measurements that use 10 ns excitation pulses.
![]() | ||
Fig. 4 Microsecond transient absorption spectra of 3 (c = 2.0 × 10−4 M) in Ar-purged (left) and O2-purged (right) CH3CN–H2O (1![]() ![]() |
The longer-lived transient in CH3CN–H2O has an absorption maximum at ≈420 nm, and decays to the baseline with the lifetime of τ = 600 ± 50 ms in Ar-purged solution and τ = 470 ± 20 ms in O2-purged solution (see ESI† Fig. S17). Based on previous reports, it is assigned to QM. Note that the transient is formed within the laser pulse (10 ns), since at longer wavelengths (500 nm) no growth of the transient absorption was detected on a short time-scale. More details on fitting results are given in the ESI.†
Results of nanosecond laser flash photolysis suggest that in fs-TA experiments at wavelengths below 425 nm transient absorption spectra probably contain contributions from the triplet excited state of 3 and the radical cation whose dynamics is too slow to be observed in fs-TA. However, the PA band at 425 nm, which remains visible in the entire time-window spanning from fs to hundreds of milliseconds, stems from the QM 4 in the ground state and not from a slowly depleting singlet excited state of 3, which should decay radiatively with nanosecond time constant.
We can give a tentative explanation for the absence of a negative signal in absorption transients. At a first glance, the apparent absence of SE of 3 seems to be in disagreement with its fluorescence quantum yield (ΦF = 0.13).30 However, it should be noted that the starting phenol derivative can adopt two distinct conformer orientations in the ground state that are energetically quite different (Fig. 5). Performed computations indicate that the conformer that forms an intramolecular bond between the two alcohol groups (3) is by 7.1 kcal mol−1 more stable than the conformer where such hydrogen bonding does not occur (3′). It is probable that the hydrogen-bonded conformer 3 is reactive and undergoes the QM formation reaction that quenches fluorescence on the 100-fs timescale. On the other hand, the less stable conformer 3′ is unreactive and therefore emissive, which would account for the observed fluorescence quantum yield. Note, however, that in an aqueous environment the conformer distribution can shift somewhat in favor of the less stable conformer since water molecules, with their hydrogen bonding network, can compete with the solute, influence the energetic stabilization and thereby increase the fluorescence efficiency, as was demonstrated previously.30 In addition, decay of fluorescence measured by TC-SPC for 3 in CH3CN–H2O (1:
1 v/v) was fitted to a sum of two exponents with decay times τ1 = 0.20 ± 0.03 ns and τ2 = 3.26 ± 0.01 ns. The detection of short and long decay times is also in agreement with the existence of both reactive and non-reactive conformers.
It is expected that the fluorescence of both conformers should be observed at wavelengths slightly higher compared to the absorption band of 3, around 300 nm. Indeed, this was confirmed by fluorescence upconversion measurements described in the next paragraph. The absence of SE signal in this wavelength range, that was covered by additional UV transient absorption measurements, can be explained by overlap of a weaker SE signal with a stronger PA signal.
To give further evidence of the rise component, we measured also the FlUC data for anisole in acetonitrile (see ESI,† Fig. S21), which is the non-photoreactive model for adamantylphenol 3. It is confirmed that 0.46 ps rise is not an instrumental artifact, since it is absent in the reference system and present only in the adamantylphenol 3. As we do not see any change in the shape, pure vibrational relaxation mechanism (energy redistribution and cooling) can be ruled out. Instead, a pure rise of a band without a change in shape would point to an indirect population mechanism, which probably stems from a conformational change, since for a simple phenol chromophore absorption and emission transition dipole moments should be co-linear. Based on the available information related to the photochemistry of this molecule, the most plausible explanation for the time-resolved rise of fluorescence is the rotation around the C–C single bond connecting the adamantane and the phenol, converting 3′ into 3 and vice versa, or the rotation of the phenolic OH group resulting in a different spatial approach of the two OH groups. Both species 3 and 3′ are present in the solution prior to the laser excitation. We assume that the excited reactive conformer 3 would undergo dehydration to QM 4 (100 fs, as indicated by fs-TA), and inefficiently convert to the emissive 3′, with the time constant of 460 fs. This difference in dynamics is also in agreement with the difference in quantum yields for dehydration (ΦR = 0.51) and fluorescence (ΦF = 0.13) from adamantylphenol.30 On the contrary, the emissive conformer 3′, which is higher in energy, may convert into 3 with 200 ps time constant, and that would account for the observed dual decay from 3′ (200 ps and 3.26 ns) detected by TC-SPC.
Moving to the interpretation of decay components, slow decay that was fixed to 2 ns in our analysis is ascribed to nonreactive emissive conformer 3′ in line with its fluorescence quantum yield (ΦF = 0.13).30 Even if we would have set the value from fluorescence decay measured with TC-SPC (3.26 ns), the only effect would have been on the relative amplitude of the two contributions but not on the values of the time constants. Interpretation of much faster decay that was clearly resolved within the available scan range of 50 ps is less trivial. Here, we cannot exclude that it reflects an anisotropy relaxation since the FLUC measurements could not have been performed at magic angle (see ESI† for more details). The fact that we found the same component in the reference system (see ESI,† Fig. S21), which does not undergo any conformational change, confirms the diffusional origin of the 22 ps component and rules out that it might describe any conformational relaxation. Details on IRF measurements, global analysis and FlUC data for the anisole reference molecule can be found in the ESI.† To summarize, temporal evolution of fluorescence can be explained by conversion between the two conformers 3′ and 3 appearing on different time scales without any need to invoke additional species.
Since the fs-TA experimental results indicate that 4(S0) is formed directly from 3(S1) without an intermediate structure, a CI pathway is the most probable one. Moreover, the failure to adiabatically optimize TS(S1) also indirectly points toward a presence of a CI since TD-DFT methods in general encounter problems when faced with CI geometries. However, a new implementation in the Orca 4.2.0 program38 enables simultaneous optimization of the ground and first excited state of the target structure in a way that facilitates approach toward the CI region, essentially confirming a reaction pathway towards a CI. In order to scan the potential energy surface of adamantylphenol 3 for such a path, we performed a CI optimization of the ground and the first excited state geometries and found the energy difference between them to amount to only 3.3 kcal mol−1, which is small enough to confirm the presence of the CI (see ESI,† Table S5).
Moreover, we indirectly resolved the passage of the wave packet formed at the Franck–Condon point in 3(S1) to the 4(S0) ground state via a conical intersection by fitting the kinetic traces at several chosen wavelengths where rise time of spectral features was longer (∼100 fs) than the time resolution of the measurement (∼50 fs). Femtosecond FlUC was employed with the goal to resolve temporal evolution of adamantylphenol 3 in its excited state revealing the existence of a nonreactive conformer 3′ with temporal evolution distinctive from the conformer involved in QM formation. Appearance of non-emissive QM could not be tracked with fluorescence measurements but lack of any changes related to photochemical reaction, also absent in transient absorption, indicates that ultrafast measurements are in line with direct formation of QM in the ground state. In addition, in the LFP experiment, transient spectra associated with absorption from QM in the ground state were tracked from nanoseconds up to hundreds of milliseconds, further supporting our interpretation of fs-TA results. These results, providing unifying picture on investigated photochemical reaction of QM formation, are summarized in Scheme 3.
The results shed new light on the photochemical reactivity of phenols. We have demonstrated that a change of the molecular structure, while not altering the chromophore, nevertheless changes the photochemical reaction pathway. This principle should be taken into account when designing molecules for particular reaction schemes.
Footnote |
† Electronic supplementary information (ESI) available: Additional information regarding the theoretical computations and analysis of experimental data. See DOI: 10.1039/d1cp05690e |
This journal is © the Owner Societies 2022 |