‘Exact’ integral equation theory and local formulation for excess thermodynamic properties of hard spheres

A hybrid hard sphere bridge function is proposed, which, in combination with the standard Ornstein–Zernike integral equation, can predict extremely accurately hard sphere compressibility, virial pressure, and correlation function. Second, a local formulation for determination of excess chemical potential is derived out, which, in combination with the present hybrid hard sphere bridge function and OZ integral equation, can predict the excess chemical potential also extremely accurately. The resultant excess entropy is in excellent agreement with that from the Carnahan–Starling equation of state. The present formalism performs excellently over the whole density range, i.e. from zero to freezing density, and is largely superior to a formalism available in the literature.     

Non-hard sphere thermodynamic perturbation theory

A non-hard sphere (HS) perturbation scheme, recently advanced by the present author, is elaborated for several technical matters, which are key mathematical details for implementation of the non-HS perturbation scheme in a coupling parameter expansion (CPE) thermodynamic perturbation framework. NVT-Monte Carlo simulation is carried out for a generalized Lennard-Jones (LJ) 2n-n potential to obtain routine thermodynamic quantities such as excess internal energy, pressure, excess chemical potential, excess Helmholtz free energy, and excess constant volume heat capacity. Then, these new simulation data, and available simulation data in literatures about a hard core attractive Yukawa fluid and a Sutherland fluid, are used to test the non-HS CPE 3rd-order thermodynamic perturbation theory (TPT) and give a comparison between the non-HS CPE 3rd-order TPT and other theoretical approaches. It is indicated that the non-HS CPE 3rd-order TPT is superior to other traditional TPT such as van der Waals/HS (vdW/HS), perturbation theory 2 (PT2)/HS, and vdW/Yukawa (vdW/Y) theory or analytical equation of state such as mean spherical approximation (MSA)-equation of state and is at least comparable to several currently the most accurate Ornstein-Zernike integral equation theories. It is discovered that three technical issues, i.e., opening up new bridge function approximation for the reference potential, choosing proper reference potential, and/or using proper thermodynamic route for calculation of f(ex-ref), chiefly decide the quality of the non-HS CPE TPT. Considering that the non-HS perturbation scheme applies for a wide variety of model fluids, and its implementation in the CPE thermodynamic perturbation framework is amenable to high-order truncation, the non-HS CPE 3rd-order or higher order TPT will be more promising once the above-mentioned three technological advances are established.     

A theoretical investigation on the honeycomb potential fluid

A local self-consistent Ornstein-Zernike (OZ) integral equation theory (IET) is proposed to provide a rapid route for obtaining thermodynamic and structural information for any thermodynamically stable or metastable state points in the bulk phase diagram without recourse to traditional thermodynamic integration, and extensive NVT-Monte Carlo simulations are performed on a recently proposed honeycomb potential in three dimensions to test the theory's reliability. The simulated quantities include radial distribution function (rdf) and excess internal energy, pressure, excess chemical potential, and excess Helmholtz free energy. It is demonstrated that (i) the theory reproduces the rdf very satisfactorily only if the bulk state does not enter deep into a two phases coexistence region; (ii) the excess internal energy is the only one of the four thermodynamic quantities investigated amenable to the most accurate prediction by the present theory, and the simulated pressure is somewhat overestimated by the theoretical calculations, but the deviation tends to vanish along with rising of the temperature; (iii) using the structural functions from the present local self-consistent OZ IET, a previously derived local expression, due to the present author, achieves even a higher accuracy in calculating for the excess chemical potential than the exact virial pressure formula for the pressure, and the resulting excess Helmholtz free energy is in surprisingly same with the simulation results due to offset of the errors. Based on the above observations, it is suggested that it may be a good procedure to integrate the theoretical excess internal energy along the isochors to get the excess Helmholtz free energy, which is then fitted to a polynomial to be used for calculation of all of other thermodynamic quantities in the framework of the OZ IET.     

Extending the simple weighted density approximation for a hard-sphere fluid to a Lennard-Jones fluid II. Application

A simple weighted density approximation (SWDA) was extended to nonuniform Lennard-Jones fluids by following the spirit of a partitioned density function theory [S. Zhou, Phys. Rev. E 68 (2003) 061201] and mapping the hard-core part onto an effective hard-sphere fluid whose higher order terms beyond the second order of the functional perturbation expansion are treated by the SWDA. The resultant DFT formalism performs well for Lennard-Jones fluids under the influence of diverse external fields. With the present DFT formalism, we investigate in detail the structure and adsorption properties of a low-density LJ gas in a spherical cavity with a wall consisting of hard-sphere or LJ particles. It was found that when the cavity wall exerts an attractive external potential on the LJ particles in the cavity, the excess adsorption decreases as the temperature increases, while when the cavity wall exerts a hard repulsive external potential on the LJ particles in the cavity, the excess adsorption increases as the temperature increases.     

Local solvent density augmentation around a solute in supercritical solvent bath: 1. A mechanism explanation and a new phenomenon

A recently proposed partitioned density functional (DF) approximation (Phys. Rev. E 2003, 68, 061201) and an adjustable parameter-free version of a Lagrangian theorem-based DF approximation (LTDFA: Phys. Lett. A 2003, 319, 279) are combined to propose a DF approximation for nonuniform Lennard-Jones (LJ) fluid. Predictions of the present DF approximation for local LJ solvent density inhomogeneity around a large LJ solute particle or hard core Yukawa particle are in good agreement with existing simulation data. An extensive investigation about the effect of solvent bath temperature, solvent-solute interaction range, solvent-solute interaction magnitude, and solute size on the local solvent density inhomogeneity is carried out with the present DF approximation. It is found that a plateau of solvent accumulation number as a function of solvent bath bulk density is due to a coupling between the solvent-solute interaction and solvent correlation whose mathematical expression is a convolution integral appearing in the density profile equation of the DF theory formalism. The coupling becomes stronger as the increasing of the whole solvent-solute interaction strength, solute size relative to solvent size, and the closeness to the critical density and temperature of the solvent bath. When the attractive solvent-solute interaction becomes large enough and the bulk state moves close enough to the critical temperature of the solvent bath, the maximum solvent accumulation number as a function of solvent bath bulk density appears near the solvent bath critical density; the appearance of this maximum is in contrast with a conclusion drawn by a previous investigation based on an inhomogeneous version of Ornstein-Zernike integral equation carried out only for a smaller parameter space than that in the present paper. Advantage of the DFT approach over the integral equation is discussed.     

Approximative "one particle" bridge function B(1)(r) for the theory of simple fluids

New properties for the one particle bridge function B(1)(r), which are necessary to the calculation of the excess chemical potential betamue), are derived for the hard sphere fluid. The method, which only requires the knowledge of the bridge function B(2)(r), is based on an investigation of the correlation function dependence on the Kirkwood charging parameter. In this framework, the unavoidable question of topological homotopy is addressed. As far as B(2)(r) is considered as exact, this work provides useful information on B(1)(r) in the well identified dynamical regimes of the hard sphere fluid. Signatures of the transitions between these regimes are identified on the trends of B(1)(r). This approach provides self-consistent results for betamue) that agree very well with simulation data.     

A consistent calculation of the chemical potential for dense simple fluids

A general method to calculate the excess chemical potential betamuex, that is based on the Kirkwood coupling parameter's dependence of the correlation functions, is presented. The expression for the one particle bridge function B(1)r is derived for simple fluids with spherical interactions. Only the knowledge of the bridge function B(2)r is required. The accuracy of our approach is illustrated for a dense hard sphere fluid. As far as B(2)r is considered as exact, B(1)r is found to be, at high densities, the normalized bridge function -B(2)rB(2)(r=0). This expression ensures a consistent calculation of the excess chemical potential by satisfying implicitly the Gibbs-Duhem constraint. Only the pressure-consistency condition is necessary to calculate the structural and thermodynamic properties of the fluid.     

Equation of state for expanded fluid mercury: variational theory with many-body interaction

A variational associating fluid theory is proposed to describe equations of state for expanded fluid mercury. The theory is based on the soft-sphere variational theory, incorporating an ab initio diatomic potential and an attractive many-body potential; the latter is evaluated with quantum chemical methods and expressed as a function of the local atomic coordination number and the nearest-neighbor distance. The resultant equation of state can reproduce the observed gas-liquid coexistence curve with good accuracy, without introducing phenomenological effective pair potentials. Various thermodynamic quantities such as pressure, isocloric thermal pressure coefficient, adiabatic sound velocity, and specific heat are calculated over a wide density-temperature range and compared with available experimental data.     

A consistent integral equation theory for hard spheres

The standard integral equation approach is used to extract the bridge function and other correlation functions of hard spheres fluid. To achieve this, we first use a recent consistent closure relation proposed by Bomont et al. [J. Chem. Phys. 119, 2188 (2003)] that has already proven to be accurate to describe the Lennard-Jones fluid properties. Second, we take advantage of the coherent scheme derived by Bomont [J. Chem. Phys. 119, 11484 (2003)] to calculate the excess chemical potential, the entropy and some relative transport properties. Very good agreement is obtained for structural quantities and thermodynamic properties as compared to exact data at densities ranging from 0.1 to 0.9.     

Density functional approach to thermodynamics of self-assembly

A local density functional approximation for predicting the surface crystallization of a thermodynamically small system under gravity is described and tested. Using the model of the classical soft-sphere fluid, the state parameters for such systems are identified. A generalized phase diagram based upon the scaling variables is obtained; systems with the same reduced-state parameters exhibit identical profiles of thermodynamic properties such as density, pressure, and intrinsic chemical potential, measured in the direction of the applied field. The point-thermodynamic approximation of Rowlinson and the local density approximation of the density functional formalism are found to be remarkably accurate. A configurational temperature is defined and shown to agree with the corresponding kinetic temperature for inhomogeneous systems at equilibrium. The structural profiles at the crystal-fluid interface are indicative of a mesolayer of lower density crystal, not seen in the field-free isobaric crystal-liquid interface.     

Structure of liquid solutions near the critical point

Specific features of the structure of the critical state of binary liquid solutions leading to an anomalous behavior of the Rayleigh line due to a dramatic increase in concentration and density fluctuations are considered. It is shown that an experimental treatment must deal with two fluctuation regions near the critical point of solvent vaporization. In the first region, one can achieve a sufficient degree of accuracy by using theories like selfconsistent field theory. In the second region, which is closer to the critical point than the first region, scaling theory of secondorder phase transitions may be applied. It is found that the anomalous behavior of the Rayleigh line associated with kinetic coefficients is determined by the equilibrium thermodynamic properties and by the radius of fluctuation correlation (rinc). A general theory is developed for calculating thermodynamic potentials, especially the chemical potential and its concentration derivative in the fluctuation region. The results of these calculations are compared with the experimental data briefly described in the paper. Translated fromZhumal Strukturnoi Khimii, Vol. 39, No. 4, pp. 655–668, July–August, 1998.     

How critical fluctuations influence adsorption properties of a van der Waals fluid onto a spherical colloidal particle

A recently proposed 3rd-order thermodynamic perturbation theory (TPT) is extended to its 5th-order version and non-uniform counterpart by supplementing with density functional theory (DFT) and a number of ansatzs for a bulk 2nd-order direct correlation function (DCF). Employment of the ansatzs DCF enables the resultant non-uniform formalism devoid of any adjustable parameter and free from numerically solving an Ornstein–Zernike integral equation theory. Density profiles calculated by the present non-uniform formalism for a hard core attractive Yukawa (HCAY) fluid near a spherical geometry are favorably compared with corresponding simulation data available in literature, and are more accurate than those based on a previous 3rd + 2nd-order perturbation DFT. The non-uniform 5th-order TPT is employed to investigate adsorption of the HCAY fluid onto a colloidal particle; it is disclosed that a depletion adsorption can be induced when the coexistence bulk fluid is situated in neighborhood of a critical point or near a bulk vapor–liquid coexistence gaseous phase or liquid phase density. A physical interpretation is given for such depletion adsorption and for its connection with parameters of the potential under consideration, which is ascribed to critical density fluctuations existing within a wide region of the bulk diagram. For a large spherical external potential inducing wetting transition, it is found that only round wetting transition is found instead of 1st-order pre-wetting transition in the case of a planar wall external potential, and the wetting transition temperature increases relative to that for the planar wall external potential. The present theoretical results for wetting transitions are supported by previous investigation based on thermodynamic considerations and a phenomenological Landau mean field theory, and are also in conformity with the present qualitative physical interpretation.     

用密度泛函理论研究Lennard-Jones 流体在狭缝中的相平衡

用改进的基础度量理论(modified fundamental measure theory, MFMT)和密度Taylor展开分别表达过剩自由能中的短程作用和色散作用. 流体分子与狭缝壁之间的相互作用以10-4-3势能函数表达. 由巨势最小原理确定Lennard-Jones (LJ)流体在狭缝中的密度分布和过剩吸附量, 所得结果与分子模拟数据吻合良好. 根据平衡时两相温度, 化学势及巨势相等, 计算了LJ流体在狭缝中的相平衡.     

Interatomic potential-based semiclassical theory for Lennard-Jones fluids

An interatomic potential based semiclassical theory is proposed to predict the concentration and potential profiles of a Lennard-Jones (LJ) fluid confined in a channel. The inputs to the semiclassical formulation are the LJ parameters of the fluid and the wall, the density of channel wall atoms, and the average concentration of the fluid inside the channel. Using the semiclassical formulation, fluid confinement in channel with widths ranging from 2sigma ff to 100sigma ff, where sigma ff is the fluid-fluid LJ distance parameter, is investigated. The concentration and potential predicted by the semiclassical formulation are found to be in good agreement with those from equilibrium molecular dynamics simulations. While atomistic simulations in large channels are computationally expensive, the proposed semiclassical formulation can rapidly and accurately predict the concentration and potential profiles. The proposed semiclassical theory is thus a robust and fast method to predict the interfacial and "bulk" fluid phenomena in channels with widths ranging from the macroscale down to the scale of a few atomic diameters.     

Nonideal gas solvation thermodynamics

General expressions are obtained for the thermodynamic properties of nonideal gases of arbitrary composition. These include the equation of state and solvation thermodynamic functions of mixtures with any number of components, expanded to first order in total number density. The results are expressed in terms of binary second virial coefficients which are in turn related to binary interaction potential energy functions. Solvation thermodynamic functions (at both constant pressure and constant volume) are decomposed into solute-solvent and solvent-reorganization contributions, thus indicating how each of the latter quantities may be experimentally measured. Moreover, the results are used to identify solute-induced changes in solvent thermodynamic functions, as well as nonideal (excess) contributions to chemical reaction equilibria and solvation heat capacities. The effects of a solute on the reorganization energy and chemical potential of a nonideal solvent are shown to differ qualitatively from both the corresponding ideal gas and high density fluid results.     

Thermodynamic quantities of adsorption described by Freundlich isotherm

Thermodynamic analysis of the Freundlich adsorption isotherm was performed. Equations describing the thermodynamic characteristics of adsorption as a function of the adsorbed quantity are presented, and relationships between the excess differential and average molar thermodynamic adsorption quantities were obtained. It was shown that the thermodynamic approach does not contradict the molecular statistical theory and provides more general equations. A model of adsorption described by the Freundlich isotherm is discussed. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 175–177, January, 2000.     

Thermodynamic quantities of adsorption described by Freundlich isotherm

Thermodynamic analysis of the Freundlich adsorption isotherm was performed. Equations describing the thermodynamic characteristics of adsorption as a function of the adsorbed quantity are presented, and relationships between the excess differential and average molar thermodynamic adsorption quantities were obtained. It was shown that the thermodynamic approach does not contradict the molecular statistical theory and provides more general equations. A model of adsorption described by the Freundlich isotherm is discussed. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 175–177, January, 2000.