Abstract
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.
Similar content being viewed by others
Introduction
The liquid-liquid interface plays an important role in many physical, chemical and biological processes. From phase transfer catalysis to liquid chromatography; liquid-liquid extraction to pharmaceutical drug delivery1. In particular, the interface between water and oil is becoming increasingly examined within the scientific community. Recent work by Mugele et al. introduces an intriguing unsolved problem: that the wettability of a surface can be manipulated by altering the ionic composition of the surrounding fluid2. Mugele et al. examined the three-phase wettability of a mineral-oil-water system, measuring contact angle variations between mica, water and alkane as a function of salt concentration and composition. Despite extensive scientific efforts, the underlying mechanisms driving this wettability change are yet to be fully discerned, especially at the molecular level3. Mugele et al. proposed several models explaining the behaviour, and further challenged the computational chemistry community to help explain the phenomenon.
Classical MD simulations of salinity induced wettability alterations are becoming more prevalent in literature. Work by Jiménez-Ángeles & Firoozabadi4 attempt to solve the problem posed by Mugele directly. The study systematically measures the contact angle formed at the three-phase interface at varying salt concentrations. The dynamic properties of wettability as a function of salinity have further been examined by Zhang et al.5, who examine the non-equilibrium flow of oil through brines in various nanopores. No single force field (the input parameters for the MD simulation) is truly universal, and most MD studies of salinity induced wettability alterations combine at least three separate force fields without prior validation of their respective interfacial properties. Indeed, in the authors previous work6,7, the CGenFF8,9,10 model of oil was combined with the TIP3Pc8 water model and the ClayFF11 force field to model the mineral surface without thorough appreciation of the interfacial properties in the system. Without prior validation of the interfacial properties at the mineral-oil, mineral-water and water-oil interface, the subsequent results of such MD simulations can be questioned. Consequently, it is becoming increasingly important to understand the properties of the noted interfaces and to validate their respective properties against experimental data. In the present study, we focus on the water-oil interface, as it possesses a variety of interesting features, and has a large set of experimental data to compare against. Using classical MD simulations, we are able to examine the properties of the oil-water interface with atomic resolution, and how these properties change upon introduction (or depletion) of inorganic salts.
Experimental studies of the atomic structuring of fluid-fluid interfaces have historically been hindered by the lack of appropriate techniques capable of probing the interface with sufficient resolution. Recently, however, surface sensitive experimental techniques, such as sum frequency generation (SFG) spectroscopy, second harmonic generation (SHG) spectroscopy, and x-ray reflectivity measurements, have been able to probe the properties of fluid-vapour and fluid-fluid interfaces12. For example, Braslau et al. measured the surface width of the water-vapour interface (the distance between bulk-like liquid region of water and bulk-like vapour region) as 3.2 Å using x-ray reflectivity and capillary wave theory13, whilst similar experiments by Sanyal et al. discerned the ethanol-vapour interfacial width as 6.9 Å14. The interface between water and alkane has also been experimentally measured using x-ray reflectivity, as in the study of Mitrinovic et al., who measured the width of the water-hexane interface as 3.5 Å15 and Tikhonov et al., who measured the interfacial width of water-docosane (C 22 H 46) as 5.7 Å16.
Such advances in experimental techniques have been mirrored in the computational realm. Classical molecular dynamics (MD) simulations have been used to help interpret the phenomena present at the interfaces of water-vapour17, electrolyte-vapour18 and water-alkane19,20. The interfacial properties of water-vapour were calculated by Matsumoto & Kataoka in the late 1980s17, whilst the role of salts at the water-vapour interface was studied critically by Jungwirth & Tobias18. Water-alkane interfaces have also been heavily studied, for example by van Buuren et al., whom looked at the interface between water and decane19, and by Rivera et al., who looked at the water-alkane interface with the addition of methanol20. More recently, computation techniques have been able to resolve the intrinsic density profile of the liquid-vapour and liquid-liquid interface. That is, the interface between two phases excluding the contribution of thermal fluctuations (capillary waves), which act as to smear the density profiles across the interface. For example, Partay et al. were able to examine the intrinsic interface of the water-vapour system21, whilst Hantal et al. examined the interface between water and CCl4 22. Using such techniques improves the clarity of the molecular structuring at the interface, and is seeing increasing use when studying the water-alkane interface23,24. The examination of the intrinsic density profile frequently leads to more insightful observations from the simulations.
As previously stated, no single force field is truly universal. Many different water models exist in the literature, and most excel at capturing certain aspects of water’s unusual physical properties. For example, the TIP4P2005 model excels at modelling the phase behaviour of water over a variety of temperatures and pressures, especially when compared to other models such as SPC/E25,26. However, the SPC/E model of water more accurately reproduces the experimental dielectric constant of water25,26. As no one water model is truly universal, it is not uncommon to observe various models being used in three-phase MD simulations. Commonly used water models applied to low-salinity EOR include the SPC27, SPC/E28, TIP3P29 and TIP4P29 parameterisations. Similar arguments also apply to the various parameterisation of oil molecules. Typical oil models used in previous work include CGenFF8, TraPPE30 and OPLS31. Both TraPPE and OPLS were parameterised to model organic solvents (with TraPPE specifically parameterised to model alkanes)30,31. OPLS, however, is frequently combined with ClayFF5,32 (a force field frequently used to model mineral surfaces), despite the fact that both force fields are parameterised with different Lennard-Jones mixing rules, and are therefore not inherently compatible. In contrast to TraPPE and OPLS, CGenFF is more general. CGenFF offers fully automated atom typing, and therefore offers large throughput.
In the present work, we investigate the utility of various force fields to accurately reproduce the interfacial behaviour between water and alkane at different NaCl concentrations. We have used decane as a representative model alkane. The aim of the article is to interpret the phenomena of the liquid-liquid interface at the molecular level, and to trial various force fields for future use in more complex three-phase simulations. Classical MD simulations have been used to calculate the surface tension of water and decane in their liquid-vapour phases. The interfacial tension of the water-decane interface has also been examined. The properties of six different force fields for water (SPC, SPC/E, TIP3P, TIP3Pc, TIP4P and TIP4P2005), and three further force fields for alkane (CGenFF, OPLS-AA and TraPPE-UA) are compared against experimental data. The salinity dependence on surface/interfacial tension has been modelled using the Smith & Dang parameterisation of NaCl. Both the intrinsic and non-intrinsic density profiles for each interface has been calculated to show the structural similarities and differences between each force field.
Results
The Alkane-Vapour Interface
Figure 1 presents the calculated surface tension values of decane over a range of temperatures. The results show that the three tested force fields accurately reproduce the overall trend of surface tension at various temperatures.
Whilst we have examined decane solely, the considered force fields are widely transferable. Both OPLS-AA and TraPPE-UA have previously been shown to accurately predict the surface tension values of ethane through to hexadecane33,34. To the author’s knowledge, there have been no surface tension calculations of alkanes using the CGenFF force field.
Both OPLS-AA and TraPPE-UA systematically underestimate the surface tension of decane at all temperatures, OPLS more-so than TraPPE-UA. CGenFF overestimates surface tension at lower temperatures, and underestimates surface tension at higher temperatures. Our results concur with the work of Ismail et al., who show that OPLS-AA consistently underestimates the surface tension of various saturated n-alkanes33 and with Mendoza et al., who present that TraPPE-UA is generally in good agreement with experimental data34. Our results are also consistent with the study of Caleman et al., who show that OPLS-AA often underestimates the surface tension values of a wide variety of organic liquids35.
Previous work by Nicolas & Smit has shown that the accurate calculation of surface tension for saturated alkanes is linearly proportional to the liquid density predicted in the simulations36. Values for the liquid and vapour density of decane have been calculated by fitting the interfacial density profile to a tanh function (see Equation 6 in Methods). The liquid density (ρ l ) of decane at 293.15 K is presented in Table 1, along with values for the interfacial width (Δ). Note that the vapour density (ρ v ) at 293.15 K is calculated to be zero, and is therefore not presented.
Our results show that both TraPPE-UA and OPLS-AA marginally underestimate the liquid density of decane at 293.15 K, and consequently, both underestimate the surface tension. Conversely, CGenFF overestimates the liquid density, and therefore the surface tension of decane. The accurate reproduction of liquid-vapour densities depends on the treatment of long-range dispersion forces37. Ismail et al. show that Lennard-Jones cutoffs greater than 16 Å are required to yield results in agreement with explicit long range methods (evaluated in reciprocal space)33, whilst the simulations of Mendoza et al. and Lopez et al. also show the utility of using Ewald-summations to treat these long range dispersion forces34,38. Note that our simulations satisfy both these criteria.
Figure 2 presents the interfacial density profiles (both intrinsic and non-intrinsic) of the decane-vapour interface at 293.15 K. The density profiles show that all three force fields act remarkably similarly, even though CGenFF and OPLS-AA are all-atom representations of decane, and TraPPE-UA is a united-atom representation. The intrinsic density profiles present four discrete peaks of decane before the density becomes bulk-like.
The Water-Vapour Interface
The calculated surface tension values of pure water are presented in Table 2. It can be seen that almost all water models underestimate the experimentally observed surface tension of water by at least 15%. The only exception to this is TIP4P2005, which underestimates the surface tension of water by 7%. Our results compare favourably to the study of Vega & De Miguel25. Like their work, we find that the accurate reproduction of surface tension follows the trend:
TIP4P2005 > SPC/E > TIP4P > SPC > TIP3P
The TIP3Pc model, which introduces LJ sites on hydrogen atoms, performs better than the original TIP3P parameterisation, but still underestimates the surface tension of water by 23%. We observe that, generally, the four-point water models perform better than the three-point models, in keeping with the observations of Vega & De Miguel25. Of all the three-point water models, SPC/E most accurately reproduces the experimental surface tension of water, whilst of the four-point water models, the TIP4P2005 parameterisation performs the best.
The density profiles of TIP4P2005 and SPC/E water are presented in Fig. 3. Analysing the intrinsic density profiles, it can be seen that water forms three separate density peaks at the water-vapour interface. This information is lost due to the capillary wave fluctuations present in the simulation when examining the non-intrinsic density profile solely. Also notable, is that most water models appear to overestimate the interfacial width of water compared to the results of Braslau et al.13, see Table 2. The TIP4P2005 and SPC/E water model most accurately predict the experimentally observed interfacial width of water. Overall, the TIP4P2005 is best placed to model water-vapour interfaces. Of the three-point water models, SPC/E is the most favourable.
The Water-Vapour Interface at Various NaCl Concentrations
Proceeding, we have chosen three water models to examine the trends in surface tension as a function of salinity. The TIP4P2005 and SPC/E water models were examined as they represent the best four-point and three-point water models respectively. The TIP3Pc water model was also examined, as it is frequently simulated in conjunction with secondary organic phases6,7. The variation of surface tension with NaCl concentration is presented in Fig. 4.
Firstly, we note that the overall trend of surface tension variation with NaCl concentration is well captured by the classical MD simulations. The simulations are able to accurately predict the experimental surface tensions of NaCl solutions up to 1.5 mol/kg, however, the simulations slightly overestimate the surface tension of NaCl electrolyte at higher concentrations, deviating more drastically as concentrations increase beyond 2.0 mol/kg. Similar behaviour was noted by D’Auria & Tobias, who observed that the surface tension of the TIP3P water model was overestimated when used in conjunction with the Smith & Dang parameterisation at NaCl concentrations of 1.2 mol/kg and 6.17 mol/kg39.
Secondly, we observe that the simulated surface tensions of all three water models agree with each other at each concentration. This implies that the accurate modelling of interfaces containing electrolytes depends more so on the parameterisation of the ion, rather than the water model. This conclusion concurs with the work of Neyt et al., who presented that the salinity dependence of the electrolyte-vapour interface varied primarily with ion parameterisation, rather than the water model, whilst the utilised water model sets the accuracy of the surface tension calculation without additional NaCl40. Notably, Neyt et al. presented that the TIP4P2005 water model used in conjunction with the OPLS parameterisation of NaCl, typically returned surface tensions within 0.2% of the experimental trend40.
The Water-Alkane Interface at Various NaCl Concentrations
Table 3 presents the interfacial tension values of the decane-water interface, without salts, for various force field combinations as predicted by the MD simulations. We have chiefly examined TraPPE-UA and CGenFF due to their compatibility (i.e. consistent use of Lennard-Jones mixing rules) with ClayFF; which makes these force fields more suitable for use in complex three-phase simulations. Priority has also been placed on testing the TIP4P2005 and SPC/E water models, as these most accurately model the interfacial properties of pure water. TIP3Pc has also been examined as it is the default three-point water model used in conjunction with CGenFF8.
In addition to the usual TIP4P2005 water model interacting with TraPPE-UA, we have examined the parameters as described by Ashbaugh et al.41, denoted TIP4P2005*. The TIP4P2005* water model uses custom Lennard-Jones interactions between the oxygen atoms of water and the TraPPE-UA beads. These custom interactions have been finely tuned to accurately reproduce the hydration free energy of alkanes41. All other parameters of the TIP4P2005* water model are identical to the original TIP4P2005 parameterisation.
We note that the MD simulations are able to predict the interfacial tension of the water-decane interface within 10% of the empirically observed values. The simulations of this interface are notably more accurate than the simulations of the water-vapour interface. The combination of different force fields appears more robust for the water-decane interface, compared to the water-vapour interface. This is likely due to the absence of a vapour phase in the simulations, which when inaccurately modelled (in particular, the vapour density), can lead to deviations in the predicted surface tension, as noted by Nicolas & Smit36.
Figure 5 presents the deviation of decane-water interfacial tension with increasing NaCl concentration for various combinations of force fields. The computed interfacial tensions accurately follow the experimental trend observed by Aveyard & Saleem42 up to 1.0 mol/L. Much like the water-vapour interface, the different parameter sets agree with each other at each NaCl concentration. This again highlights that the ion parameterisation may be more important for the accurate modelling of the water-vapour/water-decane interface.
Discussion
The presented results show the quality with which MD simulations can capture the bulk properties and structural properties of the liquid-liquid interface. One key advantage of such simulations, compared to many experiments, is the unparalleled resolution with which one can interpret the bulk phenomena on display. In this instance, we can use the MD simulations to further our understanding of why the surface tension at the water-vapour/water-alkane interface increases with increasing NaCl concentration. Typically, this behaviour is explained due to an exclusion of ions from the liquid-vapour interface3. This directly alters the surface tension between the two phases as described by the Gibbs adsorption isotherm:
where Γ i is the surface excess and μ i is the chemical potential of species i∈{ Na, Cl }. In this study, the surface excess of Na+ and Cl− ions across the electrolyte-water interface has been calculated as:
where ρ i (z) is the non-intrinsic density profile of phase i across the interface, \({\rho }_{i}^{{\rm{liq}}}\) is the bulk liquid density of phase i, and \({\rho }_{i}^{{\rm{vap}}}\) is the bulk vapour density of phase i (typically \({\rho }_{i}^{{\rm{vap}}}=0\) for ions in liquid-vapour simulations). Z Gibbs is the location of the Gibbs dividing plane, which is priorly calculated by minimising Γwater, such that:
Note that equation 3 has been used to centre the non-intrinsic density profiles presented in the results section. Figure 6 presents the negative total surface excess of ions (\({{\rm{\Gamma }}}_{{\rm{NaCl}}}={{\rm{\Gamma }}}_{{\rm{Na}}}+{{\rm{\Gamma }}}_{{\rm{Cl}}}\)) at the water-vapour interface calculated from all simulations as a function of NaCl concentration. We observe that, compared to the experimental data of Ali et al., the simulations marginally overestimate the negative total surface excess of NaCl, \(-{{\rm{\Gamma }}}_{{\rm{NaCl}}}\), especially at higher concentrations43. In turn, this elucidates why the surface tension values are overestimated in the simulations at higher salt concentrations. The surface tension values are overestimated due to the overestimation of ionic exclusion at the liquid-vapour interface. Notably, the lack of polarization may explain why our classical molecular dynamics simulations overestimate the ionic exclusion. It is well known that, whilst ion exclusion generally increases with salt concentration, certain anions (for example, bromide and iodide) can accumulate at the interface due to the polarization44. Recent computational work by Neyt et al. and Neyt et al. has examined the effects of polarization on electrolyte-vapour interfaces40,44. Both studies concluded that current models of polarizable force fields (including both polarizable water and/or ions), are not mature enough to capture the salinity dependence on surface tension. These studies tested both the Drude oscillator model and the electronic continuum correction model (whereby classical point charges are shifted by a screening factor to account for polarization). Work by Jiang & Panagiotopolous further show that theelectronic continuum correction model fails to accurately predict the interfacial properties of electrolytes with polarizable models45. More recently, work by Jiang et al. presents that the BK3 polarizable model is able to accurately capture the interfacial properties of NaCl electrolyte solutions46. However, it is worth noting that such polarizable models typically run 5–10 times slower compared to the classical parameterizations used in this study, and that the implementation of the BK3 model currently only runs on a single processor using GROMACS (the software suite used in this study).
Also noteworthy is that the variation in interfacial tension predicted at the interface between water and decane is remarkably similar to that at the water-vapour interface. This is apparent in Fig. 7, which compares the change in surface tension at the water-vapour interface against the change in interfacial tension at the water-alkane interface (using the TIP4P2005 water model and the TraPPE-UA oil model). The overall change in interfacial tension is therefore due to the properties primarily within the water phase. This result is somewhat unsurprising as the dielectric properties of decane and water vapour are much similar when compared to bulk liquid water. The surface tensions and interfacial tensions scale linearly as the ions are excluded from their respected interfaces in equal value, and the ions are not miscible in the vapour/decane phase respectively.
To conclude; in this study, the interfacial properties of water and decane have been examined at various NaCl concentrations using classical molecular dynamics simulations. By choosing an appropriate set of interaction parameters (force fields), one can obtain remarkable agreement between model and experimental observation. In particular, the TIP4P2005 water model is best placed to examine the interfacial properties of water. The SPC/E parameterisation is the best three-point water model to interpret the interfacial behaviour of water. CGenFF, OPLS-AA and TraPPE-UA all accurately reproduce the interfacial properties of decane. In combination, CGenFF and TraPPE-UA are compatible with various water models, and are able to accurately predict the interfacial tension of the water-decane interface. The salinity dependence on surface/interfacial tension is well captured using the Smith & Dang parameterisation of NaCl. We observe that the model slightly overestimates the surface/interfacial tensions at higher salinities. This is due to an overestimation of the ion exclusion at the interface. The simulations further suggest that the salinity dependence on surface/interfacial tension is dictated by the parameterisation of the salt ions. Future work will examine this hypothesis, and will model different ion parameterizations. Future work will also examine the role of divalent cations at the water-vapour/water-alkane interface, whereby, charge inversion may play a determining role on the behaviour of the liquid-liquid interface.
Methods
The datasets generated and analysed during the current study are available from the corresponding authors on reasonable request.
All simulations were calculated using GROMACS 5.1.447 with an electrostatic and van der Waals cutoff radii of 2.0 nm. Long range electrostatics were calculated using a Particle-Mesh-Ewald (PME) summation with grid spacings of 0.1 nm. The PME summation used a spline interpolation order of 4, and long-range electrostatic interactions were accurate to within 99.999%.
All simulations were initialised with an energy minimisation calculation to minimise any unphysical atomic overlaps. This was achieved using a steepest descents algorithm, which was terminated once the maximum force on any one atom was less than 100 kJ/(mol nm).
All simulations were subsequently equilibrated for 1 ns. Simulations of liquid-vapour phase surface tensions were calculated in the canonical ensemble (constant particle number - N, constant volume - V and constant temperature T) at 293.15 K, using a V-rescale thermostat set to rescale system temperatures every 0.5 ps. Liquid-liquid interfacial tension simulations were calculated in the isothermal-isobaric ensemble (constant particle number - N, constant volume - V and constant pressure - P) at 293.15 K and 1.01325 bar. NPT simulations were equilibrated using a V-rescale thermostat with a temperature coupling constant of 0.5 ps. Pressure coupling was achieved using a Berendsen barostat, with a pressure coupling constant of 1 ps.
Following equilibration, all simulations were run for a production period of 10 ns. During the production period, all simulations used a Nosé-Hoover thermostat with a temperature coupling constant of 0.5 ps and a Nosé-Hoover chain length of 1. NPT simulations used an isotropic Parrinello-Rahman barostat during the production period, with a pressure coupling constant of 1 ps. Unless otherwise stated, all production simulations have been calculated at 293.15 K and 1.01325 bar.
Three different force field parameterisations for decane have been tested in this study: TraPPE-UA30, OPLS-AA31,48 and CGenFF8,9,10. TraPPE-UA is a united atom force field, where CH3- and -CH2- groups are modelled as soft spheres. OPLS-AA and CGenFF are both all atom models, where carbon and hydrogen atoms are modelled explicitly. Six different water models were examined in the water-vapour interface: SPC27, SPC/E28, TIP3P29, TIP3Pc8, TIP4P29 and TIP4P200549. TIP3Pc refers to the CHARMM variant of the original TIP3P parameterisation, whereby hydrogen atoms contain Lennard-Jones sites. The Smith and Dang parameterisation has been used to model aqueous sodium and chloride ions50.
System Setup
Three different systems are presented in this study. In the first, a 5 × 5 × 5 nm3 box of 392 decane molecules is inserted into 5 × 5 × 20 nm3 simulation box. In the second, a 5× 5× 5 nm3 box of 4139 water molecules is inserted into 5 × 5× 20 nm3 simulation box. The number of decane and water molecules in each film is calculated to match the bulk density of each solvent at 293.15 K and 1.01325 bar, 0.727 g/cm3 for decane and 0.9982 g/cm3 for water respectively51. In the third system, the decane and water films are combined in a 5 × 5× 10 nm3 simulation box. Each system was generated using the Packmol software package52. Systems involving water were further examined at various NaCl concentrations. Systems were setup in terms of NaCl molal concentration, up to a maximum of 3.00 mol/kg, in increments of 0.50 mol/kg. This was achieved by replacing water molecules with the relevant amount of sodium and chloride ions. The number of water molecules and ions present in each simulation is presented in Table 4. Additionally, 0.20 mol/kg NaCl solution was examined. Whilst simulation results are presented in molal concentration (mol/kg), experimental results are often presented in terms of molar concentration (mol/L). The conversion from molar concentration to molal concentration is presented in Table 4, using data extrapolated from the CRC Handbook of Chemistry & Physics53.
Analyses
Thermodynamic data from each simulation were output every 1 ps. Final values for thermodynamic quantities were averaged over all 10 ns, and errors were calculated using a block-averaging method, with each block averaging over a 1 ns timeframe. In all figures, error bars are presented to ±2 standard errors of the mean (a confidence interval of 95%).
The surface tension across an interface has been calculated using the diagonal components of the local pressure tensor:
where L z is the length of the simulation in z (the direction normal to the interface) and p N(z) & p T(z) represent the normal and tangential components of the pressure tensor with respect to the interface:
The diagonal components of the local pressure tensor (p xx, p yy & p zz) have been calculated using the Irving-Kirkwood formalism54.
Density profiles across the interface have been calculated following a two-stage process. Firstly, the simulation trajectory is centred about the centre of mass of the primary solvent phase in each simulation (typically water). This reduces artefacts caused by the collective drift of the interface throughout the simulation. The primary phase is calculated by clustering all molecules in the system. The largest cluster is then selected as the primary (liquid) phase. Any molecule further than 0.35 nm from the primary phase is excluded from the centre of mass calculation, and therefore does not affect the centering of the system. A cutoff of 0.35 nm was selected as this corresponds to the first minimum in the radial distribution function (RDF) of water55,56. Consequently, water molecules in the vapour phase are excluded in the centre of mass calculation for the bulk liquid water phase, and therefore do not artificially shift the resulting density profiles normal to the interface. After the system has been centred, the density profile is calculated across the interface using a bin size of 0.01 nm. Where applicable, the density profile has been fit to the equation:
where \({\rho }_{l}\) is the liquid density of the primary phase, \({\rho }_{v}\) is the vapour density of the primary phase, z0 is the location of the interface, and Δ is the interfacial width. The density profiles of Na+ and Cl− ions have subsequently been calculated relative to the definition of the water interface.
The density profiles calculated using the above methodology are subject to capillary waves due to thermal fluctuations. Recently, computation techniques have been able to resolve the intrinsic density profile of the liquid-vapour and liquid-liquid interface21. That is, the interface between two phases excluding the contribution of thermal fluctuations (capillary waves), which act as to smear the density profiles across the interface. The amplitude of these capillary waves scales as:
where the maximum wave vector, q, depends upon the size of the simulated interface.
The calculation of the intrinsic density profile has been evaluated by offsetting the amplitude of the thermal fluctuations (ξ) from the interface:
where index i sums over all atoms of phase i, and z is the position of the local non-intrinsic interface. The intrinsic density profiles have been calculated using the ITIM method21 as presented by Sega et al.57, using a probe sphere radius of 0.2 Å. Within the ITIM method a probe sphere is moved along test lines perpendicular to the plane of the fluid-vapour or fluid-fluid interface. Once the probe sphere touches the first atom within of the phase of interest, this molecule is marked as being interfacial. This process is repeated over the entire interfacial area in the simulation. The intrinsic density profile is then calculated using the offset (ξ) of the marked interfacial atoms.
References
Benjamin, I. Molecular structure and dynamics at liquid-liquid interfaces. Annual Review of Physical Chemistry 48, 407–451 (1997).
Mugele, F. et al. Ion adsorption-induced wetting transition in oil-water-mineral systems. Scientific reports 5 (2015).
Schlossman, M. Liquid-liquid interfaces: studied by X-ray and neutron scattering. Current Opinion in Colloid & Interface Science 7, 235–243 (2002).
Jiménez-Ángeles, F. & Firoozabadi, A. Tunable substrate wettability by thin water layer. The Journal of Physical Chemistry C 120, 24688–24696 (2016).
Zhang, L., Lu, X., Liu, X., Yang, K. & Zhou, H. Surface wettability of basal surfaces of clay minerals: insights from molecular dynamics simulation. Energy & Fuels 30, 149–160 (2016).
Underwood, T., Erastova, V., Cubillas, P. & Greenwell, H. Molecular dynamic simulations of montmorillonite-organic interactions under varying salinity: an insight into enhanced oil recovery. Journal of Physical Chemistry C 119, 7282–7294 (2015).
Underwood, T., Erastova, V. & Greenwell, H. Wetting effects and molecular adsorption at hydrated kaolinite clay mineral surfaces. Journal of Physical Chemistry C 120, 11433–11449 (2016).
Vanommeslaeghe, K. et al. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. Journal of Computational Chemistry 31, 671–690 (2010).
Vanommeslaeghe, K. & MacKerell, A. Jr. Automation of the CHARMM general force field (CGenFF) i: bond perception and atom typing. Journal of Chemical Information and Modeling 52, 3144 (2012).
Vanommeslaeghe, K., Raman, E. & MacKerell, A. Jr. Automation of the CHARMM general force field (CGenFF) ii: assignment of bonded parameters and partial atomic charges. Journal of Chemical Information and Modeling 52, 3155 (2012).
Cygan, R., Liang, J.-J. & Kalinichev, A. Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field. Journal of Physical Chemistry B 108, 1255–1266 (2004).
Eisenthal, K. Liquid interfaces probed by second-harmonic and sum-frequency spectroscopy. Chemical Reviews 96, 1343–1360 (1996).
Braslau, A. et al. Surface roughness of water measured by x-ray reflectivity. Physical Review Letters 54, 114 (1985).
Sanyal, M., Sinha, S., Huang, K. & Ocko, B. X-ray-scattering study of capillary-wave fluctuations at a liquid surface. Physical Review Letters 66, 628 (1991).
Mitrinovic, D., Zhang, Z., Williams, S., Huang, Z. & Schlossman, M. X-ray reflectivity study of the water-hexane interface. Journal of Physical Chemistry B 103, 1779–1782 (1999).
Tikhonov, A., Mitrinovic, D., Li, M., Huang, Z. & Schlossman, M. An x-ray reflectivity study of the water-docosane interface. Journal of Physical Chemistry B 104, 6336–6339 (2000).
Matsumoto, M. & Kataoka, Y. Study on liquid-vapor interface of water. i. simulational results of thermodynamic properties and orientational structure. Journal of Chemical Physics 88, 3233–3245 (1988).
Jungwirth, P. & Tobias, D. Specific ion effects at the air/water interface. Chemical Reviews 106, 1259–1281 (2006).
van Buuren, A., Marrink, S. & Berendsen, H. A molecular dynamics study of the decane/water interface. Journal of Physical Chemistry 97, 9206–9212 (1993).
Rivera, J., McCabe, C. & Cummings, P. Molecular simulations of liquid-liquid interfacial properties: Water-n-alkane and water-methanol-n-alkane systems. Physical Review E 67, 011603 (2003).
Pártay, L., Hantal, G., Jedlovszky, P., Vincze, Á. & Horvai, G. A new method for determining the interfacial molecules and characterizing the surface roughness in computer simulations. application to the liquid-vapor interface of water. Journal of Computational Chemistry 29, 945–956 (2008).
Hantal, G., Darvas, M., Pártay, L., Horvai, G. & Jedlovszky, P. Molecular level properties of the free water surface and different organic liquid/water interfaces, as seen from ITIM analysis of computer simulation results. Journal of Physics: Condensed Matter 22, 284112 (2010).
Bresme, F., Chacó’n, E. & Tarazona, P. Molecular dynamics investigation of the intrinsic structure of water-fluid interfaces via the intrinsic sampling method. Physical Chemistry Chemical Physics 10, 4704–4715 (2008).
Chen, M. et al. Retardation of water reorientation at the oil/water interface. Journal of Physical Chemistry C 119, 16639–16648 (2015).
Vega, C. & De Miguel, E. Surface tension of the most popular models of water by using the test-area simulation method. Journal of Chemical Physics 126, 154707 (2007).
Kusalik, P. & Svishchev, I. The spatial structure in liquid water. Science 265, 1219–1221 (1994).
Berendsen, H., Postma, P., van Gunsteren, W. & Hermans, J. Interaction models for water in relation to protein hydration. In Intermolecular Forces, 331–342 (Springer, 1981).
Berendsen, H., Grigera, J. & Straatsma, T. The missing term in effective pair potentials. Journal of Physical Chemistry 91, 6269–6271 (1987).
Jorgensen, W., Chandrasekhar, J., Madura, J., Impey, R. & Klein, M. Comparison of simple potential functions for simulating liquid water. Journal of Chemical Physics 79, 926–935 (1983).
Martin, M. & Siepmann, I. Transferable potentials for phase equilibria. 1. united-atom description of n-alkanes. Journal of Physical Chemistry B 102, 2569–2577 (1998).
Jorgensen, W., Maxwell, D. & Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. Journal of the American Chemical Society 118, 11225–11236 (1996).
Greathouse, J. et al. Methylene blue adsorption on the basal surfaces of kaolinite: Structure and thermodynamics from quantum and classical molecular simulation. Clays and Clay Minerals 63, 185–198 (2015).
Ismail, A., Tsige, M., Veld, P. & Grest, G. Surface tension of normal and branched alkanes. Molecular Physics 105, 3155–3163 (2007).
Mendoza, F., Lopez-Rendon, R., Lopez-Lemus, J., Cruz, J. & Alejandre, J. Surface tension of hydrocarbon chains at the liquid-vapour interface. Molecular Physics 106, 1055–1059 (2008).
Caleman, C. et al. Force field benchmark of organic liquids: density, enthalpy of vaporization, heat capacities, surface tension, isothermal compressibility, volumetric expansion coefficient, and dielectric constant. Journal of Chemical Theory and Computation 8, 61–74 (2011).
Nicolas, J. & Smit, B. Molecular dynamics simulations of the surface tension of n-hexane, n-decane and n-hexadecane. Molecular Physics 100, 2471–2475 (2002).
Sega, M. & Dellago, C. Long-range dispersion effects on the water/vapor interface simulated using the most common models. Journal of Physical Chemistry B 121, 3798–3803 (2017).
López-Lemus, J., Romero-Bastida, M., Darden, T. & Alejandre, J. Liquid-vapour equilibrium of n-alkanes using interface simulations. Molecular Physics 104, 2413–2421 (2006).
D’Auria, R. & Tobias, D. Relation between surface tension and ion adsorption at the air-water interface: A molecular dynamics simulation study. Journal of Physical Chemistry A 113, 7286–7293 (2009).
Neyt, J., Wender, A., Lachet, V., Ghoufi, A. & Malfreyt, P. Prediction of the concentration dependence of the surface tension and density of salt solutions: atomistic simulations using drude oscillator polarizable and nonpolarizable models. Physical Chemistry Chemical Physics 15, 11679–11690 (2013).
Ashbaugh, H., Liu, L. & Surampudi, L. Optimization of linear and branched alkane interactions with water to simulate hydrophobic hydration. Journal of Chemical Physics 135, 054510 (2011).
Aveyard, R. & Saleem, S. Interfacial tensions at alkane-aqueous electrolyte interfaces. Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases 72, 1609–1617 (1976).
Ali, K. & Bilal, S. Surface tension, surface excess concentration, enthalpy and entropy of surface formation of aqueous salt solutions. Colloids and Surfaces A: Physicochemical and Engineering Aspects 417, 183–190 (2013).
Neyt, J. C. et al. How does the electronic continuum model perform in the prediction of the surface tension of salt solutions? Chemical Physics Letters 595, 209–213 (2014).
Jiang, H. & Panagiotopoulos, A. Z. Note: Activity coefficients and solubilities for the NaCl/ε force field. The Journal of chemical physics 145, 046101 (2016).
Jiang, H., Mester, Z., Moultos, O. A., Economou, I. G. & Panagiotopoulos, A. Z. Thermodynamic and transport properties of H2O+ NaCl from polarizable force fields. Journal of chemical theory and computation 11, 3802–3810 (2015).
Abraham, M. et al. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1, 19–25 (2015).
Jorgensen, W. & Tirado-Rives, J. The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. Journal of the American Chemical Society 110, 1657–1666 (1988).
Abascal, J. & Vega, C. A general purpose model for the condensed phases of water: TIP4P/2005. Journal of Chemical Physics 123, 234505 (2005).
Smith, D. & Dang, L. Computer simulations of NaCl association in polarizable water. Journal of Chemical Physics 100, 3757–3766 (1994).
Lemmon, E., McLinden, M., Friend, D., Linstrom, P. & Mallard, W. NIST chemistry webbook, NIST standard reference database number 69. National Institute of Standards and Technology, Gaithersburg (2011).
Martnez, L., Andrade, R., Birgin, E. & Martnez, J. PACKMOL: a package for building initial configurations for molecular dynamics simulations. Journal of Computational Chemistry 30, 2157–2164 (2009).
Haynes, W. CRC Handbook of Chemistry and Physics, 92nd Edition. 100 Key Points (CRC Press), https://books.google.co.uk/books?id=pYPRBQAAQBAJ (2011).
Irving, J. & Kirkwood, J. The statistical mechanical theory of transport processes. iv. the equations of hydrodynamics. The Journal of Chemical Physics 18, 817–829 (1950).
Mark, P. & Nilsson, L. Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. Journal of Physical Chemistry A 105, 9954–9960 (2001).
Soper, A. & Phillips, M. A new determination of the structure of water at 25 C. Chemical Physics 107, 47–60 (1986).
Sega, M., Kantorovich, S., Jedlovszky, P. & Jorge, M. The generalized identification of truly interfacial molecules (ITIM) algorithm for nonplanar interfaces. Journal of Chemical Physics 138, 044110 (2013).
Jasper, J. & Kring, E. The isobaric surface tensions and thermodynamic properties of the surfaces of a series of n-alkanes, C5 to C18, 1-Alkenes, C6 to C16, and of n-Decylcyclopentane, n-Decylcyclohexane and n-Dcylbenzene. Journal of Physical Chemistry 59, 1019–1021 (1955).
Vargaftik, N., Volkov, B. & Voljak, L. International tables of the surface tension of water. Journal of Physical and Chemical Reference Data 12, 817–820 (1983).
Ali, K., Bilal, S. & Siddiqi, S. Concentration and temperature dependence of surface parameters of some aqueous salt solutions. Colloids and Surfaces A: Physicochemical and Engineering Aspects 272, 105–110 (2006).
Zeppieri, S., Rodrguez, J. & López de Ramos, A. Interfacial tension of alkane + water systems. Journal of Chemical & Engineering Data 46, 1086–1088 (2001).
Goebel, A. & Lunkenheimer, K. Interfacial tension of the water/n-alkane interface. Langmuir 13, 369–372 (1997).
Aveyard, R. & Haydon, D. Thermodynamic properties of aliphatic hydrocarbon/water interfaces. Transactions of the Faraday Society 61, 2255–2261 (1965).
Author information
Authors and Affiliations
Contributions
T.R.U. conceived, conducted, analysed the simulations, and drafted the initial manuscript. H.C.G. helped conceive the simulations. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Underwood, T.R., Greenwell, H.C. The Water-Alkane Interface at Various NaCl Salt Concentrations: A Molecular Dynamics Study of the Readily Available Force Fields. Sci Rep 8, 352 (2018). https://doi.org/10.1038/s41598-017-18633-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-18633-y
This article is cited by
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.