Self-Accumulation of Aromatics at the Oil−Water Interface through Weak Hydrogen Bonding

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.

New insight into the role of the self-assembly of heteroatom compounds in heavy oil viscosity enhancement

Unveiling the role of heteroatom compounds in heavy oil viscosity is pivotal for finding targeted viscosity reduction methods to improve oil recovery. This research investigates the impact of heteroatoms in asphaltene molecules by utilizing quantum chemical calculations and molecular dynamics simulations to analyze their electrostatic potential characteristics, pairwise interactions, and dynamic behavior within realistic reservoirs. Heteroatom compounds can influence the molecular-level properties of asphaltenes and thus impact the macroscopic behavior of heavy oils. Research results suggest that the presence of ketone and aromatic rings in asphaltene molecules leads to the unrestricted movement of pi electrons due to their collective electronegativity. Two distinct configurations of asphaltene dimers, face-to-face, and edge-to-face, were observed. Intermolecular interactions were predominantly governed by van der Waals forces, highlighting their significant role in stabilizing asphaltene aggregates. The distribution of asphaltene molecules in the oil phase can be summarized as the "rebar-cement" theory. In the heteroatom-free system, the face-to-face peaks in the radial distribution function exhibit significantly reduced magnitudes compared to those in the heteroatom-containing system. This emphasizes the pivotal function of heteroatoms in connecting molecular components to form a more compact asphaltene structure, which may result in a higher viscosity of heavy oil. These findings give insight into the significance of heteroatoms in bridging molecular components and shaping the intricate structure of asphaltene and advance our understanding of heavy oil viscosity properties.

Experiments and Molecular Simulations to Study the Effect of Surface-Active Compounds in Mixtures of Model Oils on CO2 Corrosion during Intermittent Oil-Water Wetting

Intermittent oil-water wetting can have a significant effect on the internal corrosion of steel pipelines. This paper presents a combined experimental and molecular modeling study of several influential factors on the surface properties and corrosion behavior of mild steel in CO2 environments. The influence of different model oils (LVT-200 and Aromatic-200) and select surface-active compounds (myristic acid, cyclohexane butyric acid, and oleic acid) on the corrosion behavior of carbon steel during intermittent oil-water wetting was determined by measuring the corrosion rate after intermittent wetting cycles. The interfacial tension measurements were performed to study the incorporation of the oil phase along with surface-active molecules in the protective layer formed on the specimen surface. Results showed that the interfacial tension for an aromatic oil-water interface is lower than that for an aliphatic oil-water interface. To understand this result, molecular dynamics simulations of oil-water interfaces were performed in the presence of surface-active molecules and different oils to analyze the structure of the layer formed at the interface. The simulations supported the hypothesis that aromatic molecules are less structured at the interface, which results in the incorporation of more water molecules into the protective layer formed at the steel surface, causing a higher corrosion rate. On the other hand, the simulations revealed that myristic acid in an aliphatic oil forms a well-aligned structure at the interface, devoid of any water molecules. This is in agreement with the hypothesis that the linear molecular structure of myristic acid favors the alignment of molecules at an aliphatic oil-water interface, resulting in a lower interfacial tension and more effective corrosion mitigation as compared to the other two nonlinear compounds tested. It is concluded that an important factor controlling the corrosion behavior is the molecular structure of the oil-water interface, which is adopted by the steel surface layer through the Langmuir-Blodgett process.

Computational predictions of interfacial tension, surface tension, and surfactant adsorption isotherms

All-atom (AA) molecular dynamics (MD) simulations are employed to predict interfacial tensions (IFT) and surface tensions (ST) of both ionic and non-ionic surfactants. The general AMBER force field (GAFF) and variants are examined in terms of their performance in predicting accurate IFT/ST, γ, values for chosen water models, together with the hydration free energy, ΔGhyd, and density, ρ, predictions for organic bulk phases. A strong correlation is observed between the quality of ρ and γ predictions. Based on the results, the GAFF-LIPID force field, which provides improved ρ predictions is selected for simulating surfactant tail groups. Good γ predictions are obtained with GAFF/GAFF-LIPID parameters and the TIP3P water model for IFT simulations at a water-triolein interface, and for GAFF/GAFF-LIPID parameters together with the OPC4 water model for ST simulations at a water-vacuum interface. Using a combined molecular dynamics-molecular thermodynamics theory (MD-MTT) framework, a mole fraction of C12E6 molecule of 1.477 × 10-6 (from the experimental critical micelle concentration, CMC) gives a simulated surface excess concentration, ΓMAX, of 76 C12E6 molecules at a 36 nm2 water-vacuum surface (3.5 × 10-10 mol cm-2), which corresponds to a simulated ST of 35 mN m-1. The results compare favourably with an experimental ΓMAX of C12E6 of 3.7 × 10-10 mol cm-2 (80 surfactants for a 36 nm2 surface) and experimental ST of C12E6 of 32 mN m-1 at the CMC.

Molecular Simulation of Surfactant Displacement of Residual Oil in Nanopores: Formation of Water Channels and Electrostatic Interaction

The water-oil-rock system's surfactant and electrostatic interactions are essential for removing oil droplets from rock substrates. Our work illustrates the impact of surface charge on the oil contact angle in an ideal system comprising silica, water, and dodecane; smaller contact angles are observed for more polar substrates. Modifying the polarity of the model silica surface allows for the observation of the creation of heteromolecule channels and the process of stripping crude oil while accounting for the impacts of water flow and different types of surfactant molecules. In solutions containing ionic surfactants, the injection and diffusion of water molecules between the oil layer and the silica substrate are facilitated by the disturbance of the oil molecules by the surfactant molecules. By comparing different surfactants in water flow, the characterization of water molecular channels and the stripping process of crude oil can be observed. The disruption of oil molecules by the surfactant molecules has been found to enhance the injection and diffusion of water molecules between the oil layer and the silica substrate in solutions containing ionic surfactants. The size of the contact angle and the extension of the water channel are simultaneously greatly influenced by the surfactant's molecular characteristics and the substrate's polarity. These simulation results show that several factors influence the process of water molecule channel creation that water molecules diffuse, and the detachment of oil from the silica substrate is facilitated by the migration of surfactants to the bottom of the oil molecule and the electrostatic interactions between the water molecules and the silica substrate.

Theoretical investigation of asphaltene molecules in crude oil viscoelasticity enhancement

Understanding the mechanisms of viscosity enhancement in crude oil phases is crucial for optimizing extraction and transportation processes. The enhanced viscosity mechanism of crude oil phase can be attributed to the intricate intermolecular interactions between asphaltene molecules. However, the molecular mechanism of the viscosification of asphaltene molecules in crude oil is not yet to be fully understood. In this work, molecular dynamics simulations were employed to investigate the dynamic behavior and viscosification mechanism of asphaltene molecules in complex oil phases. Research suggests that the neutral surface of asphaltenes features abundant positive and negative electrostatic potential regions, facilitating complementary pairing between these areas. This significantly augments electrostatic interactions among asphaltene molecules. Besides, the expansive nonpolar expanse on the normal asphaltene surface facilitates interactions between asphaltenes and crude oil molecules. This leads the crude oil viscosity of the system containing normal asphaltene is higher than that of the system containing acidic asphaltene under the same mass fraction (382 μ Pa·s for AAsp and 416 μ Pa·s for NAsp). This work provides insight into the viscosity enhancement mechanisms in crude oil phases and is helpful in improving the efficiency of crude oil extraction and transportation.

Contributions of London Dispersion Forces to Solution-Phase Association Processes

ConspectusDespite their ubiquity and early discovery, London dispersion forces are often overlooked. This is due, in part, to the difficulty in assessing their contributions to molecular and polymeric structure, stability, properties, and reactivities. However, recent advances in modeling have revealed that dispersion interactions play an important role in many important chemical and biological processes. Experimental confirmation of their impact in solution has been challenging, leading to controversies about their relative importance.In the course of studying noncovalent interactions using molecular devices, our understanding and appreciation for the importance of dispersion interactions have evolved. This Account follows this intellectual journey by using examples from the literature. The goals are twofold: to describe recent advances in understanding the interaction and to provide guidance to researchers studying weak noncovalent interactions. However, first, the experimental methods for measuring the effects of dispersion interactions and the strategies for isolating their influence are described. These include the design of molecular devices to measure these weak noncovalent interactions and the strategies to disentangle the solvation, solvophobic, and dispersion components of the resulting equilibria.The literature examples are organized around five fundamental questions. (1) Do dispersion interactions have a measurable effect on solution equilibria? (2) To what extent do solvents attenuate or compensate for dispersion interactions? (3) To what extent do the solvation and solvophobic terms influence the dispersion equilibria? (4) Can we predict whether a system will form attractive dispersion or repulsive steric interactions? (5) Can the dispersion term be isolated and interrogated? We were often surprised by the answers to these questions. In each case, we describe how the systems were designed to address these questions and discuss possible interpretations of the results.While dispersion interactions in solution were weak (usually <1 kcal/mol), their influence on complexation and conformational equilibria can be observed and measured. This underscores the significance of these interactions in molecular recognition, coordination chemistry, reaction design, and catalysis. The solvent components of the dispersion equilibria can also be significant. Therefore, the isolation of the dispersion contributions from the solvation and solvophobic effects represents an ongoing challenge. The experimental studies also provide important benchmarks and offer valuable insights to help refine the next generation of computational solvent models.

The Solvation Shell of Small Solutes in Aqueous-Organic Solvent Mixtures and Its Implications for Reversed-Phase Liquid Chromatography

Reversed-phase liquid chromatography (RPLC) operates with water-organic solvent (W-OS) mobile phases where preferential solvation (PS) of solutes is likely. To investigate the relevance of the solute solvation shell in the mobile phase for RPLC retention, we combine data from molecular dynamics simulations of small, neutral solutes (six analytes and two dead time markers) in W-methanol (MeOH) and W-acetonitrile (ACN) mixtures with corresponding retention data obtained on an RPLC column over a wide range of W/OS ratios. Data derived from Kirkwood-Buff integrals show PS by the OS for analytes vs low or negative PS for dead time markers. W-ACN mixtures generate a higher amount of PS than W-MeOH mixtures, which contributes to the higher eluent strength of ACN in RPLC. Difference spatial distribution functions reveal anisotropic solvation shells with OS excess at hydrocarbon elements and W excess at functional groups, predicting that retention by the hydrophobic stationary phase is favored by hydrocarbon elements and limited by functional groups. Analysis of solute-solvent hydrogen bonds pinpoints the hydrogen-bond requirements toward W as the retention-limiting factor. The relation between the solute solvation shell and retention confirms the importance of W-OS and solute-W hydrogen bonding for RPLC retention.

Synergistic Interaction of Hyperbranched Polyglycerols and Cetyltrimethylammonium Bromide for Oil/Water Interfacial Tension Reduction: A Molecular Dynamics Study

Applying surfactants to reduce the interfacial tension (IFT) on water/oil interfaces is a proven technique. The search for new surfactants and delivery strategies is an ongoing research area with applications in many fields such as drug delivery through nanoemulsions and enhanced oil recovery. Experimentally, the combination of hyperbranched polyglycerol (HPG) with cetyltrimethylammonium bromide (CTAB) substantially reduced the observed IFT of oil/water interface, 0.9 mN/m, while HPG alone was 5.80 mN/m and CTAB alone IFT was 8.08 mN/m. Previous simulations in an aqueous solution showed that HPG is a surfactant carrier. Complementarily, in this work, we performed classical molecular dynamics simulations on combinations of CTAB and HPG with one aliphatic chain to investigate further the interaction of this pair in oil interfaces and propose the mechanism of IFT decrease. Basically, from our results, one can observe that the IFT reduction comes from a combination of effects that have not been observed for other dual systems: (i) Due to the CTAB-HPG strong interaction, a weakening of their specific and isolated interactions with the water and oil phases occurs. (ii) Aggregates enlarge the interfacial area, turning it into a less ordered interface. (iii) The spread of individual molecules charge profiles leads to the much lower interfacial tension observed with the CTAB+HPG systems.

Study on the Behavior of Saturated Cardanol-Based Surfactants at the Crude Oil/Water Interface through Molecular Dynamics Simulations

Cardanol is a green biosurfactant with broad application prospects, which is expected to be used to enhance oil recovery (EOR). This paper designed two types of surfactants (extended and nonextended), including six kinds of nonionic and anion-nonionic surfactants. The position changes of PO and EO chains and the effects of different hydrophilic groups on the interface properties were studied with molecular dynamics simulations by constructing a model of crude oil (containing four components) and water molecules. The results of interfacial tension and solvent-accessible surface area showed that the interfacial properties of sulfate were better than those of sulfonates and nonionic surfactants. Meanwhile, the interface properties of nonextended surfactants were better than those of extended surfactants. The gyration radius (Rg) and tilt angle data demonstrated that when EO chains were located between hydrophobic groups and PO chains (nonextended surfactants), the adsorption capacity of surfactants at crude oil and water interfaces could be effectively improved. The radial distribution function of the hydrophilic group and hydrophobic group of surfactants with water molecules and four components of the crude oil molecule, respectively, explained that surfactants (8EO8POSO4) had better emulsification performance when the intermolecular interactions between crude oil and water two phases were relatively balanced. This study provides a theoretical reference for the design of oil-displacement surfactants and the mechanism analysis of emulsification properties.

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.

Surface-bubble-modulated liquid chromatography: an experimental strategy for identification of molecular processes of solute retention in reversed-phase separation systems

Molecular level understanding of the chemistry at the aqueous/hydrophobe interface is crucial to separation processes in aqueous media, such as reversed-phase liquid chromatography (RPLC) and solid-phase extraction (SPE). Despite significant advances in our knowledge of the solute retention mechanism in these reversed-phase systems, direct observation of the behavior of molecules and ions at the interface in reversed-phase systems still remains a major challenge and experimental probing techniques that provide the spatial information of the distribution of molecules and ions are required. This review addresses surface-bubble-modulated liquid chromatography (SBMLC), which has a stationary gas phase in a column packed with hydrophobic porous materials and enables one to observe the molecular distribution in the heterogeneous reversed-phase systems consisting of the bulk liquid phase, the interfacial liquid layer, and the hydrophobic materials. The distribution coefficients of organic compounds referring to their accumulations onto the interface of alkyl- and phenyl-hexyl-bonded silica particles exposed to water or acetonitrile-water and into the bonded layers from the bulk liquid phase are determined by SBMLC. The experimental data obtained by SBMLC show that the water/hydrophobe interface exhibits an accumulation selectivity for organic compounds, which is quite different from that of the interior of the bonded chain layer, and the overall separation selectivity of the reversed-phase systems is determined by the relative sizes of the aqueous/hydrophobe interface and the hydrophobe. The solvent composition and the thickness of the interfacial liquid layer formed on octadecyl-bonded (C₁₈) silica surfaces are also estimated from the bulk liquid phase volume determined by the ion partition method employing small inorganic ions as probes. It is clarified that various hydrophilic organic compounds as well as inorganic ions recognize the interfacial liquid layer formed on the C₁₈-bonded silica surfaces as being different from the bulk liquid phase. The behavior of some solute compounds exhibiting substantially weak retention in RPLC or the so-called negative adsorption, such as urea, sugars, and inorganic ions, can rationally be interpreted with a partition between the bulk liquid phase and the interfacial liquid layer. The spatial distribution of solute molecules and the structural properties of the solvent layer on the C₁₈-bonded layer determined by the liquid chromatographic methods are discussed in comparison to the results obtained by other research groups using molecular simulation methods.

Tunable Hierarchically Structured Meso-Macroporous Carbon Spheres from a Solvent-Mediated Polymerization-Induced Self-Assembly

Herein, we report a versatile solvent-mediated polymerization-induced self-assembly (PISA) strategy to directly synthesize highly N-doped hierarchically porous structured carbon spheres with a tunable meso-macroporous configuration. The introduction of intermolecular hydrogen bonds is verified to enhance the interfacial interactions between block copolymers, oil droplets, and polyphenols. Moreover, the dominant hydrogen-bond-driven interactions can be systematically manipulated by selecting different cosolvent systems to generate diverse porous structures from the transformation of micellar and precursor co-assembly. Impressively, hierarchically structured meso-macroporous N-doped carbon spheres present simultaneously tunable sphere sizes and mesopores and macropores, ranging from 1.2 μm, 9/50 and 227 nm to 1.0 μm, 40, and 183 nm and 480, 24, and 95 nm. As a demonstration, dendritic-like N-doped hierarchically meso-macroporous carbon spheres manifest excellent enzyme-like activity, which is attributed to the continuous mass transport from the multiordered porosity. The current study provides a new platform for the synthesis of novel well-defined porous materials.

Hydrophilic Silica Nanoparticles in O/W Emulsion: Insights from Molecular Dynamics Simulation

Previous studies have been carried out on the effect of silica nanoparticles (SNPs) on the stability of oil-water emulsions. However, the combining configuration of SNPs and oil droplets at the molecular level and the effect of SNP content on the coalescence behavior of oil droplets cannot be obtained through experiments. In this paper, molecular dynamics (MD) simulation was performed to investigate the adsorption configuration of hydrophilic SNPs in an O/W emulsion system, and the effect of adsorption of SNPs on coalescence of oil droplets. The simulation results showed: (i) SNPs adsorbed on the surface of oil droplets, and excessive SNPs self-aggregated and connected by hydrogen bonds. (ii) Partially hydrophilic asphaltene and resin molecules formed adsorption configurations with SNPs, which changed the distribution of oil droplet components. Furthermore, compared with hydrophobic asphaltene, the hydrophilic asphaltene was easier to combine with SNPs. (iii) SNPs would extend the oil droplet coalescence time, and the π-π stacking structures were formed between asphaltene and asphaltene or resin molecules to enhance the connection between oil droplets during the oil droplet contact process. (iv) Enough SNPs tightly wrapped around the oil droplet, similar to the formation of a rigid film on the surface of an oil droplet, which hindered the contact and coalescence of components between oil droplets.

Synthesis of Polycarboxylate Viscosity Reducer and the Effect of Different Chain Lengths of Polyether on Viscosity Reduction of Heavy Oil

Since there are not many studies on the application of polymeric surfactants in viscosity reduction emulsification of heavy oil, a series of polyether carboxylic acid–sulfonate polymeric surfactants were synthesized. The viscosity reduction performance and the effect of different chain lengths on the viscosity reduction effect were also investigated. The viscosity reduction, emulsification, wetting, and foaming performance tests showed that the viscosity reduction performance of this series of polymeric surfactants was excellent, with the viscosity reduction rate exceeding 95%, and the viscosity was reduced to 97 mPa·s by the polymeric surfactant with a molecular weight of 600 polyethers. It was also concluded that among the three surfactants with different side chains, the polymeric surfactant with a polyether molecular weight of 600, which is the medium side-chain length, had the best viscosity reduction performance. The study showed that the polyether carboxylic acid–sulfonate polymer surfactant had a promising application in the viscosity reduction of heavy oil.

Molecular Dynamics Studies on Effective Surface-Active Additives: Toward Hard Water-Resistant Chemical Flooding for Enhanced Oil Recovery

Divalent ions, which are omnipresent in brine, may be detrimental to surfactant functionalities during chemical flooding in the enhanced oil recovery (EOR) process. Surfactant blending is one potential solution to overcome such an adverse effect. Herein, we report a molecular dynamics (MD) study to investigate the molecular arrangement and possible applications of surfactant blending in hard water-resistant chemical flooding for oil recovery. We chose commonly used anionic surfactants, sodium dodecyl sulfate (SDS), as primary surfactants. The non-ionic (propanol) and cationic [cetrimonium bromide (CTAB)] surfactants with a wide range of concentrations are introduced to the primary system. We demonstrate that CTAB can disaggregate the cation bridging when their concentration is above a certain threshold. This threshold value is related to the surfactant and cosurfactant surface charge in the interface region. The cation bridging density is maintained at a low level when the sum of surfactants and cosurfactant interface charges is neutral or positive. On the other hand, propanol barely disaggregates the cation bridging. When propanol concentration is above a certain value, it even facilitates the cation bridging formation. Both propanol and CTAB can further decrease the oil-brine interfacial tension (IFT) while having different efficacies (IFT decrement rate is different as their interface concentration increases). More rapid IFT decrement is observed when cation bridging is disaggregated (i.e., in systems with high CTAB concentrations). Increasing propanol concentration barely affects hydrogen bond (H-bond) formation between SDS and H₂O because of a low propanol distribution around SDS. On the other hand, the first increasing and then decreasing trend in H-bond density between SDS and H₂O is observed as CTAB concentration increases. Our work should provide important insights into designing chemical formulas in chemical flooding applications.

Molecular Dynamics Simulation for the Demulsification of O/W Emulsion under Pulsed Electric Field

A bidirectional pulsed electric field (BPEF) method is considered a simple and novel technique to demulsify O/W emulsions. In this paper, molecular dynamics simulation was used to investigate the transformation and aggregation behavior of oil droplets in O/W emulsion under BPEF. Then, the effect of surfactant (sodium dodecyl sulfate, SDS) on the demulsification of O/W emulsion was investigated. The simulation results showed that the oil droplets transformed and moved along the direction of the electric field. SDS molecules can shorten the aggregation time of oil droplets in O/W emulsion. The electrostatic potential distribution on the surface of the oil droplet, the elongation length of the oil droplets, and the mean square displacement (MSD) of SDS and asphaltene molecules under an electric field were calculated to explain the aggregation of oil droplets under the simulated pulsed electric field. The simulation also showed that the two oil droplets with opposite charges have no obvious effect on the aggregation of the oil droplets. However, van der Waals interactions between oil droplets was the main factor in the aggregation.

Characterizing Anion Adsorption to Aqueous Interfaces: Toluene-Water versus Air-Water

We continue our investigation of the behavior of simple ions at aqueous interfaces, employing the combination of two surface-sensitive nonlinear spectroscopy tools, broadband deep UV electronic sum-frequency generation and UV second harmonic generation, to characterize the adsorption of thiocyanate to the interface of water with toluene─a prototypical hydrophobe. We find that both the interfacial spectrum and the Gibbs free energy of adsorption closely match results previously reported for the air-water interface. We observe no relative spectral shift in the higher-energy CTTS transition of thiocyanate, implying similar solvation environments for the two interfaces. Similarly, the Gibbs free energies of adsorption agree within error; however, we expect the respective enthalpic and entropic contributions to differ between the two interfaces, similar to our earlier findings for the air-water versus graphene-water interfaces. Further experiments and theoretical modeling are necessary to quantify the mechanistic differences.

Effect of electric field on coalescence of an oil-in-water emulsion stabilized by surfactant: a molecular dynamics study

The microscopic mechanisms of electrocoalescence of O/W emulsions stabilized by surfactant were analyzed from the electric dipole moment of the surfactant, the interaction between surfactant and oil molecules and the deformation of the surfactant.

Water/oil interfacial tension reduction - an interfacial entropy driven process

The interfacial tension (IFT) of a fluid-fluid interface plays an important role in a wide range of applications and processes. When low IFT is desired, surface active compounds (e.g. surfactants) can be added to the system. Numerous attempts have been made to relate changes in IFT arising from such compounds to the specific nature of the interface. However, the IFT is controlled by an interplay of factors such as temperature and molecular structure of surface-active compounds, which make it difficult to predict IFT as those conditions change. In this study, we present the results from molecular dynamics simulations revealing the specific role surfactants play in IFT. We find that, in addition to reducing direct contact between the two fluids, surfactants serve to increase the disorder at the interface (related to interfacial entropy) and consequently reduce the water/oil IFT, especially when surfactants are present at high surface density. Our results suggest that surfactants that yield more disordered interfacial films (e.g. with flexible and/or unsaturated tails) reduce the water/oil IFT more effectively than surfactants which yield highly ordered interfacial films. Our results shed light on some of the factors that control IFT and could have important practical implications in industrial applications such as the design of cosmetics, food products, and detergents.

Ion Valency and Concentration Effect on the Structural and Thermodynamic Properties of Brine-Decane Interfaces with Anionic Surfactant (SDS)

Salt ion valency and concentration vary in actual oil reservoirs, which play an important role in the functionalities of surfactant formula during chemical flooding processes to enhance oil recovery. Herein, we report a molecular dynamics (MD) study to investigate the ion valency and concentration effect on the structural and thermodynamic properties of brine-decane interfaces with anionic surfactant (SDS), under typical reservoir conditions (353 K and 200 bar). We use two different cations (Na⁺ and Ca²⁺) and a wide range of ion concentrations (up to 3.96 M) to simulate reservoir conditions. We find that ion valency has a significant effect on the molecular configurations, which further influences the thermodynamic properties. Ca²⁺ ions can have a strong adsorption at the interface due to the strong electrostatic interactions between Ca²⁺ ions and SDS, which also results in the Cl^- ion enrichment at the interface. Furthermore, Ca²⁺ ions can form pentagon-like SDS-Ca²⁺ complexes through SDS-Ca²⁺-SDS cation bridging, which renders a nonuniform distribution of SDS at the interface. On the other hand, the cation bridging density monotonically increases as ion concentration increases for the systems without Ca²⁺ ions, while first increases, then decreases for the systems with Ca²⁺ ions. This is because the accumulation of Cl^- ions at the interface at high salt concentrations can melt SDS-Ca²⁺ complexes. This work should provide new insights into the structural and thermodynamic properties of brine-oil interfaces with an anionic surfactant, which can facilitate the optimization of chemical flooding processes.

Determination of Minimum Miscibility Pressure of CO2-Oil System: A Molecular Dynamics Study

CO₂ enhanced oil recovery (CO₂-EOR) has become significantly crucial to the petroleum industry, in particular, CO₂ miscible flooding can greatly improve the efficiency of EOR. Minimum miscibility pressure (MMP) is a vital factor affecting CO₂ flooding, which determines the yield and economic benefit of oil recovery. Therefore, it is important to predict this property for a successful field development plan. In this study, a novel model based on molecular dynamics to determine MMP was developed. The model characterized a miscible state by calculating the ratio of CO₂ and crude oil atoms that pass through the initial interface. The whole process was not affected by other external objective factors. We compared our model with several famous empirical correlations, and obtained satisfactory results—the relative errors were 8.53% and 13.71% for the two equations derived from our model. Furthermore, we found the MMPs predicted by different reference materials (i.e., CO₂/crude oil) were approximately linear (R² = 0.955). We also confirmed the linear relationship between MMP and reservoir temperature (T_R). The correlation coefficient was about 0.15 MPa/K in the present study.

Calculation of Solid-Fluid Interfacial Free Energy with Consideration of Solid Deformation by Molecular Dynamics Simulations

Fluid-fluid interfacial free energy can be measured accurately and can also be calculated from molecular simulations. However, it is challenging to measure solid-fluid interfacial free energy directly. Accurate computation has not yet been advanced by molecular simulations. In this study, we derive working expressions for estimating solid-fluid interfacial free energy based on the free-energy perturbation method with consideration of solid deformation. A Lennard-Jones solid-fluid system is simulated. Our derivations indicate that the effect of solid deformation is pronounced on solid-fluid interfacial free energy, and the results may be significantly different from the conventional test area method. Our results reveal that the contribution of the solid deformation highly depends on the stress conditions in the solid, which can be either positive or negative. Adsorption of fluids onto the solid surface has a significant effect on interfacial free energy. In weak adsorption, the interfacial free energy is close to the solid-vacuum surface free energy. Strong adsorption results in a significant reduction in interfacial free energy.

Water as the reaction medium in organic chemistry: from our worst enemy to our best friend

A review presenting water as the logical reaction medium for the future of organic chemistry. A discussion is offered that covers both the "on water" and "in water" phenomena, and how water is playing unique roles in each, specifically with regard to its use in organic synthesis.

Intrinsic difference between phenyl hexyl- and octadecyl-bonded silicas in the solute retention selectivity in reversed-phase liquid chromatography with aqueous mobile phase

For choosing an optimal column for a particular separation by reversed-phase liquid chromatography (RPLC), it is essential to quantitatively understand the effects of the chemical structure of hydrophobic bonded layer derived onto silica particles on the distribution equilibrium of a solute compound at the interface between the aqueous mobile phase and the packing material. However, there is still a lack of understanding of the solute distribution equilibrium in RPLC separation due to the complexities of the chemistry at the interface between the mobile phase and the bonded layer. We successfully determined the distribution coefficients of various organic compounds concerning to their accumulation onto the water/bonded layer interface and into the bonded layer from bulk water using surface-bubble-modulated liquid chromatography with octadecyl- and phenyl hexyl-bonded silica columns. The water/phenyl hexyl-bonded layer interface accumulates organic compounds much less than the water/octadecyl-bonded layer interface due to its lower interfacial tension, and this result suggests that phenyl hexyl group orient their benzene ring facing toward water. On the other hand, aromatic moiety of phenyl hexyl group enhances partitioning of the organic compounds into the bonded layer. Experimental findings in the present work demonstrated that the water/bonded layer interface and the bonded layer itself have independent contributions to the solute distribution and the water/phenyl hexyl-bonded layer interface shows quite different solute retention selectivity from the water/octadecyl-bonded layer interface.

Advances and Opportunities of Oil-in-Oil Emulsions

Emulsions are mixtures of two immiscible liquids in which droplets of one are dispersed in a continuous phase of the other. The most common emulsions are oil-water systems, which have found widespread use across a number of industries, for example, in the cosmetic and food industries, and are also of advanced scientific interest. In addition, the past decade has seen a significant increase in both the design and application of nonaqueous emulsions. This has been primarily driven by developments in understanding the mechanism of effective stabilization of oil-in-oil (o/o) systems, either using block copolymers (BCPs) or solid (Pickering) particles with appropriate surface functionality. These systems, as highlighted in this review, have enabled emergent applications in areas such as pharmaceutical delivery, energy storage, and materials design (e.g., polymerization, monolith, and porous polymer synthesis). These o/o emulsions complement traditional emulsions that utilize an aqueous phase and allow the use of materials incompatible with water. We assess recent advances in the preparation and stabilization of o/o emulsions, focusing on the identity of the stabilizer (BCP or particle), the interplay between stabilizer and oils, and highlighting applications and opportunities associated with o/o emulsions.

Transport Properties of Waxy Crude Oil: A Molecular Dynamics Simulation Study

To study the effects of paraffin on viscosity of waxy crude oil and transport properties of small molecules, light and waxy crude oil models were investigated at atmospheric pressure and 293-323 K temperature range using atomistic molecular dynamics simulations. The optimized parameters for liquid simulations all-atom (OPLS-AA) and atomistic polarizable potential for liquids, electrolytes, and polymers (APPLE&P) force fields were employed. The self-diffusion coefficients, viscosity, and paraffin configurations were compared for two oil models and between the two employed force fields. However, the behavior of paraffin molecules predicted by two force fields was quite different. Simulations using the OPLS-AA force field showed crystallization of longer paraffin molecules below 323 K, while simulations with the APPLE&P force field demonstrated a homogeneous mixture down to 293 K. To provide additional validation of the employed force fields, the density, diffusion coefficient, and crystallization of pure alkanes were compared with experimental data. The density and diffusion coefficients of n-C₆ and n-C₁₄ simulated with the APPLE&P force field were found to be in much closer agreement with the experimental data. The OPLS-AA force field was found to overestimate the crystallization temperature of pure alkanes. Therefore, simulations with the APPLE&P provide more realistic description of the waxy oil structure and transport properties. In this temperature range, the paraffin molecules are homogeneously distributed in the mixture, and viscosity of the system increased by a factor of two compared to light oil. Crystallization of paraffins requires lower temperatures or/and the presence of other components such as nanoparticles or asphaltene molecules to facilitate nucleation.

Emulsification of Surfactant on Oil Droplets by Molecular Dynamics Simulation

Heavy oil in crude oil flooding is extremely difficult to extract due to its high viscosity and poor fluidity. In this paper, molecular dynamics simulation was used to study the emulsification behavior of sodium dodecyl sulfonate (SDSn) micelles on heavy oil droplets composed of asphaltenes (ASP) at the molecular level. Some analyzed techniques were used including root mean square displacement, hydrophile-hydrophobic area of an oil droplet, potential of mean force, and the number of hydrogen bonds between oil droplet and water phase. The simulated results showed that the asphaltene with carboxylate groups significantly enhances the hydration layer on the surface of oil droplets, and SDSn molecules can change the strength of the hydration layer around the surface of the oil droplets. The water bridge structure between both polar heads of the surfactant was commonly formed around the hydration layer of the emulsified oil droplet. During the emulsification of heavy oil, the ratio of hydrophilic hydrophobic surface area around an oil droplet is essential. Molecular dynamics method can be considered as a helpful tool for experimental techniques at the molecular level.

Brine-Oil Interfacial Tension Modeling: Assessment of Machine Learning Techniques Combined with Molecular Dynamics

The physical chemistry mechanisms behind the oil-brine interface phenomena are not yet fully clarified. The knowledge of the relation between brine composition and concentration for a given oil may lead to the ionic tuning of the injected solution on geochemical and enhanced oil recovery processes. Thus, it is worth examining the parameters influencing the interfacial properties. In this context, we have combined machine learning (ML) techniques with classical molecular dynamics simulations (MD) to predict oil/brine interfacial tensions (IFT) effectively and compared this process to a linear regression (LR) method. To diversify our data set, we have introduced a new atomistic crude oil model (medium) with 36 different types of hydrocarbon molecules. The MD simulations were performed for mono- and multicomponent (toluene, heptane, Heptol, light, and medium) oil systems interfaced with sulfate and chloride brines with varying cations (Na⁺, K⁺, Ca²⁺, and Mg²⁺) and salinity concentration. Thus, a consistent IFT data set was built for the ML training and LR fitting at room temperature and pressure conditions, over the feature space considering oil density, oil composition, salinity, and ionic concentrations. On the basis of gradient boosted (GB) algorithms, we have observed that the dominant quantities affecting the IFT are related to the oil attributes and the salinity concentration, and no specific ion dominates the IFT changes. When the obtained LR model was validated against MD and experimental data from the literature, the error varied up to 2% and 9%, respectively, showing a robust and consistent transferability. The combination of MD simulations and ML techniques may provide a fast and cost-effective IFT determination over multiple and complex fluid-fluid and fluid-solid interfaces.