Molecular simulations of liquid-liquid interfacial properties: water-n-alkane and water-methanol-n-alkane systems

Molecular Dynamics Study on the Heat of Phase Transition of Chloroacetate from the Bulk to the Surface Phases

Based on the surface adsorption theory developed in our group and molecular dynamics simulations, the adsorption behaviors of methyl chloroacetate, ethyl chloroacetate, propyl chloroacetate, methyl dichloroacetate, and methyl trichloroacetate at the vapor-liquid interface have been thoroughly investigated. The surface layer thickness of these five liquids shows a significant increase as temperature increases, but it decreases significantly with the increase of intermolecular interactions of chloroacetates. It is found with our surface adsorption theory that the reversible heat of phase transition for chloroacetates from bulk phase to surface phase increases from methyl chloroacetate to propyl chloroacetate with the alkyl chain and also increases from methyl chloroacetate to methyl trichloroacetate as the number of chloro group increases. Molecular dynamics simulation is exploited to calculate the entropies of these five liquids in the surface and bulk regions at various temperatures. The variation trends of the simulated results are consistent with those of the phase transition heats determined with our surface adsorption theory.

The Effect of Spacer Chain Length on Foam Stability of Mixed Amine Gemini Surfactant/NaOl by Molecular Dynamic Simulations

Foam stability critically determines the efficiency of the mineral flotation process. Although the mixed amine Gemini surfactant/anionic surfactants exhibit excellent flotation performance, atomic-level investigations of the mechanism of their impact on foam stability remain limited. This study employs molecular dynamics simulations to investigate the self-aggregation behavior of mixed amine Gemini surfactant/sodium oleate (NaOl) systems with varying spacer chain lengths at the air/water interface. The structural parameters of self-aggregation, surface tension, synergistic energy, and diffusion coefficient of water molecules were calculated in detail. The results of molecular dynamics simulations indicated that synergistic adsorption between surfactants occurred. Compared with single amine Gemini surfactant systems, the mixed surfactant systems exhibited an enhanced interfacial activity. The spacer chain length significantly affected the adsorption configurations of the mixed surfactant at the air/water interface. For spacer chains containing fewer than five methylene groups, carboxyl groups preferentially adsorbed between two intramolecular amine groups, forming independently clustered aggregates. Conversely, longer spacer chains promoted adsorption between carboxyl groups and intermolecular amine groups, forming interconnected network-like aggregates. Both structural configurations constrained interfacial water mobility, thereby reducing the liquid flow rate between foam films, suppressing water loss and enhancing the mechanical stability of flotation foams.

Molecular Insights into the Concentration Dependent Promotion Effect of Tetrabutylammonium Bromide on Hydrate Growth: A Molecular Dynamics Simulation Study

Tetrabutylammonium bromide (TBAB) has been proven to improve the growth of hydrate via experimental methods, which may be attributed to its different concentrations. In this study, the molecular dynamics (MD) simulation is employed to investigate the concentration dependent promotion effect of TBAB on the growth of CO2 hydrate. The tetrahedral order parameter, number of cages, hydrate crystal growth trajectories, significant microconfigurations, and distribution of CO2 and TBAB are analyzed in detail. It is found that the promotion effect of TBAB is more prominent at low concentrations (5 and 10 wt %) than that at high concentrations (15 and 20 wt %). During the growth of hydrate crystal, tetrabutylammonium (TBA+) ions adsorb on the hydrate crystal and serve as guest molecules to form TBA+ semiclathrate hydrate cages. Then, the TBA+ semiclathrate hydrate cages undergo a self-adjustment process and induce the generation of CO2 hydrate cages. At high concentrations, the great accumulation of TBA+ at the hydrate-liquid interface disturbs the effective adsorption and self-adjustment processes, and the tightly packed arrangement of TBA+ at the gas-liquid interface partially inhibits the mass transfer of CO2. This study provides visible mechanisms of the concentration dependent promotion effect of TBAB from the microscopic level, which complements the vacancy in experimental studies.

Improved Emulsifying Performance of Agarose Microgels by Cross-Interfacial Diffusion of Polyphenols

Noninterface-active polysaccharides can acquire better emulsifying properties through microgelation, yet optimizing their emulsifying performance remains a significant challenge. This study introduces a novel approach to enhance the emulsifying performance of polysaccharide microgels by leveraging the cross-interfacial diffusion of polyphenols, which promotes the interfacial adsorption of microgels. Tannic acid (TA) was predispersed in oil phases and subsequently emulsified with agarose microgel (AM) suspensions, and the impacts of TA diffusion on the emulsifying performance of AMs was investigated. In addition, the transmittance profiles of oil-water biphasic systems were found to innovatively indicate the cross-interfacial diffusion of TA and the interfacial adsorption of AMs. The current results suggest that an appropriate level of TA incorporation can benefit the emulsifying performance of AMs, correlating with decreased droplet sizes and improved physical stability of the emulsion. However, excessive TA might trigger the clustering of AMs before they reach the interfacial layer, adversely affecting the emulsion stability. In conclusion, the cross-interfacial diffusion of polyphenols offers a promising strategy to overcome the stability challenges encountered in polysaccharide microgel-stabilized emulsions.

Effect of molecular branching and surface wettability on solid-liquid surface tension and line-tension of liquid alkane surface nanodroplets

HYPOTHESIS

Surface nanodroplets have important technological applications. Previous experiments and simulations have shown that their contact angle deviates from Young's equation. A modified version of Young's equation considering the three-phase line tension (τ) has been widely used in literature, and a wide range of values for τ are reported. We have recently shown that molecular branching affects the liquid-vapour surface tension γlv of liquid alkanes. Therefore, the wetting behaviour of surface nanodroplets should be affected by molecular branching. This study conducted molecular dynamics (MD) simulations to gain insight into the wetting behaviour of linear and branched alkane nanodroplets on oleophilic and oleophobic surfaces. We aim to examine the Young equation's validity and branching's effect on fundamental properties, including solid-liquid surface tension γsl and line tension τ.

SIMULATIONS

The simulations were performed on a linear alkane, triacontane (C30H62), as well as four of its branched isomers: 2,6,13,17-tetrapropyloctadecane,2,6,9,10,13,17-hexaethyloctadecane, 2,5,7,8,11,12,15-heptaethylhexadecane and 2,3,6,7,10,11-hexapropyldodecane. Nanodroplets with a diameter of approximately 15 nm were released onto the surfaces, and their contact angles were measured. Additionally, using a novel approach, the solid-liquid surface tension (γsl), the validity of Young's equation and line tension for all alkane and surface combinations are determined.

FINDINGS

It was discovered that the calculated γsl, deviated from the theoretical γsl,Young predicted from Young's equation for all alkanes on oleophilic surfaces. However, this deviation was minimal for branched alkanes on the oleophobic surfaces but more significant for the linear alkane. The findings indicated that γsl < 0 for oleophilic surfaces and γsl > 0 for oleophobic surfaces. Moreover, it was observed that |γsl| was lower for branched molecules and decreased as branching increased. Line tension values were then determined through a novel method, showing τ was positive for oleophilic surfaces ranging from 1.30 × 10-10 to 6.27 × 10-11N. On an oleophobic surface, linear alkane shows a negative line tension of -1.15 × 10-10N and branched alkanes up to two orders of magnitude lower values ranging from -2.09 × 10-12 to 2.43 × 10-11N. Line tension values between -1.15 × 10-10 and + 1.1 × 10-10N are calculated for various linear alkane and surface combinations. These findings show the dependence of line tension on the contact angle and branching, demonstrating that for linear alkanes, τ is significant, whereas, for branched alkanes, line tension is smaller or negligible for large contact angles.

Surfactant Synergistic Effect and Interfacial Properties of Microemulsions Compounded with Anionic and Nonionic Surfactants Using Dissipative Particle Dynamics

Microemulsions are one of the most promising directions in enhanced oil recovery, but conventional screening methods are time-consuming and labor-intensive and lack the means to analyze them at the microscopic level. In this paper, we used the Clint model to predict the changes in the synergistic effect of the mixed system of anionic surfactant sodium dodecyl benzenesulfonate and nonionic surfactant polyethoxylated fatty alcohols (C12E6), generated microemulsions using surfactant systems with different mole fractions, and used particle size to analyze the performance and stability of microemulsions, analyze the properties and stability of microemulsions using particle size, and analyze the interfacial behaviors and changes of microemulsions when different systems constitute microemulsions from the point of view of mesoscopic microemulsion self-assembly behaviors by combining with dissipative particle dynamics. It has been shown that microemulsion systems generated from anionic and nonanionic surfactants with a synergistic effect, based on the Clint model, exhibit excellent performance and stability at the microscopic level. The method proposed in this paper can dramatically improve the screening efficiency of microemulsions of anionic and nonanionic surfactants and accurately analyze the properties of microemulsions, so as to provide a theoretical basis for the subsequent research on microemulsions.

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.

Experiments, molecular dynamics simulations, and quantum chemistry calculations on the effect of gemini surfactants' headgroup on the oil-water interfacial tension

The effect of the Gemini surfactant headgroup on the oil-water interfacial tension has yet to be systematically revealed. In this work, anionic Gemini surfactants with different hydrophilic headgroups (carboxylic, sulfuric, and sulfonic) were designed and synthesized. The oil-water interfacial tension was tested. The essential parameters for evaluating the interface characteristics, including the oil-water interfacial layer thickness, the coordination number, and the diffusion coefficient, were calculated employing molecular dynamics simulation. The surface electrostatic potential explained the quantitative mechanism of the hydrophobicity and lipophilicity of three types of Gemini surfactants through quantum chemical calculations. The oil-water interfacial tension difference of the Gemini surfactants was revealed at the electronic level. This paper will provide theoretical guidance for designing Gemini surfactants with a high-efficiency performance to enhance oil recovery.

Molecular Dynamics Simulations on the Adsorbed Monolayers of N-Dodecyl Betaine at the Air-Water Interface

Betaine is a kind of zwitterionic surfactant with both positive and negative charge groups on the polar head, showing good surface activity and aggregation behaviors. The interfacial adsorption, structures and properties of n-dodecyl betaine (NDB) at different surface coverages at the air-water interface are studied through molecular dynamics (MD) simulations. Interactions between the polar heads and water molecules, the distribution of water molecules around polar heads, the tilt angle of the NDB molecule, polar head and tail chain with respect to the surface normal, the conformations and lengths of the tail chain, and the interfacial thickness of the NDB monolayer are analyzed. The change of surface coverage hardly affects the locations and spatial distributions of the water molecules around the polar heads. As more NDB molecules are adsorbed at the air-water interface, the number of hydrogen bonds between polar heads and water molecules slightly decreases, while the lifetimes of hydrogen bonds become larger. With the increase in surface coverage, less gauche defects along the alkyl chain and longer NDB chain are obtained. The thickness of the NDB monolayer also increases. At large surface coverages, tilted angles of the polar head, tail chain and whole NDB molecule show little change with the increase in surface area. Surface coverages can change the tendency of polar heads and the tail chain for the surface normal.

An in silico osmotic pressure approach allows characterization of pressure-area isotherms of lipid monolayers at low molecular areas

Surface pressure-area isotherms of lipid monolayers at the air-water interface provide essential information about the structure and mechanical behaviour of lipid membranes. These curves can be readily obtained through Langmuir trough measurements and, as such, have been collected for decades in the field of membrane biochemistry. However, it is still challenging to directly observe and understand nanoscopic features of monolayers through such experiments, and molecular dynamics (MD) simulations are generally used to provide a molecular view of such interfaces. In MD simulations, the surface pressure-area (Π-A) isotherms are generally computed using the Kirkwood-Irving formula, that relies on the evaluation of the pressure tensor. This approach, however, has intrinsic limitations when the molecular area in the monolayer is low (typically < 60 Ų per lipid). Recently, an alternative method to compute Π-A isotherms of surfactants, based on the calculation of the three-dimensional osmotic pressure via the implementation of semipermeable barriers was proposed. In this work, we investigate the feasibility of this approach for long-chain surfactants such as phospholipids. We identify some discrepancies between the computed values and experimental results, and we propose a semi-empirical correction based on the molecular structure of the surfactants at the monolayer interface. To validate the potential of this new approach, we simulate several phosphatidylcholine and phosphatidylethanolamine lipids at various temperatures using all-atom and coarse-grained force fields, and we compute the corresponding Π-A isotherms. Our results show that the Π-A isotherms obtained using the new method are in very good agreement with experiments and far superior to the canonical pressure tensor-based method at low molecular areas. This corrected osmotic pressure method allows for accurate characterization of the molecular packing in monolayers in various physical phases.

Detecting Liquid-Liquid Phase Separations Using Molecular Dynamics Simulations and Spectral Clustering

A stringent test of the accuracy of empirical force fields is reproducing the phase diagram of bulk phases and mixtures. Exploring the phase diagram of mixtures requires the detection of phase boundaries and critical points. In contrast to most solid-liquid transitions, in which a global order parameter (average density) can be used to discriminate between two phases, some demixing transitions entail relatively subtle changes in the local environment of each molecule. In such cases, finite sampling errors and finite-size effects make the identification of trends in local order parameters extremely challenging. Here we analyze one such example, namely a methanol/hexane mixture, and compute several local and global structural properties. We simulate the system at various temperatures and study the structural changes associated with demixing. We show that despite a seemingly continuous transformation between mixed and demixed states, the topological properties of the H-bond network change abruptly as the system crosses the demixing line. In particular, by using spectral clustering, we show that the distribution of cluster sizes develops a fat tail (as expected from percolation theory) in the vicinity of the critical point. We illustrate a simple criterion to identify this behavior, which results from the emergence of large system-spanning clusters from a collection of aggregates. We further tested the spectral clustering analysis on a Lennard-Jones system as a standard example of a system with no H-bonds, and also, in this case, we were able to detect the demixing transition.

Study on the Stability of Bio-Oil Modified Prime Coat Oil Based on Molecular Dynamics

To explore the effect of different emulsifier contents on the stability performance of biomass-emulsified asphalt, three types of emulsified asphalt with 1%, 3%, and 5% anionic emulsifiers were prepared and analyzed by molecular dynamics simulation and macroscopic experiments. Firstly, we used molecular simulation software (Material Studio, MS) to construct a model of biomass-emulsified asphalt with different emulsifier contents and analyzed the microscopic mechanism of the emulsifier to improve the stability of the emulsified asphalt by the radial distribution function, interaction energy, interfacial layer thickness, and solubility parameters of the emulsified asphalt system with different emulsifier contents. The results were validated by macro and micro tests including storage stability, particle size determination, and infrared spectroscopy. The results show that at low emulsifier contents, the emulsifier can reduce the interfacial tension between the oil-water interface and expand the transition region between the two phases (interfacial layer thickness), which will prevent interparticle agglomeration and reduce the emulsion particle size, thus reducing the settling rate and ensuring the stability of the emulsion. When the emulsifier content is further increased beyond the critical micelle concentration, the emulsifiers will agglomerate with each other and show larger peaks in the radial distribution function, and the phenomenon of emulsifier agglomeration will appear in the five-day storage stability test, resulting in a corresponding decrease in the proximity of the infrared absorption peak area ratio in the same wavelength band of the upper and lower layers of the biomass-emulsified asphalt, and the emulsion stability decreases instead.

Comprehensive review of the interfacial behavior of water/oil/surfactant systems using dissipative particle dynamics simulation

A comprehensive understanding of interfacial behavior in water/oil/surfactant systems is critical to evaluating the performance of emulsions in various industries, specifically in the oil and gas industry. To gain fundamental knowledge regarding this interfacial behavior, atomistic methods, e.g., molecular dynamics (MD) simulation, can be employed; however, MD simulation cannot handle phenomena that require more than a million atoms. The coarse-grained mesoscale methods were introduced to resolve this issue. One of the most effective mesoscale coarse-grained approaches for simulating colloidal systems is dissipative particle dynamics (DPD), which bridges the gap between macroscopic time and length scales and molecular-scale simulation. This work reviews the fundamentals of DPD simulation and its progress on colloids and interface systems, especially surfactant/water/oil mixtures. The effects of temperature, salt content, a water/oil ratio, a shear rate, and a type of surfactant on the interfacial behavior in water/oil/surfactant systems using DPD simulation are evaluated. In addition, the obtained results are also investigated through the lens of the chemistry of surfactants and emulsions. The outcome of this comprehensive review demonstrates the importance of DPD simulation in various processes with a focus on the colloidal and interfacial behavior of surfactants at water-oil interfaces.

Biphasic Behaviors of Nd3+ Bound with Cyanex272, Cyanex301, and Cyanex302: A Molecular Dynamics Simulation Study

By means of molecular dynamics simulations, this work addresses the conformational flexibility and migration of trivalent neodymium (Nd³⁺) coordinated with three or six titled (thio)phosphinic ligands and shows that the fluxionality of the complexes enables them to adapt to the solvent environment during the migration. Cyanex272 forms a more compact complex than the other two types of ligands and screens more significantly the interaction between the water solvent and the metal ion in the complex, which weakens the detainment of the aqueous environment. This results in faster motion of the Nd(C272)₃ complex both in its translation and rotation than the other complexes when migrating to the organic phase and wins over the other two ligands in transporting the metal ions from the aqueous phase to the organic phase. Depending on the solvent environment, these complexes may take two types of conformations to balance the forces from the environment benefited from their fluxionality. The migration of the M:L = 1:6 complexes, Nd[H(C272)₂]₃ and Nd[H(C301)₂]₃, was also investigated. The rich presence of the alkyl groups in the complexes screens the influence of the aqueous environment and benefits the transportation of metal ions to the interface. This work is expected to contribute to the community of inorganic chemistry interested in the coordination chemistry of metal ions and their behaviors in the condensed phase.

Compact Peptoid Molecular Brushes for Nanoparticle Stabilization

Controlling the interfaces and interactions of colloidal nanoparticles (NPs) via tethered molecular moieties is crucial for NP applications in engineered nanomaterials, optics, catalysis, and nanomedicine. Despite a broad range of molecular types explored, there is a need for a flexible approach to rationally vary the chemistry and structure of these interfacial molecules for controlling NP stability in diverse environments, while maintaining a small size of the NP molecular shell. Here, we demonstrate that low-molecular-weight, bifunctional comb-shaped, and sequence-defined peptoids can effectively stabilize gold NPs (AuNPs). The generality of this robust functionalization strategy was also demonstrated by coating of silver, platinum, and iron oxide NPs with designed peptoids. Each peptoid (PE) is designed with varied arrangements of a multivalent AuNP-binding domain and a solvation domain consisting of oligo-ethylene glycol (EG) branches. Among designs, a peptoid (PE5) with a diblock structure is demonstrated to provide a superior nanocolloidal stability in diverse aqueous solutions while forming a compact shell (∼1.5 nm) on the AuNP surface. We demonstrate by experiments and molecular dynamics simulations that PE5-coated AuNPs (PE5/AuNPs) are stable in select organic solvents owing to the strong PE5 (amine)-Au binding and solubility of the oligo-EG motifs. At the vapor-aqueous interface, we show that PE5/AuNPs remain stable and can self-assemble into ordered 2D lattices. The NP films exhibit strong near-field plasmonic coupling when transferred to solid substrates.

Effect of Oil-Displacing Agent Composition on Oil/Water Interface Stability of the Asphaltene-Rich ASP Flooding-Produced Water

In this work, the effect of oil-displacing agent composition on oil/water interface stability of the asphaltene-rich alkali-surfactant-polymer (ASP) flooding-produced water was systematically investigated, especially from the perspective of the interaction between oil displacement agents and asphaltene at the oil/water interface. Primarily, adsorption behavior of the artificial and natural interfacial substances (oil displacement surfactant and asphaltenes) on the oil/water interface was investigated by molecular dynamics simulation. The oil displacement surfactant and asphaltenes formed a cross-linked and compact interfacial film structure, which significantly enhanced the interface stability; the more the oil displacement surfactants adsorbed on the interface, the more stable is the cross-linked structure formed between them and asphaltenes. Then, the interfacial property variations that are originating from the interactions differences between oil displacement agents and asphaltenes were monitored via interfacial tension, zeta potential, and interfacial film rheology tests. Moreover, the effect of oil displacement agent concentrations on the interfacial film thinning and rupture kinetic behavior was further investigated. Finally, cream experiments were conducted to verify the effect of oil displacement agent composition on the oil/water separation efficiencies of asphaltene-rich ASP flooding-produced water. When 5% asphaltenes was added, the creaming oil removal rate reduced from 90.0 to 85.3% at 19 h. The interactions between asphaltenes and oil displacement agents immensely enhance the oil/water interfacial film strength and impede the oil/water separation process.

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.

Molecular Insight into the Adsorption Thermodynamics and Interfacial Behavior of GOs at the Liquid-Liquid Interface

The adsorption of two-dimensional (2-D) graphene oxide (GO) nanosheets at liquid-liquid interfaces has broad technological implications from functional material preparations to oil-water emulsification. Molecular-level understanding of the adsorption thermodynamics and the interfacial behavior is of great significance. Here, the adsorption free energy of GO nanosheets at the water-cyclohexane system was simulated, in which the effect of oxygen-containing groups and deprotonation has been investigated. It was observed that the neutral GO (GO-COOH) has obvious interfacial activity with a reduction of interfacial tension, while the deprotonated GO (GO-COO^-) shows a weak interface affinity. There exists an optimal oxidization degree that could cause the best interfacial stability, which is attributed to the balance of interfacial hydrophilic-hydrophobic interactions. The interaction arising from water is the main factor determining interfacial activity. The interfacial morphology and dynamics of GO nanosheets have also been simulated, in which an anisotropic 2-D translation and rotation along the interface were revealed. Our simulation results provide new insight into the adsorption mechanism and dynamics behavior of GO at the oil-water interface.

Development of coarse-grained force field for alcohols: an efficient meta-multilinear interpolation parameterization algorithm

Coarse-grained (CG) molecular dynamics are powerful tools to access a mesoscopic phenomenon and simultaneously record microscopic details, but currently the CG force fields (FFs) are still limited by low parameterization efficiency and poor accuracy especially for polar molecules. In this work, we developed a Meta-Multilinear Interpolation Parameterization (Meta-MIP) algorithm to optimize the CG FFs for alcohols. This algorithm significantly boosts parameterization efficiency by constructing on-the-fly local databases to cover the global optimal parameterization path. In specific, an alcohol molecule is mapped to a heterologous model composed of an OH bead and a hydrocarbon portion which consists of alkane beads representing two to four carbon atoms. Non-bonded potentials are described by soft Morse functions that have no tail-corrections but can still retain good continuities at truncation distance. Nearly all of the properties in terms of density, heat of vaporization, surface tension, and solvation free energy for alcohols predicted by the current FFs deviate from experimental values by less than 7%. This Meta-MIP algorithm can be readily applied to force field development for a wide variety of molecules or functional groups, in many situations including but not limited to CG FFs.

Negative Pressure within a Liquid-Fluid Interface Determines Its Thickness

The density within the interface between two fluid phases at equilibrium gradually changes from that of one phase to that of the other. The main change in density, according to experimental measurements, practically occurs over a finite distance of O [1 nm]. If we assume that the average stress difference within the interface is on the order of magnitude of ambient pressure, then the Bakker equation implies that for a liquid with surface tensions, say ∼50 mN/m, we get an interface thickness of ∼500 nm. This is certainly too big because it contradicts experimental findings. Alternatively, if the thickness is assumed to be O [10 nm] or less, as is usually believed, the average stress difference must be ∼5 × 10⁶ N/m² or bigger, which is surprisingly high. This paper shows using a few approaches that such a high average stress difference is due to negative stresses in the interface.

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.

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

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

Interfacial properties of the ionic liquid [bmim][triflate] over a wide range of temperatures

We carried out molecular dynamics simulations of the liquid/vacuum equilibrium of the ionic liquid [bmim][triflate] in a wide range of temperatures (323.15 to 573.15 K). The results showed liquid phases with high densities even at temperatures close to the decomposition temperature of the liquid. The density and surface tension behaviors are linear across this wide range of temperatures, which is an extension of the behaviors of these systems at low temperatures, where these properties have been experimentally measured. The interfacial region shows peaks of adsorption of the ions; they are ordered, with the alkyl chains of the [bmim] cations pointing out of the liquid, and the tailing angle of the chains becomes 90° at higher temperatures. The alkyl chains are part of the outermost interfacial region, where intra- and intermolecular tangential forces are in equilibrium; thus, they do not contribute to the total surface tension. Unlike simpler organic liquids, the surface tension is composed of positive normal contributions of intermolecular interactions; these are almost in equilibrium with the negative normal contributions of intramolecular interactions, which are mainly vibrations of the distance and the angle of valence. The pressure profiles show that the molecules are in 'crushed' conformations internally in the bulk liquid and even more so in the normal direction at the interface. The total pressure profiles show values with physical meaning, where the tangential peaks show higher values than normal pressures and give rise to the surface tension. Short cutoff radii for the calculation of intermolecular forces (less than 16.5 Å) produce a system that is not mechanically stable in the region of the bulk liquid (confirmed by radial distribution function calculations); this produces a difference between the normal pressure and the average of the tangential pressures, which affects the calculation of the surface tension due to overestimation by up to 20% when using the global expression, which is extensively used for the calculation of surface tension. The use of a sufficiently long cutoff radius avoids these mechanical balance problems.

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

In this study, classical molecular dynamic simulations have been used to examine the molecular properties of the water-alkane interface at various NaCl salt concentrations (up to 3.0 mol/kg). A variety of different force field combinations have been compared against experimental surface/interfacial tension values for the water-vapour, decane-vapour and water-decane interfaces. Six different force fields for water (SPC, SPC/E, TIP3P, TIP3Pcharmm, TIP4P & TIP4P2005), and three further force fields for alkane (TraPPE-UA, CGenFF & OPLS) have been compared to experimental data. CGenFF, OPLS-AA and TraPPE-UA all accurately reproduce the interfacial properties of decane. The TIP4P2005 (four-point) water model is shown to be the most accurate water model for predicting the interfacial properties of water. The SPC/E water model is the best three-point parameterisation of water for this purpose. The CGenFF and TraPPE parameterisations of oil accurately reproduce the interfacial tension with water using either the TIP4P2005 or SPC/E water model. The salinity dependence on surface/interfacial tension is accurately captured using the Smith & Dang parameterisation of NaCl. We observe that the Smith & Dang model slightly overestimates the surface/interfacial tensions at higher salinities (>1.5 mol/kg). This is ascribed to an overestimation of the ion exclusion at the interface.

The Ninth Industrial Fluid Properties Simulation Challenge

The Ninth Industrial Fluid Properties Simulation Challenge aimed to test the ability of molecular modeling approaches to predict water/oil interfacial tension (IFT) at conditions of high temperature and pressure. In particular, the challenge featured water/oil IFT where the oil was n-dodecane, toluene, or a 50:50 n-dodecane/toluene blend at 1.825 MPa and temperatures in the range of 383 K to 443 K. Seven entries were received including approaches such as molecular dynamics (MD) and Monte Carlo (MC) simulations, COSMO-RS, and iSAFT, and they were judged by comparison to pendant drop tensiometer benchmark data. The quality of predictions varied among the entries between approximately 20 % and 70 % of the total points possible with the entries based on MD and MC having the highest scores in most cases. As is often the case in molecular modeling, predictions of the relative trends tended to be reliable even if the absolute values were not.

Quantitative prediction of the position and orientation for an octahedral nanoparticle at liquid/liquid interfaces

Shape-controlled polyhedral particles and their assembled structures have important applications in plasmonics and biosensing, but the interfacial configurations that will critically determine their resultant assembled structures are not well-understood. Hence, a reliable theory is desirable to predict the position and orientation of a polyhedron at the vicinity of a liquid/liquid interface. Here we demonstrate that the free energy change theory can quantitatively predict the position and orientation of an isolated octahedral nanoparticle at a liquid/liquid interface, whose vertices and facets can play crucial roles in biosensing. We focus on two limiting orientations of an octahedral nanoparticle, vertex up and facet up. Our proposed theory indicates that the surface wettability (hydrophilic/hydrophobic ratio) of the nanoparticle determines its most stable position and the preferred orientation at a water/oil interface. The surface wettability of an octahedron is adjusted from extremely hydrophobic to extremely hydrophilic by changing the amount of charge on the Ag surface in molecular dynamics (MD) simulations. The MD simulations results are in excellent agreement with our theoretical prediction for an Ag octahedral nanoparticle at a hexane/water interface. Our proposed theory bridges the gap between molecular-level simulations and equilibrium configurations of polyhedral nanoparticles in experiments, where insights from nanoparticle intrinsic wettability details can be used to predict macroscopic superlattice formation experimentally. This work advances our ability to precisely predict the final structures of the polyhedral nanoparticle assemblies at a liquid/liquid interface.

A coarse-grain molecular dynamics study of oil-water interfaces in the presence of silica nanoparticles and nonionic surfactants

In this work, we have studied the effect of hydrophilic silica nanoparticles (NPs), in the presence of nonionic surfactants (Triethylene glycol monododecyl ether and Tween 20), on the oil-water (n-octane-water, n-dodecane-water and n-hexadecane-water) interfacial tensions (IFTs) at 300 K, using coarse-grained molecular dynamics simulations based on the MARTINI force field. Simulation results indicate that silica NPs solely do not affect the IFT. However, the silica NPs may or may not increase the IFT of oil-water containing nonionic surfactant, depending on the tendency of the surfactant to adsorb on the surface of NPs. The adsorption occurs due to the formation of hydrogen bonds, and adsorption increases with a decrease in pH, as seen in experimental studies. In this work, we found that the oil-water IFT increases with an increasing amount of adsorption of the surfactant on NPs. At a fixed amount of adsorption of the surfactant on NPs, the IFT behavior is indifferent to the change in concentration of NPs. However, the IFT decreases with an increase in surfactant concentration. We present a detailed analysis of the density profile and intrinsic width of the interface. The IFT behavior is found to correlate extremely well with the intrinsic width of the interface. The current study provides an explanation for the increase in IFT observed in a recent experiment [N. R. Biswal et al., J. Phys. Chem. B 120, 7265-7274 (2016)] for various types of NPs and nonionic surfactant systems.

Desalination by dragging water using a low-energy nano-mechanical device of porous graphene

We propose a dragging nano-structured suction system based on graphene sheets for water desalination processes.

Computer modelling of the surface tension of the gas–liquid and liquid–liquid interface

This review presents the state of the art in molecular simulations of interfacial systems and of the calculation of the surface tension from the underlying intermolecular potential.

Three Dimensional Nano "Langmuir Trough" for Lipid Studies

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