Dynamics of water molecules at liquid–vapour interfaces of aqueous ionic solutions: effects of ion concentration

Mechanism, Kinetics, and Potential of Mean Force of Evaporation of Water from Aqueous Sodium Chloride Solutions of Varying Concentrations

The mechanism, kinetics, and potential of mean force of evaporation of water from aqueous NaCl solutions are investigated through both unbiased molecular dynamics simulations and also biased simulations using the umbrella sampling method. The results are obtained for aqueous solutions of three different NaCl concentrations ranging from 0.6 to 6.0 m and also for pure water. The rate of evaporation is found to decrease in the presence of ions. It is found that the process of evaporation of a surface water molecule from ionic solutions can be triggered through its collision with another water or chloride ion. Such collisions provide the additional kinetic energy that is required for evaporation. However, when the collision takes place with a Cl^- ion, the evaporation of the escaping water also involves a collision with water in the vicinity of the ion at the same time along with the ion-water collision. These two collisions together provide the required kinetic energy for escape of the evaporating water molecule. Thus, the mechanism of evaporation process of ionic solutions can be more complex than that of pure water. The potential of mean force (PMF) of evaporation is found to be positive and it increases with increasing ion concentration. Also, no barrier in the PMF is found to be present for the condensation of water from vapor phase to the surfaces of the solutions. A detailed analysis of the unsuccessful evaporation attempts by surface water molecules is also made in the current study.

Phase separation property of a hydrophobic deep eutectic solvent-water binary mixture: A molecular dynamics simulation study

Over the past decade, deep eutectic solvents (DESs) have earned applicability in numerous fields as non-flammable, non-volatile, and greener alternatives to conventional organic solvents. In a first of its kind, a hydrophobic DES composed of a 1:1 mixture of oleic acid and lidocaine was recently reported, possessing a lower critical solution temperature in water. The thermoreversible phase property of this DES-water system was utilized to sequester out dye molecules from their aqueous solutions. In this article, we explore the phase separation phenomena for this particular DES in its aqueous solution using an all-atom molecular dynamics simulation. A 50 wt. % solution of the DES in water was studied at three different temperatures (253, 293, and 313 K) to understand the various molecular interactions that dictate the phase segregation property of these systems. In this work, we have elaborated on the importance of hydrogen bonding interactions and the non-bonding interactions between the components and the competition between the two that leads to phase separation. Overall, we observe that the increase in unfavorable interaction between the DES components and water with increasing temperature determines the phase separation behavior. We have also studied the modification in the dynamical properties of water molecules close to the phase boundary. Such molecular insights would be beneficial for designing novel solvent systems that can be used as extraction-based media in industries.

Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

A first principles molecular dynamics study of excess electron and lithium atom solvation in water-ammonia mixed clusters: structural, spectral, and dynamical behaviors of [(H2O)5NH3]- and Li(H2O)5NH3 at finite temperature

First principles molecular dynamics simulations are carried out to investigate the solvation of an excess electron and a lithium atom in mixed water-ammonia cluster (H(2)O)(5)NH(3) at a finite temperature of 150 K. Both [(H(2)O)(5)NH(3)](-) and Li(H(2)O)(5)NH(3) clusters are seen to display substantial hydrogen bond dynamics due to thermal motion leading to many different isomeric structures. Also, the structures of these two clusters are found to be very different from each other and also very different from the corresponding neutral cluster without any excess electron or the metal atom. Spontaneous ionization of Li atom occurs in the case of Li(H(2)O)(5)NH(3). The spatial distribution of the singly occupied molecular orbital shows where and how the excess (or free) electron is primarily localized in these clusters. The populations of single acceptor (A), double acceptor (AA), and free (NIL) type water and ammonia molecules are found to be significantly high. The dangling hydrogens of these type of water or ammonia molecules are found to primarily capture the free electron. It is also found that the free electron binding motifs evolve with time due to thermal fluctuations and the vertical detachment energy of [(H(2)O)(5)NH(3)](-) and vertical ionization energy of Li(H(2)O)(5)NH(3) also change with time along the simulation trajectories. Assignments of the observed peaks in the vibrational power spectra are done and we found a one to one correlation between the time-averaged populations of water and ammonia molecules at different H-bonding sites with the various peaks of power spectra. The frequency-time correlation functions of OH stretch vibrational frequencies of these clusters are also calculated and their decay profiles are analyzed in terms of the dynamics of hydrogen bonded and dangling OH modes. It is found that the hydrogen bond lifetimes in these clusters are almost five to six times longer than that of pure liquid water at room temperature.

Electrostatic properties of aqueous salt solution interfaces: a comparison of polarizable and nonpolarizable ion models

The effects of ion force field polarizability on the interfacial electrostatic properties of approximately 1 M aqueous solutions of NaCl, CsCl, and NaI are investigated using molecular dynamics simulations employing both nonpolarizable and Drude-polarizable ion sets. Differences in computed depth-dependent orientational distributions, "permanent" and induced dipole and quadrupole moment profiles, and interfacial potentials are obtained for both ion sets to further elucidate how ion polarizability affects interfacial electrostatic properties among the various salts relative to pure water. We observe that the orientations and induced dipoles of water molecules are more strongly perturbed in the presence of polarizable ions via a stronger ionic double layer effect arising from greater charge separation. Both anions and cations exhibit enhanced induced dipole moments and strong z alignment in the vicinity of the Gibbs dividing surface (GDS) with the magnitude of the anion induced dipoles being nearly an order of magnitude larger than those of the cations and directed into the vapor phase. Depth-dependent profiles for the trace and z z components of the water molecular quadrupole moment tensors reveal 40% larger quadrupole moments in the bulk phase relative to the vapor which mimics a similar observed 40% increase in the average water dipole moment. Across the GDS, the water molecular quadrupole moments increase nonmonotonically (in contrast to the water dipoles) and exhibit a locally reduced contribution just below the surface due to both orientational and polarization effects. Computed interfacial potentials for the nonpolarizable salts yield values 20-60 mV more positive than pure water and increase by an additional 30-100 mV when ion polarizability is included. A rigorous decomposition of the total interfacial potential into ion monopole, water and ion dipole, and water quadrupole components reveals that a very strong, positive ion monopole contribution is offset by negative contributions from all other potential sources. Water quadrupole components modulated by the water density contribute significantly to the observed interfacial potential increments and almost entirely explain observed differences in the interfacial potentials for the two chloride salts. By lumping all remaining nonquadrupole interfacial potential contributions into a single "effective" dipole potential, we observe that the ratio of quadrupole to "effective" dipole contributions range from 2:1 in CsCl to 1:1.5 in NaI, suggesting that both contributions are comparably important in determining the interfacial potential increments. We also find that oscillations in the quadrupole potential in the double layer region are opposite in sign and partially cancel those of the "effective" dipole potential.

Theory and computer simulation of solute effects on the surface tension of liquids

A complete description of the thermodynamics of planar mixed solute-solvent interfaces suitable for the analysis of computer simulation data is provided. The approach uses surface probability distributions to characterize the interface regions, coupled with radial distribution functions and the Kirkwood-Buff theory of solutions to characterize the bulk solution properties. The approach is then used to understand the relationship between changes in the surface tension, the degree of surface adsorption or depletion, and the bulk solution properties of two aqueous solute systems. The first, aqueous NaCl solutions, provides an example of a surface excluded solute. The second, aqueous methanol solutions, provides an example of a surface adsorbed solute. The numerical results support the theoretical relationships described here and provide a consistent picture of the thermodynamics of solution interfaces involving any number of components which can be applied to a wide variety of systems.

On the nature of ions at the liquid water surface

A qualitatively new understanding of the nature of ions at the liquid water surface is emerging. Traditionally, the characterization of liquid surfaces has been limited to macroscopic experimental techniques such as surface tension and electrostatic potential measurements, wherein the microscopic picture then has to be inferred by applying theoretical models. Because the surface tension of electrolyte solutions generally increases with ion concentration, all inorganic ions have been thought to be repelled from the air-water interface, leaving the outermost surface layer essentially devoid of ions. This oversimplified picture has recently been challenged: first by chemical kinetics measurements, then by theoretical molecular dynamics simulations using polarizable models, and most recently by new surface sensitive experimental observations. Here we present an overview of the nature of the interfacial structure of electrolyte solutions and give a detailed description of the new picture that is emerging.

Monte Carlo simulation of equilibrium reactions at modified vapor-liquid interfaces

The equilibrium conversion of a chemical reaction is known to be affected by its local environment. Various factors may alter reaction equilibria, including shifts in pressure or temperature, solvation, adsorption within porous materials, or the presence of an interface. Previously, reactive Monte Carlo simulations have been used to predict the equilibrium behavior of chemical reactions at vapor-liquid interfaces. Here, a route is tested for tuning the interfacial conversion of a Lennard-Jones dimerization reaction by adding surfactants to the vapor-liquid interface. Several temperatures are explored as well as several different surfactant models. Even with the addition of a small concentration of surfactants, the simulations predict significant shifts in the conversion at the interface. In general, the shifts in the conversion tend to be related to the values of the interfacial tension.

Electrochemistry of TEMPO in the aqueous liquid/vapor interfacial region: measurements of the lateral mobility and kinetics of surface partitioning

A new method is described to simultaneously determine the kinetics of surface partitioning and the lateral diffusion constant of redox active amphiphiles. It concerns water-soluble amphiphiles for which the surface adsorption equilibrium constant and the solution diffusion constant are measured independently. The method involves cyclic voltammetric experiments carried out at the air/water interface with microband electrodes aligned with the plane of the water surface. Typically, 100 nm wide, 1.0 cm long microband electrodes are fabricated by the vacuum vapor deposition of gold films on glass. The front face of the electrode substrates are coated with impermeable, dimensionally stable, polymer barrier films with thickness L in the range of approximately 0.1-1.0 microm. Fracturing such gold-coated glass substrates exposes gold microbands. The recorded voltammetric current sensitively depends on the barrier film thickness, the surfactant surface diffusion constant, Dsurf, and its rate constant of desorption, kdes. For a given surfactant, such as the nitroxyl piperidine free radical TEMPO featured in this report, large currents are observed with microband electrodes that do not carry a barrier film (L = 0). This is because the surfactant surface population diffusing along the air/water interface can be directly electro-oxidized at the edge of the microband. Smaller currents are measured in the presence of a barrier film, since, in those instances, the surface population may contribute to the voltammetric current only via a mechanism involving surfactant desorption from the water surface into bulk, where it contributes to the three-dimensional solution diffusion processes. The quantitative interpretation of the voltammetric experiments was made possible with finite element simulations with FEMLAB. These produce a set of calibration curves, Dsurf versus log kdes, for each value of the barrier film thickness. The intersection of the calibration curves determines the unique values of Dsurf and kdes. For TEMPO, Dsurf = 4.4 +/- 1.2 x 10(-5) cm2/s and kdes >/= 2 x 10(4) s(-1). Surfactant desorption rate constants of this magnitude have not been previously experimentally accessible. Since, in our earlier report (Wu, D. G.; Malec, A. D.; Head-Gordon, M.; Majda, M. J. Am. Chem. Soc. 2005, 27, 4490-4496), we showed that TEMPO is not immersed in water and that it diffuses along the interface hydrogen-bonded to just one or two water molecules, its Dsurf value approximates the water diffusion constant in the aqueous liquid-vapor interfacial region.

Polarization and experimental configuration analyses of sum frequency generation vibrational spectra, structure, and orientational motion of the air/water interface

Here we report a detailed study on spectroscopy, structure, and orientational distribution, as well as orientational motion, of water molecules at the air/water interface, investigated with sum frequency generation vibrational spectroscopy (SFG-VS). Quantitative polarization and experimental configuration analyses of the SFG data in different polarizations with four sets of experimental configurations can shed new light on our present understanding of the air/water interface. Firstly, we concluded that the orientational motion of the interfacial water molecules can only be in a limited angular range, instead of rapidly varying over a broad angular range in the vibrational relaxation time as suggested previously. Secondly, because different vibrational modes of different molecular species at the interface has different symmetry properties, polarization and symmetry analyses of the SFG-VS spectral features can help the assignment of the SFG-VS spectra peaks to different interfacial species. These analyses concluded that the narrow 3693 cm(-1) and broad 3550 cm(-1) peaks belong to C(infinityv) symmetry, while the broad 3250 and 3450 cm(-1) peaks belong to the symmetric stretching modes with C2v symmetry. Thus, the 3693 cm(-1) peak is assigned to the free OH, the 3550 cm(-1) peak is assigned to the singly hydrogen-bonded OH stretching mode, and the 3250 and 3450 cm(-1) peaks are assigned to interfacial water molecules as two hydrogen donors for hydrogen bonding (with C2v symmetry), respectively. Thirdly, analysis of the SFG-VS spectra concluded that the singly hydrogen-bonded water molecules at the air/water interface have their dipole vector directed almost parallel to the interface and is with a very narrow orientational distribution. The doubly hydrogen-bonded donor water molecules have their dipole vector pointing away from the liquid phase.

Analysis of chemical kinetics at the gas-aqueous interface for submicron aerosols

The effect of kinetics of chemical reactions in the gas-liquid interface between atmospheric gases and reactive solute in dilute aqueous aerosols is analysed in order to see if such processes will affect the overall uptake rate. Accordingly, a parameterization of such heterogeneous reactions was derived, taking into account interfacial reactions. Gibbs surface excess concentration of both reactive compounds and stable compounds leads to higher heterogeneous reaction rates in comparison to aqueous phase bulk reactions. An analytical formulation shows that the surface reactions may be of considerable importance for the uptake process in the case of small liquid aerosols even in the absence of organic film on the surface. In particular, we demonstrate that the uptake rate of atmospheric gas-phase oxidants (such as OH, NO(3) or O(3)) reacting with volatile organic compounds (such as ethanol or methanol) is increased by more than 10% for atmospheric aerosols with diameters lower than 0.1 microm. This effect is in addition intensified in the case of reactions of atmospheric oxidants with liquid aerosols containing organic surfactants, such as semi-volatile organic compounds, i.e., the chemical reactions at the gas-liquid interface may be dominant in the main uptake process for atmospheric submicron aerosols.

Monte Carlo Simulation of Equilibrium Reactions at Vapor−Liquid Interfaces

Chemical reactions are known to behave differently, depending upon their local environment. While the interactions with neighboring molecules may alter both the kinetics of chemical reactions and the overall equilibrium conversion, we have performed simulations of the latter. The particular environment that we address is the vapor-liquid interface, since only a few, limited studies have explored the influence of an interface on equilibrium reaction behavior. Simple dimerization reactions are modeled, as well as more complex multicomponent reactions, using the reactive Monte Carlo (RxMC) simulation technique. We find that the conversion of a reaction can be markedly different at an interface as compared to the bulk vapor and liquid phases, and these trends are analyzed with respect to specific intermolecular interactions. In conjunction, we calculate the surface tension of the reacting fluids at the interface, which is found to have unusual scaling behavior, with respect to the system temperature.

Liquid-vapor interfaces of water-acetonitrile mixtures of varying composition

Detailed molecular-dynamics simulations are carried out to investigate the equilibrium and dynamical properties of water-acetonitrile mixtures of varying composition. Altogether, we have simulated eight different systems of different concentrations of acetonitrile. The inhomogeneous density and anisotropic orientational profiles at interfaces, surface tension, and also the distribution of hydrogen bonds are calculated for both water and acetonitrile molecules. The dynamical aspects of the interfaces are investigated in terms of the anisotropic diffusion and dipole orientational relaxation of interfacial water and acetonitrile molecules. For both structural and dynamical properties, the behaviors of the interfaces are compared with those of the corresponding bulk phases. A comparison between the present theoretical results and experimental findings, wherever available, is also made to verify the usefulness of the molecular models employed in the present study for predicting interfacial properties.

Liquid-vapor interfacial properties of water-ammonia mixtures: Dependence on ammonia concentration

The equilibrium and dynamical properties of the liquid-vapor interfaces of water-ammonia mixtures are investigated by means of molecular-dynamics simulations. Altogether, we have simulated seven different systems of different concentration of ammonia. The inhomogeneous density, anisotropic orientational profiles, surface tension, and the pattern of hydrogen bonding are calculated for both water and ammonia molecules in order to characterize the location, width, thermodynamic aspects, and microscopic structure of the liquid-vapor interfaces of each of the water-ammonia systems. The dynamical aspects of the interfaces are investigated in terms of the anisotropic diffusion and dipole orientational relaxation of water and ammonia molecules. The properties of the interfaces are compared with those of the corresponding bulk phases. The present theoretical results are also compared with experimental findings wherever available.

Molecular Dynamics Study of the Liquid−Vapor Interface of Acetonitrile: Equilibrium and Dynamical Properties

The equilibrium and dynamical properties of the liquid-vapor interface of pure acetonitrile are studied by means of molecular dynamics simulations. Both nonpolarizable and polarizable models are employed in the present study. For the nonpolarizable model, the simulations are carried out for two different system sizes and at two different temperatures whereas the simulation with the polarizable model is done for a single system. The inhomogeneous density, anisotropic orientational profile, the width of the interface, and also the surface tension are calculated at room temperature and also at a lower temperature of 273 K. The dynamical aspects of the interface are investigated in terms of the single-particle dynamical properties such as the relaxation of velocity autocorrelation and the translational diffusion coefficients along the perpendicular and parallel directions and the dipole orientational relaxation of the interfacial acetonitrile molecules. The results of the interfacial dynamics are compared with those of the corresponding bulk phases at both temperatures. The convergence of the calculated results with respect to the length of simulation runs and the system size are also discussed.

Model of Uptake of OH Radicals on Nonreactive Solids

A model of adsorption and recombination of OH radicals was developed for nonreactive solid surfaces of atmospheric interest. A parametrization of this heterogeneous mechanism was carried out to determine the role of the catalytic properties of these solid surfaces, taking into account the adsorption energy, defects, surface diffusion, and chemical reactions in the gas-solid interface. The uptake process was simulated for diffusion-controlled chemical reactions on the surface on the basis of Langmuir-Hinshelwood and Eley-Rideal mechanisms. Using an analytical approach and the Monte Carlo technique, we show the dependencies of the uptake probability of the heterogeneous reactions on the OH concentration and adsorption energy. The model is employed in the analysis of the empirically derived uptake coefficient for water ice, Al(2)O(3), NaCl, NH(4)NO(3), NH(4)HSO(4), and (NH(4))(2)SO(4). We found the following values for the free energy of adsorption of OH radicals: E(ice) = 7.3-7.6 kcal/mol, E(Al)(2)(O)(3) = 11-11.7 kcal/mol, E(NH)(4)(NO)(3) = 10.2 kcal/mol, E(NaCl) = 10.2 kcal/mol, E(NH)(4)(HSO)(4) = 9.8 kcal/mol, and E((NH)(4))(2)(SO)(4) = 9.8 kcal/mol. The atmospheric implications of the catalytic reactions of OH with adsorbed reactive molecules are discussed. The results of the modeling of the uptake process showed that the heterogeneous decay rate can exceed the corresponding gas-phase reaction rate under atmospheric conditions.

Diffusion at the Liquid−Vapor Interface of an Aqueous Ionic Solution Utilizing a Dual Simulation Technique

The recently proposed dual simulation technique [J. Phys. Chem. B 2004, 108, 6595.], with slight modification, was used to determine the diffusion coefficients for a variety of regions of a 2.2 M sodium chloride aqueous solution with a vapor-liquid interface. The diffusion of all species was shown to be isotropic far away from the interface, but at different regions in the interface, the diffusion coefficients parallel and perpendicular to the interface did not agree for water and chloride. Specifically, interfacial water diffusion parallel to the interface was significantly higher than diffusion perpendicular to the interface. Chloride ions showed even larger anisotropicity in its diffusion coefficient at the interface, with its perpendicular diffusion being similar to its bulk value, but parallel diffusion being much higher, corresponding to the region of highest chloride ion concentration. The origin for this was found to be hydrogen bonds with waters which are highly oriented perpendicular to the interface, somewhat impeding chloride ion diffusion perpendicular to the interface. While sodium ion diffusion increased at the interface, its interfacial concentration is low in that region, and its diffusion was fairly isotropic throughout all regions.