Revisiting the hexane-water interface via molecular dynamics simulations using nonadditive alkane-water potentials

We present results addressing properties of a polarizable force field for hexane based on the fluctuating charge (FQ) formalism and developed in conjunction with the Chemistry at Harvard Molecular Mechanics (CHARMM) potential function. Properties of bulk neat hexane, its liquid-vapor interface, and its interface with a polarizable water model (TIP4P-FQ) are discussed. The FQ model is compared to a recently modified alkane model, C27r, also based on the CHARMM potential energy function. With respect to bulk properties, both models predict bulk density within 1%; the FQ model predicts the liquid vaporization enthalpy within 2%, while the C27r force field underestimates the property by roughly 20% (and in this sense reflects the quality of the C27r force field across the spectrum of linear and branched alkanes). The FQ hexane model realistically captures the dielectric properties of the bulk in terms of a dielectric constant of 1.94, in excellent agreement with experimental values in the range of 1.9-2.02. This behavior is also in conformity with a recent polarizable alkane model based on Drude oscillators. Furthermore, the bulk dielectric is essentially captured in the infinite frequency, or optical, dielectric contribution. The FQ model is in this respect a more realistic force field for modeling lipid bilayer interiors for which most current state-of-the-art force fields do not accurately capture the dielectric environment. The molecular polarizability of the FQ model is 11.79 A3, in good agreement with the range of experimental and ab initio values. In contrast to FQ models of polar solvents such as alcohols and water, there was no need to scale gas-phase polarizabilities in order to avoid polarization catastrophes in the pure bulk. In terms of the liquid-vapor and liquid-liquid interfaces, the FQ model displays a rich orientational structure of alkane and water in the respective interfacial systems, in general conforming with earlier simulation studies of such interfaces. The FQ force field shows a marked deviation in the interfacial dipole potentials computed from the charge densities averaged over simulation trajectories. At the liquid-vapor interface, the FQ model predicts a potential drop of -178.71 mV in contrast to the C27r estimate of -433.80 mV. For the hexane-water interface, the FQ force field predicts a dipole potential drop of -379.40 mV in contrast to the C27r value of -105.42 mV. Although the surface dipole potential predicted by the FQ model is roughly 3.5 times that predicted by the C27r potential, it is consistent with reported experimental potentials across solvated lipid bilayers in the range of 400-600 mV.     

Capillary waves at the liquid-vapor interface and the surface tension of water

Capillary waves occurring at the liquid-vapor interface of water are studied using molecular dynamics simulations. In addition, the surface tension, determined thermodynamically from the difference in the normal and tangential pressure at the liquid-vapor interface, is compared for a number of standard three- and four-point water models. We study four three-point models (SPC/E, TIP3P, TIP3P-CHARMM, and TIP3P-Ew) and two four-point models (TIP4P and TIP4P-Ew). All of the models examined underestimate the surface tension; the TIP4P-Ew model comes closest to reproducing the experimental data. The surface tension can also be determined from the amplitude of capillary waves at the liquid-vapor interface by varying the surface area of the interface. The surface tensions determined from the amplitude of the logarithmic divergence of the capillary interfacial width and from the traditional thermodynamic method agree only if the density profile is fitted to an error function instead of a hyperbolic tangent function.     

Structure, thermodynamics, and liquid-vapor equilibrium of ethanol from molecular-dynamics simulations using nonadditive interactions

We present a molecular-dynamics simulation study of the bulk and liquid-vapor interfacial properties of ethanol using a polarizable force field based on the fluctuating charge (FQ) formalism, as well as the nonpolarizable CHARMM22 force field. Both models are competitive with respect to the prediction of ambient liquid properties such as liquid density, enthalpy of vaporization, dielectric constant, and self-diffusion constants. The polarizable model predicts an average condensed-phase dipole moment of 2.2 D associated with an induced liquid-phase dipole moment of 0.6 D; though qualitatively in agreement with earlier nonadditive models as well as recent Car-Parinello calculations, the current FQ model underestimates the condensed-phase dipole moment. In terms of liquid structure, both models are in agreement with recent neutron-diffraction results of liquid ethanol structure, although the polarizable model predicts the hydroxyl-hydrogen-hydroxyl-hydrogen structure factor in closer agreement with the experimental data. In terms of interfacial properties, both models predict ambient surface tension to within 4% of the experimental value of 22.8 dyncm, while overestimating the surface excess entropy by almost a factor of 2. Both models display the characteristic preferential orientation of interfacial molecules. The polarizable model allows for a monotonic variation of the average molecular dipole moment from the bulk value to that of the vapor phase. Consequently, there is a dramatic difference in the surface potential predicted by the polarizable and nonpolarizable models. The polarizable model estimates a surface potential of -209+/-3 mV, while the nonpolarizable model yields a value of -944+/-10 mV. Finally, based on the vapor-liquid equilibrium simulation data from several temperatures, we estimate the critical properties of both models. As observed with other FQ models for associating fluids (such as water and methanol), and counter to what one would anticipate by modeling more physically the electrostatic response to local environment, the current FQ model underestimates the critical temperature and overestimates the critical density of ethanol; moreover, the FQ model is, in this respect, equivalent to the underlying fixed-charge model. These results further suggest the need to revisit polarizable models in terms of quantitative vapor-liquid equilibrium prediction.     

Molecular dynamics simulations of the surface tension of ionic liquids

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.     

Structure and dynamics of the aqueous liquid-vapor interface: a comprehensive particle-based simulation study

This research addresses a comprehensive particle-based simulation study of the structural, dynamic, and electronic properties of the liquid-vapor interface of water utilizing both ab initio (based on density functional theory) and empirical (fixed charge and polarizable) models. Numerous properties such as interfacial width, hydrogen bond populations, dipole moments, and correlation times will be characterized with identical schemes to draw useful conclusions on the strengths and weakness of the proposed models for interfacial water. Our findings indicate that all models considered in this study yield similar results for the radial distribution functions, hydrogen bond populations, and orientational relaxation times. Significant differences in the models appear when examining both the dipole moments and surface relaxation near the aqueous liquid-vapor interface. Here, the ab initio interaction potential predicts a significant decrease in the molecular dipole moment and expansion in the oxygen-oxygen distance as one approaches the interface in accordance with recent experiments. All classical polarizable interaction potentials show a less dramatic drop in the molecular dipole moment, and all empirical interaction potentials studied yield an oxygen-oxygen contraction as the interface is approached.     

Free energy partitioning analysis of the driving forces that determine ion density profiles near the water liquid-vapor interface

Free energy partitioning analysis is employed to explore the driving forces for ions interacting with the water liquid-vapor interface using recently optimized point charge models for the ions and SPC/E water. The Na(+) and I(-) ions are examined as an example kosmotrope/chaotrope pair. The absolute hydration free energy is partitioned into cavity formation, attractive van der Waals, local electrostatic, and far-field electrostatic contributions. We first compute the bulk hydration free energy of the ions, followed by the free energy to insert the ions at the center of a water slab. Shifts of the ion free energies occur in the slab geometry consistent with the SPC/E surface potential of the water liquid-vapor interface. Then the free energy profiles are examined for ion passage from the slab center to the dividing surface. The profiles show that, for the large chaotropic I(-) ion, the relatively flat total free energy profile results from the near cancellation of several large contributions. The far-field electrostatic part of the free energy, largely due to the water liquid-vapor interface potential, has an important effect on ion distributions near the surface in the classical model. We conclude, however, that the individual forms of the local and far-field electrostatic contributions are expected to be model dependent when comparing classical and quantum results. The substantial attractive cavity free energy contribution for the larger I(-) ion suggests that there is a hydrophobic component important for chaotropic ion interactions with the interface.     

Surface properties of the polarizable Baranyai-Kiss water model

The water surface properties using the Baranyai-Kiss (BK) model [A. Baranyai and P. T. Kiss, J. Chem. Phys. 133, 144109 (2010)] are studied by molecular dynamics simulation, and compared to popular rigid water potentials, namely to the extended simple point charge (SPC/E) and the transferable interaction potential with 4 points (TIP4P) models. The BK potential is a polarizable model of water with three Gaussian charges. The negative charge is connected to its field-free position by a classical harmonic spring, and mechanical equilibrium is established between this spring force and the force due to the charge distribution of the system. The aim of this study is, on the one hand, to test the surface properties of the new model, and on the other hand, to identify differences between the models listed above. The obtained results reveal that the BK model reproduces very well a number of properties corresponding to liquid-vapor equilibrium, such as the coexisting liquid and vapor densities, saturated vapor pressure or surface tension. Further, this model reproduces excellently the critical point of water even in comparison with a large number of widely used polarizable and nonpolarizable models. The structural properties of the liquid surface of BK water turns out to be very similar to that of the SPC/E model, while the surface of TIP4P water is found to be somewhat less ordered. This finding is related to the fact that the critical temperature of the TIP4P model is lower than that of either SPC/E or BK.     

The Wolf method applied to the liquid-vapor interface of water

The Wolf method for the calculation of electrostatic interactions is applied in a liquid phase and at the liquid-vapor interface of water and its results are compared with those from the Ewald sums method. Molecular dynamics simulations are performed to calculate the radial distribution functions at room temperature. The interface simulations are used to obtain the coexisting densities and surface tension along the coexistence curve. The water model is a flexible version of the extended simple point charge model. The Wolf method gives good structural results, fair coexistence densities, and poor surface tensions as compared with those obtained using the Ewald sums method.     

A study on orientation and absorption spectrum of interfacial molecules by using continuum model

In this work, a numerical procedure based on the continuum model is developed and applied to the solvation energy for ground state and the spectral shift against the position and the orientation of the interfacial molecule. The interface is described as a sharp boundary separating two bulk media. The polarizable continuum model (PCM) allows us to account for both electrostatic and nonelectrostatic solute-solvent interactions when we calculate the solvation energy. In this work we extend PCM to the interfacial system and the information about the position and orientation of the interfacial molecule can be obtained. Based on the developed expression of the electrostatic free energy of a nonequilibrium state, the numerical procedure has been implemented and used to deal with a series of test molecules. The time-dependent density functional theory (TDDFT) associated with PCM is used for the electron structure and the spectroscopy calculations of the test molecules in homogeneous solvents. With the charge distribution of the ground and excited states, the position- and orientation-dependencies of the solvation energy and the spectrum have been investigated for the interfacial systems, taking the electrostatic interaction, the cavitation energy, and the dispersion-repulsion interaction into account. The cavitation energy is paid particular attention, since the interface portion cut off by the occupation of the interfacial molecule contributes an extra part to the stabilization for the interfacial system. The embedding depth, the favorable orientational angle, and the spectral shift for the interfacial molecule have been investigated in detail. From the solvation energy calculations, an explanation has been given on why the interfacial molecule, even if symmetrical in structure, tends to take a tilting manner, rather than perpendicular to the interface.     

Surface properties of fluids of charged platelike colloids

Surface properties of mixtures of charged platelike colloids and salt in contact with a charged planar wall are studied within density functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model [J. Chem. Phys. 39, 1714 (1963)] and the charges of the particles are concentrated at their centers. The density functional applied is an extension of a recently introduced functional for charged platelike colloids. It provides a qualitative approach because it does not determine the relation between the actual and the effective charges entering into the model. Technically motivated approximations, such as using the Zwanzig model, are expected not to influence the results qualitatively. Analytically and numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large macroion charges and nematic bulk order. The asymptotic wetting and drying behaviors are investigated by means of effective interface potentials which turn out to be asymptotically the same as for a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying points as well as predrying lines and the corresponding critical points have been located numerically. A crossover from monotonic to nonmonotonic electrostatic potential profiles upon varying the surface charge density has been observed. Nonmonotonic electrostatic potential profiles are equivalent to the occurrence of charge inversion. Due to the presence of both the Coulomb interactions and the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, i.e., the point of zero charge and the isoelectric point of the surface do not coincide.     

Long-range Lennard-Jones and electrostatic interactions in interfaces: application of the isotropic periodic sum method

Molecular dynamics (MD) simulations of heptane/vapor, hexadecane/vapor, water/vapor, hexadecane/water, and dipalmitoylphosphatidylcholine (DPPC) bilayers and monolayers are analyzed to determine the accuracy of treating long-range interactions in interfaces with the isotropic periodic sum (IPS) method. The method and cutoff (rc) dependences of surface tensions, density profiles, water dipole orientation, and electrostatic potential profiles are used as metrics. The water/vapor, heptane/vapor, and hexadecane/vapor interfaces are accurately and efficiently calculated with 2D IPS (rc=10 A). It is demonstrated that 3D IPS is not practical for any of the interfacial systems studied. However, the hybrid method PME/IPS [Particle Mesh Ewald for electrostatics and 3D IPS for Lennard-Jones (LJ) interactions] provides an efficient way to include both types of long-range forces in simulations of large liquid/vacuum and all liquid/liquid interfaces, including lipid monolayers and bilayers. A previously published pressure-based long-range LJ correction yields results similar to those of PME/IPS for liquid/liquid interfaces. The contributions to surface tension of LJ terms arising from interactions beyond 10 A range from 13 dyn/cm for the hexadecane/vapor interface to approximately 3 dyn/cm for hexadecane/water and DPPC bilayers and monolayers. Surface tensions of alkane/vapor, hexadecane/water, and DPPC monolayers based on the CHARMM lipid force fields agree very well with experiment, whereas surface tensions of the TIP3P and TIP4P-Ew water models underestimate experiment by 16 and 11 dyn/cm, respectively. Dipole potential drops (DeltaPsi) are less sensitive to long-range LJ interactions than surface tensions. However, DeltaPsi for the DPPC bilayer (845+/-3 mV proceeding from water to lipid) and water (547+/-2 mV for TIP4P-Ew and 521+/-3 mV for TIP3P) overestimate experiment by factors of 3 and 5, respectively, and represent expected deficiencies in nonpolarizable force fields.     

Effects of ion solvation on phase equilibrium and interfacial tension of liquid mixtures

We study the bulk thermodynamics and interfacial properties of electrolyte solution mixtures by accounting for electrostatic interaction, ion solvation, and inhomogeneity in the dielectric medium in the mean-field framework. Difference in the solvation energy between the cations and anions is shown to give rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The ion solvation affects the phase equilibrium of the solvent mixture, depending on the dielectric constants of the solvents, reflecting the competition between the solvation energy and translation entropy of the ions. Miscibility is decreased if both solvents have low dielectric constants and is enhanced if both solvents have high dielectric constant. At the mean-field level, the ion distribution near the interface is determined by two competing effects: accumulation in the electrostatic double layer and depletion in a diffuse interface. The interfacial tension shows a nonmonotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations and decreases approximately as the square root of the salt concentration for dilute solutions, reaching a minimum near 1 mM. We also find that, for a fixed cation type, the interfacial tension decreases as the size of anion increases. These results offer qualitative explanations within one unified framework for the long-known concentration and ion size effects on the interfacial tension of electrolyte solutions.     

Effect of SDBS solution on the surface potential properties of the Zhaozhuang coal

It is of great significance to study the effect of surfactants on the coal surface potential for effective dust suppression in mining faces. The effect of different concentrations of sodium dodecyl benzene sulfonate (SDBS) solution on the surface potential of the Zhaozhuang coal was measured by atomic force microscopy (AFM). The experimental results show that the SDBS solution has significant influence on the surface potential of the Zhaozhuang coal. The electrical characteristics of the coal surface at the nanometer scale are different from those of macroscopic or the mesoscopic level. The surface potential of coal is basically a normal distribution, showing certain electrical characteristics. The mean value of the surface potential of the Zhaozhuang coal is increased with the increase in concentration of the SDBS solution; when the concentration of the SDBS solution is 0.3%, the mean value of surface potential is 5.59 mv, which is about two times of the mean value of the surface potential without SDBS solution added. With the increase of the concentration of the SDBS solution, the maximum value of the surface potential of the Zhaozhuang coal increases, and the minimum value decreases. It shows that the SDBS solution has a significant effect on the potential distribution law and the magnitude of the coal surface. Subsequently, on the basis of the constructed Zhaozhuang coal macromolecule model, xtb and Multiwfn simulation software were used to calculate the molecular surface electrostatic potential value and electron density value of the Zhaozhuang coal molecules after adding water molecules. The variation law for the electrostatic potential surface of the molecule was obtained after adding numbers of water molecules to the surface of the coal molecule. The simulation results show that the area proportion of absolute molecular surface electrostatic potential greater than 10 kcal/mol is increased with the growth in the number of water molecules, while the area proportion of absolute molecular surface electrostatic potential less than 10 kcal/mol is decreased. Because of the free state O─H bond polarity of water molecules, the charges on the molecular surface are rearranged in order to change the electron density on the surface of coal molecules, which affects the overall electrostatic potential of the configuration.     

Thermodynamic and structural properties of methanol-water solutions using nonadditive interaction models

We study bulk structural and thermodynamic properties of methanol-water solutions via molecular dynamics simulations using novel interaction potentials based on the charge equilibration (fluctuating charge) formalism to explicitly account for molecular polarization at the atomic level. The study uses the TIP4P-FQ potential for water-water interactions, and the CHARMM-based (Chemistry at HARvard Molecular Mechanics) fluctuating charge potential for methanol-methanol and methanol-water interactions. In terms of bulk solution properties, we discuss liquid densities, enthalpies of mixing, dielectric constants, self-diffusion constants, as well as structural properties related to local hydrogen bonding structure as manifested in radial distribution functions and cluster analysis. We further explore the electronic response of water and methanol in the differing local environments established by the interaction of each species predominantly with molecules of the other species. The current force field for the alcohol-water interaction performs reasonably well for most properties, with the greatest deviation from experiment observed for the excess mixing enthalpies, which are predicted to be too favorable. This is qualitatively consistent with the overestimation of the methanol-water gas-phase interaction energy for the lowest-energy conformer (methanol as proton donor). Hydration free energies for methanol in TIP4P-FQ water are predicted to be -5.6 +/- 0.2 kcal/mol, in respectable agreement with the experimental value of -5.1 kcal/mol. With respect to solution microstructure, the present cluster analysis suggests that the microscale environment for concentrations where select thermodynamic quantities reach extremal values is described by a bipercolating network structure.     

Poisson-Boltzmann theory of pH-sensitive (annealing) polyelectrolyte brush

We present a self-consistent field analytical theory of a polymer brush formed by weakly charged pH-sensitive (annealing) polyelectrolytes tethered to a solid-liquid interface and immersed in buffer solution of low molecular weight salt. We use the Poisson-Boltzmann framework, applied by us previously to polyelectrolyte (PE) brushes with quenched charge (Zhulina, E. B.; Borisov, O. V. J. Chem. Phys. 1997, 107, 5952). This approach allows for detailed analysis of the internal structure of annealing PE brush in terms of polymer density distribution, profiles of electrostatic potential and of local degree of chain ionization as a function of buffer ionic strength and pH without any assumptions on mobile ion distribution imposed in earlier scaling-type models. The presented analytical theory recovers all major asymptotic dependences for average brush properties predicted earlier. In particular, a nonmonotonic dependence of brush thickness on ionic strength and grafting density is confirmed and specified with accuracy of numerical coefficients including crossover regions. Moreover, the theory predicts qualitatively new effects, such as, e.g., disproportionation of tethered polyions into weakly charged concentrated proximal and strongly charged sparse distal brush domains at low salt and moderate grating densities. The presented results allow us to quantify responsive features of annealing PE brushes whose large-scale and local conformational properties can be manipulated by external stimuli.     

Surface tension and scaling of critical nuclei in diatomic and triatomic fluids

Density functional theory has been used to investigate surface tension and scaling of critical clusters in fluids consisting of diatomic and rigid triatomic molecules. The atomic sites are hard spheres with attractive interactions obtained from the tail part of the Lennard-Jones potential. Asymmetry in attractive interactions between the atomic sites has been introduced to cause molecular orientation and oscillatory density profiles at liquid-vapor interfaces. The radial dependence of cluster surface tension in fluids showing modest orientation in unimolecular layer at the interface or no orientation at all resembles the surface tension behavior of clusters in simple monoatomic fluids, although the surface tension maximum becomes more pronounced with increasing chain length of the molecule. Surface tension of clusters having multiple oscillatory layers at the interface shows a prominent maximum at small cluster sizes; however, the surface tension of large clusters is lower than the planar value. The scaling relation for the number of molecules in the critical cluster and the nucleation barrier height developed by McGraw and Laaksonen [Phys. Rev. Lett. 76, 2754 (1996)] are well obeyed for fluids with little structure at liquid-vapor interface. However, fluids having enhanced interfacial structure show some deviation from the particle number scaling, and the barrier height scaling breaks up seriously.     

Intermolecular interactions and electrostatic properties of the β-hydroquinone apohost: implications for supramolecular chemistry

The crystal structure of the β-polymorph of hydroquinone (β-HQ), the apohost of a large family of clathrates, is reported with a specific focus on intermolecular interactions and the electrostatic nature of its cavity. Hirshfeld surface analysis reveals subtle close contacts between two interconnecting HQ networks, and the local packing and related close contacts were examined by breakdown of the fingerprint plot. An experimental multipole model containing anisotropic thermal parameters for hydrogen atoms has been successfully refined against 15(2) K single microcrystal synchrotron X-ray diffraction data. The experimental electron density model has been compared with a theoretical electron density calculated with the molecule embedded in its own crystal field. Hirshfeld charges, interaction energies and the electrostatic potential calculated for both models are qualitatively in good agreement, but small differences in the electrostatic potential persist due to charge transfer from all hydrogen atoms to the oxygen atoms in the theoretical model. The electrostatic potential in the center of the cavity is positive, very shallow and highly symmetric, suggesting that the inclusion of polar molecules in the void will involve a balance between opposing effects. The electric field is by symmetry zero in the center of the cavity, increasing to a value of 0.0185 e/?(2) (0.27 V/?) 1 ? along the 3-fold axis and 0.0105 e/?(2) (0.15 V/?) 1 ? along the perpendicular direction. While these values are substantial in a macroscopic context, they are quite small for a molecular cavity and are not expected to strongly polarize a guest molecule.     

Computer simulation studies of surfactant monolayer mixtures at the water/oil interface: charge distribution effects

Charge distribution effects on polar head groups for a mixture of amphiphilic molecules at the water/oil interface were studied. For this purpose a model which allowed us to investigate the charge effects exclusively was created. As a molecular model we used the structure of sodium dodecyl sulfate. Then we prepared molecules with the same molecular structure but with different charge distributions in order to have one cationic and one nonionic molecule. So, in this way, we were able to focus only in the charge effects. The monolayer mixtures were composed of anionic/nonionic and cationic/nonionic surfactants. Simulations of these systems show that the location of the different surfactants at the interface is determined by the interaction and the charge distribution of the molecules. Due to the difference in the charge distribution of the surfactant monolayers, the water molecules present distinct orientations in the mixture. Finally, it was found that the electrostatic potential difference across the interface depended on the interactions (charge distribution) of the anionic, cationic, and nonionic molecules in the mixture.     

反胶束系统及萃取蛋白质的分子热力学模型1: 反胶束系统的分 子热力学模型

以反胶束系统稳定性的热力学分析为基础,综合分析了反胶束系统的三大效应,即低界面张力效应、界面弯矩效应、混合熵效应,提出了一个分子热力学模型,模型所预言的反胶束水分含量随无机盐种类、浓度、表面活性剂浓度以及助表面活性剂含量的变化与所获实验规律定量相符,还能预言反胶束内表面处电势值、表面活性剂解离度。