Molecular Simulations of the Transport of Molecules across the Liquid/Vapor Interface of Water

Unexpected concentration dependence of the mass accommodation coefficient of water on aqueous triethylene glycol droplets

The mass accommodation coefficient αM of water on aqueous triethylene glycol droplets was determined for water mole fractions in the range xmol = 0.1-0.93 and temperatures between 21 and 26 °C from modulated Mie scattering measurement on single optically-trapped droplets in combination with a kinetic multilayer model. αM reaches minimum values around 0.005 at a critical water concentration of xmol = 0.38, and increases with decreasing water content to a value of ≈0.1 for almost pure triethylene glycol droplets, essentially independent of the temperature. Above xmol = 0.38, αM first increases with increasing water content and then stabilises at a value of ≈0.1 at the lowest temperatures, while at the highest temperature its value remains around 0.005. We analysed the unexpected concentration and temperature dependence with a previously proposed two-step model for mass accommodation which provides concentration and temperature-dependent activation enthalpies and entropies. We suggest that the unexpected minimum in αM at intermediate water concentrations might arise from a more or less saturated hydrogen-bond network that forms at the droplet surface.

Molecular Insights into Chemical Reactions at Aqueous Aerosol Interfaces

Atmospheric aerosols facilitate reactions between ambient gases and dissolved species. Here, we review our efforts to interrogate the uptake of these gases and the mechanisms of their reactions both theoretically and experimentally. We highlight the fascinating behavior of N2O5 in solutions ranging from pure water to complex mixtures, chosen because its aerosol-mediated reactions significantly impact global ozone, hydroxyl, and methane concentrations. As a hydrophobic, weakly soluble, and highly reactive species, N2O5 is a sensitive probe of the chemical and physical properties of aerosol interfaces. We employ contemporary theory to disentangle the fate of N2O5 as it approaches pure and salty water, starting with adsorption and ending with hydrolysis to HNO3, chlorination to ClNO2, or evaporation. Flow reactor and gas-liquid scattering experiments probe even greater complexity as added ions, organic molecules, and surfactants alter the interfacial composition and reaction rates. Together, we reveal a new perspective on multiphase chemistry in the atmosphere.

On the Statistical Mechanics of Mass Accommodation at Liquid-Vapor Interfaces

We propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

The Role of Sum-Frequency Generation Spectroscopy in Understanding On-Surface Reactions and Dynamics in Atmospheric Model-Systems

Surfaces, both water/air and solid/water, play an important role in mediating a multitude of processes central to atmospheric chemistry, particularly in the aerosol phase. However, the study of both static and dynamic properties of surfaces is highly challenging from an experimental standpoint, leading to a lack of molecular level information about the processes that take place at these systems and how they differ from bulk. One of the few techniques that has been able to capture ultrafast surface phenomena is time-resolved sum-frequency generation (SFG) spectroscopy. Since it is both surface-specific and chemically sensitive, the extension of this spectroscopic technique to the time domain makes it possible to study dynamic processes on the femtosecond time scale. In this Perspective, we will explore recent advances made in the field both in terms of studying energy dissipation as well as chemical reactions and the role the surface geometry plays in these processes.

Early Steps in the O2 Scavenger Process in the Aqueous Phase: Hydrazine vs DEHA

We investigate the first direct proton abstraction reactions from reducing agents (RAHs) hydrazine and diethyl hydroxylamine (DEHA), toward dioxygen (O2) in the aqueous phase, spanning ambient to high-temperature conditions. Quantum chemistry methods and molecular dynamics simulations are employed in this study. Quantum chemistry methods are used to analyze the quasi-equilibrium between a reactive conformation and a transition state in the [RAH,O2] cluster. On the other hand, molecular dynamics simulations estimate the probability of observing a reactive conformation of the [RAH,O2] cluster in the solution. In this study, we assume that the energy barrier of the quasi-equilibrium is sufficiently high for the RAH/O2 association process to be at equilibrium. Our findings indicate that the first proton abstraction process from a reactive conformation cluster by DEHA is energetically favored compared to hydrazine. Conversely, the association process of hydrazine and O2 in solution is more favorable than that of DEHA. Consequently, the rate constant for the first proton abstraction process is similar for both hydrazine and DEHA, particularly at high temperatures, with activation energies of approximately 21.5 ± 1.5 kcal mol-1 for both compounds. These results align with recent experiments investigating the complete O2 scavenger process in liquid water with hydrazine and DEHA. Therefore, our findings support the assumption that first proton abstraction reactions are the rate-determining steps in O2 scavenger processes in the aqueous phase.

Phase State, Surface Tension, Water Activity, and Accommodation Coefficient of Water-Organic Clusters Near the Critical Size for Atmospheric New Particle Formation

Interactions between water and organic molecules in sub-4 nm clusters play a significant role in the formation and growth of secondary organic aerosol (SOA) particles. However, a complete understanding of the relevant water microphysics has not yet been achieved due to challenges in the experimental characterization of soft nuclei. Here, we use molecular dynamics simulations to study the phase-mixing states, surface tension, water activity, and water accommodation coefficient of organic-water clusters representative of freshly nucleated SOA particles. Our results reveal large deviations from the behavior expected based on continuum theories. In particular, the phase-mixing state has a strong dependence on cluster size; surface tension displays a minimum at a specific organic-water mass ratio (morg/mw ∼ 4.5 in this study) corresponding to a minimum inhibition of droplet nucleation associated with the Kelvin effect; and the water accommodation coefficient increases by a factor of 2 with nanocluster hygroscopic growth, in agreement with recent experimental studies. Overall, our results yield parametric relations for water microphysical properties in sub-4 nm clusters and provide insight into the role of water in the initial stages of SOA nucleation and growth.

Transport of Heptane Molecules across Water-Vapor Interfaces Laden with Surfactants

Molecular transport across liquid-vapor interfaces covered by surfactant monolayers plays a key role in applications such as fire suppression by foams. The molecular understanding of such transport, however, remains incomplete. This work uses molecular dynamics simulations to investigate the heptane transport across water-vapor interfaces populated with sodium dodecyl sulfate (SDS) surfactants. Heptane molecules' potential of mean force (PMF) and local diffusivity profiles across SDS monolayers with different SDS densities are calculated to obtain heptane's transport resistance. We show that a heptane molecule experiences a finite resistance as it crosses water-vapor interfaces covered by SDS. Such interfacial transport resistance is contributed significantly by heptane molecules' high PMF in the SDS headgroup region and their slow diffusion there. This resistance increases linearly as the SDS density rises from zero but jumps as the density approaches saturation when its value is equivalent to that afforded by a 5 nm thick layer of bulk water. These results are understood by analyzing the micro-environment experienced by a heptane molecule crossing SDS monolayers and the local perturbation it brings to the monolayers. The implications of these findings for the design of surfactants to suppress heptane transport through water-vapor interfaces are discussed.

Mass Transfer through Vapor-Liquid Interfaces Studied by Non-Stationary Molecular Dynamics Simulations

Molecular dynamics (MD) simulations are highly attractive for studying the influence of interfacial effects, such as the enrichment of components, on the mass transfer through the interface. In a recent work, we have presented a steady-state MD simulation method for investigating this phenomenon and tested it using model mixtures with and without interfacial enrichment. The present study extends this work by introducing a non-stationary MD simulation method. A rectangular simulation box that contains a mixture of two components 1 + 2 with a vapor phase in the middle and two liquid phases on both sides is used. Starting from a vapor-liquid equilibrium state, a non-stationary molar flux of component 2 is induced by inserting particles of component 2 into the center of the vapor phase in a pulse-like manner. During the isothermal relaxation process, particles of component 2 pass through the vapor phase, cross the vapor-liquid interface, and enter the liquid phase. The system thereby relaxes into a new vapor-liquid equilibrium state. During the relaxation process, spatially resolved responses for the component densities, fluxes, and pressure are sampled. To reduce the noise and provide measures for the uncertainty of the observables, a set of replicas of simulations is carried out. The new simulation method was applied to study mass transfer in two binary Lennard-Jones mixtures: one that exhibits a strong enrichment of the low-boiling component 2 at the vapor-liquid interface and one that shows no enrichment. Even though both mixtures have similar transport coefficients in the bulk phases, the results for mass transfer differ significantly, indicating that the interfacial enrichment influences the mass transfer.

Solvation, Surface Propensity, and Chemical Reactions of Solutes at Atmospheric Liquid-Vapor Interfaces

surface is covered by oceans, a large number of liquid aerosol particles fill the air, and clouds hold a tiny but critical fraction of Earth's water in the air to influence our climate and hydrology, enabling the lives of humans and ecosystems. The surfaces of these liquids provide the interface for the transfer of gases, for nucleation processes, and for catalyzing important chemical reactions. Coupling a range of spectroscopic tools to liquid microjets has become an important approach to better understanding dynamics, structure, and chemistry at liquid interfaces. Liquid microjets offer stability in vacuum and ambient pressure environments, thus also allowing X-ray photoelectron spectroscopy (XPS) with manageable efforts in terms of differential pumping. Liquid microjets are operated at speeds sufficient to allow for a locally equilibrated surface in terms of water dynamics and solute surface partitioning. XPS is based on the emission of core-level electrons, the binding energy of which is selective for the element and its chemical environment. Inelastic scattering of electrons establishes the probing depth of XPS in the nanometer range and thus its surface sensitivity.In this Account, we focus on aqueous solutions relevant to the surface of oceans, aqueous aerosols, or cloudwater. We are interested in understanding solvation and acid dissociation at the interface, interfacial aspects of reactions with gas-phase reactants, and the interplay of ions with organic molecules at the interface. The strategy is to obtain a link between the molecular-level picture and macroscopic properties and reactivity in the atmospheric context.We show consistency between surface tension and XPS for a range of surface-active organic species as an important proof for interrogating an equilibrated liquid surface. Measurements with organic acids and amines offer important insight into the question of apparent acidity or basicity at the interface. Liquid microjet XPS has settled the debate of the surface enhancement of halide ions, shown using the example of bromide and its oxidation products. Despite the absence of a strong enhancement for the bromide ion, its rate of oxidation by ozone is surface catalyzed through the stabilization of the bromide ozonide intermediate at the interface. In another reaction system, the one between Fe²⁺ and H₂O₂, a similar intermediate in the form of highly valent iron species could not be detected by XPS under the experimental conditions employed, shedding light on the abundance of this intermediate in the environment but also on the constraints within which surface species can be detected. Emphasizing the importance of electrostatic effects, we show how a cationic surfactant attracts charged bromide anions to the interface, accompanied by enhanced oxidation rates by ozone, overriding the role of surfactants as a barrier for the access of gas-phase reactants. The reactivity and structure at interfaces thus result from a subtle balance between hygroscopic and hydrophobic interactions, electrostatic effects, and the structural properties of both liquids and solutes.

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

Infrared Spectroscopy of Stepwise Hydration Motifs of Sulfur Dioxide

Experimental characterization of microscopic events and behaviors of SO₂-H₂O interactions is crucial to understanding SO₂ atmospheric chemistry but has been proven to be very challenging due to the difficulty in size selection. Here, size-dependent development of SO₂ hydrate structure and cluster growth in the SO₂(H₂O)_n (n = 1-16) complexes was probed by infrared spectroscopy based on threshold photoionization using a tunable vacuum ultraviolet free electron laser. Spectral changes with cluster size demonstrate that the sandwich structure initially formed at n = 1 develops into cycle structures with the sulfur and oxygen atoms in a two-dimensional plane (n = 2 and 3) and then into three-dimensional cage structures (n ≥ 4). SO₂ is favorably bound to the surface of larger water clusters. These stepwise features of SO₂ hydration on various sized water clusters contribute to understanding the reactive sites and electrophilicity of SO₂ on cloud droplets, which may have important atmospheric implications for studying the SO₂-containing aerosol systems.

Functionalization and exfoliation of graphite with low temperature pulse plasma in distilled water

Graphene materials exhibit extraordinary properties, but are difficult to produce. The present work describes the possibility of using a plasma process to exfoliate and functionalize graphite flakes. An impulse plasma phase is generated at a liquid surface to produce chemical species and shock waves in order to modify the reactive liquid as well as the graphite flakes. With this process, industrial graphite was treated. 20% thickness diminution was observed, and the formation of a random turbostratic structure. The exfoliation occurs with small amount of functionalization of the surface. Even after treatment, the graphite flakes present a low defect density compared with other treated graphite obtained by more conventional chemical treatments. This process is a new way to exfoliate graphite and to produce functionalized graphenic materials.

Experimental evidence that halogen bonding catalyzes the heterogeneous chlorination of alkenes in submicron liquid droplets

A key challenge in predicting the multiphase chemistry of aerosols and droplets is connecting reaction probabilities, observed in an experiment, with the kinetics of individual elementary steps that control the chemistry that occurs across a gas/liquid interface. Here we report evidence that oxygenated molecules accelerate the heterogeneous reaction rate of chlorine gas with an alkene (squalene, Sqe) in submicron droplets. The effective reaction probability for Sqe is sensitive to both the aerosol composition and gas phase environment. In binary aerosol mixtures with 2-decyl-1-tetradecanol, linoleic acid and oleic acid, Sqe reacts 12-23× more rapidly than in a pure aerosol. In contrast, the reactivity of Sqe is diminished by 3× when mixed with an alkane. Additionally, small oxygenated molecules in the gas phase (water, ethanol, acetone, and acetic acid) accelerate (up to 10×) the heterogeneous chlorination rate of Sqe. The overall reaction mechanism is not altered by the presence of these aerosol and gas phase additives, suggesting instead that they act as catalysts. Since the largest rate acceleration occurs in the presence of oxygenated molecules, we conclude that halogen bonding enhances reactivity by slowing the desorption kinetics of Cl2 at the interface, in a way that is analogous to decreasing temperature. These results highlight the importance of relatively weak interactions in controlling the speed of multiphase reactions important for atmospheric and indoor environments.

Assessing the Impact of Solvent Selection on Vibrational Sum-Frequency Scattering Spectroscopy Experiments

The development of vibrational sum-frequency scattering (S-VSF) spectroscopy has opened the door to directly probing nanoparticle surfaces with an interfacial and chemical specificity that was previously reserved for planar interfacial systems. Despite its potential, challenges remain in the application of S-VSF spectroscopy beyond simplified chemical systems. One such challenge includes infrared absorption by an absorptive continuous phase, which will alter the spectral lineshapes within S-VSF spectra. In this study, we investigate how solvent vibrational modes manifest in S-VSF spectra of surfactant stabilized nanoemulsions and demonstrate how corrections for infrared absorption can recover the spectral features of interfacial solvent molecules. We also investigate infrared absorption for systems with the absorptive phase dispersed in a nonabsorptive continuous phase to show that infrared absorption, while reduced, will still impact the S-VSF spectra. These studies are then used to provide practical recommendations for anyone wishing to use S-VSF to study nanoparticle surfaces where absorptive solvents are present.

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.

Transport of Gaseous Hydrogen Peroxide and Ozone into Bulk Water vs. Electrosprayed Aerosol

Production and transport of reactive species through plasma–liquid interactions play a significant role in multiple applications in biomedicine, environment, and agriculture. Experimental investigations of the transport mechanisms of typical air plasma species: hydrogen peroxide (H2O2) and ozone (O3) into water are presented. Solvation of gaseous H2O2 and O3 from an airflow into water bulk vs. electrosprayed microdroplets was measured, while changing the water flow rate and applied voltage, during different treatment times and gas flow rates. The solvation rate of H2O2 and O3 increased with the treatment time and the gas–liquid interface area. The total surface area of the electrosprayed microdroplets was larger than that of the bulk, but their lifetime was much shorter. We estimated that only microdroplets with diameters below ~40 µm could achieve the saturation by O3 during their lifetime, while the saturation by H2O2 was unreachable due to its depletion from air. In addition to the short-lived flying microdroplets, the longer-lived bottom microdroplets substantially contributed to H2O2 and O3 solvation in water electrospray. This study contributes to a better understanding of the gaseous H2O2 and O3 transport into water and will lead to design optimization of the water spray and plasma-liquid interaction systems.

Understanding the properties of methanesulfinic acid at the air-water interface

Methanesulfinic acid (MSIA), an organic sulfur compound, is mainly produced in the oxidation process of dimethyl sulfide in the atmosphere. The properties of MSIA at the air-water interface were studied using molecular dynamics (MD) simulations. The result shows that the lowest system free energy is located at the interface. Because the free energy difference between the interface and water phase is 3.2 kJ mol^-1, the MSIA molecule can easily get out of the free energy well and travel to water phase by the thermal motion, leading to only a 21% probability of its occurrence at the interface. The MSIA molecule tends to tilt at the interface with the sulfino group (-S(O)-OH) pointing toward the water phase. The feature of hydration status at the air-water interface may be favorable to the heterogeneous oxidation of MSIA.

Mass Transfer Thermodynamics through a Gas-Liquid Interface

Molecular level information about thermodynamic variations (enthalpy, entropy, and free energy) of a gas molecule as it crosses a gas-liquid interface is strongly lacking from an experimental perspective under equilibrium conditions. Herein, we perform in situ measurements of water interacting with the ionic liquid (IL) 1-butyl-3-methylimidazolium acetate, [C₄mim][Ace], using ambient pressure X-ray photoelectron spectroscopy in order to assess the interfacial uptake of water quantitatively as a function of temperature, pressure, and water mole fraction ( x_w). The surface spectroscopy results are compared to existing bulk water absorption experiments, showing that the amount of water in the interfacial region is consistently greater than that in the bulk. The enthalpy and entropy of water sorption vary significantly between the gas-liquid interface and the bulk as a function of x_w, with a crossover that occurs near x_w = 0.6 where the water-IL mixture converts from being homogeneous ( x_w < 0.6) to nanostructured ( x_w > 0.6). Free energy results reveal that water at the gas-IL interface is thermodynamically more favorable than that in the bulk, consistent with the enhanced water concentration in the interfacial region. The results herein show that the efficacy for an ionic liquid to absorb a gas phase molecule is not merely a function of bulk solvation parameters but also is significantly influenced by the thermodynamics occurring across the gas-IL interface during the mass transfer process.

An Overview of Dynamic Heterogeneous Oxidations in the Troposphere

Due to the adverse effect of atmospheric aerosols on public health and their ability to affect climate, extensive research has been undertaken in recent decades to understand their sources and sinks, as well as to study their physical and chemical properties. Atmospheric aerosols are important players in the Earth’s radiative budget, affecting incoming and outgoing solar radiation through absorption and scattering by direct and indirect means. While the cooling properties of pure inorganic aerosols are relatively well understood, the impact of organic aerosols on the radiative budget is unclear. Additionally, organic aerosols are transformed through chemical reactions during atmospheric transport. The resulting complex mixture of organic aerosol has variable physical and chemical properties that contribute further to the uncertainty of these species modifying the radiative budget. Correlations between oxidative processing and increased absorptivity, hygroscopicity, and cloud condensation nuclei activity have been observed, but the mechanisms behind these phenomena have remained unexplored. Herein, we review environmentally relevant heterogeneous mechanisms occurring on interfaces that contribute to the processing of aerosols. Recent laboratory studies exploring processes at the aerosol–air interface are highlighted as capable of generating the complexity observed in the environment. Furthermore, a variety of laboratory methods developed specifically to study these processes under environmentally relevant conditions are introduced. Remarkably, the heterogeneous mechanisms presented might neither be feasible in the gas phase nor in the bulk particle phase of aerosols at the fast rates enabled on interfaces. In conclusion, these surface mechanisms are important to better understand how organic aerosols are transformed in the atmosphere affecting the environment.

Chain-propagation, chain-transfer, and hydride-abstraction by cyclic carbocations on water surfaces

New mechanisms for the growth and increase in complexity of atmospheric aerosol particles are elucidated. The present findings will also be useful for interfacial polymer/oligomer synthesis.

A comparison of sodium and hydrogen halides at the air-water interface

New molecular models, parameterized to ab initio calculations, were developed to describe HBr and HI at the air-water interface. These were used to compare how the air-water interface influenced dissociation of NaX and HX, with X being Cl, Br, or I, and also their propensity for the interface. The polarizable multistate empirical valence bond method, which explicitly describes proton sharing, was used to model HX. Results showed that the air-water interface suppressed HX dissociation from a contact ion pair to a solvent separated to a greater degree than NaX dissociation. Furthermore, HX had a greater propensity for the interface than NaX, which was a consequence of the hydronium ion having a greatest interfacial activity of all species studied. As a consequence of this, the average configuration of dissociated HX, while in both contact ion and solvent separated ion pairs near the air-water interface, is with the dissociated hydrogen oriented more towards the air than the X atom.

Solvation Dynamics of CO₂(g) by Monoethanolamine at the Gas-Liquid Interface: A Molecular Mechanics Approach

A classical force field approach was used to characterize the solvation dynamics of high-density CO₂(g) by monoethanolamine (MEA) at the air–liquid interface. Intra- and intermolecular CO₂ and MEA potentials were parameterized according to the energetics calculated at the MP2 and BLYP-D2 levels of theory. The thermodynamic properties of CO₂ and MEA, such as heat capacity and melting point, were consistently predicted using this classical potential. An approximate interfacial simulation for CO₂(g)/MEA(l) was performed to monitor the depletion of the CO₂(g) phase, which was influenced by amino and hydroxyl groups of MEA. There are more intramolecular hydrogen bond interactions notably identified in the interfacial simulation than the case of bulk MEA(l) simulation. The hydroxyl group of MEA was found to more actively approach CO₂ and overpower the amino group to interact with CO₂ at the air–liquid interface. With artificially reducing the dipole moment of the hydroxyl group, CO₂–amino group interaction was enhanced and suppressed CO₂(g) depletion. The hydroxyl group of MEA was concluded to play dual but contradictory roles for CO₂ capture.

A polarizable continuum model for molecules at spherical diffuse interfaces

We present an extension of the Polarizable Continuum Model (PCM) to simulate solvent effects at diffuse interfaces with spherical symmetry, such as nanodroplets and micelles. We derive the form of the Green's function for a spatially varying dielectric permittivity with spherical symmetry and exploit the integral equation formalism of the PCM for general dielectric environments to recast the solvation problem into a continuum solvation framework. This allows the investigation of the solvation of ions and molecules in nonuniform dielectric environments, such as liquid droplets, micelles or membranes, while maintaining the computationally appealing characteristics of continuum solvation models. We describe in detail our implementation, both for the calculation of the Green's function and for its subsequent use in the PCM electrostatic problem. The model is then applied on a few test systems, mainly to analyze the effect of interface curvature on solvation energetics.

Super-Maxwellian helium evaporation from pure and salty water

Helium atoms evaporate from pure water and salty solutions in super-Maxwellian speed distributions, as observed experimentally and modeled theoretically. The experiments are performed by monitoring the velocities of dissolved He atoms that evaporate from microjets of pure water at 252 K and 4-8.5 molal LiCl and LiBr at 232-252 K. The average He atom energies exceed the flux-weighted Maxwell-Boltzmann average of 2RT by 30% for pure water and 70% for 8.5m LiBr. Classical molecular dynamics simulations closely reproduce the observed speed distributions and provide microscopic insight into the forces that eject the He atoms from solution. Comparisons of the density profile and He kinetic energies across the water-vacuum interface indicate that the He atoms are accelerated by He-water collisions within the top 1-2 layers of the liquid. We also find that the average He atom kinetic energy scales with the free energy of solvation of this sparingly soluble gas. This free-energy difference reflects the steeply decreasing potential of mean force on the He atoms in the interfacial region, whose gradient is the repulsive force that tends to expel the atoms. The accompanying sharp decrease in water density suppresses the He-water collisions that would otherwise maintain a Maxwell-Boltzmann distribution, allowing the He atom to escape at high energies. Helium is especially affected by this reduction in collisions because its weak interactions make energy transfer inefficient.

Hydrogen-Bonding Structure at Zwitterionic Lipid/Water Interface

Interface structure of water/[3-palmitoyl-2-oleoyl-D-glycero-1-phosphatidylcholine] (POPC) lipid layer is investigated with molecular dynamics (MD) simulation by analyzing the recent heterodyne-detected vibrational sum-frequency generation (HD VSFG) spectroscopy. The MD simulation clearly reproduced the experimental HD VSFG spectrum of imaginary susceptibility (Im[χ]), which exhibits two positive bands in the OH stretching vibrations of water. With the help of decomposition MD analysis, we found three kinds of interfacial water in relation to the HD VSFG spectrum. The low-frequency positive band is attributed to the water pointing toward the lipid side, whose orientation is influenced by negatively charged phosphate and positively charged choline of POPC. The high-frequency positive band is attributed to the water bonding with the carbonyl groups of the lipid. The gap between the two positive bands indicates the interfacial water pointing toward the bulk water phase in the vicinity of the choline groups.

Temperature dependence of local solubility of hydrophobic molecules in the liquid-vapor interface of water

One important aspect of the hydrophobic effect is that solubility of small, nonpolar molecules in liquid water decreases with increasing temperature. We investigate here how the characteristic temperature dependence in liquid water persists or changes in the vicinity of the liquid-vapor interface. From the molecular dynamics simulation and the test-particle insertion method, the local solubility Σ of methane in the liquid-vapor interface of water as well as Σ of nonpolar solutes in the interface of simple liquids are calculated as a function of the distance z from the interface. We then examine the temperature dependence of Σ under two conditions: variation of Σ at fixed position z and that at fixed local solvent density around the solute molecule. It is found that the temperature dependence of Σ at fixed z depends on the position z and the system, whereas Σ at fixed local density decreases with increasing temperature for all the model solutions at any fixed density between vapor and liquid phases. The monotonic decrease of Σ under the fixed-density condition in the liquid-vapor interface is in accord with what we know for the solubility of nonpolar molecules in bulk liquid water under the fixed-volume condition but it is much robust since the solvent density to be fixed can be anything between the coexisting vapor and liquid phases. A unique feature found in the water interface is that there is a minimum in the local solubility profile Σ(z) on the liquid side of the interface. We find that with decreasing temperature the minimum of Σ grows and at the same time the first peak in the oscillatory density profile of water develops. It is likely that the minimum of Σ is due to the layering structure of the free interface of water.

Ballistic Evaporation and Solvation of Helium Atoms at the Surfaces of Protic and Hydrocarbon Liquids

Atomic and molecular solutes evaporate and dissolve by traversing an atomically thin boundary separating liquid and gas. Most solutes spend only short times in this interfacial region, making them difficult to observe. Experiments that monitor the velocities of evaporating species, however, can capture their final interactions with surface solvent molecules. We find that polarizable gases such as N2 and Ar evaporate from protic and hydrocarbon liquids with Maxwell-Boltzmann speed distributions. Surprisingly, the weakly interacting helium atom emerges from these liquids at high kinetic energies, exceeding the expected energy of evaporation from salty water by 70%. This super-Maxwellian evaporation implies in reverse that He atoms preferentially dissolve when they strike the surface at high energies, as if ballistically penetrating into the solvent. The evaporation energies increase with solvent surface tension, suggesting that He atoms require extra kinetic energy to navigate increasingly tortuous paths between surface molecules.

Near-unity mass accommodation coefficient of organic molecules of varying structure

Atmospheric aerosol particles have a significant effect on global climate, air quality, and consequently human health. Condensation of organic vapors is a key process in the growth of nanometer-sized particles to climate relevant sizes. This growth is very sensitive to the mass accommodation coefficient α, a quantity describing the vapor uptake ability of the particles, but knowledge on α of atmospheric organics is lacking. In this work, we have determined α for four organic molecules with diverse structural properties: adipic acid, succinic acid, naphthalene, and nonane. The coefficients are studied using molecular dynamics simulations, complemented with expansion chamber measurements. Our results are consistent with α = 1 (indicating nearly perfect accommodation), regardless of the molecular structural properties, the phase state of the bulk condensed phase, or surface curvature. The results highlight the need for experimental techniques capable of resolving the internal structure of nanoparticles to better constrain the accommodation of atmospheric organics.

Nonequilibrium kinetic boundary condition at the vapor-liquid interface of argon

A boundary condition for the Boltzmann equation (kinetic boundary condition, KBC) at the vapor-liquid interface of argon is constructed with the help of molecular dynamics (MD) simulations. The KBC is examined at a constant liquid temperature of 85 K in a wide range of nonequilibrium states of vapor. The present investigation is an extension of a previous one by Ishiyama, Yano, and Fujikawa [Phys. Rev. Lett. 95, 084504 (2005)] and provides a more complete form of the KBC. The present KBC includes a thermal accommodation coefficient in addition to evaporation and condensation coefficients, and these coefficients are determined in MD simulations uniquely. The thermal accommodation coefficient shows an anisotropic behavior at the interface for molecular velocities normal versus tangential to the interface. It is also found that the evaporation and condensation coefficients are almost constant in a fairly wide range of nonequilibrium states. The thermal accommodation coefficient of the normal velocity component is almost unity, while that of the tangential component shows a decreasing function of the density of vapor incident on the interface, indicating that the tangential velocity distribution of molecules leaving the interface into the vapor phase may deviate from the tangential parts of the Maxwell velocity distribution at the liquid temperature. A mechanism for the deviation of the KBC from the isotropic Maxwell KBC at the liquid temperature is discussed in terms of anisotropic energy relaxation at the interface. The liquid-temperature dependence of the present KBC is also discussed.

Insight into the molecular mechanism of water evaporation via the finite temperature string method

The process of water's evaporation at its liquid/air interface has proven challenging to study experimentally and, because it constitutes a rare event on molecular time scales, presents a challenge for computer simulations as well. In this work, we simulated water's evaporation using the classical extended simple point charge model water model, and identified a minimum free energy path for this process in terms of 10 descriptive order parameters. The measured free energy change was 7.4 kcal/mol at 298 K, in reasonable agreement with the experimental value of 6.3 kcal/mol, and the mean first-passage time was 1375 ns for a single molecule, corresponding to an evaporation coefficient of 0.25. In the observed minimum free energy process, the water molecule diffuses to the surface, and tends to rotate so that its dipole and one O-H bond are oriented outward as it crosses the Gibbs dividing surface. As the water molecule moves further outward through the interfacial region, its local density is higher than the time-averaged density, indicating a local solvation shell that protrudes from the interface. The water molecule loses donor and acceptor hydrogen bonds, and then, with its dipole nearly normal to the interface, stops donating its remaining hydrogen bond. At that point, when the final, accepted hydrogen bond is broken, the water molecule is free. We also analyzed which order parameters are most important in the process and in reactive trajectories, and found that the relative orientation of water molecules near the evaporating molecule, and the number of accepted hydrogen bonds, were important variables in reactive trajectories and in kinetic descriptions of the process.