On the relation between virial coefficients and the close-packing of hard disks and hard spheres

The question of whether the known virial coefficients are enough to determine the packing fraction η(∞) at which the fluid equation of state of a hard-sphere fluid diverges is addressed. It is found that the information derived from the direct Pade? approximants to the compressibility factor constructed with the virial coefficients is inconclusive. An alternative approach is proposed which makes use of the same virial coefficients and of the equation of state in a form where the packing fraction is explicitly given as a function of the pressure. The results of this approach both for hard-disk and hard-sphere fluids, which can straightforwardly accommodate higher virial coefficients when available, lends support to the conjecture that η(∞) is equal to the maximum packing fraction corresponding to an ordered crystalline structure.     

A simplified hard-sphere equation of state meeting the high and low density limits

A simplified hard-sphere equation of state has been developed, which meets the correct limit for close-packed conditions. It is shown that the proposed equation of state for hard spheres can represent accurately the computer simulation compressibility factor data and virial coefficients over a wide density range. The comparison of the results of the calculations using this equation, the Carnahan-Starling equation, and the two equations proposed by Iglesias-Silva and Hall, shows that the equation proposed here represents the compressibility factor data and the virial coefficients with better accuracy.     

On the accurate direct computation of the isothermal compressibility for normal quantum simple fluids: application to quantum hard spheres

A systematic study of the direct computation of the isothermal compressibility of normal quantum fluids is presented by analyzing the solving of the Ornstein-Zernike integral (OZ2) equation for the pair correlations between the path-integral necklace centroids. A number of issues related to the accuracy that can be achieved via this sort of procedure have been addressed, paying particular attention to the finite-N effects and to the definition of significant error bars for the estimates of isothermal compressibilities. Extensive path-integral Monte Carlo computations for the quantum hard-sphere fluid (QHS) have been performed in the (N, V, T) ensemble under temperature and density conditions for which dispersion effects dominate the quantum behavior. These computations have served to obtain the centroid correlations, which have been processed further via the numerical solving of the OZ2 equation. To do so, Baxter-Dixon-Hutchinson's variational procedure, complemented with Baumketner-Hiwatari's grand-canonical corrections, has been used. The virial equation of state has also been obtained and several comparisons between different versions of the QHS equation of state have been made. The results show the reliability of the procedure based on isothermal compressibilities discussed herein, which can then be regarded as a useful and quick means of obtaining the equation of state for fluids under quantum conditions involving strong repulsive interactions.     

van der Waals-Tonks-type equations of state for hard-hypersphere fluids in four and five dimensions

Recently, we developed accurate van der Waals-Tonks-type equations of state for hard-disk and hard-sphere fluids by using the known virial coefficients. In this paper, we derive the van der Waals-Tonks-type equations of state. We further apply these equations of state to hard-hypersphere fluids in four and five dimensions. In the low-density fluid regime, these equations of state are in good agreement with the simulation results and existing equations of state.     

A cubic perturbed hard-chain equation of state

A cubic equation of state is developed on the basis of perturbation theory. The equation is an association of three segments: the hard-sphere, the hard-chain, and the attraction. The expression for each segment was invoked from approximations of computer simulations of rigorous molecular theories of fluids, but compromised to some extent accuracy and theory for simplicity. This model equation is shown to be potentially capable of describing the PVT behavior of real fluids. As limiting cases, the new equation is reduced to expressions for the hard-sphere and the hard-body fluids. It also represents square-well fluids when the hard-chain contribution is eliminated. The square-well equation was found satisfactory in conforming with the molecular simulation results for square-well fluids and their mixtures.     

Equation of state of a seven-dimensional hard-sphere fluid. Percus-Yevick theory and molecular-dynamics simulations

Following the work of Leutheusser [Physica A 127, 667 (1984)], the solution to the Percus-Yevick equation for a seven-dimensional hard-sphere fluid is explicitly found. This allows the derivation of the equation of state for the fluid taking both the virial and the compressibility routes. An analysis of the virial coefficients and the determination of the radius of convergence of the virial series are carried out. Molecular-dynamics simulations of the same system are also performed and a comparison between the simulation results for the compressibility factor and theoretical expressions for the same quantity is presented.     

A theoretical form of the Martin-Hou equation of state

A new equation of state is derived from the Barker-Henderson hard-sphere perturbation theory. It has the form similar to the Martin-Hou equation of state. The numerical values of the characteristic constants in the equation can be calculated by the method of Martin and Hou. The equation can be used to predict P-V-T properties accurately for fluids when the critical parameters (T_c, P_c and V_c) and one point on the vapor pressure cure are given. By using the functional relationships between the characteristic constants and the microscopic parameters, the molecular microscopic parameters of the substance can be obtained.     

Communication: inferring the equation of state of a metastable hard-sphere fluid from the equation of state of a hard-sphere mixture at high densities

A possible approximate route to obtain the equation of state of the monodisperse hard-sphere system in the metastable fluid region from the knowledge of the equation of state of a hard-sphere mixture at high densities is discussed. The proposal is illustrated by using recent Monte Carlo simulation data for the pressure of a binary mixture. It is further shown to exhibit high internal consistency.     

Magnetophoresis, sedimentation, and diffusion of particles in concentrated magnetic fluids

A dynamic mass transfer equation for describing magnetophoresis, sedimentation, and gradient diffusion of colloidal particles in concentrated magnetic fluids has been derived. This equation takes into account steric, magnetodipole, and hydrodynamic interparticle interactions. Steric interactions have been investigated using the Carnahan-Starling approximation for a hard-sphere system. In order to study the effective interparticle attraction, the free energy of the dipolar hard-sphere system is represented as a virial expansion with accuracy to the terms quadratic in particle concentration. The virial expansion gives an interpolation formula that fits well the results of computer simulation in a wide range of particle concentrations and interparticle interaction energies. The diffusion coefficient of colloidal particles is written with regard to steric, magnetodipole and hydrodynamic interactions. We thereby laid the foundation for the formulation of boundary-value problems and for calculation of concentration and magnetic fields in the devices (for example, magnetic fluid seals and acceleration sensors), which use a concentrated magnetic fluid as a working fluid. The Monte-Carlo methods and the analytical approach are employed to study the magnetic fluid stratification generated by the gravitational field in a cylinder of finite height. The coefficient of concentration stratification of the magnetic fluid is calculated in relation to the average concentration of particles and the dipolar coupling constant. It is shown that the effective particle attraction causes a many-fold increase in the concentration inhomogeneity of the fluid if the average volume fraction of particles does not exceed 30%. At high volume concentrations steric interactions play a crucial role.     

Semi-empirical insertion probability for hard-spheres and hard-sphere chains

In spite of its simplicity and a well-defined theoretical basis, the Flory–Guggenheim approach is conventionally regarded as inapplicable to off-lattice system since the insertion probability of the approach does not account for the excluded region, existing in the off-lattice system. In this work, we propose the insertion probability accounting for the excluded region of off-lattice fluids and derive a new version of equation of state (EOS) for hard-sphere chains basing on the Flory–Guggenheim approach. To investigate the behavior of the excluded regions, a Monte Carlo sampling was performed for hard disks and the various excluded regions were found to have different density dependence. On the basis of the simulation result, we formulated the insertion probability for hard-sphere and that of hard-sphere chain which accounts for the effect of chain-connectivity on the monomer insertion. The proposed insertion probability was found to correctly predict the simulation data for monomer and correctly correlate the simulation data for chain fluids. The resulting EOS was found to meet closed-packed limit and predict the simulation data of compressibility factor for monomer and chains with a reasonable degree of accuracy. When compared with other off-lattice based EOS, it shows a comparable or better result. For second virial coefficient of chain molecules, the model was found to reasonably predict the simulation data.     

Equation of state for hard-spheres and excess function calculations of binary liquid mixtures

Using our previously proposed matrix method, an equation of state for hard spheres is presented, which can reproduce the exact values of the first-eight virial coefficients. This equation meets both the low density and the close-packed limits and can predicts the first order fluid–solid phase transition of hard spheres. The results obtained show that the new equation of state can correlate the simulation data of compressibility factor up to high densities better than other equations of state.The new equation of state is extended to mixtures of hard spheres and excess functions of various binary liquid mixtures are calculated using the perturbation theory of Leonard–Henderson–Barker. The results are compared with existing theoretical and experimental data and with those calculated by other hard-sphere equations of state.It is seen that the results obtained by the new equation of state is quite satisfactory compared to other equations of state for the hard spheres and mixture of hard spheres.     

Inverse power potentials: virial coefficients and a general equation of state

Virial coefficients up to the seventh are calculated for pair potentials depending on inverse powers of separation, for inverse powers from 5 to 80. Unlike the limiting (infinite inverse power) hard-sphere potential, some virial coefficients for finite inverse power potentials are found to be negative. This makes resummation of the virial series for general inverse power potentials more difficult than that for hard spheres, and some alternative resummation methods are presented and compared. A general equation of state is proposed for fluids of particles interacting through inverse power pair potentials, for inverse powers greater than about 10. This includes the "molecular" inverse power of 12, for which the current results support and extend the results of previous studies.     

Model for the free-volume distributions of equilibrium fluids

We introduce and test via molecular simulation a simple model for predicting the manner in which interparticle interactions and thermodynamic conditions impact the single-particle free-volume distributions of equilibrium fluids. The model suggests a scaling relationship for the density-dependent behavior of the hard-sphere system. It also predicts how the second virial coefficients of fluids with short-range attractions affect their free-volume distributions.     

Excess energy and equation of state of fluids with hard-core potential models from a second-order Monte Carlo perturbation theory

Monte Carlo simulations in the NVT ensemble of the reference hard-sphere fluid have been performed to obtain the “exact” first- and second-order terms in the inverse temperature expansion of the free energy of fluids with hard-core potentials. The results have been used to obtain parametrizations of the free energy of fluids with Sutherland potentials with variable range as well as for a fluid with a hard-core Lennard–Jones potential. The results for the excess energy and the equation of state are compared with simulation data available in the literature for these fluids.     

Shape factors in equations of state. Part II. Repulsion phenomena in multicomponent chain fluids

An equation of state for the multicomponent fluid phase of nonattracting rigid particles of arbitrary shape is presented. The equation is a generalization of a previously presented equation of state for pure fluids of rigid particles; the approach describes the volumetric properties of a pure fluid in terms of a shape factor, zeta, which can be back calculated by scaling the volumetric properties of pure fluids to that of a hard sphere. The performance of the proposed equation is tested against mixtures of chain fluids immersed in a "monomeric" solvent of hard spheres of equal and different sizes. Extensive new Monte Carlo simulation data are presented for 19 binary mixtures of hard homonuclear tangent freely-jointed hard sphere chains (pearl-necklace) of various lengths (three to five segments), with spheres of several size ratios and at various compositions. The performance of the proposed equation is compared to the hard-sphere SAFT approach and found to be of comparable accuracy. The equation proposed is further tested for mixtures of spheres with spherocylinders. In all cases, the equation proved to be accurate and simple to use.     

Theory of transport in liquid metals: III. Calculation of shear viscosity coefficients of binary alloys

A theoretical procedure for the calculation of the shear viscosity coefficients of liquid binary alloys is proposed based on the simulation of an alloy by a hypothetical single-component hard-sphere liquid and treatment of the latter according to a corrected Enskog approach. The Enskog result is corrected according to a corresponding states principles developed by the authors in a previous work. The pair correlation function at contact necessary for the Enskog calculation is obtained from the solution of the Carnahan and Starling equation of state for single component hard-sphere fluids. The proposed theory is tested numerically by comparing its viscosity results for fifteen alloys with experimental data and with the corresponding results of the Eyring approach and the Tham and Gubbins theory. The latter theory is implemented numerically with pair-correlation values at contact satisfying exactly the Carnahan and Starling equation of state for binary hard-sphere fluids. These values result from analytic and closed-form expressions for the pair-correlation function developed presently through a correction of the corresponding Lebowitz expressions.     

Semiclassical statistical mechanics of two-dimensional hard-body fluids

The problem of calculating the thermodynamic properties of two-dimensional semiclassical hard-body fluids is studied. Explicit expressions are given for the first-order quantum corrections to the free energy, equation of state, and virial coefficients. The numerical results are calculated for the planar hard dumbbell fluid. Significant features are the increase in quantum corrections with increasing eta and increasing L*=L/sigma(0).     

Equation of state for hard-sphere fluid in restricted geometry

Many practical applications require the knowledge of the equation of state of fluids in restricted geometry. We study a hard-sphere fluid at equilibrium in a narrow cylindrical pore with hard walls for pore radii R<((square root 3)+2)/4 (in units of the hard sphere diameter). In this case each particle can interact only with its nearest neighbors, which makes possible the use of analytical methods to study the thermodynamics of the system. Using a transfer operator formalism and expanding in low- and high-pressure regions, we can obtain a simple analytical equation of state for almost all ranges of pressure. The results agree with Monte Carlo simulations. Additionally, it is shown that a convenient analytical representation can be chosen to accurately describe the equation of state within the error of the Monte Carlo simulation.     

Equation of state and liquid-vapor equilibrium of polarizable Stockmayer fluids

In this work we develop the concept of an effective potential to obtain the equation of state of polarizable Stockmayer (PSM) fluids. This potential consists of a Lennard-Jones function with appropriate energy and distance parameters that depend on the reduced dipolar moment μ(?) and polarizability α(?). The approach deals accurately with polarizable SM fluids with μ(?)≤2.0 and α(?)≤0.1. However, prediction of second virial coefficients is reliable up to μ(?)≤4.0. When the low-density sphericalized potential is used at moderate and large densities, the effect of the dipole-dipole attraction is overestimated in agreement with an effect previously found in the literature. This effect can be traced back to a frustration mechanism due to the interaction between three and more dipoles. We propose a model to account for this frustration effect and are able to reproduce the vapor-liquid equilibrium of polarizable SM fluids in agreement with simulated results from the literature. Molecular dynamics simulations were carried out to show that the effective SM fluid has a radial distribution function very close to that of the true SM system.     

Computational study of the melting-freezing transition in the quantum hard-sphere system for intermediate densities. I. Thermodynamic results

The points where the fluid-solid (face-centered-cubic) transition takes place in the quantum hard-sphere system, for reduced densities 0.85>rhoN*>0.5 (reduced de Broglie wavelengths lambdaB*相似文献