of 15
Low temperature growth of sub 10 nm particles by
ammonium nitrate condensation
Neil M. Donahue,
*
a
Mao Xiao,
b
Ruby Marten,
b
Mingyi Wang,
ad
Weimeng Kong,
h
Meredith Schervish,
ac
Qing Ye,
a
Victoria Hofbauer,
a
Lubna Dada,
b
Jonathan Duplissy,
g
Henning Finkenzeller,
g
Hamish Gordon,
a
Jasper Kirkby,
fj
Houssni Lamkaddam,
b
Vladimir Makhmutov,
k
Maxim Philippov,
k
Birte Rörup,
g
Rainer Volkamer,
i
Dongyu Wang,
b
Stefan K. Weber,
j
Richard C. Flagan,
h
Dominik Stolzenburg
e
and Imad El Hadad
b
Co-condensation of nitric acid and ammonia vapors to form ammonium nitrate transforms from a fully
semi-volatile behavior when it is relatively warm (273 K and above, typical of the seasonal planetary
boundary layer) into e
ff
ectively non-volatile and irreversible uptake for the limiting vapor when it is cold
(well below 273 K, typical of the upper troposphere and occasionally the wintertime boundary layer). This
causes the system to switch in character from the one governed by semi-volatile equilibrium (how it is
usually portrayed) to the one governed by irreversible reactive uptake to even the smallest particles.
Uptake involves an activation diameter, which can be as small as 1 nm for typical vapor concentrations,
and subsequent growth rates can be very high, exceeding 1000 nm h
1
. In addition to this somewhat
surprising behavior, the system provides an exemplary case for semi-volatile reactive uptake within the
context of volatility and saturation ratios.
Environmental signi
cance
Ammonium nitrate is a major constituent of submicron aerosol particles, and with rapidly decreasing sulfur emissions, it is likely to become even mor
e
important to submicron aerosol composition in the future. Though it is o
en treated with equilibrium thermodynamics in chemical transport models, when it
is very cold, the condensation of the less abundant (rate limiting) vapor, either ammonia or nitric acid, becomes irreversible. This can drive growth
of extremely
small particles and may be a prominent growth mechanism in the parts of the upper troposphere such as the Asian tropopause aerosol layer, where nitrate
is
known to be a major aerosol constituent.
1 Introduction
In atmospheric aerosols, ammonium nitrate (NH
4
$
NO
3
)isin
many ways a quintessential example of semi-volatile reactive
uptake, with strong temperature and humidity dependent par-
titioning between the gas-phase acid
base pair and the
condensed-phase (aerosol) crystalline or dissolved salt.
1
NH
3
ð
v
Þþ
HONO
2
ð
v
Þ
%
H
2
O
NH
4
$
NO
3
ð
a
Þ
The equilibrium and volatility are a
ff
ected by both water
activity (for an aqueous solution) and of course temperature (for
any aerosol phase). Both ammonia and nitric acid are by
themselves highly volatile; nitric acid has a saturation vapor
pressure of 117 Pa at a triple point of 235 K,
2
while ammonia has
a saturation vapor pressure of 6063 Pa at a triple point of 195 K.
3
Consequently, by themselves, neither would be found in the
condensed phase under any realistic conditions of Earth's
atmosphere. It is only the proton transfer reaction to make the
ionized salt solution that renders them at least semi volatile.
The system thus serves as a canonical example of a condensed-
phase reaction (proton transfer) that drives volatility.
a
Carnegie Mellon University, Department of Chemistry, Pittsburgh, PA, USA. E-mail:
nmd@andrew.cmu.edu; Tel: +1 412 268-4415
b
Laboratory for Atmospheric Chemistry, Paul Scherrer Institute, Villigen, Switzerland
c
University of California, Irvine Department of Chemistry, Irvine, CA, USA
d
University of Chicago Department of the Geophysical Sciences, Chicago, IL, USA
e
Technical University of Vienna, Vienna, Austria
f
Institute, for Atmospheric and Environmental Sciences, Goethe University Frankfurt,
Altenhöferallee 1, 60438 Frankfurt am Main, Germany
g
Institute for Atmospheric and Earth System Research/Physics, Faculty of Science,
University of Helsinki, 00014 Helsinki, Finland
h
Division of Chemistry and Chemical Engineering, California Institute of Technology,
Pasadena, CA 91125, USA
i
Department of Chemistry, University of Colorado, Boulder, CO 80305, USA
j
CERN, Geneva, CH-1211, Switzerland
k
P.N. Lebedev Physical Institute of the Russian Academy of Sciences, 119991 Moscow,
Russia
Cite this:
Environ. Sci.: Atmos.
,2025,
5
,
67
Received 17th August 2024
Accepted 19th October 2024
DOI: 10.1039/d4ea00117f
rsc.li/esatmospheres
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environ. Sci.: Atmos.
,2025,
5
,67
81 |
67
Environmental Science:
Atmospheres
PAPER
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
View Journal
| View Issue
Because the individual constituents are so volatile and the
equilibrium partitioning between the vapor and aerosol phases
is so sensitive to temperature, most treatments of nitrate par-
titioning assume that (submicron) particles equilibrate rapidly,
4
though supermicron particles (o
en with di
ff
erent cations)
require a dynamical treatment.
5
Because of this, ammonium
nitrate has rarely been regarded as important to nucleation and
growth of nanoparticles. However, when it is cold enough and
when there is a su
ffi
cient quantity (production or emission rate)
of ammonia and nitric acid, ammonium nitrate will participate
in the growth of extremely small particles and eventually
nucleation itself.
6
9
Our objective here is to explore the transi-
tion of ammonium nitrate partitioning from an essentially
semi-volatile nature to an e
ff
ectively non-volatile (irreversible in
the limiting vapor) behavior.
2 Setup
We are interested in the transition of ammonium nitrate par-
titioning from being qualitatively
semi volatile
at relatively
high temperature to qualitatively
irreversible
at lower
temperature. In this case the individual constituents remain
volatile under all relevant atmospheric conditions, but the
equilibrium reaction to an e
ff
ectively non-volatile salt becomes
e
ff
ectively irreversible and ultimately rate-limited by uptake of
the limiting vapor. Because we are interested in the low
temperature behavior, here we shall consider only the solid,
crystalline phase of ammonium nitrate. We need to consider
three fully reversible processes:
NH
3
(v)
%
NH
3
(s)
(R1)
HONO
2
(v)
%
HONO
2
(s)
(R2)
NH
3
(s) + HONO
2
(s)
%
NH
4
$
NO
3
(s)
(R3)
A full treatment, especially with other ions and for aqueous
solutions, would break reaction R3 into its elementary steps as
well.
HONO
2
(s)
%
NO
3
+H
+
NH
3
(s) + H
+
%
NH
4
+
NH
4
+
+NO
3
%
NH
4
$
NO
3
(s)
However, as we are considering pure, solid ammonium
nitrate we omit these reactions here.
The experimental evidence motivating this analysis is
observations from the CERN CLOUD experiment demonstrating
that systems containing mixtures of sulfuric acid, nitric acid,
and ammonia (and in some cases oxidized organics) show clear
evidence for nucleation followed by activation and very rapid
growth (20
2000 nm h
1
) of small particles with 1.7 nm <
d
p
<
15 nm. The experiments were designed to mimic conditions of
a wintertime megacity and were thus carried out at 263 K and
278 K, with 1 ppbv <
c
NH
3
< 3 ppbv and 10 pptv <
c
HONO
2
< 1000
pptv. The observed activation diameters and growth rates are
shown in Fig. 1.
6
Both the activation diameters and the growth
rates showed an empirical relationship to the product
c
HONO
2
-
$
c
NH
3
. This is expected for activation as it is fundamentally
related to the saturation ratio of the vapors, and for this second-
order process the saturation ratio involves the product of the
vapor concentrations (activities). However, the quadratic
growth-rate dependence is more of a surprise, as one might
expect the growth rate to depend only on the concentration of
the limiting vapor to an expected stoichiometric uptake.
3 Thermodynamics
For our calculations we need to know the saturation vapor
pressure of ammonia,
p

NH
3
;
the saturation vapor pressure,
p

i
;
of nitric acid,
p

HONO
2
;
and the equilibrium constant for
ammonium nitrate formation in a condensed phase containing
pure saturated ammonium nitrate,
K
eq
R3
. From this we can
determine activities,
a
, and the reactive activity product in the
vapor (v) and suspended (s) phases. The gas-phase activity at
pressure,
p
, and temperature,
T
, is the saturation ratio with
Fig. 1
Activation diameter (top) and activated growth rate (bottom) of
particles
via
ammonium nitrate condensation observed in the CERN
CLOUD experiment. Both are proportional to the product of nitric acid
and ammonia vapors, and both shift by roughly a factor of 10 at 263 K
(green circles)
vs.
278 K (purple squares). Growth rates can exceed
1000 nm h
1
with activation diameters as small as 1.5 nm.
68
|
Environ. Sci.: Atmos.
,2025,
5
,67
81
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environmental Science: Atmospheres
Paper
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
respect to partial pressure (
p
i
), vapor concentration (
c
v
i
), or vapor
mole fraction (mixing ratio,
x
v
i
):
a
v
i
¼
p
i
p

i
¼
c
v
i
k
B
T
p

i
¼
x
v
i
p
p

i
a
s
NH
3
a
s
HONO
2
¼
1
K
eq
R3
(1)
The condensed-phase (suspended) activity is the product of
a fraction (mole,
x
i
,
p
; mass,
w
i
,
p
, and volume) and the appro-
priate activity coe
ffi
cient (mole,
h
i
,
p
; mass,
z
i
,
p
):
a
s
i
,
p
=
h
i
,
p
x
s
i
,
p
=
z
i
,
p
w
s
i
,
p
(2)
In practice, values are known for the overall equilibrium
between the condensed phase and the vapor, giving the vapor
activity product,
1
but it is informative to separate out the gas-
particle partitioning of the unreacted vapors and the particle-
phase proton-transfer reaction as we have done in order to
illustrate the overall process. All of the processes are
func-
tionally
reversible, making this a true semi-volatile system. We
need to know the physical and thermodynamic properties of all
constituents to describe the dynamics. Ultimately the accuracy
for this system depends on the overall (vapor to salt) equilib-
rium and so we ensure that the individual vapor pressures and
condensed-phase equilibrium constant are consistent with the
more accurate overall equilibrium constant.
1
3.1 Ammonia
Ammonia is extremely volatile. The molar mass is
M
NH
3
=
17.03026
×
10
3
kg mol
1
Its triple point properties are
T
t
;
NH
3
¼
195
:
48 K
p
t
;
NH
3
¼
6063 Pa
r
l
t
;
NH
3
¼
733
:
74 kg m

3
3.1.1 NH
3
functions.
The density and saturation vapor
pressure of ammonia were reported by Haar and Gallagher at
NIST (NBS) in the 1970s.
3
r
l
NH
3
¼½
636

1
:
43
ð
T
=

C
Þ
kg m

3
ln

p

NH
3
.
atm

¼
A
=
T
þ
B
þ
CT
þ
DT
2
þ
ET
3
A
¼
3684
:
7798 K
B
¼
20
:
428787
C
¼
0
:
02893289
D
¼
3
:
4798128

10

5
E
¼
9
:
2219845

10

9
or since 1 atm
=
1.01325
×
10
5
Pa,
B
/Pa
=
31.9549. This all
gives the vapor pressure plot under atmospheric conditions
shown in Fig. 2. That is such an enormous value it is worth
considering the gas-phase activity of a reference level of
ammonia. For this we shall consider a mole fraction of 1 ppbv,
meaning a partial pressure of
p
NH
3
=
10
4
Pa. Fig. 3 shows the
gas-phase activity
ð
a
v
NH
3
¼
p

NH
3
=
p
NH
3
Þ
of 1 ppbv of ammonia.
Even at 200 K it is 10
8
and at 280 K it is just above 10
10
. The
gas-phase activity of 1 ppbv of ammonia at +5 °C is 2.0
×
10
10
while at
10 °C it is 3.5
×
10
10
.
3.2 Nitric acid
Nitric acid is also quite volatile, with a molar mass
M
HONO
2
=
63.012
×
10
3
kg mol
1
Combining values from Mozurkewich
1
and the saturation
vapor pressure from Duisman.
2
T
t
;
HONO
2
¼
235 K
p
t
;
HONO
2
¼
117 Pa
r
l
t
;
HONO
2
¼
1513 kg m

3
Fig. 2
Saturation vapor pressure of ammonia over a pure liquid.
Fig. 3
Saturation ratio (gas-phase activity) of 1 ppbv of ammonia over
a pure liquid.
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environ. Sci.: Atmos.
,2025,
5
,67
81 |
69
Paper
Environmental Science: Atmospheres
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
3.2.1 HONO
2
functions.
The density and saturation vapor
pressure of nitric acid were reported by Duisman at Union
Carbide.
2
r
l
HONO
2
¼
1510 kg m

3
ln

p

HONO
2
.
torr

¼
A
=
ð
T
þ
230

C
Þþ
B
A
¼
1486
:
238

C
B
¼
7
:
61628
ln

p

HONO
2
.
torr

¼
A
=
ð
T
þ
230

273
:
15
Þþ
B
¼
A
=
ð
T

43
:
15 K
Þþ
B
ln

p

HONO
2
.
Pa

¼
A
=
ð
T

43
:
15 K
Þþ
B
þ
4
:
8928
B
0
¼
12
:
5091
Nitric acid has a lower vapor pressure than ammonia but not
by all that much, as we see in Fig. 4. Thus, again for a reference
mixing ratio of 1 ppbv, the gas-phase activities in Fig. 5 are
about 2
×
10
7
at 280 K and 4
×
10
6
at 200 K (over a super-
cooled liquid solution, not nitric acid trihydrate or anything like
that). The gas-phase activity of 1 ppbv of ammonia at +5 °C is 2.1
×
10
7
while at
10 °C it is 3.2
×
10
7
.
3.3 Ammonium nitrate
Finally, ammonium nitrate has the following properties:
M
NH
4
$
NO
3
¼
80
:
044

10

3
kg mol

1
r
NH
4
$
NO
3
¼
1725 kg m

3
Given the density of ammonium nitrate we can easily
calculate the number of NH
4
$
NO
3
per particle
vs.
diameter, as
shown in Fig. 6.
3.3.1 Equilibrium.
Thermodynamic data are known for
ammonium nitrate equilibrium directly between the vapor and
condensed phase
1
NH
3
(v) + HONO
2
(v)
%
NH
4
$
NO
3
(s)
However, we can consider the equilibrium of reaction R3,
where in terms of activity the equilibrium can be written as
K
eq
R3
¼
a
s
NH
4
$
NO
3
a
s
NH
3
$
a
s
HONO
2
x
1
a
s
NH
3
$
a
s
HONO
2
:
Here, we make use of the negligibly small activities of the
undissociated molecules in the suspended particle phase at
equilibrium with the vapor partial pressures we have already
calculated. For our purposes, over a
at surface the gas-particle
system is at equilibrium for 100 pptv of ammonia and about 2
ppbv of nitric acid at +5 °C and about 20 pptv of nitric acid at
10 °C.
The overall equilibrium between vapor and a saturated
condensed phase of ammonium nitrate can be expressed
via
a dissociation constant for
NH
4
$
NO
3
%
NH
3
(v) + HONO
2
(v)
Fig. 4
Saturation vapor pressure of nitric acid over a pure liquid.
Fig. 5
Saturation ratio (gas-phase activity) of 1 ppbv of nitric acid over
a pure liquid.
Fig. 6
Molecules per particle of pure ammonium nitrate.
70
|
Environ. Sci.: Atmos.
,2025,
5
,67
81
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environmental Science: Atmospheres
Paper
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
K
p
¼
p
NH
3
p
NHO
3
¼

a
v
NH
3
$
a
v
HONO
2

p

NH
3
$
p

NHO
3

K
p
¼
p

NH
3
p

NHO
3
K
eq
R3
ln
ð
K
eq
R3
Þ¼
ln

K
p

þ
ln

p

NH
3

þ
ln

p

NHO
3

Mozurkewich
1
gives an overall dissociation constant for
solid-vapor (in nb
2
=
10
8
Pa
2
)of
ln

K
p

nb
2

¼
118
:
87

24 084
T

6
:
025 ln
ð
T
Þ
ln

K
p

Pa
2

¼
100
:
45

24 084
T

6
:
025 ln
ð
T
Þ
An e
ff
ective Arrhenius activation energy for the equilibrium
coe
ffi
cient is 18 150 K; it is extremely temperature dependent.
It is important to emphasize again that this overall equilib-
rium is well constrained, and the individual condensed-phase
equilibrium coe
ffi
cient and the vapor pressures only provide
useful insights. Fig. 7 shows that most of the temperature
dependence comes from the ammonium nitrate dissociation,
which spans 18 orders of magnitude over the atmospheric
temperature range (compared w
ith 2 or 3 decades for the vapor
pressures). This is why the potential for relatively rapid growth
(rather than equilibrium ammonium nitrate formation in
larger particles) is con
ned to low temperature regions of the
atmosphere.
4 Microphysics
With the thermodynamics (vapor pressures and equilibrium
coe
ffi
cients) known, we can now consider the microphysics of
collisions. We will use a microphysical formulation grounded in
the kinetic regime of molecular (and particle) collisions, rather
than the bulk
uid continuum regime of di
ff
usion to and from
particles. This is because we are most interested in condensa-
tion to particles with 2
(
d
p
(
100 nm, where the Knudsen
number Kn
T
1 and microphysics are more kinetic than not.
The net condensation
ux of a species with vapor concen-
tration
c
v
i
and saturation concentration
c

i
¼
p

i
=
ð
k
B
T
Þ
to
a suspension of identical particles
p
with number concentration
N
s
p
and (physical) diameter
d
p
depends fundamentally on the
collision cross section,
s
i
,
p
, the collision speed,
s
0
i
;
p
;
and the
thermodynamic excess concentration,
c
xs
i
,
p
.
10
This can also be
expressed in terms of the (e
ff
ective spherical) particle surface
area,
A
s
p
and the (e
ff
ective perpendicular) collision speed of
molecules with that surface,
s
t
i
,
p
, but with several modifying
factors that we de
ne below.
Note that the diameter is the physical diameter, not the
mobility diameter (which is typically
d
mob
x
d
p
+ 0.3 nm).
11
The
ux of the condensing species per unit particle surface
area is
f
v
;
s
i
;
p
¼
s
t
i
;
p
c
xs
i
;
p
¼
s
t
i
;
p
c

i
h
a
v
i

a
s
i
;
p
K
i
;
p
i
¼
s
t
i
;
p
c

i
a
v
i

1

a
s
i
;
p
a
v
i
K
i
;
p

¼
s
t
i
;
p
g
i
;
p
c
v
i
;
g
i
;
p
¼

1

a
s
i
;
p
a
v
i
K
i
;
p

(3)
This depends on the collision speed, the vapor concentra-
tion, and an uptake coe
ffi
cient,
g
i
,
p
, that in turn depends on the
activities in both phases as well as a Kelvin term,
K
i
,
p
for
curvature e
ff
ects over very small particles. The
ux is either
number per unit area or mass per unit area depending on the
concentration units.
The collision speed,
s
t
i
,
p
, is the e
ff
ective line-of-centers
velocity of each vapor toward the particle surface, accounting
for enhancements (
i.e.
van der Waals, vapor size, reduced mass)
and limitations (
i.e.
accommodation, gas-phase di
ff
usion near
the particle). It is derived from the average molecular speed,
s
i
¼
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
8
k
B
T
=
ð
p
m
i
Þ
p
:
van der Waals interactions between the
vapor and particle (charge-dipole, dipole-induced dipole, and
induced dipole-induced dipole) can be represented as an
Fig. 7
Equilibrium constant for ammonium nitrate formation in the
condensed phase.
f
v
;
s
i
;
p
¼
N
s
p

p
4

d
p
þ
d
i

2
a
i
;
p
E
m
i
;
p

|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
cross section
:
s
i
;
p
s
i
;
p
B
i
;
p
|fflfflffl{zfflfflffl}
speed
:
s
0
i
;
p
h
c
v
i

a
s
i
;
p
K
i
;
p
c

i
i
|fflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflffl}
excess conc
:
:
c
xs
i
;
p
¼
N
s
p
p
d
p
2
|fflfflfflffl{zfflfflfflffl}
area
:
A
s
p
E
m
i
;
p
3
i
;
p
e
i
;
p
a
i
;
p
s
i
4
B
i
;
p
|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
collision speed
:
s
t
i
;
p
h
c
v
i

a
s
i
;
p
K
i
;
p
c

i
i
|fflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflffl}
excess conc
:
:
c
xs
i
;
p
|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
flux per unit surface area
:
f
v
;
s
i
;
p
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environ. Sci.: Atmos.
,2025,
5
,67
81 |
71
Paper
Environmental Science: Atmospheres
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
enhancement over hard-sphere collisions,
E
m
i
,
p
=
E
m
(
A
i
,
p
), where
A
i
,
p
is the Hamaker constant in joule.
12,13
The non-zero size of
the vapor is accounted for by
3
i
,
p
=
(
d
i
2
+
d
p
2
)/
d
p
2
, and the non-
in
nite mass of small clusters gives a reduced-mass correction
of
e
i
;
p
¼
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ð
m
i
þ
m
p
Þ
=
m
p
p
:
10
The mass accommodation coe
ffi
-
cient,
a
i
,
p
is likely 1.
14
Gas-phase di
ff
usion limitations emerging
for small Knudsen numbers are given by
B
i
,
p
.
10,15
The excess concentration, sometimes also known as the
thermodynamic driving force, can be written as
c
xs
i
;
p
¼
F
v
;
s
i
;
p
¼
c

i
h
a
v
i

a
s
i
;
p
K
i
;
p
i
(4)
The Kelvin term
K
i
,
p
scales with a
Kelvin diameter
for
decadal change,
d
K10
. We use base-10 because it is a little easier
to think in decades; for particles with
d
p
=
d
K10
, condensed-
phase activity is increased by a factor of 10.
log
10
K
i
;
p
¼
d
K10
d
p
þ
O
d
K10
d
p
2
þ
.
d
K10
¼
0
:
434
d
K
;
d
K
¼
4
s
p
m
i
k
B
T
r
i
(5)
This relates the bulk properties (
s
i
,
r
p
) to the molecular
properties of clusters
via
the correspondence principle
(potentially with higher-order terms). It is perhaps not widely
recognized how small
d
K10
actually is; it is the diameter where
the e
ff
ective saturation concentration increases by 1 decade,
and for typical surface tensions and densities,
d
K10
x
4
5nm.
Still, for
d
p
=
1nm,
K
i
,
p
=
10
4
10
5
, which is why, in classical
nucleation theory, very high supersaturations can be required
to surmount the critical cluster. This is crucial to the activa-
tion behavior observed in ammonium-nitrate condensation in
Fig. 1.
We can express the excess concentration as the product of
the vapor concentration and an uptake coe
ffi
cient,
g
i
,
p
(the
fraction of collisions resulting in uptake), de
ned in eqn (3). We
can re-write this in terms of the saturation ratio (including
curvature),
S
i
,
p
, and the excess saturation ratio,
S
xs
i
,
p
=
S
i
,
p
1.
S
i
;
p
¼
a
v
i
a
s
i
;
p
K
i
;
p
g
i
;
p
¼

1

1
S
i
;
p

¼
S
i
;
p

1
S
i
;
p
¼
S
xs
i
;
p
S
i
;
p
(6)
The uptake coe
ffi
cient is simply the fractional excess satu-
ration ratio over a small, curved particle: the ratio of the
Kelvin
adjusted
suspended-phase activity,
a
s
i
K
i
,
p
, to the equilibrium
suspended-phase activity, which would be
a
v
i
. For example, for
water vapor during cloud droplet activation, the fractional
excess saturation ratio is typically 0.002
0.01, so between 1%
and 0.2% of water vapor collisions are taken up (assuming
a mass accommodation coe
ffi
cient,
a
=
1). If the system is at
equilibrium,
g
i
,
p
=
0, whereas for kinetic uptake on every
collision with the surface,
g
i
,
p
=
1. The uptake coe
ffi
cient is thus
just the fractional excess saturation ratio (at the surface
we
neglect any condensed phase di
ff
usion limitations, which in
any event are irrelevant for this pure model system).
The growth rate,
R
gr
i
,
p
, is proportional to the
ux per unit
surface area,
f
v,s
i
,
p
. Assuming a spherical particle (or a spherical
equivalent
d
p
), the growth rate due to condensation of a given
species is
R
gr
i
;
p
¼
d
d
p
d
t
i
¼
2
V
i
;
p
f
v
;
s
i
;
p

½
c
i
#
volume

3

¼
2
r
i
;
p
f
v
;
s
i
;
p

½
c
i
mass volume

3

(7)
The factor of 2 is because growth extends the diameter by
extending the radius symmetrically about the sphere. The
molecular volume in the particle is
V
i
,
p
=
m
i
,
p
/
r
i
,
p
; we retain the
{
i
,
p
}
species within particle
designation to account for the
form of the species within the particle. Speci
cally in the case
we are considering here, one mole of either ammonia or nitric
acid condensation will result in one mole of ammonium nitrate,
and so the appropriate mass and density are those of ammo-
nium nitrate. More generally, we could write
d
V
i
,
p
as the change
in particle volume arising from condensation of this species;
here we can also treat condensation of water vapor for deli-
quesced and/or aqueous particles.
We can solve eqn (7) for any given parameter, but for the
most part we will either want to know the growth rate for a given
vapor concentration
c
v
i
or mixing ratio
x
v
i
or conversely
nd
either of those in terms of a given growth rate:
R
gr
i
,
p
=
2
V
i
,
p
s
t
i
,
p
g
i
,
p
c
v
i
=
2
V
i
,
p
s
t
i
,
p
g
i
,
p
x
v
i
c
c
M
(8)
We also may want to identify the maximum possible growth
rate, when
g
i
,
p
=
1:
Fig. 8
Collision speed
s
t
i
,
p
for NH
3
and HONO
2
vs.
the particle diam-
eter with
A
HONO
2
=
6
×
10
20
joule and
A
NH
3
=
0 joule. Dashed hori-
zontal lines are asymptotic values equal to one-quarter of the average
gas speed. This is the so-called
kinetic
limit of nominally constant
growth; for
d
p
x
10 nm the NH
3
speed exhibits a modest plateau near
the
kinetic
value. For very small particles, the size of the molecule,
(small) mass of the particle, and the interaction potential between
them all enhance the speed in the free-molecular or
collision
limit.
72
|
Environ. Sci.: Atmos.
,2025,
5
,67
81
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environmental Science: Atmospheres
Paper
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
R
max
i
,
p
=
2
V
i
,
p
s
t
i
,
p
c
v
i
=
2
V
i
,
p
s
t
i
,
p
x
v
i
c
v
M
(9)
The concentration of all (air) molecules is
c
v
M
.
Fig. 8 shows the collision speeds for both NH
3
and HONO
2
,
including the
kinetic regime
asymptote
s
t
;
kin
i
¼
s
i
=
4
:
This
assumes Hamaker constants of
A
HONO
2
=
6
×
10
20
joule and
A
NH
3
=
0 joule, because HONO
2
has a larger dipole moment
than NH
3
. Values are based on growth rate observations for
H
2
SO
4
/
(NH
4
)
2
$
SO
4
(ref. 15) because the dipole moment of
HONO
2
vapor is similar to the dipole moment of H
2
SO
4
vapor
and the relative values appropriate for (NH
4
)
$
NO
3
vs.
(NH
4
)
2
$
SO
4
are not known. The exact values are not overly
important to this example, but the di
ff
erence between a highly
polar but relatively heavy acid and a non-polar but light base is
interesting and illustrative. As Fig. 8 shows, NH
3
has a speed
advantage above 5 nm, but for very small particles the van der
Waals e
ff
ect (as well as the molecular size and reduced mass
e
ff
ects) give HONO
2
a speed advantage. Growth by ammonium
nitrate formation ultimately means a 1 : 1 uptake frequency for
NH
3
and HONO
2
, and this collision speed transition means that
in some circumstances a system can go from being NH
3
limited
for the smallest particles to HONO
2
limited for larger particles.
5 Semivolatile condensation
Let us consider a
ux balance for ammonium nitrate uptake.
For the sake of simplicity and generality of (mono)acid
base
condensation we shall refer to HONO
2
as A and NH
3
as B.
5.1 Setup
Because of the 1 : 1 stoichiometry enforced by the reactive
uptake, there will be a stoichiometric, molar,
ux balance:
f
v
;
s
A
;
p
¼
f
v
;
s
B
;
p
s
t
B
;
p
c
v
B
g
B
;
p
¼
s
t
A
;
p
c
v
A
g
A
;
p
g
B
;
p
¼
s
t
A
;
p
c
v
A
s
t
B
;
p
c
v
B
g
A
;
p
g
B
;
p
¼
r
n
AB
g
A
;
p
r
n
AB
¼
s
t
A
;
p
c
v
A
s
t
B
;
p
c
v
B
(10)
Here we de
ne a collision ratio of A and B with the surface,
r
n
AB
. This is the concentration ratio of the two species scaled by
the collision speed for each molecule seen in Fig. 8. Eqn (10) is
symmetrical, but it is easier to think of the case where
r
n
AB
$
1so
that A is in excess and B is the rate-limiting contributor to
condensational growth. In the limit of
r
n
AB
[
1, there is a large
excess of A and we expect A to be equilibrated (
a
s
A,
p
K
A,
p
/
a
v
A
), so
g
A,
p
/
0.
When the collision ratio
r
n
AB
=
1, the system is at the
equi-
collision
point where with stoichiometric uptake we must have
g
A
=
g
B
. For large excess saturation,
g
=
1. In the real case we
are considering, NH
3
(B) is about twice as fast as HONO
2
(A), so
the equi-collision point occurs when
c
HONO
2
x
2
c
NH
3
(A : B
x
2).
The maximal growth rates, for
g
=
1, form a chevron meeting
along this
equi-collision
line, as shown in Fig. 9 for
d
p
=
10 nm. This shows that a few 100 pptv of either NH
3
or HONO
2
have the potential to drive particle growth at
s
gr
> 100 nm h
1
,
provided that there is enough driving force (supersaturation) for
this to actually occur. As a point of reference, growth rates
between 1 and 10 nm h
1
are typical of the continental
boundary layer when conditions permit accurate determination
of aerosol growth.
16,17
However, we also have the ammonium nitrate equilibrium to
consider. Here we assume rapid equilibration of the salt with
the undissociated molecules in the suspended particle phase
(slow particle phase kinetics is a completely di
ff
erent limit,
mostly controlled by volatility). This imposes a balance on the
activity product in the particles.
a
s
A
a
s
B
¼
1
K
eq
AB
or
a
s
A
¼
1
a
s
B
K
eq
AB
(11)
The overall saturation of vapors toward the particle phase
can be written as
S
AB
;
p
¼
a
v
A
a
v
B
a
s
A
a
s
B
K
A
;
p
K
B
;
p
¼
a
v
A
a
v
B
K
eq
AB
K
A
;
p
K
B
;
p
¼
a
v
A
a
v
B
K
eq
AB
K
2
AB
;
p
(12)
This is the product of the vapor activities (saturation ratios)
and suspended particle phase activities (including curvature)
and it is the fundamental imposed thermodynamic driver for
growth. The two Kelvin terms multiply, which we represent as
the square based on the average of the two Kelvin diameters,
d
K10,AB
=
0.5
×
(
d
K10,A
+
d
K10,B
).
The overall equilibrium between vapor concentrations and
the salt is what is best constrained by thermodynamic data and
calculations (without the Kelvin e
ff
ect).
1
As Fig. 10 shows, for
Fig. 9
Maximum growth rates of 10 nm particles at 263 K for
ammonium nitrate formation. Growth rates are shown as contour
labels. The maximum depends sharply on the limiting vapor, with
ammonia limitation causing horizontal contours at relatively high nitric
acid and nitric acid limitation causing vertical contours at relatively
high ammonia. The spine is o
ff
set from 1 : 1 based on the relative vapor
speeds and depends on particle size.
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environ. Sci.: Atmos.
,2025,
5
,67
81 |
73
Paper
Environmental Science: Atmospheres
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online
sub-freezing conditions (263 K) the saturation isopleths over
a
at surface cut through typical ambient (urban) mixing ratios
of NH
3
and HONO
2
.
We have two constraints on the actual growth rate. One is the
maximum possible growth rate shown in Fig. 9, and the other is
the saturation ratio,
S
AB,
p
, shown in Fig. 10, both at 263 K. These
two combine to give the overall constraints to growth, as shown
in Fig. 11. We expect the actual growth rate to be at the
maximum value for
S
AB,
p
[
1, but obviously much lower for
S
AB,
p
T
1, where the natural quantity is the excess saturation,
S
AB,
p
1. Of course, for
S
AB,
p
< 1 there is no growth, and the
dynamics turn to evaporation.
The collision rates and thus the maximum growth rates are
only a weak function of temperature, whereas the saturation
ratios (at a
xed concentration or mixing ratio) are extremely
strong functions of temperature. The saturation limits will thus
sweep from the lower le
of Fig. 11 toward the upper right as
temperature increases.
5.2 Activation diameter
Before anything grows, the growth must start. This is activation,
and it is characterized by an activation diameter,
d
act
. Unac-
tivated particles are made of something, and here we assume
that ammonium nitrate (acid and base) does not interact with
the unactivated particles beyond the Kelvin e
ff
ect. The activa-
tion diameter simply depends on supersaturation and so strictly
scales with the square of the gas-phase concentrations (and
thus activities) of the acid and base.
S
AB
;
p
¼
a
v
A
a
v
B
K
eq
AB
K
2
AB
;
p
¼
1
1
¼
S

AB
K
2
AB
;
p
;
S

AB
¼
a
v
A
a
v
B
K
eq
AB
S

AB
¼
10
ð
2
d
K10
;
AB
=
d
act
Þ
d
act
¼
2
d
K10
;
AB
log
10

S

AB

(13)
As stated above we treat the Kelvin diameter(s) as an
adjustable parameter given th
e very small number of mole-
cules in the activating clusters. However, we expect the values
to lie somewhere between those of water (
d
K10,H
2
O
x
1nm)
and organics (
d
K10,org
x
6nm).Ifthegrowingparticleis
a quasi aqueous solution with
s
x
75 mN m
1
,wemight
Fig. 10
Saturation ratio over a pure ammonium nitrate crystalline
surface at 263 K
vs.
nitric acid and ammonia mixing ratios. The solid
green diagonal shows
S
=
1. Regions to the lower left are undersatu-
rated and particles will evaporate; those to the upper right, with
saturation ratios indicated on the dashed contours, are super saturated
and particles will grow. The saturation threshold moves toward the
lower left as temperature falls and towards the upper right as it
increases.
Fig. 11
Phase space combining maximum growth rates and the
saturation ratio over an ammonium nitrate surface at 263 K. Only
systems in the saturated region to the upper right of the solid green
saturation contour (
S
=
1) will grow, and growth near
S
=
1 will be
(much) slower than the maximum possible values.
Fig. 12
Activation diameters at 263 K for solid ammonium nitrate with
d
K10
=
3 nm. Diameters are shown with contour labels, and
d
act
=
N
at
the saturation line
S
=
1. To the lower left of this saturation line
particles will not activate. Unlike the growth rates, the activation
diameters depend exponentially on the saturation ratios; for example,
the 5 nm activation diagonal from upper right to lower left corre-
sponds to
S
x
20, whereas the smallest diameter shown (2 nm)
corresponds to
S
x
1000.
74
|
Environ. Sci.: Atmos.
,2025,
5
,67
81
© 2025 The Author(s). Published by the Royal Society of Chemistry
Environmental Science: Atmospheres
Paper
Open Access Article. Published on 21 October 2024. Downloaded on 2/21/2025 12:34:07 AM.
This article is licensed under a
Creative Commons Attribution 3.0 Unported Licence.
View Article Online