Aerosol Science and Technology
ISSN: 0278-6826 (Print) 1521-7388 (Online) Journal homepage: https://www.tandfonline.com/loi/uast20
Stochastic effects in H2SO4-H2O cluster growth
Christoph Köhn, Martin Bødker Enghoff & Henrik Svensmark
To cite this article: Christoph Köhn, Martin Bødker Enghoff & Henrik Svensmark (2020): Stochastic effects in H2SO4-H2O cluster growth, Aerosol Science and Technology, DOI: 10.1080/02786826.2020.1755012 To link to this article: https://doi.org/10.1080/02786826.2020.1755012
Accepted author version posted online: 15 Apr 2020. Published online: 29 Apr 2020.
Submit your article to this journal
Article views: 10
View related articles
View Crossmark data
Full Terms & Conditions of access and use can be found at https://www.tandfonline.com/action/journalInformation?journalCode=uast20 AEROSOL SCIENCE AND TECHNOLOGY https://doi.org/10.1080/02786826.2020.1755012
Stochastic effects in H2SO4-H2O cluster growth
Christoph Kohn€ , Martin Bødker Enghoff , and Henrik Svensmark National Space Institute (DTU Space), Technical University of Denmark, Kgs Lyngby, Denmark
ABSTRACT ARTICLE HISTORY The nucleation of sulfuric acid-water clusters plays a significant role in the formation of Received 8 January 2020 aerosols. Based on a recently developed particle Monte Carlo (MC) Code, we analyze how Accepted 6 April 2020 the growth of sulfuric acid-water clusters is influenced by stochastic fluctuations. We here EDITOR consider samples of H2SO4-H2O clusters at T ¼ 200 K with a relative humidity of 50%, with – – – Jim Smith particle concentrations between 105 and 107 cm 3 in volumes between 10 6 and 10 2 cm3. We present the temporal evolution of the formation rate and of the size distribution as well as growth rates and the onset time of the nucleation above a given cluster size with and without constant production of new monomers. Clear evidence is revealed by the MC code that fluctuations result in a faster growth rate of the smallest clusters compared to deter- ministic continuum models that do not contain the stochastic effects. The faster growth of small clusters in turn influences the growth of larger clusters. Depending on the volume size, the onset time for clusters larger than 0.85 nm varies between 1000 s and 20,000 s for n ¼ 105 cm–3 and between 10 s and 100 s for n ¼ 107 cm–3.
1. Introduction respectively, and C the monomer concentration. Recently it was shown that the inclusion of stochastic On a global scale approximately half of cloud condensa- fluctuations to the size distribution can significantly tion nuclei (CCN)—the seeds for cloud droplet forma- increase the calculated growth rates of sub 5- tion—originate from nucleation and condensation of nanometer sized aerosols (Olenius et al. 2018). low-volatility gases in the modern day atmosphere Stochastic fluctuations can be understood as the ran- (Merikanto et al. 2009;YuandLuo2009); this number dom collisions between molecular clusters (and corre- increases to about two thirds for the pre-industrial atmos- sponding evaporations) which happen outside of the phere (Gordon et al. 2017). Since newly formed aerosols average growth of the system given by the term bC are small in size ( 1 nm) and thus have a high mobility, c: These events are not included in the continuous they have a high risk of being scavenged by larger particles GDE (Equation (1)), which forms the foundation for (Pierce and Adams 2007) before they can grow to sizes many aerosol growth models (Olenius et al. 2013; relevant for cloud formation. The faster the aerosols can Trostl€ et al. 2016). grow to larger sizes the higher their probability of survival Not including the stochastic term can thus have so the dynamics of both the formation and growth of new the risk of underestimating the potential of nucleated aerosols are therefore important to understand the global aerosols to form CCN, but this is still a relatively CCN concentration (Pierce and Adams 2009). unexplored area. This work is the first application of a Traditionally the dynamics of aerosol populations recently developed 3D particle Monte Carlo model, are described by the general dynamic equation (GDE). tracking individual clusters. Such a model is capable The continuous equation only treating condensation of capturing the individual behavior of each cluster and evaporation (thus neglecting coagulation, losses to and is therefore suitable to investigate the effects of other processes, and production of new monomers) is: stochastic fluctuations on a developing aerosol popula- @cði tÞ @ , ¼ ½ ðbðiÞC cðiÞÞcði tÞ tion whilst traditional models such as the classical @t @i , (1) thermodynamic nucleation theory (Hamill et al. 1982), where c is the concentration as a function of size (i), kinetic models (Pirjola and Kulmala 1998; Lovejoy, b and c the condensation and evaporation constants Curtius, and Froyd 2004;Yu2006), nucleation theory
CONTACT Christoph Kohn€ [email protected] National Space Institute (DTU Space), Technical University of Denmark, Elektrovej 328, 2800 Kgs Lyngby, Denmark. ß 2020 American Association for Aerosol Research 2 C. KÖHN ET AL.
– (Vehkam€aki et al. 2002), or Atmospheric Cluster 107 cm 3 lie in the order of 10 2 10 3 s. With these Dynamics (McGrath et al. 2012) only capture macro- time steps, the diffusive length scale varies from 42 scopic or averaged quantities and are therefore not and 134 lm. Since the mean free path of sulfuric acid suited to deal with the stochastics of particle growth. molecules in air is approximately 10–60 nm (Seinfeld The model used was previously described in Kohn,€ and Pandis 2006), hence significantly smaller than Enghoff, and Svensmark (2018) with modifications as 1 lm, the justification of the diffusion approach in noted in Section 2 where we also describe the setup this study is justified. and simulation runs. In Section 3, we present results After every time step, we check whether clusters on the temporal evolution of the size distribution and evaporate or whether two clusters collide, and if so, its dependence on simulation parameters such as determine the mass and volume of the merged cluster. monomer density and simulation volume. Modeled Throughout this article, we connect the particle i growth rates are compared with those expected from size to the number of H2SO4 molecules in the cluster the continuous GDE and we analyze the condensation (Yu 2005). In this sense, the particle radius of clusters rates in an attempt to determine the critical cluster with i monomers is sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi size. Finally, we summarize and conclude our study in mðiÞ Section 4. In addition, we describe details about the Rði TÞ¼ 3 3 , p.ði TÞ (6) parallelization of our computational model, which is 4 , encoded in Python, in Appendix A. Implicitly including the size dependent amount of
H2O molecules found at a relative humidity of approx. 2. Modeling 50%. The mass and the density for temperatures between 233 and 323 K are given as in Vehkam€aki We model the motion and clustering of hydrated sul- et al. (2002): furic acid clusters as described in detail in Kohn,€ kg kg Enghoff, and Svensmark (2018) where a cluster 0:098 0:018 i i MFðiÞ mðiÞ¼ mol i þ mol kg (H2SO4)i-(H2O)j consists of H2SO4 molecules and NA NA MFðiÞ implicitly includes a size dependent amount of H2O (7) molecules. Here we recapitulate the main features of kg our simulations: We trace individual particles of given .ði, TÞ ¼½0:7681724 þ 2:184714 MFðiÞ cm3 size and update their position through pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi þ 7:163002 MFðiÞ2 44:31447 MFðiÞ3 rðt þ DtÞ¼rðtÞþ 2DðRÞDtG (2) þ 88:75606 MFðiÞ4 with the size-dependent diffusion coefficient D(R) and 75:73729 MFðiÞ5 þ 23:43228 MFðiÞ6 the time step Dt þ T½ K ð0:001808255 0:009294656 MFðiÞ 2 3 2 0:03742147 MFðiÞ þ 0:2565321 MFðiÞ Rð1Þ DðRÞ¼D 4 5 0 R (3) 0:5362872 MFðiÞ þ 0:4857736 MFðiÞ sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 0:1629592 MFðiÞ6ÞþT½ K 2 3 3 2 kBT 1 ð : þ : MFðiÞ D ¼ (4) 0 000003478524 0 00001335867 0 p3mð Þ 2 3 1 4PRð1Þ þ 0:00005195706 MFðiÞ2 0:0003717636 1 MFðiÞ3 þ 0:0007990811 MFðiÞ4 DtՇ pffiffiffiffiffi (5) 3 2 : MFðiÞ5 þ : 2D0 n 0 000745806 0 000258139 1 for molecules of size R where Rð1Þ¼0:329 nm and MFðiÞ6Þ 1000 mð Þ¼ : 25 1 2 0033 10 kg are the initial size and the (8) mass of single sulfuric acid molecules without any 23 –1 23 with the Avogadro constant NA 6:022 10 mol attached water molecules. kB 1:38 10 J/K is the Boltzmann constant and P ¼ 1 bar the ambient pres- and the mole fraction 0:2097 sure. For T ¼ 200 K the diffusion coefficient is D0 MFðiÞ¼0:4505 i (9) : 7 2 G 8 96 10 m /s. is the Gaussian random number which is a powerlaw fit to Fig. 2b of Yu (2005) for a G ¼ð. / h . / h . hÞ . ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffifficos sin , sin sin , cos with temperature of approx. 300 K for a humidity of 50%. 2 log ðr1Þ, / ¼ 2pr2, h ¼ arccosð2r3 1Þ, ri 2½0, 1Þ: Since we here investigate the nucleation of clusters at The typical time steps for densities between 105 and 200 K only, we multiply Equation (9) with a factor of AEROSOL SCIENCE AND TECHNOLOGY 3