Molecular simulations of the n -alkane liquid-vapor interface: interfacial properties and their long range corrections

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.

A perspective on the microscopic pressure (stress) tensor: History, current understanding, and future challenges

The pressure tensor (equivalent to the negative stress tensor) at both microscopic and macroscopic levels is fundamental to many aspects of engineering and science, including fluid dynamics, solid mechanics, biophysics, and thermodynamics. In this Perspective, we review methods to calculate the microscopic pressure tensor. Connections between different pressure forms for equilibrium and nonequilibrium systems are established. We also point out several challenges in the field, including the historical controversies over the definition of the microscopic pressure tensor; the difficulties with many-body and long-range potentials; the insufficiency of software and computational tools; and the lack of experimental routes to probe the pressure tensor at the nanoscale. Possible future directions are suggested.

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.

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.

Microscopic Pressure Tensor in Cylindrical Geometry: Pressure of Water in a Carbon Nanotube

The microscopic pressure tensor plays an important role in understanding the mechanical stability, transport, and high-pressure phenomena of confined phases. The lack of an exact formulation to account for the long-range Coulombic contribution to the local pressure tensor in cylindrical geometries prevents the characterization of molecular fluids confined in cylindrical pores. To address this problem, we first derive the local cylindrical pressure tensor for Lennard-Jones fluids based on the Harasima (H) definition, which is expected to be compatible with the Ewald summation method. The test of the H-definition pressure equations in a homogeneous system shows that the radial and azimuthal pressure have unphysical radial dependence near the origin, while the axial pressure gives physically meaningful values. We propose an alternative contour definition that is more appropriate for cylindrical geometry and show that it leads to physically realistic results for all three pressure tensor components. With this definition, the radial and azimuthal pressures are of Irving-Kirkwood (IK) type, and the axial pressure is of Harasima type. Because of the practical interest in the axial pressure, we develop a Harasima/Ewald (H/E) method for calculating the long-range Coulombic contribution to the local axial pressure for rigid molecules. As an application, the axial pressure profile of water inside and outside a (20, 20) single-wall carbon nanotube is determined. The H/E method is compared to the IK method, which assumes a spherically truncated Coulombic potential. Detailed analysis of the pressure profile by both methods shows that the water confined in the nanotube is in a stretched state overall in the axial direction.

The pair distribution function in the planar gas-liquid interface: Application to the calculation of the surface tension

A Monte Carlo simulation is used to calculate the pair distribution function g⁽²⁾r₁,r₂ for a planar gas-liquid interface. Due to the cylindrical symmetry of the system, g⁽²⁾ can be stored as a three-dimensional array that can be readily manipulated and used to calculate the surface tension and the single atom density profile directly. The consistency and accuracy of our calculation of g⁽²⁾(r₁, r₂) is demonstrated by a calculation of the single atom density through the first Born-Green-Yvon equation. We show that the surface tension calculated directly from the pair distribution function and from other well-established routes is completely consistent. In the case of the gas-liquid interface for argon modeled with an explicit inclusion of the three-body forces, an accurate pair distribution can be used to estimate the long-range contribution to the three-body part of the surface tension. A detailed analysis of this correction, its dependence on the three-body cutoff, and its overall contribution to the surface tension are presented.

Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields

The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.

Calculation of the interfacial tension of the graphene-water interaction by molecular simulations

We report the calculation of the solid-liquid interface tension of the graphene-water interaction by using molecular simulations. Local profiles of the interfacial tension are given through the mechanical and thermodynamic definitions. The dependence of the interfacial tension on the graphene area is investigated by applying both reaction field and Ewald summation techniques. The structure of the interfacial region close to the graphene sheet is analyzed through the profiles of the density and hydrogen bond number and the orientation of the water molecules. We complete this study by plotting the profiles of the components of the pressure tensor calculated by the Ewald summation and reaction field methods. We also investigate the case of a reaction field version consisting in applying a damped shifted force in the case of the calculation of the pressure components.

On the equilibrium contact angle of sessile liquid drops from molecular dynamics simulations

We present a new methodology to estimate the contact angles of sessile drops from molecular simulations by using the Gaussian convolution method of Willard and Chandler [J. Phys. Chem. B 114, 1954-1958 (2010)] to calculate the coarse-grained density from atomic coordinates. The iso-density contour with average coarse-grained density value equal to half of the bulk liquid density is identified as the average liquid-vapor (LV) interface. Angles between the unit normal vectors to the average LV interface and unit normal vector to the solid surface, as a function of the distance normal to the solid surface, are calculated. The cosines of these angles are extrapolated to the three-phase contact line to estimate the sessile drop contact angle. The proposed methodology, which is relatively easy to implement, is systematically applied to three systems: (i) a Lennard-Jones (LJ) drop on a featureless LJ 9-3 surface; (ii) an SPC/E water drop on a featureless LJ 9-3 surface; and (iii) an SPC/E water drop on a graphite surface. The sessile drop contact angles estimated with our methodology for the first two systems are shown to be in good agreement with the angles predicted from Young's equation. The interfacial tensions required for this equation are computed by employing the test-area perturbation method for the corresponding planar interfaces. Our findings suggest that the widely adopted spherical-cap approximation should be used with caution, as it could take a long time for a sessile drop to relax to a spherical shape, of the order of 100 ns, especially for water molecules initiated in a lattice configuration on a solid surface. But even though a water drop can take a long time to reach the spherical shape, we find that the contact angle is well established much faster and the drop evolves toward the spherical shape following a constant-contact-angle relaxation dynamics. Making use of this observation, our methodology allows a good estimation of the sessile drop contact angle values even for moderate system sizes (with, e.g., 4000 molecules), without the need for long simulation times to reach the spherical shape.

Importance of the tail corrections on surface tension of curved liquid-vapor interfaces

We report molecular simulations of the liquid-vapor cylindrical interface of methane. We apply the truncated Lennard-Jones potential and specific long-range corrections for the surface tension developed especially for cylindrical interfaces. We investigate the impact of the cutoff on the radial density profile, the intrinsic and long-range correction parts to the surface tension, and Tolman length. We also study the curvature dependence of the surface tension as a function of the cutoff used. In this work we shed light that both density and Tolman length are cutoff-dependent whereas the total surface tension is slightly curvature and cutoff dependent.

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.

Ultrafast diffusion of Ionic Liquids Confined in Carbon Nanotubes

Over the past decade many works have focused on various aspects of the dynamics of liquids confined at the nanoscale such as e.g. water flow enhancement through carbon nanotubes (CNTs). Transport of room temperature ionic liquids (RTILs) through various nanochannels has also been explored and some conflicting findings about their translational dynamics have been reported. In this work, we focus on translational dynamics of RTILs confined in various CNTs. By means of molecular dynamics simulations we highlight a substantially enhanced diffusion of confined RTILs with an increase up to two orders of magnitude with respect to bulk-phase properties. This ultrafast diffusion of RTILs inside CNTs is shown to result from the combination of various factors such as low friction, molecular stacking, size, helicity, curvature and cooperative dynamics effects.

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.

Local pressure components and interfacial tensions of a liquid film in the vicinity of a solid surface with a nanometer-scale slit pore obtained by the perturbative method

A classical molecular dynamics simulation was conducted for a liquid-solid interfacial system with a nanometer-scale slit pore in order to reveal local thermodynamic states: local pressure components and interfacial tensions of a liquid film in the vicinity of the slit. The simulation also examined the transition mechanism between the two states of the liquid film: (a) liquid film on the slit and (b) liquid film in the slit, based on the local thermodynamic quantities from a molecular point of view. An instantaneous expression of the local pressure components and interfacial tensions, which is based on a volume perturbation, was presented to investigate time-dependent phenomena in molecular dynamics simulations. The interactions between the particles were described by the 12-6 Lennard-Jones potential, and effects of the fluid-solid interaction intensity on the local pressure components and interfacial tensions of the fluid in the vicinity of the slit were examined in detail by the presented perturbative method. The results revealed that the local pressure components tangential to the solid surface in the vicinity of the 1st fluid layer from the solid surface are different in a two dimensional plane, and the difference became pronounced in the vicinity of the corner of the slit, for cases where the fluid-solid interaction intensities are relatively strong. The results for the local interfacial tensions of the fluid inside the slit suggested that the local interfacial tensions in the vicinity of the 2nd and 3rd layers of the solid atoms from the entrance of the slit act as a trigger for the transition between the two states under the influence of a varying fluid-solid interaction.

Temperature dependence of local solubility of hydrophobic molecules in the liquid-vapor interface of water

One important aspect of the hydrophobic effect is that solubility of small, nonpolar molecules in liquid water decreases with increasing temperature. We investigate here how the characteristic temperature dependence in liquid water persists or changes in the vicinity of the liquid-vapor interface. From the molecular dynamics simulation and the test-particle insertion method, the local solubility Σ of methane in the liquid-vapor interface of water as well as Σ of nonpolar solutes in the interface of simple liquids are calculated as a function of the distance z from the interface. We then examine the temperature dependence of Σ under two conditions: variation of Σ at fixed position z and that at fixed local solvent density around the solute molecule. It is found that the temperature dependence of Σ at fixed z depends on the position z and the system, whereas Σ at fixed local density decreases with increasing temperature for all the model solutions at any fixed density between vapor and liquid phases. The monotonic decrease of Σ under the fixed-density condition in the liquid-vapor interface is in accord with what we know for the solubility of nonpolar molecules in bulk liquid water under the fixed-volume condition but it is much robust since the solvent density to be fixed can be anything between the coexisting vapor and liquid phases. A unique feature found in the water interface is that there is a minimum in the local solubility profile Σ(z) on the liquid side of the interface. We find that with decreasing temperature the minimum of Σ grows and at the same time the first peak in the oscillatory density profile of water develops. It is likely that the minimum of Σ is due to the layering structure of the free interface of water.

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

We report molecular simulations of oil-water liquid-liquid interfaces by using atomistic and coarse grained (CG) MARTINI force fields. We also apply the electronic continuum (EC) model to the MARTINI force field for the calculation of the interfacial tension of oil/water-salt systems. In a first step, we propose to calculate the interfacial tensions using thermodynamic and mechanical definitions of hydrocarbon-water interfacial systems modified by the addition of salts and alcohol. We also establish here the order of magnitude of the long-range corrections to the interfacial tension in fluid-fluid interfaces. Whereas the atomistic models are able to reproduce quantitatively the interfacial tension and the coexisting densities of oil-water systems, the coarse-description shows some deviations in the prediction of the interfacial tensions. Nevertheless, the physical features of these liquid-liquid interfaces are well-captured by this CG description. The CG force field offers then a very challenging alternative that will require however a more developed calibration of the parameters on the basis of liquid-liquid properties.

Recent advances in Many Body Dissipative Particles Dynamics simulations of liquid-vapor interfaces

Many Body Dissipative Particles Dynamics (MDPD) simulation is a novel promising mesoscopic method to model the liquid-vapor interfaces. Based upon works of Paganobarraga and Frenkel (J. Chem. Phys. 15, 5015 (2001)) and Trofimov (J. Chem. Phys. 117, 9383 (2002)) and of Warren (Phys. Rev. E 68, 066702 (2003)) this method has been critically reviewed during this last decade. We propose here to give an overview of the Many Body Dissipative Particles Dynamic simulation within the framework of the liquid-vapor interfaces. We recall the theoretical background of MDPD and we present some recent results of systems of interest such as water liquid-vapor interfaces and salt effect on water surface tension. Additionally we discuss the ability of MDPD to capture the mechanisms at the mesoscopic scale through the formation of micelles and the coalescence of a nanodroplet water on water surface.

Semi-infinite boundary conditions for the simulation of interfaces: the Ar/CO2(s) model revisited

We propose a method to account for the long tail corrections of dispersive forces in inhomogeneous systems. This method deals separately with the two interfaces that are usually present in a simulation setup, effectively establishing semi-infinite boundary conditions that are appropriate for the study of the interface between two infinite bulk phases. Using the wandering interface method, we calculate surface free energies of vapor-liquid, wall-liquid, and wall-vapor interfaces for a model of Lennard-Jones argon adsorbed on solid carbon dioxide. The results are employed as input to Young's equation, and the wetting temperature located. This estimate is compared with predictions from the method of effective interface potentials and good agreement is found. Our results show that truncating Ar-Ar interactions at two and a half molecular diameters results in a dramatic decrease of the wetting temperature of about 40%.

COMPUTER SIMULATION OF LIQUID–VAPOR INTERFACIAL TENSION: LENNARD-JONES FLUID AND WATER REVISITED

We review several commonly used simulation methods for computing liquid–vapor surface tension and associated theoretical treatments of the long-range correction for inhomogeneous systems. Prototype model systems considered in this review are the Lennard-Jones (LJ) fluid and the SPC/E model water. In addition, we examine a variety of factors that can affect calculation of the surface tension γ via the mechanical approach (i.e. using either KB or IK method). It is found that for the LJ fluid, the size of simulation box and the number of particles in the system can have notable effects on the computed surface tension. For SPC/E water, the Ewald parameters can influence computed surface tensions (γ) as well, e.g., very small Ewald parameters tend to overestimate γ. It is also found that the IK method consistently gives γ that are 0.6 - 0.9 mN/m greater than γ computed based on the KB method. When computing the first reciprocal–space contribution to the surface tension, the Ghoufi's strategy gives rise to more sensible profile of pressure difference PN(z)-PT(z) than the Alejandre's strategy although both strategies result in nearly the same average surface tension through the integration of PN(z)-PT(z).

Comparison of united-atom potentials for the simulation of vapor-liquid equilibria and interfacial properties of long-chain n-alkanes up to n-C100

Canonical ensemble molecular dynamics (MD) simulations are reported which compute both the vapor-liquid equilibrium properties (vapor pressure and liquid and vapor densities) and the interfacial properties (density profiles, interfacial tensions, entropy and enthalpy of surface formation) of four long-chained n-alkanes: n-decane (n-C(10)), n-eicosane (n-C(20)), n-hexacontane (n-C(60)), and n-decacontane (n-C(100)). Three of the most commonly employed united-atom (UA) force fields for alkanes (SKS: Smit, B.; Karaborni, S.; Siepmann, J. I. J. Chem. Phys. 1995,102, 2126-2140; J. Chem. Phys. 1998,109, 352; NERD: Nath, S. K.; Escobedo, F. A.; de Pablo, J. J. J. Chem. Phys. 1998, 108, 9905-9911; and TraPPE: Martin M. G.; Siepmann, J. I. J. Phys. Chem. B1998, 102, 2569-2577.) are critically appraised. The computed results have been compared to the available experimental data and those fitted using the square gradient theory (SGT). In the latter approach, the Lennard-Jones chain equation of state (EoS), appropriately parametrized for long hydrocarbons, is used to model the homogeneous bulk phase Helmholtz energy. The MD results for phase equilibria of n-decane and n-eicosane exhibit sensible agreement both to the experimental data and EoS correlation for all potentials tested, with the TraPPE potential showing the lowest deviations. However, as the molecular chain increases to n-hexacontane and n-decacontane, the reliability of the UA potentials decreases, showing notorious subpredictions of both saturated liquid density and vapor pressure. Based on the recommended data and EoS results for the heaviest hydrocarbons, it is possible to attest, that in this extreme, the TraPPE potential shows the lowest liquid density deviations. The low absolute values of the vapor pressure preclude the discrimination among the three UA potentials studied. On the other hand, interfacial properties are very sensitive to the type of UA potential thus allowing a differentiation of the potentials. Comparing the interfacial tension MD results to the available experimental data and SGT results, the TraPPE model exhibits the lowest deviations for all hydrocarbons.