Limononic acid oxidation by hydroxyl radicals and ozone in the

after elimination of formaldehyde from the primary alkoxy radical a C9 α-hydroxyalkyl radical is. 249 formed. The C9 α-hydroxyalkyl radical can reac...
0 downloads 8 Views 2MB Size
Subscriber access provided by UNIV OF DURHAM

Article

Limononic acid oxidation by hydroxyl radicals and ozone in the aqueous phase Bartlomiej Witkowski, sara jurdana, and Tomasz Gierczak Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04867 • Publication Date (Web): 14 Feb 2018 Downloaded from http://pubs.acs.org on February 20, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 20

Environmental Science & Technology

1

Limononic acid oxidation by hydroxyl radicals and ozone in the

2

aqueous phase

3

Bartłomiej Witkowski, 1* Sara Jurdana, 1,2 Tomasz Gierczak 1 1

4 5

University of Warsaw, Faculty of Chemistry,

Al. Żwirki i Wigury 101, 02-089 Warsaw, Poland

6

7 8 9 10 11 12 13 14 15 16 17 18 19

Keywords: Limononic acid, limonene, secondary organic aerosol, hydroxyl radicals, ozone, relative

20

rate, tandem mass spectrometry

21

*Corresponding author e-mail: [email protected]

22

2

23

Rijeka, Croatia

Permanent address: University of Rijeka, Department of Biotechnology, Radmile Matejčić 2, 51000

1

ACS Paragon Plus Environment

Environmental Science & Technology

24

Page 2 of 20

Abstract

25

Kinetics and mechanism of limononic acid, (3-isopropenyl-6-oxoheptanoic acid, LA) oxidation

26

by hydroxyl radicals (OH) and ozone (O3) were studied in the aqueous phase at 298±2 K. These

27

reactions were investigated using liquid chromatography coupled to the electrospray ionization and

28

quadrupole tandem mass spectrometry (LC-ESI/MS/MS). The rate coefficients determined for LA +

29

OH reaction were: 1.3 ± 0.3 × 1010 M-1s-1 at pH=2 and 5.7 ± 0.6 × 109 M-1s-1 at pH=10. The rate

30

coefficient determined for LA ozonolysis was 4.2 ± 0.2 × 104 M-1s-1 at pH=2. Calculated Henry’s Law

31

constant (H) for LA was ca. 6.3 × 106 M × atm-1 thereby indicating that in fogs and clouds with LWC =

32

0.3-0.5 g × m-3 LA will reside entirely in the aqueous-phase. Calculated atmospheric lifetimes due to

33

reaction with OH and O3 strongly indicate that aqueous-phase oxidation can be important for LA

34

under realistic atmospheric conditions. Under acidic conditions, the aqueous-phase oxidation of LA

35

by OH will dominate over reaction with O3, whereas the opposite is more likely when pH ≥ 4.5. The

36

aqueous-phase oxidation of LA produced keto-limononic acid and a number of low-volatility

37

products, such as hydroperoxy-LA and α-hydroxyhydroperoxides.

38 39 40 41 42 43 44 45 46 47 48 49 50 51 2

ACS Paragon Plus Environment

Page 3 of 20

Environmental Science & Technology

52

1. Introduction

53

Secondary organic aerosols (SOA) are formed due to nucleation or partitioning of the semi1, 2

54

volatile oxidation products of VOCs in accordance with the partitioning theory.

55

methane VOCs budget is dominated by biogenic emission of isoprene and monoterpenes: α, β-

56

pinene, and limonene.

57

(BVOC) contributes significantly to the secondary organic aerosols (SOA) formation.

58

because of our limited understanding of the SOA formation mechanisms, impact of the terpene

59

biogenic emission on the climate is still uncertain. 5-10

3, 4

Globally, non-

Gas-phase oxidation of these atmospherically abundant biogenic VOC 5

However,

60

Limonene is a more efficient secondary organic aerosol (SOA) precusor than α and β-pinene

61

because it possess two carbon-carbon double bonds. 11-15 A significant fraction of the first generation

62

oxidation products of limonene will retain the less reactive,

63

Consequently, these unsaturated, semi-volatile compounds can undergo further oxidation in the gas,

64

particle or aqueous phase thereby increasing the overall SOA yield. 14, 15

65

11, 12, 16

exocyclic double bond.

17-20

Recently, it has been recognized that the aqueous-phase processing of the semi-volatile 7, 21-23

66

monoterpene derivatives can contribute to SOA

67

characterized.

68

phase processing of limonene oxidation products.

69

products of limonene can react with ozone (O3) and hydroxyl radicals (OH) in the aqueous-phase. If

70

relevant under atmospheric conditions, these aqueous-phase reaction may lead to aqueous-phase

71

SOA (aqSOA).

72

been focused on the less reactive, saturated cis-pinonic acid (CPA) and other saturated derivatives of

73

α-pinene. 26-29

74

10, 21, 22, 24

7, 14, 21-23

but these processes are very poorly

Also, little attention has been directed to the investigation of the aqueous25

The unsaturated, first-generation oxidation

However, to-date, most studies of the terpene derived aqSOA precursors have

Limononic acid (3-isopropenyl-6-oxoheptanoic acid, LA) is a first-generation product of limonene 12, 17, 30

28

75

oxidation in the gas-phase

76

double bond and the estimated Henry's law constant (H) for this compound is 6.3 × 106 M × atm-1. 31

77

In fog and clouds with LWC = 0.3-0.5 g × m-3 LA will reside entirely in the aqueous-phase. Thus, LA is a

78

good model compound for studying

79

volatile derivatives of monoterpenes. In the gas-phase LA reacts with O3 (k≈5 × 10-19 cm3 × molecule-1

80

s-1)

81

reactions in the aqueous-phase are currently unknown.

11, 12, 15, 16

and also a major product of CPA photolysis.

aqSOA

LA has a exocyclic

formation due to oxidation of the unsaturated, semi-

and OH (k=7 × 1011 cm3 × molecule-1 s-1)

32, 33

but the rate coefficients for these two

82

The aim of this work is to investigate the kinetics and mechanism of LA reaction with OH

83

(reaction 1) and with O3 (reaction 2) in the aqueous-phase. These reactions were studied in a batch

84

reactor using liquid chromatography coupled to the electrospray ionization tandem mass

85

spectrometry (LC-ESI/MS/MS). 3

ACS Paragon Plus Environment

Environmental Science & Technology

86 87 88

Page 4 of 20

2. Experimental List of reagents is provided in section S1 of the Supplementary Information (SI). 2.1. Limononic acid synthesis and purification

89

LA was synthesized by irradiating CPA solution in deionized (DI) water/MeOH with a 254 nm

90

mercury Pen-Ray lamp. To confirm the structure of LA as the major product of CPA Norrish-II

91

isomerization,

92

agreement with the previously published data. 28 Details of this synthesis are provided in section S2.

93

2.2. Kinetics of limononic acid reaction with OH and O3

94 95

28 1

H NMR spectrum was acquired. LA 1H NMR spectrum obtained is in a good

The rate coefficients for reaction (1) reaction (2) in the aqueous-phase were obtained using competition kinetics by monitoring the relative loss of LA and the reference compound. Limononic acid + OH → products 1 Limononic acid + O → products 2

96 97

If LA and the reference compound are lost only due to reaction with OH or O3, the kLA can be calculated using eq. (I) Ln 

98 99

[Limononic acid] k ! [Reference] Ln  =  I [Limononic acid] k "#$ [Reference]

[Limononic acid] and [Reference] are the initial (0) and intermediate (t) concentrations, kLA and kref are the bimolecular rate coefficients for LA and reference compound.

100

Caffeine (kOH=6.9 ± 0.7 × 109 M−1 s−1 at pH=2), phenylalanine (kOH=9 ± 1 × 109 M−1 s−1 at pH=10

101

and 5.7 ± 0.5 × 109 M−1 s−1 at pH=2), p-toluic acid (kOH= 8 ± 0.8 × 109 M−1 s−1 at pH=10) 34 were chosen

102

as the reference compounds to study kinetics of reaction (1) – k1. Gallic acid (kozone= 9.7 × 104 M−1 s−1

103

at pH =2) 27, 38 was used to study kinetics of reaction (2) – k2. It was assumed that one molecule of O3

104

reacts with one molecule of gallic acid and LA. 35

105 106

2.3. Oxidation of limononic acid by OH The aqueous-phase photo-reactor was described in our previous study.

36

Briefly, the photo-

107

reactor consisted of a 50 ml Pyrex vessel and the reaction mixture was stirred with magnetic stirrer.

108

The reaction solution was irradiated with a UVAHAND 250 GS H1/BL lamp (Honle UV technology, 310

109

W) to generate OH by H2O2 photolysis.

4

ACS Paragon Plus Environment

Page 5 of 20

Environmental Science & Technology

110

The reaction mixture consisted of 25 ml of 17 mM H2O2 solution in DI water. Concentrations of

111

LA and the reference compounds were ca. 30 µM. Solution pH was adjusted by adding HCl or NaOH.

112

OH concentration in the absence of other reactants was ca. 5 × 10-12 M.

113

concentrations were used in some experiments to study the mechanism and identify reaction (1)

114

products. Aliquots of the reaction mixture (100 µL) were sampled from the reactor, incubated with

115

the equal volume of the catalase solution (≈ 0.1 mg/ml) and injected into LC/MS.

17

Slightly higher reactants

116

2.4. Ozonolysis of limononic acid

117

Kinetics of reaction (2) was studied by adding different volumes of the ozone solution (ca. 6

118

µM, determined with the indigo method 37) to the aqueous solution of LA and gallic acid (2 µM each).

119

38

120

(pH=2). 38, 39 The samples always contained excess reactants relative to O3 and the tert-butyl alcohol

121

(t-BuOH) was used as OH scavenger (10-20 mM). 38, 39

Ozone solution was prepared by bubbling the ozone through phosphoric acid solution in DI water

122

Mechanism of reaction (2) was studied by slowly bubbling the ozone (≈ 5 ml/min) through

123

solution of LA (50 µM) in DI water. 150 µL aliquots of the reaction mixture were sampled from the

124

bubbler and subjected to the LC/MS analysis.

125

2.5. Liquid chromatography coupled to the mass spectrometry

126

LC/MS analysis were carried out with LC20A liquid chromatograph (Shimadzu) coupled to the 17

127

QTRAP 3200 (AB Sciex) triple quadrupole mass spectrometer as previously described.

128

products identification a 55 min gradient elution program was used. Shorter gradient elution

129

program was used for the kinetic analysis (section S3).

For the

130

Mass spectra were acquired in the scan mode in the mass range 50-700 m/z and in the

131

multiple reaction monitoring (MRM) mode. The ESI conditions were as follows: capillary voltage was -

132

4.5 kV and 5.5 kV in the negative and positive ionization modes. Ion source temperature was 450 °C,

133

nitrogen was used as curtain (3 × 105 Pa), auxiliary (3 × 105 Pa) and collision gas. Fragmentation

134

spectra of reactions (1) and (2) products were acquired during chromatographic run, after the

135

analytes were separated by LC. Ion lenses voltages for the MRM mode were adjusted by directly

136

injecting the analytes into the ESI ion source.17

137

2.6.

138

Control experiments and error analysis Control experiments were carried out to search for systematic errors.

17

Briefly, UV-vis

139

irradiation, H2O2, phosphoric acid and t-BuOH by themselves did not caused any measureable

140

decomposition of LA and the reference compounds within the time-scale of the experiment as 5

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 20

141

described in section S4. Unfortunately, reaction (2) kinetics was only investigated under acidic

142

conditions due to some experimental problems (data not shown).

143

It was also impossible to avoid some minor repartitioning of LA into the gas-phase under

144

acidic pH conditions. To account for the non-OH losses of LA under acidic pH conditions, two Pyrex

145

vessels were placed in front of the lamp after carefully adjusting their positions. The reaction mixture

146

volumes and compositions were identical in both bottles but H2O2 was added only to one. During

147

these experiments, the reaction solutions were sampled simultaneously from both bottles. In the

148

second set of control experiments, the reaction vessel without the H2O2 was place at exactly the

149

same distance from the lamp immediately after the photo-oxidation experiment (with the H2O2

150

added to the reaction mixture) was completed. The results from both sets of these control

151

experiments were averaged. During the experiments under acidic conditions, LA was lost due to

152

reaction (1) and also due to repartitioning into the gas phase. Using the first order loss rates (min-1),

153

it was possible to account for the non-OH losses of LA as shown in Fig. S8.

154

The total uncertainty in the bimolecular rate coefficients listed in Table 1 is a sum of

155

precision of the kLA/kref data fits obtained from linear regression analysis (standard error) as well as

156

uncertainties of the kref values (±10%). 34, 36 The uncertainties of the analytes measurements was ±5%

157

as estimated from the reproducibility of the integrated MRM peak areas (2σ). Thus, the combined

158

uncertainty for k1 value calculated for basic pH conditions was ±10% and ± 20 % for reaction (1)

159

under acidic conditions.

160

3. Results and discussion

161

3.1.

162

Determination of the reaction rate coefficients with OH and O3 Relative kinetic plots for reaction (1) are shown in Fig. 1.

6

ACS Paragon Plus Environment

Page 7 of 20

Environmental Science & Technology

163 164

Figure 1 Relative kinetic data for reaction (1)

165

Linear regression analysis (straight lines in Fig. 1) was used to calculate k1/kref ratios listed in

166 167

Table 1. Table 1 Measured k1 and k2

Reference

168

R2

Slope, k1/krefa Caffeine 0.996 1.45 ± 0.04 Phenylalanine 0.994 3.00 ± 0.09 p-Toluic acid 0.993 0.72 ± 0.02 Phenylalanine 0.995 0.62 ± 0.02 a ± standard uncertainty; b listed by Buxton et al. 34

kref (109 M−1 s−1)b 6.9 ± 0.7 5.7 ± 0.5 8.0 ± 0.8 9.0 ± 0.9

pH 2 10

Average k1 (1010 M−1 s−1) k1 for pH=2 1.3 ± 0.3 k1 for pH=10 0.57 ± 0.06

169

As shown in Fig. 1, the kinetic plots were linear: R2 > 0.99. Calculated k1 values listed in Table

170

1 for pH=10 were very consistent for both reference compounds, and the average value of 5.70 ±

171

0.06 × 109 M−1 s−1 was obtained. Because there was a need to correct for a non-OH losses of LA under

172

acidic conditions (section 2.6), the uncertainty of k1 obtained for pH=2 was higher. For the

173

undissociated LA, an average k1 value of 1.3 ± 0.3 × 1010 M−1 s−1 was obtained.

174

For CPA, a saturated compound, considerably lower k1 values of 3.6 ± 0.3 and 3.3 ± 0.5 × 109

175

M−1 s−1 for pH=2 and 3.0 ± 0.3 × 109 M−1 s−1 for pH = 10 were recently reported.

176

increased reactivity is most likely due to the OH addition to the double bond. Using aqueous-phase

177

structure–activity relationship (SAR)

178

will react with OH via addition mechanism (section S5).

17, 18

29, 36

Thus, this

it can be roughly estimated that ca. 70 % of LA molecules

7

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 20

179

It is difficult to discuss the pH-dependence of k1, especially since the current literature

180

database for keto-acids oxidation by OH is rather limited. 34, 36 Similar results were presented in our

181

previous study for OH reaction with CPA, where the apparent lack of pH dependence on the kOH was

182

associated with the unexpectedly high acidity of the hydrogen atom in α position to the carbonyl

183

group. 36

184

For reaction (2), k2 = 4.2 ± 0.4 × 104 M−1 s−1 was obtained as shown in Fig. S12. As expected,

185

ozonolysis of the exocyclic double bond of LA was significantly slower than the ozonolysis of

186

endocyclic double bonds. 40, 41

187 188 189

The data acquired was used to calculate atmospheric lifetimes of LA - section 3.3. 3.2. Mechanisms of limononic acid oxidation by OH and O3 The products of reaction (1) and (2) detected using LC/MS as shown in Fig 2. 2

190 191 192

Figure 2 Chromatograms (total ion current): red - reaction mixtures at t=0 and, blue – 60 min after reactions initiation

193

Fragmentation spectra of the products detected with LC/MS were used to propose the

194

mechanisms of reactions (1) and (2). Molar yields (section S7) were estimated using integrated peak

195

areas for the [M-H]- ions. Since the lack of standards the yields calculated using this method are most

196

likely underestimated. The ESI response increases together with the percentage of volatile organic

197

modifier (ACN), 42 hence it was lower in the beginning of the gradient run - see section S3. Thus, the

8

ACS Paragon Plus Environment

Page 9 of 20

Environmental Science & Technology

198

highly oxygenated compounds that were eluted earlier than limononic acid (tr =27.5 min) were

199

ionized less efficiently. 3.2.1. Aqueous phase oxidation of limononic acid by OH

200 201 202

Extracted ion chromatograms (EICs) and fragmentation spectra of the reaction (1) products are shown in Fig. 3.

203 204 205

Figure 3 EICs and fragmentation spectra of reaction (1) products

206

As it is shown in Fig. 3, the major product of reaction (1) was keto-LA – see section S7. Keto-

207

LA was also identified using the data reported in our previous study of limonene SOA composition. 17

208

The small quantities of compounds with MW 234 Da were most likely produced due to the non-

209

radical addition of H2O2 to LA as it could be inferred from the retention time of these peaks – see

210

section S8.

211

To further investigate the non-radical addition of H2O2 to LA, a series of experiments was carried

212

out, using linear carboxylic acids, carboxylic acids containing hydroxyl and carbonyl groups as well as

213

dicarboxylic acids – see section S8. These carboxylic acids were dissolved in the H2O/H2O2 mixture,

214

evaporated to dryness and re-dissolved in pure water. H2O2 only reacted with the keto-acids (Table

215

S3), including LA. As shown in Fig. 3, LA + H2O2 reaction generated products detected as negatively 9

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 20

216

charged ions with m/z 233. Retention times of some of these ions overlapped with the products of

217

reaction (1), as shown in Fig.3. Interestingly, H2O2 did not cause any measurable decay of LA during

218

control experiments but small peaks as m/z 233 ions were detected in these experiments. Perhaps a

219

substantial amount of H2O2 was produced during the photooxidation experiments, thereby increasing

220

the yields of these non-radical products.

221

The other chromatographic peaks for the m/z 185 ion are most likely a result of the in-

222

source fragmentation as m/z 233 ions (the same retention times of the chromatographic peaks). As

223

shown in Fig. 3, compounds detected as m/z 233 ions readily produced m/z 185 fragmentation ion

224

following collision induced dissociation (CID). A number of minor products of reaction (1) were also

225

detected as shown in Fig. 4.

226 227

Figure 4 EICs and fragmentation spectra of the additional products of reaction (1)

228

As shown in Fig. 4, minor products of reaction (1) were detected as m/z 215 (two isomers),

229

201 (three isomers) and 197 (two isomers). Compounds detected as m/z 215 and 197 ions are most

230

likely carboxylic acids as indicated by the neutral loss of 44 Da (CO2). Neutral losses of 18 Da (H2O)

231

and 62 Da (CO2 + H2O) from m/z 215, 201 as well as from 197 ions indicate that these compounds

232

could be hydroxy carboxylic acids. These compounds were most likely produced by the abstraction

233

mechanism that is shown Fig. 5

234

The proposed mechanism for reaction (1) is shown in Fig. 5.

235 10

ACS Paragon Plus Environment

Page 11 of 20

Environmental Science & Technology

236 237 238 239

Figure 5 Proposed mechanism of LA oxidation by OH, products detected with LC/MS are shown in blue As shown in Fig.5 addition mechanism is expected to dominate over H atom abstraction as 17, 18

240

already discussed in section 3.1 and section S5.

241

experimental data, since keto-LA (major product detected) is formed by conversion of the terminal

242

double bond to a carbonyl.

This assumption is consistent with the

243

As shown in Fig. 5, the two peroxy radicals (RO2) are produced due to OH addition to the

244

double bond and reaction of the hydroxylalkyl radicals with O2. The formation of peroxy radical (1)

245

should be favored because it is produced from a tertiary alkyl radical.

246

with each other to form the unstable tetroxides (the water-cage effect).

247

decompose to primary or tertiary alkoxy radicals. Keto-LA can be produced following •CH2OH

248

elimination from the tertiary alkoxy radical – this pathway is most likely favored. 43, 44 Alternatively,

249

after elimination of formaldehyde from the primary alkoxy radical a C9 α-hydroxyalkyl radical is

250

formed. The C9 α-hydroxyalkyl radical can react with O2 thereby yielding keto-LA by eliminating HO2

43

The two RO2 radicals react 44

Both tetroxides can

11

ACS Paragon Plus Environment

Environmental Science & Technology

251 252

45

or to produce α-hydroxyalkyl alkoxy radical.

alkoxy radical by OH elimination.

44

Page 12 of 20

Keto-LA is then produced from the α-hydroxyalkyl

44

253

Minor products of reaction (1) are most likely produced by addition pathway and subsequent

254

oxidation of the terminal double-bonds of the first generation products. Hydroperoxy (MW = 216

255

Da), hydroxyl (MW = 200 Da) and carbonyl (MW = 198 Da) – substituted LA derivatives are produced

256

due to reaction of the peroxy radicals with HO2 or RO2. 36 Formation mechanism of hydroperoxy-LA

257

was most likely similar to the formation of hydroperoxy-CPA that was previously described. 36, 46, 47

258

Estimated LA conversion yield to keto-LA is rather low (ca. 10%). Clearly other products of

259

reaction (1) are formed, and these additional compounds were most likely not detected with LC/MS

260

(e.g. carbonyl compounds and lower-MW carboxylic acids). This might be a result of the

261

fragmentation of the alkoxy radicals. In case of carboxylic acids, it has been recently demonstrated

262

that fragmentation of alkoxy radicals is enhanced for the carboxylate ions. 48 Clearly, there is a need

263

to quantify the lower-MW products of reaction (1) for LA and similar compounds to further

264

investigate the impact of the carboxylic acid/carboxylate equilibrium on the alkoxy radicals

265

fragmentation.

266 267

3.2.2. Aqueous phase ozonolysis of limononic acid EICs and fragmentation spectra of the reaction (2) products are shown in Fig. 6.

268 269

Figure 6 EIC chromatograms and fragmentation spectra of the reaction (2) products

12

ACS Paragon Plus Environment

Page 13 of 20

Environmental Science & Technology

270

As shown in Fig. 6, the major product of reaction (2) is keto-LA; only a single chromatographic

271

peak was detected for m/z 185 ion. Fragmentation spectra and retention time of keto-LA were

272

consistent with data reported in our previous study.

273

and 34 Da (H2O2) from the m/z 201 ion indicated that the second product most likely contains

274

hydroperoxide moiety. This characteristic fragmentation mechanism confirms that at least the

275

compounds with MW 202 Da might be acyloxyhydroperoxy aldehydes. Fragmentation of the second

276

m/z 201 ion was very difficult (low intensity of the fragmentation ions), which might be consistent

277

with the formation of secondary ozonide (data not shown).

278

17

As shown in Fig.6 neutral losses of 18 (H2O)

A proposed mechanism for reaction (2) is shown in Fig. 7.

279 Figure 7 Proposed mechanism of LA ozonolysis, products detected with LC/MS are shown in

280 281 282

blue As it is shown in Fig. 7, when O3 adds to the terminal double bond, the unstable primary 40, 41, 49

283

ozonide (POZ) is formed.

284

formaldehyde (strongly favored) 49 or C1 CI and keto-LA (minor, not shown). 40 In the aqueous-phase,

285

the C9 CI is efficiently stabilized and a stabilized CI (SCI) is formed. 47, 57 In the laboratory reactor, the

286

C9 SCI can react with water to produce the α-hydroxyhydroperoxide (MW 220 Da). The α-

287

hydroxyhydroperoxides are unstable as confirmed by the very low intensity of the chromatographic

288

peak for the m/z 219 ion (data not shown) and decompose to keto-LA and H2O2. 40 The experimental

289

results also indicated that small fraction of the C9 SCIs can self-react, thus two products with MW 202

290

Da are observed as shown in Fig 7.

291

3.3. Atmospheric implications

POZ can decompose to C9 Criegee intermediate (CI) and

13

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 20

292

The following data was used to calculate the gas-phase lifetimes of LA that are shown in Fig. 8:

293

kOH = 7 × 10-11 cm3 × molecule-1 s-1, (calculated with SAR) 32, 33 and kozone = 5 × 10-19 cm3 × molecule-1 s-1

294

(ozonolysis reaction rate coefficient for the limonene exocyclic double bond), 11, 12, 16 [OH] = 2 × 106

295

molecule × cm-3 and [O3] = 2 × 1012 molecule × cm-3.

296

[OH]= 2 × 10

-14

-9

M, [O3] =2 × 10 M,

25

32

In the aqueous-phase, it was assumed that

and pKa of LA ≈ 4.76. 50

297 298

Figure 8 Calculated LA atmospheric lifetimes in the gas and aqueous phase

299

As shown in Fig. 8, in the gas-phase LA primarily reacts with OH due to very slow ozonolysis

300

of the terminal double bond. 11, 12, 16 The data acquired also indicate that LA lifetime in the aqueous-

301

phase is similar to the estimated lifetime of this compound in the gas-phase. Moreover, since both k1

302

(Fig. S16) and k2 are pH-dependent, reaction (2) may dominate over reaction (1) for LA and related

303

compounds when pH ≥ 4.5 due to dissociation of the terpene-derived carboxylic acids under mildly

304

acidic and basic conditions.

305

carboxylate ion of LA,

306

marine aerosols but not in the acidic aerosol particles (Fig. 8).

307

ozonolysis of monoterepene-derived, unsaturated carboxylic acids are needed.

25

50, 51

Given a reasonable assumption that k2 is factor of 2-4 higher for

reaction (2) can compete with reaction (1) in e.g. clouds, fogs and some 25

Therefore, further studies of the

308

When LWC ≥ 0.05 g × m-3 LA will reside entirely in the aqueous-phase (Fig. S17), therefore

309

reactions (1) and (2) can be relevant under realistic atmospheric conditions. However, in clouds the

14

ACS Paragon Plus Environment

Page 15 of 20

Environmental Science & Technology

310

mechanisms of these reactions and products distribution will be different than in the laboratory

311

reactor.

312

The RO2 radicals that are formed due to reaction (1) will most likely react with the other RO2

313

radicals that are more abundant in clouds rather than self react. Consequently, the distribution of

314

products will depend on the branching ratios of the “mixed” tetroxides. Moreover, inside the water-

315

containing particles, atmospherically abundant carboxylic acids (formic, acetic etc.) 52, 53 will rapidly 54

316

scavenge the C9 SCIs that are produced folloging reaction (2) to yield the relatively stable,55 non-

317

volatile α-acyloxyhydroperoxy aldehydes.

318

likely remain in the particle phase after water evaporation thereby contributing to SOA. 56 Previously,

319

it was demonstrated that even when a low amounts of a non-volatile products are formed (such as

320

methyl-1,2,3-butanetricarboxylic acid, MBTCA from the CPA + OH reaction) the aqSOA relative yield

321

can be as high as 50%.

322

atmosphere. Both reactions (1) and (2) also produce keto-LA, a well-known limonene SOA tracer that

323

was identified in the ambient particles. 58-60

324 325

29

5, 56, 57

These α-acyloxyhydroperoxy aldehydes will most

Thus, reactions (1) and (2) can most likely contribute to

aqSOA

in the

Acknowledgments This

work

was

founded

by

the

Polish

National

Science

Centre:

grant

number

326

2014/13/B/ST4/04500. Sara Jurdana was supported by the EU Foundation for the Development of

327

the Education System Erasmus+ program 2014-2020. We thank the Structural Research Laboratory

328

(SRL) at the Department of Chemistry of University of Warsaw for making LC/MS measurements

329

possible. SRL has been established with financial support from European Regional Development

330

Found in the Sectoral Operational Programme “Improvement of the Competitiveness of Enterprises,

331

years 2004–2005” project no: WPK_1/ 1.4.3./1/2004/72/72/165/2005/U. The study was carried out

332

at the Biological and Chemical Research Centre, University of Warsaw, established within the project

333

co-financed by European Union from the European Regional Development Fund under the

334

Operational Programme Innovative Economy, 2007 – 2013. We thank dr Dagmara Tymecka for LA

335

purification with a semi-preparative HPLC. We thank dr Marcin Wilczek for the

336

measurements. We thank the anonymous reviewers for very helpful, constructive and insightful

337

comments.

338 339

1

H NMR

Supporting Information available, this information is available free of charge via the Internet at http://pubs.acs.org.

340 341

Literature

15

ACS Paragon Plus Environment

Environmental Science & Technology

342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391

Page 16 of 20

1. Odum, J. R.; Hoffmann, T.; Bowman, F.; Collins, D.; Flagan, R. C.; Seinfeld, J. H., Gas/Particle Partitioning and Secondary Organic Aerosol Yields. Environ. Sci. Technol. 1996, 30, (8), 2580-2585. 2. Seinfeld, J. H.; Pankow, J. F., Organic atmospheric particulate material. Annu. Rev. Phys. Chem. 2003, 54, 121-140. 3. Messina, P.; Lathière, J.; Sindelarova, K.; Vuichard, N.; Granier, C.; Ghattas, J.; Cozic, A.; Hauglustaine, D. A., Global biogenic volatile organic compound emissions in the ORCHIDEE and MEGAN models and sensitivity to key parameters. Atmos. Chem. Phys. 2016, 16, (22), 14169-14202. 4. Atkinson, R.; Arey, J., Gas-phase tropospheric chemistry of biogenic volatile organic compounds: a review. Atmos. Environ. 2003, 37, (Supplement 2), 197-219. 5. Hallquist, M.; Wenger, J. C.; Baltensperger, U.; Rudich, Y.; Simpson, D.; Claeys, M.; Dommen, J.; Donahue, N. M.; George, C.; Goldstein, A. H.; Hamilton, J. F.; Herrmann, H.; Hoffmann, T.; Iinuma, Y.; Jang, M.; Jenkin, M. E.; Jimenez, J. L.; Kiendler-Scharr, A.; Maenhaut, W.; McFiggans, G.; Mentel, T. F.; Monod, A.; Prévôt, A. S. H.; Seinfeld, J. H.; Surratt, J. D.; Szmigielski, R.; Wildt, J., The formation, properties and impact of secondary organic aerosol: current and emerging issues. Atmos. Chem. Phys. 2009, 9, (14), 5155-5236. 6. Calvo, A. I.; Alves, C.; Castro, A.; Pont, V.; Vicente, A. M.; Fraile, R., Research on aerosol sources and chemical composition: Past, current and emerging issues. Atmospheric Res. 2013, 120– 121, 1-28. 7. Gilardoni, S.; Massoli, P.; Paglione, M.; Giulianelli, L.; Carbone, C.; Rinaldi, M.; Decesari, S.; Sandrini, S.; Costabile, F.; Gobbi, G. P.; Pietrogrande, M. C.; Visentin, M.; Scotto, F.; Fuzzi, S.; Facchini, M. C., Direct observation of aqueous secondary organic aerosol from biomass-burning emissions. Proc. Natl. Acad. Sci. U.S.A. 2016, 113, (36), 10013-10018. 8. Chung, C. E.; Ramanathan, V.; Kim, D.; Podgorny, I. A., Global anthropogenic aerosol direct forcing derived from satellite and ground-based observations. J. Geophys. Res. 2005, 110, (D24), D24207. 9. Goldstein, A. H.; Galbally, I. E., Known and Unexplored Organic Constituents in the Earth's Atmosphere. Environ. Sci. Technol. 2007, 41, (5), 1514-1521. 10. McNeill, V. F., Aqueous Organic Chemistry in the Atmosphere: Sources and Chemical Processing of Organic Aerosols. Environ. Sci. Technol. 2015, 49, (3), 1237-1244. 11. Jonsson, Å. M.; Hallquist, M.; Ljungström, E., Impact of Humidity on the Ozone Initiated Oxidation of Limonene, Δ3-Carene, and α-Pinene. Environ. Sci. Technol. 2006, 40, (1), 188-194. 12. Leungsakul, S.; Jaoui, M.; Kamens, R. M., Kinetic Mechanism for Predicting Secondary Organic Aerosol Formation from the Reaction of d-Limonene with Ozone. Environ. Sci. Technol. 2005, 39, (24), 9583-9594. 13. Youssefi, S.; Waring, M. S., Indoor transient SOA formation from ozone+α-pinene reactions: Impacts of air exchange and initial product concentrations, and comparison to limonene ozonolysis. Atmos. Environ. 2015, 112, 106-115. 14. Ng, N. L.; Kroll, J. H.; Keywood, M. D.; Bahreini, R.; Varutbangkul, V.; Flagan, R. C.; Seinfeld, J. H.; Lee, A.; Goldstein, A. H., Contribution of First- versus Second-Generation Products to Secondary Organic Aerosols Formed in the Oxidation of Biogenic Hydrocarbons. Environ. Sci. Technol. 2006, 40, (7), 2283-2297. 15. Gallimore, P. J.; Mahon, B. M.; Wragg, F. P. H.; Fuller, S. J.; Giorio, C.; Kourtchev, I.; Kalberer, M., Multiphase composition changes and reactive oxygen species formation during limonene oxidation in the new Cambridge Atmospheric Simulation Chamber (CASC). Atmos. Chem. Phys. 2017, 17, (16), 9853-9868. 16. Zhang, J.; Huff Hartz, K. E.; Pandis, S. N.; Donahue, N. M., Secondary Organic Aerosol Formation from Limonene Ozonolysis:  Homogeneous and Heterogeneous Influences as a Function of NOx. J. Phys. Chem. A. 2006, 110, (38), 11053-11063. 17. Witkowski, B.; Gierczak, T., Characterization of the limonene oxidation products with liquid chromatography coupled to the tandem mass spectrometry. Atmos. Environ. 2017, 154, 297-307.

16

ACS Paragon Plus Environment

Page 17 of 20

392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443

Environmental Science & Technology

18. Nørgaard, A. W.; Vibenholt, A.; Benassi, M.; Clausen, P. A.; Wolkoff, P., Study of OzoneInitiated Limonene Reaction Products by Low Temperature Plasma Ionization Mass Spectrometry. J. Am. Soc. Mass Spectrom. 2013, 24, (7), 1090-1096. 19. Pathak, R. K.; Salo, K.; Emanuelsson, E. U.; Cai, C.; Lutz, A.; Hallquist, Å. M.; Hallquist, M., Influence of Ozone and Radical Chemistry on Limonene Organic Aerosol Production and Thermal Characteristics. Environ. Sci. Technol. 2012, 46, (21), 11660-11669. 20. Kundu, S.; Fisseha, R.; Putman, A. L.; Rahn, T. A.; Mazzoleni, L. R., High molecular weight SOA formation during limonene ozonolysis: insights from ultrahigh-resolution FT-ICR mass spectrometry characterization. Atmos. Chem. Phys. 2012, 12, (12), 5523-5536. 21. Blando, J. D.; Turpin, B. J., Secondary organic aerosol formation in cloud and fog droplets: a literature evaluation of plausibility. Atmos. Environ. 2000, 34, (10), 1623-1632. 22. Ervens, B.; Turpin, B. J.; Weber, R. J., Secondary organic aerosol formation in cloud droplets and aqueous particles (aqSOA): a review of laboratory, field and model studies. Atmos. Chem. Phys. 2011, 11, (21), 11069-11102. 23. Sullivan, A. P.; Hodas, N.; Turpin, B. J.; Skog, K.; Keutsch, F. N.; Gilardoni, S.; Paglione, M.; Rinaldi, M.; Decesari, S.; Facchini, M. C.; Poulain, L.; Herrmann, H.; Wiedensohler, A.; Nemitz, E.; Twigg, M. M.; Collett Jr, J. L., Evidence for ambient dark aqueous SOA formation in the Po Valley, Italy. Atmos. Chem. Phys. 2016, 16, (13), 8095-8108. 24. Burkholder, J. B.; Abbatt, J. P. D.; Barnes, I.; Roberts, J. M.; Melamed, M. L.; Ammann, M.; Bertram, A. K.; Cappa, C. D.; Carlton, A. G.; Carpenter, L. J.; Crowley, J. N.; Dubowski, Y.; George, C.; Heard, D. E.; Herrmann, H.; Keutsch, F. N.; Kroll, J. H.; McNeill, V. F.; Ng, N. L.; Nizkorodov, S. A.; Orlando, J. J.; Percival, C. J.; Picquet-Varrault, B.; Rudich, Y.; Seakins, P. W.; Surratt, J. D.; Tanimoto, H.; Thornton, J. A.; Tong, Z.; Tyndall, G. S.; Wahner, A.; Weschler, C. J.; Wilson, K. R.; Ziemann, P. J., The Essential Role for Laboratory Studies in Atmospheric Chemistry. Environ. Sci. Technol. 2017, 51, (5), 2519-2528. 25. Herrmann, H.; Schaefer, T.; Tilgner, A.; Styler, S. A.; Weller, C.; Teich, M.; Otto, T., Tropospheric Aqueous-Phase Chemistry: Kinetics, Mechanisms, and Its Coupling to a Changing Gas Phase. Chem. Rev. 2015, 115, (10), 4259-4334. 26. Enami, S.; Sakamoto, Y., OH-Radical Oxidation of Surface-Active cis-Pinonic Acid at the Air– Water Interface. J. Phys. Chem. A. 2016, 120, (20), 3578-3587. 27. Lai, C.; Liu, Y.; Ma, J.; Ma, Q.; Chu, B.; He, H., Heterogeneous Kinetics of cis-Pinonic Acid with Hydroxyl Radical under Different Environmental Conditions. J. Phys. Chem. A. 2015, 119, (25), 65836593. 28. Lignell, H.; Epstein, S. A.; Marvin, M. R.; Shemesh, D.; Gerber, B.; Nizkorodov, S., Experimental and Theoretical Study of Aqueous cis-Pinonic Acid Photolysis. J. Phys. Chem. A. 2013, 117, (48), 12930-12945. 29. Aljawhary, D.; Zhao, R.; Lee, A. K. Y.; Wang, C.; Abbatt, J. P. D., Kinetics, Mechanism, and Secondary Organic Aerosol Yield of Aqueous Phase Photo-oxidation of α-Pinene Oxidation Products. J. Phys. Chem. A. 2016, 120, (9), 1395-1407. 30. Chen, J. J.; Griffin, R. J., Modeling secondary organic aerosol formation from oxidation of alpha-pinene, beta-pinene, and d-limonene. Atmos. Environ. 2005, 39, (40), 7731-7744. 31. Henry's law constants were calcualted with HenryWin v3.20 Website: https://www.epa.gov/tsca-screening-tools/epi-suitetm-estimation-program-interface. 32. Percival, C.; McGillen, M., Overview Of Structure-Activity Relationship Methods For Predicting Gas-Phase Rate Coefficients. In Simulation and Assessment of Chemical Processes in a Multiphase Environment; Barnes, I.; Kharytonov, M. M., Eds. Springer Netherlands: Dordrecht 2008; pp 47. 33. Atkinson, R., A structure-activity relationship for the estimation of rate constants for the gasphase reactions of OH radicals with organic compounds. Int J Chem Kinet. 1987, 19, (9), 799-828. 34. Buxton, G. V.; Greenstock, C. L.; Helman, W. P.; Ross, A. B., Critical Review of rate constants for reactions of hydrated electrons, hydrogen atoms and hydroxyl radicals (⋅OH/⋅O− in Aqueous Solution. J. Phys. Chem. Ref. Data. 1988, 17, (2), 513-886. 17

ACS Paragon Plus Environment

Environmental Science & Technology

444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493

Page 18 of 20

35. Beltrán, F. J.; García-Araya, J. F.; Rivas, F. J.; Álvarez, P.; Rodríguez, E., Kinetics Of Competitive Ozonation Of Some Phenolic Compounds Present In Wastewater From Food Processing Industries. Ozone: Sci. Eng. 2000, 22, (2), 167-183. 36. Witkowski, B.; Gierczak, T., cis-Pinonic Acid Oxidation by Hydroxyl Radicals in the Aqueous Phase under Acidic and Basic Conditions: Kinetics and Mechanism. Environ. Sci. Technol. 2017, 51, (17), 9765-9773. 37. Bader, H.; Hoigné, J., Determination of ozone in water by the indigo method. Water Res. 1981, 15, (4), 449-456. 38. Peter, A.; Von Gunten, U., Oxidation Kinetics of Selected Taste and Odor Compounds During Ozonation of Drinking Water. Environ. Sci. Technol. 2007, 41, (2), 626-631. 39. Schöne, L.; Herrmann, H., Kinetic measurements of the reactivity of hydrogen peroxide and ozone towards small atmospherically relevant aldehydes, ketones and organic acids in aqueous solutions. Atmos. Chem. Phys. 2014, 14, (9), 4503-4514. 40. Zhang, X.; Chen, Z.; Wang, H.; He, S.; Huang, D., An important pathway for ozonolysis of alpha-pinene and beta-pinene in aqueous phase and its atmospheric implications. Atmos. Environ. 2009, 43, (29), 4465-4471. 41. Leviss, D. H.; Van Ry, D. A.; Hinrichs, R. Z., Multiphase Ozonolysis of Aqueous α-Terpineol. Environ. Sci. Technol. 2016, 50, (21), 11698-11705. 42. Kostiainen, R.; Kauppila, T. J., Effect of eluent on the ionization process in liquid chromatography–mass spectrometry. J. Chromatogr. A. 2009, 1216, (4), 685-699. 43. Schone, L.; Schindelka, J.; Szeremeta, E.; Schaefer, T.; Hoffmann, D.; Rudzinski, K. J.; Szmigielski, R.; Herrmann, H., Atmospheric aqueous phase radical chemistry of the isoprene oxidation products methacrolein, methyl vinyl ketone, methacrylic acid and acrylic acid - kinetics and product studies. Phys. Chem. Chem. Phys. 2014, 16, (13), 6257-6272. 44. Zhang, X.; Chen, Z. M.; Zhao, Y., Laboratory simulation for the aqueous OH-oxidation of methyl vinyl ketone and methacrolein: significance to the in-cloud SOA production. Atmos. Chem. Phys. 2010, 10, (19), 9551-9561. 45. Atkinson, R., Atmospheric chemistry of VOCs and NOx. Atmos. Environ. 2000, 34, (12), 20632101. 46. Enami, S.; Hoffmann, M. R.; Colussi, A. J., In Situ Mass Spectrometric Detection of Interfacial Intermediates in the Oxidation of RCOOH(aq) by Gas-Phase OH-Radicals. J. Phys. Chem. A. 2014, 118, (23), 4130-4137. 47. Praplan, A. P.; Barmet, P.; Dommen, J.; Baltensperger, U., Cyclobutyl methyl ketone as a model compound for pinonic acid to elucidate oxidation mechanisms. Atmos. Chem. Phys. 2012, 12, (22), 10749-10758. 48. Liu, M. J.; Wiegel, A. A.; Wilson, K. R.; Houle, F. A., Aerosol Fragmentation Driven by Coupling of Acid–Base and Free-Radical Chemistry in the Heterogeneous Oxidation of Aqueous Citric Acid by OH Radicals. J. Phys. Chem. A. 2017, 121, (31), 5856-5870. 49. Wang, H. L.; Huang, D.; Zhang, X.; Zhao, Y.; Chen, Z. M., Understanding the aqueous phase ozonolysis of isoprene: distinct product distribution and mechanism from the gas phase reaction. Atmos. Chem. Phys. 2012, 12, (15), 7187-7198. 50. Marvin v17.28.0 developed by ChemAxon was used to calculate pKa of LA, Website: https://chemaxon.com/products/marvin. 51. Howell, H.; Fisher, G. S., The Dissociation Constants of Some of the Terpene Acids. J. Am. Chem. Soc 1958, 80, (23), 6316-6319. 52. Khare, P.; Kumar, N.; Kumari, K. M.; Srivastava, S. S., Atmospheric formic and acetic acids: An overview. Rev. Geophys. 1999, 37, (2), 227-248. 53. Löflund, M.; Kasper-Giebl, A.; Schuster, B.; Giebl, H.; Hitzenberger, R.; Puxbaum, H., Formic, acetic, oxalic, malonic and succinic acid concentrations and their contribution to organic carbon in cloud water. Atmos. Environ. 2002, 36, (9), 1553-1558.

18

ACS Paragon Plus Environment

Page 19 of 20

494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514

Environmental Science & Technology

54. Tobias, H. J.; Ziemann, P. J., Kinetics of the Gas-Phase Reactions of Alcohols, Aldehydes, Carboxylic Acids, and Water with the C13 Stabilized Criegee Intermediate Formed from Ozonolysis of 1-Tetradecene. J. Phys. Chem. A. 2001, 105, (25), 6129-6135. 55. Witkowski, B.; Gierczak, T., Analysis of α-acyloxyhydroperoxy aldehydes with electrospray ionization–tandem mass spectrometry (ESI-MSn). J. Mass Spectrom. 2013, 48, (1), 79-88. 56. Lee, S.; Kamens, R. M., Particle nucleation from the reaction of α-pinene and O3. Atmos. Environ. 2005, 39, (36), 6822-6832. 57. Kamens, R. M.; Jaoui, M., Modeling Aerosol Formation from α-Pinene + NOx in the Presence of Natural Sunlight Using Gas-Phase Kinetics and Gas-Particle Partitioning Theory. Environ. Sci. Technol. 2001, 35, (7), 1394-1405. 58. Jaoui, M.; Corse, E.; Kleindienst, T. E.; Offenberg, J. H.; Lewandowski, M.; Edney, E. O., Analysis of Secondary Organic Aerosol Compounds from the Photooxidation of d-Limonene in the Presence of NOX and their Detection in Ambient PM2.5. Environ. Sci. Technol. 2006, 40, (12), 38193828. 59. Yasmeen, F.; Szmigielski, R.; Vermeylen, R.; Gómez-González, Y.; Surratt, J. D.; Chan, A. W. H.; Seinfeld, J. H.; Maenhaut, W.; Claeys, M., Mass spectrometric characterization of isomeric terpenoic acids from the oxidation of α-pinene, β-pinene, d-limonene, and Δ3-carene in fine forest aerosol. J. Mass Spectrom. 2011, 46, (4), 425-442. 60. Kourtchev, I.; Fuller, S.; Aalto, J.; Ruuskanen, T. M.; McLeod, M. W.; Maenhaut, W.; Jones, R.; Kulmala, M.; Kalberer, M., Molecular Composition of Boreal Forest Aerosol from Hyytiälä, Finland, Using Ultrahigh Resolution Mass Spectrometry. Environ. Sci. Technol. 2013, 47, (9), 4069-4079.

515

19

ACS Paragon Plus Environment

Environmental Science & Technology

Table of Contents Art 47x26mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 20 of 20