Götz
Bucher‡
*a,
Mukul
Lal
b,
Anup
Rana
b and
Michael
Schmittel
*b
aLehrstuhl für Organische Chemie II, Ruhr-Universität Bochum, Universitätsstr. 150, 44801 Bochum, Germany
bFB 8 (Chemie-Biologie), Organische Chemie I, Universität Siegen, Adolf-Reichwein-Str. 2, 57068 Siegen, Germany. E-mail: schmittel@chemie.uni-siegen.de
First published on 18th July 2018
The photochemistry of two Barton esters, one derived from a dioxolane carboxylic acid and the other from pivalic acid, was investigated by product analysis and nanosecond laser flash photolysis (LFP). As expected, photolysis of the pivalate ester resulted in formation of the pyridine-2-thiyl and the t-butyl radical. Photolysis of the Barton ester of 2,2-dimethyl-1,3-dioxolane-4-carboxylic acid, on the other hand, revealed a complex multi-step fragmentation. In addition to the pyridine-2-thiyl and dioxolanyl radical, we gained evidence for the formation of the vinyloxy radical, CH2CHO˙. The latter was identified in the LFP by its π-complexes with benzene and diphenylether, its rapid quenching by electron-rich arenes and tri-n-butyl tin hydride, and its oxidative power in presence of trifluoroacetic acid as demonstrated by the oxidation of ferrocene to ferrocenium. Formation of CH2CHO˙ can be rationalized via fragmentation of the dioxolanyl radical. As the calculated barriers are too high for the reaction sequence to occur on the LFP time scale, we investigated the fragmentation of the photoexcited Barton ester via Born–Oppenheimer molecular dynamics simulations. In one trajectory, we could observe all reaction steps including ring opening of the dioxolanyl radical, suggesting that the excess energy gained in the ester cleavage and decarboxylation may lead to fragmentation of the hot dioxolanyl radical.
The present work was undertaken with the aim to establish an alternative photochemical route to enol radical cations. It is well known that acetals are acid-sensitive functionalities.14 We therefore reasoned that the acid-catalyzed ring opening of a radical derived from a cyclic acetonide should represent a potential route to obtain the radical cation of a free enol and/or its conjugated base, the vinyloxy radical. Our approach aimed at the generation of ethenol+˙ (1), the radical cation of the simplest enol. Explicitly, the intention was to generate 1 by acid-induced fragmentation of the dioxolanyl radical 3, which in turn should be readily available from decarboxylation of the carbonyloxy radical 4 (Scheme 1). The latter species was intended to be generated by photolysis of the N-hydroxypyridine-2-thione ester (“Barton-ester”) 7. As side products, bis-(2-pyridyl)disulfide (6) and acetone (2) are expected to be formed.
For the present study, we employed product analysis, UV-vis spectroscopic monitoring during steady-state irradiation, and laser flash photolysis (LFP) with UV-vis detection. In addition, we applied density functional theory (DFT) including Born–Oppenheimer molecular dynamic (BOMD) computations to elucidate the reaction pathways occurring in the photolysis. For comparison, we investigated the photochemical decomposition of the t-butyl derivative 8, which should not undergo fragmentation beyond formation of the t-butyl radical (10) (Scheme 2).
The UV-vis spectral changes in the acid addition experiments indicate that protonation of 7 predominantly occurs at the chromophore of the molecule. At low concentration of TFA, the most likely protonation should afford the cation 17 (Scheme 4).
For preparative purposes the photolyses were conducted at higher concentrations of 7 (1.57 mM). Product formation was studied by IR, 1H NMR, 13C NMR and HPLC (for quantification). Notably formation of product 6 was increased, most likely because at higher concentrations it arises from the reaction of radical 5 with 7.15 HPLC quantification showed that the overall conversion of 7 → 18 in the absence of TFA was 89% (by 1H NMR no further product was detectable) and 59% in the presence of acid (see ESI†). Formation of 6 was the same for both acid-free and acidic conditions.
A control experiment was conducted with Barton ester 8. Photolysis of 8 (568 μM) in the presence of an equimolar concentration of ferrocene and TFA showed a decrease of the absorbance at 370 nm and the formation of 6 and 19. Again, no formation of any long absorbing species above 400 nm was detectable. Photolysis of 8 (568 μM) in the presence of an equimolar concentration of ferrocene and TFA showed a decrease in the 370 nm band but no formation of any long absorbing species above 400 nm. Even with 9 equiv. of TFA (5.19 mM) no formation of any long absorbing species was observed (see ESI†).
The above experiments with 8 in presence of TFA and ferrocene suggest that the putative oxidants 6+˙ or protonated 5·H+ are either not formed during the photolytic conditions or have very short transient lifetimes insufficient for intermolecular electron transfer with ferrocene. Laser flash photolysis of 6 has been carried out by Osamu Ito et al. in various solvents,16 but without any evidence for 6+˙ or 5·H+. Thus, the finding of ferrocenium formation with 7, TFA, ferrocene and light suggests that another, yet unknown strong oxidant is formed in the photolysis that is a result of the dioxolane's follow-up reactions.
As Barton esters are known to be prone to decomposition in radical chain reactions, a plausible mechanism for the formation of 18 (or 19) appears to be the attack of radical 3 (or 10) onto the sulfur atom in 7 (or 8), which would eventually release another equivalent of radical 3 (or 10) to continue the radical chain process (vide infra). Disulfide 6 in turn could be formed by dimerization of radical 5 or by radical-induced decomposition of 7 (or 8),15 both mechanisms explaining increased formation of 6 at higher concentrations of 7 (or 8).
In summary, the observations indicate that in the presence of TFA and light an oxidizing species is formed in the solution from 7 that is derived from the dioxolane unit. While the data would be consistent with formation of ethenol radical cation (1) as hypothesized in Scheme 1, we note that the oxidation of ferrocene could also be due to reactions of other short-lived oxidants.
LFP of pivalate 8 (c = 0.1 mM) in benzene solution gave the transient spectrum shown in Fig. 2. We observed a single transient with λmax ∼ 480 and ca. 570 nm (second order decay, half-life in the μs range), which we assign to the pyridine 2-thiyl radical (5), complexed to a benzene molecule. Our spectrum agrees well with the spectrum reported by Aveline, Kochevar and Redmond,15 using benzene as solvent, whereas in THF solution, radical 5 has been reported to show an absorption maximum with λmax = 490 nm.16 π-Complex formation with benzene is predicted (TD-B3LYP/6-311++G(d,p)(pcm,benzene)//M06-2X/cc-pVTZ)(pcm,benzene) to result in a broadening of the spectrum with two absorption maxima at λmax = 425 and 556 nm, in reasonable agreement with the experimental observation (for the calculated absorption spectrum of a complex of 5 and benzene, see the ESI†). Adding 0.5 mM TFA did not change the transient spectrum. In agreement with the absorption maxima of the precursor, there is a bleach at λ = 295 and 375 nm. Interestingly, at both wavelengths we observed an apparent recovery of the signal, which occurred with a first-order rate constant k = 7.4 × 104 s−1, corresponding to τ = 13.5 μs. Usually, bleach and recovery like this is observed, if a long-lived triplet excited state returns to the ground state and the precursor concentration thus is replenished. This certainly cannot be the case here. We therefore suggest that the bleach and recovery reflects the reaction of the t-butyl radical (which we do not observe directly) with the precursor 8, yielding a species that absorbs at a similar place as 8, yet stronger.
LFP of acetonide 7 gave significantly weaker transient signals. A larger number of different transient species observed conveyed a much more complex mechanistic picture. Fig. 3 shows the transient spectrum observed upon LFP of a 0.1 mM solution of 7 in benzene, purged with argon. We observed a transient A decaying according to 2nd order kinetics (λmax = 445 and 485 nm) that is again assigned to a complex of thiyl radical 5 with benzene.16 Additionally, a transient B with λmax = 315 nm and a much shorter lifetime (1st order: τ = 170 ns) was observed (Fig. 3). If the laser power was reduced to ca. 20 mJ per pulse, we again observed growth or recovery at the absorption maxima of the precursor. This reaction took place with a lifetime τ = 170 ns. On a yet shorter time-scale, we detected another transient C, showing a λmax = 500 nm that appears to be formed in a monophotonic process (the plot of transient intensity vs. laser power is linear, see Fig. S6, ESI†). Its lifetime is extremely short (τ = 34 ns) (Fig. 4).
Addition of 0.03% TFA to the solution did not affect the transient spectra of A, B or C. The lifetimes of the three transients also did not change. We could, however, observe certain changes in the reactivity of transient C. This species turned out to be quenched by electron-rich aromatic quenchers such as anisole at a rate that was close to or at the diffusion-controlled limit (Fig. 5). Adding TFA further increased the reactivity slightly. Table 1 shows quenching rate constants measured both in the presence and in the absence of TFA.
Quencher | Rate constant in benzene | Rate constant in benzene + 0.03% TFA |
---|---|---|
1,2,3-Trimethoxybenzene | (8.2 ± 0.8) × 109 l mol−1 s−1 | (1.1 ± 0.1) × 1010 l mol−1 s−1 |
Anisole | (4.9 ± 0.5) × 109 l mol−1 s−1 | (8.9 ± 0.9) × 109 l mol−1 s−1 |
Diphenylether | Lifetime increases | (1.4 ± 0.1) × 1010 l mol−1 s−1 |
Toluene | No quenching observed | No quenching observed |
Tri-n-butylstannane | (1.6 ± 0.2) × 109 l mol−1 s−1 | Not measured |
Cyclohexane | No quenching observed | No quenching observed |
In the presence of TFA, the quenching rate constants turned out to be consistently slightly higher than in the absence of TFA. This effect was particularly obvious in the case of anisole as quencher. When diphenylether (DPE) was added, the transient lifetime measured at λ = 500 nm increased in the absence of TFA (τ = 90 ns in the presence of 0.08% diphenylether), and decreased in the presence of TFA. In the absence of TFA, we also observed a bathochromic shift in λmax (from λmax = 500 nm in the absence of DPE to a double maximum at ca. 500 and 630 nm in the presence of DPE), and the transient absorption spectrum became more intense and much broader (Fig. 6). Whether the growth (τ = 7 ns) that is seemingly resolved during the initial stages of the transient trace (Fig. 7, upper trace) is real or not is uncertain, as this growth lifetime is very close to the rise time of our LFP detector.17
Fig. 7 Transient traces, monitored at λ = 500 nm upon LFP (355 nm) of 7 in benzene. Bottom trace: Without diphenylether. Top trace: In the presence of 0.08% diphenylether. |
Considering the three different transient species observed, the assignment of A to thiyl radical 5 is straightforward based on a comparison with literature data,16 as well as a comparison with a calculated spectrum of the complex 5·benzene. The identity of transients B (λmax = 315 nm, τ = 170 ns) and C (λmax = 500 nm, τ = 34 ns), however, still needs to be discussed. In the experiments with pivalate 8, we observed a growth or recovery in the absorption bands of the precursor that occurred with a lifetime τ = 13.5 μs. As there is no transient observed decaying with this lifetime, it appears reasonable to attribute this growth to a reaction of the t-butyl radical, which is not expected to be detectable in our experiments owing to the lack of a good chromophore absorbing at λ > 300 nm.
While we can only speculate about the nature of this reaction, our possible candidate for transient B is the dioxolanyl radical 3. TD-DFT calculations (TD-B3LYP/6-31 + G*//B3LYP/6-31G*) reveal that this radical in the gas phase should have a longest-wavelength transition at λ = 302 nm, which should not be particularly intense (f = 0.0016). Thus there is qualitative agreement between the calculated absorption maximum of 3 and the absorption maximum of B. Given the fact, that formation of thioether 18 as a major product is strongly indicative of the intermediacy of 3, the assignment appears reasonable. Its rather short lifetime (τ = 170 ns) likely stems from a rapid reaction with the precursor. Given a diffusion-controlled reaction between 3 and 7, a sub-μs lifetime would be expected for 3 in the sub-millimolar concentration range employed for precursor 7.
Assigning the short-lived transient C is less straightforward. Possible candidates include a triplet excited state of precursor 7, carbonyloxy radical 4, the ring-opened radicals 20 or 21, and the vinyloxy radical 22 (Scheme 5 and Fig. 8).
Scheme 5 Possible intermediates in the photochemical fragmentation of Barton ester 7 in the absence of TFA. |
As far as the triplet-excited Barton ester 37 or carbonyloxy radical 4 are concerned, the barriers for fragmentation are expected to be very small. In fact, in both cases DFT optimizations resulted in either spontaneous N–O bond cleavage (37, M06-2X/cc-pVDZ) or in a very small barrier to decarboxylation (4, M06-2X/cc-pVQZ: ΔU‡ = 0.8 kcal mol−1). Born–Oppenheimer molecular dynamics simulations (M06-2X/cc-pVDZ) also in both cases resulted in C–C cleavage (4) or N–O cleavage followed by C–C cleavage (37), within a few femtoseconds each. Hence, an assignment of transient C to either species does not appear to represent a viable option. Among the remaining radical species 20–22, the vinyloxy radical (or formylmethyl radical) 22 represents the most likely assignment, as DFT calculations indicate that the barrier for decay of both 20 and 21 is smaller than the barrier for their formation, meaning that 20 or 21 should normally only be present in quasistationary concentration (Fig. 8). We note that the results obtained using the M06-2X/cc-pVQZ method are in excellent agreement with the results from the CCSD(T) single point energy calculations. Over all, the fragmentation of 3 is predicted to be endothermic, but exergonic (at T = 298 K), hence entropy-driven.
In the gas phase, the longest-wavelength absorption maximum of 22 is calculated (TD-B3LYP/6-31+G*) at λmax = 271 nm. However, if π-complexation to benzene or DPE is taken into account, additional weak long-wavelength CT bands appear at λmax = 468 nm (π-complex with benzene) or λmax = 464 and 689 nm (π-complex with DPE) (all values TD-UB3LYP/aug-cc-pVTZ//UM062X/cc-pVQZ). Hence, with these CT bands, the calculated UV-vis spectra of the vinyloxy/arene complexes qualitatively match the experimental transient spectra shown in Fig. 4 and 6, including the double maximum observed for the DPE complex. Fig. 9 shows a juxtaposition of the calculated vs. the experimental spectrum of 22·DPE, recorded in the presence of 0.08% DPE.
As far as the quenching reactions of the short-lived transient C are concerned, we note that both the high reactivity towards electron-rich arenes and the increase in reactivity upon addition of TFA would also be consistent with an assignment to the vinyloxy radical (22). A hydrogen-bridged complex of 22 with TFA would be very similar to a free ethenol radical cation in both structure and likely also oxidative power. Such complex would explain the pronounced reactivity of 22 in acidic medium, in particular its capability to rapidly react with DPE or Fc, whereas no reaction (or an increase in lifetime by complex formation) was observed in the absence of TFA.
In the LFP of 7 and 8, at low laser powers, an apparent recovery was observed at the wavelength of the precursor absorption band. The kinetics of this recovery in case of 7 matched the kinetics of the decay of radical 3, monitored at λ = 315 nm. The fact that the recovery in case of the LFP of 7 was only observed when the laser power was reduced would be consistent with an assignment of this transient process to a reaction of the precursor Barton esters 7 and 8, as use of high laser power would result in complete bleaching of the precursor, so that its subsequent radical chain reactions would be suppressed. It cannot be the growth of thioethers 18 or 19 that is resolved here (Scheme 6), as the UV-vis spectrum of the very similar 2-methylthiopyridine, which may serve as a model system, has neither been reported18a nor is expected to show absorption around λmax = 390 nm. Likewise, the UV-vis spectrum of di-2-pyridyldisulfide (6) has not been reported to show absorption at this wavelength.18b This hypothesis would also be consistent with the observation that the lifetime of the t-butyl radical (10) in the presence of 0.1 mmol L−1 of 7 is τ = 13.5 μs, whereas the lifetime of 3 under analogous conditions ([8] = 0.1 mmol L−1) is only τ = 170 ns. The rate constant of the reaction of 8 with the sterically hindered radical 10 can be expected to be smaller than the rate constant of the analogous reaction of 7 with the more accessible radical 3. Reactions of free radicals like 3 or 10 with Barton esters are expected to proceed via addition to the sulfur atom of the reactive CS bond, followed by cleavage of the rather weak N–O bond in adducts such as 23 (Scheme 6).
It is noted that the barrier for cleavage of the radical adduct 23 is expected to be very small. At the M06-2X/6-311++G(2df,p) (pcm,benzene) level of theory, a relaxed surface scan systematically stretching the N–O bond in 23 resulted in an increase in energy of only 0.8 kcal mol−1 required to reach a point where fragmentation into 19, 10 and CO2 occurred. Hence, any free radical intermediate formed along the reaction coordinate linking 8/10 with 10/19 will be extremely short-lived. We will therefore leave the assignment of the species responsible for the growth at λ = 390 nm open at present, merely noting that kinetic data indicate that this species results from reaction of radicals 3/10 with Barton esters 7/8. The results of an in-depth computational study on reactions of Barton esters with free radicals will be published in due course.
Notably, we have been able to detect radical 3 in the photolysis of 7 in benzene by time resolved EPR (see ESI†). Radical 3, generated by hydrogen abstraction or photoionization, had been detected earlier by EPR spectroscopy and had not demonstrated any fragmentation at ambient temperature.19,20 The barriers for its fragmentation are significant – even the more favourable pathway via radical 20 bears a 21.7 kcal mol−1 free energy of activation in the initial step (see Fig. 8), which would normally make a fragmentation on the ps timescale completely impossible. Our observation of the vinyloxy radical (22), however, necessitates such a rapid fragmentation. The BDE of the N–O bond in Barton esters is small, estimates range between 16 and 28 kcal mol−1.21,22 Given the fact that a 355 nm photon carries an energy of 82 kcal mol−1, and that the decarboxylation of 4 is exothermic by 35 kcal mol−1 (M06-2X/cc-pVQZ), a sizeable excess energy of at least 90 kcal mol−1 remains for the complete photoinduced fragmentation of 7 into thiyl radical 5, dioxolanyl radical 3, and CO2. Even if this excess energy is evenly spread over all vibrational degrees of freedom in all three fragments, sufficient energy (ca. 66 kcal mol−1) will remain in 3 to potentially result in rapid fragmentation. In order to gain further evidence for potential nonstatistical reaction dynamics23 in the fragmentation of 3, we performed Born–Oppenheimer molecular dynamic (BOMD) simulations at the UB3LYP-D3(BJ)/6-31G* level. Four possible conformers of 7 (Fig. 10) were considered with 7c being the most stable one. Quasiclassical trajectories were initialized from normal mode sampling of 7c in the T1 state to approximate a quantum mechanical Boltzmann distribution of the vibrational levels at 298 K. The overlay of the sampled geometries is shown in Fig. 11. A total of 98 trajectories were initiated from the sampled structures.
Fig. 10 Conformers of 7. Numerical values are relative free energies (in kcal mol−1) at 25 °C (B3LYP-D3(BJ)/6-31G* level). |
For all the trajectories, N1–O2 and C3–C4 (Fig. 12) bonds cleave within the fs time scale. The average time to break the N1–O2 bond, as defined by lengthening of the same beyond 1.60 Å, is 196 fs whereas the C3–C4 bond cleaves at 276 fs, considering a bond distance larger than 1.90 Å. Accordingly, the average time lag between the two bond cleavage events is 80 fs. The minimum time required to cleave the N1–O2 bond is 23 fs while the maximum time is 787 fs. Interestingly, the sluggish trajectory took only 45 fs to cleave the second bond. Although all the trajectories ended up eliminating CO2 within 1121 fs, only one went further to cleave the third bond (O8–C7) dynamically. It took a total 630 fs to disrupt all three bonds. Thereafter, the dihedral angle C4–C5–O6–C7 started to rotate but the C5–O6 bond did not break within the simulation time limit of 1352 fs. Thus, the finding indicates that the C5–O6 bond breaking is the rate determining step, in contrast to the result obtained via the calculated barrier heights (cf.Fig. 8).
Footnotes |
† Electronic supplementary information (ESI) available: Synthesis, characterization data, photolysis, EPR, computational data. See DOI: 10.1039/c8cp03311k |
‡ Current address: WestCHEM, School of Chemistry, University of Glasgow, Joseph-Black-Building, University Avenue, Glasgow G12 8QQ, UK. E-mail: goetz.bucher@glasgow.ac.uk |
This journal is © the Owner Societies 2018 |