The line tension and the generalized Young equation: the choice of dividing surface

Using Gibbs method of dividing surfaces, the condition of equilibrium of a sessile drop on a flat non-deformable solid substrate is investigated. The dependence of the line tension on the curvature radius of the dividing three-phase contact line is found. It has been derived a relationship between the partial derivative of the line tension with respect to the curvature radius of the three-phase contact line (which stands in the generalized Young equation) and the total derivative of the line tension with respect to the same radius along the equilibrium states. Various approximated formulas of the generalized Young equation used in the literature are analyzed.     

Liquid-vapor equilibrium in the dimethyl sulfoxide-methanol system

The liquid-vapor phase equilibrium in the dimethyl sulfoxide-methanol binary system was studied with the use of a statistical method. Partial pressures of dimethyl sulfoxide and methanol were calculated by integrating the Gibbs-Duhem equation. Molar excess Gibbs energies were described by the Redlich-Kister equation, and correlation parameters were calculated. It was found that molar excess Gibbs energies are negative, and the deviation from ideality increases as temperature increases.     

Dynamic interfacial effect of electroosmotic slip flow with a moving capillary front in hydrophobic circular microchannels

Miniaturization of chemical analysis using microfabrication is an emerging technology. The use of polymeric materials as opposed to conventional glass substrate is also a promising alternative. As most polymeric materials are hydrophobic relative to glass, we describe here the implication for the loading process of electroosmotic flow (EOF) when a three-phase (solid-liquid-vapor) contact line exists. The presence of these interfaces can result in a large Laplace pressure that resists EOF and hence hinders its flow performance. This effect depends on the phenomenological contact angle at the solid-liquid interface. In our model for EOF, we considered simultaneously the presence of an electric double layer, liquid slips via a weaker solid-liquid interaction and Laplace pressure across a liquid-vapor interface.     

Monte Carlo simulation methodology of the ghost interface theory for the planar surface tension

A novel "ghost interface" expression for the surface tension of a planar liquid-vapor interface is derived in detail from consideration of the free energy of the system, and a methodology for utilization of this new technique is given. An augmented Monte Carlo computer simulation procedure is developed specifically for the ghost interface, including derivation of long-range corrections resulting from potential truncation and a modified Gibbs ensemble technique for the simulation of adjacent coexisting phases. Results generated from the ghost interface theory for the surface tension are presented and found to be in good quantitative agreement with those resulting from the Kirkwood-Buff equation. Applications of this new approach to curved and to supersaturated systems are also discussed.     

Generalized theory of capillarity for axisymmetric menisci

A high-curvature generalization of the Laplace equation of capillarity and the Young equation of capillarity (including line tension) is developed for an axisymmetric solid-liquid-fluid system. The most general expressions for the Laplace and Young equations do not assume a particular form for the specific surface free energy. However, when a particular form, i.e., ω^(A) = γ^(A)_∞+ C_JJ+ C_kK, which is related to Gibbs' expression for a highly curved menisci,¹ is assumed to hold for the specific surface free energy then we are able to recover the expected simplified form of the Laplace equation. The corresponding high-curvature Young equation includes a couple which balances the surface moments at the contact line. Unfortunately, the effect of this couple could be confused with the effect of line tension in experiments which attempt to measure the dependence of the contact angle on the contact line radius.     

Thermodynamic approach to calculating the surface tension of single-component liquids by computer simulations

A thermodynamic method for computing the surface tension at a flat liquid-vapor interface by the Monte Carlo or molecular dynamics methods over a wide temperature range was proposed. The approach is based on the Gibbs separating surface method; it does not require information on the mechanical state of the surface layer.     

The surface tension of a solid at the solid-vacuum interface, an evaluation from adsorption and wall potential calculations

A method for the evaluation of quantities that are experimentally inaccessible such as the surface tension at the solid-vacuum interface and the superficial tension of the fluid in contact with the solid is presented. The approach is based on consideration of an equilibrium of a fluid in solid capillary wherein a balance between surface and capillary forces has been replaced by conceptual alternative interfacial and centrifugal forces. This approach involves the simultaneous numerical solution one the special forms of the Gibbs equation for solid-fluid interface and a generalized Kelvin equation derived earlier. The latter equation takes into account interactions between the solid thick cylindrical wall and confined fluid, this body-body interaction potential has been primarily calculated using the Lennard-Jones (6-12) expression for the atom-atom pair potentials and expressed by hypergeometrical functions having good convergences. All numerical calculations shown here have been performed for the model graphite-argon system at 90 K. Finally, an analysis of the accuracy of the proposed method is considered.     

The development of the fundamental concepts of surface thermodynamics

Author’s results concerning the most fundamental problems of the thermodynamics of surface phenomena are reviewed. The generalized Laplace-Young-Kelvin equation, phase rules, and Gibbs adsorption equations are presented. Analogs of Konovalov’s laws are describes as applied to surface phenomena. The surface tension dualism, the Gibbs equation for adsorption on solid surfaces, and the phase equilibrium condition for a soluble nanoparticle are explained. The general mechanochemical approach, chemical affinity tensor, and the discovery of the mechanochemical dissolution effect are characterized. A novel approach to the monolayer state equation is formulated based on an excluded area. The problems of nucleation and the theory of surface separation are reperted.     

Metastable extension of the liquid-vapor phase equilibrium curve and surface tension

The method of molecular dynamics has been used to calculate the parameters of liquid-vapor phase equilibrium and the surface tension in a two-phase system of 4096 Lennard-Jones particles. Calculations have been made in a range from the triple point to near-critical temperature and also at temperatures below the triple point corresponding to the metastable equilibrium of a supercooled liquid and supersaturated vapor. To determine the surface tension, along with a mechanical approach a thermodynamic one has been used as well. The latter was based on calculation of the excess internal energy of an interfacial layer. It has been shown that in accuracy the thermodynamic approach is as good as the more sophisticated mechanical one. Low-temperature asymptotics of the phase-equilibrium curve and also of liquid and vapor spinodals have been considered in the Lennard-Jones and the van der Waals models. The behavior of the surface tension and the excess internal energy of an interfacial layer at T-->0 is discussed.     

Adhesion and wetting in an aqueous environment: theoretical assessment of sensitivity to the solid surface energy

The sensitivity of contact angles in an aqueous environment to the surface energy of the solid is discussed. It is demonstrated that this sensitivity is much higher in an aqueous than in a vapor environment. The Girifalco-Good equation, in combination with the Owens-Wendt equation, is used for the approximate demonstrations. It is shown that the transition from complete wetting to complete dewetting by the aqueous phase in a solid-liquid-liquid system occurs over a much narrower range of the surface energy of the solid than in a solid-liquid-vapor system. It is also demonstrated that the contact angle may be extremely sensitive to small variations in the relationship between surface tensions and the corresponding interfacial tension.     

Nanoparticle adsorption at liquid-vapor surfaces: influence of nanoparticle thermodynamics, wettability, and line tension

We developed a statistical mechanical theory that describes the adsorption of nanoparticles (NPs) at liquid-vapor surfaces. This theory accounts for the surface to bulk NP thermodynamic equilibrium, as well as the NP mechanical equilibrium, wettability, and line tension at liquid-vapor surfaces. The theory is tested by examining the adsorption of 5 nm diameter dodecanethiol-ligated gold NPs at the liquid-vapor surface of a homologous series of n-alkane solvents, from n-nonane to n-octadecane, where the NP wettability decreases with an increasing n-alkane chain length.     

On the applicability of Young–Laplace equation for nanoscale liquid drops

Debates continue on the applicability of the Young–Laplace equation for droplets, vapor bubbles and gas bubbles in nanoscale. It is more meaningful to find the error range of the Young–Laplace equation in nanoscale instead of making the judgement of its applicability. To do this, for seven liquid argon drops (containing 800, 1000, 1200, 1400, 1600, 1800, or 2000 particles, respectively) at T = 78 K we determined the radius of surface of tension R_s and the corresponding surface tension γ_s by molecular dynamics simulation based on the expressions of R_s and γ_s in terms of the pressure distribution for droplets. Compared with the two-phase pressure difference directly obtained by MD simulation, the results show that the absolute values of relative error of two-phase pressure difference given by the Young–Laplace equation are between 0.0008 and 0.027, and the surface tension of the argon droplet increases with increasing radius of surface of tension, which supports that the Tolman length of Lennard-Jones droplets is positive and that Lennard-Jones vapor bubbles is negative. Besides, the logic error in the deduction of the expressions of the radius and the surface tension of surface of tension, and in terms of the pressure distribution for liquid drops in a certain literature is corrected.     

Liquid-vapor and liquid-liquid interfaces in Ising fluids: an integral equation approach

The microscopic structure and thermodynamic properties of liquid-vapor and liquid-liquid interfaces in Ising fluids are studied using an integral equation approach. The calculations are performed in the absence and presence of an external magnetic field by solving the corresponding set of Lovett-Mou-Buff-Wertheim integrodifferential equations for the one-particle density distribution functions. The two-particle inhomogeneous direct correlation functions are consistently constructed by nonlinear interpolation between the bulk ones. The bulk correlation functions of the coexisting phases are obtained from the Ornstein-Zernike equations with a modified soft mean spherical approximation for the closure relation. As a result, the density and magnetization profiles at liquid-vapor and liquid-liquid interfaces as well as the surface tension and adsorption coefficients are evaluated in a wide temperature range including subcritical regions. The influence of an external magnetic field on the liquid-vapor interfaces is also considered.     

Dynamic surface tension,critical micelle concentration,and activity coefficients of aqueous solutions of nonyl phenol ethoxylates

Dynamic surface tensions, σ(t) for aqueous solutions of nonyl phenol ethoxylates (NPEOs) at the temperature 298.15 K were measured using a Lauda drop volume tensiometer. The non-ionic surfactants analyzed in this work were Tergitol NP-9, NP-35 and NP-40. By using the classical Ward and Torday equation, the diffusion coefficient for each bulk surfactant concentration was calculated. The equilibrium surface tension values were determined by extrapolating the dynamic surface tension to t → ∞ on the σ(t) vs. t^−1/2 curves. These values were used to determine the critical micelle concentrations (CMC) of the surfactant aqueous solutions as well as to calculate the infinite dilution activity coefficient of the surfactant, following a model that combines the Volmer surface equation of state and the Gibbs adsorption equation.     

Thermodynamic properties of binary mixtures containing aldehydes. II. Liquid-vapor equilibria in normal or branched alkanal + normal alkane mixtures: Analysis in terms of a quasichemical group-contribution model

Recent liquid-vapor equilibrium data for three alkanal + n-alkane mixtures are examined on the basis of the surface-interaction version of the quasichemical group-contribution theory used in Part I to correlate and predict excess enthalpies and excess Gibbs energies for such mixtures. The predictions prove to be accurate to better than 10%. Using the new data, revised interaction parameters are proposed for the estimation of liquid-vapor equilibrium for normal or branched alkanal + normal alkane mixtures.     

Constitutive modeling of contact angle hysteresis

We introduce a phase field model of wetting of surfaces by sessile drops. The theory uses a two-dimensional non-conserved phase field variable to parametrize the Gibbs free energy of the three-dimensional system. Contact line tension and contact angle hysteresis arise from the gradient term in the free energy and the kinetic coefficient respectively. A significant advantage of this approach is in the constitutive specification of hysteresis. The advancing and receding angles of a surface, the liquid-vapor interfacial energy and three-phase line tension are the only required constitutive inputs to the model. We first simulate hysteresis on a smooth chemically homogeneous surface using this theory. Next we show that it is possible to study heterogeneous surfaces whose component surfaces are themselves hysteretic. We use this theory to examine the wetting of a surface containing a circular heterogeneous island. The contact angle for this case is found to be determined solely by the material properties at the contact line in accord with recent experimental data.     

Adsorption hysteresis for a slit-like pore model

The Frenkel-Halsey-Hill equation is used to describe the adsorption branch of a hysteresis loop upon polylayer adsorption with an H3 loop according to IUPAC nomenclature. The equation for the desorption branch of a hysteresis loop is derived from a combined solution to the equation for the Gibbs potential change, given the adsorbent swelling and pore connectivity function, and the Laplace equation taken for the conditions of infinitely elongated meniscus. This equation is shown to connect the adsorbate relative pressure in a bulk phase for the desorption branch of a hysteresis loop with the key parameters of the adsorption system. The equation obtained was verified by a water adsorption isotherm on natural mineral schungite.     

A molecular dynamics simulation study of nanoscale liquid threads

The properties of liquid threads in angular direction are studied. On the basis of Gibbs theory, a thermodynamic model is constructed and a formula of surface of tension is derived. For the determination of surface tension, different methods are presented. We investigate seven different systems (the numbers of molecules N are 1600, 2240, 2880, 3360, 4000, 4800, and 5280) by molecular dynamics simulations. We analyze the surface tension obtained by various routes, then whether the classical Laplace equation is applicable in nanoscale can be determined. The results obtained by molecular dynamics simulations support that surface tension is dependent on the dividing surface curvature.     

Lateral capillary forces between solid bodies on liquid surface: a lattice Boltzmann study

When two solid bodies are placed on the surface of a dense liquid under gravitation, they deform the liquid surface to experience a lateral capillary force between themselves that can be attractive and repulsive, depending on the wettabilities and weights of the bodies. In the present study, the lateral capillary force between two square bodies at a liquid-vapor interface has been examined using numerical simulations based on a two-dimensional two-phase lattice Boltzmann (LB) method. The particular situations were simulated, where every body was vertically constrained and had the fixed triple points at its upper or lower corners. Here, the triple point indicates the place at which vapor, liquid, and solid phases meet. The interaction force between these two bodies was calculated as a function of the separation distance, the interfacial tension, and the gravitational acceleration. The simulation results agree well with the analytical expression of the lateral capillary interaction, indicating that our LB method can reproduce the interaction force between two bodies of various wettabilities at a liquid-vapor interface in mechanical equilibrium.