How well can we predict cluster fragmentation inside a mass spectrometer?

Fragmentation of molecular clusters inside mass spectrometers is a significant source of uncertainty in a wide range of chemical applications. We have measured the fragmentation of sulfuric acid clusters driving atmospheric new-particle formation, and developed a novel model, based on first principles calculations, capable of quantitatively predicting the extent of fragmentation.

Modeling on Fragmentation of Clusters inside a Mass Spectrometer.

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.

Deviation from equilibrium conditions in molecular dynamic simulations of homogeneous nucleation.

We present a comparison between Monte Carlo (MC) results for homogeneous vapour-liquid nucleation of Lennard-Jones clusters and previously published values from molecular dynamics (MD) simulations. Both the MC and MD methods sample real cluster configuration distributions. In the MD simulations, the extent of the temperature fluctuation is usually controlled with an artificial thermostat rather than with more realistic carrier gas. In this study, not only a primarily velocity scaling thermostat is considered, but also Nosé-Hoover, Berendsen, and stochastic Langevin thermostat methods are covered. The nucleation rates based on a kinetic scheme and the canonical MC calculation serve as a point of reference since they by definition describe an equilibrated system. The studied temperature range is from T = 0.3 to 0.65 ϵ/k. The kinetic scheme reproduces well the isothermal nucleation rates obtained by Wedekind et al. [J. Chem. Phys. 127, 064501 (2007)] using MD simulations with carrier gas. The nucleation rates obtained by artificially thermostatted MD simulations are consistently lower than the reference nucleation rates based on MC calculations. The discrepancy increases up to several orders of magnitude when the density of the nucleating vapour decreases. At low temperatures, the difference to the MC-based reference nucleation rates in some cases exceeds the maximal nonisothermal effect predicted by classical theory of Feder et al. [Adv. Phys. 15, 111 (1966)].

Effect of Conformers on Free Energies of Atmospheric Complexes.

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 (H2SO4)3(NH3)3(H2O)4 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.

Evaporation rate of nucleating clusters.

The Becker-Döring kinetic scheme is the most frequently used approach to vapor liquid nucleation. In the present study it has been extended so that master equations for all cluster configurations are included into consideration. In the Becker-Döring kinetic scheme the nucleation rate is calculated through comparison of the balanced steady state and unbalanced steady state solutions of the set of kinetic equations. It is usually assumed that the balanced steady state produces equilibrium cluster distribution, and the evaporation rates are identical in the balanced and unbalanced steady state cases. In the present study we have shown that the evaporation rates are not identical in the equilibrium and unbalanced steady state cases. The evaporation rate depends on the number of clusters at the limit of the cluster definition. We have shown that the ratio of the number of n-clusters at the limit of the cluster definition to the total number of n-clusters is different in equilibrium and unbalanced steady state cases. This causes difference in evaporation rates for these cases and results in a correction factor to the nucleation rate. According to rough estimation it is 10(-1) by the order of magnitude and can be lower if carrier gas effectively equilibrates the clusters. The developed approach allows one to refine the correction factor with Monte Carlo and molecular dynamic simulations.

Connection between the virial equation of state and physical clusters in a low density vapor.

We carry out Monte Carlo simulations of physical Lennard-Jones and water clusters and show that the number of physical clusters in vapor is directly related to the virial equation of state. This relation holds at temperatures clearly below the critical temperatures, in other words, as long as the cluster-cluster interactions can be neglected--a typical assumption used in theories of nucleation. Above a certain threshold cluster size depending on temperature and interaction potential, the change in cluster work of formation can be calculated analytically with the recently proposed scaling law. The breakdown of the scaling law below the threshold sizes is accurately modeled with the low order virial coefficients. Our results indicate that high order virial coefficients can be analytically calculated from the lower order coefficients when the scaling law for cluster work of formation is valid. The scaling law also allows the calculation of the surface tension and equilibrium vapor density with computationally efficient simulations of physical clusters. Our calculated values are in good agreement with those obtained with other methods. We also present our results for the curvature dependent surface tension of water clusters.

Origin of the failure of classical nucleation theory: incorrect description of the smallest clusters.

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.

Comparison between the classical theory predictions and molecular simulation results for heterogeneous nucleation of argon.

We have performed Monte Carlo simulations of homogeneous and heterogeneous nucleations of Lennard-Jones argon clusters. The simulation results were interpreted using the major concept posing a difference between the homogeneous and heterogeneous classical nucleation theories-the contact parameter. Our results show that the multiplication concept of the classical heterogeneous nucleation theory describes the cluster-substrate interaction surprisingly well even for small molecular clusters. However, in the case of argon nucleating on a rigid monolayer of fcc(111) substrate at T=60 K, the argon-substrate atom interaction being approximately one-third as strong as the argon-argon interaction, the use of the classical theory concept results in an underestimation of the heterogeneous nucleation rate by two to three orders of magnitude even for large clusters. The main contribution to this discrepancy is induced by the failure of the classical theory of homogeneous nucleation to predict the energy involved in bringing one molecule from the vapor to the cluster for clusters containing less than approximately 15 molecules.

Analysis of nucleation ability of cluster configurations with Monte Carlo simulations of argon.

We determine the nucleation ability of argon clusters from Monte Carlo simulations. The nucleation rate appears to be defined by a sole characteristic of the clusters, namely, the stability. The stability is calculated as the ratio of grand canonical growth and decay rates and can be assigned to individual cluster configurations. We study the connection between the stability of the cluster configurations and their volume and total potential energy. Neither the potential energy nor the volume of a cluster configuration has a clear relation to its stability, and thus to the nucleation ability. On the other hand, we show that it is possible to use a specific volume for each cluster size to calculate the work of the cluster formation. These clusters with a unique volume have the same average stability as the full set of clusters. Our simulation method allows us to study the effect of possible deviations from equilibrium in the cluster configuration distributions. We argue that the nucleation process itself can produce a source for such a deviation. We show that even a small deviation from equilibrium in the cluster configuration distribution can lead to a dramatic deceleration of the nucleation rate. Although our simulations may overestimate the magnitude of the effect, they give qualitative estimates for its importance.

The molecular approach to heterogeneous nucleation.

A molecular approach to heterogeneous nucleation has been developed. The expressions for the equilibrium cluster distribution, the reversible work of the cluster formation, and the nucleation rate have been derived. Two separate statements for the work of formation were formulated. If the equilibrium cluster distribution is normalized on the monomer concentration near the substrate surface, the reversible work of formation is expressed by DeltaG(het) (I) = (F(n) (het)-F(n) (hom))-(F(1) (het)-F(1) (hom)) + DeltaG(hom) where F(n) (het) and F(n) (hom) are the Helmholtz free energies of a cluster interacting with a substrate and a cluster not interacting with the substrate, respectively. If the equilibrium cluster distribution is normalized on the monomer concentration far from the substrate surface, the work of cluster formation is given by DeltaG(het) (II) = (F(n) (het)-F(n) (hom)) + DeltaG(hom). The former expression corresponds to the approach of the classical heterogeneous nucleation theory. The cluster partition function appears to be dependent on the location of a virtual plane, which separates the volume, where the interaction of the clusters with the substrate is effective from the one where interaction is negligible. Our Monte Carlo simulations have shown that the dependence is rather weak and thus the location of the plane is not very important. According to the simulations the variation of the plane position in the range from 20 to 50 Angstroms does not lead to a considerable change of the heterogeneous nucleation rate.

Monte Carlo simulations of critical cluster sizes and nucleation rates of water.

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.