Quantum corrections to the radial distribution function of a fluid

The perturbation theory developed for liquids is used to derive an expression for the first-order quantum correction to the radial distribution function of a fluid. The result is given in terms of grand canonical ensemble distribution functions for the classical fluid. The equations giving the thermodynamic functions in terms of the radial distribution function are discussed, and differences in the quantum and classical cases emphasized. An equation relating the two- and three-body classical distribution functions, derived recently by Singh and Ram using an indirect method, is shown to be simply related to the second equation of the BGY hierarchy.     

Integral equation theory of a one dimensional hard-ellipse fluid between hard walls

Density functional theory is used to study the structure of a one dimensional fluid model of hard-ellipse molecules with their axes freely rotating in a plane, confined between hard walls. A simple Hypernetted chain (HNC) approximation is used for the density functional of the fluid and the integral equation for the density is obtained from the grand potential. The only required input is the direct correlation function of the one dimensional hard-ellipse fluid. For this model, the pressure, sum rule and the density at the walls are obtained. The Percus Yevick (PY), for lower density, and HNC, for higher density, integral equations are also solved to obtain the direct correlation function of hard-ellipse model introduced here. We obtain the average density at the wall as well as the radial density profile. We compare these with Monte Carlo simulations of the same model and find reasonable agreement.     

Density functional theory for nonuniform polymer melts: Based on the universality of the free energy density functional

A density functional theory is proposed for nonuniform freely jointed tangential hard sphere polymer melts in which the bonding interaction is treated on the basis of the properties of the Dirac δ-function, thus avoiding the use of the single chain simulation in the theory. The excess free energy is treated by making use of the universality of the free energy density functional and the Verlet-modified (VM) bridge function. To proceed numerically, one of the input parameters, the second-order direct correlation function of a uniform polymer melt is obtained by solving numerically the Polymer-RISM integral equation with the Percus-Yevick (PY) closure. The predictions of the present theory for the site density distribution, the partition coefficient and the adsorption isotherm, near a hard wall or between two hard walls are compared with computer simulation results and with those of previous theories. Comparison indicates that the present approach is more accurate than the previous integral equation theory and the most accurate Monte Carlo density functional theories. The predicted oscillations of the medium-induced force between two hard walls immersed in polymer melts are consistent with the experimental results available in the literature. Received 18 April 2000     

Universal cubic equation of state and contact values of the radial distribution functions for multi-component additive hard-sphere mixtures

A universal cubic equation of state (UC EOS) is proposed based on a modification of the virial Percus-Yevick (PY) integral equation EOS for hard-sphere fluid. The UC EOS is extended to multi-component hard-sphere mixtures based on a modification of Lebowitz solution of PY equation for hard-sphere mixtures. And expressions of the radial distribution functions at contact (RDFC) are improved with the form as simple as the original one. The numerical results for the compressibility factor and RDFC are in good agreement with the simulation results. The average errors of the compressibility factor relative to MC data are 3.40%, 1.84% and 0.92% for CP3P, BMCSL equations and UC EOS, respectively. The UC EOS is a unique cubic one with satisfactory precision among many EOSs in the literature both for pure and mixture fluids of hard spheres.     

Nonexistence of Kirkwood's instability in a uniform fluid

Kirkwood's instability in the theory of fluid-solid transitions is proved to be impossible. Fluctuation of the one-particle distribution function in the first equation of the BGY hierarchy is investigated beyond Kunkin and Frisch's treatment. The second equation of the BGY hierarchy is utilized to eliminate the three-particle distribution function left in the Kunkin-Frisch result. The final expression for the first-order fluctuation of the one-particle distribution function under the presence of an external field is written in a form including only the pair correlation function and agrees identically with the one obtained from the direct expansion of the oneparticle distribution function in terms of the external field.     

广义Morse流体解析状态方程及对氮的应用

 根据Ross变分微扰理论以及硬球流体Percus-Yevick（PY） 径向分布函数表达式,建立了广义Morse势流体的解析状态方程。与模拟结果的比较一方面证实了广义Morse 势模型的合理性;另一方面表明了解析Ross变分微扰理论的精度相当或略好于非解析的Weeks-Chandler-Anderson （mWCA）理论,而优于复杂的优化超网络链积分方程理论（RHNC）。该解析状态方程被应用于拟合处于环境温度和压强小于1 GPa情形流体氮的实验数据,所得到的势能参数被用于预测高温高密度情形氮流体的压强,预测结果证实,该解析状态方程可以很好地适用于较宽的压强和温度范围。     

Structural properties of effective potential model by liquid state theories

This paper investigates the structural properties of a model fluid dictated by an effective inter-particle oscillatory potential by grand canonical ensemble Monte Carlo (GCEMC) simulation and classical liquid state theories.The chosen oscillatory potential incorporates basic interaction terms used in modeling of various complex fluids which is composed of mesoscopic particles dispersed in a solvent bath,the studied structural properties include radial distribution function in bulk and inhomogeneous density distribution profile due to influence of several external fields.The GCEMC results are employed to test the validity of two recently proposed theoretical approaches in the field of atomic fluids.One is an Ornstein-Zernike integral equation theory approach;the other is a third order + second order perturbation density functional theory.Satisfactory agreement between the GCEMC simulation and the pure theories fully indicates the ready adaptability of the atomic fluid theories to effective model potentials in complex fluids,and classifies the proposed theoretical approaches as convenient tools for the investigation of complex fluids under the single component macro-fluid approximation.     

Analytical expressions for the correlation function of a hard sphere dimer fluid

A closed form expression is given for the correlation function of a hard sphere dimer fluid. A set of integral equations is obtained from Wertheim's multidensity Ornstein-Zernike integral equation theory with Percus-Yevick approximation. Applying the Laplace transformation method to the integral equations and then solving the resulting equations algebraically, the Laplace transforms of the individual correlation functions are obtained. By the inverse Laplace transformation, the radial distribution function (RDF) is obtained in closed form out to 3D (D is the segment diameter). The analytical expression for the RDF of the hard dimer should be useful in developing the perturbation theory of dimer fluids.     

Classical multicomponent fluid structure near solid substrates: Born-Green-Yvon equation versus density-functional theory

We report a study of adsorption of binary mixtures of hard spheres of different sizes on a hard wall by using a version of density-functional theory, the Born-Green-Yvon (BGY) equation and Monte Carlo simulations. Following the BGY approach introduced by Fischer and Methfessel for single-component fluids, the proposed extension uses coarse-grained densities to approximate the contact values of pair distribution function of hard spheres. A procedure for evaluation of the coarse-grained densities, leading to an exact theory in one dimension, is proposed. The density-functional theory employed here, however, uses the Meister-Kroll and Groot approach. Comparisons of theoretical calculations with Monte Carlo simulations, as well as with previous theoretical predictions, have shown that density-functional theory reproduces the pseudo-experimental data accurately, even for extremely large size ratios of molecules of both species. The accuracy of the predictions of the BGY approach is less satisfactory, and for higher bulk fluid densities discrepancies with numerical simulations have been found.     

Chandler—Silbey—Ladanyi integral equation theory for semiflexible molecules

Extension of Chandler—Silbey—Ladanyi (CSL) integral equation theory for the fluid of semi-flexible site—site molecules is proposed. The Percus—Yevick type of the closure is used to describe the structural properties of the fluid consisting of semiflexible linear chain triatomic molecules. Results for the site—site intramolecular and intermolecular radial distribution functions (RDFs) are compared with the corresponding computer simulation results and results of self-consistent reference interaction site model (RISM) theory. In general both theories give reasonably good agreement with corresponding Monte Carlo simulation data. The exception is the RDF between the terminal sites, for which none of the theories is satisfactory. The present version of CSL theory appears to be slightly more accurate than the self-consistent RISM approach.     

Equilibrium theories of simple liquids

In this article, the theory of equilibrium properties of simple classical fluids is reviewed. The various relationships between the pair potential φ(r) and the pair correlation function g(r) are explored, from the usual integral equations and the perturbation theories to the generalized random phase approximation proposed recently. Particular attention is devoted to the extraction of the intermolecular forces from a given experimental data on the structure and thermodynamics of fluids. In particular, the propagation of errors in these calculations arising due to uncertainties in the input data is discussed. Finally, the recent use of BGY integral equation and the vacancy-cell model in the study of solid-liquid transition and melting is discussed.     

An efficient method for accurate evaluation of the site–site direct correlation functions of molecular fluids

The direct correlation function (DCF) plays an important role in liquid integral-equation theories and non-mean-field applications of the classical density functional theory (DFT). While for a simple fluid the DCF can easily be calculated from the radial distribution function via the Fourier transform and/or, for special cases, can be derived from analytical solutions of the Ornstein–Zernike equation, computation of the site–site DCFs of a molecular fluid is more challenging because of numerical issues associated with solving the matrix integral equations. This paper describes a new theoretical method for accurate evaluation of the site–site DCFs of molecular fluids by combination of molecular simulation and analytical asymptotic analysis. The computational procedure entails four steps: (1) molecular simulation is used to calculate the site–site total correlation functions (TCFs) in real space; (2) the reference-interaction-site model (RISM) is used to calculate the site–site DCFs in Fourier space at large wavenumbers; (3) asymptotic expressions are derived for the TCFs and DCFs in the limit of small wavenumbers; and (4) site–site DCFs over the entire range are obtained by interpolation of the asymptotic results. The numerical procedure has been illustrated by application to bulk SPC/E water. Accurate evaluation of the site–site DCFs for water lays a foundation for future applications of the DFT to aqueous systems with atomic details.     

Phase diagram of two-dimensional binary Yukawa mixtures

We have used Ramakrishnan–Yussouff (RY) density functional theory (DFT) to explore the topology of the phase diagram of two-component charge stabilised colloidal suspensions confined to a two-dimensional plane. The particles of the system interact via purely repulsive soft core Yukawa potential. Pair correlation functions (PCFs) used as input informations in DFT were calculated by solving both the hypernetted chain (HNC) and Percus–Yevick (PY) integral equation theories. To test the relative performance of the HNC and PY theories in the context of phase transitions, we have also studied the corresponding one-component systems. We found that RY DFT with HNC PCFs does not stabilise solid in both the one- and two-component cases, whereas the PY theory does. By considering the freezing into the substitutionally disordered triangular solid, we found that the temperature-composition phase diagrams of the binary mixture are narrow spindles whose thickness depends on the symmetry of the mixture components and the value of the screening constant of the Yukawa potential. Although the phase diagram obtained by RY DFT with structural inputs calculated by the PY theory is found to be shifted to higher temperature region in the temperature-composition plane, however, it captures qualitatively all the essential features of the phase diagram. Our results are in principle verifiable through computer simulations and experiments.     

Binary fluid mixture of hard ellipses: Integral equation and weighted density functional theory

We study a two-dimensional (2D) classical fluid mixture of hard convex shapes. The components of the mixture are two kinds of hard ellipses with different aspect ratios. Two different approaches are used to calculate the direct, pair and total correlation functions of this fluid and results are compared. We first use a formalism based on the weighted density functional theory (WDFT), introduced by Chamoux and Perera [Phys. Rev. E 58 (1998) 1933]. Second, in general the Percus-Yevick (PY) and the hypernetted chain (HNC) integral equations are solved numerically for the 2D fluid mixtures of hard noncircular particles. Explicit results are obtained for the fluid mixtures of hard ellipses and comparisons are made by the two approaches. Also, the results are compared with the recent Monte Carlo simulation for the one-component fluids of hard ellipses. Finally we obtained the equation of state of hard ellipses for the aspect ratio sufficiently close to 1 and compared our results with the simulations of the fluid mixtures of hard disks.     

Analytic solution of Percus-Yevick equation for fluid of hard spheres

A new method of analytic solution of the Percus-Yevick equation for the radial distribution functiong(r) of hard-sphere fluid is proposed. The original non-linear integral equation is reduced to non-homogeneous linear integral equation of Volterra's type of the second order. The kernel of this new equation has a polynomial form which allows to find analytic expression forg(r) itself without using the Laplace transformation. In addition, the first three moments of the total correlation function can be found.     

Comportement asymptotique de la fonction de distribution radiale dans les liquides

We study the influence of pair potential on the asymptotic behaviour of the radial distribution function in liquids. We use a perturbation expansion from a model of rigid spherical core inKirkwood's integral equation. It is pointed out that in metals as well as in insulators the asymptotic behavior is connected to that of the pair potential. However, the results are numerically different from those obtained from thePercus-Yevick approximation.     

Variants of the optimised random phase approximation: from one direct correlation function to many different radial distribution functions

It is shown that there are an infinite number of radial distribution functions (RDFs), corresponding to only one direct correlation function (DCF) of the optimised random phase approximation (ORPA). This observation in the thermodynamic perturbation theory is in sharp contrast to that of integral equation theories in which they uniquely correspond. By devising a new method we will be able to introduce various perturbation theories of simple liquids all coming from one DCF. Among all, we will only present analytically seven variants of the ORPA in the thermodynamic perturbation theory of liquids. The DCF of hard-core potential for all variants is assumed to be the same as the ORPA. However, interestingly enough the resulted expressions for the Helmholtz free energies and the RDF are obtained very differently. Furthermore, the resulted thermodynamic properties come out somehow the same, whereas the structural functions of some variants are found to behave much better than the standard ORPA.     

Theoretical Investigation of Uniform and Non-uniform Penetrable Sphere Fluid

A bridge function approximation is proposed for a single-component fluid consisting of penetrable sphere interacting via a potential that remains finite and constant for center-center distance smaller than the particle diameter and is zero otherwise. The radial distribution function from the Ornstein-Zernike integral equation combined with the present bridge function approximation is in satisfactory agreement with the corresponding simulation data for all of the investigated state points. The presently calculated excess Helmholtz free energy respectively based on virial route and compressibility route is highly self-consistent, and is in very good agreement with simulational results for the case of low temperatures. The present bridge function approximation, combined with the bridge density functional approximation,can reproduce very accurately density profiles of the penetrable sphere fluid confined in a hard spherical cavity for all the cases where simulational results are available.     

Theoretical Investigation of Uniform and Non-uniform Penetrable Sphere Fluid

A bridge function approximation is proposed for a single-component fluid consisting of penetrable sphere interacting via a potential that remains finite and constant for center-center distance smaller than the particle diameter and is zero otherwise. The radial distribution function from the Ornstein-Zernike integral equation combined with the present bridge function approximation is in satisfactory agreement with the corresponding simulation data for all of the investigated state points. The presently calculated excess Helmholtz free energy respectively based on virial route and compressibility route is highly self-consistent, and is in very good agreement with simulational results for the case of low temperatures. The present bridge function approximation, combined with the bridge density functional approximation, can reproduce very accurately density profiles of the penetrable sphere fluid confined in a hard spherical cavity for all the cases where simulational results are available.     

An equation of state contribution for dipolar and quadrupolar square-well fluids

A dipolar–quadrupolar contribution to the residual Helmholtz energy for a polar square well (a square well plus either a point dipole or a point quadrupole) fluid is developed based on the Padé approximation. Taking the square well system as reference, the contribution is formulated using an expansion for radial distribution function of the reference system. In addition to square well potential parameters the contribution depends only on dipole and quadrupole moments. This term is added as perturbation to a generalized equation of state for square well fluids. The results are then compared with the available simulation data in the literature. With the new equation obtained, it was possible to predict liquid–vapour equilibrium properties and critical properties of polar square well fluids more accurately than with available perturbation theories for multipolar square well systems. Application of the equation of state to a real dipolar (water) and a real quadrupolar (carbon dioxide) fluid indicated that the polar contribution greatly improved the predictions of saturation properties. Accurate prediction of critical properties for polar square well fluids remains as a challenge. This work can be useful in the development of better equations of state.