The Water-Alkane Interface at Various NaCl Salt Concentrations: A Molecular Dynamics Study of the Readily Available Force Fields

Experimental and Molecular Dynamics Simulation to Investigate Oil Adsorption and Detachment from Sandstone/Quartz Surface by Low-Salinity Surfactant Brines

In this study, we explore the impact of monovalent (NaCl) and divalent (CaCl2) brines, coupled with sodium dodecyl sulfate (SDS) surfactant at varying low concentrations, on the detachment and displacement of oil from sandstone rock surfaces. Employing the sessile drop method and molecular dynamics simulations, we scrutinize the behavior of the brine solutions. Our findings reveal that both low salinity and low-salinity surfactant solutions induce a gradual shift in rock wettability toward a more water-wet state. This wettability transformation is not instantaneous but evolves over time, as observed through meticulous molecular motion analyses. Through contact angle measurements and molecular dynamics simulations, we delve into the molecular motion at subpore and micropore scales on sandstone/quartz surfaces. The adsorption of surface-active agents from the oil to the oil-brine interface results in a reduced interfacial tension, significantly contributing to oil displacement. Notably, low salinity concentrations ranging from 1000 to 10,000 ppm exhibit the lowest contact angles within 30 min across all solutions. However, higher concentrations deviate from this declining trend, especially with divalent ions like Ca2+, which bridge polar molecules onto the rock surface, resulting in an increased oil-wetting state. This research unveils the intricate molecular motions involved in employing low-salinity surfactant solutions for oil detachment from surfaces. Furthermore, it provides valuable insights into the underlying forces driving oil detachment and wettability alteration.

Hydrocarbon Transportation in Heterogeneous Shale Pores by Molecular Dynamic Simulation

Shale oil in China is widely distributed and has enormous resource potential. The pores of shale are at the nanoscale, and traditional research methods encounter difficulty in accurately describing the fluid flow mechanism, which has become a bottleneck restricting the industrial development of shale oil in China. To clarify the distribution and migration laws of fluid microstructure in shale nanopores, we constructed a heterogeneous inorganic composite shale model and explored the fluid behavior in different regions of heterogeneous surfaces. The results revealed the adsorption capacity for alkanes in the quartz region was stronger than that in the illite region. When the aperture was small, solid-liquid interactions dominated; as the aperture increased, the bulk fluid achieved a more uniform and higher flow rate. Under conditions of small aperture/low temperature/low pressure gradient, the quartz region maintained a negative slip boundary. Illite was more hydrophilic than quartz; when the water content was low, water molecules formed a "liquid film" on the illite surface, and the oil flux percentages in the illite and quartz regions were 87% and 99%, respectively. At 50% water content, the adsorbed water in the illite region reached saturation, the quartz region remained unsaturated, and the difference in the oil flux percentage of the two regions decreased. At 70% water content, the adsorbed water in the two regions reached a fully saturated state, and a layered structure of "water-two-phase region-water" was formed in the heterogeneous nanopore. This study is of great significance for understanding the occurrence characteristics and flow mechanism of shale oil within inorganic nanopores.

Revealing innovative JAK1 and JAK3 inhibitors: a comprehensive study utilizing QSAR, 3D-Pharmacophore screening, molecular docking, molecular dynamics, and MM/GBSA analyses

The heterocycle compounds, with their diverse functionalities, are particularly effective in inhibiting Janus kinases (JAKs). Therefore, it is crucial to identify the correlation between their complex structures and biological activities for the development of new drugs for the treatment of rheumatoid arthritis (RA) and cancer. In this study, a diverse set of 28 heterocyclic compounds selective for JAK1 and JAK3 was employed to construct quantitative structure-activity relationship (QSAR) models using multiple linear regression (MLR). Artificial neural network (ANN) models were employed in the development of QSAR models. The robustness and stability of the models were assessed through internal and external methodologies, including the domain of applicability (DoA). The molecular descriptors incorporated into the model exhibited a satisfactory correlation with the receptor-ligand complex structures of JAKs observed in X-ray crystallography, making the model interpretable and predictive. Furthermore, pharmacophore models ADRRR and ADHRR were designed for each JAK1 and JAK3, proving effective in discriminating between active compounds and decoys. Both models demonstrated good performance in identifying new compounds, with an ROC of 0.83 for the ADRRR model and an ROC of 0.75 for the ADHRR model. Using a pharmacophore model, the most promising compounds were selected based on their strong affinity compared to the most active compounds in the studied series each JAK1 and JAK3. Notably, the pharmacokinetic, physicochemical properties, and biological activities of the selected compounds (As compounds ZINC79189223 and ZINC66252348) were found to be consistent with their therapeutic effects in RA, owing to their non-toxic, cholinergic nature, absence of P-glycoprotein, high gastrointestinal absorption, and ability to penetrate the blood-brain barrier. Furthermore, ADMET properties were assessed, and molecular dynamics and MM/GBSA analysis revealed stability in these molecules.

Effects of molecular size and orientation on the interfacial properties and wetting behavior of water/n-alkane systems: a molecular-dynamics study

Molecular dynamics simulations (MD) are performed to study the interfacial structure/tension and wetting behavior of water/n-alkane systems (water/nC5 to water/nC16 where nCx = CxH(2x + 2)). In particular, we study complete-to-partial wetting transitions by changing the n-alkane chain length (N_C) at a constant temperature, T = 295 K. Simulations are carried out with a united-atom TraPPE model for n-alkanes and the TIP4P-2005 model of water. Simulation results are in excellent agreement with the initial spreading coefficients and contact angles calculated using experimental values of the surface and interfacial tensions. In addition, it has been determined that water/(nC5-nC7) and water/(nC8-nC16), respectively, exhibit complete and partial initial wetting modes. Simulations show that the interfacial structures of water/(nC5-nC7) are different from water/(nC8-nC16) systems. In the latter, water preferentially orients near the interface to increase the number of hydrogen bonds and the charge and mass densities. Moreover, the orientation of n-alkane molecules at water/(nC8-nC16) interfaces has a long-range persistence, resulting in layered structures that increase with N_C. In addition, simulation results of the orientational order parameter S_z show alignment behavior of the n-alkane molecules with respect to the interfaces. Simulations predict that the central segments of n-alkane are strongly packed in the interfaces while the end segments (methyl groups) form smaller peaks in the outer edge of the layer. This observation confirms the "horseshoe" or "C-shaped" structure of n-alkane molecules in the water/n-alkane interfaces. At constant temperature, the interface widths of both water and the n-alkanes decrease with increasing n-alkane molecular length. These results suggest that increasing the n-alkane chain length affects the water/n-alkane interfacial properties in a manner similar to that of cooling.

Quantifying the Effect of Spatial Confinement on the Diffusion and Nuclear Magnetic Resonance Relaxation of Water/Hydrocarbon Mixtures: A Molecular Dynamics Simulation Study

We implement molecular dynamics (MD) simulations to explore the relaxation mechanisms involved in the description of translational diffusion and rotational dynamics in water/hydrocarbon interfaces. The analysis of density profiles, self-diffusion coefficients, and nuclear magnetic resonance (NMR) relaxation properties as a function of the confinement layer width and type of hydrocarbons improves the understanding of confined water properties at water/oil interfaces. Density profile fluctuations reveal the presence of water-oil interactions close to the interface. MD results show that self-diffusion coefficients and NMR relaxation times of planar and cylindrical water/oil interfaces are strongly influenced by layer thickness and geometry. Shorter (between 20 and 60%) self-diffusion coefficients and ¹H NMR relaxation times were obtained for water/n-pentane, water/n-decane, and water/n-hexadecane systems than bulk diffusion coefficients. An increase in Larmor frequency from 2.3 MHz to 400 MHz shows that longitudinal relaxation time (T₁) of confined oil has slightly larger differences at higher frequencies than the transverse relaxation time (T₂). At 400 MHz, n-alkanes (n-pentane, n-decane, and n-hexadecane) exhibit longer relaxation times than at smaller frequency values (2.3 and 22 MHz). Analysis of spin-spin and spin-lattice times provides relevant information about inter- and intramolecular relaxation mechanisms of water and oil as a function of geometry and width of the interface layer. These MD results suggest that the strength of confinement and geometry play a vital role in the diffusion and NMR relaxation properties of water/oil interfaces.

Flotation surface chemistry of water-soluble salt minerals: from experimental results to new perspectives

The flotation separation of water-soluble salt minerals has to be conducted under the condition of saturation in brines which represents a challenging but exciting topic of colloid and surface chemistry. Despite several proposals on explaining the success of this industrial application for many decades, our understanding of the flotation separation is still far from complete yet, owing to the complexity of the highly selective collection of salt crystals by air bubbles in brines. Here, we thoroughly review the experimental results for halogen, oxyanion, and double salts and match them with the proposed theories on the flotation of soluble salts to identify the agreed and disagreed cases. The experimental results show that the flotation of these salts varies from collectors (surfactants applied to control the crystal hydrophobicity) to collectors and is strongly affected by the brine ion composition and pH conditions. We find some exceptional flotation results that cannot be simply explained by the crystal surface charge and wettability. Furthermore, we outline several disputes and discrepancies between the experiments and the theories when different collectors are applied. Apart from the extensive consideration of surface hydration, the presence of external ion species exhibits ubiquitous effects on the surface properties of salt crystals and the colloidal properties of collectors. We conclude that the interactions between salt ions, water molecules, collectors, and salt crystals must be considered more thoroughly, and the activity of collectors at the air-liquid interface should also be the focus. Advanced techniques such as molecular dynamics simulation, atomic force microscopy, X-ray photoelectron spectroscopy, and sum-frequency generation spectroscopy are expected to be promising research tools for future studies.

Coarse-grained molecular models of the surface of hair

We present a coarse-grained molecular model of the surface of human hair, which consists of a supported lipid monolayer, in the MARTINI framework. Using coarse-grained molecular dynamics (MD) simulations, we identify a lipid grafting distance that yields a monolayer thickness consistent with both atomistic MD simulations and experimental measurements of the hair surface. Coarse-grained models for fully-functionalised, partially damaged, and fully damaged hair surfaces are created by randomly replacing neutral thioesters with anionic sulfonate groups. This mimics the progressive removal of fatty acids from the hair surface by bleaching and leads to chemically heterogeneous surfaces. Using molecular dynamics (MD) simulations, we study the island structures formed by the lipid monolayers at different degrees of damage in vacuum and in the presence of polar (water) and non-polar (n-hexadecane) solvents. We also use MD simulations to compare the wetting behaviour of water and n-hexadecane droplets on the model surfaces through contact angle measurements, which are compared to experiments using virgin and bleached hair. The model surfaces capture the experimentally-observed transition of the hair surface from hydrophobic (and oleophilic) to hydrophilic (and oleophobic) as the level of bleaching damage increases. By selecting surfaces with specific damage ratios, we obtain contact angles from the MD simulations that are in good agreement with experiments for both solvents on virgin and bleached human hairs. To negate the possible effects of microscale curvature and roughness of real hairs on wetting, we also conduct additional experiments using biomimetic surfaces that are co-functionalised with fatty acids and sulfonate groups. In both the MD simulations and experiments, the cosine of the water contact angle increases linearly with the sulfonate group surface coverage with a similar slope. We expect that the proposed systems will be useful for future molecular dynamics simulations of the adsorption and tribological behaviour of hair, as well as other chemically heterogeneous surfaces.

Interfacial behavior of the decane + brine + surfactant system in the presence of carbon dioxide, methane, and their mixture

Molecular dynamics simulations are carried out to get insights into the interfacial behavior of the decane + brine + surfactant + CH₄ + CO₂ system at reservoir conditions. Our results show that the addition of CH₄, CO₂, and sodium dodecyl sulfate (SDS) surfactant at the interface reduces the IFTs of the decane + water and decane + brine (NaCl) systems. Here the influence of methane was found to be less pronounced than that of carbon dioxide. As expected, the addition of salt increases the IFTs of the decane + water + surfactant and decane + water + surfactant + CH₄/CO₂ systems. The IFTs of these surfactant-containing systems decrease with temperature and the influence of pressure is found to be less pronounced. The atomic density profiles show that the sulfate head groups of the SDS molecules penetrate the water-rich phase and their alkyl tails are stretched into the decane-rich phase. The sodium counterions of the surfactant molecules are located very close to their head groups. Furthermore, the density profiles of water and salt ions are hardly affected by the presence of the SDS molecules. However, the interfacial thickness between water and decane/CH₄/CO₂ molecules increases with increasing surfactant concentration. An important result is that the enrichment of CH₄ and/or CO₂ in the interfacial region decreases with increasing surfactant concentration. These results may be useful in the context of the water-alternating-gas approach that has been utilized during CO₂-enhanced oil recovery operations.

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.

Temperature and Pressure Dependence of Salt-Brine Dihedral Angles in the Subsurface

Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na⁺ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.

A hydrogen bond-modulated soft nanoscale water channel for ion transport through liquid-liquid interfaces

Ion transport through interfaces is of ubiquitous importance in many fields such as electrochemistry, emulsion stabilization, phase transfer catalysis, liquid-liquid extraction and enhanced oil recovery. However, the knowledge of interfacial structures that significantly affect ion transport through liquid-liquid interfaces is still lacking due to the difficulty of observing nanoscale interfaces. We studied here the evolution of interfacial structures during ion transport through the decane-water interface under different ionic concentrations and external forces using molecular dynamics simulations. The roles of hydrogen bonds in ion transport through interfaces are revealed. We identified a soft nanoscale channel during ion transport through liquid-liquid interfaces and the decane phase under specific external force. The stability of the water channel and the ion transport velocity both increase with ionic concentration due to the layered ordering structures of the water near the channel surface. We observed that the stability and connectivity of the water channel in the decane phase are remarkably improved both by the high increase of the number of hydrogen bonds in the water channel with increasing ionic concentration, and by the conformational change in water molecules near the water channel surface. Our discovery of a soft nanoscale water channel by molecular simulations implies that there is a potential stable passage for ion transport through liquid-liquid interfaces.

Molecular dynamics study of wetting of alkanes on water: from high temperature to the supercooled region and the influence of second inflection points of interfacial tensions

To explore the wetting behavior of alkanes on bulk water interfaces, molecular dynamics (MD) simulations were carried out for united-atom PYS alkane models, and for SPC/E and TIP4P/2005 water models over a wide temperature range. The MD results at each temperature were used to find (1) the surface tension of the alkanes (octane, nonane) and water, and (2) the interfacial tensions of the alkane-water systems. These quantities were then used to calculate the spreading coefficient (S) and contact angle (θc) for each alkane on water. At higher temperatures, the contact angle of octane and nonane on water is found to behave in accord with conventional expectations, i.e., it decreases with increasing temperature for both water models as each system approaches the usual high-temperature transition to perfect wetting. At lower temperatures, we found an unusual temperature dependence of S and θc for each PYS alkane on SPC/E water. In contrast to conventional expectations, θc decreases with a decrease in the temperature. For octane-SPC/E water, this unusual behavior of θc occurs due to the presence of second inflection points (SIP) in the vapor-water and the octane-water interfacial tensions, whereas the SIP effect is much less important for the nonane-water system. The unusual temperature dependence of θc observed for nonane on SPC/E water is also found for nonane on TIP4P/2005 water. On the other hand, such unusual wetting behavior has not been observed in the PYS octane-TIP4P/2005 water system, except possibly for the two lowest temperatures studied.

Atomistic insight into salinity dependent preferential binding of polar aromatics to calcite/brine interface: implications to low salinity waterflooding

This paper resolve the salinity-dependent interactions of polar components of crude oil at calcite-brine interface in atomic resolution. Molecular dynamics simulations carried out on the present study showed that ordered water monolayers develop immediate to a calcite substrate in contact with a saline solution. Carboxylic compounds, herein represented by benzoic acid (BA), penetrate into those hydration layers and directly linking to the calcite surface. Through a mechanism termed screening effect, development of hydrogen bonding between –COOH functional groups of BA and carbonate groups is inhibited by formation of a positively-charged Na⁺ layer over CaCO₃ surface. Contrary to the common perception, a sodium-depleted solution potentially intensifies surface adsorption of polar hydrocarbons onto carbonate substrates; thus, shifting wetting characteristic to hydrophobic condition. In the context of enhanced oil recovery, an ion-engineered waterflooding would be more effective than injecting a solely diluted saltwater.

Discovery of Novel Coumarin Analogs against the α-Glucosidase Protein Target of Diabetes Mellitus: Pharmacophore-Based QSAR, Docking, and Molecular Dynamics Simulation Studies

Diabetes mellitus (DM) is a chronic metabolic disease, the third killer of mankind. The finding of potent drugs against diabetes remains challenging. In the present study, coumarin derivatives with known biological activity against diabetic protein have been used to predict functional groups' positions on coumarin derivatives. α-Glucosidase is a brush border membrane-bound lysosomal enzyme from the hydrolase enzyme family. It plays an important role in the metabolism of glycoproteins. Inhibitors of lysosomal α-glucosidase can reduce postprandial hyperglycemia. Due to this, lysosomal α-glucosidase is a good therapeutic target for drugs. A total of 116 coumarin derivatives with IC50 values against lysosomal α-glucosidase were selected for a CADD (computer-aided drug design) approach to identify more potent drugs. Pharmacophore modeling and atom-based 3-QSAR of 116 active compounds against lysosomal α-glucosidase were performed and identified positions and types of groups to increase activity. We performed molecular docking of 116 coumarin derivatives against the lysosomal α-glucosidase enzyme, and three compounds (isorutarine, 10_, and 36) resulted in a docking score of -7.64, -7.12, and -6.86 kcal/mol. The molecular dynamics simulation of the above three molecules and protein complex performed for 100 ns supported the interaction stability of isorutarine, 10_, and 36 with the lysosomal binding site α-glucosidase.

Interfacial and bulk properties of concentrated solutions of ammonium nitrate

We conducted molecular dynamics (MD) simulations to calculate the density and surface tension of concentrated ammonium nitrate (AN) solutions up to the solubility limit of ammonium nitrate in water, by combining the SPC/E, SPCE/F and TIP4P/2005 water models with OPLS model for ammonium and nitrate ions. This is the first time that the properties of concentrated solutions of nitrates, especially AN, have been studied by molecular dynamics. We effectively account for the polarisation effects by the electronic continuum correction (ECC), practically realised via rescaling of the ionic charges. We found that, the full-charge force field MD simulations overestimate the experimental results, as the ions experience repulsion from the interface and prefer to remain in the subsurface layer and the bulk solution. In contrast, reducing the ionic charges results in the behaviour that fits well with the experimental data. The nitrate anions display a greater propensity for the interface than the ammonium cations. We accurately predict both the density and the rise in the surface tension of concentrated solutions of AN, recommending TIP4P/2005 for water and the scaled-charge OPLS model (OPLS/ECC) for the ions in the solutions. We observe that, the adsorption of anions to the interface accompanies their depletion in the subsurface layer, which is preferentially occupied by cations, resulting in the formation of the electric double layer. We demonstrate the ion deficiency for up to 3 Å below the surface and establish the requirement to include the polarisability effects in the OPLS model for AN. While these results confirmed the findings of the previous studies for dilute solutions, they are new in the solubility limit. Concentrated solutions exhibit a strong effect of the abundance of solute on the coordination numbers of ions and on the degree of ion pairing. Surprisingly, ion pairing decreases significantly at the interface compared with the bulk. The present study identifies OPLS/ECC, along with TIP4P/2005, to yield accurate predictions of physical properties of concentrated AN, with precision required for industrial applications, such as a formulation of emulsion and fuel-oil explosives that now predominate the civilian use of AN. An application of this model will allow one to predict the surface properties of supersaturated solutions of AN which fall outside the capability of the present laboratory experiments but are important industrially.

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)-block-poly(ε-caprolactone) (PEG-b-PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer-air and polymer-water interfaces, indicating the accuracy of the NB interactions. The structural properties (density, radius of gyration (R_g), and end-to-end distance (h)) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG₂₃-b-PCL₉ diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

Molecular dynamics simulation of CO2-switchable surfactant regulated reversible emulsification/demulsification processes of a dodecane-saline system

CO2-Switchable surfactants are of great potential in a wide range of industrial applications related to their ability to stabilize and destabilize emulsions upon command. Molecular dynamics simulations have been performed to reveal the fundamental mechanism of the reversible emulsification/demulsification processes of a dodecane-saline system by a CO2-switchable surfactant that switches between active (i.e., N'-dodecyl-N,N-dimethylacetamidinium (DMAAH+)) and inactive (i.e., N'-dodecyl-N,N-dimethylacetamidine (DMAA)) forms. The density profiles indicate that DMAAH+ could increase the oil-water interfacial thickness to a greater extent compared to DMAA. DMAAH+ could sharply reduce the interfacial tension of the dodecane-saline system, while DMAA only exhibits a limited decrease, which is in accordance with the experimental observation that DMAAH+/DMAA can reversibly emulsify/demulsify alkane-water systems. Our simulations showed that both the number and lifetime of hydrogen bonds (HBs) between DMAA and water are almost equal to those between DMAAH+ and water. In DMAA, the N atom connecting with the alkyl tail acted as a HB acceptor, while the N atom attached by a proton in DMAAH+ acted as a HB donor. Furthermore, the HBs between DMAAH+ and HCO3- at the interfaces are relatively limited. Hence, it is deduced that the HBs are insufficient to achieve the CO2-switchability of DMAA/DMAAH+. The Lennard Jones and coulombic potentials between DMAA/DMAAH+ and other species show that the coulombic potentials between DMAAH+ and water or anions (i.e., Cl- and HCO3-) sharply decrease with the increase of DMAAH+ and are much lower than those in models with DMAA. The enhanced coulombic interactions between DMAAH+ and anions lead to a remarkable reduction in interfacial tension and the emulsification of the alkane-saline system. Therefore, coulombic interactions are of crucial importance to the reversible emulsification/demulsification processes regulated by CO2-switchable surfactants, namely DMAAH+/DMAA.

Learning Coarse-Grained Potentials for Binary Fluids

For a multiple-fluid system, CG models capable of accurately predicting the interfacial properties as a function of curvature are still lacking. In this work, we propose a new probabilistic machine learning (ML) model for learning CG potentials for binary fluids. The water-hexane mixture is selected as a typical immiscible binary liquid-liquid system. We develop a new CG force field (FF) using the Shinoda-DeVane-Klein (SDK) FF framework and compute parameters in this CG FF using the proposed probabilistic ML method. It is shown that a standard response-surface approach does not provide a unique set of parameters, as it results in a loss function with multiple shallow minima. To address this challenge, we develop a probabilistic ML approach where we compute the probability density function (PDF) of parameters that minimize the loss function. The PDF has a well-defined peak corresponding to a unique set of parameters in the CG FF that reproduces the desired properties of a liquid-liquid interface. We compare the performance of the new CG FF with several existing FFs for the water-hexane mixture, including two atomistic and three CG FFs with respect to modeling the interface structure and thermodynamic properties. It is demonstrated that the new FF significantly improves the CG model prediction of both the interfacial tension and structure for the water-hexane mixture.

Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS2 nanopores

Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based on low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based on nanoporous materials have been suggested in recent years. One of the materials suggested is MoS₂, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct force fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were built based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, and the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through nanopores of two sizes were used - one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter - and two rigid water models from the TIP4P family - TIP4P/2005 and TIP4P/ε. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only can the water model influence the outcomes, but the ion model plays a crucial role when the pore is small enough.

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.

Surface Energies and Structure of Salt-Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt-brine and salt-salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. Here, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt-brine interface to provide insight into the stability of salt systems. Mixed NaCl-KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt-brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2-3 times (5.0 M systems). Size differences in Na⁺ and K⁺ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt-brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

Evidence for water ridges at oil-water interfaces: implications for ion transport

Understanding ion transport across interfaces is of fundamental importance in many processes such as liquid-liquid extraction, phase transfer catalysis, enhanced oil recovery and emulsion stabilisation. However, the factors that control ion transport across interfaces are poorly known due to a lack of knowledge of structural changes at interfaces. We studied here the effects of ionic concentration and external force on the transport of ions across the decane-water interface using classical molecular dynamics simulations. The results show that the evolution of interfacial structures during ion transfer across the interface is controlled by hydrogen bonding and ionic interactions at the interface. We also identified a new mode of ion transfer across the interface at low ionic concentrations, involving a 'water ridge', rather that the classical 'water finger'. In the water ridge mode, hydrogen bonds are not broken due to low ion levels, and the water ridge induces gradual interface deformation. Whereas, at high ionic concentrations, hydrogen bonds are broken by the strong ion electrostatic repulsion, thus inducing the formation of a water finger. We also found that the variation of the Gibbs free energy during ion transfer is directly relevant to the ionic concentration. The water ridge at low ionic concentrations, which displaces more water molecules towards the decane phase, induces less free energy variation than the water finger at high ionic concentrations.

On the transferability of ion parameters to the TIP4P/2005 water model using molecular dynamics simulations

Countless molecular dynamics studies have relied on available ion and water force field parameters to model aqueous electrolyte solutions. The TIP4P/2005 model has proven itself to be among the best rigid water force fields, whereas many of the most successful ion parameters were optimized in combination with SPC/E, TIP3P, or TIP4P/Ew water. Many researchers have combined these ions with TIP4P/2005, hoping to leverage the strengths of both parameter sets. To assess if this widely used approach is justified and to provide a guide in selecting ion parameters, we investigated the transferability of various commonly used monovalent and multivalent ion parameters to the TIP4P/2005 water model. The transferability is evaluated in terms of ion hydration free energy, hydration radius, coordination number, and self-diffusion coefficient at infinite dilution. For selected ion parameters, we also investigated density, ion pairing, chemical potential, and mean ionic activity coefficients at finite concentrations. We found that not all ions are equally transferable to TIP4P/2005 without compromising their performance. In particular, ions optimized for TIP3P water were found to be poorly transferable to TIP4P/2005, whereas ions optimized for TIP4P/Ew water provided nearly perfect transferability. The latter ions also showed good overall agreement with experimental values. The one exception is that no combination of ion parameters and water model considered here was found to accurately reproduce experimental self-diffusion coefficients. Additionally, we found that cations optimized for SPC/E and TIP3P water displayed consistent underpredictions in the hydration free energy, whereas anions consistently overpredicted the hydration free energy.

Challenges and future of chemical assisted heavy oil recovery processes

The primary method for heavy oil and bitumen production across the world is still in-situ steam-based technology. There are some drawbacks associated with steam-driven heavy oil recovery methods such as cyclic steam stimulation (CSS), steam flooding, and steam-assisted gravity drainage (SAGD). These cons include the high greenhouse gas footprint, low heavy oil/bitumen recovery, and difficulty in stop operation in emergency conditions. There exists a need for an improved method for recovering residual oils after applying steam injection. One of the potential technologies for doing this is chemical assisted heavy oil recovery, especially alkaline and surfactant additives. But the challenging question is how to develop a chemical-based oil recovery method considering long-term steam-rock interactions. Several associated issues of chemical additives, including adsorption behavior of surfactant at reservoir conditions and thermal stability of surfactant at steam chamber temperature, make this question more complex. This paper addresses all these concerns and provides solid knowledge regarding this technology. We delve into newly formulated chemicals for coupling with thermal oil recovery techniques that are still limited to lab-scale research, with the need for further studies. This critical review also provides the opportunities and challenges associated with chemical assisted heavy oil/bitumen production in a post-steam injection scenario. Finally, different aspects of such a method are covered in this review, along with practical information on field trials and best practices across the world.

Reducing uncertainty in simulation estimates of the surface tension through a two-scale finite-size analysis: thicker is better

Recent simulation studies of the surface tension γ, and other properties of thin free-standing films, have revealed unexpected finite size effects in which the variance of the properties vary monotonically with the in-plane width of the films, complicating the extrapolation of estimates of film properties to the thermodynamic limit. We carried out molecular dynamics simulations to determine the origin of this phenomenon, and to address the practical problem of developing a more reliable methodology for estimating γ in the thermodynamic limit. We find that there are two distinct finite size effects that must be addressed in a finite size analysis of γ in thin films. The first finite size scale is the in-plane width of the films and the second scale is the simulation cell size in the transverse direction. Increasing the first scale enhances fluctuations in γ, measured by the standard deviation of their distribution, while increasing the second reduces γ fluctuations due to a corresponding increased ‘freedom’ of the film to fluctuate out of plane. We find that using progressively large simulation cells in the transverse direction, while keeping the film width fixed to an extent in which the full bulk liquid zone is developed, allows us to obtain a smooth extrapolation to the thermodynamic limit, enabling a reduction of the γ uncertainty to a magnitude on the order of 1% for systems having a reasonably large size, i.e., O (1 μm).

The variance in the surface tension of systems under vapor/liquid equilibrium is strongly affected by the size of the interfacial area. Wider layers increase the variance, but these increments disappear as the interfacial area grows.

Effect of confinement on the adsorption behavior of inorganic and organic ions at aqueous-cyclohexane interfaces: a molecular dynamics study

Molecular dynamics simulation was used to study the adsorption behavior of inorganic and organic ions at aqueous-cyclohexane interfaces, and that with such systems confined in a kaolinite nanopore. Four aqueous solutions were used, each containing one of the four solutes (NaCl, NaOH, CaCl₂ and Ca(OH)₂) at the concentration of 1.0 M. At the interface of each of the solutions with neat cyclohexane, there existed an ion depletion zone. The more strongly hydrated ions, such as calcium and hydroxide, were more intensely depleted from the interface as compared to sodium and chloride. Such surface exclusion led to interfacial tension increases, with greater increments for solutions containing calcium or hydroxide. Upon addition of sodium decanoate to the cyclohexane phase, they partly migrated to the aqueous-organic interface, and the remaining formed inverted micelle complexes. Also, a small fraction of the solvated cations in the aqueous phase drifted to the depletion zone to interact with the organic anions, with their affinity towards the interface still being controlled by their hydration-strength. When such systems were confined in a kaolinite nanopore, behavior of decanoate anions - which were found to be prone to strong adsorption to solid surfaces in the absence of any aqueous solution - was determined by the nature of the solvated ions in water. In these systems, the more weakly-hydrated ionic species exhibited preferential adsorption to the solid surface, while affinity of the more strongly-hydrated ones was towards remaining within the water layers or to the aqueous-organic interface. With the calcium chloride solution, almost all of the organic ions were detached from the surface and adsorbed at the aqueous-cyclohexane interface. This was caused by the release of double amount of inorganic anions by this 1 : 2 salt, the inner-sphere adsorption of majority of chloride anions to the octahedral surface, and the higher charge density of calcium cations.