Surface tension of hydrocarbon chains at the liquid–vapour interface

Molecular dynamics simulations for interfacial structure and affinity between carboxylic acid-modified Al2O3 and polymer melts

Controlling the dispersion state of nanoparticles in a polymer matrix is necessary to produce polymer nanocomposites. The surface modification of nanoparticles is used to enable their dispersion in polymers. Moreover, molecular dynamics (MD) simulations are useful for revealing the interfacial properties between nanoparticles and polymers to aid in the design of materials. In this study, the effect of surface coverage, modifier length, and polymer species on the interfacial structure and affinity between surface-modified Al2O3 and polymer melts were investigated using all-atom MD simulations. Hexanoic, decanoic, and tetradecanoic acids were used as surface modifiers, and polypropylene (PP), polystyrene (PS), and poly (methyl methacrylate) (PMMA) were used as polymers. The work of adhesion Wadh and the work of immersion Wimm were selected as quantitative measures of affinity. Wadh was calculated using the phantom-wall approach, and Wimm was calculated by simply subtracting the surface tension of polymers γL from Wadh. The results showed that Wadh and Wimm were improved by surface modification with low coverage, owing to a good penetration of the polymer. The effect of modifier length on Wadh and Wimm was small. Whereas Wadh increased in the following order: PP < PS < PMMA, Wimm increased as follows: PMMA < PS < PP. Finally, the trend of Wadh and Wimm was organized using the Flory-Huggins interaction parameter χ between the modifier and the polymer. This study demonstrates that the interfacial affinity can be improved by tuning the surface coverage and modifier species depending on the polymer matrix.

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.

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.

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

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

Structure and Interfacial Tension of a Hard-Rod Fluid in Planar Confinement

The structural properties and interfacial tension of a fluid of rodlike hard-spherocylinder particles in contact with hard structureless flat walls are studied by means of Monte Carlo simulation. The calculated surface tension between the rod fluid and the substrate is characterized by a nonmonotonic trend as a function of the bulk concentration (density) over the range of isotropic bulk concentrations. As suggested by earlier theoretical studies, a surface-ordering scenario is confirmed by our simulations: the local orientational order close to the wall changes from uniaxial to biaxial nematic when the bulk concentration reaches about 85% of the value at the onset of the isotropic-nematic phase transition. The surface ordering coincides with a wetting transition whereby the hard wall is wetted by a nematic film. Accurate values of the fluid-solid surface tension, the adsorption, and the average particle-wall contact distance are reported (over a broad range of densities into the dense nematic region for the first time), which can serve as a useful benchmark for future theoretical and experimental studies on confined rod fluids. The simulation data are supplemented with predictions from second-virial density functional theory, which are in good qualitative agreement with the simulation results.

Molecular dynamics simulations of 2-(dimethylamino)ethanol (DMEA)

We develop a multipurpose force field to investigate the properties of the condensed phases of 2-(dimethylamino)ethanol (DMEA). We use ab initio computations at the HF/6-311++G(2d,2p) level to derive partial charges, obtain force constants, and compute the electrostatic potential of the DMEA. We find that the HF predictions for the dipole moment are in excellent agreement with the experimental result (2.6 D). The computations also show the strong preference of DMEA to form intramolecular hydrogen bonds between the hydrogen in the alcohol group and nitrogen. We have tested the accuracy of our force field by computing coexistence and interfacial properties as well as thermal conductivities in wide range of thermodynamic states. In all these instances we find excellent agreement with the available experimental data. We have further investigated the structure of the liquid by computing pair correlations. Our results indicate a clear preference for DMEA to form low-dimensional structures, such as linear and bifurcated chains, which are driven by the association of the DMEA molecules via the alcohol group. Overall, our force field provides a good basis to compute the bulk and interfacial properties of DMEA.

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.