Quantitative Predictions of the Interfacial Tensions of Liquid-Liquid Interfaces through Atomistic and Coarse Grained Models

Molecular Insights into the Wettability and Adsorption of Acid Gas-Water Mixture

Sequestration of acid gas in geological formations is a disposal method with potential economic and environmental benefits. The process is governed by variables such as gas-water interfacial tension, wetting transition, and gas adsorption into water, among other things. However, the influence of the pressure and temperature on these parameters is poorly understood. This study investigates these parameters using coarse-grained molecular dynamics (CG-MD) simulations and density gradient theory (DGT). Simulations were carried out at 313.15 K and a pressure range of 0-15 MPa. A comparison was made against H2S-water systems to clarify the effects of adsorption on interfacial tension due to vapor-liquid-liquid equilibrium. The predicted H2S-water interfacial tension and phase densities by CG-MD and DGT matched the experimental values well. The adsorption can be quantified via the Gibbs Adsorption function Γ12, which correlated well with the three-phase transition. On the one hand, pressure increments below the three-phase transition revealed a significant adsorption of H2S. On the other hand, above the three-phase transition, the Gibbs Adsorption capacity remained constant, which indicated a saturation of H2S at the water surface due to liquid-liquid equilibrium. Finally, H2S behaves markedly differently in wetting transition, rather than the involved for CO2 to different molecular layers beneath the surface of aqueous solutions. In this respect, H2S is represented by a first-order wetting transition while CO2 presents a critical wetting. Finally, it has also been found that the preferential adsorption of H2S over the H2O interface is greater if compared to that of CO2, due to its strong interaction with water. In fact, we have also demonstrated that CO2 under triphasic conditions strongly influences the wetting of the ternary system.

Liquid-Liquid Interface-Assisted Self-Assembly of Ag-Doped ZnO Nanosheets for Atomic Switch Application

Developing facile and inexpensive methods for obtaining large-area two-dimensional semiconducting nanosheets is highly desirable for mass-scale device application. Here, we report a method for producing uniform and large-area films of a Ag-doped ZnO (AZO) nanosheet network via self-assembly at the hexane-water interface by controlling the solute/solvent ratio. The self-assembled film comprises of uniformly tiled nanosheets with size ∼1 μm and thicknesses∼60-100 nm. Using these films in a Pt/AZO/Ag structure, an atomic switch operation is realized. The switching mechanism is found to be governed by electrochemical metallization with nucleation as the rate-limiting step. Our results establish the protocol for large-scale device applications of AZO nanosheets for exploring advanced atomic switch-based neuromorphic systems.

Thermodynamics and morphology of linear multiblock copolymers at homopolymer interfaces

Block copolymers at homopolymer interfaces are poised to play a critical role in the compatibilization of mixed plastic waste, an area of growing importance as the rate of plastic accumulation rapidly increases. Using molecular dynamics simulations of Kremer-Grest polymer chains, we have investigated how the number of blocks and block degree of polymerization in a linear multiblock copolymer impacts the interface thermodynamics of strongly segregated homopolymer blends, which is key to effective compatibilization. The second virial coefficient reveals that interface thermodynamics are more sensitive to block degree of polymerization than to the number of blocks. Moreover, we identify a strong correlation between surface pressure (reduction of interfacial tension) and the spatial uniformity of block junctions on the interface, yielding a morphological framework for interpreting the role of compatibilizer architecture (number of blocks) and block degree of polymerization. These results imply that, especially at high interfacial loading, the choice of architecture of a linear multiblock copolymer compatibilizing surfactant does not greatly affect the modification of interfacial tension.

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 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.

Density functional theory for the microscopic structure of nanoparticles at the liquid-liquid interface

We propose an extension of the density functional approach to study the structure and thermodynamic properties of a system comprising a certain amount of nanoparticles at the interface between two partially miscible liquids. Model calculations have been carried out for a binary symmetric mixture of Yukawa fluids and hard-sphere nanoparticles. Despite its simplicity, the model captures the principal features of this type of system. The results indicate that nanoparticles form layers and the number of layers depends on the amount of nanoparticles and on their diameters. For the systems studied the formation of the layers evidences strong localization of the nanoparticles at the interface.

Modelling realistic microgels in an explicit solvent

Thermoresponsive microgels are polymeric colloidal networks that can change their size in response to a temperature variation. This peculiar feature is driven by the nature of the solvent-polymer interactions, which triggers the so-called volume phase transition from a swollen to a collapsed state above a characteristic temperature. Recently, an advanced modelling protocol to assemble realistic, disordered microgels has been shown to reproduce experimental swelling behavior and form factors. In the original framework, the solvent was taken into account in an implicit way, condensing solvent-polymer interactions in an effective attraction between monomers. To go one step further, in this work we perform simulations of realistic microgels in an explicit solvent. We identify a suitable model which fully captures the main features of the implicit model and further provides information on the solvent uptake by the interior of the microgel network and on its role in the collapse kinetics. These results pave the way for addressing problems where solvent effects are dominant, such as the case of microgels at liquid-liquid interfaces.

Simulations of Interfacial Tension of Liquid-Liquid Ternary Mixtures Using Optimized Parametrization for Coarse-Grained Models

In this work, liquid-liquid systems are studied by means of coarse-grained Monte Carlo simulations (CG-MC) and Dissipative Particle Dynamics (DPD). A methodology is proposed to reproduce liquid-liquid equilibrium (LLE) and to provide variation of interfacial tension (IFT), as a function of the solute concentration. A key step is the parametrization method based on the use of the Flory-Huggins parameter between DPD beads to calculate solute/solvent interactions. Parameters are determined using a set of experimental compositional data of LLE, following four different approaches. These approaches are evaluated, and the results obtained are compared to analyze advantages/disadvantages of each one. These methodologies have been compared through their application on six systems: water/benzene/1,4-dioxane,water/chloroform/acetone, water/benzene/acetic acid, water/benzene/2-propanol, water/hexane/acetone, and water/hexane/2-propanol. CG-MC simulations in the Gibbs (NVT) ensemble have been used to check the validity of parametrization approaches for LLE reproduction. Then, CG-MC simulations in the osmotic (μ_soluteN_solventP _zzT) ensemble were carried out considering the two liquid phases with an explicit interface. This step allows one to work at the same bulk concentrations as the experimental data by imposing the precise bulk phase compositions and predicting the interface composition. Finally, DPD simulations were used to predict IFT values for different solute concentrations. Our results on variation of IFT with solute concentration in bulk phases are in good agreement with experimental data, but some deviations can be observed for systems containing hexane molecules.

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.

Coarse-Graining the Liquid-Liquid Interfaces with the MARTINI Force Field: How Is the Interfacial Tension Reproduced?

We report two-phase coarse-grained (CG) simulations of organic-water liquid-liquid interfaces with the MARTINI force field. We discuss the ability of the CG force field to predict quantitatively the interfacial tension of alkanes-water, benzene-water, chloroform-water, and alcohol-water systems. The performance of the prediction of the interfacial tension is evaluated through its dependence on temperature and alkane length. This study contributes to the challenging discussion about the robustness and the transferability of the MARTINI force field to interfacial properties. We have also used the distributions of the molecules along the direction normal to the interface to investigate the composition of the interfacial region and to compare the simulated densities of the coexisting phases with experiments.

Understanding the interfacial behavior in isopycnic Lennard-Jones mixtures by computer simulations

The physical characterization of the singular interfacial behavior of heterogeneous fluid systems is a very important step in preliminary stages of the design process, and also in the subsequent procedures for the determination of the optimal operating conditions. Molar isopycnicity or molar density inversion is a special case of phase equilibrium behavior that directly affects the relative position of phases in heterogeneous mixtures, without being affected by gravitational fields. This work is dedicated to characterize the impact of molar density inversion on the interfacial properties of Lennard-Jones binary mixtures. The results and specific trends of the molar density inversion phenomena on the peculiar calculated composition profiles across the interface and interfacial tensions are explored by using canonical molecular dynamics simulations of the Lennard-Jones binary mixtures. Our results show that the density inversion causes drastic changes in the density profiles of the mixtures. In particular, symmetrical and equal-sized Lennard-Jones mixtures always exhibit desorption along the interfacial zone, i.e. the interfacial concentration profiles show a relative minimum at the interface of the total density profiles that increases when the dispersive energy parameter (ε(ij)) between unlike species decreases. However, as the asymmetry of the Lennard-Jones mixtures increases (σ(i) ≠ σ(j)), the concentration profiles display a relative maximum at the interface, which implies the adsorption of the total density profiles along the interfacial zone.