A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

Tuning density profiles and mobility of inhomogeneous fluids

Density profiles are the most common measure of inhomogeneous structure in confined fluids, but their connection to transport coefficients is poorly understood. We explore via simulation how tuning particle-wall interactions to flatten or enhance the particle layering of a model confined fluid impacts its self-diffusivity, viscosity, and entropy. Interestingly, interactions that eliminate particle layering significantly reduce confined fluid mobility, whereas those that enhance layering can have the opposite effect. Excess entropy helps to understand and predict these trends.     

Surface phase transitions of multiple-site associating fluids

Surface phase transitions of Lennard–Jones (LJ) based two- and four-site associating fluids have been studied for various associating strengths using grand-canonical transition matrix Monte Carlo simulations. Our results suggest that, in the case of a smooth surface, represented by a LJ 9-3-type potential, multiple-site associating fluids display a prewetting transition within a certain temperature range. However, the range of the prewetting transition decreases with increasing associating strength and increasing number of sites on the fluid molecules. With the addition of associating sites on the surface, a quasi-2D vapor–liquid transition may appear, which is observed at a higher surface site density for weaker associating fluids. The prewetting transition at lower associating strength is found to shift towards the quasi-2D vapor–liquid transition with increasing surface site density. However, for highly associating fluids, the prewetting transition is still intact, but shifts slightly towards the lower temperature range. Adsorption isotherms, chemical potentials and density profiles are used to characterize surface phase transitions.     

Continuum percolation of the four-bonding-site associating fluids

Wertheim’s integral equation theory for associating fluids is reformulated for the study of the connectedness properties of associating hard spheres with four bonding sites. The association interaction is described as a square-well saturable attraction between these sites. The connectedness version of the Ornstein-Zernike (OZ) integral equation is supplemented by the PY-like closure relation and solved analytically within an ideal network approximation in which the network is represented as resulting from the crossing of ideal polymer chains. The pair connectedness functions and the mean cluster size are calculated and discussed. The condition for the percolation transition and the analytical form of the percolation threshold are derived. The connection of the percolation with the gas-liquid phase transition is discussed.     

Theory of the liquid-vapour interface of a binary mixture of Lennard-Jones fluids

We present results of calculations of the equilibrium surface tension and density profiles for the liquid-vapour interface of a binary mixture of Lennard-Jones 12-6 fluids. The calculations are based on a density-functional theory for the Helmholtz free energy of the inhomogeneous mixture. This is a ‘microscopic’ generalization of the van der Waals-Cahn-Hilliard theory for the interface of a binary mixture.

Our calculations cover the full range of liquid-vapour coexistence and the whole range of concentration. We find a correlation between the excess surface tension of the mixture and the surface segregation (adsorption) of the species with the lower surface tension. The ways in which segregation and excess surface tension depend on the Lennard-Jones parameters of the pure components are briefly discussed. Our results for the excess surface tension of mixtures of Ar and N₂ and Ar and CH₄ are compared with experiment; the agreement is reasonable.     

Generalized density-functional theory: extended weighted density approaches

A third-order density-functional theory is introduced by an approach that may be used to find density-functional theories to any higher-order accuracy provided only that known homogeneous state correlation functions are utilized as input. It is constructed from the required knowledge of a single weight function at each order. By way of application results are presented for the melting of classical hard spheres using functionals accurate to a third-order functional Taylor series in the homogeneous limit. Within the framework of the modified weighted density approximation, there is a uniform improvement in the solid phase-free energies, pressures and melting parameters, and further improvement also results when these functionals are optimized in a way that utilizes the close packing limit. The sensitivity of the results to existing and proposed models of the third-order direct correlation function is discussed.     

The vicinity of an equilibrium three-phase contact line using density-functional theory: density profiles normal to the fluid interface

The paper by Nold et al. [Phys. Fluids 26 (7), 072001 (2014)] examined density profiles and the micro-scale structure of an equilibrium three-phase (liquid–vapour–solid) contact line in the immediate vicinity of the wall using elements from the statistical mechanics of classical fluids, namely density-functional theory. The present research note, building on the above work, further contributes to our understanding of the nanoscale structure of a contact line by quantifying the strong dependence of the liquid–vapour density profile on the normal distance to the interface, when compared to the dependence on the vertical distance to the substrate. A recent study by Benet et al. [J. Phys. Chem. C 118 (38), 22079 (2014)] has shown that this could explain the emergence of a film-height-dependent surface tension close to the wall, with implications for the Frumkin–Derjaguin theory.     

Thermodynamic and transport properties of simple fluids using lattice sums: bulk phases and liquid-vapour interface

Molecular dynamics simulations in the canonical ensemble have been performed to obtain the thermodynamic and transport properties of the Lennard-Jones fluid. The dispersion interactions were calculated using lattice sums. This method makes it possible to simulate the full potential avoiding the inclusion of the long range corrections (LRC) during or at the end of simulations. In the calculation of dynamic properties in bulk phases and thermodynamic quantities of inhomogeneous systems where the interface is physically present, in general the LRC cannot easily be included. By using the lattice sums method, the results are independent of the truncation of the potential. In the liquid-vapour interface simulations it is not necessary to make any pre-judgments about the form of the LRC formula to calculate coexisting properties such as the surface tension. The lattice sums method has been applied to evaluate how well the full interaction can be calculated in the liquid phase and in the liquid-vapour interface. In the liquid phase the pressure, configurational energy, diffusion coefficient and shear viscosity were obtained. The results of the thermodynamic properties are compared with those obtained using the spherically truncated and shifted (STS) potential with the LRC added at the end of simulations, and excellent agreement is found. The transport properties are calculated on different system sizes for a state near the triple point. The diffusion coefficient using the lattice sums method increases with the number of molecules, and the results are higher than those of the STS model truncated at 2.5σ (STS2.5). The shear viscosity does not show any system size dependence for systems with more than 256 molecules, and the lattice sums results are essentially the same as those for the STS2.5. In the liquid-vapour equilibria the coexisting densities and vapour pressures for the full potential agree well with those obtained using the Gibbs ensemble and the NPT + test particle methods. The surface tension using lattice sums and truncation of forces at 2.5σ agrees well with STS results using large system sizes and cutoff distances.     

Classical charged fluids at equilibrium near an interface: Exact analytical density profiles and surface tension

The structure of equilibrium density profiles in an electrolyte in the vicinity of an interface with an insulating or conductive medium is of crucial importance in chemical physics and colloidal science. The Coulomb interaction is responsible for screening effects, and in dilute solutions the latter effects give rise to universal leading corrections to nonideality, which distinguish electrolyte from nonelectrolyte solutions. An example is provided by the excess surface tension for an air-water interface, which is determined by the excess particle density, and which was first calculated by Onsager and Samaras. Because of the discrepancy between the dielectric constants on both sides of the interface, every charge in the electrolyte interacts with an electrostatic image, and the Boltzmann factor associated with the corresponding self-energy has an essential singularity over the length scalel from the wall. Besides Coulomb interactions, short-range repulsions must be taken into account in order to prevent the collapse between charges with opposite signs or between each charge and its image when the solvent dielectric constant is lower than that of the continuous medium on the other side of the interface. For a dilute and weaklycoupled electrolyte,l is negligible with respect to the bulk Debye screening length ξ_D. In the framework of the grand-canonical ensemble, systematic partial resummations in Mayer diagrammatics allow one to exhibit that, in this regime, the exact density profiles at leading order are the same as if they were calculated in a partially-linearized mean-field theory, where the screened pair interaction obeys an inhomogeneous Debye equation. In the latter equation the effective screening length depends on the distancex from the interface: it varies very fast over the lengthl and tends to its bulk value over a few ξDs. The equation can be solved iteratively at any distancex, and the exact density profiles are calculated analytically up to first order in the coupling parameter l/ξ_D. They show the interplay between three effects: (1) the geometric repulsion from the interface associated with the deformation of screening clouds, (2) the polarization effects described by the images on the other side of the interface, (3) the interaction between each charge and the potential drop created by the electric layer which appears as soon as the fluid has not a charge-symmetric composition. Moreover, the expressions allow us to go beyond Onsager-Samaras theory: the surface tension is calculated for charge-asymmetric electrolytes and for any value of the ratio between the dielectric constants on both sides of the interface. Similar diagrammatic techniques also allow one to investigate the charge renormalization in the dipolar effective pair interaction along the interface with an insulating medium.     

A self-consistent density-functional approach to the structure of electric double layer: charge-asymmetric electrolytes

A self-consistent density-functional approach has been employed to study the structure of an electric double layer formed from a charge-asymmetric (2:l) electrolyte within the restricted primitive model which corresponds to charged hard sphere ions and a continuum solvent. The particle correlation due to hard-core exclusions is evaluated by making use of the universality of the density functionals and the correlation function of the uniform hard sphere fluid obtained through the integral equation theory with an accurate closure relation whereas mean spherical approximation is employed for the electrical contribution. Numerical results on the diffuse layer potential drop, ionic density profile, and the mean electrostatic potential near the electrode surface at several surface charge densities are found to be in quantitative agreement with the available simulation data.     

Local elastic properties of nano-confined fluids: A density functional study

The understanding of mechanical properties of confined fluids is essential for modeling and manipulating of nano-scaled systems. Unlike the uniform phase, the confined fluids usually display different features in structure and related properties. Due to the presence of the confining geometry, the density profile and many physical and chemical properties may be position-dependent. The aim of our research is to derive an expression for the local elastic property by using the classical elastic theory. Both the bulk and shear moduli are expressed as functional of density of particle. The theoretical result derived is applied to the Lennard-Jones fluids confined in nano-cavity. Comparison of our numerical result and the simulation result is made and qualitative agreement is observed. Further, influence of bulk density, temperature and external potential on moduli is calculated and the physical mechanism is analyzed. Relationship between contact modulus and the interfacial tension is also calculated. Their opposite trend with temperature is observed.     

Structure and phase transitions in a network-forming associating Lennard-Jones fluid in a slit-like pore: a density functional approach

A study is reported of adsorption of an associating Lennard-Jones fluid with four associative sites per molecule in a slit-like pore. The density distribution of particles in the pore and thermodynamics properties are evaluated by using a density functional method. It is found that at low temperatures the fluid exhibits a set of layering transitions, followed by capillary condensation. Transitions are localized by analysing the grand canonical potential. The density profiles of particles and the distribution of unbound and differently bonded particles demonstrate changes in the structure of the fluid in the pore along the phase coexistence. The critical temperature is lower for a confined fluid, compared with the bulk counterpart. However, an increase in the energy of association increases the critical temperature. The envelope of the capillary condensation is narrower than the bulk liquid-vapour phase diagram. The dependence of the solvation force on the energy of association and on the bulk density is discussed.     

Domain growth in ternary fluids: A level set approach

We analyze phase separation in ternary systems in the asymptotic hydrodynamic regime when the volume fractions and concentrations are constant. The multiphase Navier-Stokes equations are solved using a level set method. A new projection method was developed to treat multiple junctions for systems with more than two phases. It is found that surface tension ratios can alter the growth mechanism of a minority phase in the presence of two majority phases. When the minority phase wets the interface of the majority phases the domain growth rate of all three phases is initially similar to that of a symmetric binary fluid but slows down at later times.     

Bridging dielectric fluids by light: A ray optics approach

Rayleigh-Plateau instability is known to impose a stability limit for the length of a liquid bridge in weightless conditions. This fundamental limit may be exceeded by using a light field to form and stabilize dielectric fluid bridges (A. Casner, J.P. Delville, Europhys. Lett. 65, 337 (2004)). Using both new experimental data as well as a new theoretical approach, we show that both the size and the stability of such light-sustained dielectric bridge can be qualitatively explained. We present a ray optics model that encompasses the competition between surface tension effects and optical radiation pressure arising from total internal reflection inside the bridge. A critical power below which a liquid bridge can no longer be sustained by light is predicted and confirmed experimentally. The observed power dependence of the bridge diameter also agrees with the proposed stabilization mechanism.