Monte Carlo Methods for Estimating Interfacial Free Energies and Line Tensions

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given

Evaporation of sulfate aerosols at low relative humidity

Evaporation of sulfuric acid from particles can be important in the atmospheres of Earth and Venus. However, the equilibrium constant for the dissociation of H2SO4 to bisulfate ions, which is the one of the fundamental parameters controlling the evaporation of sulfur particles, is not well constrained. In this study we explore the volatility of sulfate particles at very low relative humidity. We measured the evaporation of sulfur particles versus temperature and relative humidity in the CLOUD chamber at CERN. We modelled the observed sulfur particle shrinkage with the ADCHAM model. Based on our model results, we conclude that the sulfur particle shrinkage is mainly governed by H2SO4 and potentially to some extent by SO3 evaporation. We found that the equilibrium constants for the dissociation of H2SO4 to HSO4-(KH2SO4) and the dehydration of H2SO4 to SO3 ((x) K-SO3) are K H2SO4 = 2-4 x 10(9) mol kg(-1) and (x) K SO3 >= 1.4 x 10(10) at 288.8 +/- 5 K.Peer reviewe

New particle formation in the sulfuric acid-dimethylamine-water system: reevaluation of CLOUD chamber measurements and comparison to an aerosol nucleation and growth model

A recent CLOUD (Cosmics Leaving OUtdoor Droplets) chamber study showed that sulfuric acid and dimethylamine produce new aerosols very efficiently and yield particle formation rates that are compatible with boundary layer observations. These previously published new particle formation (NPF) rates are reanalyzed in the present study with an advanced method. The results show that the NPF rates at 1.7 nm are more than a factor of 10 faster than previously published due to earlier approximations in correcting particle measurements made at a larger detection threshold. The revised NPF rates agree almost perfectly with calculated rates from a kinetic aerosol model at different sizes (1.7 and 4.3 nm mobility diameter). In addition, modeled and measured size distributions show good agreement over a wide range of sizes (up to ca. 30 nm). Furthermore, the aerosol model is modified such that evaporation rates for some clusters can be taken into account; these evaporation rates were previously published from a flow tube study. Using this model, the findings from the present study and the flow tube experiment can be brought into good agreement for the high base-to-acid ratios (similar to 100) relevant for this study. This confirms that nucleation proceeds at rates that are compatible with collision-controlled (a.k.a. kinetically controlled) NPF for the conditions during the CLOUD7 experiment (278 K, 38% relative humidity, sulfuric acid concentration between 1 x 10(6) and 3 x 10(7) cm(-3), and dimethylamine mixing ratio of similar to 40 pptv, i.e., 1 x 10(9) cm(-3)).Peer reviewe

The role of low-volatility organic compounds in initial particle growth in the atmosphere

About half of present-day cloud condensation nuclei originate from atmospheric nucleation, frequently appearing as a burst of new particles near midday1. Atmospheric observations show that the growth rate of new particles often accelerates when the diameter of the particles is between one and ten nanometres2, 3. In this critical size range, new particles are most likely to be lost by coagulation with pre-existing particles4, thereby failing to form new cloud condensation nuclei that are typically 50 to 100 nanometres across. Sulfuric acid vapour is often involved in nucleation but is too scarce to explain most subsequent growth5, 6, leaving organic vapours as the most plausible alternative, at least in the planetary boundary layer7, 8, 9, 10. Although recent studies11, 12, 13 predict that low-volatility organic vapours contribute during initial growth, direct evidence has been lacking. The accelerating growth may result from increased photolytic production of condensable organic species in the afternoon2, and the presence of a possible Kelvin (curvature) effect, which inhibits organic vapour condensation on the smallest particles (the nano-Köhler theory)2, 14, has so far remained ambiguous. Here we present experiments performed in a large chamber under atmospheric conditions that investigate the role of organic vapours in the initial growth of nucleated organic particles in the absence of inorganic acids and bases such as sulfuric acid or ammonia and amines, respectively. Using data from the same set of experiments, it has been shown15 that organic vapours alone can drive nucleation. We focus on the growth of nucleated particles and find that the organic vapours that drive initial growth have extremely low volatilities (saturation concentration less than 10−4.5 micrograms per cubic metre). As the particles increase in size and the Kelvin barrier falls, subsequent growth is primarily due to more abundant organic vapours of slightly higher volatility (saturation concentrations of 10−4.5 to 10−0.5 micrograms per cubic metre). We present a particle growth model that quantitatively reproduces our measurements. Furthermore, we implement a parameterization of the first steps of growth in a global aerosol model and find that concentrations of atmospheric cloud concentration nuclei can change substantially in response, that is, by up to 50 per cent in comparison with previously assumed growth rate parameterizations

Global atmospheric particle formation from CERN CLOUD measurements

Fundamental questions remain about the origin of newly formed atmospheric aerosol particles because data from laboratory measurements have been insufficient to build global models. In contrast, gas-phase chemistry models have been based on laboratory kinetics measurements for decades. Here we build a global model of aerosol formation using extensive laboratory-measured nucleation rates involving sulfuric acid, ammonia, ions and organic compounds. The simulations and a comparison with atmospheric observations show that nearly all nucleation throughout the present-day atmosphere involves ammonia or biogenic organic compounds in addition to sulfuric acid. A significant fraction of nucleation involves ions, but the relatively weak dependence on ion concentrations indicates that for the processes studied variations in cosmic ray intensity do not significantly affect climate via nucleation in the present-day atmosphere