The surface tension of Martini 3 water mixtures

The Martini model, a coarse-grained forcefield for biomolecular simulations, has experienced a vast increase in popularity in the past decade. Its building-block approach balances computational efficiency with high chemical specificity, enabling the simulation of organic and inorganic molecules. The modeling of coarse-grained beads as Lennard-Jones particles poses challenges for the accurate reproduction of liquid-vapor interfacial properties, which are crucial in various applications, especially in the case of water. The latest version of the forcefield introduces refined interaction parameters for water beads, tackling the well-known artifact of Martini water freezing at room temperature. In addition, multiple sizes of water beads are available for simulating the solvation of small cavities, including the smallest pockets of proteins. This work focuses on studying the interfacial properties of Martini water, including surface tension and surface thickness. Employing the test-area method, we systematically compute the liquid-vapor surface tension across various combinations of water bead sizes and for temperatures from 300 to 350 K. These findings are of interest to the Martini community as they allow users to account for the low interfacial tension of Martini water by properly adjusting observables computed via coarse-grained simulations to allow for accurate matching against all-atom or experimental results. Surface tension data are also interpreted in terms of local enrichment of the various mixture components at the liquid-vapor interface by means of Gibbs' adsorption formalism. Finally, the critical scaling of the Martini surface tension with temperature is reported to be consistent with the critical exponent of the 3D Ising universality class.

Influence of molecular parameters on the representativeness of interfacial properties of simple fluids

New parameterizations for the Lennard-Jones 12/6 potential capable of reproducing the vapor pressure and surface tension with sufficient precision, but not the liquid-vapor equilibrium densities for the case of simple fluids that include Ar, Kr, Xe, Ne, and CH4 are presented in this work. These results are compared with those derived from the family of Mie(n, 6) potentials, which adequately reproduce the coexistence curve and the vapor pressure, leaving aside the surface tension. In addition, a detailed analysis is presented on different parameterizations and methodologies, which have been developed in recent decades to estimate the interfacial properties of interest here for simple fluids, such as argon, which is a molecule that is, in principle, "simple" to study but that clearly reveals the enormous discrepancy between the results reported in the literature throughout these years. These facts undoubtedly reveal one of the fundamental problems in the context of molecular thermodynamics of fluids: reproducing different thermodynamic properties with sufficient precision from a single set of free parameters for some interaction potential. In order to show the scope of the parameterizations presented for the Lennard-Jones model, they were successfully applied to the case of binary mixtures, which included Ar-Kr, Ar-CH4, and Xe-Kr. Finally, and with the aim of showing a possible solution to the problem posed in this research, results of the same interfacial properties above mentioned for argon and methane were presented in this work by using a set of molecular interactions, called ANC2s, whose flexibility allowed to reproduce the experimental evidence with just one parameterization. The results reported in this work were generated using molecular dynamics simulations.

Self-Amputating and Interfusing Machines

Biological organisms exhibit phenomenal adaptation through morphology-shifting mechanisms including self-amputation, regeneration, and collective behavior. For example, reptiles, crustaceans, and insects amputate their own appendages in response to threats. Temporary fusion between individuals enables collective behaviors, such as in ants that temporarily fuse to build bridges. The concept of morphological editing often involves the addition and subtraction of mass and can be linked to modular robotics, wherein synthetic body morphology may be revised by rearranging parts. This work describes a reversible cohesive interface made of thermoplastic elastomer that allows for strong attachment and easy detachment of distributed soft robot modules without direct human handling. The reversible joint boasts a modulus similar to materials commonly used in soft robotics, and can thus be distributed throughout soft robot bodies without introducing mechanical incongruities. To demonstrate utility, the reversible joint is implemented in two embodiments: a soft quadruped robot that self-amputates a limb when stuck, and a cluster of three soft-crawling robots that fuse to cross a land gap. This work points toward future robots capable of radical shape-shifting via changes in mass through autotomy and interfusion, as well as highlights the crucial role that interfacial stiffness change plays in autotomizable biological and artificial systems.

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.

Fundamental Aspects of Phase-Separated Biomolecular Condensates

Biomolecular condensates, formed through phase separation, are upending our understanding in much of molecular, cell, and developmental biology. There is an urgent need to elucidate the physicochemical foundations of the behaviors and properties of biomolecular condensates. Here we aim to fill this need by writing a comprehensive, critical, and accessible review on the fundamental aspects of phase-separated biomolecular condensates. We introduce the relevant theoretical background, present the theoretical basis for the computation and experimental measurement of condensate properties, and give mechanistic interpretations of condensate behaviors and properties in terms of interactions at the molecular and residue levels.

Importance of the Electrostatic Correlations in Surface Tension of Hydrated Reline Deep Eutectic Solvent from Combined Experiments and Molecular Dynamics Simulations

In this study, the surface tension and the structure of hydrated reline are investigated by using diverse methods. Initially, the surface tension displays a nonlinear pattern as water content increases, decreasing until reaching 45 wt %, then gradually matching that of pure water. This fluctuation is associated with strong electrostatic correlations present in pure reline, which decrease as more water is added. Changes in surface tension reflect a shift from charge layering in pure reline to an increased interfacial hydrogen bonding as the water content rises. This shift causes the segregation of urea molecules into the bulk phase and a gradual anchoring of water molecules to the air-reline interface. An interesting observation is the antisurfactant effect, where heightened interfacial anchoring results in an unexpected increase in real contribution of surface tension. This, along with weakened electrostatic correlations beyond 45 wt % due to reinforced interfacial hydrogen bonding, contributes to the complex behavior of surface tension observed in this study.

Multiscale simulation of fluids: coupling molecular and continuum

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.

Thermodynamics of liquid and fluid mixtures from the kinetic Monte Carlo viewpoint

In this study a binary mixture is modelled in a uniform simulation cell at various temperatures using an extended version of the grand canonical kinetic Monte Carlo (GC-kMC) method. The main goal of this study is to consider the thermodynamic properties of binary liquids, gases, and gas-liquid mixtures from a more general point of view than that applied to the particular case of vapour-liquid equilibrium when the pressure and partial chemical potentials in coexisting phases are the same. Particular attention is paid to thermodynamic functions such as chemical potentials, Gibbs free energy and entropy. For the pair potential of unlike molecules a more universal scheme is proposed in comparison with the Lorentz-Berthelot combining rule. The approach is tested on an Ar-Kr mixture in a wide range of temperatures. In all cases, the obtained values of chemical potentials, pressure and internal energy for the entire set of component densities and temperature fully satisfy the Gibbs-Duhem equation with a high degree of accuracy. For the case of vapour-liquid equilibrium, the developed approach made it possible to reproduce the experimental pressure-composition diagrams with the highest accuracy ever achieved in the literature. Despite the fact that the Ar-Kr mixture, according to Raoult's law, is close to an ideal system, it was found that the partial pressures in the liquid phase or in a dense supercritical gas mixture are non-linear functions of the composition, and the partial pressure of the heavier component (Kr) can even be negative.

Phase behaviour of coarse-grained fluids

Soft condensed matter structures often challenge us with complex many-body phenomena governed by collective modes spanning wide spatial and temporal domains. In order to successfully tackle such problems, mesoscopic coarse-grained (CG) statistical models are being developed, providing a dramatic reduction in computational complexity. CG models provide an intermediate step in the complex statistical framework of linking the thermodynamics of condensed phases with the properties of their constituent atoms and molecules. These allow us to offload part of the problem to the CG model itself and reformulate the remainder in terms of reduced CG phase space. However, such exchange of pawns to chess pieces, or 'Hamiltonian renormalization', is a radical step and the thermodynamics of the primary atomic and CG models could be quite distinct. Here, we present a comprehensive study of the phase diagram including binodal and interfacial properties of a dissipative particle dynamics (DPD) model, extended to include finite-range attraction to support the liquid-gas equilibrium. Despite the similarities with the atomic model potentials, its phase envelope is markedly different featuring several anomalies such as an unusually broad liquid range, change in concavity of the liquid coexistence branch with variation of the model parameters, volume contraction on fusion, temperature of maximum density in the liquid phase and negative thermal expansion in the solid phase. These results provide new insight into the connection between simple potential models and complex emergent condensed matter phenomena.

Bubble-Mediated Large-Scale Hierarchical Assembly of Ultrathin Pt Nanowire Network Monolayer at Gas/Liquid Interfaces

Extensive macroscale two-dimensional (2-D) platinum (Pt) nanowire network (NWN) sheets are created through a hierarchical self-assembly process with the aid of biomolecular ligands. The Pt NWN sheet is assembled from the attachment growth of 1.9 nm-sized 0-D nanocrystals into 1-D nanowires featuring a high density of grain boundaries, which then interconnect to form monolayer network structures extending into centimeter-scale size. Further investigation into the formation mechanism reveals that the initial emergence of NWN sheets occurs at the gas/liquid interfaces of the bubbles produced by sodium borohydride (NaBH₄) during the synthesis process. Upon the rupture of these bubbles, an exocytosis-like process releases the Pt NWN sheets at the gas/liquid surface, which subsequently merge into a continuous monolayer Pt NWN sheet. The Pt NWN sheets exhibit outstanding oxygen reduction reaction (ORR) activities, with specific and mass activities 12.0 times and 21.2 times greater, respectively, than those of current state-of-the-art commercial Pt/C electrocatalysts.

Impact of the Force Field on the Calculation of Density and Surface Tension of Epoxy-Resins

The molecular simulation of interfacial systems is a matter of debate because of the choice of many input parameters that can affect significantly the performance of the force field of reproducing the surface tension and the coexisting densities. After developing a robust methodology for the calculation of the surface tension on a Lennard-Jones fluid, we apply it with different force fields to calculate the density and surface tension of pure constituents of epoxy resins. By using the model that best reproduces the experimental density and surface tension, we investigate the impact of composition in mass fraction on uncured epoxy resins and the effects of degree of cross-linking on cured resins.

An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation

Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.

An extended Stokes-Einstein model for condensed ionic water structures with topological complexity

'What is the structure of water?' This has been a perplexing question for a long time and water structure with various phases is a great topic of research interest. Topological complexity generally occurs because hydrophilic ions strongly influence the size and shape of condensed water structures owing to their kosmotropic and chaotropic transitions. In this study, an extended Stokes-Einstein model incorporating Flory-Huggins free energy equation is proposed to describe the constitutive relationship between dynamic diffusion and condensed water structure with a topological complexity. The newly developed model provides a geometrical strategy of end-to-end distance and explores the constitutive relationship between condensed ionic water structures and their dynamic diffusion behaviors. A free-energy function is then formulated to study thermodynamics in electrolyte aqueous solution, in which the condensed ionic water structures undergo topologically complex changes. Finally, effectiveness of the proposed model is verified using both molecular dynamics simulations and experimental results reported in literature.

A simple theory for interfacial properties of dilute solutions

Recent studies suggest that cosolute mixtures may exert significant non-additive effects upon protein stability. The corresponding liquid–vapor interfaces may provide useful insight into these non-additive effects. Accordingly, in this work, we relate the interfacial properties of dilute multicomponent solutions to the interactions between solutes. We first derive a simple model for the surface excess of solutes in terms of thermodynamic observables. We then develop a lattice-based statistical mechanical perturbation theory to derive these observables from microscopic interactions. Rather than adopting a random mixing approximation, this dilute solution theory (DST) exactly treats solute–solute interactions to lowest order in perturbation theory. Although it cannot treat concentrated solutions, Monte Carlo (MC) simulations demonstrate that DST describes the interactions in dilute solutions with much greater accuracy than regular solution theory. Importantly, DST emphasizes a fundamental distinction between the “intrinsic” and “effective” preferences of solutes for interfaces. DST predicts that three classes of solutes can be distinguished by their intrinsic preference for interfaces. While the surface preference of strong depletants is relatively insensitive to interactions, the surface preference of strong surfactants can be modulated by interactions at the interface. Moreover, DST predicts that the surface preference of weak depletants and weak surfactants can be qualitatively inverted by interactions in the bulk. We also demonstrate that DST can be extended to treat surface polarization effects and to model experimental data. MC simulations validate the accuracy of DST predictions for lattice systems that correspond to molar concentrations.

Molecular interactions at the metal-liquid interfaces

We reported molecular simulations of the interactions among water, an epoxy prepolymer diglycidic ether of bisphenol A (DGEBA), and a hardener isophorone diamine (IPDA) on an aluminum surface. This work proposes a comprehensive thermodynamic characterization of the adhesion process from the calculation of different interfacial tensions. The cross-interactions between the atoms of the metal surface and different molecules are adjusted so as to reproduce the experimental work of adhesion. Water nanodroplets on the metal surface are then simulated to predict their contact angle. Liquid-vapor surface tensions of the epoxy prepolymer (DGEBA) and hardener (IPDA) and the solid-vapor surface tension of the aluminum surface are also calculated to provide the solid-liquid interfacial tension that remains very difficult to obtain from the mechanical definition.

Insulin Adsorption onto PE and PVC Tubings

We investigate the adsorption of insulin onto PE and PVC materials by using HPLC measurements and computer simulations. We interpret the experiments by calculating the Gibbs free energy profiles during the adsorption process. The values of free energy of adsorption show a good agreement with the experimental measurements. The adsorption of insulin onto the different materials is characterized through the conformational changes with respect to its conformation in water and the interfacial regions, which are described by specific arrangements of polymer chains, water, insulin, and plasticizer molecules.

Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?

In principle, the answer to the posed titular question is undoubtedly "yes." But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the "cut-and-shift" version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I_h. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.

The Intrinsic Fragility of the Liquid-Vapor Interface: A Stress Network Perspective

The evolution of the liquid-vapor interface of a Lennard-Jones fluid is examined with molecular dynamics simulations using the intrinsic sampling method. Results suggest clear damping of the intrinsic profiles with increasing temperature. Investigating the surface stress distribution, we have identified a linear variation of the space-filling nature (fractal dimension) of the stress clusters at the intrinsic surface with increasing surface tension or, equivalently, with decreasing temperature. A percolation analysis of these stress networks indicates that the stress field is more disjointed at higher temperatures. This leads to more fragile (or poorly connected) interfaces which result in a reduction in surface tension.

Osmotic Method for Calculating Surface Pressure of Monolayers in Molecular Dynamics Simulations

Surface pressure is a fundamental thermodynamic property related to the activity of molecules at interfaces. In molecular simulations, it is typically calculated from its definition: the difference between the surface tension of the air-water and air-surfactant interfaces. In this Letter, we show how to connect the surface pressure with a two-dimensional osmotic pressure and how to take advantage of this analogy to obtain a practical method of calculating surface pressure-area isotherms in molecular simulation. As a proof-of-concept, compression curves of zwitterionic and ionic surfactant monolayers were obtained using the osmotic approach and the curves were compared with the ones from the traditional pressure tensor-based scheme. The results shown an excellent agreement between both alternatives. Advantageously, the osmotic approach is simple to use and allows to obtain the surface pressure-area isotherm on the fly with a single simulation using equilibration stages.

Characterization of axial dispersion in vertical helical coil for gas-liquid-liquid flow at low Reynolds numbers

Gas-Liquid-Liquid (GLL) slug flow reactors offer several advantages like higher interfacial area, excellent mass transfer, and lower backmixing. The mesoscale (diameter ~ few mm) helical coiled reactors operating in slug...

Mesoscale perspective on the Tolman length

We demonstrate that the multiphase Shan-Chen lattice Boltzmann method (LBM) yields a curvature dependent surface tension σ as computed from three-dimensional hydrostatic droplets and bubbles simulations. Such curvature dependence is routinely characterized, at first order, by the so-called Tolman length δ. LBM allows one to precisely compute σ at the surface of tension R_{s} and determine the Tolman length from the coefficient of the first order correction. The corresponding values of δ display universality for different equations of state, following a power-law scaling near the critical temperature. The Tolman length has been studied so far mainly via computationally demanding Molecular Dynamics simulations or by means of Density Functional Theory approaches playing a pivotal role in extending Classical Nucleation Theory. The present results open a hydrodynamic-compliant mesoscale arena, in which the fundamental role of the Tolman length, alongside real-world applications to cavitation phenomena, can be effectively tackled. All the results can be independently reproduced through the "idea.deploy" framework.

Surface morphology effects on clathrate hydrate wettability

HYPOTHESIS

Clathrate hydrates preferentially form at interfaces; hence, wetting properties play an important role in their formation, growth, and agglomeration. Experimental evidence suggests that the hydrate preparation process can strongly affect contact angle measurements, leading to the different results reported in the literature. These differences hamper technological progress. We hypothesize that changes in hydrate surface morphologies are responsible for the wide variation of contact angles reported in the literature.

EXPERIMENTS

Experimental testing of our hypothesis is problematic due to the preparation history of hydrates on their surface properties, and the difficulties in advanced surface characterization. Thus, we employ molecular dynamics simulations, which allow us to systematically change the interfacial features and the system composition. Implementing advanced algorithms, we quantify fundamental thermodynamic properties to validate our observations.

FINDINGS

We achieve excellent agreement with experimental observations for both atomically smooth and rough hydrate surfaces. Our results suggest that contact line pinning forces, enhanced by surface heterogeneity, are accountable for altering water contact angles, thus explaining the differences among reported experimental data. Our analysis and molecular level insights help interpret adhesion force measurements and yield a better understanding of the agglomeration between hydrate particles, providing a microscopic tool for advancing flow assurance applications.

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.

Cassie State Stability and Gas Restoration Capability of Superhydrophobic Surfaces with Truncated Cone-Shaped Pillars

The gas layer stability on superhydrophobic surfaces and gas restoration on the immersed superhydrophobic surfaces have been great challenges for their practical applications in recent years. Inspired by the naturally existing mushroom-like super-repellent superhydrophobic patterns, we choose superhydrophobic surfaces with truncated cone-shaped pillars as our research objects to tackle such challenges by tuning their geometrical parameters. We perform molecular dynamics simulations to investigate the Cassie-Wenzel transition under external pressure and the Wenzel-Cassie transition due to underwater spreading of compressed bubbles. Theories based on the Young-Laplace equation and total free-energy variation are developed to explore the influence of geometrical parameters of pillars on the pressure resistance and underwater gas restoration, which is in good agreement with simulation results. These simulation results and theoretical analysis suggest that cork-shaped pillars, analogous to the surface structures of natural organisms like springtails and Salvinia leaves, can be super-repellent to the liquid and favorable for the gas spreading process. Our study provides theoretical guidance for the design of superhydrophobic surfaces with both Cassie state stability and gas restoration capability.

Structure and Contact Angle in Sessile Droplets of Binary Mixtures of Lennard-Jones Chains: A Molecular Dynamics Study

Molecular dynamics simulations were carried out to investigate cylindrical droplets consisting of binary mixtures of Lennard-Jones (LJ) fluids in contact with a solid substrate. The droplets are composed of mixtures of the monomeric LJ fluid plus linear-tangent chains of 2, 10, 20, and 30 segments per chain that interact through a harmonic potential and the spherically truncated and shifted potential Lennard-Jones. The solid surface was modeled as a semi-infinite platinum substrate with an FCC structure that interacts with the fluid by means of a LJ 9-3 potential. We place emphasis on the effect of mixing a monomeric LJ fluid with heavy components on the contact angle and on the droplet structure, especially in the liquid-solid region. The density profiles of the droplets reveal a strong discrete layering of the fluid in the vicinity of the solid-liquid interface. The layering is more pronounced at low temperatures and for mixtures of short chains (symmetric mixtures). The ordering of the fluid was much less intense for fluids of long chains (asymmetric mixtures), and some cases even show gas enrichment at the solid-liquid interface. Enrichment at the vapor-liquid interfaces and density inversion can also be observed. However, these effects are not as marked as in planar interfaces. The contact angle between the droplet and the substrate is calculated by fitting an ellipse to the vapor-liquid interface defined by the Gibbs dividing surface. In general, an increment in the concentration of the heavy component and a reduction of the temperature resulted in an increase of the contact angle, which in turn disfavored the wetting of the droplet.

A novel hybrid method for the calculation of methane hydrate-water interfacial tension along the three-phase (hydrate-liquid water-vapor) equilibrium line

We use a novel hybrid method to explore the temperature dependence of the solid-liquid interfacial tension of a system that consists of solid methane hydrate and liquid water. The calculated values along the three-phase (hydrate-liquid water-vapor) equilibrium line are obtained through the combination of available experimental measurements and computational results that are based on approaches at the atomistic scale, including molecular dynamics and Monte Carlo. An extensive comparison with available experimental and computational studies is performed, and a critical assessment and re-evaluation of previously reported data is presented.

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.

Macromolecular regulators have matching effects on the phase equilibrium and interfacial tension of biomolecular condensates

The interfacial tension of phase-separated biomolecular condensates affects their fusion and multiphase organization, and yet how this important property depends on the composition and interactions of the constituent macromolecules is poorly understood. Here we use molecular dynamics simulations to determine the interfacial tension and phase equilibrium of model condensate-forming systems. The model systems consist of binary mixtures of Lennard-Jones particles or chains of such particles. We refer to the two components as drivers and regulators; the former has stronger self-interactions and hence a higher critical temperature (T_c ) for phase separation. In previous work, we have shown that, depending on the relative strengths of driver-regulator and driver-driver interactions, regulators can either promote or suppress phase separation (i.e., increase or decrease T_c ). Here we find that the effects of regulators on T_c quantitatively match the effects on interfacial tension (γ). This important finding means that, when a condensate-forming system experiences a change in macromolecular composition or a change in intermolecular interactions (e.g., by mutation or posttranslational modification, or by variation in solvent conditions such as temperature, pH, or salt), the resulting change in T_c can be used to predict the change in γ and vice versa. We also report initial results showing that disparity in intermolecular interactions drives multiphase coexistence. These findings provide much needed guidance for understanding how biomolecular condensates mediate cellular functions.

On the Gibbs-Thomson equation for the crystallization of confined fluids

The Gibbs-Thomson (GT) equation describes the shift of the crystallization temperature for a confined fluid with respect to the bulk as a function of pore size. While this century old relation is successfully used to analyze experiments, its derivations found in the literature often rely on nucleation theory arguments (i.e., kinetics instead of thermodynamics) or fail to state their assumptions, therefore leading to similar but different expressions. Here, we revisit the derivation of the GT equation to clarify the system definition, corresponding thermodynamic ensemble, and assumptions made along the way. We also discuss the role of the thermodynamic conditions in the external reservoir on the final result. We then turn to numerical simulations of a model system to compute independently the various terms entering in the GT equation and compare the predictions of the latter with the melting temperatures determined under confinement by means of hyper-parallel tempering grand canonical Monte Carlo simulations. We highlight some difficulties related to the sampling of crystallization under confinement in simulations. Overall, despite its limitations, the GT equation may provide an interesting alternative route to predict the melting temperature in large pores using molecular simulations to evaluate the relevant quantities entering in this equation. This approach could, for example, be used to investigate the nanoscale capillary freezing of ionic liquids recently observed experimentally between the tip of an atomic force microscope and a substrate.

Associated molecular liquids at the graphene monolayer interface

We report molecular simulations of the interaction between a graphene sheet and different liquids such as water, ethanol, and ethylene glycol. We describe the structural arrangements at the graphene interface in terms of density profiles, number of hydrogen bonds (HBs), and local structuration in neighboring layers close to the surface. We establish the formation of a two-dimensional HB network in the layer closest to the graphene. We also calculate the interfacial tension of liquids with a graphene monolayer and its profile along the direction normal to the graphene to rationalize and quantify the strengthening of the intermolecular interactions in the liquid due to the presence of the surface.

Interfacial properties of binary mixtures of Lennard-Jones chains in planar interfaces by molecular dynamics simulation

Binary mixtures of fully flexible linear tangent chains composed of bonded Lennard-Jones interaction sites (monomers) were studied using the molecular dynamics simulation in the NVT ensemble. Their interfacial properties were investigated in planar interfaces by direct simulation of an explicit liquid film in equilibrium with its vapor. A method for the calculation of long-range interactions in inhomogeneous fluids was implemented to take into account the potential truncation effects. Surface tension and the pressure tensor were calculated via the classical Irving-Kirkwood method; vapor pressure, orthobaric densities, density profiles, and Gibbs relative adsorption of the volatile component with respect to the heavy component were also obtained. The properties were studied as a function of the temperature, molar concentration of the heavy component, and the asymmetry of the mixture. According to the results of this work, the temperature loses influence on the surface tension, vapor pressure, and Gibbs relative adsorption curves as the molecular length of the heavy component increases. This suggests that the universal behavior observed in pure fluids of Lennard-Jones chains also holds for binary mixtures. The contribution of the long-range interactions turned out to account for about 60%, 20%, and 10% of the surface tension, vapor pressure, and orthobaric density final values, respectively. This contribution was even larger at high temperatures and for large molecules. Strong enrichment of the volatile component at the interface was observed in the asymmetric mixtures. One of these mixtures even showed a barotropic effect at elevated pressures and a class III phase behavior.

Liquid/liquid interface in periodic boundary condition

We study how surface phenomena can change the interface geometry in liquid-liquid two-phase systems with periodic boundary conditions. Without any curvature effect on surface tension, planar (slab), cylindrical, and spherical structures are successively obtained as a function of the total composition and elongation of the box, in accordance with molecular dynamics simulations for a water/heptane system. The curvature effects described by Tolman relationship desymmetrize the phase diagram by stabilizing a concavity but it leads to inconsistencies with high curvature. Helfrich model partially resolves this and predicts the possible presence of shells reflecting a frustrated system.

Bottom-up coarse-grained models for external fields and interfaces

Bottom-up coarse-grained (CG) models accurately describe the structure of homogeneous systems but sometimes provide limited transferability and a poor description of thermodynamic properties. Consequently, inhomogeneous systems present a severe challenge for bottom-up models. In this work, we examine bottom-up CG models for interfaces and inhomogeneous systems. We first analyze the effect of external fields upon the many-body potential of mean force. We also demonstrate that the multiscale CG (MS-CG) variational principle for modeling the external field corresponds to a generalization of the first Yvon-Born-Green equation. This provides an important connection with liquid state theory, as well as physical insight into the structure of interfaces and the resulting MS-CG models. We then develop and assess MS-CG models for a film of liquid methanol that is adsorbed on an attractive wall and in coexistence with its vapor phase. While pair-additive potentials provide unsatisfactory accuracy and transferability, the inclusion of local-density (LD) potentials dramatically improves the accuracy and transferability of the MS-CG model. The MS-CG model with LD potentials quite accurately describes the wall-liquid interface, the bulk liquid density, and the liquid-vapor interface while simultaneously providing a much improved description of the vapor phase. This model also provides an excellent description of the pair structure and pressure-density equation of state for the bulk liquid. Thus, LD potentials hold considerable promise for transferable bottom-up models that accurately describe the structure and thermodynamic properties of both bulk and interfacial systems.

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.