Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations

The Atmospheric Cluster Dynamics Code (ACDC) is presented and explored. This program was created to study the first steps of atmospheric new particle formation by examining the formation of molecular clusters from atmospherically relevant molecules. The program models the cluster kinetics by explici...

Full description

Bibliographic Details
Main Authors: M. J. McGrath, T. Olenius, I. K. Ortega, V. Loukonen, P. Paasonen, T. Kurtén, M. Kulmala, H. Vehkamäki
Format: Article
Language:English
Published: Copernicus Publications 2012-03-01
Series:Atmospheric Chemistry and Physics
Online Access:http://www.atmos-chem-phys.net/12/2345/2012/acp-12-2345-2012.pdf
_version_ 1818199998683152384
author M. J. McGrath
T. Olenius
I. K. Ortega
V. Loukonen
P. Paasonen
T. Kurtén
M. Kulmala
H. Vehkamäki
author_facet M. J. McGrath
T. Olenius
I. K. Ortega
V. Loukonen
P. Paasonen
T. Kurtén
M. Kulmala
H. Vehkamäki
author_sort M. J. McGrath
collection DOAJ
description The Atmospheric Cluster Dynamics Code (ACDC) is presented and explored. This program was created to study the first steps of atmospheric new particle formation by examining the formation of molecular clusters from atmospherically relevant molecules. The program models the cluster kinetics by explicit solution of the birth–death equations, using an efficient computer script for their generation and the MATLAB <tt>ode15s</tt> routine for their solution. Through the use of evaporation rate coefficients derived from formation free energies calculated by quantum chemical methods for clusters containing dimethylamine or ammonia and sulphuric acid, we have explored the effect of changing various parameters at atmospherically relevant monomer concentrations. We have included in our model clusters with 0–4 base molecules and 0–4 sulfuric acid molecules for which we have commensurable quantum chemical data. The tests demonstrate that large effects can be seen for even small changes in different parameters, due to the non-linearity of the system. In particular, changing the temperature had a significant impact on the steady-state concentrations of all clusters, while the boundary effects (allowing clusters to grow to sizes beyond the largest cluster that the code keeps track of, or forbidding such processes), coagulation sink terms, non-monomer collisions, sticking probabilities and monomer concentrations did not show as large effects under the conditions studied. Removal of coagulation sink terms prevented the system from reaching the steady state when all the initial cluster concentrations were set to the default value of 1 m<sup>−3</sup>, which is probably an effect caused by studying only relatively small cluster sizes.
first_indexed 2024-12-12T02:30:40Z
format Article
id doaj.art-8327aad53c6147d0bf2ace4bcc6bacf7
institution Directory Open Access Journal
issn 1680-7316
1680-7324
language English
last_indexed 2024-12-12T02:30:40Z
publishDate 2012-03-01
publisher Copernicus Publications
record_format Article
series Atmospheric Chemistry and Physics
spelling doaj.art-8327aad53c6147d0bf2ace4bcc6bacf72022-12-22T00:41:25ZengCopernicus PublicationsAtmospheric Chemistry and Physics1680-73161680-73242012-03-011252345235510.5194/acp-12-2345-2012Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equationsM. J. McGrathT. OleniusI. K. OrtegaV. LoukonenP. PaasonenT. KurténM. KulmalaH. VehkamäkiThe Atmospheric Cluster Dynamics Code (ACDC) is presented and explored. This program was created to study the first steps of atmospheric new particle formation by examining the formation of molecular clusters from atmospherically relevant molecules. The program models the cluster kinetics by explicit solution of the birth–death equations, using an efficient computer script for their generation and the MATLAB <tt>ode15s</tt> routine for their solution. Through the use of evaporation rate coefficients derived from formation free energies calculated by quantum chemical methods for clusters containing dimethylamine or ammonia and sulphuric acid, we have explored the effect of changing various parameters at atmospherically relevant monomer concentrations. We have included in our model clusters with 0–4 base molecules and 0–4 sulfuric acid molecules for which we have commensurable quantum chemical data. The tests demonstrate that large effects can be seen for even small changes in different parameters, due to the non-linearity of the system. In particular, changing the temperature had a significant impact on the steady-state concentrations of all clusters, while the boundary effects (allowing clusters to grow to sizes beyond the largest cluster that the code keeps track of, or forbidding such processes), coagulation sink terms, non-monomer collisions, sticking probabilities and monomer concentrations did not show as large effects under the conditions studied. Removal of coagulation sink terms prevented the system from reaching the steady state when all the initial cluster concentrations were set to the default value of 1 m<sup>−3</sup>, which is probably an effect caused by studying only relatively small cluster sizes.http://www.atmos-chem-phys.net/12/2345/2012/acp-12-2345-2012.pdf
spellingShingle M. J. McGrath
T. Olenius
I. K. Ortega
V. Loukonen
P. Paasonen
T. Kurtén
M. Kulmala
H. Vehkamäki
Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations
Atmospheric Chemistry and Physics
title Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations
title_full Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations
title_fullStr Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations
title_full_unstemmed Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations
title_short Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations
title_sort atmospheric cluster dynamics code a flexible method for solution of the birth death equations
url http://www.atmos-chem-phys.net/12/2345/2012/acp-12-2345-2012.pdf
work_keys_str_mv AT mjmcgrath atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT tolenius atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT ikortega atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT vloukonen atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT ppaasonen atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT tkurten atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT mkulmala atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations
AT hvehkamaki atmosphericclusterdynamicscodeaflexiblemethodforsolutionofthebirthdeathequations