Subscriber access provided by NEW YORK MED COLL
Article
A Quantum Chemical Study on •Cl-initiated Atmospheric Degradation of Monoethanolamine Hong-bin Xie, Fangfang Ma, Yuanfang Wang, Ning He, Qi Yu, and Jingwen Chen Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b03324 • Publication Date (Web): 23 Oct 2015 Downloaded from http://pubs.acs.org on October 24, 2015
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 28
Environmental Science & Technology
1
A Quantum Chemical Study on •Cl-initiated Atmospheric
2
Degradation of Monoethanolamine
3
Hong-Bin Xie†, Fangfang Ma†, Yuanfang Wang†, Ning He‡, Qi Yu† and Jingwen
4
Chen†*
5
†
6
Education), School of Environmental Science and Technology, Dalian University of
7
Technology, Dalian 116024, China
8
‡
9
116024, China
10
Key Laboratory of Industrial Ecology and Environmental Engineering (Ministry of
State Key Laboratory of Fine Chemicals, Dalian University of Technology, Dalian
Table of Contents (TOC)
11 12
1
ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 28
13
ABSTRACT: Recent findings on the formation of •Cl in continental urban areas
14
necessitate the consideration of •Cl initiated degradation when assessing the fate of
15
volatile organic pollutants. Monoethanolamine (MEA) is considered as a potential
16
atmospheric pollutant since it is a benchmark and widely utilized solvent in a leading
17
CO2 capture technology. Especially, •Cl may have specific interactions with the
18
N-atom of MEA, which could make the MEA+•Cl reaction have different pathways
19
and products from those of the MEA+•OH reaction. Hence, •Cl initiated reactions
20
with
21
[CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p)] and kinetics modeling. Results show
22
that the overall rate constant for •Cl initiated H-abstraction of MEA is 5 times faster
23
than that initiated by •OH, and the tropospheric lifetimes of MEA will be
24
overestimated by 6%-46% when assuming that [•Cl]/[•OH] = 1%-10% if the role of
25
•Cl is ignored. The MEA+•Cl reaction exclusively produces MEA-N that finally
26
transforms into several products including mutagenic nitramine and carcinogenic
27
nitrosamine via further reactions with O2/NOx, and the contribution of •Cl to their
28
formation is about 25%-250% of that of •OH. Thus, it is necessary to consider •Cl
29
initiated tropospheric degradation of MEA for its risk assessment.
30
INTRODUCTION
MEA
were
investigated
by
a
quantum
chemical
method
31
Chlorine atoms (•Cl) are historically regarded to be produced primarily from
32
heterogeneous reaction cycles involving sea salt.1-4 A significant production of the •Cl
33
precursor ClNO2 via the reaction of NOx with chloride was observed at night or in the
34
early morning near Boulder, Colorado, an urban location in the middle of North
35
America.5 Accordingly, it was estimated that a total annual ClNO2 production rate for
36
the contiguous US lies in the range of 3.2-8.2 Tg yr-1, providing a photolytic •Cl
37
source of 1.4-3.6 Tg Cl yr-1. This US ClNO2 source is far larger than the first global 2
ACS Paragon Plus Environment
Page 3 of 28
Environmental Science & Technology
38
estimate of 0.06 Tg Cl yr-1, and is similar to the recent 3.2 Tg Cl yr-1 estimated for
39
global coastal and marine regions.5 Therefore, the potential importance of •Cl on the
40
oxidation of tropospheric organic pollutants has been expanded from coastal areas to
41
continental urban areas.6 Previous studies have demonstrated the high oxidation
42
ability of •Cl towards volatile organic pollutants, with rate constants that are, with
43
some exceptions, an order of magnitude larger than those of hydroxyl radicals
44
(•OH).3,6-15 Thus, although the concentration of •Cl in the marine boundary layer is ca.
45
1-10% of [•OH],16,17 the reaction with •Cl can be a significant pathway governing the
46
fate of tropospheric organic pollutants.
47
Similar to the reaction with •OH, the reaction of organic pollutants with •Cl may
48
proceed via either addition of •Cl to unsaturated bonds or H-abstraction.9 However,
49
the favorable reaction pathway for the reaction of one specific compound with •Cl is
50
not always the same as that with •OH. For example, •OH prefers to add to the
51
unsaturated bond of polycyclic aromatic hydrocarbons,18 while •Cl prefers to abstract
52
their hydrogen atoms.19,20 Therefore, the transformation initiated by •Cl may lead to a
53
fate of tropospheric organic pollutants that differs from that by •OH.
54
Monoethanolamine (MEA) is a benchmark and widely utilized solvent in
55
amine-based post-combustion CO2 capture (PCCC), which is a promising technology
56
for reducing CO2 emission from fossil fuel power plants.21-28 Given the possible
57
large-scale implementation of amine-based PCCC, it is likely that there will be
58
relatively significant discharges of MEA or other alcoholamines to the atmosphere
59
from PCCC units because of their relatively high vapor pressure.29 It has been
60
estimated that a CO2 capture plant using MEA as solvent, which removes 1 million
61
tons CO2 per year from flue gas, could potentially emit 80 tons MEA into the
62
atmosphere.30-32
Therefore,
concern
about
the
atmospheric
3
ACS Paragon Plus Environment
chemistry
and
Environmental Science & Technology
63
environmental risk of the use of amines associated with the PCCC has been
64
increasing.28,33,34
65
Previous studies have investigated the reaction kinetics and pathways of MEA with
66
•OH.35-38 These studies revealed that through •OH abstracting H on the α and the β
67
site of -NH2 in MEA, the MEA+•OH reaction favorably produces MEA-α and
68
MEA-β radicals. This finally leads to the formation of NHCHCH2OH and
69
NH2CH2CHO by their further reaction with O2.35-38 NH2CH2CHO may photolyze just
70
like other aliphatic aldehydes to form NH2CH2•, which can further react with O2/NO
71
to produce formamide or imine.37,39-42 Formamide could react with ·OH to produce
72
isocyanic acid HNCO,43,44 a potential hazardous compound.45 There are no
73
experimental data available for the gas phase reactions of NHCHCH2OH or other
74
imines.33 However, imines are water-soluble and may be taken up by aqueous aerosols
75
in the troposphere. In aqueous solution, NHCHCH2OH could hydrolyze to produce
76
ammonia and glycolaldehyde.33,46
77
Similar to the MEA+•OH reaction, the MEA+•Cl reaction may also follow the
78
H-abstraction mechanism because MEA has no unsaturated bonds. However, it is
79
unknown which H-atom of MEA can be abstracted favorably. We noted that there
80
could be a particular interaction (2-center-3-electron bond) between the nitrogen lone
81
pair electrons of the amine group and the single electron occupied p-orbital of •Cl,
82
which could make the H-abstraction at the N site become the most favorable,17
83
different from the reaction MEA+•OH. If H-abstraction at the N site of MEA, leading
84
to the formation of N-center radicals, is the most favorable pathway, then the •Cl
85
initiated MEA reaction may finally form carcinogenic nitrosamine via the reaction of
86
N-center radicals with NO.47-49 In addition, the N-center radicals may react with NO2
87
to form nitramine,33,50-53 a potential mutagenic compound.54 Therefore, the formation 4
ACS Paragon Plus Environment
Page 4 of 28
Page 5 of 28
Environmental Science & Technology
88
of N-center radicals will increase the environmental risk of MEA emission from
89
PCCC units. However, to the best of our knowledge, there are no previous studies
90
concerning the atmospheric chemical reaction pathways and kinetics for the
91
transformation of MEA initiated by •Cl.
92
Therefore, it became our purpose to investigate the reaction pathways and product
93
branching ratio of the critical steps in the •Cl-initiated oxidation of MEA using a
94
combination of highly cost-expensive coupled-cluster theory (CCSD(T))55 and
95
kinetics modeling. We further investigated the subsequent reactions of the
96
MEA-radicals formed in the initial steps with O2 and NOx (NO and NO2) as they are
97
important oxidative agents in radical-initiated atmospheric chemical reactions.55
98
COMPUTATIONAL DETAILS
99
Ab Initio Electronic Structure and Kinetics Calculations. All the structures and
100
energy calculations were performed using the GAUSSIAN 09 program package.56 The
101
geometry optimizations and harmonic frequency calculations for the reactants,
102
products,
103
MP2/6-31+G(3df,2p) level.55 The connections of the transition states between
104
designated local minima were confirmed by intrinsic reaction coordinate (IRC)
105
calculation at the MP2/6-31+G(3df,2p) level. Single-point energy calculation was
106
performed at the CCSD(T)/aug-cc-pVTZ level55 based on the geometries at the
107
MP2/6-31+G(3df,2p) level. In the application to atmospheric chemistry problems, the
108
CCSD(T) method is likely the most popular ab initio method in use today.55 It is
109
highly accurate and expensive for energy calculation, and errors of relative energies
110
and bond energies can often be calculated to within 1 kcal mol-1.55 A similar scheme
111
involving MP2 geometry optimization and CCSD(T) single point energy calculation
112
has been successfully used to predict the reaction kinetics of simple mono-, di- and
intermediates
and
transition
states
were
5
ACS Paragon Plus Environment
performed
at
the
Environmental Science & Technology
Page 6 of 28
113
trimethylamine with •Cl.17 Zero point energy at the MP2/6-31+G(3df,2p) level was
114
used to correct the single-point energy. The most stable conformer (Figure 1) of
115
MEA57 was selected as the starting reactant.
116
We employed the MultiWell-2014.1 master equation code58-61 to calculate the
117
reaction rate constants. The master equation method is a powerful tool in calculating
118
the time-dependent, temperature-dependent and pressure-dependent kinetics of a
119
multi-channel and multi-well chemical reaction system. Reaction rate constants for
120
tight transition states were calculated from the RRKM theory,62 based on sums and
121
densities
122
CCSD(T)/aug-cc-pVTZ barrier heights. The energy-grained master equation was
123
solved over 2000 grains of 10 cm−1 each, carried on to 85000 cm−1 for the continuum
124
component of the master equation. N2 was employed as buffer gas. The collision
125
transfer probability between reactive intermediates and N2 was described by the single
126
exponential-down model63 with the average transfer energy ∆Ed = 200 cm-1. The
127
Lennard-Jones parameters for intermediates were calculated from an empirical
128
method proposed by Gilbert et al.64 The overall rate constants k(T) were calculated
129
using
130
of
states
for
the
MP2/6-31+G(3df,2p)
k(T) = k∞(T)(1 - Гreactants)
structures
and
the
(1)
131
where k∞(T) is the high-pressure-limit rate constant for the combination reaction and
132
Гreactants is the fractional yield of reactants.65
133
For both reaction MEA+•Cl and reaction N-center radicals (MEA-N) + NO, the
134
initial combination was found to be barrierless. The calculation for the reaction rate
135
constant of barrierless reactions is always challenging. Several methods can be
136
employed to calculate rate constants for barrierless reactions with respective strengths
137
and weaknesses.66 When the rate constant for a target system or its similar system is 6
ACS Paragon Plus Environment
Page 7 of 28
Environmental Science & Technology
138
known, it is convenient and reliable to use the restricted Gorin model to simulate the
139
kinetics.66 For the reaction of MEA-N radicals with NO, the capture rate for the
140
related •NH2 + •NO association (7×10−11 cm3 molecule−1 s−1) is available.49 The
141
barrierless combination of MEA-N radicals with NO was therefore treated with the
142
restricted Gorin model.67 However, the restricted Gorin model cannot be used for the
143
reaction of MEA with •Cl since no related experimental data are available. We noted
144
that the long-range transition state theory with a dispersion force potential68 had been
145
successfully used in calculating the formation rate constant of pre-complexes for the
146
reactions of mono-, di- and trimethylamine with •Cl.17 Similarly, we also used the
147
long-range transition state theory to calculate the formation rate constant of
148
pre-reactive complexes for the reaction MEA+•Cl. Additional details for the reaction
149
rate constant calculations are presented in the Supporting Information (SI).
150 151
Figure 1. MP2/6-31+G(3df,2p)-optimized geometries for MEA and some important
152
complexes and transition states involved in the reaction of MEA+•Cl.
7
ACS Paragon Plus Environment
Environmental Science & Technology
153
RESULTS AND DISCUSSION
154
Initial Reactions with •Cl. In principle, MEA should have seven H-abstraction
155
pathways, as it is constituted by seven H-atoms with different chemical environments
156
(Figure 2). The thermodynamic calculation results (Table S2) indicated that the
157
H-abstraction pathway from the -OH group is endothermic, while the other six
158
pathways are exothermic. Therefore, we excluded the contribution of the
159
H-abstraction pathway from the -OH group to the final products. A schematic
160
potential energy surface for the six remaining H-abstraction pathways is presented in
161
Figure 3, and the optimized geometries for the important species including MEA,
162
complexes and transition states are presented in Figure 1.
163 164
Figure 2. Possible pathways for the reaction of MEA with •Cl. The symbols “TS1-m
165
and P1-m” denote the transition states and products involved in the reaction and m
166
denotes the different species.
167
As shown in Figure 3, every reaction pathway proceeds through a pre-reactive and 8
ACS Paragon Plus Environment
Page 8 of 28
Page 9 of 28
Environmental Science & Technology
168
post-reactive complex. It can be concluded from the activation energies (Ea) that
169
H-abstractions occurring at the N site of MEA via the transition state TS1-1’ and TS1-1
170
are much more favorable than those occurring at the α site and the β site. Therefore,
171
the formation of MEA-N radicals is the most favorable in terms of kinetics. In
172
addition, the pathway via TS1-1’ should be slightly more favorable than that via TS1-1
173
due to its slightly lower reaction energy barrier (-8.3 kcal mol-1) than that (-7.2 kcal
174
mol-1) of the pathway via TS1-1. There is no regularity in reaction barriers occurring at
175
the α and the β site because both the local chemical environment and the
176
σ-electron-withdrawing group (-NH2 and -OH adjacent to the C-H bond) have a great
177
influence on the reaction barrier for the H-abstraction occurring at the α and the β
178
site.69 It deserves mentioning that the energy of the complex (-14.0 kcal mol-1) and Ea
179
(-8.3 kcal mol-1) in the most favorable pathway for the reaction MEA+•Cl are both
180
lower than those of the corresponding process for the reaction of the simplest primary
181
amine CH3NH2 with •Cl.17 This may result from the stabilizing role of the interaction
182
of -OH…•Cl in the pre-reactive complex and the transition state in the reaction
183
MEA+•Cl. As no well-defined transition state was located for the H-abstraction
184
occurring at the β-site (3’), an approximated transition state (TS1-3’) was shown
185
instead in Figure 3. A detailed discussion about the approximation of TS1-3’ is
186
presented in the SI.
187
We observed four pre-reactive complexes involved in the reaction MEA+•Cl
188
(Figure 1). The two complexes RC1-1 and RC1-1’ associated with the interaction
189
between the N-atom and •Cl are the most stable. This should result from the formation
190
of a 2-center-3-electron bond between the nitrogen lone pair electrons and the single
191
electron occupied p-orbital of •Cl. Similar complexes were also observed in the
192
reactions of methylamine, dimethylamine and trimethylamine with •Cl.17 The 9
ACS Paragon Plus Environment
Environmental Science & Technology
193
formation of such stable pre-reactive complexes (RC1-1 and RC1-1’) is a reason why
194
the Ea values for the H-abstractions occurring at the N-site are much lower than those
195
at the other sites. A detailed discussion on how such complexes lead to the
196
H-abstraction occurring at the N site is presented in the SI. RC1-1 is 0.3 kcal mol-1
197
lower than RC1-1’ in energy. Structurally, the N-atom in RC1-1 rotates a little bit
198
compared with MEA and RC1-1’(Figure 1). Thus, when MEA rigidly collides with •Cl,
199
RC1-1’ forms first.
200 201
Figure 3. Schematic potential energy surface for the reaction MEA +•Cl calculated at
202
the CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p) level. The total energy of the
203
reactant MEA+•Cl is set as zero (reference state). The symbols “R1, RC1-m, PC1-m,
204
TS1-m and P1-m” stand for reactants, pre-reactive complexes, post-reactive complexes,
205
transition states and products involved in the reaction, respectively; m denotes
206
different species. ∆E was calculated at 0 K. The green, blue and red pathways 10
ACS Paragon Plus Environment
Page 10 of 28
Page 11 of 28
Environmental Science & Technology
207
correspond to H-abstraction occurring at the N, α and β sites of MEA, respectively.
208
With the master equation method, the overall rate constant (kCl) was calculated to
209
be 3.6×10-10 cm3 molecule-1 s-1 for the reaction MEA+•Cl, which is close to the values
210
for the reactions of methylamine (2.9×10-10), dimethylamine (3.9×10-10) and
211
trimethylamine (3.7×10-10 cm3 molecule-1 s-1) with •Cl.17 It deserves mentioning that
212
although the canonical variational transition state theory within the POLYRATE
213
2010-A program70 usually works well for the barrierless reaction with a pre-reactive
214
complex, a reasonable value could not be calculated for the reaction rate constant of
215
MEA with •Cl, as detailed in the SI. Over the temperature range 200-376 K, the
216
calculated kCl values decrease with increasing temperature (Figure 4). As the
217
temperature range is narrow, there seems an apparent linear relationship between kCl
218
and temperature. For a wide temperature range, the relationship become nonlinear as
219
shown in Figure S3.
A
B
220
The calculated branching ratio Г values for the formation of MEA+•Cl,
221
RC1-1’(RC1-1), MEA-N+HCl, MEA-β+HCl and MEA-α+HCl are 23.37%, 0, 72.38%,
222
0.87% and 3.38% at 298 K and 1 atm, respectively, suggesting a strong preference for
223
the N site H-abstraction. The Г values of the parent MEA+•Cl increase with
224
temperature and the Г values of the main product MEA-N+HCl decrease with
225
temperature. However, the Г values for the other species are still small even though
226
they increase with temperature. This variation of the Г values with temperature could
227
result from the fact that the difference in the transformation rate constants of CP1-1’
228
towards PC1-1(finally to MEA-N+HCl) and other species decreases with increasing
229
temperature, e.g. the difference (k1/k2=2.0×10-2) in rate constants between k1 (for
230
transformation from CP1-1’ to MEA+•Cl) and k2 (from CP1-1’ to PC1-1) at 376 K is
231
smaller than that (k1/k2=2.0×10-6) at 200 K. In addition, the branching ratio of 11
ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 28
RC1-1’(RC1-1) is zero although it seems that RC1-1’(RC1-1) lies in a relatively deep well.
233
This can be due to that the excess internal energy (~14 kcal mol-1) along the formation
234
of RC1-1’(RC1-1) is high enough to make RC1-1’(RC1-1) overcome the 6-7 kcal mol-1
235
energy barrier to form other species, but cannot make it deactivated. The similar case
236
was also found for the reactions of CH3NH2, (CH3)2NH, and (CH3)3N with •Cl.17 10 3 -1 -1 kCl×10 (cm molecule s )
232
A
4.0
1.0
MEA-α+HCl MEA+Cl
0.8
3.8
RC1-1' MEA-N+HCl MEA-β+HCl
B
Г 0.6
3.6
0.4
3.4
0.2
3.2
0.0
200
250
300 T(K)
350
200 225 250 275 300 325 350 375 400 T (K)
400
237
Figure 4. Calculated reaction kinetics for the reaction MEA+•Cl. A: Variation of rate
238
constants (kCl) at 1 atm; B: Variation of branching ratios (Г) for the formation of the
239
parent MEA+•Cl, RC1-1’, MEA-N+HCl, MEA-β+HCl and MEA-α+HCl at 1 atm.
240
These results confirm our hypothesis that the reaction of •Cl with MEA leads to the
241
formation of N-center radicals (MEA-N) as the primary activated products, in contrast
242
to the reaction of •OH with MEA that forms MEA-β and MEA-α radicals.37 To the
243
best of our knowledge, this is the first report that the reaction of •Cl with MEA, a
244
benchmark and widely utilized solvent in amine-based PCCC, can produce MEA-N
245
radicals. We further investigated the atmospheric transformation of MEA-N radicals
246
as the atmospheric transformation of the minor products MEA-β and MEA-α radicals
247
has previously been investigated.37,46
248
Subsequent Reactions of MEA-N radicals. The chemically activated MEA-N
249
radicals can subsequently self-isomerize, self-dissociate or react with the main
250
atmospheric oxidants (O2 and NOx). Schematic potential energy surfaces of
251
isomerization and dissociation of the MEA-N radicals are depicted in Figure 5A, and 12
ACS Paragon Plus Environment
Page 13 of 28
Environmental Science & Technology
252
the corresponding reaction enthalpy (∆H) values are listed in Table S2. The Ea values
253
for all the pathways are > 24 kcal mol-1, indicating that the isomerization and
254
dissociation reactions proceed slowly, and consequently the MEA-N radicals have a
255
great chance to react with atmospheric O2 and NOx.
256
The calculated schematic potential energy surface for the MEA-N+O2 reaction is
257
shown in Figure 5B with the corresponding ∆H values listed in Table S2. It can be
258
observed from Figure 5B that there are two kinds of pathways for the MEA-N+O2
259
reaction: one is the H-abstraction pathway in which ground state O2 (at triplet 3O2)
260
abstracts H-atom at the α site of –NH2 to form imine and O2H (P3-1); the other is the
261
addition pathway in which O2 adds to the N-site of MEA-N radicals. The results show
262
that the overall reaction energy barrier (6.4 kcal mol-1) of H-abstraction pathway is
263
much lower than those of the addition pathways (at least 18.4 kcal mol-1). Therefore,
264
the H-abstraction pathway leading to the formation of imine and O2H is more
265
favorable. To the best of our knowledge, this is the first time to reveal that the
266
reaction of N-centered radicals with O2 follows the H-abstraction reaction mechanism.
267
With the H-abstraction pathway, the reaction rate constant kO2 was calculated to be
268
2.9×10-18 cm3 molecule-1 s-1 and the lifetime of MEA-N is 0.07 s. The slow reaction
269
MEA-N+O2 could make the MEA-N radicals have a chance to react with NOx. A
B
13
ACS Paragon Plus Environment
Environmental Science & Technology
B
A
270 271
Figure 5. Schematic potential energy surfaces for the isomerization and dissociation
272
of MEA-N radicals (A) and the reaction MEA-N+O2 (B) calculated at the
273
CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p) level. The total energy of the reactants
274
MEA-N and MEA-N+O2 is set as zero (reference state). The symbols “Rn, TSn-m and
275
Pn-m” stand for reactants, transition states and products involved in the reaction,
276
respectively; m denotes different species, n = 2 is for the isomerization and
277
dissociation reaction and n = 3 is for the reaction MEA-N+O2. ∆E was calculated at 0
278
K.
279
The calculated schematic potential energy surface for the reaction MEA-N+NO is
280
shown in Figure 6, the optimized geometries for the corresponding intermediates and
281
transition states are shown in Figure S4, and the corresponding ∆H values are listed in
282
Table S2. It can be observed from Figure 6 that the MEA-N+NO reaction can
283
barrierlessly proceed to form MEA-N-NO adducts (IM4-1, nitrosamine) with six
284
conformations, depending on the attacking direction of NO. The six conformations of
285
IM4-1 can interconvert into each other (Figure S5).
286
The excess vibrational energy along the formation of nitrosamine IM4-1 could 14
ACS Paragon Plus Environment
Page 14 of 28
Page 15 of 28
Environmental Science & Technology
287
promote IM4-1 to isomerize or dissociate to other species. A recent study also indicates
288
that hot nitrosamine CH3NHNO formed from CH3NH+NO can partly isomerize via a
289
H-shift to CH3NNOH or CH2NHNOH.49 From IM4-1, eight H-shift processes (e.g. 1,3
290
H-shift from the N-atom to the O-atom of -NO, 1,3 H-shift from the Cα-atom to the
291
N-atom of -NO, and 1,4 H-shift from the Cα-atom of -NH2 to the O-atom of –NO)
292
were identified. Among these processes, the 1,3 H-shifts from the N-atom to the
293
O-atom of -NO for three IM4-1 conformations, via respectively TS4-1, TS4-2 and TS4-3,
294
to form three different conformations of IM4-2 are the most favorable. The three
295
formed conformations of IM4-2 are in the trans-form in view of the O-site of -NO and
296
the Cα-site and can interconvert into each other as shown in Figure S6. The trans-form
297
IM4-2g can convert to the cis-form IM4-2h via TS4-9. The cis-form HOCH2CH2NNOH
298
(IM4-2h) can dehydrate via TS4-10 to form diazoethanol (OHCH2CHNN) and water
299
(product set P4-1). This dehydration transition state (TS4-10) is comparable in energy to
300
TS4-9 and TS4-1. Importantly, the dissociation of the chemically activated IM4-1 to
301
OHCH2CHNN+H2O proceeds via the pathways for which the transition states are at
302
least 13 kcal mol-1 below the reactants. Thus, the reaction pathway via the
303
intermediate IM4-1, finally leading to the formation of OHCH2CHNN+H2O, is the
304
most favorable for the reaction MEA-N+NO. It deserves mentioning that the relative
305
energies of IM4-1 and all the transition states in the most favorable channel for the
306
reaction of MEA-N+NO are in general 1.0 kcal mol-1 higher than those of the
307
corresponding species for the reaction CH3NH+NO.49
15
ACS Paragon Plus Environment
Environmental Science & Technology
20.0
Page 16 of 28
∆Ε (kcal mol-1) TS 4-8 14.6
10.0
TS4-7 6.9 MEA-N+NO
0.0
IM4-3 7.8
TS4-5 -2.8 TS4-6 -1.4
R4 0.0 TS4-4 -3.1
IM4-2b -13.3 IM4-2e -12.4
-10.0
TS4-9 -13.6
TS4-3 -15.7 TS4-2 -16.2
-20.0
TS4-10 -14.2
IM4-2d -15.7 TS4-1-16.4 OHCH CHNH+HNO 2 P4-2 -19.9
-30.0 TS4-3a -39.4 -40.0 IM4-2f-47.2 IM4-2a -47.2 -50.0
-60.0
IM4-1a -44.5 IM4-1b -45.6 IM4-1c -44.5 IM4-1d -45.8 IM4-1e -45.6 IM4-1f -44.5
IM4-2c -47.3
IM4-2g-45.8
IM4-2h -47.2
P4-1 -54.9 NNCHCH2OH+H2O
308 309
Figure 6. Schematic potential energy surface for the MEA-N+NO reaction calculated
310
at the CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p) level. The total energy of the
311
reactant MEA-N+NO is set as zero (reference state). The symbols “R4, IM4-m, TS4-m
312
and P4-m” stand for reactants, intermediates, transition states and products,
313
respectively; m denotes different species. ∆E was calculated at 0 K.
314
To evaluate the branching ratio of important species involved, the formation rate
315
constant of IM4-1 and the overall rate constants (kNO) for the reaction MEA-N+NO, we
316
performed a master equation simulation on the favorable pathway, that is, the pathway
317
via the intermediate IM4-1 to finally form P4-1. The calculated kNO at 298 K and 1 atm
318
is 6.99×10-11 cm3 molecule-1 s-1 and the phenomenological reaction rate constant for
319
the formation of IM4-1 is 6.06 ×10-11 cm3 molecule-1 s-1. As can be seen from Figure 7,
320
in the temperature range of 200-376 K and 1 atm, the calculated kNO values increase
321
with temperature. Therefore, in some cases where temperature and concentration of 16
ACS Paragon Plus Environment
Page 17 of 28
Environmental Science & Technology
322
NO are decreased, e.g. from lower-troposphere to mid-troposphere,71 the reaction rate
323
of MEA-N with NO is decreased if we suppose that the concentration of MEA-N
324
remains constant. In the studied temperature range, the reactants predominantly form IM4-1f
326
(nitrosamine) with Г > 70% and the Г values for the other species (IM4-2f, IM4-2g,
327
IM4-2h and P4-1) are small (Figure 7). The product branching ratio of nitrosamine in the
328
reaction MEA-N+NO is about 10% higher than that in the reaction CH3NH+NO at
329
298 K and 1 atm.49 Therefore, under the tropospheric conditions, the reaction of •Cl
330
with MEA can finally lead to the formation of nitrosamine that has been proven to
331
have strong carcinogenic activities.72-74 The formed nitrosamine may also photolyze
332
in daytime.38,50 A recent modelling study predicted that the concentration of
333
nitrosamine is 0.3 pg m-3 in •OH initiated MEA-oxidation based on the assumption
334
that the Г value of MEA-N is 15% in the reaction MEA+•OH, nitrosamine can be
335
formed from the reaction MEA-N+NO and nitrosamine can photolyze in daytime.50
336
Thus, due to the high Г value (72%) of MEA-N in the reaction MEA+•Cl, the
337
concentration of nitrosamine could be higher than that estimated without the
338
consideration of •Cl. 10
3
7 6
IM4-1f IM4-2g
1.0 0.8
8 Г
IM4-2f IM4-2h P4-1
B
0.6 0.4 0.2
5
0.0
4 200
339
1.2
A
9
11
-1 -1
kNO×10 (cm molecule s )
325
250
300 T(K)
350
200 225 250 275 300 325 350 375 400 T (K)
400
Figure 7. Calculated reaction kinetics for the reaction MEA-N+NO. A: Variation of 17
ACS Paragon Plus Environment
Environmental Science & Technology
340
reaction rate constants (kNO) at 1 atm; B: Variation of branching ratios (Г) for the
341
formation of IM4-1f, IM4-2f, IM4-2g, IM4-2h and P4-1 at 1 atm.
342
With the calculated values of kO2 and kNO, we evaluated the reaction competition
343
between MEA-N+O2 and MEA-N+NO. The kO2/kNO ratio is 4.1×10-6, which is close
344
to the corresponding experimental value (1.5×10-6) for the (CH3)2N radicals.52 If we
345
assume that the atmospheric concentration of NO is 5 ppb, which is commonly
346
encountered in urban atmospheres and could be achieved in polluted air masses where
347
MEA emissions are expected to be significant (e.g., coal-fired power stations),46,52
348
then the lifetime of MEA-N with respect to the NO reaction is calculated to be 0.11 s.
349
The lifetime of MEA-N with respect to the NO reaction is comparable to that (0.07 s)
350
with respect to reaction with O2. Therefore, the reaction of MEA with NO is an
351
important removal pathway for MEA-N radicals at NO levels of around 5 ppb or
352
above.
353
The reaction with NO2 can also be a removal pathway for the N-center radicals to
354
form nitramine and imine.33,50-53 We further investigated the potential energy surface
355
(Figure 8) of MEA-N with NO2. The initial combination of MEA-N and NO2
356
resulting in two different adducts (nitramine IM5-1 and N-nitrosooxy amine IM5-2) is
357
barrierless, which means that the reaction rate constant of MEA-N with NO2 should
358
be high, just like the MEA-N+NO reaction. IM5-1 and IM5-2 have three and eight
359
conformations, respectively, depending on the attacking direction of NO2 on MEA-N.
360
The IM5-2 can undergo MEA-NO-NO bond rupture to produce NO and nitroxide
361
radicals without any transition state. However, IM5-1 needs to overcome a high
18
ACS Paragon Plus Environment
Page 18 of 28
Page 19 of 28
Environmental Science & Technology
362
reaction barrier to form imine and HONO(P5-1). By comparing the potential energy
363
surfaces, it was found that nitramine IM5-1 in the reaction MEA-N+NO2 is trapped in
364
a deeper well than nitroamine in the reaction MEA-N+NO. Therefore, nitramine IM5-1,
365
a potentially mutagenic compound,54 can be formed by collisional deactivation, just
366
like the formation of nitroamine in the reaction of MEA-N with NO. We can conclude
367
that the reaction of MEA-N with NO2 will mainy lead to the formation of nitramine
368
and nitroxide radicals with high reaction rate constants like in case of other N-center
369
radicals.52-53
370 371
Figure 8. Schematic potential energy surface for the MEA-N+NO2 reaction calculated
372
at the CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p) level. The total energy of the
373
reactants MEA-N+NO2 is set as zero (reference state). The symbols “R5, IM5-m, TS5-m
374
and P5-m” stand for reactants, intermediates, transition states and products involved in
375
the reaction, respectively; m denotes different species. ∆E was calculated at 0 K.
376
Implications. The calculated kCl of the MEA+•Cl reaction is higher than that (kOH) of
377
the MEA+•OH reaction36,37 by a factor of 5 (at 298K) and is close to the rate
378
constants for the reactions of methylamine (2.9×10-10), dimethylamine (3.9×10-10) and 19
ACS Paragon Plus Environment
Environmental Science & Technology
379
trimethylamine (3.7×10-10 cm3 molecule-1 s-1) with •Cl. As in the marine boundary
380
layer, •Cl concentrations [•Cl] are estimated to be as much as 1-10% of the [•OH],16,17
381
the contribution of •Cl to the transformation of MEA is about 5-50% (estimated by
382
kCl[•Cl]/kOH[•OH]) of the contribution of •OH. Based on the atmospheric
383
concentrations of •OH (9.7 × 105 molecules cm-3)75 and •Cl (9.7 × 103 - 9.7 × 104
384
molecules cm-3),16,17 the tropospheric lifetimes of MEA with respect to the reaction
385
with •OH (τOH) and a total tropospheric lifetimes with respect to both reaction with
386
•OH and •Cl (τOH,Cl) were calculated to be τOH = 3.8 hours and τOH,Cl = 2.6-3.6 hours.
387
Thus, the tropospheric lifetimes of MEA will be overestimated by 6-46 % if the role
388
of •Cl is ignored, which further proves the importance of •Cl in the transformation of
389
MEA, especially when considering the huge area of oceans (71% of the earth’s
390
surface) as well as the potential source of •Cl from continental urban areas.5
391
This study for the first time reveals that the N-center radicals MEA-N are produced
392
most favorably from the MEA+•Cl reaction. The yield of the MEA-N radicals in the
393
reaction of MEA with •Cl is about 5 times higher than that of the reaction
394
MEA+•OH.36 The subsequent atmospheric reaction of the MEA-N radicals with O2
395
and NOx finally leads to the formation of imine, nitroxide radicals, nitramine and
396
nitrosamine, and the contribution of •Cl to their formation is about 25-250%
397
(estimated by kCl[•Cl]ГCl,MEA-N/kOH[•OH]ГOH,MEA-N, where ГCl,MEA-N and ГOH,MEA-N are
398
the branching ratios of MEA-N in the reaction MEA+•Cl and MEA+•OH,
399
respectively) of that of •OH in marine environments. Thus, if •Cl is not considered
400
when assessing the risk of MEA, the risk will be underestimated significantly.
401
This study further proves that the reaction pathways and products for some organic
402
pollutants with •Cl can be different from those with •OH, implying that more studies
403
should be performed on •Cl initiated tropospheric degradation of volatile organic 20
ACS Paragon Plus Environment
Page 20 of 28
Page 21 of 28
Environmental Science & Technology
404
pollutants for their fate assessment, although it actually has already been taken into
405
consideration for decades.2
406 407
ASSOCIATED CONTENT
408
Supporting Information. Texts, figures, and tables giving detailed discussion about
409
an approximated transition state for the H-abstraction occurring at the β site, master
410
equation calculations and long-range transition state theory treatment, reaction
411
process of the H-abstraction occurring at the N site of MEA, variational transition
412
state theory calculation for the reaction MEA+•Cl, geometries involved in the reaction
413
MEA-N+NO, potential energy surfaces for the interconversion of different
414
conformers, and energetic values for the reactions. This material is available free of
415
charge via the Internet at http://pubs.acs.org.
416 417
AUTHOR INFORMATION
418
Corresponding Author
419
*
420
ACKNOWLEDGEMENTS
421
We thank Prof. Willie Peijnenburg (Leiden University) for improving the contents of
422
the manuscript and Prof. John R. Barker (University of Michigan) for providing the
423
MultiWell-2014.1 program and instructions on the calculation. The study was
424
supported by the National Natural Science Foundation of China (21207016,
425
21325729), and Program for Changjiang Scholars and Innovative Research Team in
426
University (IRT_13R05).
Phone/fax: +86-411-84706269; e-mail:
[email protected].
427
21
ACS Paragon Plus Environment
Environmental Science & Technology
428
REFERENCES
429
(1) Keene, W. C.; Khalil, M. A. K.; Erickson, D. J.; McCulloch, A.; Graedel, T. E.;
430
Lobert, J. M.; Aucott, M. L.; Gong, S. L.; Harper, D. B.; Kleiman, G.; Midgley, P.;
431
Moore, R. M.; Seuzaret, C.; Sturges, W. T.; Benkovitz, C. M.; Koropalov, V.; Barrie,
432
L. A.; Li, Y. F. Composite Global Emissions of Reactive Chlorine from
433
Anthropogenic and Natural Sources: Reactive Chlorine Emissions Inventory. J.
434
Geophys. Res. 1999, 104, 8429-8440.
435
(2) Finlayson-Pitts, B. J., Chlorine Chronicles. Nat. Chem. 2013, 5, 724-724.
436
(3) Faxon, C. B.; Allen, D. T. Chlorine Chemistry in Urban Atmospheres: A Review.
437
Environ. Chem. 2013, 10, 221-233.
438
(4) Knipping, E. M.; Lakin, M. J.; Foster, K. L.; Jungwirth, P.; Tobias, D. J.; Gerber,
439
R. B.; Dabdub, D.; Finlayson-Pitts, B. J. Experiments and Simulations of
440
Ion-Enhanced Interfacial Chemistry on Aqueous NaCl Aerosols. Science 2000, 288,
441
301-306.
442
(5) Thornton, J. A.; Kercher, J. P.; Riedel, T. P.; Wagner, N. L.; Cozic, J.; Holloway,
443
J. S.; Dubé, W. P.; Wolfe, G. M.; Quinn, P. K.; Middlebrook, A. M.; Alexander, B.;
444
Brown, S. S. A Large Atomic Chlorine Source Inferred from Mid-continental
445
Reactive Nitrogen Chemistry. Nature 2010, 464, 271-274.
446
(6) Young, C. J.; Washenfelder, R. A.; Edwards, P. M.; Parrish, D. D.; Gilman, J. B.;
447
Kuster, W. C.; Mielke, L. H.; Osthoff, H. D.; Tsai, C.; Pikelnaya, O.; Stutz, J.; Veres,
448
P. R.; Roberts, J. M.; Griffith, S.; Dusanter, S.; Stevens, P. S.; Flynn, J.; Grossberg, N.;
449
Lefer, B.; Holloway, J. S.; Peischl, J.; Ryerson, T. B.; Atlas, E. L.; Blake, D. R.;
450
Brown, S. S. Chlorine as a Primary Radical: Evaluation of Methods to Understand Its
451
Role in Initiation of Oxidative Cycles. Atmos. Chem. Phys. 2014, 14, 3427-3440.
452
(7) Ceacero-Vega, A. A.; Ballesteros, B.; Albaladejo, J.; Bejan, I.; Barnes, I.
453
Temperature Dependence of the Gas-phase Reactions of Cl Atoms with Propene and
454
1-butene between 285< T< 315K. Chem. Phys. Lett. 2009, 484, 10-13.
455
(8) Atkinson, R.; Baulch, D. L.; Cox, R. A.; Hampson, R. F.; Kerr , J. A.; Troe, J.
456
Evaluated Kinetic and Photochemical Data for Atmospheric Chemistry: Supplement
457
III. IUPAC Subcommittee on Gas Kinetic Data Evaluation for Atmospheric
458
Chemistry. J. Phys. Chem. Ref. Data. 1989, 18, 881-1097.
459
(9) Sun, C.; Xu, B.; Zhang, S. Atmospheric Reaction of Cl + Methacrolein: A
460
Theoretical Study on the Mechanism, and Pressure- and Temperature-Dependent Rate
461
Constants. J. Phys. Chem. A. 2014, 118, 3541-3551.
462
(10) Wang, W.; Ezell, M. J.; Ezell, A. A.; Soskin, G.; Finlayson-Pitts, B. J. Rate 22
ACS Paragon Plus Environment
Page 22 of 28
Page 23 of 28
Environmental Science & Technology
463
Constants for the Reactions of Chlorine Atoms with a Series of Unsaturated
464
Aldehydes and Ketones at 298 K: Structure and Reactivity. Phys. Chem. Chem. Phys.
465
2002, 4, 1824-1831.
466
(11) Canosa-Mas, C.; Hutton-Squire, H.; King, M.; Stewart, D.; Thompson, K.;
467
Wayne, R. Laboratory Kinetic Studies of the Reactions of Cl Atoms with Species of
468
Biogenic Origin: δ3-Carene, Isoprene, Methacrolein and Methyl Vinyl Ketone. J.
469
Atmos. Chem. 1999, 34, 163-170.
470
(12) Finlayson-Pitts, B. J.; Keoshian, C. J.; Buehler, B.; Ezell, A. A. Kinetics of
471
Reaction of Chlorine Atoms with Some Biogenic Organics. Int. J. Chem. Kinet. 1999,
472
31, 491-499.
473
(13) Orlando, J. J.; Tyndall, G. S.; Apel, E. C.; Riemer, D. D.; Paulson, S. E. Rate
474
Coefficients and Mechanisms of the Reaction of Cl-Atoms with a Series of
475
Unsaturated Hydrocarbons under Atmospheric Conditions. Int. J. Chem. Kinet. 2003,
476
35, 334-353.
477
(14) Li, J.; Cao, H.; Han, D.; Li, M.; Li, X.; He, M.; Ma, S. Computational Study on
478
the Mechanism and Kinetics of Cl-Initiated Oxidation of Vinyl Acetate. Atmos.
479
Environ. 2014, 94, 63-73.
480
(15) Han, D.; Cao, H.; Li, M.; Li, X.; Zhang, S.; He, M.; Hu, J. Computational Study
481
on the Mechanisms and Rate Constants of the Cl-Initiated Oxidation of Methyl Vinyl
482
Ether in the Atmosphere. J. Phys. Chem. A. 2015, 119, 719-727.
483
(16) Wingenter, O. W.; Sive, B. C.; Blake, N. J.; Blake, D. R.; Rowland, F. S. Atomic
484
Chlorine Concentrations Derived from Ethane and Hydroxyl Measurements over the
485
Equatorial Pacific Ocean: Implication for Dimethyl Sulfide and Bromine Monoxide. J.
486
Geophys. Res. 2005, 110, 1-10.
487
(17) Nicovich, J. M.; Mazumder, S.; Laine, P. L.; Wine, P. H.; Tang, Y.; Bunkan, A. J.
488
C.; Nielsen, C. J. An Experimental and Theoretical Study of the Gas Phase Kinetics of
489
Atomic Chlorine Reactions with CH3NH2, (CH3)2NH, and (CH3)3N. Phys. Chem.
490
Chem. Phys. 2015, 17, 911-917.
491
(18) Atkinson, R.; Arey, J.; Zielinska, B.; Aschmann, S. M. Kinetics and Products of
492
the Gas-Phase Reactions of OH Radicals and N2O5 with Naphthalene and Biphenyl.
493
Environ. Sci. Technol. 1987, 21, 1014-1022.
494
(19) Riva, M.; Healy, R. M.; Flaud, P.-M.; Perraudin, E.; Wenger, J. C.; Villenave, E.
495
Kinetics of the Gas-Phase Reactions of Chlorine Atoms with Naphthalene,
496
Acenaphthene, and Acenaphthylene. J. Phys. Chem. A. 2014, 118, 3535-3540.
497
(20) Wang, L.; Arey, J.; Atkinson, R. Reactions of Chlorine Atoms with a Series of 23
ACS Paragon Plus Environment
Environmental Science & Technology
498
Aromatic Hydrocarbons. Environ. Sci. Technol. 2005, 39, 5302-5310.
499
(21) Veawab, A.; Tontiwachwuthikul, P.; Chakma, A. Corrosion Behavior of Carbon
500
Steel in the CO2 Absorption Process Using Aqueous Amine Solutions. Ind. Eng.
501
Chem. Res. 1999, 38, 3917-3924.
502
(22) Liu, Y.; Zhang, L.; Watanasiri, S. Representing Vapor−Liquid Equilibrium for an
503
Aqueous MEA−CO2 System Using the Electrolyte Nonrandom-Two-Liquid Model.
504
Ind. Eng. Chem. Res. 1999, 38, 2080-2090.
505
(23) Puxty, G.; Rowland, R.; Allport, A.; Yang, Q.; Bown, M.; Burns, R.; Maeder, M.;
506
Attalla, M. Carbon Dioxide Postcombustion Capture: A Novel Screening Study of the
507
Carbon Dioxide Absorption Performance of 76 Amines. Environ. Sci. Technol. 2009,
508
43, 6427-6433.
509
(24) Xie, H. B.; Zhou, Y.; Zhang, Y.; Johnson, J. K. Reaction Mechanism of
510
Monoethanolamine with CO2 in Aqueous Solution from Molecular Modeling. J. Phys.
511
Chem. A. 2010, 114, 11844-11852.
512
(25) Xie, H. B.; Johnson, J. K.; Perry, R. J.; Genovese, S.; Wood, B. R. A
513
Computational Study of the Heats of Reaction of Substituted Monoethanolamine with
514
CO2. J. Phys. Chem. A. 2011, 115, 342-350.
515
(26) Xie, H. B.; He, N.; Song, Z.; Chen, J.; Li, X. Theoretical Investigation on the
516
Different Reaction Mechanisms of Aqueous 2-Amino-2-methyl-1-propanol and
517
Monoethanolamine with CO2. Ind. Eng. Chem. Res. 2014, 53, 3363-3372.
518
(27) Xie, H. b.; Wei, X.; Wang, P.; He, N.; Chen, J. CO2 Absorption in an Alcoholic
519
Solution of Heavily Hindered Alkanolamine: The Reaction Mechanism of
520
2-(Tert-Butylamino)- Ethanol with CO2 Revisited. J. Phys. Chem. A. 2015, 119,
521
6346–6353.
522
(28) da Silva, E. F.; Booth, A. M. Emissions from Postcombustion CO2 Capture
523
Plants. Environ. Sci. Technol. 2013, 47, 659-660.
524
(29) Kapteina, S.; Slowik, K.; Verevkin, S. P.; Heintz, A. Vapor Pressures and
525
Vaporization Enthalpies of a Series of Ethanolamines. J. Chem. Eng. Data, 2005, 50,
526
398-402.
527
(30) Veltman, K.; Singh, B.; Hertwich, E. G. Human and Environmental Impact
528
Assessment of Postcombustion CO2 Capture Focusing on Emissions from
529
Amine-Based Scrubbing Solvents to Air. Environ. Sci. Technol. 2010, 44, 1496-1502.
530
(31) Karl, M.; Wright, R. F.; Berglen, T. F.; Denby, B. Worst Case Scenario Study to
531
Assess the Environmental Impact of Amine Emissions from a CO2 Capture Plant. Int.
532
J. Greenhouse Gas Control, 2011, 5, 439-447. 24
ACS Paragon Plus Environment
Page 24 of 28
Page 25 of 28
Environmental Science & Technology
533
(32) Shao, R.; Stangeland, A. Amines Used in CO2 Capture-Health and
534
Environmental Impacts. Report: The Bellona Foundation: Oslo, Norway, 2009.
535
(33) Nielsen, C. J.; Herrmann, H.; Weller, C. Atmospheric Chemistry and
536
Environmental Impact of the Use of Amines in Carbon Capture and Storage (CCS).
537
Chem. Soc. Rev. 2012, 41, 6684-6704.
538
(34) Dai, N.; Shah, A. D.; Hu, L.; Plewa, M. J.; McKague, B.; Mitch, W. A.
539
Measurement of Nitrosamine and Nitramine Formation from NOx Reactions with
540
Amines during Amine-Based Carbon Dioxide Capture for Postcombustion Carbon
541
Sequestration. Environ. Sci. Technol. 2012, 46, 9793-9801.
542
(35) Borduas, N.; Abbatt, J. P. D.; Murphy, J. G. Gas Phase Oxidation of
543
Monoethanolamine (MEA) with OH Radical and Ozone: Kinetics, Products, and
544
Particles. Environ. Sci. Technol. 2013, 47, 6377-6383.
545
(36) Onel, L.; Blitz, M. A.; Seakins, P. W. Direct Determination of the Rate
546
Coefficient for the Reaction of OH Radicals with Monoethanol Amine (MEA) from
547
296 to 510 K. J. Phys. Chem. Lett. 2012, 3, 853-856.
548
(37) Xie, H. B.; Li, C.; He, N.; Wang, C.; Zhang, S.; Chen, J. Atmospheric Chemical
549
Reactions of Monoethanolamine Initiated by OH Radical: Mechanistic and Kinetic
550
Study. Environ. Sci. Technol. 2014, 48, 1700-1706.
551
(38) Karl, M.; Dye, C.; Schmidbauer, N.; Wisthaler, A.; Mikoviny, T.; D'Anna, B.;
552
Müller, M.; Borrás, E.; Clemente, E.; Muñoz, A.; Porras, R.; Ródenas, M.; Vázquez,
553
M.; Brauers, T. Study of OH-Initiated Degradation of 2-Aminoethanol. Atmos. Chem.
554
Phys. 2012, 12, 1881-1901.
555
(39) Carlier, P.; Hannachi, H.; Mouvier, G., The chemistry of carbonyl compounds in
556
the atmosphere—A review. Atmos. Environ., 1986, 20, (11), 2079-2099.
557
(40) Claus Jørgen Nielsen, Barbara D’Anna, Christian Dye, Christian George, Martin
558
Graus, Armin Hansel, Matthias Karl, Stephanie King, Mihayo Musabila, Marcus
559
Müller, Norbert Schmidbauer, Yngve Stenstrøm, Armin Wisthaler, Atmospheric
560
Degradation of Amines (ADA), Summary Report: Gas phase photo-oxidation of
561
2-aminoethanol (MEA)
562
(41) Nielsen, C. J.; D’Anna, B.; Dye, C.; Graus, M.; Karl, M.; King, S.; Maguto, M.
563
M.; Müller, M.; Schmidbauer, N.; Stenstrøm, Y.; Wisthaler, A.; Pedersen, S.,
564
Atmospheric chemistry of 2-aminoethanol (MEA). Energy Procedia 2011, 4,
565
2245-2252.
566
(42) da Silva, G.; Kirk, B. B.; Lloyd, C.; Trevitt, A. J.; Blanksby, S. J., Concerted
567
HO2 Elimination from α-Aminoalkylperoxyl Free Radicals: Experimental and 25
ACS Paragon Plus Environment
Environmental Science & Technology
568
Theoretical Evidence from the Gas-Phase NH2•CHCO2– + O2 Reaction. J. Phys.
569
Chem. Lett. 2012, 3, 805-811.
570
(43) Barnes, I.; Solignac, G.; Mellouki, A.; Becker, K. H. Aspects of the Atmospheric
571
Chemistry of Amides. Chem. Phys. Chem. 2010, 11,3844−3857.
572
(44) Borduas, N.; da Silva, G.; Murphy, J. G.; Abbatt, J. P. D., Experimental and
573
Theoretical Understanding of the Gas Phase Oxidation of Atmospheric Amides with
574
OH Radicals: Kinetics, Products, and Mechanisms. J. Phys. Chem. A 2015, 119, (19),
575
4298-4308.
576
(45) Roberts, J. M.; Veres, P. R.; Cochran, A. K.; Warneke, C.; Burling, I. R.;
577
Yokelson, R. J.; Lerner, B.; Gilman, J. B.; Kuster, W. C.; Fall, R.; de Gouw, J.
578
Isocyanic Acid in the Atmosphere and Its Possible Link to Smoke-Related Health
579
Effects. Proc. Natl. Acad. Sci. 2011, 108,10236-10242.
580
(46) da Silva, G. Atmospheric Chemistry of 2-Aminoethanol (MEA): Reaction of the
581
NH2•CHCH2OH Radical with O2. J. Phys. Chem. A. 2012, 116, 10980-10986.
582
(47) Tang, Y.; Hanrath, M.; Nielsen, C. J., Do primary nitrosamines form and exist in
583
the gas phase? A computational study of CH3NHNO and (CH3)2NNO. Phys.Chem.
584
Chem. Phys. 2012, 14, (47), 16365-16370.
585
(48) Tang, Y.; Nielsen, C. J., Theoretical Study on the Formation and Photolysis of
586
Nitrosamines (CH3CH2NHNO and (CH3CH2)2NNO) under Atmospheric Conditions.
587
J. Phys. Chem. A 2013, 117, (1), 126-132.
588
(49) da Silva, G. Formation of Nitrosamines and Alkyldiazohydroxides in the Gas
589
Phase: The CH3NH + NO Reaction Revisited. Environ. Sci. Technol. 2013, 47,
590
7766-7772.
591
(50) Karl, M.; Svendby, T.; Walker, S. E.; Velken, A. S.; Castell, N.; Solberg, S.,
592
Modelling atmospheric oxidation of 2-aminoethanol (MEA) emitted from
593
post-combustion capture using WRF–Chem. Sci. Total Environ. 2015, 527-528,
594
185-202.
595
(51) Pitts, J. N.; Grosjean, D.; Vancauwenberghe, K.; Schmid, J. P.; Fitz, D. R.
596
Photooxidation of Aliphatic Amines Under Simulated Atmospheric Conditions:
597
Formation of Nitrosamines, Nitramines, Amides, and Photochemical Oxidant.
598
Environ. Sci. Technol. 1978, 12, 946−953.
599
(52) Lindley, C. R. C.; Calvert, J. G.; Shaw, J. H. Rate Studies of the Reactions of the
600
(CH3)2N radical with O2, NO, and NO2. Chem. Phys. Lett. 1979, 67, 57-62.
601
(53) Tang, Y.; Nielsen, C. J., A systematic theoretical study of imines formation from
602
the atmospheric reactions of RnNH2−n with O2 and NO2 (R = CH3 and CH3CH2; n = 1 26
ACS Paragon Plus Environment
Page 26 of 28
Page 27 of 28
Environmental Science & Technology
603
and 2). Atmos. Environ. 2012, 55, 185-189.
604
(54) Fjellsbø, L. M.; Verstraelen, S.; Kazimirova, A.; Van Rompay, A. R.;
605
Magdolenova, Z.; Dusinska, M., Genotoxic and mutagenic potential of nitramines.
606
Environ. Res. 2014, 134, 39-45.
607
(55) Vereecken, L.; Francisco, J. S. Theoretical studies of atmospheric reaction
608
mechanisms in the troposphere, Chem. Soc. Rev., 2012, 41, 6259–6293
609
(56) Frisch, M. J.; Trucks, G. W.; H.B, S.; Scuseria, G. E.; Robb, M. A.; Cheeseman,
610
J. R., et.al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, 2009.
611
(57) Vorobyov, I.; Yappert, M. C.; DuPré, D. B. Hydrogen Bonding in Monomers and
612
Dimers of 2-Aminoethanol. J. Phys. Chem. A. 2002, 106, 668-679.
613
(58) MultiWell-2014.1 Software. Designed and Maintained by Barker, J. R. with
614
Contributors Ortiz, N. F.; Preses, J. M.; Lohr, L. L.; Maranzana, A.; Stimac, P. J.;
615
Nguyen, T. L.; Kumar, T. J. D. Universityof Michigan: Ann Arbor, MI, 2014.
616
http://aoss.engin.umich.edu/multiwell/.
617
(59) Barker, J. R. Multiple-Well, Multiple-Path Unimolecular Reaction Systems. I.
618
MultiWell Computer Program Suite. Int. J. Chem. Kinet. 2001, 33, 232-245.
619
(60) Barker, J. R.; Ortiz, N. F. Multiple-Well, Multiple-Path Unimolecular Reaction
620
Systems. II. 2-Methylhexyl Free Radicals. Int. J. Chem. Kinet. 2001, 33, 246-261.
621
(61) Barker, J. R. Energy Transfer in Master Equation Simulations: A New Approach.
622
Int. J. Chem. Kinet. 2009, 41, 748-763.
623
(62) RRKM Robinson, P. J.; Holbrook, K. A.; Unimolecular reactions. John Wiley &
624
Sons; 1972
625
(63) Barker, J. R.; Yoder, L. M.; and King, K. D. Vibrational Energy Transfer
626
Modeling of Nonequilibrium Polyatomic Reaction Systems, J. Phys. Chem. A 2001,
627
105, 796-809
628
(64) Gilbert, R. G.; Smith, S. C. Theory of Unimolecular and Recombination
629
Reactions, Blackwell Scientific: Carlton, Australia, 1990.
630
(65) Akbar Ali, M.; Barker, J. R., Comparison of Three Isoelectronic Multiple-Well
631
Reaction Systems: OH + CH2O, OH + CH2CH2, and OH + CH2NH. J. Phys. Chem. A
632
2015, 119, (28), 7578-7592.
633
(66) Barker, J. R.; Ortiz, N. F.; Preses, J. M.; Lohr, L. L.; Maranzana, A.; Stimac, P. J.;
634
Nguyen, T. L.; and Dhilip Kumar, T. J. MultiWell Program Suite User Manual,
635
(MultiWell-2014.1), http://aoss-research.engin.umich.edu/multiwell/
636
(67) Smith, G. P.; Golden, D. M. Application of RRKM Theory to the Reactions OH
637
+ NO2 + N2 → HONO2 + N2 (1) and ClO + NO2 + N2 → ClONO2 + N2 (2): A 27
ACS Paragon Plus Environment
Environmental Science & Technology
638
Modified Gorin Model Transition State. Int. J. Chem. Kinet. 1978, 10, 489-501.
639
(68) Georgievskii, Y.; Klippenstein, S. J. Long-Range Transition State Theory. J.
640
Chem. Phys. 2005, 122, 194103.
641
(69) Chan, B.; O’Reilly, R. J.; Easton, C. J.; Radom, L., Reactivities of Amino Acid
642
Derivatives Toward Hydrogen Abstraction by Cl• and OH•. J. Org. Chem. 2012, 77,
643
(21), 9807-9812.
644
(70) Zheng, J.; Zhang, S.; Lynch, B. J.; Corchado, J. C.; Chuang, Y. Y.; Fast, P. L.;
645
Hu, W. P.; Liu, Y. P.; Lynch, G. C.; Nguyen, K. A.; Jackels, C. F.; Ramos, A. F.;
646
Ellingson, B. A.; Melissas, V. S.; Villà, J.; Rossi, I.; Coitiño, E. L.; Pu, J. Z.; Albu, T.
647
V.; Steckler, R.; Garrett, B. C.; Truhlar, D. G. POLYRATE, version 2010-A;
648
Software Manager, University of Minnesota: Minneapolis, MN, 2010.
649
(71) Hudman, R. C.; Jacob, D. J.; Turquety, S. et al. Surface and lightning sources of
650
nitrogen oxides over the United States: Magnitudes, chemical evolution, and outflow.
651
Journal of Geophysical Research: Atmospheres 2007, 112, D12S05.
652
(72) Magee, P. N. Toxicity of nitrosamines: their possible human health hazards. Food
653
Cosmet. Toxicol. 1971, 9, 207-218.
654
(73) Scherf,
655
N-nitrodimethylamine and N-nitromethylamine in the rat. Carcinogenesis, 1989, 10,
656
1977-1981.
657
(74) Hassel, M.; Frei, E.; Streeter, A. J.; Wiessler, M. Pharmacokinetics of
658
N-nitrodimethylamine and N-nitromethylamine in the rat. Carcinogenesis, 1990, 11,
659
357-360.
660
(75) Prinn, R. G.; Weiss, R. F.; Miller, B. R.; Huang, J.; Alyea, F. N.; Cunnold, D. M.;
661
Fraser, P. J.; Hartley, D. E.; Simmonds, P. G. Atmospheric Trends and Lifetime of
662
CH3CCl3 and Global OH concentrations. Science, 1995, 269, 187-192.
H. R.; Frei,
E.; Wiessler, M.
Carcinogenic properties of
28
ACS Paragon Plus Environment
Page 28 of 28