Mary Alice
Upshur‡
a,
Marvin M.
Vega‡
a,
Ariana Gray
Bé‡
a,
Hilary M.
Chase
a,
Yue
Zhang
b,
Aashish
Tuladhar
c,
Zizwe A.
Chase
c,
Li
Fu
c,
Carlena J.
Ebben
a,
Zheming
Wang
c,
Scot T.
Martin
bd,
Franz M.
Geiger
*a and
Regan J.
Thomson
*a
aDepartment of Chemistry, Northwestern University, Evanston, IL 60208, USA. E-mail: r-thomson@northwestern.edu; f-geiger@northwestern.edu
bJohn A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
cWilliam R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA 99352, USA
dDepartment of Earth and Planetary Sciences, Harvard University, Cambridge, MA 02138, USA
First published on 31st July 2019
Atmospheric aerosol–cloud interactions remain among the least understood processes within the climate system, leaving large uncertainties in the prediction of future climates. In particular, the nature of the surfaces of aerosol particles formed from biogenic terpenes, such as α-pinene, is poorly understood despite the importance of surface phenomena in their formation, growth, radiative properties, and ultimate fate. Herein we report the coupling of a site-specific deuterium labeling strategy with vibrational sum frequency generation (SFG) spectroscopy to probe the surface C–H oscillators in α-pinene-derived secondary organic aerosol material (SOM) generated in an atmospheric flow tube reactor. Three α-pinene isotopologues with methylene bridge, bridgehead methine, allylic, and vinyl deuteration were synthesized and their vapor phase SFG spectra were compared to that of unlabeled α-pinene. Subsequent analysis of the SFG spectra of their corresponding SOM revealed that deuteration of the bridge methylene C–H oscillators present on the cyclobutane ring in α-pinene leads to a considerable signal intensity decrease (ca. 30–40%), meriting speculation that the cyclobutane moiety remains largely intact within the surface bound species present in the SOM formed upon α-pinene oxidation. These insights provide further clues as to the complexity of aerosol particle surfaces, and establish a framework for future investigations of the heterogeneous interactions between precursor terpenes and particle surfaces that lead to aerosol particle growth under dynamically changing conditions in the atmosphere.
The formation, composition, and physicochemical properties of biogenic SOM-containing particles are reasonably well known to depend on the structure and reactivity of terpene volatile organic compound (VOC) precursors, tropospheric conditions of temperature and relative humidity (RH), and the extent of particle phase state heterogeneity.4–7 In contrast, heterogeneous phenomena occurring at particle surfaces remain elusive, despite their intrinsic importance in key formation and growth processes, i.e. condensational growth (gas:particle conversion),8 coagulation (particle:particle merging),9 and cloud droplet activation (cloud condensation nucleation activity).10,11 Scientific insights into the surface-specific molecular origins of SOM-containing particle formation, growth, and cloud activation are therefore critical to the study of such constituents in the climate system, yet direct studies of the particle interface are notoriously challenging to carry out nondestructively with chemical specificity and under tropospherically relevant conditions of temperature, pressure, and RH.4,12,13
Our previous work has demonstrated the applicability of vibrational sum frequency generation (SFG) spectroscopy as a nondestructive nonlinear optical probe to analyze surface-localized oscillators of particle-phase SOM.8,13–15 SFG spectroscopy has revealed molecular information regarding the structure, dynamics, and orientations of terpenes and their oxidation products at atmospheric model interfaces in the C–H and CO vibrational stretching regions.16–19 Notably, the SFG intensity obtained from α-pinene-derived SOM in the C–H region is independent of particle size8,20 and can be readily tracked under varying conditions of RH and temperature in real time.
These initial studies have unlocked the prospect of using SFG spectroscopy to probe the dynamics of VOC uptake onto or into the particle phase under varying RH, temperature, and consequently particle viscosity.13,16 Surface-selectively investigating such early-stage heterogeneous interaction events leading to particle generation would provide critical new insights useful to solving the outstanding uncertainties within aerosol particle mechanisms. Unfortunately, the SFG spectra of most terpenes and oxidation products relevant to biogenic SOM remain complex and congested, with overlapping spectral lineshapes that make vibrational mode assignments challenging.16–19 Moreover, the signal frequencies and intensities obtained from SOM derived from a particular terpene overlap significantly with those of the precursor molecule, making the aforementioned dynamic studies nearly impossible.8,13 For instance, the highest SFG signal intensity obtained from α-pinene-derived SOM collected on a filter sample occurs at 2950 cm−1, shifted only by ∼15 cm−1 from the otherwise similarly shaped SFG spectrum of α-pinene itself.8,21–24 Site-specific deuterium labeling offers the opportunity to address these limitations by “silencing” signal intensity generating C–H oscillators in the SFG spectra of both the precursor terpene and its corresponding SOM, thereby enabling vibrational mode elucidation of the highly congested SFG spectra. We also envision probing the interactions between unlabeled VOCs and labeled particle surfaces using deuterated precursors and preparing labeled SOM from them to distinguish the unlabeled gas-phase molecules from those in the particle phase.
With such goals in mind, we present here results from a detailed isotope-labeling study in which we synthesized novel α-pinene isotopologues and their corresponding SOM, and characterized them using SFG spectroscopy. We expand upon our previous work on methyl C–D α-pinene isotopologues22 by presenting the synthesis and SFG spectroscopy of three new analogues with methylene bridge, bridgehead methine, allylic, and vinyl deuteration. Specifically, we report on the development of a strategy for the de novo synthesis of α-pinene isotopologues adapted from a previous total synthesis completed by Thomas and Fallis.25,26 We then provide SFG spectra of the isotopologues at vapor:solid interfaces to investigate contributions from such C–H stretches on the complex SFG spectrum of α-pinene in the pursuit of identifying high signal intensity generating oscillators. Spectra of the ring-substituted isotopologues revealed up to 80% suppression of the dominant SFG signal intensity (around 2930 cm−1) observed in the unlabeled α-pinene spectrum. Following the preparation of SOM from all accessible isotopologues, the ring-substituted SOM SFG spectra reveal a considerable, albeit far from complete reduction of the 2950 cm−1 signal intensity, indicating that individual oscillators on the SOM surface can be chemically accessed and thus identified. These first attempts toward identifying the surface oscillators in α-pinene-derived SOM will enable future investigations of adsorption, desorption, and reactive uptake of the gas-phase species onto or into aerosol particles, with direct implications for elucidating aerosol formation and growth mechanisms. Beyond such surface spectroscopic studies, our general synthetic route offers the opportunity to utilize site-specific isotopically labeled α-pinene analogues in new mass spectroscopy-coupled chamber studies on SOM formation, growth, and properties within the greater atmospheric community.
Fig. 1 α-Pinene isotopologues synthesized in this study. Atom numbering throughout this manuscript follows the convention indicated by compound 1. |
Synthetic efforts commenced with targeting deuteration of the C6 methylene bridge of α-pinene (Scheme 1), referred to here as α-pinene-d2 (5). Reduction of Hagemann's ester (8) with LiAlD4 generated diol 9 (72% yield), wherein the labeled hydroxyl methyl appendage represented incorporation of deuterium at the site that would eventually become the C6 methylene bridge within the target. Selective oxidation of the allylic hydroxyl group within 9 by DDQ provided enone 10, which was converted to ketone 11 following acetylation and methyl cuprate 1,4-addition (52% yield over three steps from 9). Ketone 11 was then transformed into tosylate 12 in 73% yield by first conducting a base-promoted aldol condensation with benzaldehyde followed by tosylation of the liberated primary alcohol. As demonstrated by Thomas and Fallis, installation of the benzylidene group on the less hindered C3 side of the ketone within 11 was necessary to enforce production of the desired regioisomer during the subsequent formation of the cyclobutane ring. Accordingly, exposure of tosylate 12 to sodium hydride in dimethoxyethane at 80 °C afforded methylene bridge labeled bicycloheptane 13 in 88% yield. Removal of the benzylidene blocking group was accomplished using γ-aminobutyric acid (GABA) as a retro–aldol catalyst, as described by Thomas and Fallis.25 Somewhat disappointingly, the best reproducible yield obtained for nopinone-d2 (14) under optimized conditions was 41% yield (71% yield based on recovered starting material), in part due to product volatility. With advanced material in hand, however, the synthesis of α-pinene-d2 (5) was completed via triflate 15, which underwent methylation with methyl cuprate under the procedure we had developed previously (51% yield over two steps from 14).22
Having established the viability of modifying the Thomas and Fallis synthesis to generate methylene bridge labeled analogue, α-pinene-d2 (5), we turned our attention to installing deuterium at the two bridgehead positions (Scheme 2). We reasoned that deprotonation of ester 8 to yield an extended enolate might allow for direct labeling at the desired position(s) upon quenching with an appropriate source of deuterium. To that end, Hagemann's ester (8) was treated with excess sodium hydride to generate a solution of the requisite enolate, which was quenched by the slow addition of D2O and then carefully acidified with deuterated sulfuric acid. Under optimized conditions, we found that ester 16 could be generated in 88% yield with complete incorporation of two deuterium atoms at both C1 and C5. We were initially surprised to see labeling of the enone vinyl position, but this can be rationalized by considering the highly delocalized nature of the enolate formed from Hagemann's ester. It was important that the reaction mixture be kept at a pH of ∼2.0 during the quench in order to ensure complete deuteration at both positions, lest an inseparable mixture of partially labeled species form. From this point deuterated Hagemann's ester (8) was treated similarly to the previous modified synthesis presented. Reduction with LiAlH4 followed by allylic oxidation with DDQ provided enone 17 (85% over two steps). Acetylation of the free hydroxyl and conjugate addition with the cuprate reagent produced acetate 19. When benzylidene formation was conducted under the original conditions, hydrogen–deuterium exchange had occurred adjacent to the ketone, removing the C1 deuterium label previously inserted. The exchange could be easily avoided by using 10% NaOD during the aldol condensation, allowing ketone 20 to be isolated with complete deuterium installation at C1. Subsequent treatment of ketone 20 with sodium hydride induced cyclization as before to generate cyclobutane 21 in 89% yield, which now possessed deuterium labels at both C1 and C5 bridgehead carbons.
Removing the benzylidene group within 21 proved to be problematic due to unforeseen loss of the bridgehead deuterium atoms. GABA was found to facilitate this process, as was the use of NaOH and regular DMSO. To circumvent this issue, we conducted the retro–aldol reaction in the absence of GABA using NaOD and deuterated DMSO as the solvent. Under our best conditions found, the benzilidene group could be cleaved without loss of bridgehead deuterium atoms in 20% yield (41% brsm), however, we now observed additional deuteration α- and β- to the ketone at C3 and C4, respectively (i.e., 22 was the major product). While deuteration next to the ketone was not unexpected, we were surprised to see deuteration at the β-position. We reason that deprotonation of 21 at the C4 allylic position must be facilitating this exchange process before the desired retro–aldol reaction takes place. Unfortunately, despite much experimentation, the efficiency of this allylic exchange could not be rendered 100% effective; 1H NMR spectroscopy of 22 shows 50–70% of double to single deuteration at C4. This inefficiency, however, was not expected to be problematic in the planned SOM synthesis and spectroscopic analysis, and so we pressed forward to generate the desired α-pinene-d5 (6) as in the prior analogue synthesis. We also prepared α-pinene-d7 (7), substituting LiAlD4 for the reduction of deuterated Hagemann's ester (8) to generate alcohol 18, which was subjected to the same sequence as described for α-pinene-d5 (6, see ESI† for complete details).
Fig. 2 shows the 0.6 cm−1 resolution SFG spectra of vapors of isotopologues 5–7 and of unlabeled α-pinene in contact with fused silica (Fig. S2† shows that spectra obtained on CaF2 are comparable, while Fig. S3† shows standard-resolution spectra obtained at NU). The SFG response of unlabeled α-pinene is predominantly observed within two vibrational frequency regions centered around 2930 cm−1 and 2880 cm−1. Previous atomistic simulation studies accounting for Fermi resonances and nonlocal modes have reported that the bridge methylene group on α-pinene contributes considerably to the dominant spectral feature around 2930–2940 cm−1.27 Specifically, static DFT calculations indicate that the 2930–2940 cm−1 peak is 77% composed of the bridgehead methylene symmetric stretch and a Fermi resonance (FR), with the remaining 23% of the SFG response at this frequency due to overtone bands of the gem-dimethyl CH3 symmetric stretches (19%), the asymmetric stretch of the vinyl CH3 (1%), and other stretches (3%).27 MD simulations indicate that the bridgehead methylene symmetric stretch and FR contribute 79%, with the remaining contributions from the overtone band FR of the gem-dimethyl CH3 symmetric stretches (12%), the gem-dimethyl CH3 asymmetric stretch (3%), vinyl CH3 asymmetric stretch (1%), bridging methine CHs (4%), and other modes (1%).27 At 2880 cm−1, static DFT calculations indicate that ∼75% of the SFG response is due to the gem-dimethyl CH3 symmetric stretches and their Fermi resonance, with the remainder being due to the allylic CH2 symmetric and asymmetric stretching modes (∼20%) and other modes.27 The C–H symmetric stretch on the vinyl CH3 group was found to contribute to the 2880 cm−1 peak by a somewhat greater extent (ca. 2%) in the ensemble-averaged spectrum from MD simulations than in the static DFT calculations.27 Isotopologues 5–7 each exhibit a marked decrease in the SFG response in the 2930 cm−1 region relative to unlabeled α-pinene, while their methyl symmetric stretching modes between ∼2850 cm−1 and ∼2880 cm−1 are largely unchanged, as expected. The latter outcome suggests that the molecular orientation distributions of the isotopologues 5–7 are comparable to that of the unlabeled compound, which is an important result given the dependence of SFG responses on molecular orientation.
Fig. 3 shows that contributions to the 2930 cm−1 signal from the methyl groups seem to be minor if not negligible, as deuterating one or two of those groups (isotopologues 2–4) still produces the highest peak intensity at that frequency.22 The overall reduction of signal intensity in the 2930 cm−1 region exhibited in the spectra of isotopologues 5 and 7 is consistent with the previously reported contributions from deuteration of the bridge methylene group,27 but is rather unexpected for isotopologue 6, given its bridge methylene group is not deuterated. The results shown in Fig. 2 and 3 then suggest that the two bridging methine groups on the ring (deuterated in isotopologues 6 and 7) also contribute somewhat to the overall SFG response at ∼2930 cm−1, perhaps to a larger extent than the ca. 4% previously reported by MD simulations.27 Similar interactions between C–H oscillators in isotopologues of isoprene epoxydiol have already been described earlier.28 Additionally, given the aforementioned coherent nature of the SFG process, this outcome may reveal how the vector field addition that produces the observed SFG signal may be ultimately influenced by deuteration pattern. The band centered at 2840 cm−1 disappears completely in the spectra of isotopologues 6 and 7, likely due to silencing of the allylic methylene group, which agrees with predictions from the discussed computations,27 though we cannot unambiguously assign this band given that there exists partial deuteration at that position (ca. 50–70% of double to single deuteration at C4). Synthetic efforts are underway to selectively deuterate the allylic and vinyl positions of α-pinene, though such studies are beyond the scope of this work. Nevertheless, we caution that the unambiguous assignment of the 2930 cm−1 peak will continue to be challenging.
Fig. 4 shows the standard resolution SFG spectra of SOM derived from isotopologues 2–7 compared to that of SOM derived from unlabeled α-pinene. As perhaps anticipated, the SFG spectral lineshapes for SOM derived from isotopologues 2–4 are largely invariant with methyl group deuteration (Fig. 4A). In contrast, and surprisingly to us, none of the complex lineshape changes observed for the ring-labeled isotopologues (5–7) are observed for the SOM prepared from them (Fig. 4B). Yet, we note a 30 to 40 percent signal intensity reduction in the 2950 cm−1 region for SOM prepared from isotopologues 5 and 7, which are the only isotopologues we studied that feature a deuterated methylene group on the four-membered ring. It appears that this reduction originates from deuterating the CH2 group on the four-membered ring, since isotope labeling of the methine positions (isotopologue 6) on the same ring do not result in any changes in the surface SFG spectrum of the corresponding SOM. We conclude that the high-intensity peak seen in the SFG spectra of α-pinene-derived SOM at ∼2950 cm−1 contains considerable contributions (ca. 30%) from the methylene group on the four-membered ring, while the methyl groups do not contribute.
Together, our findings support the hypothesis that the surface-active constituents of α-pinene-derived SOM conserve the cyclobutane functionality, as the bridge methylene group was found to contribute ca. 30% of spectral intensity in the SFG response at ∼2950 cm−1. While isotope labeling of these positions in isotopologues 5 and 7 resulted in significant changes in the signal intensities and lineshapes within the vapor phase SFG spectra, it appears as though the bridge methylene oscillators in the SOM do not contribute nearly as much intensity as initially expected to the ssp-polarized SFG spectra. It is possible that this outcome may arise from differences in the orientation distributions of the cyclobutane moieties present in the SOM surface-active oxidation products compared to those present in the precursor isotopologues. However, we expect it is more likely that the remaining signal intensity around 2950 cm−1 in the SFG spectra of SOM prepared from the bridge methylene deuterated isotopologues 5 and 7 is attributed to C–H oscillators of α-pinene and its oxidation products that have not yet been identified through deuteration.
With this goal in mind, we provided SFG spectra of the ring- and methyl-substituted isotopologues at vapor:solid interfaces to observe the relative contributions from silenced C–H oscillators on the SFG spectrum of α-pinene. The ring-substituted isotopologues revealed a significant (50–80%) suppression of the SFG signal intensity in the broad band, centered around 2930 cm−1, that is emblematic of unlabeled α-pinene. In contrast, the methyl-deuterated isotopologues we previously examined22 do not show such a signal intensity reduction, emphasizing the coherent nature of the SFG signal generation process from the various C–H oscillators in the highly strained bicyclic system of α-pinene, and the need for complementary experimental and computational studies for further spectral deconvolution. Future work will include accessing the C–D region, carrying out internal heterodyne phase-resolved measurements, improving SFG calculations, and synthesizing other α-pinene isotopologues.
By subsequently preparing and collecting SFG spectra of SOM derived from the α-pinene isotopologues, we then made the first attempts thus far toward identifying the surface oscillators in α-pinene-derived SOM. From the SOM SFG data, we speculate the presence of the dimethylcyclobutane ring at the SOM surface, as SOM prepared from the ozone oxidation of two isotopologues featuring a deuterated bridge methylene group shows up to 40% less signal intensity than SOM prepared from the methyl-substituted isotopologues to which we had access. Yet, given the substantially higher (up to 80%) SFG signal intensity reductions observed in the SOM precursor isotopologues that are bridge-deuterated, we conclude that SOM formation leads to other SFG active modes in the C–H stretching region at the particle surface.
Ultimately, our synthesis of isotopically labeled SOM and analysis of their surface-active constituents opens the door to future investigations of the dynamics of gas-phase VOC adsorption on and/or uptake into the particle phase under varying tropospherically realistic conditions.13,16 Beyond the SFG studies presented here, our syntheses of α-pinene isotopologues have been used in mechanistic chamber-based studies, including the demonstration of particle–particle chemical exchange using single-particle mass spectrometry,34 as well as identification of intramolecular hydrogen migration in autoxidation leading to highly oxygenated multifunctional organic compounds.35 More broadly speaking, synthetically accessing isotopologues of terpene hydrocarbons that are not commercially available offers the scientific community at large new opportunities to elucidate mechanistic questions in aerosol formation and growth, terpenoid biosynthesis, sustainable biofuels, and food and flavor chemistry, among other areas of study.34,36–39
The standard resolution spectra reported are an average of 7 individual spectra each taken for 2 minutes each. For the standard resolution spectra, the IR energy distribution was accounted for by normalizing the data to the ppp-polarized nonresonant SFG signal of gold deposited onto fused silica. Frequencies were calibrated using a polystyrene film.13,14,45,48 The high-resolution spectra are reported as an average of 3–5 individual acquisitions each recorded for 5 minutes. The data points were binned by 5 points, or by 1.73 cm−1, in Igor Pro Version 6.11 (WaveMetrics, Lake Oswego, OR, USA). SFG intensities in the 0.6 cm−1 resolution spectra were normalized to the ppp-polarized nonresonant SFG response of clean z-cut α-quartz, and frequencies were calibrated again using polystyrene.
Footnotes |
† Electronic supplementary information (ESI) available: Detailed information for the synthesis of isotopologues, SOM collection, and supplementary SFG data. See DOI: 10.1039/c9sc02399b |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2019 |