Structures and Properties of Newton Black Films Characterized Using Molecular Dynamics Simulations

Review of the Interfacial Structure and Properties of Surfactants in Petroleum Production and Geological Storage Systems from a Molecular Scale Perspective

Surfactants play a crucial role in tertiary oil recovery by reducing the interfacial tension between immiscible phases, altering surface wettability, and improving foam film stability. Oil reservoirs have high temperatures and high pressures, making it difficult and hazardous to conduct lab experiments. In this context, molecular dynamics (MD) simulation is a valuable tool for complementing experiments. It can effectively study the microscopic behaviors (such as diffusion, adsorption, and aggregation) of the surfactant molecules in the pore fluids and predict the thermodynamics and kinetics of these systems with a high degree of accuracy. MD simulation also overcomes the limitations of traditional experiments, which often lack the necessary temporal-spatial resolution. Comparing simulated results with experimental data can provide a comprehensive explanation from a microscopic standpoint. This article reviews the state-of-the-art MD simulations of surfactant adsorption and resulting interfacial properties at gas/oil-water interfaces. Initially, the article discusses interfacial properties and methods for evaluating surfactant-formed monolayers, considering variations in interfacial concentration, molecular structure of the surfactants, and synergistic effect of surfactant mixtures. Then, it covers methods for characterizing microstructure at various interfaces and the evolution process of the monolayers' packing state as a function of interfacial concentration and the surfactants' molecular structure. Next, it examines the interactions between surfactants and the aqueous phase, focusing on headgroup solvation and counterion condensation. Finally, it analyzes the influence of hydrophobic phase molecular composition on interactions between surfactants and the hydrophobic phase. This review deepened our understanding of the micro-level mechanisms of oil displacement by surfactants and is beneficial for screening and designing surfactants for oil field applications.

Molecular Dynamics Simulation of the Effects of Complex Surfactants on Oil-Water Interaction and Aggregation Characteristics at the Interface

In response to the problem of complex interaction between oil and water in the oil-water interface, especially heavy oil and water, this study investigated the effects of complex surfactants on the interaction of two phases and their aggregation characteristics by molecular dynamics simulation. The results showed that increasing the content of sodium lauryl polyether carboxylate (AEC-9Na) was beneficial to the coordination between it and alkyl glycoside (APG-10), improved the interfacial activity, and enhanced the interfacial stability of the composite system, and the best effect was achieved when AEC-9Na:APG-10 = 8:2. The thickness of the oil and water film on the oil-water interface was irregular. When the concentration of AEC-9Na was lower than that of APG-10, the total thickness of the interfacial film (ttotal) first increased. When the content of AEC-9Na is higher, a large number of sodium ions were adsorbed near the -COO- group of AEC-9Na, which will polarize out of the hydration layer structure and attract water molecules from the second hydration layer on the heavy oil surface to the first hydration layer through electrostatic interaction. Then, the thickness of the interface film was compressed, and the interface film was reduced. When the ratio increased to 10:0, the oil and water phase competed to adsorb surfactant molecules, and the headgroup tended to lay on the interface. Moreover, the hydrophilicity of the surfactant layer was weakened, and the thickness of the water film decreased. The distribution of surfactant was looser than 8:2, the light components of heavy oil molecules (saturated and aromatic hydrocarbons) entered the gap between surfactants in large quantities, and the hydrophobic tail chain tended to be laid on the oil-water interface. The oleophilicity of the surfactant layer increased, and the thickness of the oil film remarkably increased, so the total thickness of the interface film increased again.

Contribution of air-water interface in removing PFAS from drinking water: Adsorption, stability, interaction and machine learning studies

As a class of synthetic persistent organic pollutants, contamination of Per-and poly-fluoroalkyl substances (PFAS) in drinking water has attracted widespread concern. Aeration has been confirmed to enhance the removal of PFAS in drinking water by activated carbon (AC). However, the contribution of the air-water interface in removing PFAS is not yet to be fully understood at the molecular level. In this work, molecular dynamics (MD) simulations were employed to investigate the role of nanobubble in removing PFAS in the aqueous environment. The result suggests that the free energies of the air-water interface are about 3-7 kcal mol^-1 lower than that of the bulk water region, indicating that the transformation of PFAS from the water phase into the air-water interface is favorable from the viewpoint of thermodynamics. The interface-water partition coefficients (P_sur/wat) of PFAS are in the order of PFOS > PFOA > PFHxS > PFBS. On the air-water-AC three-phase interface, PFBS can not only move along the interface region but also leave the interface region into water phase, while PFOS tended to move along the interface region until it was captured by AC. Finally, the ΔG_{water-interface} quantitative structure-activity relationships (QSAR) models were developed to predict the removal efficiencies of PFAS enhanced by aeration in aquatic systems. The proposed mechanism promotes the understanding of the contribution of air-water interface in removing PFAS from drinking water by activated carbon.

Exploring the effects of approach velocity on depletion force and coalescence in oil-in-water emulsions

An emulsion is a thermodynamically unstable system consisting of at least two immiscible liquid phases, one of which is dispersed in the other in the form of droplets of varying size. Most studies on emulsions have focused on the behaviour of emulsion droplets with diameter from ∼50 μm and upwards. However, the properties of smaller droplets may be highly relevant in order to understand the behaviour of emulsions, including their performance in numerous applications within the fields of food, industry, and medical science. The relatively long life-time and small size of these droplets compared to other emulsion droplets, make them suited for optical trapping and micromanipulation technologies. Optical tweezers have previously shown potential in the study of stabilized emulsions. Here we employ optical tweezers to examine unstable oil-in-water emulsions to determine the effects of system parameters on depletion force and coalescence times.

Effect of Sodium Dodecyl Sulfonate on the Foam Stability and Adsorption Configuration of Dodecylamine at the Gas-Liquid Interface

In this study, the effect of sodium dodecyl sulfonate (SDS) on the foam stability of dodecylamine (DDA) and on its adsorption configuration at the gas-liquid interface was investigated. Froth stability experiments, surface tension measurements, time-of-flight secondary-ion mass spectrometry measurements, and molecular dynamics simulation calculations were performed in this investigation. The results revealed that the foam stability of DDA solution was extremely strong, and the addition of SDS could decrease the foam stability when the concentration of DDA was less than a certain value. The decrease in foam stability could be ascribed to several reasons, namely, the big cross-sectional area of SDS at the gas-liquid interface and low adsorption capacity of surfactants at the gas-liquid interface, the high surface tension, the change in the double-layer structure, the small thickness of the gas-liquid interfacial layer, the weak interaction intensity between the head groups of the surfactants and the water molecules, the strong movement ability of the water molecules around the head groups, and the sparse and less upright arrangement configuration of molecules at the gas-liquid interface. These findings can greatly help in solving the strong foam stability problem in DDA flotation and provide a method for reducing foam stability.

Revealing the Origin of the Specificity of Calcium and Sodium Cations Binding to Adsorption Monolayers of Two Anionic Surfactants

The studied anionic surfactants linear alkyl benzene sulfonate (LAS) and sodium lauryl ether sulfate (SLES) are widely used key ingredients in many home and personal care products. These two surfactants are known to react very differently with multivalent counterions, including Ca²⁺. This is explained by a stronger interaction of the calcium cation with the LAS molecules, compared to SLES. The molecular origin of this difference in the interactions remains unclear. In the current study, we conduct classical atomistic molecular dynamics simulations to compare the ion interactions with the adsorption layers of these two surfactants, formed at the vacuum-water interface. Trajectories of 150 ns are generated to characterize the adsorption layer structure and the binding of Na⁺ and Ca²⁺ ions. We found that both surfactants behave similarly in the presence of Na⁺ ions. However, when Ca²⁺ is added, Na⁺ ions are completely displaced from the surface with adsorbed LAS molecules, while this displacement occurs only partially for SLES. The simulations show that the preference of Ca²⁺ to the LAS molecules is due to a strong specific attraction with the sulfonate head-group, besides the electrostatic one. This specific attraction involves significant reduction of the hydration shells of the interacting calcium cation and sulfonate group, which couple directly and form surface clusters of LAS molecules, coordinated around the adsorbed Ca²⁺ ions. In contrast, SLES molecules do not exhibit such specific interaction because the hydration shell around the sulfate anion is more stable, due to the extra oxygen atom in the sulfate group, thus precluding substantial dehydration and direct coupling with any of the cations studied.

Foam-in-vein: A review of rheological properties and characterization methods for optimization of sclerosing foams

Varicose veins are chronic venous defects that affect >20% of the population in developed countries. Among potential treatments, sclerotherapy is one of the most commonly used. It involves endovenous injection of a surfactant solution (or foam) in varicose veins, inducing damage to the endothelial layer and subsequent vessel sclerosis. Treatments have proven to be effective in the short-term, however recurrence is reported at rates of up to 64% 5-year post-treatment. Thus, once diagnosed with varicosities there is a high probability of a permanently reduced quality of life. Recently, foam sclerotherapy has become increasingly popular over its liquid counterpart, since foams can treat larger and longer varicosities more effectively, they can be imaged using ultrasound, and require lower amounts of sclerosing agent. In order to minimize recurrence rates however, an investigation of current treatment methods should lead to more effective and long-lasting effects. The literature is populated with studies aimed at characterizing the fundamental physics of aqueous foams; nevertheless, there is a significant need for appropriate product development platforms. Despite successfully capturing the microstructural evolution of aqueous foams, the complexity of current models renders them inadequate for pharmaceutical development. This review article will focus on the physics of foams and the attempts at optimizing them for sclerotherapy. This takes the form of a discussion of the most recent numerical and experimental models, as well as an overview of clinically relevant parameters. This holistic approach could contribute to better foam characterization methods that patients may eventually derive long term benefit from.

Unraveling the Differential Aggregation of Anionic and Nonionic Monorhamnolipids at Air-Water and Oil-Water Interfaces: A Classical Molecular Dynamics Simulation Study

The molecular structure of a surfactant molecule is known to have a great effect on the interfacial properties. We employ molecular dynamics simulations for a detailed atomistic study of monolayers of the nonionic and anionic form of the most common congener of monorhamnolipids, α-rhamnopyranosyl-β-hydroxydecanoyl-β-hydroxydecanoate (( R, R)-Rha-C10-C10), at the air-water and oil-water interfaces. An atomistic-level understanding of monolayer aggregation is necessary to explain a recent experimental observation indicating that nonionic and anionic Rha-C10-C10 show surprisingly different surface area per molecule at the critical micelle concentration. Surface-pressure analysis, interface formation energy calculations, and mass density profiles of the monolayers at the air-water interface show similar properties between nonionic and anionic Rha-C10-C10 aggregation. It is found that there is a significant difference in the headgroup conformations of Rha-C10-C10 in the nonionic and anionic monolayers. Hydrogen bonding interactions between the Rha-C10-C10 molecules in the monolayers is also significantly different between nonionic and anionic forms. Representative snapshots of the simulated system at different surface concentrations show the segregation of molecular aggregates from the interface into the bulk water in the anionic Rha-C10-C10 monolayer at higher concentrations, whereas in the nonionic Rha-C10-C10 monolayer, the molecules are still located at the interface. The present work provides insight into the different aggregation properties of nonionic and anionic Rha-C10-C10 at the air-water interface. Further analyses were carried out to understand the aggregation behavior of nonionic and anionic Rha-C10-C10 at the oil-water interface. It is observed that the presence of oil molecules does not significantly influence the aggregation properties of Rha-C10-C10 as compared to those of the air-water interface.

Achieving foaming control smartly: pre-solubilized flavor oil serves as an in situ homogeneous defoamer

In the wide application of aqueous foam, creating abundant foam and processing appropriate foaming control are both essential, depending upon the actual situation; the latter process is not only harder to achieve, but also more complicated to comprehensively understand on the molecular level. In this paper, a type of natural flavor oil, carvone, was solubilized in a micelle solution of sodium dodecyl sulfate (SDS) to study the effect on the foaming properties. The foamability and foam stability of the swollen micelle solutions were experimentally characterized, and the molecular behavior of the surfactant and oil molecules before, during and after the foaming process were investigated. It was found that the solubilized carvone co-adsorbed with SDS at the gas/water interface and caused a prominent effect on the foam film stability in several approaches, thereby making the flavor oil a possible foam controller that would not inhibit foam formation, but could eliminate foam efficiently once foam was undesired. Interestingly, it was found that the release of flavor in the foaming process was promoted. Detailed discussion of the interfacial behavior of carvone and the effect on the foaming properties of surfactants in different stages of foam may provide a theoretical foundation for exploring green and smart approaches in achieving foaming control.

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.

Surfactant films in lyotropic lamellar (and related) phases: Fluctuations and interactions

The analogy between soap films thinning under border capillary suction and lamellar stacks of surfactant bilayers dehydrated by osmotic stress is explored, in particular in the highly dehydrated limit where the soap film becomes a Newton black film. The nature of short-range repulsive interactions between surfactant-covered interfaces and acting across water channels in both cases will be discussed.

Understanding about How Different Foaming Gases Effect the Interfacial Array Behaviors of Surfactants and the Foam Properties

In this paper, the detailed behaviors of all the molecules, especially the interfacial array behaviors of surfactants and diffusion behaviors of gas molecules, in foam systems with different gases (N2, O2, and CO2) being used as foaming agents were investigated by combining molecular dynamics simulation and experimental approaches for the purpose of interpreting how the molecular behaviors effect the properties of the foam and find out the key factors which fundamentally determine the foam stability. Sodium dodecyl sulfate SDS was used as the foam stabilizer. The foam decay and the drainage process were determined by Foamscan. A texture analyzer (TA) was utilized to measure the stiffness and viscoelasticity of the foam films. The experimental results agreed very well with the simulation results by which how the different gas components affect the interfacial behaviors of surfactant molecules and thereby bring influence on foam properties was described.

Benchmarking the Self-Assembly of Surfactin Biosurfactant at the Liquid-Air Interface to those of Synthetic Surfactants

The adsorption of surfactin, a lipopeptide biosurfactant, at the liquid–air interface has been investigated in this work. The maximum adsorption density and the nature and the extent of lateral interaction between the adsorbed surfactin molecules at the interface were estimated from surface tension data using the Frumkin model. The quantitative information obtained using the Frumkin model was also compared to those obtained using the Gibbs equation and the Langmuir–Szyszkowski model. Error analysis showed a better agreement between the experimental and the calculated values using the Frumkin model relative to the other two models. The adsorption of surfactin at the liquid–air interface was also compared to those of synthetic anionic, sodium dodecylbenzenesulphonate (SDBS), and nonionic, octaethylene glycol monotetradecyl ether (C₁₄E₈), surfactants. It has been estimated that the area occupied by a surfactin molecule at the interface is about 3- and 2.5-fold higher than those occupied by SDBS and C₁₄E₈ molecules, respectively. The interaction between the adsorbed molecules of the anionic biosurfactant (surfactin) was estimated to be attractive, unlike the mild repulsive interaction between the adsorbed SDBS molecules.

Slow dynamics of water confined in Newton black films

Slowdown of translational and reorientational dynamics of water confined in Newton black films (NBFs) is revealed by molecular dynamics simulations. As a film becomes thinner, both translational and reorientational dynamics become slower. The polarization of water molecules in the macroscopic electrostatic field across the NBF and the coordination of Na(+) ions and surfactant anionic groups around water molecules concertedly lead to slowdown of water dynamics. The polarization effect is obvious for water not coordinated by Na(+) ions, which exhibits reorientational dynamics depending on initial dipole orientations. Na(+) ions and surfactant anionic groups retard dynamics of surrounding water by decreasing the hydrogen bond exchange probability and increasing the viscosity of water. The dependences of translational and reorientational dynamics on coordination environments of water are similar. Dynamics of water in positions close to the interfaces of NBFs are mainly retarded by Na(+) ions and surfactant anionic groups, while the macroscopic polarization effect plays the main role in influencing water dynamics in positions far from the interfaces. This study sheds light on the improvement of knowledge about the water dynamics slowdown mechanism in similar environments like reverse micelles and lamellar structures.

Surface force at the nano-scale: observation of non-monotonic surface tension and disjoining pressure

The disjoining pressures of thin aqueous salt films at different salt concentrations and temperatures were calculated using MD simulations.

Film and membrane-model thermodynamics of free thin liquid films

In spite of over 7 decades of effort, the thermodynamics of thin free liquid films (as in emulsions and foams) lacks clarity. Following a brief review of the meaning and measurement of thin-film forces (i.e., conjoining/disjoining pressures), we offer a consistent analysis of thin-film thermodynamics. By carefully defining film reversible work, two distinct thermodynamic formalisms emerge: a film model with two zero-volume membranes each of film tension γ(f) and a membrane model with a single zero-volume membrane of membrane tension 2γ(m). In both models, detailed thermodynamic analysis gives rise to thin-film Gibbs adsorption equations that allow calculation of film and membrane tensions from measurements of disjoining-pressure isotherms. A modified Young-Laplace equation arises in the film model to calculate film-thickness profiles from the film center to the surrounding bulk meniscus. No corresponding relation exists in the membrane model. Illustrative calculations of disjoining-pressure isotherms for water are presented using square-gradient theory. We report considerable deviations from Hamaker theory for films less than about 3 nm in thickness. Such thin films are considerably more attractive than in classical Hamaker theory. Available molecular simulations reinforce this finding.

Interaction of Refractory Dibenzothiophenes and Polymerizable Structures

We carried out first principles calculations to show that polymerizable structures containing hydroxyl (alcoholic chain) and amino groups are suitable to form stable complexes with dibenzothiophene (DBT) and its alkyl derivates. These sulfur pollutants are very difficult to eliminate through traditional catalytic processes. Spontaneous and exothermic interactions at 0 K primarily occur through the formation of stable complexes of organosulfur molecules with monomeric structures by hydrogen bonds. The bonds are formed between the sulfur atom and the hydrogen of the hydroxyl group; additional hydrogen bonds are formed between the hydrogen atoms of the organosulfur molecule and the nitrogen atoms of the monomers. We vary the number of methylene groups in the alcoholic chain containing the hydroxyl group of the monomer and find that the monomeric structure with four methylene groups has the best selectivity towards the interaction with the methyl derivates with reference to the interaction with DBT. Even this study does not consider solvent and competitive adsorption effects; our results show that monomeric structures containing amino and hydroxyl groups can be used to develop adsorbents to eliminate organosulfur pollutants from oil and its derivates.

Stability of Newton black films under mechanical stretch--a molecular dynamics study

The stability of Newton black films (NBFs) under lateral mechanical stretch was investigated by nonequilibrium molecular dynamics (NEMD) simulations using force field parameters validated by accurate prediction of surface tensions. The applied strains accelerated film ruptures, enabling efficient measurements of the critical thicknesses of the films. Two representative surfactants, sodium dodecyl sulfate (SDS) for ionic surfactant and pentaethylene glycol monododecyl ether (C12EO5H) for nonionic surfactant, were investigated and compared. The predicted critical thickness of C12EO5H-coated film is smaller than that of the SDS-coated film, which is consistent with previously reported experimental observations. Our simulation results show that while the two surfactant-coated films exhibit similar dynamic properties attributed to the Marangoni-Gibbs effect, their surface structural characteristics are quite different. Consequently the two films demonstrate distinct rupture mechanisms in which rupture starts at uncovered water domains in the SDS-coated film, but at lateral surfactant/water interfaces in the C12EO5H-coated film. Our findings provide new insights into the stabilization mechanisms of NBFs and will facilitate the design and development of new films with improved properties.

Study of the microcharacter of ultrastable aqueous foam stabilized by a kind of flexible connecting bipolar-headed surfactant with existence of magnesium ion

In this paper, ultrastable aqueous foam stabilized by a kind of flexible connecting bipolar-headed surfactant alkyl polyoxyethylene sulfate (AE3S) with coexisting Mg(2+) was reported. Detailed molecular behaviors of AE3S in foam film with coexisting divalent cationic Ca(2+) or Mg(2+) were investigated by molecular dynamic simulation, comparing with the traditional surfactant sodium dodecyl sulfate (SDS), to find out how the microcharacter and array behavior of molecules in the foam film determined by molecular interaction effect the foam stability. It was found that the ultrastable foam film obtained by the cooperation of magnesium ions and AE3S was driven from two aspects: one is the favorable arrangement of surfactant molecules, and the other is the increase of capacity of foam films for resolutely holding water molecules deduced by a dipolar pair formed by the flexible connecting head groups of AE3S and hydrated Mg(2+) via intermolecular coactions, both related to the presence of magnesium ions. Foam lamella stability measurement and foam decay method were both used to evaluate the stability of foam. Fourier transform infrared (FT-IR) was used to detect the composition variation of foam film in the drainage process; the vibration peak of OH for water molecule shifted from the 3390 cm(-1) (being assigned to the bulk water integrated by hydrogen bonds) to 3685 cm(-1) (being assigned to the vibration of isolated water molecules) for the ultrastable foam film after complete drainage, which agreed very well with the molecular simulation results.