We carry out molecular Monte Carlo simulations of clusters in an imperfect vapor. We show that down to very small cluster sizes, classical nucleation theory built on the liquid drop model can be used very accurately to describe the work required to add a monomer to the cluster. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range. We calculate factors needed to correct the cluster formation work given by the liquid drop model. The corrected work of formation results in nucleation rates in good agreement with recent nucleation experiments on argon and water.

In this article we show how to calculate free energies for atmospherically relevant complexes when multiple conformers and/or isomers are present. We explain why the thermal averaging methods used in several published works are incorrect. On the basis of our two sample cases, the sulfuric acid-pinic acid complex and the (HSO)(NH)(HO) cluster, we provide numerical evidence that the use of these incorrect formulas can result in errors larger than 1 kcal/mol. We recommend that if vibrational frequencies and thus Gibbs free energies of the individual conformers are unavailable, one should not attempt to correct for the presence of multiple conformers and instead use only the global minimum conformers for both reactants and products. On the contrary, if the free energies for the conformers are calculated for both reactants and products, their effect can be accounted for by the statistical mechanical methods presented in this article.

Atmospheric clusters are weakly bound and can fragment inside the measuring instruments, in particular, mass spectrometers. Since the clusters accelerate under electric fields, the fragmentation cannot be described in terms of rate constants under equilibrium conditions. Using basic statistical principles, we have developed a model for fragmentation of clusters moving under an external force. The model describes an energy transfer to the cluster internal modes caused by collisions with residual carrier gas molecules. As soon as enough energy is accumulated in the cluster internal modes, it can fragment. The model can be used for interpreting experimental measurements by atmospheric pressure interface mass spectrometers.

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.