Monte Carlo simulation of mixed lennard-jones nonionic surfactant adsorption at the liquid/vapor interface

New Monte Carlo simulations are presented for nonionic surfactant adsorption at the liquid/vapor interface of a monatomic solvent specifically investigating the roles of tail attraction and binary mixtures of different tail lengths. Surfactant molecules consist of an amphiphilic chain with a solvophilic head and a solvophobic tail. All molecules in the system, solvent and surfactant, are characterized by the Lennard-Jones (LJ) potential. Adjacent atoms along the surfactant chain are connected by finitely extensible harmonic springs. Solvent molecules move via the Metropolis random-walk algorithm, whereas surfactant molecules move according to the continuum configurational bias Monte Carlo (CBMC) method. We generate thermodynamic adsorption and surface-tension isotherms and compare results quantitatively to single-surfactant adsorption (Langmuir, 2007, 23, 1835). Surfactant tail groups with attractive interaction lead to cooperative adsorption at high surface coverage and higher maximum adsorption at the interface than those without. Moreover, adsorption and surface-tension isotherms with and without tail attraction are identical at low concentrations, deviating only near maximum coverage. Simulated binary mixtures of surfactants with differing lengths give intermediate behavior between that of the corresponding single-surfactant adsorption and surface-tension isotherms both with and without tail attraction. We successfully predict simulated mixture results with the thermodynamically consistent ideal adsorbed solution (IAS) theory for binary mixtures of unequal-sized surfactants using only the simulations from the single surfactants. Ultimately, we establish that a coarse-grained LJ surfactant system is useful for understanding actual surfactant systems when tail attraction is important and for unequal-sized mixtures of amphiphiles.     

Ellipsometry studies of nonionic surfactant adsorption at the oil--water interface

In the presented study we have developed and implemented a methodology for ellipsometry measurements at liquid interfaces that makes it possible to determine the amount adsorbed without assumptions of refractive index or thickness of the adsorbed layer. It was demonstrated that this is possible by combined measurements from different aqueous phases, H(2)O and D(2)O, which were shown to have sufficiently different refractive indices. The methodology was tested by studying adsorption of two types of nonionic poly(ethylene glycol) alkyl ether surfactants, C(n)H(2)(n)(+1)(OC(2)H(4))(m)OH or C(n)E(m) at the decane--aqueous interface, where C(12)E(5) was adsorbed from the oil phase and C(18)E(50) from the aqueous phase. The observed plateau values of the adsorbed amounts were 1.38 and 0.93 mg/m(2) for C(12)E(5) and C(18)E(50), respectively, which is in agreement with the corresponding values of 1.49 and 1.15 mg/m(2) obtained from applying the Gibbs equation to interfacial tension data for the same systems. We will briefly discuss the adsorption behavior in relation to the molecular structure of the surfactant and the phase behavior of the oil--surfactant--aqueous systems in relation to our experimental results.     

Self-assembly of a nonionic surfactant at the graphite/ionic liquid interface

In this study, we demonstrate by AFM imaging that nonionic surfactants self-assemble into hemicylindrical aggregates at the interface between graphite and the room temperature ionic liquid ethylammonium nitrate. Like aqueous systems, surfactant first adsorbs in a tail-to-tail monolayer arrangement along one of the three symmetry axes of graphite, templating subsequent self-assembly into adsorbed hemicylinders. Longer surfactant tails and higher concentrations are required to produce hemicylindrical aggregates in the ionic liquid than in aqueous solutions.     

Optimized ensemble Monte Carlo simulations of dense Lennard-Jones fluids

We apply the recently developed adaptive ensemble optimization technique to simulate dense Lennard-Jones fluids and a particle-solvent model by broad-histogram Monte Carlo techniques. Equilibration of the simulated fluid is improved by sampling an optimized histogram in radial coordinates that shifts statistical weight towards the entropic barriers between the shells of the liquid. Interstitial states in the vicinity of these barriers are identified with unprecedented accuracy by sharp signatures in the quickly converging histogram and measurements of the local diffusivity. The radial distribution function and potential of mean force are calculated to high precision.     

Hydration state of nonionic surfactant monolayers at the liquid/vapor interface: structure determination by vibrational sum frequency spectroscopy

The OH stretching region of water molecules in the vicinity of nonionic surfactant monolayers has been investigated using vibrational sum frequency spectroscopy (VSFS) under the polarization combinations ssp, ppp, and sps. The surface sensitivity of the VSFS technique has allowed targeting the few water molecules present at the surface with a net orientation and, in particular, the hydration shell around alcohol, sugar, and poly(ethylene oxide) headgroups. Dramatic differences in the hydration shell of the uncharged headgroups were observed, both in comparison to each another and in comparison to the pure water surface. The water molecules around the rigid glucoside and maltoside sugar rings were found to form strong hydrogen bonds, similar to those observed in tetrahedrally coordinated water in ice. In the case of the poly(ethylene oxide) surfactant monolayer a significant ordering of both strongly and weakly hydrogen bonded water was observed. Moreover, a band common to all the surfactants studied, clearly detected at relatively high frequencies in the polarization combinations ppp and sps, was assigned to water species located in proximity to the surfactant hydrocarbon tail phase, with both hydrogen atoms free from hydrogen bonds. An orientational analysis provided additional information on the water species responsible for this band.     

A Monte Carlo study of proton adsorption at the heterogeneous oxide/electrolyte interface

Adsorption of protons on a heterogeneous solid surface is modeled using the Monte Carlo (MC) simulation method. The surface of an oxide is assumed to consist of adsorption sites with pK assigned according to a quasi-Gaussian distribution. The influence of the electrostatic interactions combined with the energetic heterogeneity of the surface is examined, and the MC results are compared with the predictions of the mean field theory (MFT). It is demonstrated that the heterogeneity affects strongly the shape of the isotherms while it does not change the location of the common intersection point of the isotherms. On the other hand, introduction of repulsive interactions into the system is found to shift the CIP toward lower values of pH. It is also shown that the MFT, in general, describes correctly the behavior of the system. On the contrary the condensation approximation, used to derive relatively simple expressions for the adsorption isotherms, introduces serious errors unless the surface is strongly heterogeneous. Some practical remarks how to eliminate the errors associated both with the MC simulations and with the theory are also presented.     

Monte Carlo simulations of vapor–liquid–liquid equilibrium of some ternary petrochemical mixtures

The aim of this study was to determine the capability and accuracy of Monte Carlo simulations to predict ternary vapor–liquid–liquid equilibrium (VLLE) for some industrial systems. Hence, Gibbs ensemble Monte Carlo simulations in the isobaric–isothermal (NpT) and isochoric–isothermal (NVT) ensembles were performed to determine vapor–liquid–liquid equilibrium state points for three ternary petrochemical mixtures: methane/n-heptane/water (1), n-butane/1-butene/water (2) and n-hexane/ethanol/water (3). Since mixture (1) exhibits a high degree of mutual insolubility amongst its components, and hence has a large three-phase composition region, simulations in the NpT ensemble were successful in yielding three distinct and stable phases at equilibrium. The results were in very good agreement with experimental data at 120 kPa, but minor deviations were observed at 2000 kPa. Obtaining three phases for mixture (2) with the NpT ensemble is very difficult since it has an extremely narrow three-phase region at equilibrium, and hence the NVT ensemble was used to simulate this mixture. The simulated results were, once again, in excellent agreement with experimental data. We succeeded in obtaining three-phase equilibrium in the NpT ensemble only after knowing, a priori, the correct three-phase pressure (corresponding to the force fields that were implemented) from NVT simulations. The success of the NVT simulation, compared to NpT, is due to the fact that the total volume can spontaneously partition itself favorably amongst the three boxes and only one intensive variable (T) is fixed, whereas the pressure and the temperature are fixed in an NpT simulation. NpT simulations yielded three distinct phases for mixture (3), but quantitative agreement with experimental data was obtained at very low ethanol concentrations only.     

Tunable synergism/antagonism in a mixed nonionic/anionic surfactant layer at the solid/liquid interface

The use of mixed surfactants for modification of solid surfaces is important for many applications, since beneficial synergism often occurs depending on the surfactant type and mixing conditions. Systematical information on the properties of surfactant mixtures at the solid/liquid interface can be helpful for optimizing the interactions between the surfactants and then their corresponding performance. In this work, a nonionic/anionic surfactant combination, n-dodecyl beta-d-maltoside (DM) and sodium dodecyl sulfonate (SDS), was selected for the study of adsorption on an oxide solid, alumina. Interestingly, the mixture of the two surfactants with opposite pH-dependence of adsorption on alumina exhibits some unique synergistic or antagonistic features that were found to be tunable in the region of pH 4-10. In addition, the DM/SDS molar ratio in the adsorbed layer was found to decrease with concentration in the saturated region at all the pH and mixing ratios tested. The decrease is attributed to the monomer concentration changes in solution due to the difference in surface activities of the two surfactants. The tunable features of this mixture at the solid/liquid interface provide a way to optimize the properties by changing the mixing conditions. This can be valuable in many applications, such as enhanced oil recovery, flotation, and solubilization.     

Study of proton adsorption at heterogeneous oxide/electrolyte interface. Prediction of the surface potential using Monte Carlo simulations and 1-pK approach

Adsorption of protons on a heterogeneous solid surface is modeled using the Monte Carlo (MC) simulation method. The surface of an oxide is assumed to consist of adsorption sites with pK assigned according to a quasi-Gaussian distribution. The influence of the electrostatic interactions combined with the energetic heterogeneity of the surface is examined and the MC results are compared with the predictions of the analytical 1-pK approach. The surface potential behavior is examined using both "experimental" MC results and "theoretical" results obtained from the application of 1-pK model. The results are compared qualitatively with experimental determination of the surface potential of metal oxide surfaces. They confirm that the relation between the surface potential and the pH of bulk solution should not be described by the Nernst equation but by the equation with the parameter linearly reducing Nerstian potential. The values of this parameter are examined with respect to degree of surface energetic heterogeneity and site density of the surface.     

Monte Carlo study of surfactant adsorption on heterogeneous solid surfaces

The equilibrium between free surfactant molecules in aqueous solution and adsorbed layers on structured solid surfaces is investigated by lattice Monte Carlo simulation. The solid surfaces are composed of hydrophilic and hydrophobic surface regions. The structures of the surfactant adsorbate above isolated surface domains and domains arranged in a checkerboard-like pattern are characterized. At the domain boundary, the adsorption layers display a different behavior for hydrophilic and hydrophobic surface domains. For the checkerboard-like surfaces, additional adsorption takes place at the boundaries between surface domains.     

Adsorption of nonionic surfactant mixtures at the hydrophilic solid-solution interface

The adsorption of the mixed nonionic surfactants, monododecyl triethylene glycol (C2EO3) and monododecyl octaethylene glycol (C12EO8), at the hydrophilic silica-solution interface has been studied by specular neutron reflectivity. The adsorption at the solid-solution interface is compared with that previously measured at the air-solution interface. The marked differences that are observed are explained in terms of the different packing constraints or preferred curvature arising from the disparity in the respective headgroup dimensions.     

Monte Carlo simulations of liquid crystals near rough walls

The effect of surface roughness on the structure of liquid crystalline fluids near solid substrates is studied by Monte Carlo simulations. The liquid crystal is modeled as a fluid of soft ellipsoidal molecules and the substrate is modeled as a hard wall that excludes the centers of mass of the fluid molecules. Surface roughness is introduced by embedding a number of molecules with random positions and orientations within the wall. It is found that the density and order near the wall are reduced as the wall becomes rougher, i.e., the number of embedded molecules is increased). Anchoring coefficients are determined from fluctuations in the reciprocal space order tensor. It is found that the anchoring strength decreases with increasing surface roughness.     

Monte Carlo simulations of liquid crystals between microstructured substrates

The structure of a model liquid crystalline fluid confined between two microstructured substrates is studied through Monte Carlo simulations. A simple model for a structured substrate, similar in spirit to those used for rough walls of walls with grafted polymers, is introduced. It is found that varying the structure of the substrate, a transition in the alignment of the confined fluid, from parallel to perpendicular, is induced. For particular substrate structures, it is possible to induce tilted alignment in the confined fluid, the tilt angle being temperature dependent.     

Two-dimensional Monte Carlo simulations of ionic and nonionic silane self-assembly on hydrophilic surfaces

Aqueous chemistries have recently been shown to be useful for the deposition of hydrophobic films of nonionic and cationic silanes on hydrophilic substrates for the prevention of stiction in MEMS. The Monte Carlo method is used to simulate in two dimensions the self-assembly of silane films on a hydrophilic surface. We investigate the impact of charged group in cationic silane on the overall structure of the films. We characterize the film structure with spatial pair correlations at each molecular layer of the deposited films. The simulations reveal long-range correlations for the film of cationic silanes. Based on our two-dimensional simulations, we report an average "most probable" structure for the films of nonionic and cationic silanes.     

Solute rotational dynamics at the water liquid/vapor interface

The rotational dynamics of a number of diatomic molecules adsorbed at different locations at the interface between water and its own vapors are studied using classical molecular dynamics computer simulations. Both equilibrium orientational and energy correlations and nonequilibrium orientational and energy relaxation correlations are calculated. By varying the dipole moment of the molecule and its location, and by comparing the results with those in bulk water, the effects of dielectric and mechanical frictions on reorientation dynamics and on rotational energy relaxation can be studied. It is shown that for nonpolar and weekly polar solutes, the equilibrium orientational relaxation is much slower in the bulk than at the interface. As the solute becomes more polar, the rotation slows down and the surface and bulk dynamics become similar. The energy relaxation (both equilibrium and nonequilibrium) has the opposite trend with the solute dipole (larger dipoles relax faster), but here again the bulk and surface results converge as the solute dipole is increased. It is shown that these behaviors correlate with the peak value of the solvent-solute radial distribution function, which demonstrates the importance of the first hydration shell structure in determining the rotational dynamics and dependence of these dynamics on the solute dipole and location.     

Structure of the acetone liquid-vapor interface as seen from Monte Carlo simulations

The orientational order of the molecules at the liquid-vapor interface of acetone has been investigated by computer simulation. To fully describe the orientational preferences of the acetone molecules, the bivariate joint distribution of two independent orientational parameters has been determined at different layers of the interface. The strength of the orientational ordering of the interfacial molecules has been found to be liquid-like rather than crystal-like. The obtained results have revealed that the interfacial acetone molecules have dual orientational preferences. The main symmetry axis of the molecules declines by about 50-70 degrees from the interface normal axis, pointing toward the liquid phase in both of the preferred orientations. However, the plane of the molecules in the orientation preferred on the liquid side of the interface is perpendicular to the interfacial plane, whereas the other preferred orientation, which is present on the vapor side of the interface, corresponds to the alignment obtained from this orientation by an almost 90 degrees rotation around the main symmetry axis. Because the population of the liquid side is higher than that of the vapor side of the interface, the first of the two preferred orientations is the dominant alignment over the entire interface, in good agreement with recent experimental findings (Chen, H.; Gan, W.; Wu, B. H.; Wu, D.; Zhang, Z.; Wang, H. F. Chem. Phys. Lett. 2005, 408, 284).     

Monte Carlo simulations of hydrogen adsorption in alkali-doped single-walled carbon nanotubes

Monte Carlo simulations and Widom's test particle insertion method have been used to calculate the solubility coefficients (S) and the adsorption equilibrium constants (K) in single-walled (10,10) armchair carbon nanotubes including single nanotubes, and nanotube bundles with various configurations with and without alkali dopants. The hydrogen adsorption isotherms at room temperature were predicted by following the Langmuir adsorption model using the calculated constants S and K. The simulation results were in good agreement with experimental data as well as the grand canonical Monte Carlo simulation results reported in the literature. The simulations of nanotube bundle configurations suggest that the gravimetric hydrogen adsorption increases with internanotube gap size. It may be attributed to favorable hydrogen-nanotube interactions outside the nanotubes. The effect of alkali doping on hydrogen adsorption was studied by incorporating K+ or Li+ ions into nanotube arrays using a Monte Carlo simulation. The results on hydrogen adsorption isotherms indicate hydrogen adsorption of 3.95 wt% for K-doping, and 4.21 wt% for Li-doping, in reasonable agreement with the experimental results obtained at 100 atm and room temperature.