Partition Constant of Binary Mixtures for the Equilibrium between a Bulk and a Confined Phase

It is well-known that the thermodynamic, kinetic and structural properties of fluids, and in particular of water and its solutions, can be drastically affected in nanospaces. A possible consequence of nanoscale confinement of a solution is the partial segregation of its components. Thereby, confinement in nanoporous materials (NPM) has been proposed as a means for the separation of mixtures. In fact, separation science can take great advantage of NPM due to the tunability of their properties as a function of nanostructure, morphology, pore size, and surface chemistry. Alcohol-water mixtures are in this context among the most relevant systems. However, a quantitative thermodynamic description allowing for the prediction of the segregation capabilities as a function of the material-solution characteristics is missing. In the present study we attempt to fill this vacancy, by contributing a thermodynamic treatment for the calculation of the partition coefficient in confinement. Combining the multilayer adsorption model for binary mixtures with the Young equation, we conclude that the liquid-vapor surface tension and the contact angle of the pure substances can be used to predict the separation ability of a particular material for a given mixture to a semiquantitative extent. Moreover, we develop a Kelvin-type equation that relates the partition coefficient to the radius of the pore, the contact angle, and the liquid-vapor surface tensions of the constituents. To assess the validity of our thermodynamic formulation, coarse grained molecular dynamics simulations were performed on models of alcohol-water mixtures confined in cylindrical pores. To this end, a coarse-grained amphiphilic molecule was parametrized to be used in conjunction with the mW potential for water. This amphiphilic model reproduces some of the properties of methanol such as enthalpy of vaporization and liquid-vapor surface tension, and the minimum of the excess enthalpy for the aqueous solution. The partition coefficient turns out to be highly dependent on the molar fraction, on the interaction between the components and the confining matrix, and on the radius of the pore. A remarkable agreement between the theory and the simulations is found for pores of radius larger than 15 Å.

Grand Canonical Investigation of the Quasi Liquid Layer of Ice: Is It Liquid?

In this study, the solid-vapor equilibrium and the quasi liquid layer (QLL) of ice Ih exposing the basal and primary prismatic faces were explored by means of grand canonical molecular dynamics simulations with the monatomic mW potential. For this model, the solid-vapor equilibrium was found to follow the Clausius-Clapeyron relation in the range examined, from 250 to 270 K, with a Δ H_sub of 50 kJ/mol in excellent agreement with the experimental value. The phase diagram of the mW model was constructed for the low pressure region around the triple point. The analysis of the crystallization dynamics during condensation and evaporation revealed that, for the basal face, both processes are highly activated, and in particular cubic ice is formed during condensation, producing stacking-disordered ice. The basal and primary prismatic surfaces of ice Ih were investigated at different temperatures and at their corresponding equilibrium vapor pressures. Our results show that the region known as QLL can be interpreted as the outermost layers of the solid where a partial melting takes place. Solid islands in the nanometer length scale are surrounded by interconnected liquid areas, generating a bidimensional nanophase segregation that spans throughout the entire width of the outermost layer even at 250 K. Two approaches were adopted to quantify the QLL and discussed in light of their ability to reflect this nanophase segregation phenomena. Our results in the μVT ensemble were compared with NPT and NVT simulations for two system sizes. No significant differences were found between the results as a consequence of model system size or of the working ensemble. Nevertheless, certain advantages of performing μVT simulations in order to reproduce the experimental situation are highlighted. On the one hand, the QLL thickness measured out of equilibrium might be affected because of crystallization being slower than condensation. On the other, preliminary simulations of AFM indentation experiments show that the tip can induce capillary condensation over the ice surface, enlarging the apparent QLL.

Stability and Vapor Pressure of Aqueous Aggregates and Aerosols Containing a Monovalent Ion

The incidence of charged particles on the nucleation and the stability of aqueous aggregates and aerosols was reported more than a century ago. Many studies have been conducted ever since to characterize the stability, structure, and nucleation barrier of ion-water droplets. Most of these studies have focused on the free-energy surface as a function of cluster size, with an emphasis on the role of ionic charge and radius. This knowledge is fundamental to go beyond the rudimentary ion-induced classical nucleation theory. In the present article, we address this problem from a different perspective, by computing the vapor pressures of (H₂O)_nLi⁺ and (H₂O)_nCl^- aggregates using molecular simulations. Our calculations shed light on the structure, the critical size, the range of stability, and the role of ion-water interactions in aqueous clusters. Moreover, they allow one to assess the accuracy of the classical thermodynamic model, highlighting its strengths and weaknesses.

Hydrogen-Bond Heterogeneity Boosts Hydrophobicity of Solid Interfaces

Experimental and theoretical studies suggest that the hydrophobicity of chemically heterogeneous surfaces may present important nonlinearities as a function of composition. In this article, this issue is systematically explored using molecular simulations. The hydrophobicity is characterized by computing the contact angle of water on flat interfaces and the desorption pressure of water from cylindrical nanopores. The studied interfaces are binary mixtures of hydrophilic and hydrophobic sites, with and without the ability to form hydrogen bonds with water, intercalated at different scales. Water is described with the mW coarse-grained potential, where hydrogen-bonds are modeled in the absence of explicit hydrogen atoms, via a three-body term that favors tetrahedral coordination. We found that the combination of particles exhibiting the same kind of coordination with water gives rise to a linear dependence of contact angle with respect to composition, in agreement with the Cassie model. However, when only the hydrophilic component can form hydrogen bonds, unprecedented deviations from linearity are observed, increasing the contact angle and the vapor pressure above their values in the purely hydrophobic interface. In particular, the maximum enhancement is seen when a 35% of hydrogen bonding molecules is randomly scattered on a hydrophobic background. This effect is very sensitive to the heterogeneity length-scale, being significantly attenuated when the hydrophilic domains reach a size of 2 nm. The observed behavior may be qualitatively rationalized via a simple modification of the Cassie model, by assuming a different microrugosity for hydrogen bonding and non-hydrogen bonding interfaces.