Structure, Interfacial Properties, and Dynamics of the Sodium Alkyl Sulfate Type Surfactant Monolayer at the Water/Trichloroethylene Interface: A Molecular Dynamics Simulation Study

Experimental and simulation-based characterization of surfactant adsorption layers at fluid interfaces

Adsorption of surfactants to fluid interfaces occurs in numerous technological and daily-life contexts. The coverage at the interface and other properties of the formed adsorption layers determine the performance of a surfactant with regard to the desired application. Given the importance of these applications, there is a great demand for the comprehensive characterization and understanding of surfactant adsorption layers. In this review, we provide an overview of suitable experimental and simulation-based techniques and review the literature in which they were used for the investigation of surfactant adsorption layers. We come to the conclusion that, while these techniques have been successfully applied to investigate Langmuir monolayers of water-insoluble surfactants, their application to the study of Gibbs adsorption layers of water-soluble surfactants has not been fully exploited. Finally, we emphasize the great potential of these methods in providing a deeper understanding of the behavior of soluble surfactants at interfaces, which is crucial for optimizing their performance in various applications.

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.

Membrane stabilization versus perturbation by aromatic monoamine-modified γ-PGA for cryopreservation of human RBCs with high intracellular trehalose

As a nonreducing disaccharide, trehalose can be used as a biocompatible cryoprotectant for solvent-free cell cryopreservation, but the membrane-impermeability limits its cryoprotective efficiency. Herein, a series of aromatic monoamines with a 1-4 methylene spacer were grafted onto γ-poly(glutamic acid) (γ-PGA) for promoting intracellular trehalose uptake in human red blood cells (hRBCs) via membrane perturbation. The self-assembled nanoparticles of the obtained amphiphilic γ-PGA could be adsorbed on the cell membrane by the hydrophobic interaction to disturb the lipid arrangement and increase the membrane permeability of trehalose under hypertonic conditions. Results suggested that the intracellular trehalose could be enhanced progressively with the methylene spacer length, significantly increasing to 75.1 ± 0.7 mM by incubating hRBCs in 0.8 M trehalose containing phenylbutylamine-grafted γ-PGA at 4 °C for 24 h. Meanwhile, the other three polymers exhibited membrane stabilization in addition to improved intracellular trehalose, maintaining the membrane integrity during cryopreservation to achieve high cryosurvival. Molecular dynamics simulation further confirmed that defects could be formed by interaction of the above four amphiphilic polymers on the modeled phospholipid bilayer. It was believed that glycerol-free cryopreservation of human cells could be realized by using trehalose as the biocompatible cryoprotectant, and membrane stabilization can be a compensatory approach to membrane perturbation during impermeable biomolecule delivery.

Improving 1-propanol force field: a new methodology

A new force field for 1-propanol, in the united and all atom models, has been obtained by combining two different empirical methodologies. The first was developed by scaling atom charges and Lennard-Jones parameters to fit the dielectric constant, surface tension, and density; this methodology is named three steps systematic parameterization procedure (3SSPP), as reported by Pérez de la Luz et al. (J Chem Theory Comput 14:5949-5958, 2018). The second methodology consists of moving these parameters and together with the bond distance to obtain the liquid-vapor phase diagram of the CO₂ molecule as discussed by Harris and Yung (J Phys Chem 99:12021-12024, 1995). The last methodology is used to obtain the self-diffusion coefficient, which was not consider in the 3SSPP. The 3SSPP/bond methodology is the 3SSPP plus the bond distance scaling. With this new methodology, the experimental density, dielectric constant, surface tension, and self-diffusion coefficient at ambient temperature could be achieved. Furthermore, we show the temperature dependence of the aforementioned properties. The static structure factors are in accordance with the experimental spectrum. Solubility is increased to the experimental value for the united atom (UA) model after applying this methodology and for all atom (AA) scheme, the experimental solubility value is maintained. Graphical abstract The reduction in bond distance of the 1-propanol molecule does not modify the structure factor.

Reversible assembly of silica nanoparticles at water-hydrocarbon interfaces controlled by SDS surfactant

Achieving reversible and tunable assembly of silica nanoparticles at liquid-liquid interfaces is vital for a wide range of scientific and technological applications including sustainable subsurface energy applications, catalysis, drug delivery and material synthesis. In this study, we report the mechanisms controlling the assembly of silica nanoparticles (dia. 50 nm and 100 nm) at water-heptane and water-toluene interfaces using sodium dodecyl sulfate (SDS) surfactant with concentrations ranging from 0.001-0.1 wt% using operando ultrasmall/small-angle X-ray scattering, cryogenic scanning electron microscopy imaging and classical molecular dynamics simulations. The results show that the assembly of silica nanoparticles at water-hydrocarbon interfaces can be tuned by controlling the concentrations of SDS. Silica nanoparticles are found to: (a) dominate the interfaces in the absence of interfacial SDS molecules, (b) coexist with SDS at the interfaces at low surfactant concentration of 0.001 wt% and (c) migrate toward the aqueous phase at a high SDS concentration of 0.1 wt%. Energetic analyses suggest that the van der Waals and electrostatic interactions between silica nanoparticles and SDS surfactants increase with SDS concentration. However, the favorable van der Waals and electrostatic interactions between the silica nanoparticles and toluene or heptane decrease with increasing SDS concentration. As a result, the silica nanoparticles migrate away from the water-hydrocarbon interface and towards bulk water at higher SDS concentrations. These calibrated investigations reveal the mechanistic basis for tuning silica nanoparticle assembly at complex interfaces.

Self-Assembly of NaOL-DDA Mixtures in Aqueous Solution: A Molecular Dynamics Simulation Study

The self-assembly behaviors of sodium oleate (NaOL), dodecylamine (DDA), and their mixtures in aqueous solution were systematically investigated by large-scale molecular dynamics simulations, respectively. The interaction mechanisms between the surfactants, as well as the surfactants and solvent, were revealed via the radial distribution function (RDF), cluster size, solvent-accessible surface area (SASA), hydrogen bond, and non-bond interaction energy. Results showed that the molecules more easily formed aggregates in mixed systems compared to pure systems, indicating higher surface activity. The SASA values of DDA and NaOL decreased significantly after mixing, indicating a tighter aggregation of the mixed surfactants. The RDF results indicated that DDA and NaOL strongly interacted with each other, especially in the mixed system with a 1:1 molar ratio. Compared to van der Waals interactions, electrostatic interactions between the surfactant molecules were the main contributors to the improved aggregation in the mixed systems. Besides, hydrogen bonds were found between NaOL and DDA in the mixed systems. Therefore, the aggregates in the mixed systems were much more compact in comparison with pure systems, which contributed to the reduction of the repulsive force between same molecules. These findings indicated that the mixed NaOL/DDA surfactants had a great potential in application of mineral flotation.

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.

Stability, deformation and rupture of Janus oligomer enabled self-emulsifying water-in-oil microemulsion droplets

Microemulsions exist widely in nature, daily life and industrial manufacturing processes, including petroleum production, food processing, drug delivery, new material fabrication, sewage treatment, etc. The mechanical properties of microemulsion droplets and a correlation to their molecular structures are of vital importance to those applications. Despite studies on their physicochemical determinants, there are lots of challenges of exploring the mechanical properties of microemulsions by experimental studies. Herein, atomistic modelling was utilized to study the stability, deformation, and rupture of Janus oligomer enabled water-in-oil microemulsion droplets, aiming at revealing their intrinsic relationship with Janus oligomer based surfactants and oil structures. The self-emulsifying process from a water, oil and surfactant mixture to a single microemulsion droplet was modulated by the amphiphilicity and structure of the surfactants. Four microemulsion systems with an interfacial thickness in the range of 7.4-17.3 Å were self-assembled to explore the effect of the surfactant on the droplet morphology. By applying counter forces on the water core and the surfactant shell, the mechanical stability of the microemulsion droplets was probed at different ambient temperatures. A strengthening response and a softening regime before and after a temperature-dependent peak force were identified followed by the final rupture. This work demonstrates a practical strategy to precisely tune the mechanical properties of a single microemulsion droplet, which can be applied in the formation, de-emulsification, and design of microemulsions in oil recovery and production, drug delivery and many other applications.

Molecular Dynamics Simulation Insight into Interfacial Stability and Fluidity Properties of Microemulsions

Although the interfacial properties of microemulsions have been extensively studied in both experimental and simulation research studies, the molecular mechanisms of stability and fluidity about microemulsion are still poorly understood. Herein, we report a molecular dynamics simulation study to elaborate the motion of an emulsion droplet involving dichain surfactant Aerosol OT (AOT) and its dynamics evolution at the oil-water interface. By varying the concentrations of AOT, we show that the interfacial thickness and emulsification rate display a piecewise change as the interfacial coverage increases and the W/O emulsion is more stable than the O/W one while O/W emulsion presents better fluidity. In addition, the dispersed system combined with water/AOT/n-heptane tends to form a W/O microemulsion instead of an O/W microemulsion due to the structural collapse of the latter. This work provides a molecular understanding of microemulsion interfacial stability and fluidity.

Surfactant-enhanced heterogeneity of the aqueous interface drives water extraction into organic solvents

Liquid/liquid extraction (LLE) is one of the most industrially relevant separations methods. Adsorbed surfactant is demonstrated to enhance interfacial heterogeneity and lead to water protrusions that form the basis for transport into the organic phase.

Mesoscopic simulations of temperature-dependent anchoring and wetting behavior at aqueous-liquid crystal interfaces in the presence of a rod-coil amphiphilic monolayer

Dissipative particle dynamics simulations have been applied to study the temperature dependent anchoring and wetting behavior of thermotropic liquid crystals (LCs) in the presence of a rod-coil amphiphilic monolayer at the aqueous-LC interface. Upon cooling in the nematic phase, a thermally-induced anchoring transition from homeotropic through tilted to planar has been observed. The growth and propagation of smectic order from the interfaces to the bulk nematic LCs are demonstrated to be mainly responsible for this novel transition sequence. In particular, when a complete smectic layer in the amphiphile monolayer is induced around the bulk transition of nematic-smectic-A, the propagation of homeotropic alignment fails instantly and a unique planar anchoring configuration is formed instead. While heating towards the isotropic phase, simulation results show that the nematic-isotropic transition of confined LCs is slightly shifted to a higher temperature, and a nematic wetting layer with homeotropic alignment appears in the rod block monolayer when the bulk LCs is isotropic. Our systematic simulations throughout the whole phase regimes of LCs provide important molecular-level insight into how the coupling between the ordering of LCs and adsorbents and their temperature dependencies affect the anchoring behavior in this complex system, which should be instrumental in the rational design and application of advanced LC-based biosensors with optimal operating temperature range.

Creating two self-assembly micro-environments to achieve supercrystals with dual structures using polyhedral nanoparticles

Organizing nanoparticles into supercrystals comprising multiple structures remains challenging. Here, we achieve one assembly with dual structures for Ag polyhedral building blocks, comprising truncated cubes, cuboctahedra, truncated octahedra, and octahedra. We create two micro-environments in a solvent evaporation-driven assembly system: one at the drying front and one at the air/water interface. Dynamic solvent flow concentrates the polyhedra at the drying front, generating hard particle behaviors and leading to morphology-dependent densest-packed bulk supercrystals. In addition, monolayers of nanoparticles adsorb at the air/liquid interface to minimize the air/liquid interfacial energy. Subsequent solvent evaporation gives rise to various structurally diverse dual-structure supercrystals. The topmost monolayers feature distinct open crystal structures with significantly lower packing densities than their densest-packed supercrystals. We further highlight a 3.3-fold synergistic enhancement of surface-enhanced Raman scattering efficiency arising from these dual-structure supercrystals as compared to a uniform one.

Crystals with multiple structures often perform special functions in nature, inspiring the creation of synthetic analogues. Here, the authors subject polyhedral nanoparticles to two self-assembly micro-environments to realize supercrystals with dual structures, in which the order of the surface layer differs from the bulk structure.

Interfacial Properties and Monolayer Collapse of Alkyl Benzenesulfonate Surfactant Monolayers at the Decane-Water Interface from Molecular Dynamics Simulations

The molecular structure of a surfactant molecule is known to have a great effect on the interfacial properties and the type of nanostructures formed. In this work, we have performed molecular dynamics simulations on six isomers of an alkyl benzenesulfonate surfactant to investigate the effect of the degree and position of aromatic substitution on the interfacial properties and on the collapse of the surfactant monolayer at a decane-water interface. The surface pressure of the monolayers was shown to increase with increasing surface coverage, until some of the monolayers become mechanically unstable and form large undulations. Shifting the primary alkyl chain of the surfactant from the para to the meta position was found to significantly affect the orientation of the surfactant head groups, while the attachment position of the benzene ring along the primary alkyl chain plays a greater role in the orientation of the surfactant tails. In general, to the extent considered in this work, our results suggest that additional alkyl substitution and meta substitution of the primary alkyl chains increase both the effectiveness and efficiency of the surfactants, and accelerate the onset of monolayer collapse. The interface was found to consist of an inner Helmholtz layer of partially dehydrated counterions in contact with the surfactant head groups, an outer Helmholtz layer of hydrated counterions, and a diffuse layer. The di- and trisubstituted surfactants formed nearly spherical swollen micelles encapsulating pure decane, which effectively solubilizes decane in water as a microemulsion. The monosubstituted surfactants formed elongated buds that protrude from the interface, but did not detach from the monolayer. To our knowledge, the role of aromatic substitution on interfacial properties has not been investigated by molecular simulations previously. The results from this work could provide insights to design improved surfactants by exploiting aromatic substitution to encapsulate material for drug delivery and other applications.

The structure and dynamics of Nano Particles encapsulated by the SDS monolayer collapse at the water/TCE interface

The super-saturated surfactant monolayer collapses with the nanoparticles (NPs) at the water/trichloroethylene (TCE) interface are investigated using molecular dynamics (MD) simulations. The results show that sodium alkyl sulfate (SDS) monolayer collapse is initiated by buckling and followed primarily by budding and the bud encapsulating the NPs and oil molecules. The developed bud detaches from the monolayer into a water phase and forms the swollen micelle emulsion with NPs and oil molecules. We investigate the wavelength of the initial budding and the theoretical description of the budding process. The wavelength of the monolayer increases with bending modulus. The energy barrier of the budding can be easily overcome by thermal fluctuation energy, which indicates that budding process proceeds rapidly.

Surface Electrostatic Potential and Water Orientation in the presence of Sodium Octanoate Dilute Monolayers Studied by Means of Molecular Dynamics Simulations

A series of atomistic molecular dynamics simulations were performed in the present investigation to assess the spontaneous formation of surfactant monolayers of sodium octanoate at the water-vacuum interface. The surfactant surface coverage increased until a saturation threshold was achieved, after which any further surfactant addition led to the formation of micellar aggregates within the solution. The saturated films were not densely packed, as might be expected for short-chained surfactants, and all films regardless of the surface coverage presented surfactant molecules with the same ordering pattern, namely, with the ionic heads toward the aqueous solution and the tails lying nearly parallel to the interface. The major contributions to the electrostatic surface potential came from the charged heads and the counterion distribution, which nearly canceled out each other. The balance between the oppositely charged ions rendered the electrostatic contributions from water meaningful, amounting to ca. 10% of the contributions arising from the ionic species. And even the aliphatic tails, whose atoms bear relatively small partial atomic charges as compared to the polar molecules and molecular fragments, contributed with ca. 20% of the total electrostatic surface potential of the systems under investigation. Although the aliphatic tails were not so orderly arranged as in a compact film, the C-H bonds assumed a preferential orientation, leading to an increased contribution to the electrostatic properties of the interface. The most prominent feature arising from the partitioning of the electrostatic potential into individual contributions was the long-range ordering of the water molecules. This ordering of the water molecules produced a repulsive dipole-dipole interaction between the two interfaces, which increased with the surface coverage. Only for a water layer wider than 10 nm was true bulk behavior observed, and the repulsive dipole-dipole interaction faded away.

Transferability of the coarse-grained potentials for trans-1,4-polybutadiene

The CG torsion potential has a dual effect on the reproduction of crystallization of trans-1,4-polybutadiene in CG simulation.

Determining the structure of interfacial peptide films: comparing neutron reflectometry and molecular dynamics simulations

The peptides AM1 and Lac21E self-organize into switchable films at an air-water interface. In an earlier study, it was proposed that both AM1 and Lac21E formed monolayers of α-helical peptides based on consistency with neutron reflectivity data. In this article, molecular dynamics simulations of assemblies of helical and nonhelical AM1 and Lac21E at an air-water interface suggest some tendency for the peptides to spontaneously adopt an α-helical conformation. However, irrespective of the structure of the peptides, the simulations reproduced not only the structural properties of the films (thickness and distribution of the hydrophobic and hydrophilic amino acids) but also the experimental neutron reflectivity measurements at different contrast variations. This suggests that neutron reflectometry alone cannot be used to determine the structure of the peptides in this case. However, together with molecular dynamics simulations, it is possible to obtain a detailed understanding of peptide films at an atomic level.

Dynamic morphologies of microscale droplet interface bilayers

Droplet interface bilayers (DIBs) are a powerful platform for studying the dynamics of synthetic cellular membranes; however, very little has been done to exploit the unique dynamical features of DIBs. Here, we generate microscale droplet interface bilayers (μDIBs) by bringing together femtoliter-volume water droplets in a microfluidic oil channel, and characterize morphological changes of the μDIBs as the droplets shrink due to evaporation. By varying the initial conditions of the system, we identify three distinct classes of dynamic morphology. (1) Buckling and fission: when forming μDIBs using the lipid-out method (lipids in oil phase), lipids in the shrinking monolayers continually pair together and slide into the bilayer to conserve their mass. As the bilayer continues to grow, it becomes confined, buckles, and eventually fissions one or more vesicles. (2) Uniform shrinking: when using the lipid-in method (lipids in water phase) to form μDIBs, lipids uniformly transfer from the monolayers and bilayer into vesicles contained inside the water droplets. (3) Stretching and unzipping: finally, when the droplets are pinned to the wall(s) of the microfluidic channel, the droplets become stretched during evaporation, culminating in the unzipping of the bilayer and droplet separation. These findings offer a better understanding of the dynamics of coupled lipid interfaces.