Atmospheric Chemistry and Physics

An Interactive Open Access Journal of the European Geosciences Union

| Contact |
Atmos. Chem. Phys. Discuss., 11, 25263-25295, 2011
www.atmos-chem-phys-discuss.net/11/25263/2011/
doi:10.5194/acpd-11-25263-2011
© Author(s) 2011. This work is distributed
under the Creative Commons Attribution 3.0 License.


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

M. J. McGrath1,2, T. Olenius1, I. K. Ortega1, V. Loukonen1, P. Paasonen1, T. Kurtén1, M. Kulmala1, and H. Vehkamäki1
1Department of Physics, P.O. Box 64, University of Helsinki, 00014, Helsinki, Finland
2Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan

Abstract. 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 ode15s 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, the temperature and sticking probabilities both have a large impact on 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, and monomer concentrations can all have significant effects. 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−3, which is probably an effect caused by studying only relatively small cluster sizes.

Citation: McGrath, M. J., Olenius, T., Ortega, I. K., Loukonen, V., Paasonen, P., Kurtén, T., Kulmala, M., and Vehkamäki, H.: Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations, Atmos. Chem. Phys. Discuss., 11, 25263-25295, doi:10.5194/acpd-11-25263-2011, 2011.
 
Search ACPD
Discussion Paper
    XML
    Citation
    Please cite the
    Final Revised Paper
    Final Revised Paper
    Published on 2012-03-02
    Share