Molecular simulation of binary phase diagrams from the osmotic equilibrium method: vapour pressure and activity in water–ethanol mixtures

ABSTRACT

Herein, an approach for simulating phase diagrams of binary mixtures is presented, where a bulk liquid and its corresponding vapour phase are simulated by means of molecular dynamics using explicit polarisation. Time-averaged density profiles for the pure compounds and mixtures at different mole fractions provide information about the spatial distribution in the bulk liquid and the amount of evaporated species in the adjacent vapour phase. The activities in the liquid phase are calculated from the mean vapour phase densities at a given composition, providing a good qualitative agreement compared to experimental data and the precision of the method follows a previously developed Poisson model of evaporation. With the Redlich–Kister approach for the activities in a binary mixture, the directly obtained activities are fitted providing corrected activity coefficients of the two species. This method is applied to ethanol water mixtures at different mole fractions. The obtained structural data are in good agreement with experimental data and time-averaged density profiles provide a detailed insight into the composition of the liquid–vapour interface. An azeotropic point is obtained for an excess concentration of ethanol at 87% as percentage by mass compared to the experimental value of 95%.     

Surface tension of hydrocarbon chains at the liquid–vapour interface

Molecular dynamics simulations were performed at constant temperature to obtain the surface tension of hydrocarbon chains at the liquid–vapour interface. The Ewald sum was used to calculate the dispersion forces of the Lennard–Jones potential to take into account the full interaction. The NERD and TraPPE_UA flexible force field models were used to simulate molecules from ethane to hexadecane along the coexistence curve. The simulation results for the TraPPE_UA model are in good agreement with experimental data, whereas the NERD model predicts slightly higher values.     

Structure and dynamics of water at liquid–vapour interfaces covered by surfactant monolayers of neutral stearic acid and charged stearate ions

The behaviour of water molecules at liquid–vapour interfaces with a surfactant monolayer of either stearic acid molecules or anionic stearate ions is investigated by means of molecular dynamics simulations. The density and dipolar orientational profiles and also the dynamics of translational and rotational motion of interfacial water molecules are calculated in the present work and the results are compared with the bulk liquid water and also of liquid–vapour interface of surfactant-free water. The present simulation results are also compared with available experimental results of similar interfacial systems with a monolayer of either neutral or ionic surfactants.     

Nanopore wall effect on surface tension of methane

Local pressure is known to be anisotropic across the interfaces separating fluids in equilibrium. Tangential pressure profiles show characteristic negative peaks as a result of surface tension forces parallel to the interface. Nearby attractive forces parallel to the interface are larger than the repulsive forces and, hence, constitute the surface tension. In this work, using molecular dynamics simulations of methane inside nano-scale pores, we show this surface tension behaviour could be significantly influenced by confinement effects. The layering structure, characterised by damped oscillations in local liquid density and tangential pressures, extends deep into the pore and can be a few nanometers thick. The surface tension is measured numerically using local pressures across the interface. Results show that the tension is smaller under confinement and becomes a variable in small pores, mainly controlled by the thickness of the liquid density layering (or liquid saturation) and the pore width. If the liquid saturation inside the pore is high enough, the vapour–liquid interface is not interfered by the pore wall and the surface tension remains the same as the bulk values. The results are important for understanding phase change and multi-phase transport phenomena in nanoporous materials.     

Perspectives on water science: transport and application of confined water

The confinements of water can be divided into two main categories,namely,the confinements on surface or interface and the confinements in bulk water.By adding ions or applying electric field,the intensity and distribution of the hydrogen bonds can be greatly affected.These are collectively known as confinement on water surface or interface,which has potential applications in life science and industries involving evaporation control.Confined bulk water could be found everywhere in nature,such as in granular and porous materials,macromolecules and gels,etc.The investigation of the physical properties and the transports of the confined bulk water will contribute to understanding certain types of life activities such as the water transport in plant and in new application of extracting the shale oil and water.     

The dipole moment distribution on water is improved by using large flexibility and large bending angle

A highly flexible model of water with fixed charges is used to study properties of water. The bending angle of an isolated molecule is 125^○ that was chosen to match the experimental dipole moment. The geometry of water in the liquid phase is made closer to that of the rigid SPC/E model by decreasing the bending angle spring constant, k _Θ. The new model, called SPCE-FHΘ, is a modified version of the recently proposed SPCE-FH [J. Alejandre, G.A. Chapela, F. Bresme and J.-P. Hansen, J. Chem. Phys. 130, 174505 (2009)] to simulate ionic solutions which includes short ranged interactions on the hydrogen atoms. By increasing angle flexibility it is possible to obtain, in the liquid phase at ambient conditions, bending angles ?Θ(HOH)? ～ 109^○, dipole moment ?μ? ～ 2.5 D and dielectric constant ?ε? ～ 80. The dipole moment distribution at room temperature goes from 1.5 to 3.5 D due to large fluctuations in bending angle and has the same trend found in ab initio simulations of liquid water. The dipole moment profile at the interface of water varies from 1.9 D in the vapour phase to 2.5 D in the liquid region at 400 K. The SPCE-FHΘ gives dipole moment, dielectric constant, coexisting densities and surface tension along the liquid–vapour coexistence line closer to the experimental values than those obtained for the SPC/E force field.     

Force-field dependence on the interfacial structure of oil–water interfaces

We investigate the performance of different force-fields for alkanes, united (TraPPE) and all atom (OPLS-AA) models, and water (SPC/E and TIP4P-2005), in the prediction of the interfacial structure of alkane (n-octane, and n-dodecane)–water interfaces. We report an extensive comparison of the interfacial thermodynamic properties as well as the interfacial structure (translational and orientational). We use the recently introduced intrinsic sampling method, which removes the averaging effect of the interfacial capillary waves and provides a clear view of the interface structure. The alkane interfacial structure is sensitive to the environment, i.e. alkane–vapour or alkane–water interfaces, showing a stronger structure when it is in contact with the water phase. We find that this structure is fairly independent of the level of detail, full or united atom, employed to describe the alkane phase. The water surface properties show a small dependence on the water model. The dipole moment of the SPC/E model shows asymmetric fluctuations, with a tendency to point both towards the alkane and water phases. On the other hand the dipole moment of the TIP4P-2005 model shows a tendency to point towards the water phase only. Analysis of the intrinsic electrostatic field indicates that the surface water potential is confined to an interfacial region of about 8 Å. Overall we find that the intrinsic structure of alkane–water interfaces is a robust interfacial property, which is independent of the details of the force-field employed. Hence, it should provide a good reference to interpret experimental data.     

Surface triple points and multiple-layer transitions observed by tuning the surface field at smectic liquid-crystal-water interfaces

We present an ellipsometric study of the interface between a smectic liquid crystal and water in the presence of a nonionic surfactant. The surfactant concentration serves as a handle to tune the surface field. For sufficiently large surfactant concentrations, a smectic phase is present at the interface in the temperature range above the smectic-A-isotropic bulk transition; when the bulk transition is approached, the thickness of this surface phase grows via a series of layer-by-layer transitions at which single smectic layers are formed. At lower surfactant concentrations, transitions appear at which the thickness of the surface phase jumps by multiple smectic layers, thereby implying the existence of triple points at which surface phases with different smectic layer numbers coexist. This is the first experimental demonstration of such surface triple points which are predicted by theoretical models.     

Structure of fluid interfaces: an integral equation study

We develop a method for treating an interface between coexisting fluid phases to study its structure and thermodynamical properties. Our approach includes calculated correlation functions. Approximations are intrinsically optimized. The method is successfully applied to the liquid vapour surface of a Lennard-Jones fluid, to the liquid–liquid interface of a demixing system and to the surface of a Stockmayer liquid towards vapour. In all cases we compare our results to simulations. The agreement with simulation results demonstrates the reliability of our approximations and shows that the technique we applied provides a powerful and robust method for studying inhomogeneous fluids.     

Structure and fluctuations of the premelted liquid film of ice at the triple point

ABSTRACT

In this paper, we study the structure of the ice/vapour interface in the neighbourhood of the triple point for the TIP4P/2005 model. We probe the fluctuations of the ice/film and film/vapour surfaces that separate the liquid film from the coexisting bulk phases at basal, primary prismatic and secondary prismatic planes. The results are interpreted using a coupled sine Gordon plus Interface Hamiltonian model. At large length scales, the two bounding surfaces are correlated and behave as a single complex ice/vapour interface. For small length, on the contrary, the ice/film and film/vapour surfaces behave very much like independent ice/water and water/vapour interfaces. The study suggests that the basal facet of the TIP4P/2005 model is smooth, the prismatic facet is close to a roughening transition, and the secondary prismatic facet is rough. For the faceted basal face, our fluctuation analysis allows us to estimate the step free energy in good agreement with experiment. Our results allow for a quantitative characterisation of the extent to which the adsorbed quasi-liquid layer behaves as water and explains experimental observations which reveal similar activation energies for crystals grown in bulk vapour or bulk water.     

Liquid–vapour equilibrium of n -alkanes using interface simulations

Molecular Dynamics simulations were performed to calculate liquid–vapour coexisting properties of n-alkane chains up to 16 carbon atoms using interface simulations. The lattice sum or Ewald method on the dispersion forces of the Lennard–Jones potential was applied to calculate the full interaction. The liquid and vapour coexisting densities were obtained for two flexible force field models, NERD and TraPPE-UA, where the intermolecular interactions are of the Lennard–Jones type. We have recently shown [P. Orea, J. López-Lemus, and J. Alejandre, J. Chem. Phys. 123, 114702 (2005)] that the liquid–vapour densities for simple fluids do not depend on interfacial area and therefore it is possible to use a small number of molecules in a simulation. We show that the same trend is found on the simulation of these hydrocarbon molecules. The phase diagram of ethane/n-decane binary mixtures is also obtained at 410.95 K for the NERD model. The simulation results from this work were compared with those obtained using methods with interfaces using large cut-off distances and with methods without interfaces for the same potential model. In both comparisons, excellent agreement was found. The results of liquid density from the TraPPE-UA model are in good agreement with experimental data while those from the NERD model are underestimated at low temperatures. Our findings are consistent with results published by other authors for small hydrocarbons.     

Probing interfacial dynamics of water in confined nanoporous systems by NMRD

The confined dynamics of water molecules inside a pore involves an intermittence between adsorption steps near the interface and surface diffusion and excursions in the pore network. Depending on the strength of the interaction in the layer(s) close to the surface and the dynamical confinement of the distal bulk liquid, exchange dynamics can vary significantly. The average time spent in the surface proximal region (also called the adsorption layer) between a first entry and a consecutive exit allows estimating the level of ‘nanowettablity’ of water. As shown in several seminal works, NMRD is an efficient experimental method to follow such intermittent dynamics close to an interface. In this paper, the intermittent dynamics of a confined fluid inside nanoporous materials is discussed. Special attention is devoted to the interplay between bulk diffusion, adsorption and surface diffusion on curved pore interfaces. Considering the nano or meso length scale confinement of the pore network, an analytical model for calculating the inter-dipolar spin–lattice relaxation dispersion curves is proposed. In the low-frequency regime (50?KHz–100?MHz), this model is successfully compared with numerical simulations performed using a 3D-off lattice reconstruction of Vycor glass. Comparison with experimental data available in the literature is finally discussed.     

Molecular Dynamics Simulation of 4-N-Hexyl-4'-Cyanobiphenyl Adsorbed at the Air-Water Interface

The interfacial behavior of 4-n-hexyl-4'-cyanobiphenyl(6CB) molecules at the air-water interface is investigated by full atomistic molecular dynamics simulations. To understand the morphology and the structure of adsorbed 6CB molecules in detail, the snapshots and mass density profiles of the simulation system are generated. The average tilt angles between the interface normal and various vectors defined in the rigid and alkyl parts of 6CB are in good agreement with the experimental data available. The interfacial thickness and monolayer width are obtained from the mass density profiles of water and 6CB phase, respectively. The second and fourth rank orientational order parameters of cyanobiphenyl core are found to be larger than those of an elastic alkyl chain.Bond order parameters for 6CB are also calculated. The calculated oxygen-oxygen radial distribution function and hydrogen bonding statistics for bulk water are compared with those for the interfacial region. The surface tensions of the systems are calculated. All simulation results are compared with the available literature data.     

Gas—liquid interfaces of room temperature ionic liquids

The structure of the gas-liquid surface of dimethylimidazolium chloride has been studied using atomistic simulation. We find that there is a region of enhanced density immediately below the interface in which the cations are oriented with their planes perpendicular to the surface and their dipoles in the surface plane. There is negligible segregation of cations and anions. The temperature dependence of the surface tension is predicted to be anomalously low or be reversed in sign. The vapour-liquid interfaces between mixtures of water and dimethylimidazolium chloride show similar regions of enhanced density and preferential orientation of the cations. Water molecules also show preferential orientation in the interface region and are preferentially adsorbed on the vapour side of the interface. The surface tension decreases with increase in the mole fraction of water.     

Calculation of the surface tension and the surface energy of Lennard–Jones fluids from the radial distribution function in the interface zone

A detailed study is presented of the calculation of the surface tension and the surface energy of Lennard–Jones fluids from the radial distribution function and the density profile. To do so, a modification is made to Lekner and Henderson's statistical mechanics approach by introducing two simple analytical expressions for the radial distribution function of the interface zone. In these expressions the radial distribution functions of the liquid and vapour phases are weighted via step or exponential variations. The well- known exponential model for the density profile in the interface zone is considered. Finally, results are compared with values from experiment, from computer simulation and from relevant theoretical developments. It is shown that the use of the proposed radial distribution function in the interface zone represents a significant improvement in applying Lekner and Henderson's approach.     

Near-field optical study of photorefractive surface waves in BaTiO(3)

Near-field optical microscopy has been used to study photorefractive surface waves in BaTiO(3). The field distribution of the photorefractive surface wave near the crystal-air interface has been measured and compared with theory. Experimental data indicate that a micrometer-wide transition layer with dielectric and photorefractive properties that are different from the properties of the bulk BaTiO(3) exists near the interface.     

Surface behaviour and undercooling in the liquid Ga—Bi binary system detected by optical second harmonic generation

We employ an optical second harmonic generation(SHG) technique to investigate the surface behaviours at the liquid(solid)/vapour interface of the Ga-Bi binary metallic system. In a heating and cooling cycle between 280℃ and room temperature, there is no change of the SH-intensity in the heating process, whereas there exists an abrupt and abnormal change of the SH-intensity in the cooling process. It is interesting to find that a macroscopic Bi-rich solid layer is floating on the surface of the Ga-rich liquid phase just below the monotectic temperature (222℃±2℃) in the cooling process, in spite of the Bi-rich phase being heavier than the Ga-rich phase. On the other hand, different undercooling behaviours are observed at the surface and in the bulk. The behaviours of surface solidification and surface melting are different from those in the bulk.     

Electric field dependence of phase equilibria of binary Stockmayer fluid mixtures

A perturbation theory based study of the effect of an external electric field on the phase equilibrium properties of binary Stockmayer fluids is presented. The dipole–dipole interaction and the applied field are treated as independent perturbations to a Lennard–Jones mixture, and the reference fluid is treated by the van der Waals one-fluid approximation. A third-order free energy expression in the electric field strength is established, and the dielectric constant is calculated for a needle-shaped sample parallel to the field direction. We present and discuss vapour–liquid and liquid–liquid equilibrium curves at a given temperature for some dipolar mixtures exposed to an electric field, including chlorodifluoromethane +?difluoromethane and acetonitrile +?methanol. A sufficiently high electric field may result in massive shifts of vapour pressures and critical or azeotropic points, and can considerably alter the properties of coexisting phases. The vapour pressure decreases with increasing field strength.     

The vicinity of an equilibrium three-phase contact line using density-functional theory: density profiles normal to the fluid interface

The paper by Nold et al. [Phys. Fluids 26 (7), 072001 (2014)] examined density profiles and the micro-scale structure of an equilibrium three-phase (liquid–vapour–solid) contact line in the immediate vicinity of the wall using elements from the statistical mechanics of classical fluids, namely density-functional theory. The present research note, building on the above work, further contributes to our understanding of the nanoscale structure of a contact line by quantifying the strong dependence of the liquid–vapour density profile on the normal distance to the interface, when compared to the dependence on the vertical distance to the substrate. A recent study by Benet et al. [J. Phys. Chem. C 118 (38), 22079 (2014)] has shown that this could explain the emergence of a film-height-dependent surface tension close to the wall, with implications for the Frumkin–Derjaguin theory.     

Investigation of vapour--liquid nucleation properties for spherical and chain-like fluids by density functional theory

The excess Helmholtz free energy functional for nonpolar chain-like molecules is formulated in terms of a weighted density approximation （WDA） for short-range interactions and a Weaks Chandler Andersen （WCA） approximation and a Barker Henderson （BH） theory for long-range attraction. Within the framework of density functional theory （DFT）, vapour liquid interracial properties including density profile and surface tension, and vapour-liquid nucleation properties including density profile, work of formation and number of particles are investigated for spherical and chain- like molecules. The obtained vapour liquid surface tension and the number of particles in critical nucleus for Lennard- Jones （L J） fluids are consistent with the simulation results. The influences of supersaturation, temperature and chain length on vapour liquid nucleation properties are discussed.[第一段]