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

Contribution of Different Molecules and Moieties to the Surface Tension in Aqueous Surfactant Solutions. II: Role of the Size and Charge Sign of the Counterions

Understanding the role of the counterion species in surfactant solutions is a complicated task, made harder by the fact that, experimentally, it is not possible to vary independently bulk and surface quantities. Here, we perform molecular dynamics simulations at constant surface coverage of the liquid/vapor interface of lithium, sodium, potassium, rubidium, and cesium dodecyl sulfate aqueous solutions. We investigate the effect of counterion type and charge sign on the surface tension of the solution, analyzing the contribution of different species and moieties to the lateral pressure profile. The observed trends are qualitatively compatible with the Hofmeister series, with the notable exception of sodium. We point out a possible shortcoming of what is at the moment, in our experience, the most realistic nonpolarizable force field (CHARMM36) that includes the parametrization for the whole series of alkali counterions. In the artificial system where the counterion and surfactant charges are inverted in sign, the counterions become considerably harder. This charge inversion changes considerably the surface tension contributions of the counterions, surfactant headgroups, and water molecules, stressing the key role of the hardness of the counterions in this respect. However, the hydration free energy gain of the counterions, occurring upon charge inversion, is compensated by the concomitant free energy loss of the headgroups and water molecules, leading to a negligible change in the surface tension of the entire system.

Congener-dependent conformations of isolated rhamnolipids at the vacuum-water interface: A molecular dynamics simulation

HYPOTHESIS

Molecular dynamics simulation can be used to differentiate between the adsorption properties of rhamnolipid congeners at a vacuum-water interface.

EXPERIMENTS

Adsorption of five congeners with differing alkyl chains (two C10 chains, two C14 chains or mixed C14C10 and C10C14), number of rhamnose rings (mono- or di-) and carboxyl group charge (non-ionic or anionic) are simulated at the vacuum-water interface.

FINDINGS

All rhamnolipids adsorb in the interfacial region with rhamnose and carboxyl groups closer to the water phase, and alkyl chains closer to the vacuum phase, but with differing adsorbed conformations. Headgroups of uncharged congeners show two preferred conformations, closed and partially open. Di-rhamnolipid has a low proportion of closed conformation, due to the steric constraints of the second pyranose ring. Charged congeners show strong preference for closed headgroup conformations. For rhamnolipids with equal alkyl chains lengths (C10C10, C14C14) the distribution of alkyl chain tilt angles is similar for both. Where chain lengths are unequal (C14C10, C10C14) one chain has a greater tendency to tilt towards the water phase (>90°). The order parameter of the alkyl chains shows they are disordered at the interface. Together, these results show congener-dependent adsorbed conformation differences suggesting they will have differing surface-active properties at vacuum-water and oil-water interfaces.

Preferential Orientations and Anomalous Interfacial Tensions in Aqueous Solutions of Alcohols

Literature studies on interfacial tension versus temperature between normal alcohols and water show that it increases with temperature and exhibits a maximum value at a given temperature depending on the molecular weight of the alcohol. This very unusual behavior is supposedly accompanied by the formation of monolayers of alcohol molecules oriented preferentially at the interface, a structural issue not confirmed until now. We use molecular-based models for water and alcohols in combination with molecular dynamics simulations to provide physical insights, from a molecular perspective, into the structural and thermodynamic behavior at the liquid-liquid interfaces of aqueous solutions of alcohols.

Role of the Counterions in the Surface Tension of Aqueous Surfactant Solutions. A Computer Simulation Study of Alkali Dodecyl Sulfate Systems

We have investigated the surface tension contributions of the counterions, surfactant headgroups and tails, and water molecules in aqueous alkali dodecyl sulfate (DS) solutions close to the saturated surface concentration by analyzing the lateral pressure profile contribution of these components using molecular dynamics simulations. For this purpose, we have used the combination of two popular force fields, namely KBFF for the counterions and GROMOS96 for the surfactant, which are both parameterized for the SPC/E water model. Except for the system containing Na+ counterions, the surface tension of the surfactant solutions has turned out to be larger rather than smaller than that of neat water, showing a severe shortcoming of the combination of the two force fields. We have traced back this failure of the potential model combination to the unphysically strong attraction of the KBFF counterions, except for Na+, to the anionic head of the surfactants. Despite this failure of the model, we have observed a clear relation between the soft/hard character (in the sense of the Hofmeister series) and the surface tension contribution of the counterions, which, given the above limitations of the model, can only be regarded as an indicative result. We emphasize that the obtained results, although in a twisted way, clearly stress the crucial role the counterions of ionic surfactants play in determining the surface tension of the aqueous surfactant solutions.

Adsorption of Sodium Alkyl Sulfate Homologues at the Air/Solution Interface

Experimental results are presented on the adsorption of sodium alkyl sulfate homologues (nC = 8-14) at the air/solution interface. The adsorption isotherms calculated from equilibrium surface-tension vs concentration data and the critical micelle concentration change regularly with the length of the alkyl chain; the odd/even effect was not observed. The isotherms were analyzed using a model-independent approach. The analysis indicates that the total driving force of adsorption reaches a plateau value and becomes constant in the function of the adsorbed amount in the case of each homologue. With the use of different electrostatic models, it was demonstrated that this behavior is consistent with a saturation-type hydrophobic driving-force contribution, which can be interpreted by the development of a liquidlike alkane environment in the adsorbed layer above a "critical" adsorbed amount.

Adsorption of Octyl Cyanide at the Free Water Surface as Studied by Monte Carlo Simulation

Monte Carlo simulations of the adsorption layer of octyl cyanide have been performed on the canonical (N, V, T) ensemble at 300 K. The systems simulated cover the range of octyl cyanide surface densities from 0.27 to 7.83 mumol/m2. The surface density value at which the saturation of the adsorption layer occurs is estimated to be 1.7 mumol/m2. At low surface densities, the main driving force of the adsorption is found to be the formation of hydrogen bonds between the water and octyl cyanide molecules, whereas at higher surface concentrations, the dipole-dipole attraction between the neighboring adsorbed octyl cyanide molecules becomes more important. At low surface concentrations, the water-octyl cyanide hydrogen bonds prefer tilted alignments relative to the interface; however, in the case of the saturated adsorption layer, the number of such hydrogen bonds is maximized, leading to the preference of these bonds for the orientation perpendicular to the interface. Contrary to nonionic surfactants of multiple hydrogen bonding abilities (e.g., 1-octanol, C8E3), the increasing surface concentration of octyl cyanide was not found to lead to considerable competition of the molecules for positions of optimal arrangement. As a consequence, the energy and geometry of the water-octyl cyanide hydrogen bonds are found to be insensitive to the octyl cyanide surface concentration.

Counterion and Surface Density Dependence of the Adsorption Layer of Ionic Surfactants at the Vapor−Aqueous Solution Interface: A Computer Simulation Study

To test the validity of currently used adsorption theories and understand the origin of the lack of their ability of adequately describing existing surface tension measurement data, we have performed a series of molecular dynamics simulations of the adsorption layer of alkali decyl sulfate at the vapor/aqueous solution interface. The simulations have been performed with five different cations (i.e., Li+, Na+, K+, Rb+, and Cs+) at two different surface concentrations (i.e., 2 micromol/m2 and 4 micromol/m2). The obtained results clearly show that the thickness of the outer Helmholtz plate, a key quantity of the various adsorption theories, depends on two parameters, that is, the size of the cations and the surface density of the anionic surfactant. Namely, with increasing surface concentration, the electrostatic attraction between the two, oppositely charged, layers becomes stronger, leading to a considerable shrinking of the outer Helmholtz plate. Furthermore, this layer is found to be thicker in the presence of larger cations. The former effect could be important in understanding the anomalous shape of the adsorption isotherms of alkali alkyl sulfate surfactants, while the second effect seems to be essential in explaining the cation specificity of these isotherms.

Determination of the Adsorption Isotherm of Methanol on the Surface of Ice. An Experimental and Grand Canonical Monte Carlo Simulation Study

The adsorption isotherm of methanol on ice at 200 K has been determined both experimentally and by using the Grand Canonical Monte Carlo computer simulation method. The experimental and simulated isotherms agree well with each other; their deviations can be explained by a small (about 5 K) temperature shift in the simulation data and, possibly, by the non-ideality of the ice surface in the experimental situation. The analysis of the results has revealed that the saturated adsorption layer is monomolecular. At low surface coverage, the adsorption is driven by the methanol-ice interaction; however, at full coverage, methanol-methanol interactions become equally important. Under these conditions, about half of the adsorbed methanol molecules have one hydrogen-bonded water neighbor, and the other half have two hydrogen-bonded water neighbors. The vast majority of the methanols have a hydrogen-bonded methanol neighbor, as well.

Molecular Dynamics Study of 2-Nitrophenyl Octyl Ether and Nitrobenzene

The pure organic liquids nitrobenzene (NB) and 2-nitrophenyl octyl ether (NPOE) have been studied by means of molecular dynamics simulations. Both solvents are extremely important in various interfacial processes, mainly connected with ion transfer taking place across the interface with water. Thermodynamic (mass density, enthalpy of vaporization, isothermal compressibility, dipole moment) and dynamic (viscosities and self-diffusion coefficients) properties of both liquids have been calculated and are in very good agreement with the experimental data. In the case of NB, several potentials have been tested and the obtained results compared and discussed. In most cases, the OPLS all-atom potential gives results that are in better agreement with available experimental values. Atomic radial distribution functions, dihedral and angle distributions, as well as dipole-orientation correlation functions are used to probe the structure and interactions of the bulk molecules of both organic solvents. These were seen to be very similar in terms of structure and thermodynamics, but quite distinct in terms of dynamic behavior, with NPOE showing a much slower dynamic response than NB. A simulation study of the simple Cl- and K+ ions dissolved in both solvents has been also undertaken, revealing details about the diffusion and solvation mechanisms of these ions. It was found that in both liquids the positive potassium ion is solvated by the negative end of the molecular dipole, whereas the negative chloride ion is solvated by the positive end of the dipole.

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.

Structure of the nonionic surfactant triethoxy monooctylether C8E3 adsorbed at the free water surface, as seen from surface tension measurements and Monte Carlo simulations

The adsorption layer of the nonionic surfactant triethoxy monooctylether C8E3 has been investigated at the free water surface by means of both experimental and computer simulation methods. The surface tension of the aqueous solution of C8E3 has been measured by pendant drop shape analysis in the entire concentration range in which C8E3 is soluble in water. The data obtained from these measurements are used to derive the adsorption isotherm. The critical micellar concentration and the surface excess concentration of the saturated adsorption layer are found to be 7.48 mM and 4.03 micromol/m2, respectively, the latter value corresponding to the average area per molecule of 41 A2. In order to analyze the molecular level structure of the unsaturated adsorption layer, Monte Carlo simulations have been performed at four different surface concentration values, i.e., 0.68, 1.36, 2.04, and 2.72 micromol/m2, respectively. It has been found that the water surface is already almost fully covered at the lowest surface density value investigated, and the adsorbed molecules show a strong preference for lying parallel with the interface in elongated conformations. No sign of the penetration of the hydrophilic triethoxy headgroups into the aqueous phase to any extent has been observed. With increasing surface densities the preferential orientation of the apolar octyl tails gradually turns from lying parallel with the interface to pointing toward the vapor phase by their CH3 end, whereas the conformation of the adsorbed molecules becomes gradually less elongated. Both of these changes lead to the increase of the number of C8E3 molecules being in a direct contact (i.e., forming hydrogen bonds) with water. However, the increasing number of the C8E3 molecules hydrogen bonded to water is found to be accompanied by the weakening of this binding, i.e., the decrease of both the number of hydrogen bonds a bound C8E3 molecule forms with water and the magnitude of the average binding energy of the adsorbed C8E3 molecules.

Orientational order of the water molecules at the vicinity of the water–benzene interface in a broad range of thermodynamic states, as seen from Monte Carlo simulations

Monte Carlo simulation of the water/benzene liquid-liquid interfacial system has been performed at six different thermodynamic state points, ranging from ambient conditions up to the vicinity of the critical point of water. The system has been found to consist of two immiscible liquid phases at every state point studied. The orientational preferences of the interfacial water molecules have been analysed in detail using the simulated configurations. The results obtained at ambient conditions are in agreement with previous results on various different water/apolar interfaces. Thus, interfacial water molecules have been found to have dual orientational preferences: the molecules located nearest to the organic phase prefer to stay perpendicular to the interface, pointing flatly toward the apolar phase by their dipole vectors, whereas the waters located somewhat farther from the organic phase prefer the parallel alignment with the interface. The observed orientational preferences are found to be rather stable with changing thermodynamic conditions: although the increase of the temperature has led, due to the increasing thermal motion of the molecules, to a gradual weakening of the orientational preferences, both preferences are found to exist up to at least 450 K, and found to be completely washed out at 575 K only. The pressure has not been found to influence the orientation of the water molecules noticeably.