Sodium Dodecyl Sulfate at Water–Hydrophobic Interfaces: A Simulation Study

Exploration of the Adsorption Kinetics of Surfactants at the Water-Oil Interface via Grand-Canonical Molecular Dynamics Simulations

It is well-known that surfactants tend to aggregate into clusters or micelles in aqueous solutions due to their special structures, and it is difficult for the surfactant molecules involved in the aggregation to move spontaneously to the oil-water interface. In this article, we developed a new grand-canonical molecular dynamics (GCMD) model to predict the saturated adsorption amount of surfactant with constant concentration of surfactant molecules in the bulk phase, which can prevent surfactants aggregating in the bulk phase and get the atomic details of the interfacial structural change with increase of the adsorption amount through a single GCMD run. The adsorption of anionic surfactant sodium dodecyl sulfate (SDS) at the heptane-water interface was studied to validate the model. The saturated adsorption amount obtained from the GCMD simulation is consistent with the experimental results. The adsorption kinetics of SDS molecules during the simulation can be divided into three stages: linear adsorption stage, transition adsorption stage, and dynamic equilibrium stage. We also carried out equilibrium molecular dynamics (EMD) simulations to compare with GCMD simulation. This GCMD model can effectively reduce the simulation time with correct prediction of the interfacial saturation adsorption. We believe the GCMD method could be especially helpful for the computational study of surfactant adsorption under complex environments or emulsion systems with the adsorption of multiple types of surfactants.

Nanoscale Effects on the Surfactant Adsorption and Interface Charging in Hexadecane/Water Emulsions

Nanoscale properties at interfaces play a key role in the colloidal stability of emulsions and other soft matter materials where physical properties need to be controlled from the nano to macroscopically visible length scales. Our molecular level understanding of oil-water interfaces arises mostly from results at extended interfaces and the common view that emulsions are stabilized by a large number of surfactant molecules at the droplet's interface which, however, has been recently challenged. In this work, we show that the particle size and the curvature of oil droplets at the nanoscale is of great importance for the interface adsorption of dodecyl sulfate surfactants and possible counterion condensation at the charged hexadecane-water interface. Using second-harmonic scattering, we have studied the surface charge of oil droplets in nanoemulsions where we systematically varied the particle size R between 80 and 270 nm and demonstrate that the surface charge density σ changes drastically with size: For sizes >200 nm, σ is similar to what can be expected at flat extended interfaces, while σ is dramatically reduced by almost an order of magnitude when the particle size of the oil droplet is 80 nm. Using a theoretical approach that considers counterion condensation, we quantify the nanoscale effects on the change in surface charge with particle size and find excellent agreement with our experimental result. Modeling of the experimental results also implies that the charge per particle remains constant and depends on a critical balance of surfactant adsorption and ion condensation.

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.

Chemical computational approaches for optimization of effective surfactants in enhanced oil recovery

The surfactant flooding becomes an attractive method among several Enhanced Oil Recovery (EOR) processes to improve the recovery of residual oil left behind in the reservoir after secondary oil recovery process. The designing of a new effective surfactant is a comparatively complex and often time consuming process as well as cost-effective due to its dependency on the crude oil and reservoir properties. An alternative chemical computational approach is focused in this article to optimize the performance of effective surfactant system for EOR. The molecular dynamics (MD), dissipative particle dynamics (DPD) and density functional theory (DFT) simulations are mostly used chemical computational approaches to study the behaviour in multiple phase systems like surfactant/oil/brine. This article highlighted a review on the impact of surfactant head group structure on oil/water interfacial property like interfacial tensions, interface formation energy, interfacial thickness by MD simulation. The effect of entropy in micelle formation has also discussed through MD simulation. The polarity, dipole moment, charge distribution and molecular structure optimization have been illustrated by DFT. A relatively new coarse-grained method, DPD is also emphasized the phase behaviour of surfactant/oil/brine as well as polymer-surfactant complex system.

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.

Lyophobicity may not be the main driving force for long chain surfactants from the bulk phase to the interface

According to the Traube rule, a surfactant with a longer alkane chain is more hydrophobic so its tendency to be driven from a polar solvent to a less polar interface is higher. In this work, we revisited this topic by studying the adsorption of quaternary ammonium salts and carboxylic acids with various alkane chain lengths at the hexadecane-water interface. The adsorption free energies of the surfactants at this oil-water interface from the polar (aqueous solution) or nonpolar phase (hexadecane) were estimated from second harmonic generation measurements. The variation of the free energies per methylene group in the bulk phase, at the oil-water interface and at the air-water interface revealed that there are different interactions between the alkane chains of the surfactants in different environments. The chain-chain interaction at the hexadecane-water interface is lower than that at the air-water interface. The driving force for the alkane chains to adsorb at the oil-water interface from the oil phase is close to that from the aqueous phase. This observation reveals that the chain-chain interaction rather than the lyophobicity of the solute with respect to the solvent is the main contributor to the adsorption free energy. This is the first experimental comparison of the free energies of the alkane chains in oil, in water, at the air-water interface and at the oil-water interface. These results provide information for studying the interactions of hydrophobic species in different environments. This work also provides a method for estimating the solvation energy of some head groups in surfactants.

The Molecular Mechanism of Nanodroplet Stability

Mixtures of nano- and microscopic oil droplets in water have recently been rediscovered as miniature reaction vessels in microfluidic environments and are important constituents of many environmental systems, food, personal care, and medical products. The oil nanodroplet/water interface stabilized by surfactants determines the physicochemical properties of the droplets. Surfactants are thought to stabilize nanodroplets by forming densely packed monolayers that shield the oil phase from the water. This idea has been inferred from droplet stability measurements in combination with molecular structural data obtained from extended planar interfaces. Here, we present a molecular level investigation of the surface structure and stability of nanodroplets and show that the surface structure of nanodroplets is significantly different from that of extended planar interfaces. Charged surfactants form monolayers that are more than 1 order of magnitude more dilute than geometrically packed ones, and there is no experimental correlation between stability and surfactant surface density. Moreover, dilute negatively charged surfactant monolayers produce more stable nanodroplets than dilute positively charged and dense geometrically packed neutral surfactant monolayers. Droplet stability is found to depend on the relative cooperativity between charge-charge, charge-dipole, and hydrogen-bonding interactions. The difference between extended planar interfaces and nanoscale interfaces stems from a difference in the thermally averaged total charge-charge interactions in the two systems. Low dielectric oil droplets with a size smaller than the Debye length in oil permit repulsive interactions between like charges from opposing interfaces in small droplets. This behavior is generic and extends up to the micrometer length scale.

Electrokinetic transport in liquid foams

Investigating electrokinetic transport in a liquid foam is at the confluence of two well developed research areas. On one hand, the study of electrokinetic flows (i.e. surface-driven flows generated close to a charged interface) is fairly well understood in regards the solid/liquid interface. On the other hand, the flow of liquid in a 3D deformable network, i.e a foam, under a volume force such as gravity has been thoroughly studied over the past decade. The overlapping zone of these two frameworks is of great interest for both communities as it gives rise to challenging new questions such as: what is the importance of the nature of the charged interface, created by mobile and soluble surfactants in the case of foam, on electrokinetic transport? How does a foam behave when submitted to a surface-driven flow? Can we compensate a volume-driven flow, i.e. gravity, by a surface-driven flow, i.e. electroosmosis? In this review, we will explore these questions on three different scales: a surfactant laden interface, a foam film and a macroscopic foam.

Effects of Frothers and Oil at Saltwater-Air Interfaces for Oil Separation: Molecular Dynamics Simulations and Experimental Measurements

Separating oil from saltwater is a process relevant to some industries and may be aided by bubble and froth generation. Simulating saltwater-air interfaces adsorbed with surfactants and oil molecules can assist in understanding froth stability to improve separation. Combining with surface tension experimental measurements, in this work we employ molecular dynamics with a united-atom force field to linear alkane oil and three surfactant frothers, methyl isobutyl carbinol (MIBC), terpineol, and ethyl glycol butyl ether (EGBE), to investigate their synergistic behaviors for oil separation. The interfacial phenomena were measured for a range of frother surface coverages on saltwater. Density profiles of the hydrophilic and hydrophobic portions of the frothers show an expected orientation of alcohol groups adsorbing to the polar water. A decrease in surface tension with increasing surface coverage of MIBC and terpineol was observed and reflected in experiments where the frother concentration increased. Relations between surface coverage and bulk concentration were observed by comparing the surface tension decreases. Additionally, a range of oil surface coverages was explored when the interface has a thin layer of adsorbed frother molecules. The obtained results indicate that an increase in surface coverage of oil molecules led to an increase in surface tension for all frother types and the pair correlation functions depicted MIBC and terpineol as having higher distributions with water at closer distances than with oil.

Electro-osmosis at surfactant-laden liquid-gas interfaces: beyond standard models

Electro-osmosis (EO) is a powerful tool to manipulate liquids in micro and nanofluidic systems. While EO has been studied extensively at liquid-solid interfaces, the case of liquid-vapor interfaces, found e.g. in foam films and bubbles, remains to be explored. Here we perform molecular dynamics (MD) simulations of EO in a film of aqueous electrolyte covered with fluid layers of ionic surfactants and surrounded by gas. Following the experimental procedure, we compute the zeta potential from the EO velocity, defined as the velocity difference between the middle of the liquid film and the surrounding gas. We show that the zeta potential can be smaller or larger than existing predictions depending on the surfactant coverage. We explain the failure of previous descriptions by the fact that surfactants and bound ions move as rigid bodies and do not transmit the electric driving force to the liquid locally. Considering the reciprocal streaming current effect, we then develop an extended model, which can be used to predict the experimental zeta potential of surfactant-laden liquid-gas interfaces.

Exploring Fluorescent Dyes at Biomimetic Interfaces with Second Harmonic Generation and Molecular Dynamics

The adsorption of a DNA fluorescent probe belonging to the thiazole orange family at the dodecane/water and dodecane/phospholipid/water interfaces has been investigated using a combination of surface second harmonic generation (SSHG) and all-atomistic molecular dynamics (MD) simulations. Both approaches point to a high affinity of the cationic dye for the dodecane/water interface with a Gibbs free energy of adsorption on the order of -45 kJ/mol. Similar affinity was observed with a monolayer of negatively charged DPPG (1,2-dipalmitoyl-sn-glycero-3-phospho-rac-(1-glycerol)) lipids. On the other hand, no significant adsorption could be found with the zwitterionic DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphocholine) lipids. This was rationalized in terms of Coulombic interactions between the monolayer surface and the cationic dye. The similar affinity for the interface with and without DPPG, despite the favorable Coulombic attraction in the latter case, could be explained after investigating the interfacial orientation of the dye. In the absence of a monolayer, the dye adsorbs with its molecular plane almost flat at the interface, whereas in the presence of DPPG it has to intercalate into the monolayer and adopt a significantly different orientation to benefit from the electrostatic stabilization.

Molecular Dynamics Study of the Bulk and Interface Properties of Frother and Oil with Saltwater and Air

For water treatment purposes, the separation processes involving surfactants and crude oil at seawater-air interfaces are of importance for the chemical and energy industries. Little progress has been made in understanding the nanoscale phenomena of surfactants on oily saltwater-air interfaces. This work focuses on using molecular dynamics with a united-atom force field to simulate the interface of linear alkane oil, saltwater, and air with three surfactant frothers: methyl isobutyl carbinol (MIBC), terpineol, and ethyl glycol butyl ether. For each frother, although the calculated diffusivities and viscosities are lower than the expected experimental values, our results show that diffusivity trends between each frother agree with experiments but the method cannot be applied for viscosity. Binary combinations of liquid (frother or saltwater)-air and liquid-liquid interfaces are equilibrated to study the density profiles and interfacial tensions. The calculated surface tensions of the frother-air interfaces are like that of oil-air, but lower than that of saltwater-air. Only the MIBC-air and terpineol-air interfaces agreed with our experimental measurements. For the frother-saltwater interfaces, the calculated results showed that terpineol has interfacial tensions higher than those of MIBC-saltwater. The simulated results indicate that the frother-oil systems underwent mixing such that the density profiles depicted large interfacial thicknesses.

Assessment of micro-polarity anisotropy as a function of surfactant packing in sodium dodecyl sulphonate-hexane reverse micelles

The micro-polarity anisotropy behaviour across the aqueous phase of a SDS (sodium dodecyl sulphonate)-hexane reverse micelle (RM) relies on the SDS packing in the oil-water interfacial self-assembled surfactant structure of the RM.

Three Dimensional Nano "Langmuir Trough" for Lipid Studies

A three-dimensional-phospholipid monolayer with tunable molecular structure was created on the surface of oil nanodroplets from a mixture of phospholipids, oil, and water. This simple nanoemulsion preparation technique generates an in situ prepared membrane model system with controllable molecular surface properties that resembles a lipid droplet. The molecular interfacial structure of such a nanoscopic system composed of hexadecane, 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), and water was determined using vibrational sum frequency scattering and second harmonic scattering techniques. The droplet surface structure of DPPC can be tuned from a tightly packed liquid condensed phase like monolayer to a more dilute one that resembles the liquid condensed/liquid expanded coexistence phase by varying the DPPC/oil/water ratio. The tunability of the chemical structure, the high surface-to-volume ratio, and the small sample volume make this system an ideal model membrane for biochemical research.

Temperature-dependent phase transition and desorption free energy of sodium dodecyl sulfate at the water/vapor interface: approaches from molecular dynamics simulations

Adsorption of surfactants at the water/vapor interface depends upon their chemical potential at the interface, which is generally temperature-dependent. Molecular dynamics simulations have been performed to reveal temperature influences on the microstructure of sodium dodecyl sulfate (SDS) molecule adsorption layer. At room temperature, SDS molecules aggregate at the interface, being in a liquid-expanded phase, whereas they tend to spread out and probably transit to a gaseous phase as the temperature increases to above 318 K. This phase transition has been confirmed by the temperature-dependent changes in two-dimensional array, tilt angles, and immersion depths to the aqueous phase of SDS molecules. The aggregation of SDS molecules accompanies with larger immersion depths, more coordination of Na(+) ions, and less coordination of water. Desorption free energy profiles show that higher desorption free energy appears for SDS molecules at the aggregate state at low temperatures, but no energy barrier is observed. The shapes of desorption free energy profiles depend upon the distribution of SDS at the interface, which, in turn, is related to the phase state of SDS. Our study sheds light on the development of adsorption thermodynamics and kinetics theories.