Kinetics of Limonene Secondary Organic Aerosol Oxidation in the

Sep 12, 2018 - University of Warsaw , Faculty of Chemistry, Al. Żwirki i Wigury 101, 02-089 Warsaw , Poland. Environ. ... Environmental Science & Tec...
0 downloads 0 Views 818KB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Environmental Processes

Kinetics of limonene SOA oxidation in the aqueous phase Bartlomiej Witkowski, Mohammed Al-sharafi, and Tomasz Gierczak Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b02516 • Publication Date (Web): 12 Sep 2018 Downloaded from http://pubs.acs.org on September 13, 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.

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

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

ACS Paragon Plus Environment

Environmental Science & Technology

1

Kinetics of limonene SOA oxidation in the aqueous phase

2

Bartłomiej Witkowski, 1* Mohammed Al-sharafi, 1 Tomasz Gierczak 1 1

3 4

Page 2 of 20

University of Warsaw, Faculty of Chemistry,

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

5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

Keywords: Limonene, secondary organic aerosol, aqueous-phase, relative rates, hydroxyl radicals,

22

ozone

23

*Corresponding author e-mail: [email protected]

1 ACS Paragon Plus Environment

Page 3 of 20

Environmental Science & Technology

24

Abstract

25

Twenty semi-volatile organic compounds that contribute to limonene secondary organic

26

aerosol (SOA) were synthesized in the flow-tube reactor. Kinetics of the aqueous-phase oxidation of

27

the synthesized compounds by hydroxyl radicals (OH) and ozone (O3) were investigated at 298±2 K

28

using the relative rate method. Oxidized organic compounds identified as the major components of

29

limonene SOA were quantified with liquid chromatography coupled to the electrospray ionization

30

and quadrupole tandem mass spectrometry (LC-ESI/MS/MS). The bimolecular rate coefficients

31

measured for the oxidation products of limonene are: kOH = 2-5 × 109 M-1s-1 for saturated and kOH = 1-

32

2 × 1010 M-1s-1 for unsaturated compounds. Ozonolysis reaction bimolecular rate coefficients

33

obtained for the unsaturated compounds in the aqueous phase are between 2-6 × 104 M-1s-1. The

34

results obtained in this work also indicate that oxidation of limonene carboxylic acids by OH was

35

about factor of 2 slower for the carboxylate ions than for the corresponding carboxylic acids while

36

the opposite was true for the ozonolysis. The data acquired provided new insights into kinetics of the

37

limonene SOA processing in the aqueous-phase. Ozonolysis of limonene SOA also increased the

38

concentration of dimers, most likely due to reactions of the stabilized Criegee intermediates with the

39

other, stable products. These results indicate that aqueous-phase oxidation of limonene SOA by OH

40

and O3 will be relevant in clouds, fogs and wet aerosols.

41 42 43 44 45 46 47 48 49 50 51

2 ACS Paragon Plus Environment

Environmental Science & Technology

52

1.

Page 4 of 20

Introduction

53

Secondary organic aerosols (SOA) are produced by oxidation of the volatile organic

54

compounds (VOC) that are abundant in the atmosphere. 1 Because of high impact of aerosols on the

55

climate and human health,

56

atmospheric scientists. However, the incomplete understanding of SOA formation is still a source of

57

significant uncertainty in the atmospheric models. 3 Recently, more attention has been directed to

58

the aqueous-phase processes 4, 5 in order to improve the performance of the traditional atmospheric

59

models that tend to significantly underestimate the amount of SOA in the atmosphere. 6

60

2, 3

SOA sources and formation mechanisms are of great interest to the

Monoterpenes, mainly α and β-pinene as well as limonene are the second largest group of 7

61

VOC in the atmosphere.

62

monoterpene SOA, including functionalized carboxylic acids, carbonyls as well as non-volatile

63

oligomers. 8-13 In fogs and clouds with LWC = 0.3-0.5 g × m-3 these compounds will reside entirely in

64

the aqueous-phase as inferred from the estimated Henry's law solubility coefficients (H) for such

65

molecules.

66

monoterpene-derived compounds can lead to aqSOA. 15-19

67

14

A large number of semi-volatile organic compounds contribute to

More importantly, recent studies also indicate that the aqueous-phase oxidation of

However, despite some recent advantages,

15, 17, 19-21

the current database of the aqueous22

68

phase reaction of the monoterpene oxidation products is still exceedingly sparse.

69

these reactions are still not represented in the current 3D models. 6 Thus, further studies are needed

70

to expand the current understating of the aqueous-phase processing of monoterpene SOAs.

71

Moreover, the kinetic data previously acquired strongly indicate that the aqueous-phase oxidation of

72

semi-volatile organic compounds produced following gas-phase oxidation of α-pinene and limonene

73

should be relevant under realistic atmospheric conditions. 15, 17, 20, 21

74

Consequently,

In this work we direct our attention to the aqueous-phase oxidation of limonene SOA by two 22

75

tropospheric oxidants: ozone (O3) and hydroxyl radicals (OH).

76

reactions of the limonene oxidation products since some of these molecules retain the less reactive,

77

exocyclic double bond of the precursor.

78

products of limonene will be more reactive towards OH and O3 than the oxidation products of α-

79

pinene that are mostly saturated molecules.

80

unsaturated oxidation products of limonene can also yield low-volatility compounds thereby

81

contributing to aqSOA formation. 20 However, to date, the aqueous-phase oxidation of limonene SOA

82

wasn’t investigated.

12, 23

We investigate the aqueous-phase

Consequently, first generation, unsaturated oxidation

12, 17, 20

Following oxidation in the aqueous-phase, the

83

Previously, limononic acid (3-isopropenyl-6-oxoheptanoic acid, LA) was synthesized by

84

photolysis of the commercially available cis-pinonic acid and purified with the semi-preparative liquid 3 ACS Paragon Plus Environment

Page 5 of 20

Environmental Science & Technology

20, 24

85

chromatography.

86

ozonolysis of limonene.

87

products since a separate synthesis and purification of each of these compounds would be very

88

complicated and time-consuming. Afterwards, the oxidation products of limonene (LimSOA) that were

89

produced in the flow-tube reactor were extracted into water and oxidized by OH (reaction 1) and O3

90

(reaction 2) in two batch reactors. The aim of this work was to provide detailed insights into kinetics

91

of limonene SOA aging in clouds, fogs and wet aerosols. Reactions 1 and 2 were investigated with

92

liquid chromatography coupled to the electrospray ionization tandem mass spectrometry (LC-

93

ESI/MS/MS). The chromatographic (LC) separation combined with the selective and sensitive

94

detection method (MS/MS) allowed for the following a behavior of the individual LimSOA during the

95

aqueous-phase oxidation experiment. Consequently, it was possible to investigate the kinetics of

96

reactions (1) and (2) for each LimSOA while working with the mixture of these compounds; SOA

97

sample from the flow-tube reactor.

98

2.

99

In this work, a flow-tube reactor was used to generate SOA by gas-phase 25

Flow-tube reactor was used to generate twenty limonene oxidation

Experimental Materials and reagents are listed in section S1 of the Supplementary Information (SI). LA was

100

synthesized as previously described.20

101

2.1.

102

Synthesis of the limonene oxidation products The flow reactor was used to synthesize the compounds that were later used in the aqueous25

103

phase experiments.

104

mixing plunger where it was mixed with limonene to initiate the reaction. The reaction conditions

105

were: [O3] = 3 ppm, [limonene] = 10 ppm, pressure = 1 atm, T = 298 K (room temperature), relative

106

humidity ≈ 5%, total flow through the tube = 1.3 L/min, Reynolds number ≈ 14, 26 average residence

107

time = 10 min; no OH scavenger was used. Under these conditions both first and second/third

108

generation products of limonene oxidation were formed.

109

(TX40H120WW, Pall) and extracted into a buffered solution of the reference compounds (see Table

110

1) by mechanical agitation to avoid artifacts formation that may be formed during ultrasound

111

extraction. This solution was passed through a 0.22 µm PTFE syringe filter. The pH of the extract was

112

adjusted and checked with pH-meter (HI 221, Hanna Instruments) before and after each aqueous-

113

phase oxidation experiment (see sections 2.3 and 2.4). The aqueous-phase oxidation was carried out

114

immediately after the extraction and the estimated concentrations of the individual LimSOA were

115

several µM.

116

2.2. Relative rate method

Ozone was added to the 1.5 m x 12 cm Pyrex glass tube through a movable

25

SOA was collected on a filter

4 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 20

117

Relative rate method was used to study kinetics of reactions (1) and (2) for LimSOA. Reference

118

compounds used to measure the bimolecular rate coefficients (kOH and kOzone) for reactions (1) and

119

(2) are listed in Table 1.

120

Table 1 List of the reference compounds and rate coefficients, kref, for their reactions with OH

121

and ozone

Reference compounds

kref for oxidation by OH (M−1 s−1) × 10-9 pH =2

pH=5

pH=10

6.9

6.9

-

Phenylalanine

5.7

6.2

a

p-Toluic acid

-

-

Caffeine

Reference compounds

122

9.0 8.0

kref for oxidation by O3 (M−1 s−1) × 10-5 pH =2

pH=7 and 8

Gallic acid

0.97

4.7

Cinnamic acid

0.31b

3.8

Phenylalanine

-

0.19

a

b

calculated using kOH for caffeine (probably pH-independent) calculated using kOzone for gallic acid-see SI

The kref values listed in Table 1 were taken from the literature.

123

22, 30

22, 27-29

124

discrepancies in kozone values are reported for cinnamic acid.

125

104 M−1 s−1 listed in Table 1 for pH=2 was measured in this work (section S3.3).

126

Note that significant

For this reason, the value of 3.1 ×

Reaction (2) was studied at pH =7 and 8 to avoid dissociation of the hydroxyl groups of gallic 31

and also to avoid decomposition of this compound under strongly basic pH conditions. 29 Note

127

acid

128

that pH = 7 and 8 was sufficient to convert the monoterpene acids into carboxylate ions. 32, 33

129 130

Assuming that a given LimSOA and ref are removed only by reacting with a single oxidant, the unknown bimolecular rate coefficient is calculated using eq. (I). 17, 21 Ln 

k  [Lim ] [Ref] Ln 

=

(I) [Lim ] k  [Ref]

131

[LimSOA] and [Ref] are the initial (0) and intermediate (t) concentrations, kLimSOA and kref are the

132

bimolecular rate coefficients for LimSOA and the reference compound. Note that ozone reacts with

133

gallic and cinnamic acids as well as with phenylalanine in 1:1 molar ratio. 29, 30, 34

134

2.3. Oxidation by hydroxyl radicals

5 ACS Paragon Plus Environment

Page 7 of 20

Environmental Science & Technology

135

The photooxidation experiments were carried out as previously described.

20, 21

pH of the

136

filter extract was adjusted by adding small amount of HCl (pH=2), NaOH (pH=10) or a diluted

137

phosphate buffer (pH=5). The filtered reaction mixture (ca. 15 ml) was placed in a small Pyrex glass

138

bottle and OH was generated by H2O2 (15 mM) photolysis with a 310 W lamp (UVAHAND 250 GS

139

H1/BL, Honle) equipped with a Pyrex glass UVA filter. The solution was mixed with a magnetic stirrer.

140

Aliquots of the reaction mixture (100 µl) were quenched with the catalase solution to stabilize the

141

sample before the LC/MS analysis.

142

2.4. Oxidation by ozone

143

Aqueous-phase ozonolysis was carried in a miniaturized bubble column reactor shown in Fig.

144

S3. The reaction mixture (ca. 15 ml) consisted of limonene SOA solution in DI water with the addition

145

of H3PO4 (pH=2) or H3PO4/Na3PO4 buffer (pH=7) and (pH=8) and ca. 15 mM of t-butyl alcohol (OH

146

scavenger). 35 Ozone was generated by photolysis of pure O2 using a UV generator (UVP SOG-2, 3ppm

147

at 1 L/min version). O2/O3 mixture was bubbled though the solution at a flow rate of 10 ml/min. 100

148

µl aliquots were sampled from the bubbler and incubated with 50 µl of the quenching buffer: 2 µM

149

of indigotrisulfonate and a solution of 0.001 mg/ml of catalase to stabilize the sample. Afterwards, 30

150

µl of ACN was added and the sample was analyzed with LC/MS.

151

2.5. Liquid chromatography coupled to the mass spectrometry

152

Chromatographic analyses were carried out with LC20A liquid chromatograph (Shimadzu) 25

153

coupled to the QTRAP 3200 (AB Sciex) triple quadrupole mass spectrometer.

154

C18 column (100 mm × 2.1 mm, 3 µm, 100 Å) with the security guard cartridge with a 2 mm ID C18

155

pre-column was used for analytes separation. Eluent A was 0.03 % formic acid solution in water (pH =

156

2.8) and eluent B was acetonitrile (ACN), flow rate of the mobile phase was 0.2 mL/min. Gradient

157

elution programs are provided in section S2.

Luna (Phenomenex)

158

The mass spectrometer was equipped with the electrospray (ESI) ion source and was operating

159

in the multiple reaction monitoring (MRM) and scan modes: 50-700 m/z. The ESI conditions were as

160

follows: spray voltage: 4.5 kV (negative ionization mode) and 5.5 kV (positive ionization mode), N2

161

was curtain (3 × 105 Pa), auxiliary (3 × 105 Pa) and collision gas. The MRMs for the LimSOA were

162

published previously. 25 MRM conditions and the sample chromatogram of the reference compounds

163

are provided in section S2, Fig. S1 and S2 and Table S1.

164

2.6. Control experiments and uncertainty

6 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 20

165

Validation of the bubble reactor included: calculations of the ozone mass transfer through the

166

gas-liquid interface to ensure that these experiments were carried out in the kinetic regime (section

167

S3.2) and comparison of the kOzone values obtained for limononic and cinnamic acid under static and

168

dynamic conditions (section S3.3). Control experiments were carried out by turning off the ozone

169

generator to check if the analytes are purged from the solution by the gas flow – Fig. S4-S6.

170

LimSOA and the reference compounds decomposition in the aqueous solution of H2O2 and due to

171

UV-Vis irradiation was checked by carrying out appropriate control experiments – Fig. S8-S10. Some

172

repartitioning of LimSOA into the gas-phase was observed under acidic pH as shown in Fig. S8. To

173

correct for these non-OH losses a control experiment was always carried out simultaneously with the

174

photooxidation as previously described.

175

using a smaller bottle with a reduced headspace volume.

20

However, it was possible to minimize this problem by

176

The uncertainties of the kOH and kOzone values obtained depends on the: precision of the

177

kLimSOA/kref (slope of the linear fits to the experimental data), uncertainty of the kref values that are

178

reported in the literature and also uncertainty of the LC/MS measurements. All of these factors were

179

taken into account and the cumulative uncertainties obtained for kOH and kOzone values were between

180

10-20%.

181

3.

182

3.1. Kinetics of the aqueous-phase oxidation of limonene SOA

Results and discussion

183

In our previous study, 25 limonene SOA composition was investigated with LC-ESI/MS using both

184

low and high-resolution mass spectrometers. The high-resolution fragmentation spectra together

185

with the data reported in the literature was used for structural identification of the semi-volatile

186

organic compounds that contribute to limonene SOA. These tentatively identified structures of the

187

individual LimSOA are shown in Table 2. Detailed discussion regarding structural identification and

188

possible formation pathways of the compounds listed in Table 2 is provided elsewhere. 25

25

As

189

previously demonstrated,

dimers can be unambiguously distinguished from the “monomers” due

190

to chromatographic separation of these higher and lower-MW molecules before introducing the

191

sample into the ESI ion-source of the mass spectrometer.

192

Sample relative kinetic plots for limonic acid (compound 5) and dimer identified as ester of

193

limonic acid and 7-hydroxylimononic (compound 18) are shown in Fig. 1; temporal profiles for the

194

compounds 1-20 and the reference compounds are provided in Fig. S11-S16.

7 ACS Paragon Plus Environment

Page 9 of 20

Environmental Science & Technology

195 196 197

Figure 1 The plots for relative kinetic experiments for the selected LimSOA, red lines are linear fits to the experimental data (blue squares)

198

As shown in Fig. 1, the relative kinetic plots were linear. Therefore the kOH and kOzone values were

199

calculated under different pH conditions using the experimental data acquired. Similar results were

200

obtained for the rest of the LimSOA using the reference compounds listed in Table 1. Average kOH and

201

kOzone values that were calculated with eq. I are listed in Table 2.

8 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 20

Table 2 Bimolecular rate coefficients measured for the LimSOA

202

Compound

Namea

number

1

Structureb

MRM (Q1/Q3

kOH (M−1 s−1) × 10-10

m/z, Da)

Limononic acid

183/139

kOzone (M−1 s−1) × 10-4

Neutral

Carboxylate

Neutral

Carboxylate

1.3 ± 0.1

0.54 ± 0.1

2.6 ± 0.2

5.7 ± 0.3

0.20

0.21

(SAR)

(SAR)

O OH

O

2

Isomers of keto-limononic acid

185/115 O

3

OH

Stable

O O

4

Keto-limononic acid

5

Limonic acid

Increase O

185/141

OH

1.3 ± 0.2

0.39 ± 0.06

0.11

0.17

(SAR)

(SAR)

1.3 ± 0.2

0.50 ± 0.10

2.3 ± 0.2

4.6 ± 0.3

O

OH

6

Keto-limonic acid

O

187/143

OH

O O OH

7

Keto-limonic acid isomer

187/115

O

Increase

OH O HO OH HO O O

O

O

8

Dicarbonyl derivative of limononic acid

197/153

O O

4.0 ± 0.2

6.4 ± 0.3

O OH

9

ACS Paragon Plus Environment

Page 11 of 20

Environmental Science & Technology

9

Dicarbonyl derivative of

197/153

O O

limononic acid

10

1.1 ± 0.1

1.1 ± 0.1

4.1 ± 0.2

6.4 ± 0.3

1.2 ± 0.3

0.8 ± 0.1

2.9 ± 0.3

5.3 ± 0.4

0.35

0.48 ± 0.05

O OH

OH

7-hydroxy limononic acid

O

O

OH

11

199/181 Carbonyl-substituted keto-

O

O O

12 13

O

O

O

O

OH

O

limononic acids

Hydroxy keto – keto-limononic

O

215/97

acid

(SAR)

OH

O OH O

Stable

0.20 ± 0.05

0.81 ± 0.1

0.27

0.30

(SAR)

(SAR)

1.3 ± 0.2

0.69 ± 0.14

3.2 ± 0.3

7.1 ± 0.4

2.2 ± 0.4

0.69 ± 0.06

10 ± 1

11 ± 1

1.2 ± 0.2

0.67 ± 0.2

HO

14

Dihydroperoxy limononic acid

O

233/183

OH O

OH O OH

15

Aldol reaction producs or a

337/319

O

O O

hemiacetal

OH O O O

HO OH O

16

339/115

O

O

Esters/aldol reaction products

O

Increase

OH O OH

17

339/115

O

O

1.3 ± 0.1

0.90 ± 0.1

OH O O O

10

ACS Paragon Plus Environment

Environmental Science & Technology

18

367/185

Page 12 of 20

1.2 ± 0.2

0.70 ± 0.04

3.8 ± 0.5

5.9 ± 0.7

1.6 ± 0.4

1.2 ± 0.2

10 ± 2

10 ± 1

1.4 ± 0.2

1.2 ± 0.1

5.1 ± 1

6.0 ± 1

O OH

O

19 20

Esters of limonic and 7hydroxylimononic acid

O

HO

367/185

O

O OH

367/185

O

O O

HO

O O

203

a

tentatively identified in our previous work and/or based on the literature data available b If known or one of the possible isomers– see below

204

11

ACS Paragon Plus Environment

Page 13 of 20

205

Environmental Science & Technology

As expected, the kOH values that are listed in Table 2 are within diffusion limit 22

36

due to high

206

reactivity of unsaturated compounds towards OH.

207

good agreement with the kinetic data reported for the ozonolysis of unsaturated compounds with

208

the terminal C=C bonds. 22, 35, 37 The pH-dependence of the kOH and kOzone is discussed in section 3.2.

The kOzone values listed in Table 2 are also in a

209

Concentrations of compounds 4, 6 and 7 increased during reactions (1) and (2) since their

210

formation (oxidation of the first generation products) was faster than decomposition (oxidation of

211

saturated molecules by OH). Therefore, compounds 4, 6 and 7 were identified as second-generation

212

products (with the terminal C=C bonds converted into C=O). This conclusion is in a good agreement

213

with the previously proposed structures of these compounds. 9, 25, 38, 39

214

Concentrations of compounds 2, 3, 11-13 decreased during reactions (1) and (2) in some

215

experiments but were stable or had a local maximum (initial increase followed by decomposition) in

216

other as shown in Fig. S11-12 and S14-15. Therefore, compounds 2, 3, 11-13 were identified as the

217

third-generation products. Previously, compounds 11-13 were incorrectly identified as the first

218

generation products and their revised structures are presented in Table 2. A detailed discussion

219

about elucidation of the structures of compounds 11-13 is provided in section S6.

220

When it was possible (concentration decreased exponentially in some but not in all experiments),

221

kOH values for the saturated LimSOA were estimated with eq. (I). The rest of kOH values for the

222

saturated LimSOA was estimated with the structure-activity relationship (SAR) parameters.

223

reactivity of the second and third generation products towards OH is most likely caused by the

224

deactivating effects of the carbonyl and carboxylic groups 49, 50 and also by the relatively low number

225

of aliphatic hydrogen atoms in these molecules. 40, 41

226

40, 41

Low

Ozonolysis of the saturated compounds is too slow to compete with their oxidation by OH,

22

227

therefore, in this case, we did not estimate the kOzone values for the second and third generation

228

products.

229

Reaction (1) also leads to the degradation of dimers, which is in a good agreement with the 15

230

results reported by Zhao et al.

231

neither OH nor O3 is selective towards dimers. More importantly, compounds 16 and 17 were

232

generated during reaction (2) as listed in Table 2 and shown in Fig. S14-S15. In the aqueous solution,

233

the Criegee intermediates are efficiently stabilized 47, 57 thus they can react with the stable molecules

234

(mainly carboxylic acids) to generate higher-MW products. 42, 43 Note that aqueous-phase ozonolysis

235

can dominate over oxidation by OH as a result of the pH-dependence of kOH and kOzone values – see

236

section 3.3.

However, the kinetic data summarized in Table 2 indicate that

237 12 ACS Paragon Plus Environment

Environmental Science & Technology

238

Page 14 of 20

3.2. pH-dependence of kOH and kOzone

239

For the LimSOA that are shown in Table 2, the carboxylic groups are most likely not adjacent to

240

the terminal double bond. Consequently, the rate coefficient of reaction (2) is only moderately

241

enhanced for the carboxylate ions due to the electrophilic character of ozone, as previously

242

concluded.

243

literature database for the carboxylic acid/carboxylate equilibrium impact on the reaction (1)

244

mechanism, especially for the unsaturated carboxylic acids. 22, 27, 45

22, 44

The pH-dependence of the kOH is however difficult to explain due to very limited

245

Kinetics of reaction (1) were also studied at pH=5 (see Table 1) to compare the

246

experimentally obtained and predicted kOH values. The kOH values listed in Table 2 were used to

247

calculate the cumulative rate coefficients at the intermediate pH with eq. (II) and eq. (III). 17 α=!

[H # ] ' (II) [H # ] + Ka

k ( )**+, - = k ( × α + k / × (1 − α) (III) 248

α –fraction of the non-dissociated molecules: AH

249

Ka- estimated carboxylic acid/carboxylate equilibrium constant, 10-pKa - see Table 3

250

[H+] – concentration of the H+ ions (M), 10-pH

251

k OH AH = rate coefficients measured for the protonated acids (AH)

252

k OH A- = rate coefficient measured for the carboxylate ions (A-)

253

For the reaction (1) the calculated pH-dependence of the kOH values were compared with the

254

experimental data - the kOH values measured at pH=5 - as listed in Table 3. Only compounds with the

255

kOH values obtained experimentally (first generation products - see Table 2) are included in Table 3.

256

Table 3 Estimated and experimentally obtained kOH values at pH=5

Compound number pKaa

kOH (M−1 s−1) × 10-9 Calculated Experimental

1

4.76

8.2

6.2 ± 1.0

5

4.17

6.0

5.2 ± 1.0

4.91

8

4.49

6.9

10.0 ± 1.4

9

4.49

9.9

11.0 ± 1.0 13

ACS Paragon Plus Environment

Page 15 of 20

Environmental Science & Technology

10

4.10

14

4.35

b

7.3 ± 1.0

8.0

6.6 ± 1.1

-

15

16 17

4.3

18

8.8 ± 1.2 7.6

5.8 ± 0.8

7.6

5.4 ± 0.9

5.4

6.7 ± 1.0

19

3.90

9.0

7.7 ± 1.1

20

4.50

7.0

6.1 ± 0.8

a

257

8.4

b

estimated with Marvin software Not an acid

258

As listed in Table 3, the measured and predicted kOH values at pH=5 are in a reasonably good

259

agreement. These results indicate that eq. II and III can be used to calculate the kOH values between

260

pH 2 and 10 using the experimental data acquired. Consequently, the rate coefficients obtained were

261

used to estimate the atmospheric lifetimes of LimSOA as a function of pH of the aqueous medium as

262

presented in section 3.3.

263

3.3. Atmospheric implications

264

Average lifetimes calculated for saturated and unsaturated LimSOA are presented in Fig. 2.

265

Values of kOH and kOzone under intermediate pH conditions were estimated as described in section 3.2.

266

It was assumed that only unsaturated compounds will react with O3 (kOzone for the terminal C=C bond

267

≈ 5 × 10-19 cm3 × molecule-1 s-1 in the gas-phase).

268

estimate kOH values for the reaction of LimSOA with OH – see Table S4. 23, 46, 47

269

concentrations were assumed:

270

molecule × cm-3 and 2 × 10-9 M. 22, 32, 50

23, 46, 47

Gas-phase SAR parameters were used to

6

48, 49

The following oxidants

-3

[OH] ≈ 2 × 10 molecule × cm and 2 × 10-14 M, [O3] = 2 × 1012

271 14 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 20

272

Figure 2 The average lifetimes for the saturated and unsaturated LimSOA as a function of pH,

273

uncertainty bars represent the data range, data points are the averaged results for saturated,

274

unsaturated compounds and dimers (na rysunku słowo oxidation jest napisane niepoprawnie

275

Estimated lifetimes for unsaturated LimSOA (as well as for the dimers) due to oxidation by OH in

276

the gas and aqueous-phase are of similar order as shown in Fig. 2. The unsaturated oxidation

277

products of limonene can also react with O3 in addition to the reaction with OH whereas aqueous-

278

phase ozonolysis is not important for the saturated molecules, as expected.

279

unsaturated limonene acids can also compete with their oxidation by OH under mildly acidic and

280

basic pH-conditions due to higher reactivity of the carboxylate ions towards ozone as discussed in

281

section 3.2.

22

Ozonolysis of the

282

The kinetic data acquired, strongly indicate that the reactions studied here can occur in clouds

283

and fogs which agrees well with the previously published data about the aqueous-phase processing

284

of monoterpene SOA.

285

present entirely in the aqueous-phase as inferred from their H values. 14 Consequently, whether or

286

not the LimSOA will be oxidized in aqueous-phase depends on the partitioning of the individual

287

compounds contributing to limonene SOA.

17, 20, 24

When LWC ≈ 0.3-0.5 g × m-3 the molecules listed in Table 2 will be

288

As previously demonstrated, aside from limonene, gas-phase oxidation of the other

289

atmospherically abundant terpenes yields semi-volatile products and these molecules will often

290

retain the less reactive double bond of the parent hydrocarbon. 51-53 Results presented here strongly

291

indicate that such molecules will be highly susceptible to the aqueous-phase oxidation by OH and O3

292

under realistic atmospheric conditions. Aqueous-phase ozonolysis of the compounds with the

293

terminal double bonds appears to be less important than their oxidation by OH. However,

294

generation of the carboxylate ions can sufficiently enhance the rate of aqueous-phase ozonolysis

295

thus it should still be considered as the potential oxidation mechanisms of the unsaturated terpene

296

acids.

297

Both types of reactions that were studied here can yield low-volatility products as previously

298

demonstrated, 17, 19-21 thereby potentially contributing to aqSOA formation under certain conditions.

299

However, it is very difficult to speculate about the aqSOA yields following aqueous-phase oxidation of

300

limonene SOA using the data obtained here. For this reason, a more quantitative approach is

301

necessary to estimate the branching ratios of the unstable by-products as well as yields of highly

302

oxidized products of reactions (1) and (2).

303

Supplementary data

15 ACS Paragon Plus Environment

Page 17 of 20

304

Environmental Science & Technology

Supporting Information available, this information is available free of charge via the Internet at

305

http://pubs.acs.org.

306

Acknowledgments

307

This project was founded by the polish National Science Centre: grant number

308

2014/13/B/ST4/04500. We thank the Structural Research Laboratory (SRL) at the Department of

309

Chemistry of University of Warsaw for the LC/MS measurements. SRL was established with financial

310

support

311

1.4.3./1/2004/72/72/165/2005/U. We thank dr Dagmara Tymecka for LA purification with a semi-

312

preparative HPLC. The study was carried out at the Biological and Chemical Research Centre,

313

University of Warsaw, established within the project co-financed by European Union from the

314

European Regional Development Fund under the Operational Programme Innovative Economy, 2007

315

– 2013. We thank the anonymous reviewers for very helpful, constructive and insightful comments

316

and suggestions.

317

References

318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345

1. 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. 2. Poschl, U., Atmospheric aerosols: Composition, transformation, climate and health effects. Angew. Chem.-Int. Edit. 2005, 44, (46), 7520-7540. 3. Carslaw, K. S.; Lee, L. A.; Reddington, C. L.; Pringle, K. J.; Rap, A.; Forster, P. M.; Mann, G. W.; Spracklen, D. V.; Woodhouse, M. T.; Regayre, L. A.; Pierce, J. R., Large contribution of natural aerosols to uncertainty in indirect forcing. Nature 2013, 503, 67. 4. Ervens, B., Modeling the Processing of Aerosol and Trace Gases in Clouds and Fogs. Chem. Rev 2015, 115, (10), 4157-4198. 5. McNeill, V. F., Aqueous Organic Chemistry in the Atmosphere: Sources and Chemical Processing of Organic Aerosols. Environ. Sci. Technol 2015, 49, (3), 1237-1244. 6. Carlton, A. G.; Turpin, B. J.; Altieri, K. E.; Seitzinger, S. P.; Mathur, R.; Roselle, S. J.; Weber, R. J., CMAQ Model Performance Enhanced When In-Cloud Secondary Organic Aerosol is Included: Comparisons of Organic Carbon Predictions with Measurements. Environ. Sci. Technol 2008, 42, (23), 8798-8802. 7. Sindelarova, K.; Granier, C.; Bouarar, I.; Guenther, A.; Tilmes, S.; Stavrakou, T.; Müller, J. F.; Kuhn, U.; Stefani, P.; Knorr, W., Global data set of biogenic VOC emissions calculated by the MEGAN model over the last 30 years. Atmos. Chem. Phys 2014, 14, (17), 9317-9341. 8. Bateman, A. P.; Nizkorodov, S. A.; Laskin, J.; Laskin, A., Time-resolved molecular characterization of limonene/ozone aerosol using high-resolution electrospray ionization mass spectrometry. Phys. Chem. Chem. Phys. 2009, 11, (36), 7931-7942. 9. Walser, M. L.; Desyaterik, Y.; Laskin, J.; Laskin, A.; Nizkorodov, S. A., High-resolution mass spectrometric analysis of secondary organic aerosol produced by ozonation of limonene. Phys. Chem. Chem. Phys 2008, 10, (7), 1009-1022.

from

European

Regional

Development

Found,

project

no:

WPK_1/

16 ACS Paragon Plus Environment

Environmental Science & Technology

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 392 393 394 395 396

Page 18 of 20

10. Farhat, Y.; Rafal, S.; Reinhilde, V.; Yadian, G. G.; D., S. J.; H., C. A. W.; H., S. J.; Willy, M.; Magda, C., 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. 11. 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. 12. Fang, W.; Gong, L.; Sheng, L., Online analysis of secondary organic aerosols from OH-initiated photooxidation and ozonolysis of α-pinene, β-pinene, Δ3-carene and d-limonene by thermal desorption–photoionisation aerosol mass spectrometry. Environ. Chem. 2017, 14, (2), 75-90. 13. Jianzhen, Y.; J., G. R.; R., C. D.; C., F. R.; H., S. J.; Pierrette, B., Observation of gaseous and particulate products of monoterpene oxidation in forest atmospheres. Geophys Res Lett. 1999, 26, (8), 1145-1148. 14. Brook, R. D.; Rajagopalan, S.; Pope, C. A.; Brook, J. R.; Bhatnagar, A.; Diez-Roux, A. V.; Holguin, F.; Hong, Y.; Luepker, R. V.; Mittleman, M. A.; Peters, A.; Siscovick, D.; Smith, S. C.; Whitsel, L.; Kaufman, J. D., Particulate Matter Air Pollution and Cardiovascular Disease. Circulation. 2010, 121, (21), 2331. 15. Zhao, R.; Aljawhary, D.; Lee, A. K. Y.; Abbatt, J. P. D., Rapid Aqueous-Phase Photooxidation of Dimers in the α-Pinene Secondary Organic Aerosol. Environ. Sci. Technol. Lett. 2017, 4, (6), 205-210. 16. Cortés, D. A.; Elrod, M. J., Kinetics of the Aqueous Phase Reactions of Atmospherically Relevant Monoterpene Epoxides. J. Phys. Chem. A 2017, 121, (48), 9297-9305. 17. 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. 18. 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. 19. 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. 20. Witkowski, B.; Jurdana, S.; Gierczak, T., Limononic Acid Oxidation by Hydroxyl Radicals and Ozone in the Aqueous Phase. Environ. Sci. Technol 2018, 52, (6), 3402-3411. 21. 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. 22. 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. 23. 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. 24. 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. 25. 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. 26. Young, D. F.; Munson, B. R.; Okiishi, T. H.; Huebsch, W. W., A Brief Introduction to Fluid Mechanics. John Wiley & Sons: 2011; p 523. 27. 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

Page 19 of 20

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 444 445 446 447

Environmental Science & Technology

28. Sharma, V. K.; Graham, N. J. D., Oxidation of Amino Acids, Peptides and Proteins by Ozone: A Review. Ozone Sci. Eng. 2010, 32, (2), 81-90. 29. 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. 30. Leitzke, A.; Reisz, E.; Flyunt, R.; von Sonntag, C., The reactions of ozone with cinnamic acids: formation and decay of 2-hydroperoxy-2-hydroxyacetic acid. J. Chem. Soc., Perkin Trans. 2. 2001, (5), 793-797. 31. Erdemgil, F. Z.; Şanli, S.; Şanli, N.; Özkan, G.; Barbosa, J.; Guiteras, J.; Beltrán, J. L., Determination of pKa values of some hydroxylated benzoic acids in methanol–water binary mixtures by LC methodology and potentiometry. Talanta 2007, 72, (2), 489-496. 32. Pilewskie, P., Aerosols heat up. Nature 2007, 448, 541. 33. Howell, H.; Fisher, G. S., The Dissociation Constants of Some of the Terpene Acids. J. Am. Chem. Soc 1958, 80, (23), 6316-6319. 34. Pryor, W. A.; Giamalva, D. H.; Church, D. F., Kinetics of ozonation. 2. Amino acids and model compounds in water and comparisons to rates in nonpolar solvents. J. Am. Chem. Soc 1984, 106, (23), 7094-7100. 35. 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. 36. 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. 37. Chen, Z. M.; Wang, H. L.; Zhu, L. H.; Wang, C. X.; Jie, C. Y.; Hua, W., Aqueous-phase ozonolysis of methacrolein and methyl vinyl ketone: a potentially important source of atmospheric aqueous oxidants. Atmos. Chem. Phys 2008, 8, (8), 2255-2265. 38. 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. nviron. Sci. Technol 2012, 46, (21), 11660-11669. 39. Walser, M. L.; Park, J.; Gomez, A. L.; Russell, A. R.; Nizkorodov, S. A., Photochemical Aging of Secondary Organic Aerosol Particles Generated from the Oxidation of d-Limonene. J. Phys. Chem. A 2007, 111, (10), 1907-1913. 40. Doussin, J. F.; Monod, A., Structure–activity relationship for the estimation of OH-oxidation rate constants of carbonyl compounds in the aqueous phase. Atmos. Chem. Phys 2013, 13, (23), 11625-11641. 41. Monod, A.; Doussin, J. F., Structure-activity relationship for the estimation of OH-oxidation rate constants of aliphatic organic compounds in the aqueous phase: alkanes, alcohols, organic acids and bases. Atmos. Environ. 2008, 42, (33), 7611-7622. 42. 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. 43. Bartłomiej, W.; Tomasz, G., Analysis of α-acyloxyhydroperoxy aldehydes with electrospray ionization–tandem mass spectrometry (ESI-MSn). J. Mass Spectrom 2013, 48, (1), 79-88. 44. Leitzke, A.; Sonntag, C. v., Ozonolysis of Unsaturated Acids in Aqueous Solution: Acrylic, Methacrylic, Maleic, Fumaric and Muconic Acids. Ozone Sci. Eng. 2009, 31, (4), 301-308. 45. Ervens, B.; Gligorovski, S.; Herrmann, H., Temperature-dependent rate constants for hydroxyl radical reactions with organic compounds in aqueous solutions. Phys. Chem. Chem. Phys. 2003, 5, (9), 1811-1824. 46. 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. 18 ACS Paragon Plus Environment

Environmental Science & Technology

448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467

Page 20 of 20

47. 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. 48. 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. 49. 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. 50. Percival, C.; McGillen, M. In Overview Of Structure-Activity Relationship Methods For Predicting Gas-Phase Rate Coefficients, Dordrecht, 2008; Springer Netherlands: Dordrecht, 2008; pp 47-59. 51. 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. 52. Mohammed, J.; Michael, L.; E., K. T.; H., O. J.; O., E. E., β-caryophyllinic acid: An atmospheric tracer for β-caryophyllene secondary organic aerosol. Geophys. Res. Lett. 2007, 34, (5). 53. Li, Y. J.; Chen, Q.; Guzman, M. I.; Chan, C. K.; Martin, S. T., Second-generation products contribute substantially to the particle-phase organic material produced by β-caryophyllene ozonolysis. Atmos. Chem. Phys 2011, 11, (1), 121-132.

468

19 ACS Paragon Plus Environment