Particle Partitioning ... - ACS Publications

Aug 31, 2017 - Accommodation Coefficients in Environmental Chambers ... from which we quantify the aerosol mass accommodation coefficient (α) as 0.7 ...
0 downloads 0 Views 1MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Article

Direct measurements of gas/particle partitioning and mass accommodation coefficients in environmental chambers Jordan Edward Krechmer, Douglas A. Day, Paul J. Ziemann, and Jose L. Jimenez Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b02144 • Publication Date (Web): 31 Aug 2017 Downloaded from http://pubs.acs.org on September 3, 2017

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 free 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 accessible to all readers and 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 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36

Environmental Science & Technology

Direct measurements of gas/particle partitioning and mass accommodation coefficients in environmental chambers Jordan E. Krechmer%, Douglas A. Day, Paul J. Ziemann, Jose L. Jimenez* Cooperative Institute for Research in Environmental Sciences (CIRES) and Department of Chemistry and Biochemistry, Boulder, Colorado 80309, United States Corresponding author: * Jose-Luis Jimenez University of Colorado; UCB 216, Boulder, CO 80309-0216 Phone: 303-492-3557 Fax: 303-492-1149 [email protected] Keywords: Aerosols, SOA, environmental chambers, gas/particle partitioning Abstract Secondary organic aerosols (SOA) are a major contributor to fine particulate mass and wield substantial influences on the Earth’s climate and human health. Despite extensive research in recent years, many of the fundamental processes of SOA formation and evolution remain poorly understood. Most atmospheric aerosol models use gas/particle equilibrium partitioning theory as a default treatment of gas-aerosol transfer, despite questions about potentially large kinetic effects. We have conducted fundamental SOA formation experiments in a Teflon environmental chamber using a novel method. A simple chemical system produces a very fast burst of lowvolatility gas-phase products, which are competitively taken up by liquid organic seed particles and Teflon chamber walls. Clear changes in the species time evolution with differing amounts of seed allow us to quantify the particle uptake processes. We reproduce gas- and aerosol-phase observations using a kinetic box model, from which we quantify the aerosol mass accommodation coefficient (α) as 0.7 on average, with values near unity especially for low volatility species. α appears to decrease as volatility increases. α has historically been a very difficult parameter to measure with reported values varying over three orders of magnitude. We use the experimentally-constrained model to evaluate the correction factor (Φ) needed for chamber SOA mass yields due to losses of vapors to walls as a function of species volatility and particle condensational sink. Φ ranges from 1 – 4.

1 ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 31

37

Introduction

38

Gas/particle partitioning (GPP) is a fundamental physical process describing the interaction

39

between atmospheric organic gases and particles. Originally treated as an adsorptive process,1,2

40

most atmospheric models currently treat GPP of organic compounds as an absorptive process.3–5

41

In this interpretation, semivolatile organic gases condense into and evaporate from particles until

42

quickly reaching a volatility- and available organic aerosol mass-dependent equilibrium.3,4

43

Semivolatile equilibrium partitioning is often parameterized using a volatility basis set, which

44

lumps the amount of condensable material in the air in multiple bins spaced typically by a decade

45

of the saturation mass concentration of the compound.5 GPP has been incorporated into

46

numerous box, regional, and global models6–10 as a default treatment of organic aerosol

47

formation and evolution.

48

Still, modeling of GPP by rapid absorptive equilibrium remains controversial. There are

49

now several real-time analysis techniques capable of measuring the gas/particle distribution of

50

trace organic compounds at relatively high time resolution in the field.11–14 Many of these

51

measurements have shown substantial disagreements with semivolatile equilibrium partitioning

52

theory and one another, however.15,16 Some recent works have even suggested that semivolatile

53

equilibrium partitioning is not reached under certain conditions due to strong kinetic limitations,

54

and more detailed treatments are required to match atmospheric observations.17–21 GPP model treatment is expected to have a large influence on secondary organic aerosol

55 56

(SOA)

mass

concentrations

predicted

by

atmospheric

models.

Most

models

use

57

parameterizations of the SOA yield, i.e. the mass of SOA formed per unit mass of volatile

58

organic compound (VOC) precursor reacted.22 Almost all experiments quantifying SOA mass

59

yields that have been used in models were performed in Teflon-walled environmental

2 ACS Paragon Plus Environment

Page 3 of 31

Environmental Science & Technology

60

chambers.22–26 However, it has recently been reported that Teflon chambers suffer from losses of

61

semivolatile gaseous compounds to chamber walls.27–30 Wall losses can have strong effects on

62

aerosol mass yields,30 and could also affect conclusions about chemical composition and

63

processes gleaned from chamber experiments. The question even arises: if wall losses are

64

substantial for a certain range of compounds, and fast enough to be on similar time scales as the

65

approach to equilibrium GPP, is it possible to observe and accurately characterize GPP kinetics

66

and equilibrium in environmental chambers?

67

In mathematical representations of the kinetics of GPP, the mass accommodation

68

coefficient (α; also known as the sticking coefficient; and equal to the evaporation coefficient at

69

equilibrium) is a critical parameter that defines the fraction of gas/particle collisions that result in

70

a colliding species being taken up by a particle. α for organic species was historically assumed to

71

be 1, mostly from liquid evaporation measurements.31–33 More recent measurements of mass

72

accommodation coefficients for aerosol-phase compounds34–36 range across two orders of

73

magnitude. A very wide range of values (as low as 0.001) has been used in models and data

74

fitting, sometimes as a tuning parameter which might obscure other model structural

75

limitations.30,37–40 Comparing these literature values is complicated by the different methods used

76

to derive α. Some models34,41 and measurements34,35 conflate particle diffusion kinetics and

77

surface accommodation as one effective mass accommodation coefficient. Other models have

78

attempted to specifically separate mass accommodation and particle-phase diffusion

79

coefficients.36,42,43

80

Given the discrepancies that persist between ambient SOA measurements and

81

models,44,45 the substantial fraction of organic material in the atmosphere which is thought to be

82

semivolatile,46 and the difficulty of quantifying GPP in field studies,16 it is important to closely

3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 31

83

examine in detail the kinetics of approach to equilibrium GPP. . It is also important to further

84

evaluate the impact that losses of vapors to Teflon chamber walls may have in SOA yield

85

measurements. Herein, we conduct fundamental gas/particle partitioning experiments in an

86

environmental chamber with a particular focus on GPP kinetics. A well-characterized simple

87

chemical system is used to produce low-volatility organic compounds very rapidly, which are

88

taken up by liquid organic seed particles and/or the Teflon chamber walls. We use liquid organic

89

seed aerosol specifically to establish an experimental method under which particle-phase kinetic

90

limitations are not present. Both gas-phase products and total aerosol volume and surface area

91

concentrations are continuously monitored. Analysis of these experiments allow us to thoroughly

92

analyze the gas-phase and accommodation kinetics of GPP in the presence of walls for semi- and

93

low-volatility gas-phase compounds. A box model is used to quantify c*and α. We use liquid

94

aerosol particles as seeds, which removes particle-phase diffusion as a limiting parameter and

95

allows us to infer the true value of α. We discuss the implications for quantifying gas/particle

96

partitioning under more complex conditions.

97 98

MATERIALS AND METHODS

99

Environmental chamber experiments

100

A schematic outlining the chamber and instruments used in this work is shown in Figure

101

1. We conducted all of our experiments in the one of the two chambers in the CU Environmental

102

Chamber (CUEC) Facility. The CUEC bag was constructed of FEP Teflon with a maximum

103

volume of 20 m3. The enclosure temperature was maintained at 26 ± 1° C for all experiments.

104

The chambers were filled with clean, dry (< 1% RH) air from two AADCO Model 737-15 clean

105

air generators (Cleaves, OH, USA). The chambers were equipped with an automated flushing

4 ACS Paragon Plus Environment

Page 5 of 31

Environmental Science & Technology

106

system that alternately evacuated and refilled the bag when it was not in use (~5% of full volume

107

every 5 minutes), while maintaining the bag volume at ~95% of maximum volume using

108

differential pressure between the bag and climate-controlled enclosure. Running the automated

109

flushing sequence overnight before each experiment ensured that each experiment started with

110

aerosol concentrations of 0 particles cm-3 as measured with an ultrafine condensation particle

111

counter (UCPC, TSI model 3776 with a diameter cutoff of 2.5 nm). The experimental procedure is similar to the one detailed in Krechmer et al.27 to measure

112 113

vapor wall losses to Teflon walls, but with the addition of seed aerosol. Dioctyl sebacate (DOS)

114

seed was chosen for this set of experiments because it has been shown19 to be reliably liquid

115

under the experimental conditions used here and should not impose particle-phase diffusion

116

limitations on time scales relevant to the experiment. DOS also has a low vapor pressure (2 × 10-

117

6

118

reaction, thus not affecting GPP through particle-phase reactions.

119

Pa)47 and, as a diester, should be unreactive with any of the products formed in the OH˙

First, we injected dioctyl sebacate (DOS; >97%; Sigma-Aldrich; saturation concentration

120

at 298 K (c*) < 0.1 µg m-3)48 liquid seed aerosol into the chamber using an evaporation-

121

condensation generator modeled after Sinclair and Lamer.49 The aerosol generator was

122

configured to seed the chamber to 100 µg m-3 of DOS in ~15 minutes. Particle-phase wall losses

123

occurred on a ~10 hr timescale, but it was not necessary to account for them explicitly in this

124

model since the time-dependent aerosol size distribution measurements were used to constrain

125

the box model. If needed, we diluted the chamber with clean air until the DOS seed surface area

126

reached a desired level. '-condensation generator was tuned to produce aerosol with a lognormal

127

surface area distribution centered at 200 nm. A scanning mobility particle sizer (SMPS)

128

comprised of a TSI model 3080 electrostatic classifier and model 3775 condensation particle

5 ACS Paragon Plus Environment

Environmental Science & Technology

129

counter monitored aerosol number, surface area, and volume concentrations. The SMPS scan

130

was programmed to last 120 s, with 15 s to return to the starting voltage, resulting in a time

131

resolution of 135 s. We provide a summary of different experiments and seed concentrations in

132

Table S1.

133

Page 6 of 31

The chemical system used here was the oxidation of long-chain alcohols by the OH

134

radical under high-NO conditions. This system is the same used by Krechmer et al.27 because the

135

oxidation products are simple and well-characterized (and with consistent chemical functionality

136

in a progression of carbon numbers) as opposed to those of other gaseous precursors such as α-

137

pinene or isoprene.

138

After the seed aerosol concentration stabilized, we injected 50 µL of 1-hexanol, 40 µL of

139

1-octanol, 25 µL of 1-nonanol, 18 µL of 1-decanol, and 12 µL of 1-dodecanol (all >99% purity;

140

Sigma-Aldrich) into the chamber by gently heating the liquids into a clean nitrogen stream (UHP

141

liquid nitrogen evaporate) flowing into the chamber. Injection took ~5 minutes. This resulted in a

142

total VOC concentration in the chamber of 1.4 ppm. We fully mixed the bag contents by turning

143

on a Teflon-coated fan inside the bag for 1 min to mix the contents. Finally, we injected 4 ppm

144

of NO (99%, Matheson Tri-gas) and 4 ppm of methyl nitrite (MeONO, synthesized after the

145

method of Taylor et al.50) into the chamber as pure gases using a small glass bulb and mixed the

146

contents again for 1 min.

147

To initiate chemical production, we turned on the chamber UV black lights (~300-400

148

nm) at 100% intensity (JNO2 = 0.01 s-1) for precisely 10 s via computer-control. MeONO formed

149

OH radicals, which then oxidized the alkanols forming ppt-levels of several oxidized products,27

150

including hydroxynitrates (HN), dihydroxynitrates (DHN), trihydroxynitrates (THN), and

151

carbonyl dihydroxynitrates (CDHN). A list of the compounds, along with their exact masses and

6 ACS Paragon Plus Environment

Page 7 of 31

Environmental Science & Technology

152

estimated saturation concentrations (c*) is provided in Table S2. We did not observe new particle

153

formation. We were able to continuously monitor the latter three types of products with a

154

chemical ionization mass spectrometer (CIMS). This “rapid burst” method, developed

155

previously,27 is critical to these experiments because it produces compounds over a wide range of

156

volatilities over the entire volume of the chamber much faster than the vapors can equilibrate

157

with the particles and walls. After each burst, the system was left to come to equilibrium for at

158

least 1 hr.

159

The amount of aerosol produced by each photooxidation burst comprised an average of

160

7% of the total aerosol volume in the chamber (Figure S1). Therefore, the system can be

161

modeled assuming that the gas-phase compounds are partitioning into pure DOS particles and

162

not into a product compound mixture. The specific molecules studied here are estimated to form

163

~ 20 ppt and account for < ~5% of the SOA. The majority of the aerosol formed is likely to be

164

semivolatile hydroxynitrates, which are not detectable by the nitrate CIMS.27 The commercial

165

iodide CIMS, which would likely detect the majority of aerosol-forming compounds, suffers

166

from a significant time delay for semivolatile products due to gas-wall partitioning in the

167

instrument source.27 This makes it unsuitable for a study such as this where it is critical to

168

measure the gas-phase decays as close to real time as possible.

169

NO3-CIMS

170

A chemical ionization time-of-flight mass spectrometer equipped with a nitrate ion

171

ionization source (NO3-CIMS) monitored gas-phase product concentrations. The instrument51,52

172

and nitrate anion source53 have been described extensively in previous works. We used the

173

instrument in a configuration similar to the one in Krechmer et al.27. Notably here, we placed the

174

instrument inside the CUEC enclosure directly adjacent to the reactor bag. This allowed us to use

7 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 31

175

a 0.6 m long electropolished stainless steel inlet to isothermally bring 10 standard L min-1

176

(defined as 25° C, 1 Atm) of sample air directly from the bag to the instrument without dilution.

177

The short inlet with a sub-sampled center flow minimizes wall losses and results in a residence

178

time of sampled air in the instrument inlet and source of < 4 s. The NO3-CIMS acquired spectra

179

at 1 Hz. We then processed instrument data using Tofware (Tofwerk, AG and Aerodyne, version

180

2.5.8) toolkit within IGOR Pro 6 (Wavemetrics, Inc.) at the same time resolution. Throughout

181

this work we report CIMS results in units of ions s-1 because we are only concerned with relative

182

changes in gas-phase concentration and the CIMS response has previously been shown to be

183

linear.54,55 The NO3-CIMS detected gaseous product molecules as clusters with the nitrate ion,

184

but the NO3- prefix has been removed from reported formulas throughout this paper for clarity.

185

An example mass spectrum of some of the product compounds is provided in Figure S2.

186

Box model

187

A chemical kinetic box model was used to simulate the behavior of gas- and aerosol-

188

phase compounds, accounting for partitioning and evaporation of each oxidized products to walls

189

and particles. A schematic of the model is shown in Figure S3. The model was solved using

190

KinSim v3.24 within IGOR Pro.56

191

The vapor wall loss or condensation rate coefficient, kc was determined using the

192

procedure of Krechmer et al.,27 the results of which are shown in Figure S4. Under these

193

conditions (c*