Solvation enthalpies of neutral solutes in water and octanol

The balance between electrostatic and non-electrostatic enthalpic contributions to the free energy of solvation of a series of neutral solutes in water and n-octanol is examined by means of continuum solvation calculations based on the Miertus–Scrocco–Tomasi (MST) method. The experimental data indicate that the solvation enthalpy of hydrocarbons is very similar in water and n-octanol, and that the enthalpic contribution measured for polar compounds is larger in water than in n-octanol. According to MST calculations, the different magnitude of the solvation enthalpy found for polar compounds in the two solvents can be largely attributed to the electrostatic contribution. Moreover, the results point out that there is close resemblance between the non-electrostatic components for both hydrocarbons and polar compounds in the two solvents. Finally, the results show the power of current continuum models like MST to dissect the total free energy of solvation in entropic and enthalpic contributions and suggest that new refinements of continuum solvation models should include not only the fitting to solvation free energies, but also their enthalpic components.     

Hydrophobic and ionic interactions in nanosized water droplets

A number of situations such as protein folding in confined spaces, lubrication in tight spaces, and chemical reactions in confined spaces require an understanding of water-mediated interactions. As an illustration of the profound effects of confinement on hydrophobic and ionic interactions, we investigate the solvation of methane and methane decorated with charges in spherically confined water droplets. Free energy profiles for a single methane molecule in droplets, ranging in diameter (D) from 1 to 4 nm, show that the droplet surfaces are strongly favorable as compared to the interior. From the temperature dependence of the free energy in D = 3 nm, we show that this effect is entropically driven. The potentials of mean force (PMFs) between two methane molecules show that the solvent separated minimum in the bulk is completely absent in confined water, independent of the droplet size since the solute particles are primarily associated with the droplet surface. The tendency of methanes with charges (M(q+) and M(q-) with q(+) = |q(-)| = 0.4e, where e is the electronic charge) to be pinned at the surface depends dramatically on the size of the water droplet. When D = 4 nm, the ions prefer the interior whereas for D < 4 nm the ions are localized at the surface, but with much less tendency than for methanes. Increasing the ion charge to e makes the surface strongly unfavorable. Reflecting the charge asymmetry of the water molecule, negative ions have a stronger preference for the surface compared to positive ions of the same charge magnitude. With increasing droplet size, the PMFs between M(q+) and M(q-) show decreasing influence of the boundary owing to the reduced tendency for surface solvation. We also show that as the solute charge density decreases the surface becomes less unfavorable. The implications of our results for the folding of proteins in confined spaces are outlined.     

Molecular origin of anticooperativity in hydrophobic association

The structuring of water molecules in the vicinity of nonpolar solutes is responsible for hydrophobic hydration and association thermodynamics in aqueous solutions. Here, we studied the potential of mean force (PMF) for the formation of a dimer and trimers of methane molecules in three specific configurations in explicit water to explain multibody effects in hydrophobic association on a molecular level. We analyzed the packing and orientation of water molecules in the vicinity of the solute to explain the effect of ordering of the water around nonpolar solutes on many-body interactions. Consistent with previous theoretical studies, we observed cooperativity, manifested as a reduction of the height of the desolvation barrier for the trimer in an isosceles triangle geometry, but for linear trimers, we observed only anticooperativity. A simple mechanistic picture of hydrophobic association is drawn. The free energy of hydrophobic association depends primarily on the difference in the number of water molecules in the first solvation shell of a cluster and that in the monomers of a cluster; this can be approximated by the molecular surface area. However, there are unfavorable electrostatic interactions between the water molecules from different parts of the solvation shell of a trimer because of their increased orientation induced by the nonpolar solute. These electrostatic interactions make an anticooperative contribution to the PMF, which is clearly manifested for the linear trimer where the multibody contribution due to changes in the molecular surface area is equal to zero. The information theory model of hydrophobic interactions of Hummer et al. also explains the anticooperativity of hydrophobic association of the linear trimers; however, it predicts anticooperativity with a qualitatively identical distance dependence for nonlinear trimers, which disagrees with the results of simulations.     

Homogeneous ice nucleation from aqueous inorganic/organic particles representative of biomass burning: water activity, freezing temperatures, nucleation rates

Homogeneous ice nucleation plays an important role in the formation of cirrus clouds with subsequent effects on the global radiative budget. Here we report on homogeneous ice nucleation temperatures and corresponding nucleation rate coefficients of aqueous droplets serving as surrogates of biomass burning aerosol. Micrometer-sized (NH(4))(2)SO(4)/levoglucosan droplets with mass ratios of 10:1, 1:1, 1:5, and 1:10 and aqueous multicomponent organic droplets with and without (NH(4))(2)SO(4) under typical tropospheric temperatures and relative humidities are investigated experimentally using a droplet conditioning and ice nucleation apparatus coupled to an optical microscope with image analysis. Homogeneous freezing was determined as a function of temperature and water activity, a(w), which was set at droplet preparation conditions. The ice nucleation data indicate that minor addition of (NH(4))(2)SO(4) to the aqueous organic droplets renders the temperature dependency of water activity negligible in contrast to the case of aqueous organic solution droplets. The mean homogeneous ice nucleation rate coefficient derived from 8 different aqueous droplet compositions with average diameters of ～60 μm for temperatures as low as 195 K and a(w) of 0.82-1 is 2.18 × 10(6) cm(-3) s(-1). The experimentally derived freezing temperatures and homogeneous ice nucleation rate coefficients are in agreement with predictions of the water activity-based homogeneous ice nucleation theory when taking predictive uncertainties into account. However, the presented ice nucleation data indicate that the water activity-based homogeneous ice nucleation theory overpredicts the freezing temperatures by up to 3 K and corresponding ice nucleation rate coefficients by up to ～2 orders of magnitude. A shift of 0.01 in a(w), which is well within the uncertainty of typical field and laboratory relative humidity measurements, brings experimental and predicted freezing temperatures and homogeneous ice nucleation rate coefficients into agreement. The experimentally derived ice nucleation data are applied to constrain the water activity-based homogeneous ice nucleation theory to smaller than ±1 order of magnitude compared to the predictive uncertainty of larger than ±6 orders of magnitude. The atmospheric implications of these findings are discussed.     

Charge asymmetries in hydration of polar solutes

We study the solvation of polar molecules in water. The center of water's dipole moment is offset from its steric center. In common water models, the Lennard-Jones center is closer to the negatively charged oxygen than to the positively charged hydrogens. This asymmetry of water's charge sites leads to different hydration free energies of positive versus negative ions of the same size. Here, we explore these hydration effects for some hypothetical neutral solutes, and two real solutes, with molecular dynamics simulations using several different water models. We find that, like ions, polar solutes are solvated differently in water depending on the sign of the partial charges. Solutes having a large negative charge balancing diffuse positive charges are preferentially solvated relative to those having a large positive charge balancing diffuse negative charges. Asymmetries in hydration free energies can be as large as 10 kcal/mol for neutral benzene-sized solutes. These asymmetries are mainly enthalpic, arising primarily from the first solvation shell water structure. Such effects are not readily captured by implicit solvent models, which respond symmetrically with respect to charge.     

Influence of oil polarity on droplet growth in oil-in-water emulsions stabilized by a weakly adsorbing biopolymer or a nonionic surfactant

The influence of oil type (n-hexadecane, 1-decanol, n-decane), droplet composition (hexadecane:decanol), and emulsifier type (Tween 20, gum arabic) on droplet growth in oil-in-water emulsions was studied. Droplet size distributions of emulsions were measured over time (0-120 h) by laser diffraction and ultrasonic spectroscopy. Emulsions containing oil molecules of low polarity and low water solubility (hexadecane) were stable to droplet growth, irrespective of the emulsifier used to stabilize the droplets. Emulsions containing oil molecules of low polarity and relatively high water solubility (decane) were stable to coalescence, but unstable to Ostwald ripening, irrespective of emulsifier. Droplet growth in emulsions containing oil molecules of relatively high polarity and high water solubility (decanol) depended on emulsifier type. Decanol droplets stabilized by Tween 20 were stable to droplet growth in concentrated emulsions but unstable when the emulsions were diluted. Decanol droplets stabilized by gum arabic exhibited rapid and extensive droplet growth, probably due to a combination of Ostwald ripening and coalescence. We proposed that coalescence was caused by the relatively low interfacial tension at the decanol-water boundary, which meant that the gum arabic did not absorb strongly to the droplet surfaces and therefore did not prevent the droplets from coming into close proximity.     

Interfacial water on hydrophobic surfaces recognized by ions and molecules

Recent spectrophotometric and molecular dynamics simulation studies have shown that the physicochemical properties and structures of water in the vicinity of hydrophobic surfaces differ from those of the bulk water. However, the interfacial water acting as a separation medium on hydrophobic surfaces has never been detected and quantified experimentally. In this study, we show that small inorganic ions and organic molecules differentiate the interfacial water formed on the surfaces of octadecyl-bonded (C(18)) silica particles from the bulk water and the chemical separation of these solutes in aqueous media with hydrophobic materials can be interpreted with a consistent mechanism, partition between the bulk water phase and the interfacial water formed on the hydrophobic surface. Thermal transition behaviour of the interfacial water incorporated in the nanopores of the C(18) silica materials and the solubility parameter of the water calculated from the distribution coefficients of organic compounds have indicated that the interfacial water may have a structure of disrupted hydrogen bonding. The thickness of the interfacial water or the limit of distance from the hydrophobic surface at which molecules and ions can sense the surface was estimated to be 1.25 ± 0.13 nm from the volume of the interfacial water obtained by a liquid chromatographic method and the surface area, suggesting that the hydrophobic effect may extend beyond the first solvation shell of water molecules directly surrounding the surfaces.     

Mobile phase effects in reversed-phase liquid chromatography: a comparison of acetonitrile/water and methanol/water solvents as studied by molecular simulation

Molecular simulations of water/acetonitrile and water/methanol mobile phases in contact with a C(18) stationary phase were carried out to examine the molecular-level effects of mobile phase composition on structure and retention in reversed-phase liquid chromatography. The simulations indicate that increases in the fraction of organic modifier increase the amount of solvent penetration into the stationary phase and that this intercalated solvent increases chain alignment. This effect is slightly more apparent for acetonitrile containing solvents. The retention mechanism of alkane solutes showed contributions from both partitioning and adsorption. Despite changes in chain structure and solvation, the molecular mechanism of retention for alkane solutes was not affected by solvent composition. The mechanism of retention for alcohol solutes was primarily adsorption at the interface between the mobile and stationary phase, but there were also contributions from interactions with surface silanols. The interaction between the solute and surface silanols become very important at high concentrations of acetonitrile.     

Solvation model based on order parameters and a fast sampling method for the calculation of the solvation free energies of peptides

An analytical solvation model is proposed as a function of an order parameter, which represents the local arrangement of water molecules in the first solvation shell of peptide atoms. The model is combined with a fast sampling method, rotational isomeric state Monte Carlo, to sample efficiently the torsional degrees of freedom on a peptide backbone. This order parameter solvation model is shown to reproduce without ad hoc fitting parameters the solvation free energies of single amino acids and tripeptides with slightly better accuracy than the generalized Born model but with several orders of magnitude improvement in efficiency. This method is a potential candidate for efficiently and accurately tackling some important issues in biophysical chemistry that are related to solvation, for example, protein folding, ligand binding, etc. Our results also present fundamental new insights into solvation. Specifically, the local water geometry, represented in this work by a properly defined order parameter, carries the majority, if not all, of the energetic information of solvation, including solute-solvent interactions and solvent reorganization in the presence of the solute.     

Evaporation of sessile droplets of dilute aqueous solutions containing sodium n-alkylates from polymer surfaces: influences of alkyl length and concentration of solute

The evaporation of sessile droplets placed on polymer surfaces was studied by microscopic observation of the changes in shape of aqueous solution droplets in which the alkyl lengths and the initial concentrations of sodium n-alkylates were varied. Although the initial contact angles of the droplets were not significantly different, the evaporation process varied significantly with the alkyl length of the sodium n-alkylate employed. For the sodium dodecanoate (C 12), showing the highest surface activity, the concentration was found to have a significant effect on the evaporation process of the droplets. In the evaporation of water droplets, variations in the three distinct stages were caused by the different concentration of solutes distributed near or at the air/water interface. It is revealed that the concentration of droplet solute near the air/water interface requires not only solvent evaporation but also some affinity of the solute for the interface. The initial C 12 concentration-dependence of the evaporation of C 12 solution droplets is discussed with particular emphasis on the sudden spreading or sudden contraction of the contact area near the end of evaporation. It is suggested that the cluster formation by C 12 molecules at the air/liquid interface during the evaporation causes Marangoni instability in an evaporating droplet, and the clusters are expected to move dynamically, depending on the droplet concentration of C 12, from the droplet center to the contact line and vice versa, showing Marangoni flow along the air/water interface.     

Enthalpy-entropy contributions to the potential of mean force of nanoscopic hydrophobic solutes

Entropic and enthalpic contributions to the hydrophobic interaction between nanoscopic hydrophobic solutes, modeled as graphene plates in water, have been calculated using molecular dynamics simulations in the isothermal-isobaric (NPT) ensemble with free energy perturbation methodology. We find the stabilizing contribution to the free energy of association (contact pair formation) to be the favorable entropic part, the enthalpic contribution being highly unfavorable. The desolvation barrier is dominated by the unfavorable enthalpic contribution, despite a fairly large favorable entropic compensation. The enthalpic contributions, incorporating the Lennard-Jones solute-solvent terms, largely determine the stability of the solvent separated configuration. We decompose the enthalpy into a direct solute-solute term, the solute-solvent interactions, and the remainder that contains pressure-volume work as well as contributions due to solvent reorganization. The enthalpic contribution due to changes in water-water interactions arising from solvent reorganization around the solute molecules is shown to have major contribution in the solvent induced enthalpy change.     

Effects of atmospherically important solvated ions on organic acid adsorption at the surface of aqueous solutions

The effects of salts on the solubility of amphiphilic organic molecules are of importance to numerous atmospheric, environmental, and biological systems. A detailed picture of the influence of dissolved atmospheric salts, NaCl and Na(2)SO(4), on the adsorption of hexanoic acid at the vapor/water interface is developed using vibrational sum-frequency spectroscopy and surface tension measurements as a function of time, organic concentration, and solution pH. We have found that for hexanoic acid adsorption at the vapor/water interface, a fast initial adsorption is followed by two considerably slower processes: a reorientation of the polar headgroup and a restructuring of the headgroup solvation shell. The addition of salts affects this restructuring by reducing the range of water--headgroup interactions immediately upon surface adsorption for ion containing solutions. Reorientation of the organic headgroup with time occurs at the surface of both salt-containing and salt-free solutions, but the most stable orientation differs with the added ions. The dissolved salts also enhance the interfacial concentration of hexanoic acid, consistent with the known salting-out behavior of Cl(-) and SO(4)(2-) anions.     

Adsorption of polycyclic aromatic hydrocarbons at the air-water interface: molecular dynamics simulations and experimental atmospheric observations

Adsorption of benzene, naphthalene, anthracene, and phenanthrene at the aqueous surface is investigated by means of molecular dynamics simulations. Potentials of mean force, i.e., free energy profiles obtained when moving the studied molecules across an aqueous slab were evaluated. In all cases, deep surface free energy minima, corresponding to orders of magnitude of surface enhancement of the aromatic molecule, were located. This enhancement, which increases with the size of the solute, points to the importance of the aqueous surface for the chemistry of polycyclic aromatic hydrocarbons (PAHs). Supporting evidence in the atmospheric environment related to the heterogeneous chemistry of PAHs on water droplets and planar surfaces is summarized. There is good agreement between the hydration free energies computed from MD calculations and the experimentally determined values. Data pertaining to the importance of the air-water interface in the adsorption and transport of PAHs on micron sized water droplets are described. The relevant data on adsorption and reaction (ozonation and photochemical) at the air-water interface of planar surfaces and droplets are also summarized.     

Solvent screening for a hard-to-dissolve molecular crystal

Materials with a high-degree of inter- and intra-molecular hydrogen bonding generally have limited solubility in conventional organic solvents. This presents a problem for the dissolution, manipulation and purification of these materials. Using a state-of-the-art density-functional-theory based quantum chemical solvation model we systematically evaluated solvents for a known hydrogen-bonded molecular crystal. This, coupled with direct solubility measurements, uncovered a class of ionic liquids involving fluoride anions that possess more than two orders of magnitude higher solvation power as compared with the best conventional solvents. The crystal structure of one such ionic liquid, determined by X-ray diffraction spectroscopy, indicates that F(-) ions are stabilized through H-bonded chains with water. The presence of coordinating water in such ionic liquids seems to facilitate the dissolution process by keeping the chemical activity of the F(-) ions in check.     

Multiply Charged Helium Droplet Anions

The detection of multiply charged helium droplet anions is reported for the first time. By ionizing droplets of superfluid helium with low energy electrons (up to 25 eV), it was possible to produce droplets containing up to five negative charges, which remain intact on the timescale of the experiment. The appearance sizes for different charge states are determined and are found to be orders of magnitude larger than for the equivalent cationic droplets, starting at 4 million He atoms for dianions. Droplets with He*⁻ as charge carriers show signs of being metastable, but this effect is quenched by the pickup of water molecules.     

Free energy of solvation from molecular dynamics simulation applying Voronoi-Delaunay triangulation to the cavity creation

The free energy of solvation for a large number of representative solutes in various solvents has been calculated from the polarizable continuum model coupled to molecular dynamics computer simulation. A new algorithm based on the Voronoi-Delaunay triangulation of atom-atom contact points between the solute and the solvent molecules is presented for the estimation of the solvent-accessible surface surrounding the solute. The volume of the inscribed cavity is used to rescale the cavitational contribution to the solvation free energy for each atom of the solute atom within scaled particle theory. The computation of the electrostatic free energy of solvation is performed using the Voronoi-Delaunay surface around the solute as the boundary for the polarizable continuum model. Additional short-range contributions to the solvation free energy are included directly from the solute-solvent force field for the van der Waals-type interactions. Calculated solvation free energies for neutral molecules dissolved in benzene, water, CCl4, and octanol are compared with experimental data. We found an excellent correlation between the experimental and computed free energies of solvation for all the solvents. In addition, the employed algorithm for the cavity creation by Voronoi-Delaunay triangulation is compared with the GEPOL algorithm and is shown to predict more accurate free energies of solvation, especially in solvents composed by molecules with nonspherical molecular shapes.     

Laser tweezers Raman study of optically trapped aerosol droplets of seawater and oleic acid reacting with ozone: implications for cloud-droplet properties

In this communication we report the first use of the Raman laser tweezers technique to trap and hold a mixed droplet of oleic acid and water at atmospheric pressure for 30 min, oxidize the oleic acid on the droplet, follow the decay of reactants and the growth of chemical products using Raman spectroscopy, and monitor the growth in size of the droplet as it becomes more hydrophilic. We demonstrate that the oxidation of organic films on water droplets could have large climatic effects in the atmosphere. We show that cloud-droplet growth and activation of cloud condensation nuclei (to become cloud droplets) is retarded by the presence of an organic film and that chemical oxidation of this film would allow a cloud droplet to grow, reducing cloud albedo by inducing precipitation, and would allow a cloud condensation nucleus to grow to a cloud droplet, thus forming a cloud and increasing the albedo of the Earth.     

SURFACE ENTHALPY AND ENTROPY AND THE PHYSICO-CHEMICAL NATURE OF HYDROPHOBIC AND HYDROPHILIC INTERACTIONS

The attractive Interactions between typically hydrophobic molecules such as hexane or CCl₄, and the repulsive Interactions between extremely hydrophilic molecules such as poly(ethylene oxide) (PEO), when immersed in water, as well as the interactions between these molecules and water, have been examined from a surface thermodynamic viewpoint, taking the changes in surface free energy into account, as a function of temperature. It was found that attractive hydrophobic Interactions are not, as vas generally believed up to now, invariably entropic. Hydrophobic Interactions can be mainly enthalpic or mainly entropic, or more or less equal mixtures of both, depending on each individual case; however, all hydrophobic interactions are polar (in the sense of Lewis acid-base) in nature. Repulsive hydrophilic interactions are enthalpic, and also polar in nature. The interaction between hydrophobic solutes and water is mainly enthalpic, and is apolar in nature.     

Controlled shrinkage and re-expansion of a single aqueous droplet inside an optical vortex trap

This paper describes the shrinkage and re-expansion of individual femtoliter-volume aqueous droplets that were suspended in an organic medium and held in an optical vortex trap. To elucidate the mechanism behind this phenomenon, we constructed a heat- and mass-transfer model and carried out experimental verifications of our model. From these studies, we conclude that an evaporation mechanism sufficiently describes the shrinkage of aqueous droplets held in a vortex trap, whereas a mechanism based on the supersaturation of the organic phase by water that surrounds the droplet adequately explains the re-expansion of the shrunk droplet. The proposed mechanisms correlated well with experimental observations using different organic media, when H2O was replaced with D2O and when an optical tweezer was used to induce droplet shrinkage rather than an optical vortex trap. For H2O droplets, the temperature rise within the droplet during shrinkage was on the order of 1 K or less, owing to the rapid thermal conduction of heat away from the droplet at the microscale and the sharp increase in solubility for water by the organic phase with slight elevations in temperature. Because most chemical species confined to droplets can be made impenetrable to the aqueous/organic interface, a change in the volume of aqueous droplets translates into a change in concentration of the dissolved species within the droplets. Therefore, this phenomenon should find use in the study of fundamental chemical processes that are sensitive to concentration, such as macromolecular crowding and protein nucleation and crystallization.     

Molecular dynamic simulation of dicarboxylic acid coated aqueous aerosol: structure and processing of water vapor

Organic monolayers at the surfaces of aqueous aerosols play an important role in determining the mass, heat transfer rate and surface reactivity of atmospheric aerosols. They can potentially contribute to the formation of cloud condensation nuclei (CCN) and are involved in a series of chemical reactions occurring in atmosphere. Recent studies even suggest that organic-coated interfaces could have played some role in prebiotic biochemistry and the origin of life. However, creating reproducible, well-characterized aqueous aerosol particles coated with organic films is an experimental challenge. This opens the opportunity for computer simulations and modeling of these complex structures. In this work, molecular dynamics simulation was used to probe the structure and the interfacial properties of the dicarboxylic acid coated aqueous aerosol. Low molecular weight dicarboxylic acids of various chain lengths and water solubility were chosen to coat a water droplet consisting of 2440 water molecules. For malonic acid coated aerosol, the surface acid molecules dissolved into the water core and formed an ordered structure due to the hydrophobic interactions. The acid and the water are separated inside the aerosol. For other nanoaerosols coated with low solubility acids, phase separation between water and acid molecules was observed on the surface of the particle. To study the water processing of the coated aerosols, the water vapor accommodation factors were calculated.