Structure of the liquid-vapor interface of water-methanol mixtures as seen from Monte Carlo simulations

Monte Carlo simulation of the vapor-liquid interface of water-methanol mixtures of different compositions, ranging from pure water to pure methanol, have been performed on the canonical (N, V, T) ensemble at 298 K. The analysis of the systems simulated has revealed that the interface is characterized by a double layer structure: methanol is strongly adsorbed at the vapor side of the interface, whereas this adsorption layer is followed at its liquid side by a depletion layer of methanol of lower concentration than in the bulk liquid phase of the system. The dominant feature of the interface has been found to be the adsorption layer in systems of methanol mole fractions below 0.2, and the depletion layer in systems of methanol mole fractions between 0.25 and 0.5. The orientation of the molecules located at the depletion layer is found to be already uncorrelated with the interface, whereas the methanol molecules of the adsorption layer prefer to align perpendicular to the interface, pointing straight toward the vapor phase by their methyl group. Although both the preference of the molecular plane for a perpendicular alignment with the interface and the preference of the methyl group for pointing straight to the vapor phase are found to be rather weak, the preference of the methyl group for pointing as straight toward the vapor phase as possible within the constraint imposed by the orientation of the molecular plane is found to be fairly strong. One of the two preferred orientations of the interfacial water molecules present in the neat system is found to disappear in the presence of methanol, because methanol molecules aligned in their preferred orientation can replace these water molecules in the hydrogen-bonding pattern of the interface.     

Computer simulation investigation of the water-benzene interface in a broad range of thermodynamic States from ambient to supercritical conditions

The dependence of the properties of the water-benzene system on the thermodynamic conditions in a broad range of temperatures and pressures has been investigated by computer simulation methods. For this purpose, Monte Carlo simulations have been performed at 23 different thermodynamic states, ranging from ambient to supercritical conditions. The density profiles of the water and benzene molecules have been determined at each of the thermodynamic states investigated. Information on the dependence of the mutual solubility of the two components in each other as well as of the width of the interface on the temperature and pressure has been extracted from these profiles. The width of the interface has been found to increase with increasing temperature up to a certain point, where it diverges. The temperature of this divergence corresponds to the mixing of the two phases. The determination of the critical mixing temperature at various pressures allowed us to estimate the upper critical curve, separating the two-phase and one-phase liquid systems, of the phase diagram of the simulated water-benzene system. In analyzing the preferential orientation of the interfacial molecules relative to the interface, it has been found that the main orientational preference of the benzene molecules is to lie parallel with the plane of the interface, and the water molecules penetrated deepest into the benzene phase prefer to stay perpendicular to the interface, pointing by one of their O-H bonds almost straight toward the benzene phase, whereas the waters located at the aqueous side of the interface are preferentially aligned parallel with the interfacial plane. Although the strength of the observed orientational preferences decreases rapidly with increasing temperature, the preferred orientations themselves are found to be independent of the thermodynamic conditions. Remains of the orientational preferences of the molecules are found to be present up to temperatures as high as 650 K. The analysis of the relative orientation of the neighboring water-benzene pairs has revealed that the radius of the first hydration shell of the benzene molecules is independent of the thermodynamic conditions, even if the system consists of one single phase. It has been found that the nearest water neighbors of the benzene molecules are preferentially located above and below the benzene ring, whereas more distant water neighbors, belonging still to the first hydration shell, prefer to stay within the plane of the benzene molecule. In the two-phase systems the dipole vector of the nearest waters has been found to be preferentially perpendicular to the vector pointing from the center of the benzene molecule to the water O atom.     

Molecular level structure of the liquid/liquid interface. Molecular dynamics simulation and ITIM analysis of the water-CCl4 system

The molecular level properties of the liquid/liquid interface between water and CCl(4) are analysed in detail on the basis of molecular dynamics computer simulation. This analysis requires a full list of the molecules that are right at the interface in both phases. Such a list can be provided by the novel method for identifying truly interfacial molecules (ITIM). The full list of the truly interfacial molecules various properties (e.g., width, molecular level roughness) of the interface can be meaningfully analysed. The residence time of the molecules at the interface, the percolation of the water molecules at the interfacial layer as well as in the second layer beneath the surface, the preferred orientations of the interfacial water molecules and the dependence of these orientational preferences on the local curvature of the interface are also analysed and discussed in detail.     

Orientation of the 3-methylpyridine molecules at the liquid–vapour interface of their aqueous solution

The orientation of the 3-methylpyridine (3MPy) molecules at the liquid–vapour interface of their 3 mol% solution is analysed in detail at 298 K on the basis of a Monte Carlo simulation. The 3MPy molecules are found to have a strong tendency for being adsorbed at the interface, exhibiting a dual orientational preference. At the vapour side of the interface they are preferentially aligned perpendicular to the interface, in such a way that the para carbon atom of the pyridine ring points straight to the vapour phase, whereas the C–Me bond declines by about 60° from the interface normal. On the other hand, at the liquid side of the interface the preferred orientation of the 3MPy molecules is close to the parallel alignment with the plane of the interface.     

Adsorption of 1-octanol at the free water surface as studied by Monte Carlo simulation

The adsorption of 1-octanol at the free water surface has been investigated by Monte Carlo computer simulation. Six different systems, built up by an aqueous and a vapor phase, the latter also containing various number of octanol molecules, have been simulated. The number of the octanol molecules has been chosen in such a way that the octanol surface density varies in a broad range, between 0.27 and 7.83 micromol/m(2) in the six systems simulated. For reference, the interfacial system containing bulk liquid octanol in the apolar phase has also been simulated. The results have shown that the formation of hydrogen bonds between the interfacial water and adsorbed octanol molecules is of key importance in determining the properties of the adsorbed layer. At low octanol surface concentration values all the octanol molecules are strongly (i.e., by hydrogen bonds) bound to the aqueous phase, whereas their interaction with each other is negligibly small. Hence, they are preferentially oriented in such a way that their own binding energy (and thus their own free energy) is minimized. In this preferred orientation the O-H bond of the octanol molecule points flatly toward the aqueous phase, declining by about 30 degrees from the interfacial plane, irrespectively from whether the octanol molecule is the H-donor or the H-acceptor partner in the hydrogen bond. Hence, in its preferred orientation the octanol molecule can form at least two low energy hydrogen bonds with water: one as a H-donor and another one as a H-acceptor. Moreover, the preferred orientation of the hydrogen bonded water partners is close to one of the two preferred interfacial water alignments, in which the plane of the water molecule is parallel with the interface. When increasing the octanol surface density, the water surface gets saturated with hydrogen bonded octanols, and hence any further octanol molecule can just simply condense to the layer of the adsorbed octanols. The surface density value at which this saturation occurs is estimated to be about 1.7 micromol/m(2). Above this surface density value the hydrogen bonded octanols and their water partners are oriented in such a way that the number of the water-octanol hydrogen bonds is maximized. Hence, the preferred alignment of the O...O axes of these hydrogen bonds is perpendicular to the interface. This orientation is far from the optimal alignment of the individual octanol molecules, which is also reflected in the observed fact that, unlike in the case of many other adsorbents, the average molecular binding energy of the adsorbed octanol molecules increases (i.e., becomes less negative) with increasing octanol surface density.     

Orientational motions of vibrational chromophores in molecules at the air/water interface with time-resolved sum frequency generation

The first time-resolved experiments in which interfacial molecules are pumped to excited electronic states and probed by vibrational sum frequency generation (SFG) are reported. This method was used to measure the out-of-plane rotation dynamics, i.e. time dependent changes in the polar angle, of a vibrational chromophore of an interfacial molecule. The chromophore is the carbonyl group, the rotation observed is that of the -C=O bond axis, with respect to the interfacial normal, and the interfacial molecule is coumarin 314 (C314) at the air/water interface. The orientational relaxation time was found to be 220+/-20 ps, which is much faster than the orientational relaxation time of the permanent dipole moment axis of C314 at the same interface, as obtained from pump-second harmonic probe experiments. Possible effects on the rotation of the -C=O bond axis due to the carbonyl group hydrogen bonding with interfacial water are discussed. From the measured equilibrium orientation of the permanent dipole moment axis and the carbonyl axis, and knowledge of their relative orientation in the molecule, the absolute orientation of C314 at the air/water interface is obtained.     

A systematic molecular simulation study of ionic liquid surfaces using intrinsic analysis methods

In this paper, we apply novel intrinsic analysis methods, coupled with bivariate orientation analysis, to obtain a detailed picture of the molecular-level structure of ionic liquid surfaces. We observe pronounced layering at the interface, alternating non-polar with ionic regions. The outermost regions of the surface are populated by alkyl chains, which are followed by a dense and tightly packed layer formed of oppositely charged ionic moieties. We then systematically change the cation chain length, the anion size, the temperature and the molecular model, to examine the effect of each of these parameters on the interfacial structure. Increasing the cation chain length promotes orientations in which the chain is pointing into the vapor, thus increasing the coverage of the surface with alkyl groups. Larger anions promote a disruption of the dense ionic layer, increasing the orientational freedom of cations and increasing the amount of free space. The temperature had a relatively small effect on the surface structure, while the effect of the choice of molecular model was clearly significant, particularly on the orientational preferences at the interface. Our study demonstrates the usefulness of molecular simulation methods in the design of ionic liquids to suit particular applications.     

Orientational correlations in liquid acetone and dimethyl sulfoxide: a comparative study

The structure of acetone and dimethyl sulfoxide in the liquid state is investigated using a combination of neutron diffraction measurements and empirical potential structure refinement (EPSR) modeling. By extracting the orientational correlations from the EPSR model, the alignment of dipoles in both fluids is identified. At short distances the dipoles or neighboring molecules are found to be in antiparallel configurations, but further out the molecules tend to be aligned predominately as head to tail in the manner of dipolar ordering. The distribution of these orientations in space around a central molecule is strongly influenced by the underlying symmetry of the central molecule. In both liquids there is evidence for weak methyl hydrogen to oxygen intermolecular contacts, though these probably do not constitute hydrogen bonds as such.     

Absolute orientation of molecules at interfaces

A method to determine the absolute orientation of molecules at liquid interfaces by sum frequency generation (SFG) is reported. It is based on measurements of the orientations of two nonparallel vibrationally active chromophores in the molecule of interest combined with a rotation matrix formulation to obtain the absolute molecular orientation. We chose m-tolunitrile, a planar molecule adsorbed to the air/water interface, as a proof-of-method experiment. Quantitative analysis of different polarization sum frequency intensities facilitate unique peak assignments of the methyl and nitrile groups of m-tolunitrile. The SFG analysis of the measurement yields a nitrile group tilting at 53 degrees to the surface normal, and the C3 axis of the methyl group is almost upright at 23 degrees with respect to the surface normal. Using a rotation matrix formulation, we found that the angle between the surface plane and the m-tolunitrile molecular plane is 70 degrees.     

Imaging the binding ability of proteins immobilized on surfaces with different orientations by using liquid crystals

We report an investigation of the binding ability of a protein immobilized on surfaces with different orientations but in identical interfacial microenvironments. The surfaces present mixed self-assembled monolayers (SAMs) of 11-[19-carboxymethylhexa(ethylene glycol)]undecyl-1-thiol, 1, and 11-tetra(ethylene glycol) undecyl-1-thiol, 2. Whereas 2 is used to define an interfacial microenvironment that prevents nonspecific adsorption of proteins, 1 was activated by two different schemes to immobilize ribonuclease A (RNase A) in either a preferred orientation or random orientations. The binding of the ribonuclease inhibitor protein (RI) to RNase A on these surfaces was characterized by using ellipsometry and the orientational behavior of liquid crystals. Ellipsometric measurements indicate identical extents of immobilization of RNase A via the two schemes. Following incubation of both surfaces with RI, however, ellipsometric measurements indicate a 4-fold higher binding ability of the RNase A immobilized with a preferred orientation over RNase A immobilized with a random orientation. The higher binding ability of the oriented RNase A over the randomly oriented RNase A was also apparent in the orientational behavior of nematic liquid crystals of 4-cyano-4'-pentylcyanobiphenyl (5CB) overlayed on these surfaces. These results demonstrate that the orientations of proteins covalently immobilized in controlled interfacial microenvironments can influence the binding activities of the immobilized proteins. Results reported in this article also demonstrate that the orientational states of proteins immobilized at surfaces can be distinguished by examining the optical appearances of liquid crystals.     

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.     

The extent of molecular orientation at liquid/vapor interface of pyridine and its alkyl derivatives by molecular dynamics simulation

In this paper, molecular dynamics simulation was performed to investigate the liquid∕vapor interfacial structure of neat polar liquids. Large-scale ensembles of liquid pyridine and its alkyl derivatives, 4-methylpyridine and 4-ethylpyridine, were simulated by classical molecular dynamics at 298 K. For the liquid system of low polarity, the surface density profile of the atoms meet exactly at the middle of interfacial region, and atoms of hydrophobic nature can be hardly discriminated from hydrophilic ones in either vapor or liquid sides. For a liquid system of high polarity, the density profile of atoms with different nature is highly discriminated all over the interfacial region, and as the polarity increases, a dense region of atomic density is clearly developed in the subsurface region. The recognized bivariate method was also used to study the molecular orientational distribution quantitatively. Orientational analysis of the three liquid systems indicates that the pyridine ring plane in the outmost surface tends to be vertical. Its tendency in the innermost interfacial region is parallel. The orientational states available to 4-ethylpyridine and pyridine are discriminated by predicting the possibility of a bisector-wise tumbling for the ring plane in pyridine and a side-wise tumbling in 4-ethylpyridine. The orientational distribution maps explain the trend of experimental surface tension and surface entropy. As the dipole moment of these liquids increases with the alkyl chain length, the surface structural profile changes from a regular definite one to a surface of complex atomic structure involving a dense phase near the interface. The development of dense region in alkyl derivatives is the result of segregation of molecules due to the alkyl group, which is captured and discriminated by molecular dynamics simulation even when the length of a short alkyl chain is increased by one carbon atom.     

Hydrogen bonded structure and dynamics of liquid-vapor interface of water-ammonia mixture: an ab initio molecular dynamics study

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations.     

Determination of structure and energetics for gibbs surface adsorption layers of binary liquid mixture 1. Acetone + water

The orientation, structure, and energetics of the vapor/acetone-water interface are studied with sum frequency generation vibrational spectroscopy (SFG-VS). We used the polarization null angle (PNA) method in SFG-VS to accurately determine the interfacial acetone molecule orientation, and we found that the acetone molecule has its C=O group pointing into bulk phase, one CH3 group pointing up from the bulk, and the other CH3 group pointing into the bulk phase. This well-ordered interface layer induces an antiparallel structure in the second layer through dimer formation from either dipolar or hydrogen bond interactions. With a double-layer adsorption model (DAM) and Langmuir isotherm, the adsorption free energies for the first and second layer are determined as deltaG degrees (ads,1) = - 1.9 +/- 0.2 kcal /mol and deltaG degrees (ads,2) = - 0.9 +/- 0.2 kcal /mol, respectively. Since deltaG degrees (ads,1) is much larger than the thermal energy kT = 0.59 kcal /mol, and deltaG degrees (ads,2) is close to kT, the second layer has to be less ordered. Without either strong dipolar or hydrogen bonding interactions between the second and the third layer, the third layer should be randomly thermalized as in the bulk liquid. Therefore, the thickness of the interface is not more than two layers thick. These results are consistent with previous MD simulations for the vapor/pure acetone interface, and undoubtedly provide direct microscopic structural evidences and new insight for the understanding of liquid and liquid mixture interfaces. The experimental techniques and quantitative analysis methodology used for detailed measurement of the liquid mixture interfaces in this report can also be applied to liquid interfaces, as well as other molecular interfaces in general.     

A novel mesophase formed by top-shaped molecules in the bulk and unsupported thin films: a molecular dynamics study

We have used molecular dynamics simulations to investigate the ordering of top-shaped molecules in bulk phases and in unsupported thin films. Each rigid anisotropic molecule was composed of 11 Lennard-Jones interaction centers (beads). In an attempt to enhance the nematic stability in preference to smectic, the three central beads were assigned a larger Lennard-Jones diameter than the tail beads, giving the molecule a shape resembling a top. The molecular model was found to exhibit an unusual bulk mesophase with long-range orientational order and with molecular center-of-mass positions arranged in parallel interdigitated layers, with layer spacing smaller than half the length of the long axis of a molecule. However, despite the toplike molecular shape, no nematic phase was observed in the pressure range studied. Unsupported films of the isotropic liquid were cooled in order to locate a triple point between the novel mesophase, vapor, and isotropic liquid. At temperatures slightly above the triple point, enhanced surface ordering of molecules was found to occur in the unsupported film. At temperatures slightly below the triple point, the preferred molecular alignment in the unsupported film was parallel to the interface, in violation of arguments that have been proposed based on the relative enthalpies of various cleavage planes for close-packed structures.     

Molecular dynamics simulations of ionic liquid-vapour interfaces: effect of cation symmetry on structure at the interface

The influence of alkyl chain symmetry of the imidazolium cation on the structure and properties of the ionic liquid-vapour interface has been addressed through molecular dynamics simulations. The anion chosen is bis(trifluoromethylsulfonyl)imide (NTf(2)). Profiles of number densities, orientation of cations, charge density, electrostatic potential, and surface tension have been obtained. At the interface, both cations and anions were present, and the alkyl chains of the former preferred to orient out into the vapour phase. A large fraction of cations preferred to be oriented with their ring-normal parallel to the surface and alkyl chains perpendicular to it. These orientational preferences are reduced in ionic liquids with symmetric cations. Although the charge densities at the interface were largely negative, an additional small positive charge density has been observed for systems with longer alkyl chains. The electrostatic potential difference developed between the liquid and the vapour phases were positive and decreased with increasing length of the alkyl group. The calculated surface tension of the liquids also decreased with increasing alkyl chain length, in agreement with experiment. The surface tension of an ionic liquid with symmetric cation was marginally higher than that of one with an asymmetric, isomeric cation.     

Classical density functional theory of orientational order at interfaces: application to water

A classical density functional formalism has been developed to predict the position-orientation number density of structured fluids. It is applied to the liquid-vapor interface of pure water, where it consists of a classical term, a gradient correction, and an anisotropic term that yields order through density gradients. The model is calibrated to predict that water molecules have their dipole moments almost parallel to a planar interface, while the molecular plane is parallel to it on the liquid side and perpendicular to it on the vapor side. For a planar interface, the surface tension obtained is twice its experimental value, while the surface potential is in qualitative agreement with that calculated by others. The model is also used to predict the orientation of water molecules near the surface of droplets, as well as the dependence of equilibrium vapor pressure around them on their size.     

Structural characterization of interfacial n-octanol and 3-octanol using molecular dynamic simulations

Structurally isomeric octanol interfacial systems, water/vapor, 3-octanol/vapor, n-octanol/vapor, 3-octanol/water, and n-octanol/water are investigated at 298 K using molecular dynamics simulation techniques. The present study is intended to investigate strongly associated liquid/liquid interfaces and probe the atomistic structure of these interfaces. The octanol and water molecules were initially placed randomly into a box and were equilibrated using constant pressure techniques to minimize bias within the initial conditions as well as to fully sample the structural conformations of the interface. An interface formed via phase separation during equilibration and resulted in a slab geometry with a molecularly sharp interface. However, some water molecules remained within the octanol phase with a mole fraction of 0.12 after equilibration. The resulting "wet" octanol interfaces were analyzed using density profiles and orientational order parameters. Our results support the hypothesis of an ordered interface only 1 or 2 molecular layers deep before bulk properties are reached for both the 3-octanol and water systems. However, in contrast to most other interfacial systems studied by molecular dynamics simulations, the n-octanol interface extends for several molecular layers. The octanol hydroxyl groups form a hydrogen-bonding network with water which orders the surface molecules toward a preferred direction and produces a hydrophilic/hydrophobic layering. The ordered n-octanol produces an oscillating low-high density of oxygen atoms out of phase with a high-low density of carbon atoms, consistent with an oscillating dielectric. In contrast, the isomeric 3-octanol has only a single carbon-rich layer directly proximal to the interface, which is a result of the different molecular topology. Both 3-octanol and n-octanol roughen the water interface with respect to the water/vapor interface. The "wet" octanol phases, in the octanol/water systems reach bulk properties in a shorter distance than the "dry" octanol/vapor interfaces.     

Structure and dynamics of surfactant and hydrocarbon aggregates on graphite: a molecular dynamics simulation study

We have examined the structure and dynamics of sodium dodecyl sulfate (SDS) and dodecane (C12) molecular aggregates at varying surface coverages on the basal plane of graphite via classical molecular dynamics simulations. Our results suggest that graphite-hydrocarbon chain interactions favor specific molecular orientations at the single-molecule level via alignment of the tail along the crystallographic directions. This orientational bias is reduced greatly upon increasing the surface coverage for both molecules due to intermolecular interactions, leading to very weak bias at intermediate surface coverages. Interestingly, for complete monolayers, we find a re-emergent orientational bias. Furthermore, by comparing the SDS behavior with C12, we demonstrate that the charged head group plays a key role in the aggregate structures: SDS molecules display a tendency to form linear file-like aggregates while C12 forms tightly bound planar ones. The observed orientational bias for SDS molecules is in agreement with experimental observations of hemimicelle orientation and provides support for the belief that an initial oriented layer governs the orientation of hemimicellar aggregates.