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.     

Principles of Surface Potential Estimation in Mixed Electrolyte Solutions:Taking into Account Dielectric Saturation

The dielectric properties between in-particle/water interface and bulk solution are significantly different, which are ignored in the theories of surface potential estimation. The analytical expressions of surface potential considering the dielectric saturation were derived in mixed electrolytes based on the nonlinear Poisson-Boltzmann equation. The surface potentials calculated from the approximate analytical and exact numerical solutions agreed with each other for a wide range of surface charge densities and ion concentrations. The effects of dielectric saturation became important for surface charge densities larger than 0.30 C/m\begin{document}$ ^2 $\end{document}. The analytical models of surface potential in different mixed electrolytes were valid based on original Poisson-Boltzmann equation for surface charge densities smaller than 0.30 C/m\begin{document}$ ^2 $\end{document}. The analytical model of surface potential considering the dielectric saturation for low surface charge density can return to the result of classical Poisson-Boltzmann theory. The obtained surface potential in this study can correctly predict the adsorption selectivity between monovalent and bivalent counterions.     

Determination of stagnant layer conductivity in polystyrene suspensions: temperature effects

In the classical theory of electrokinetic phenomena, it is admitted that the whole electrokinetic behavior of any colloidal system is fully determined by the zeta potential, zeta, of the interface. However, both experimental data and theoretical models have shown that this is an incomplete picture, as ions in the stagnant layer (the region between the solid surface and the slip plane--the plane where the equilibrium potential equals zeta) may respond to the field. In this paper, we aim at the evaluation of this contribution by the estimation of both K(SL)(sigma) (the surface conductivity of the stagnant layer) and K(d)(sigma) (the conductivity associated with the diffuse layer). This will be done by measuring the high-frequency dielectric dispersion (HFDD) in polystyrene suspensions; here "high-frequency" means the frequency interval where Maxwell-Wagner-O'Konski relaxation takes place (typically at MHz frequencies). Prior to any conclusions, a treatment of electrode polarization effects in the measurements was needed: we used two methods, and both led to similar results. Simulating the existence of surface conductivity by bulk conductivity, we reached the conclusion that no consistent explanation can be given for our HFDD and additional electrophoresis data based on the existence of diffuse-layer conductivity alone. We thus show how K(SL)(sigma) can be estimated and demonstrate that it can be explained by an ionic mobility very close to that characteristic of ions in the bulk solution. Such mobility, and hence also the values of K(SL)(sigma), increases with temperature as expected on simple physical grounds.     

INTERFACIAL FREE ENERGIES,ZETA POTENTIALS AND THROMBORESI STANCE IN BOVINE AND OVINE CAROTID ARTERIES

The physical properties of the intimal surfaces of both natural and artificial arteries are important In the prevention of thrombosis and clotting. It Is expected thati a) the interface between the blood plasma and the intima should be highly negatively charged, to repel negatively charged platelets and other blood constituents, and b) the interfacial free energy at the plasma-intima interface should be as low as possible, to prevent adhesion If the charge barrier is breached. Few measurements of either the zeta potential or the interfacial free energy at the intimal surface can be found in the literature. Zeta potentials of natural and cross-linked bovine carotid intima, and of an expanded Teflon prosthesis were determined by a streaming potential method in 0.9% (w/v) NaC1. Interfacial free energies of both bovine and ovine carotid intima, and of an expanded Teflon prosthesis were also determined by a contact-angle method using water immiscible liquid drops on the intimal surface Immersed in 0.9% (w/v) NaC1.     

Structure of the surface layer of mercury and adhesion of Langmuir monolayers on it

A model crystalline approximant of the surface layer (SL) structure of liquid mercury is formed using the algorithm, previously applied to form a crystalline approximant of the SL structure of water, the experimental data on the SL structure of mercury, and the stereochemical data for Hg atoms in the α-Hg crystalline phase and bulk liquid mercury. The correspondence of the metrics of the crystalline approximant proposed for the SL structure of mercury and Langmuir monolayers composed of organic molecules with hydrocarbon tails and different hydrophilic groups on liquid mercury to the position of molecules along mercury SL is shown.     

Sensitivity analysis of thermodynamic properties of liquid water: a general approach to improve empirical potentials

A sensitivity analysis of bulk water thermodynamics is presented in an effort to understand the relation between qualitative features of molecular potentials and properties that they predict. The analysis is incorporated in molecular dynamics simulations and investigates the sensitivity of the Helmholtz free energy, internal energy, entropy, heat capacity, pressure, thermal pressure coefficient, and static dielectric constant to components of the potential rather than the parameters of a given functional form. The sensitivities of the properties are calculated with respect to the van der Waals repulsive and the attractive parts, plus short- and long-range Coulomb parts of three four site empirical water potentials: TIP4P, Dang-Chang and TTM2R. The polarization sensitivity is calculated for the polarizable Dang-Chang and TTM2R potentials. This new type of analysis allows direct comparisons of the sensitivities for different potentials that use different functional forms. The analysis indicates that all investigated properties are most sensitive to the van der Waals repulsive, the short-range Coulomb and the polarization components of the potentials. When polarization is included in the potentials, the magnitude of the sensitivity of the Helmholtz free energy, internal energy, and entropy with respect to this part of the potential is comparable in magnitude to the other electrostatic components. In addition similarities in trends of observed sensitivities for nonpolarizable and polarizable potentials lead to the conclusion that the complexity of the model is not of critical importance for the calculation of these thermodynamic properties for bulk water. The van der Waals attractive and the long-range Coulomb sensitivities are relatively small for the entropy, heat capacity, thermal pressure coefficient and the static dielectric constant, while small changes in any of the potential contributions will significantly affect the pressure. The analysis suggests a procedure for modification of the potentials to improve predictions of thermodynamic properties and we demonstrate this general approach for modifying potentials for one of the potentials.     

Surface electronic structure of periodically δ-doped superlattice

Surface electronic structure of a semi-infinite superlattice (SL) intentionally δ-doped in each well- or barrier-layer is studied using a terminated Kronig-Penney-type of model interspersed with a periodic array of δ-function potentials representing the deep-center defect sheets. It is shown that SL minibands can be tailored by varying the properties of inserted δ-defects, e.g., the lowest energetic miniband can be shifted below the bottom of SL quantum wells (as a result, the main SL energy gap becomes smaller than those of the host semiconductors) and miniband alignement can occur (as a result, the bandwith is enhanced by an order of magnitude). Special attention is paid to the effect of δ-doping on the properties of surface states (SS's), i.e., the states appearing within energetic minigaps and localized at the SL/substrate interface. It is found that by an appropriate choice of the δ-defect weight and location within the SL period, SS with a required energy position relative to minigap edges and a desired extension into the SL bulk can be achieved. Moreover, in contrast to standard (i.e., undoped) SL's, a possibility of SS existence for the substrate identical to SL barriers is shown in periodically δ-doped SL's.     

Influence of membrane layer properties on the electrophoretic behavior of a soft particle

The influence of the physical properties of the membrane layer of a soft particle, which comprises a rigid core and a porous membrane layer, on its electrophoretic behavior, is investigated. Because that influence was almost always neglected in the previous studies, the corresponding results can be unrealistic. The applicability of the model proposed is verified by the available theoretical and experimental results. The electrophoretic mobility of the particle under various conditions is simulated through varying the dielectric constant, the thickness, and the drag coefficient of the membrane layer, and the bulk ionic concentration. We show that under typical conditions, the deviation in the electrophoretic mobility arising from assuming that the dielectric constant of the membrane layer is the same as that of the bulk liquid phase can be in the order of 50%. In addition, the thicker the membrane layer and/or the higher the bulk ionic concentration, the larger the deviation. If the surface of the core of the particle is charged, as in the case of inorganic particles covered by an artificial membrane layer, the deviation at constant core surface potential is larger than that under other types of charged conditions. However, if the surface of the core is uncharged, as in the case of biocolloids, then that deviation becomes negligible. These findings are of fundamental significance to theoreticians in their analysis on the electrokinetic behaviors of soft particles, and to experimentalists in the interpretation of their data.     

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.     

Electrowetting-based control of static droplet states on rough surfaces

Electrowetting (EW) is a powerful tool to control fluid motion at the microscale and has promising applications in the field of microfluidics. The present work analyzes the influence of an electrowetting voltage in determining and altering the state of a static droplet resting on a rough surface. An energy-minimization-based modeling approach is used to analyze the influence of interfacial energies, surface roughness parameters, and electric fields in determining the apparent contact angle of a droplet in the Cassie and Wenzel states under the influence of an EW voltage. The energy-minimization-based approach is also used to analyze the Cassie-Wenzel transition under the influence of an EW voltage and estimate the energy barrier to transition. The results obtained show that EW is a powerful tool to alter the relative stabilities of the Cassie and Wenzel states and enable dynamic control of droplet morphology on rough surfaces. The versatility and generalized nature of the present modeling approach is highlighted by application to the prediction of the contact angle of a droplet on an electrowetted rough surface consisting of a dielectric layer of nonuniform thickness.     

Coarse-graining in interaction space: a systematic approach for replacing long-range electrostatics with short-range potentials

A short-range effective potential for long-range electrostatic interactions in homogeneously disordered condensed phase systems has been determined with a novel approach to coarse-graining in interaction space. As opposed to coarse-graining the system resolution, this approach "coarsens" the system's interactions by mapping multiple configurations of an accurate long-range atomistic potential onto a more efficient, short-range effective potential with a force-matching (FM) method. Developing an empirical potential in this manner is fundamentally different from existing strategies because it utilizes condensed-phase (as opposed to gas-phase) atomistic interactions to determine general pair potentials defined on distance meshes (as opposed to fitting predetermined functional forms). The resulting short-range ( approximately 10 A) effective potential reproduces structural, dynamical, and many thermodynamic properties of liquid water, ions in water, and hydrophobes in water, with unprecedented accuracy. The effective potential is also shown to be transferable to a nonaqueous molten salt system. With continued development, such effective potentials may provide an accurate and highly efficient alternative to Ewald-based long-range electrostatics methods.     

Electrokinetic actuation of an uncharged polarizable dielectric droplet in charged hydrogel medium

Electrokinetic transport of an uncharged nonconducting microsized liquid droplet in a charged hydrogel medium is studied. Dielectric polarization of the liquid drop under the action of an externally imposed electric field induces a non-homogeneous charge density at the droplet surface. The interactions of the induced surface charge of the droplet with the immobile charges of the hydrogel medium generates an electric force to the droplet, which actuates the drop through the charged hydrogel medium. A numerical study based on the first principle of electrokinetics is adopted. Dependence of the droplet velocity on its dielectric permittivity, bulk ionic concentration, and immobile charge density of the gel is analyzed. The surface conduction is significant in presence of charged gel, which creates a concentration polarization. The impact of the counterion saturation in the Debye layer due to the dielectric decrement of the medium is addressed. The modified Nernst–Planck equation for ion transport and the Poisson equation for the electric field is considered to take into account the dielectric polarization. A quadrupolar vortex around the uncharged droplet is observed when the gel medium is considered to be uncharged, which is similar to the induced charge electroosmosis around an uncharged dielectric colloid in free-solution. We find that the induced charge electrokinetic mechanism creates a strong recirculation of liquid within the droplet and the translational velocity of the droplet strongly depends on its size for the dielectric droplet embedded in a charged gel medium.     

A new realistic potential for neon based on dilute gas bulk properties

A potential based on the simple Maitland-Smith generalized n-6 form and fit solely to three dilute gas bulk properties is proposed. The potential predicts not only dilute gas bulk properties but also spectroscopic and molecular beam data as favorably as much more complex realistic potentials.     

Effective potential approach to bulk thermodynamic properties and surface tension of molecular fluids I. Single component n-alkane systems

The aim of this work is to develop spherically symmetric effective potentials allowing bulk thermodynamic properties and surface tension of molecular fluids to be predicted semiempirically by the use of statistical mechanical methods. Application is made to the straight chain alkane fluids from methane to decane. An effective Lennard-Jones potential is generated with temperature-dependent parameters fitted to the critical temperature and pressure and to Pitzer's acentric factor. Insertion of this potential into the generalised van der Waals (GvdW) density functional theory yields bulk properties in good agreement with experiments. The surface tension is overestimated for the longer alkane chains. In order to account for the surface tension, an independently adjustable attractive range of interaction is required and obtained through the use of square-well potentials chosen so as to leave the bulk thermodynamics unaltered while the attractive range is fitted to the surface tension at a single temperature. The GvdW theory, which includes binding energy, entropic and profile shape contributions, then generates surface tension estimates that are of good accuracy over the full range of available experimental data. It appears that, given a sufficiently flexible form, effective potentials combined with simple statistical mechanical theory can reproduce both bulk and non-uniform fluid data of great variety in an insighful and practically useful way.     

A polarizable force field for water using an artificial neural network

A force field for liquid water including polarization effects has been constructed using an artificial neural network (ANN). It is essential to include a many-body polarization effect explicitly into a potential energy function in order to treat liquid water which is dense and highly polar. The new potential energy function is a combination of empirical and nonempirical potentials. The TIP4P model was used for the empirical part of the potential. For the nonempirical part, an ANN with a back-propagation of error algorithm (BPNN) was introduced to reproduce the complicated many-body interaction energy surface from ab initio quantum mechanical calculations. BPNN, described in terms of a matrix, provides enough flexibility to describe the complex potential energy surface (PES). The structural and thermodynamic properties, calculated by isobaric-isothermal (constant-NPT) Monte Carlo simulations with the new polarizable force field for water, are compatible with experimental results. Thus, the simulation establishes the validity of using our estimated PES with a polarization effect for accurate predictions of liquid state properties. Applications of this approach are simple and systematic so that it can easily be applied to the development of other force fields besides the water-water system.     

On the behaviour of molecules of the quinoline group at the water—Mercury interface: Part III. The kinetics of the isoquinoline transition in the presence of an “inert” electrolyte

The reorientation transition, which involves two distinct monolayers of isoquinoline molecules, has been investigated by using the single-potential step and the double-potential step methods under various experimental conditions.The cathodic transients correspond to a liquid → solid transition. Their morphology and half-times have been measured at various initial and final potentials, and for several isoquinoline concentrations close to the saturation value. At small overvoltages, the reorientation process is determined by heterogeneous (“instantaneous”) nucleation, while progressive nucleation and growth prevail at higher potentials. Recourse to the double-potential step method affords an efficient way of assessing separately the influence of the overvoltage on the rate constants for nucleation and growth. Determination of the critical nucleus size and the activation energy for each of the two processes has been based on an adaptation of the Brandes theory developed for two-dimensional crystallization from a supersaturated vapour phase. The kinetics are markedly dependent on the initial state, which can be easily controlled by the potential and the bulk concentration. When the potential is stepped from the region where there is superadsorption, progressive nucleation is particularly fast, while starting from a partially depleted layer gives much slower transients with extensive tailing due to diffusional hindrance.The anodic transients show consonant characteristics which indicate that the initial film can be considered as a two-dimensional solid, and that boundary defects are acting as nucleation centres, during the solid → liquid transition triggered by the potential step.     

Tangential streaming potential as a tool in modeling of ion transport through nanoporous membranes

Tangential streaming potential (TSP) measurements have been carried out so as to assess the electrokinetic properties of the active layer of organic nanofiltration (NF) membranes. Due to the porous structure of NF membranes, cares must be taken to convert the experimental data into zeta potential. Indeed, an assumption that is implicitly made in Smoluchowski's theory (or in related approaches accounting for the surface conduction phenomenon) is that both streaming and conduction currents involved in the streaming potential process flow through an identical path. Such an assumption does not hold with porous membranes since the conduction current is expected to flow wherever the electric conductivity differs from zero. Consequently, a non-negligible share of the conduction current is likely to flow through the membrane body filled with the electrolyte solution. This phenomenon has been taken into account by carrying out a series of TSP measurements at various channel heights. Experiments have been conducted with various electrolyte solutions. The inferred zeta potentials have been further converted into membrane volume charge densities which have been used to predict the membrane performances in terms of rejection rates. The conventional NF theory, i.e. based on a steric/Donnan exclusion mechanism, has been found to be unable to describe the experimental rejection rates. Using the volume charge density of the membrane as an adjustable parameter, it has been shown that the conventional theory even predicts the opposite sign for the membrane charge. On the other hand, the experimental rejection rates have been well described by including dielectric effects in the exclusion mechanism. In this case, a noticeable lowering of the effective dielectric constant of the electrolyte solution inside pores has been predicted (with respect to the bulk value).     

Static electrification of solid oxide in liquid metal and electrical double layer at the interface

Static electrification of a solid oxide, say a semiconducting oxide in liquid metal, is mainly due to electron transfer between two phases. Excess electrons in the liquid metal phase provided by the oxide give rise to an electrical double layer at the interface. The electrical double layer may be divided into three parts, an immobile inner layer, a compressed diffuse layer, and a flat layer extending into the bulk liquid metal. Differential potential analysis and the induced emf method were used to measure the potential of the compressed diffuse layer and the excess electron density of the flat layer, respectively. Results show that most oxides in liquid metals carry positive charges on their surfaces and the potentials of the compressed diffuse layer are in the range of 3 to 42 microV. Such a low potential implies that the diffuse layer is considerably compressed. The excess electron densities of the flat layer are on the order of 10(22) electrons/m(3) of Hg and their contributions to surface charges of oxide are in the range of 10(17) to 10(18) charges/m(2) for the oxide/mercury systems with a solid density of 0.3 wt% at room temperature.