Surfactant Partitioning Dynamics in Freshly Generated Aerosol Droplets.

Aerosol droplets are unique microcompartments with relevance to areas as diverse as materials and chemical synthesis, atmospheric chemistry, and cloud formation. Observations of highly accelerated and unusual chemistry taking place in such droplets have challenged our understanding of chemical kinetics in these microscopic systems. Due to their large surface-area-to-volume ratios, interfacial processes can play a dominant role in governing chemical reactivity and other processes in droplets. Quantitative knowledge about droplet surface properties is required to explain reaction mechanisms and product yields. However, our understanding of the compositions and properties of these dynamic, microscopic interfaces is poor compared to our understanding of bulk processes. Here, we measure the dynamic surface tensions of 14-25 µm radius (11-65 pL) droplets containing a strong surfactant (either sodium dodecyl sulfate or octyl-ß-D-thioglucopyranoside) using a stroboscopic imaging approach, enabling observation of the dynamics of surfactant partitioning to the droplet-air interface on time scales of 10s to 100s of microseconds after droplet generation. The experimental results are interpreted with a state-of-the-art kinetic model accounting for the unique high surface-area-to-volume ratio inherent to aerosol droplets, providing insights into both the surfactant diffusion and adsorption kinetics as well as the time-dependence of the interfacial surfactant concentration. This study demonstrates that microscopic droplet interfaces can take up to many milliseconds to reach equilibrium. Such time scales should be considered when attempting to explain observations of accelerated chemistry in microcompartments.

Effect of the Addition of Diblock Copolymer Nanoparticles on the Evaporation Kinetics and Final Particle Morphology for Drying Aqueous Aerosol Droplets.

A deeper understanding of the key processes that determine the particle morphologies generated during aerosol droplet drying is highly desirable for spray-drying of powdered pharmaceuticals and foods, predicting the properties of atmospheric particles, and monitoring disease transmission. Particle morphologies are affected by the drying kinetics of the evaporating droplets, which are in turn influenced by the composition of the initial droplet as well as the drying conditions. Herein, we use polymerization-induced self-assembly (PISA) to prepare three types of sterically stabilized diblock copolymer nanoparticles comprising the same steric stabilizer block and differing core blocks with z-average diameters ranging from 32 to 238 nm. These well-defined nanoparticles enable a systematic investigation of the effect of the nanoparticle size and composition on the drying kinetics of aqueous aerosol droplets (20-28 µm radius) and the final morphology of the resulting microparticles. A comparative kinetics electrodynamic balance was used to obtain evaporation profiles for 10 examples of nanoparticles at a relative humidity (RH) of 0, 45, or 65%. Nanoparticles comprising the same core block with mean diameters of 32, 79, and 214 nm were used to produce microparticles, which were dried under different RH conditions in a falling droplet column. Scanning electron microscopy was used to examine how the drying kinetics influenced the final microparticle morphology. For dilute droplets, the chemical composition of the nanoparticles had no effect on the evaporation rate. However, employing smaller nanoparticles led to the formation of dried microparticles with a greater degree of buckling.

Dynamics and outcomes of binary collisions of equi-diameter picolitre droplets with identical viscosities.

The dynamics of binary collisions of equi-diameter picolitre droplets with identical viscosities, varying impact speeds and impact angles have been investigated experimentally and compared to collision outcome prediction models. Collisions between pairs of pure water droplets with a viscosity of 0.89 mPa s and pairs of aqueous-sucrose (40% w/w) droplets with a viscosity of 5.17 mPa s were examined. The colliding droplets were ∼38 µm in diameter, which is around ten times smaller than those previously investigated when examining the effect of viscosity on the outcome of binary droplet collisions. Varying the impact speed and angle resulted in different collision outcomes, including coalescence, reflexive separation and stretching separation. The collision outcomes were plotted on two viscosity dependent regime maps. The regime boundaries are generally in agreement with earlier literature for both high and low viscosity cases. The agreement between experiment and theory, for both fluids, gives more confidence in the models tested here to predict collision outcomes for droplets of this size and these viscosities.

A purely kinetic description of the evaporation of water droplets.

The process of water evaporation, although deeply studied, does not enjoy a kinetic description that captures known physics and can be integrated with other detailed processes such as drying of catalytic membranes embedded in vapor-fed devices and chemical reactions in aerosol whose volumes are changing dynamically. In this work, we present a simple, three-step kinetic model for water evaporation that is based on theory and validated by using well-established thermodynamic models of droplet size as a function of time, temperature, and relative humidity as well as data from time-resolved measurements of evaporating droplet size. The kinetic mechanism for evaporation is a combination of two limiting processes occurring in the highly dynamic liquid-vapor interfacial region: direct first order desorption of a single water molecule and desorption resulting from a local fluctuation, described using third order kinetics. The model reproduces data over a range of relative humidities and temperatures only if the interface that separates bulk water from gas phase water has a finite width, consistent with previous experimental and theoretical studies. The influence of droplet cooling during rapid evaporation on the kinetics is discussed; discrepancies between the various models point to the need for additional experimental data to identify their origin.

Drying kinetics and nucleation in evaporating sodium nitrate aerosols.

A quantitative understanding of the evaporative drying kinetics and nucleation rates of aqueous based aerosol droplets is important for a wide range of applications, from atmospheric aerosols to industrial processes such as spray drying. Here, we introduce a numerical model for interpreting measurements of the evaporation rate and phase change of drying free droplets made using a single particle approach. We explore the evaporation of aqueous sodium chloride and sodium nitrate solution droplets. Although the chloride salt is observed to reproducibly crystallize at all drying rates, the nitrate salt solution can lose virtually all of its water content without crystallizing. The latter phenomenon has implications for our understanding of the competition between the drying rate and nucleation kinetics in these two systems. The nucleation model is used in combination with the measurements of crystallization events to infer nucleation rates at varying equilibrium state points, showing that classical nucleation theory provides a good description of the crystallization of the chloride salt but not the nitrate salt solution droplets. The reasons for this difference are considered.

Multi-Solvent Microdroplet Evaporation: Modeling and Measurement of Spray-Drying Kinetics with Inhalable Pharmaceutics.

PURPOSE: Evaporation and particle formation from multi-solvent microdroplets containing solid excipients pertaining to spray-drying of therapeutic agents intended for lung delivery were studied. Various water and ethanol co-solvent systems containing a variety of actives and excipients (beclomethasone, budesonide, leucine, and trehalose) were considered. METHODS: Numerical methods were used to predict the droplet evaporation rates and internal solute transfers, and their results verified and compared with results from two separate experimental setups. In particular, an electrodynamic balance was used to measure the evaporation rates of multicomponent droplets and a monodisperse droplet chain setup collected dried microparticles for further analytical investigations and ultramicroscopy. RESULTS: The numerical results are used to explain the different particle morphologies dried from solutions at different co-solvent compositions. The obtained numerical data clearly show that the two parameters controlling the general morphology of a dried particle, namely the Péclet number and the degree of saturation, can change with time in a multi-solvent droplet. This fact complicates product development for such systems. However, this additional complexity vanishes at what we define as the iso-compositional point, which occurs when the solvent ratios and other composition-dependent properties of the droplet remain constant during evaporation, similar to the azeotrope of such systems during distillation. CONCLUSIONS: Numerical and experimental analysis of multi-solvent systems indicate that spray-drying near the iso-compositional ratio simplifies the design and process development of such systems.

Correction to: Multi-Solvent Microdroplet Evaporation: Modeling and Measurement of Spray-Drying Kinetics with Inhalable Pharmaceutics.

The Publisher regrets having introduced the following errors into the article when performing proof corrections.

Complexity of Measuring and Representing the Hygroscopicity of Mixed Component Aerosol.

The validation of approaches to predict the hygroscopicity of complex mixtures of organic components in aerosol is important for understanding the hygroscopic response of organic aerosol in the atmosphere. We report new measurements of the hygroscopicity of mixtures of dicarboxylic acids and amino acids using a comparative kinetic electrodynamic balance (CK-EDB) approach, inferring the equilibrium water content of the aerosol from close to a saturation relative humidity (100%) down to 80%. We show that the solution densities and refractive indices of the mixtures can be estimated with an accuracy of better than ±2% using the molar refractive index mixing rule and densities and refractive indices for the individual binary organic-aqueous solutions. Further, we show that the often-used mass-, volume-, and mole-weighted mixing rules to estimate the hygroscopicity parameter κ can overestimate the hygroscopic parameter by a factor of as much as 3, highlighting the need to understand the specific nonideal interactions that may arise synergistically in mixtures and cannot be represented by simple models. Indeed, in some extreme cases the hygroscopicity of a multicomponent mixture can be very close to that for the least hygroscopic component. For mixtures of similar components for which no additional synergistic interactions need be considered, the hygroscopicity of the mixed component aerosol can be estimated with high accuracy from the hygroscopic response of the binary aqueous-organic aerosol. In conclusion, we suggest that the hygroscopicity of multicomponent organic aerosol can be highly nonadditive and that simple correlations of hygroscopicity with composition may often misrepresent the level of complexity essential to interpreting aerosol hygroscopicity.

Surface Tensions of Picoliter Droplets with Sub-Millisecond Surface Age.

Aerosols are key components of the atmosphere and play important roles in many industrial processes. Because aerosol particles have high surface-to-volume ratios, their surface properties are especially important. However, direct measurement of the surface properties of aerosol particles is challenging. In this work, we describe an approach to measure the surface tension of picoliter volume droplets with surface age <1 ms by resolving their dynamic oscillations in shape immediately after ejection from a microdroplet dispenser. Droplet shape oscillations are monitored by highly time-resolved (500 ns) stroboscopic imaging, and droplet surface tension is accurately retrieved across a wide range of droplet sizes (10-25 µm radius) and surface ages (down to ∼100 µs). The approach is validated for droplets containing sodium chloride, glutaric acid, and water, which all show no variation in surface tension with surface age. Experimental results from the microdroplet dispenser approach are compared to complementary surface tension measurements of 5-10 µm radius droplets with aged surfaces using a holographic optical tweezers approach and predictions of surface tension using a statistical thermodynamic model. These approaches combined will allow investigation of droplet surface tension across a wide range of droplet sizes, compositions, and surface ages.

The evaporation kinetics of pure water droplets at varying drying rates and the use of evaporation rates to infer the gas phase relative humidity.

Numerous analytical models have been applied to describe the evaporation/condensation kinetics of volatile components from aerosol particles for use in many applications. However, the applicability of these models for treating cases that lead to substantial and rapid changes in particle temperature due to, for example, evaporative cooling remain to be compared with measurements. We consider three typical treatments, comparing predictions of the evaporation rates of pure water droplets over a wide range in gas phase relative humidity (RH) and exploring the sensitivity of the predictions to uncertainties in the thermophysical gas and condensed-phase parameters. We also compare predictions from the three treatments to measurements of the evaporation rates of pure water droplets with varying RH using an electrodynamic balance (EDB), concluding that only two of the model treatments are sufficiently able to account for the level of evaporative cooling (typically as high as 12 K). Finally, we show that the RH can be inferred accurately from the evaporation rate of pure water droplets over the full range in accessible RH and comparison with the model predictions (within absolute uncertainties of 2.5% RH over the range 20% to 95% RH), considering the level of agreement with independent measurements made through determining the equilibrated size of aqueous sodium chloride and sodium nitrate droplets.

Coalescence Sampling and Analysis of Aerosols using Aerosol Optical Tweezers.

We present a first exploratory study to assess the use of aerosol optical tweezers as an instrument for sampling and detecting accumulation- and coarse-mode aerosol. A subpicoliter aqueous aerosol droplet is captured in the optical trap and used as a sampling volume, accreting mass from a free-flowing aerosol generated by a medical nebulizer or atomizer. Real-time measurements of the initial stability in size, refractive index, and composition of the sampling droplet inferred from Raman spectroscopy confirm that these quantities can be measured with high accuracy and low noise. Typical standard deviations in size and refractive index of the sampling droplet over a period of 200 s are <±2 nm and <±0.0005, respectively, equivalent to <±0.04% in both measured quantities. A standard deviation of <±1% over a 200 s period is achieved in the spontaneous Raman intensity measurement. When sampling coarse-mode aerosol, mass changes of <10 pg can be detected by the sampling droplet as discrete coalescence events. With accumulation-mode aerosol, we show that fluxes as low as 0.068 pg s-1 can be detected over a 50 s period, equivalent to ∼3 pg of sampled material.

The influence of the surface composition of mixed monolayer films on the evaporation coefficient of water.

We explore the dependence of the evaporation coefficient of water from aqueous droplets on the composition of a surface film, considering in particular the influence of monolayer mixed component films on the evaporative mass flux. Measurements with binary component films formed from long chain alcohols, specifically tridecanol (C13H27OH) and pentadecanol (C15H31OH), and tetradecanol (C14H29OH) and hexadecanol (C16H33OH), show that the evaporation coefficient is dependent on the mole fractions of the two components forming the monolayer film. Immediately at the point of film formation and commensurate reduction in droplet evaporation rate, the evaporation coefficient is equal to a mole fraction weighted average of the evaporation coefficients through the equivalent single component films. As a droplet continues to diminish in surface area with continued loss of water, the more-soluble, shorter alkyl chain component preferentially partitions into the droplet bulk with the evaporation coefficient tending towards that through a single component film formed simply from the less-soluble, longer chain alcohol. We also show that the addition of a long chain alcohol to an aqueous-sucrose droplet can facilitate control over the degree of dehydration achieved during evaporation. After undergoing rapid gas-phase diffusion limited water evaporation, binary aqueous-sucrose droplets show a continued slow evaporative flux that is limited by slow diffusional mass transport within the particle bulk due to the rapidly increasing particle viscosity and strong concentration gradients that are established. The addition of a long chain alcohol to the droplet is shown to slow the initial rate of water loss, leading to a droplet composition that remains more homogeneous for a longer period of time. When the sucrose concentration has achieved a sufficiently high value, and the diffusion constant of water has decreased accordingly so that bulk phase diffusion arrest occurs in the monolayer coated particle, the droplet is found to have lost a greater proportion of its initial water content. A greater degree of slowing in the evaporative flux can be achieved by increasing the chain length of the surface active alcohol, leading to a greater degree of dehydration.

Accurate Measurements of Aerosol Hygroscopic Growth over a Wide Range in Relative Humidity.

Using a comparative evaporation kinetics approach, we describe a new and accurate method for determining the equilibrium hygroscopic growth of aerosol droplets. The time-evolving size of an aqueous droplet, as it evaporates to a steady size and composition that is in equilibrium with the gas phase relative humidity, is used to determine the time-dependent mass flux of water, yielding information on the vapor pressure of water above the droplet surface at every instant in time. Accurate characterization of the gas phase relative humidity is provided from a control measurement of the evaporation profile of a droplet of know equilibrium properties, either a pure water droplet or a sodium chloride droplet. In combination, and by comparison with simulations that account for both the heat and mass transport governing the droplet evaporation kinetics, these measurements allow accurate retrieval of the equilibrium properties of the solution droplet (i.e., the variations with water activity in the mass fraction of solute, diameter growth factor, osmotic coefficient or number of water molecules per solute molecule). Hygroscopicity measurements can be made over a wide range in water activity (from >0.99 to, in principle, <0.05) on time scales of <10 s for droplets containing involatile or volatile solutes. The approach is benchmarked for binary and ternary inorganic solution aerosols with typical uncertainties in water activity of <±0.2% at water activities >0.9 and ∼±1% below 80% RH, and maximum uncertainties in diameter growth factor of ±0.7%. For all of the inorganic systems examined, the time-dependent data are consistent with large values of the mass accommodation (or evaporation) coefficient (>0.1).

Comparison of Methods for Predicting the Compositional Dependence of the Density and Refractive Index of Organic-Aqueous Aerosols.

Representing the physicochemical properties of aerosol particles of complex composition is of crucial importance for understanding and predicting aerosol thermodynamic, kinetic, and optical properties and processes and for interpreting and comparing analysis methods. Here, we consider the representations of the density and refractive index of aqueous-organic aerosol with a particular focus on the dependence of these properties on relative humidity and water content, including an examination of the properties of solution aerosol droplets existing at supersaturated solute concentrations. Using bulk phase measurements of density and refractive index for typical organic aerosol components, we provide robust approaches for the estimation of these properties for aerosol at any intermediate composition between pure water and pure solute. Approximately 70 compounds are considered, including mono-, di- and tricarboxylic acids, alcohols, diols, nitriles, sulfoxides, amides, ethers, sugars, amino acids, aminium sulfates, and polyols. We conclude that the molar refraction mixing rule should be used to predict the refractive index of the solution using a density treatment that assumes ideal mixing or, preferably, a polynomial dependence on the square root of the mass fraction of solute, depending on the solubility limit of the organic component. Although the uncertainties in the density and refractive index predictions depend on the range of subsaturated compositional data available for each compound, typical errors for estimating the solution density and refractive index are less than ±0.1% and ±0.05%, respectively. Owing to the direct connection between molar refraction and the molecular polarizability, along with the availability of group contribution models for predicting molecular polarizability for organic species, our rigorous testing of the molar refraction mixing rule provides a route to predicting refractive indices for aqueous solutions containing organic molecules of arbitrary structure.

Influence of organic films on the evaporation and condensation of water in aerosol.

Uncertainties in quantifying the kinetics of evaporation and condensation of water from atmospheric aerosol are a significant contributor to the uncertainty in predicting cloud droplet number and the indirect effect of aerosols on climate. The influence of aerosol particle surface composition, particularly the impact of surface active organic films, on the condensation and evaporation coefficients remains ambiguous. Here, we report measurements of the influence of organic films on the evaporation and condensation of water from aerosol particles. Significant reductions in the evaporation coefficient are shown to result when condensed films are formed by monolayers of long-chain alcohols [C(n)H(2n+1)OH], with the value decreasing from 2.4 × 10(-3) to 1.7 × 10(-5) as n increases from 12 to 17. Temperature-dependent measurements confirm that a condensed film of long-range order must be formed to suppress the evaporation coefficient below 0.05. The condensation of water on a droplet coated in a condensed film is shown to be fast, with strong coherence of the long-chain alcohol molecules leading to islanding as the water droplet grows, opening up broad areas of uncoated surface on which water can condense rapidly. We conclude that multicomponent composition of organic films on the surface of atmospheric aerosol particles is likely to preclude the formation of condensed films and that the kinetics of water condensation during the activation of aerosol to form cloud droplets is likely to remain rapid.

Variabilities and uncertainties in characterising water transport kinetics in glassy and ultraviscous aerosol.

We present a comprehensive evaluation of the variabilities and uncertainties present in determining the kinetics of water transport in ultraviscous aerosol droplets, alongside new measurements of the water transport timescale in sucrose aerosol. Measurements are performed on individual droplets captured using aerosol optical tweezers and the change in particle size during water evaporation or condensation is inferred from shifts in the wavelength of the whispering gallery mode peaks at which spontaneous Raman scattering is enhanced. The characteristic relaxation timescale (τ) for condensation or evaporation of water from viscous droplets following a change in gas phase relative humidity can be described by the Kohlrausch-Williams-Watts function. To adequately characterise the water transport kinetics and determine τ, sufficient time must be allowed for the particle to progress towards the final state. However, instabilities in the environmental conditions can prevent an accurate characterisation of the kinetics over such long time frames. Comparison with established thermodynamic and diffusional water transport models suggests the determination of τ is insensitive to the choice of thermodynamic treatment. We report excellent agreement between experimental and simulated evaporation timescales, and investigate the scaling of τ with droplet radius. A clear increase in τ is observed for condensation with increase in drying (wait) time. This trend is qualitatively supported by model simulations.

Probing the Evaporation Dynamics of Ethanol/Gasoline Biofuel Blends Using Single Droplet Manipulation Techniques.

Using blends of bioethanol and gasoline as automotive fuel leads to a net decrease in the production of harmful emission compared to the use of pure fossil fuel. However, fuel droplet evaporation dynamics change depending on the mixing ratio. Here we use single particle manipulation techniques to study the evaporation dynamics of ethanol/gasoline blend microdroplets. The use of an electrodynamic balance enables measurements of the evaporation of individual droplets in a controlled environment, while optical tweezers facilitate studies of the behavior of droplets inside a spray. Hence, the combination of both methods is perfectly suited to obtain a complete picture of the evaporation process. The influence of adding varied amounts of ethanol to gasoline is investigated, and we observe that droplets with a greater fraction of ethanol take longer to evaporate. Furthermore, we find that our methods are sensitive enough to observe the presence of trace amounts of water in the droplets. A theoretical model, predicting the evaporation of ethanol and gasoline droplets in dry nitrogen gas, is used to explain the experimental results. Also a theoretical estimation of the saturation of the environment, with other aerosols, in the tweezers is carried out.

Timescales of water transport in viscous aerosol: measurements on sub-micron particles and dependence on conditioning history.

Evaporation studies of single aqueous sucrose aerosol particles as a function of relative humidity (RH) are presented for coarse and fine mode particles down into the submicron size range (600 nm < r < 3.0 µm). These sucrose particles serve as a proxy for biogenic secondary organic aerosols that have been shown to exist, under ambient conditions, in an ultraviscous glassy state, which can affect the kinetics of water mass transport within the bulk phase and hinder particle response to changes in the gas phase water content. A counter-propagating Bessel beams (CPBBs) optical trapping setup is employed to monitor the real-time change in the particle radius with RH decreasing from 75% to 5%. The slow-down of the size change upon each RH step and the deviation from the theoretical equilibrium hygroscopic growth curve indicate the onset of glassy behavior in the RH range of 10-40%. Size-dependent effects were not observed within the uncertainty of the measurements. The influence of the drying time below the glass transition RH on the timescale of subsequent water condensation and re-equilibration for sucrose particles is explored by optical tweezers measurements of micron-sized particles (3 µm < r < 6 µm). The timescale for water condensation and re-equilibration is shown to increase with increasing drying time, i.e. the time over which a viscous particle is dried below 5% RH. These studies demonstrate the importance of the history of the particle conditioning on subsequent water condensation and re-equilibration dynamics of ultraviscous and glassy aerosol particles.

Volatility and oxidative aging of aqueous maleic acid aerosol droplets and the dependence on relative humidity.

The microphysical structure and heterogeneous oxidation by ozone of single aerosol particles containing maleic acid (MA) has been studied using aerosol optical tweezers and cavity enhanced Raman spectroscopy. The evaporation rate of MA from aqueous droplets has been measured over a range of relative humidities and the pure component vapor pressure determined to be (1.7 ± 0.2) × 10(-3) Pa. Variation in the refractive index (RI) of an aqueous MA droplet with relative humidity (RH) allowed the subcooled liquid RI of MA to be estimated as 1.481 ± 0.001. Measurements of the hygroscopic growth are shown to be consistent with equilibrium model predictions from previous studies. Simultaneous measurements of the droplet composition, size, and refractive index have been made during ozonolysis at RHs in the range 50-80%, providing insight into the volatility of organic products, changes in the droplet hygroscopicity, and optical properties. Exposure of the aqueous droplets to ozone leads to the formation of products with a wide range of volatilities spanning from involatile to volatile. Reactive uptake coefficients show a weak dependence on ozone concentration, but no dependence on RH or salt concentration. The time evolving RI depends significantly on the RH at which the oxidation proceeds and can even show opposing trends; while the RI increases with ozone exposure at low relative humidity, the RI decreases when the oxidation proceeds at high relative humidity. The variations in RI are broadly consistent with a framework for predicting RIs for organic components published by Cappa et al. ( J. Geophys. Res. 2011 , 116 , D15204 ). Once oxidized, particles are shown to form amorphous phases on drying rather than crystallization, with slow evaporation kinetics of residual water.

Mass accommodation of water: bridging the gap between molecular dynamics simulations and kinetic condensation models.

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268-300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [Atmos. Chem. Phys. 2012, 12, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99.