Enrique Chacón

The intrinsic structure of the water surface

An operational procedure to obtain the intrinsic structure of liquid surfaces is applied here to a molecular dynamics simulation of water, with a model of point charges for the molecular interactions. The method, which had been recently proposed and used for simple fluids, is successfully extended to a molecular liquid with the complex bond structure of water. The elimination of the capillary wave fluctuations, in the intrinsic density and orientation profiles, gives a new overall view of the water surface, at the sharpest molecular level, and without the size-dependent broadening observed in the mean profiles. The molecules belonging to the outer liquid layer are clearly identified, and we find that only these molecules exhibit a clear preferential orientation to lie flat on the surface. Moreover, there is a strong correlation between the dipolar structure and the local curvatures of the intrinsic surface, so that at the extrusions of the intrinsic surface the molecular dipoles point preferentially toward the vapor side of the interface. Finally, we have found an intrinsic density layering structure, although the inner structure is strongly damped beyond the second layer.

A model for membranes, vesicles and micelles in amphiphilic systems

We present a microscopic model for the aggregates of amphiphilic molecules, based on a simple density functional approximation for the free energy. The different molecular aggregates are described as self-structured density distributions at the relative minima of the grand potential energy. We search for these structures with planar and spherical geometries, and obtain the phase diagram for bilayer membranes, and the curvature energies for vesicles and different types of micelles. The study of a global phase diagram, to get the density of micelles and isolated amphiphilic molecules, at equilibrium with free membranes, requires me link between two description levels of micelles: as self-structured density distributions, or as molecular clusters in the solution of amphiphilic molecules in water. This is done with the help of a simple harmonic model which provides an appropriate choice of the configurational unit cell for micelles.

Molecular dynamics investigation of the intrinsic structure of water–fluid interfaces via the intrinsic sampling method

Capillary wave fluctuations smooth out the structure of fluid interfaces, making difficult the detailed analysis of the interfacial structure. Most computer simulation investigations performed to date have focused on the computation of average density profiles, ignoring the characterization of the intrinsic structure of the interface. Recent theoretical developments have reversed this situation, making possible the detailed investigation of the interfacial intrinsic structure at an unprecedented level of detail. In this article we investigate via molecular dynamics simulations the intrinsic structure of water-alkane (hexane and dodecane) interfaces. The implementation of the recently introduced, intrinsic sampling method to compute the intrinsic surface of water-fluid interfaces is discussed. We provide quantitative molecular information on the structure, corrugation, and stiffness of the liquid surfaces. The intrinsic structure of water at alkane interfaces is shown to be insensitive to the alkane-chain length, and can be very accurately described by the intrinsic structure of the water free surface.

Characterization of the intrinsic density profiles for liquid surfaces

This paper presents recent advances in the characterization of the intrinsic structures in computer simulations of liquid surfaces. The use of operational definitions for the intrinsic surface, associated with each molecular configuration of a liquid slab, gives direct access to the intrinsic profile and to the wavevector dependent surface tension. However, the characteristics of these functions depend on the definition used for the intrinsic surface. We discuss the pathologies associated with a local Gibbs dividing surface definition, and consider the alternative definition of a minimal area surface, going though a set of surface pivots, self-consistently chosen to represent the first liquid layer.

Low melting temperature and liquid surface layering for pair potential models

We have recently proposed [Phys. Rev. Lett. 87, 166101 (2001)] that any isotropic fluid should exhibit surface layering at its liquid–vapor interface above the triple temperature provided that the system has a low triple temperature. In this article we present an extensive study of systems with different isotropic pair interactions, some of which present a very low triple temperature. We have confirmed that surface layering is a general characteristic of very cold liquids, independent of the specific shape of the potential, and that only pair potentials presenting a low triple-point temperature do exhibit surface oscillations; in other cases layering is preempted by solidification. Finally, we study the damping of surface oscillations due to capillary waves and conclude that for any model pair potential the temperature threshold below which layering would be observed for the typical experimental transverse sampling sizes is 15% of the critical temperature.

Atomically resolved three-dimensional structures of electrolyte aqueous solutions near a solid surface

Interfacial liquid layers play a central role in a variety of phenomena ranging from friction to molecular recognition. Liquids near a solid surface form an interfacial layer where the molecular structure is different from that of the bulk. Here we report atomic resolution three-dimensional images of electrolyte solutions near a mica surface that demonstrate the existence of three types of interfacial structures. At low concentrations (0.01–1 M), cations are adsorbed onto the mica. The cation layer is topped by a few hydration layers. At higher concentrations, the interfacial layer extends several nanometres into the liquid. It involves the alternation of cation and anion planes. Fluid Density Functional calculations show that water molecules are a critical factor for stabilizing the structure of the interfacial layer. The interfacial layer stabilizes a crystal-like structure compatible with liquid-like ion and solvent mobilities. At saturation, some ions precipitate and small crystals are formed on the mica.

The intrinsic interfacial structure of ionic surfactant monolayers at water-oil and water-vapour interfaces

Using computer simulations, we investigate the interfacial structure of sodium dodecyl sulphate (SDS) monolayers adsorbed at the water surface and water–oil interfaces. Using an algorithm that removes the averaging effect of the capillary waves, we obtain a detailed view of the solvation structure of water around the monolayer. We investigate surface concentrations between 45 and 33 Å2 per surfactant, which are near experimental conditions corresponding to the critical micellar concentration and the formation of Newton black films. The surfactants induce a layering structure in water, which disappears at approximately 1 nm from the monolayer plane. The water molecules exhibit a preferred orientation with the dipoles pointing towards the monolayer. The orientational order decays slowly, but it does not influence the hydrogen bond structure of water, which is significantly disrupted in the interfacial region only. These structural changes are qualitatively the same in SDS–water and oil–SDS–water interfaces. In the latter case, we find a small degree of penetration of oil in the monolayer (between 0.2 and 0.25 molecules per SDS). This small penetration has a measurable effect on the monolayer, which increases its thickness by approximately 10 per cent. The bending modulus of the SDS monolayers is of the order of the thermal energy, kBT.

Thermal fluctuations and bending rigidity of bilayer membranes

We present a new scheme to estimate the elastic properties of biological membranes in computer simulations. The method analyzes the thermal fluctuations in terms of a coupled undulatory mode, which disentangle the mixing of the mesoscopic undulations and the high-q protrusions. This approach makes possible the accurate estimation of the bending modulus both for membranes under stress and in tensionless conditions; it also extends the applicability of the fluctuation analysis to the small membrane areas normally used in atomistic simulations. Also we clarify the difference between the surface tension imposed in simulations through a pressure coupling barostat, and the surface tension that can be extracted from the analysis of the low wave vector dependence of the coupled undulatory fluctuation spectrum. The physical analysis of the peristaltic mode is also refined, by separating the bulk and protrusions contributions. We illustrate the procedure by analyzing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayers. The bending moduli obtained from our analysis, shows good agreement with available experiments.

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.

A functional perturbation theory for nonuniform molecular fluids. Effect of a weak dipole in a liquid–vapor interphase

A perturbation theory for nonuniform molecular fluids using the functional formalism is developed. Approximate expressions are written down for axial molecules using the spherical harmonic expansion. Under certain approximations Tarazona–Navascues and Haile et al. perturbation theories are recovered. Numerical calculation for the molecular orientation in liquid–vapor interphase of an axial fluid with a weak dipole is included.