Density expansion of the magnetic susceptibility

The pressure of a system may be expanded as a power series in the density, whose coefficients are the virial coefficients. In this paper, the magnetic susceptibility of a $\text{spin}\text{-}\text{1}\text{2}$ fermion system is also expanded in powers of density. This process explicitly separates the temperature and density dependence of the magnetic susceptibility. The coefficients of this series are shown to be related to certain virial coefficients. The first (previously established) and second corrections to Curie's law are explicitly expressed in terms of second and third virial coefficients. These corrections to Curie's law are small for temperatures above 4K, but become important below that temperature. The first correction has been previously measured. Given a set of second and third virial coefficients, the importance of the second correction can be calculated immediately at any density of interest.     

Integral-equation theories and Mayer-sampling Monte Carlo: a tandem approach for computing virial coefficients of simple fluids

Mayer-sampling Monte Carlo (MSMC) has enabled computation of higher-order virial coefficients than previously possible for a variety of potential models, but it is not required for computation of the entire virial coefficient for models that are spherically symmetric: approximations that result from the hypernetted-chain (HNC) or Percus–Yevick (PY) integral-equation theories in conjunction with the compressibility equation (c) or virial equation (v) can be computed quickly by fast Fourier transforms. For the fourth and fifth virial coefficients of the Lennard–Jones potential (with parameters σ and ε), we demonstrate that the corrections to each of the four approximations (HNC(c), HNC(v), PY(c), and PY(v)) are faster to compute to a desired precision by MSMC than the full coefficient itself, with the exception of the PY(v) correction at fifth order, and that the optimal decomposition with regard to precision can be identified using a fraction of the steps required to obtain precise virial coefficients. At reduced temperatures kT/ε greater than 4, the PY(c) correction is fastest to compute by MSMC at both fourth and fifth orders. For lower temperatures, the HNC(v) decomposition is most efficient at fourth order, while the HNC(c) decomposition is most efficient at fifth order. These results are specific to the Lennard–Jones potential, but the method for determining the optimal decomposition is applicable to any spherically symmetric potential.     

The influence of intermolecular interactions on the Kerr effect in gases

Expressions are derived for the second and third Kerr virial coefficients, B _K and C _K, of spherical top molecules in terms of irreducible cluster polarizabilities, and values are calculated using the dipole-induced dipole model for argon, krypton, xenon, methane, tetrafluoromethane, neopentane and sulphur hexafluoride. For mixtures of rare gases it is shown that the collision-induced dipole moment makes a negligible contribution to B _K. The effect of the choice of intermolecular potential function on the calculated second Kerr virial coefficients is also demonstrated. It is found that the predominant contributions to C _K arise from the pair polarizability, and that the triplet polarizability is only of minor importance.     

Quantum virial coefficients of molecular nitrogen

We report virial coefficients up to third order in density for molecular nitrogen, investigating 103 temperatures in the range (15 K, 3000 K). All calculations are based on an ab initio-based potential taken from the literature. Path-integral Monte Carlo (PIMC) is applied to account for nuclear quantum effects, and these results are compared to a more approximate but faster semiclassical treatment. Additionally, we examine a PIMC approach that employs semiclassical beads for the path-integral images, but find that it offers marginal advantage. A recently developed orientation sampling algorithm is used in conjunction with Mayer sampling to compute precise virial coefficients. We find that, within the precision of our calculations of the second-order coefficient (B₂), semiclassical methods are adequate for temperatures greater than 250 K, and are needed to correct classical behaviour for temperatures as high as 800 K. For the third-order coefficient (B₃), the semiclassical methods are adequate above 150 K, and are required up to the highest temperature examined (3000 K) in order to correct the classical treatment within the precision of the calculations. However, three-body contributions to the potential are much more significant than nuclear quantum effects for the evaluation of B₃.     

Virial coefficients expressed by heat kernel coefficients

In this paper, we generally expressed the virial expansion of ideal quantum gases by the heat kernel coefficients for the corresponding Laplace type operator. As examples, we give the virial coefficients for quantum gases in d-dimensional confined space and spheres, respectively. Our results show that, the relative correction from the boundary to the second virial coefficient is independent of the dimension and it always enhances the quantum exchange interaction. In d-dimensional spheres, however, the influence of the curvature enhances the quantum exchange interaction in two dimensions, but weakens it in higher dimensions ($d>3$).     

On the second-order virial coefficient and the relationships between quantum and classical statistics and scattering

We clarify some aspects of the relationship between quantum statistical mechanics and scattering theory, which show up in their simplest form in the behaviour of the second-order virial coefficient b₂. For this purpose we derive a new representation for b₂. The relationship with classical statistical mechanics is also illuminated by a recently obtained formula for the difference between the classical and quantum second order virial coefficients which allows the determination of the leading quantum correction to b₂. The several approaches represent alternative ways of implementing the cancellation of divergences in b₂.     

Atomic interactions in neon and helium

The first two quantum corrections to the second virial coefficients of the Smith-Thakkar potential are calculated. Parameters for neon and helium, gases in which quantum effects are important, are then determined by fitting to semiempirical dispersion coefficients and experimental second virial coefficients. Viscosity coefficients for both gases and vibrational energy level spacings for the neon dimer are calculated as independent tests of the potentials. Overall agreement with experiment is excellent for neon and moderate for helium. The previously determined parameters for argon are found to be only very slightly perturbed by the inclusion of quantum corrections in the calculated second virial coefficients.     

Fluids of general hard triatomic molecules

The second, third, and fourth virial coefficients, B_i , of a fluid of general symmetric hard triatomic molecules (fused hard spheres) have been calculated both numerically and theoretically for a variety of potential parameters. It has been found that: (i) for B₂ a valency angle ω_c exists such that for ω>w_c, B ₂ is independent of ω, (ii) B ₃ is very flat for ω>w_c, and (iii) B ₄ exhibits a maximum at ω～π/2. Theoretical calculations employing an assigned convex body fit very well the second and fairly well the third virial coefficients, but fail for the fourth except in the case of a linear molecule.     

Equilibrium properties of hard non-sphere fluids

Derivation of the thermodynamic properties of fluids of hard non-spherical molecules of arbitrary symmetry is based on the decoupling approximation. Theoretical expressions are given and calculations made for the equation of state and virial coefficients for hard ellipsoids. These results are compared with Monte Carlo values and show fair agreement in all cases. The theoretical predictions for the equation of state for binary mixtures are compared with the Monte Carlo results for hard spheres and hard prolate spherocylinders. Theoretical expressions for the first order quantum correction to the free energy, pressure and virial coefficients are also given. The quantum effects increase with increase of density and with increase of anisotropy parameter.     

The Quantum Excess Free Energy for Two Component Plasma

The aim of this paper is to calculate the quantum excess free energy until the third virial coefficient for the two component plasma (TCP). We consider only the thermal equilibrium plasma in the case of nλ³_ab ≪ 1 where λ²_ab = (ħ²/(m_abKT)) is the thermal De Broglie wave‐length. The second and third virial coefficients are represented in forms of a convergent series expansion in terms of the interaction Parameter ξ²_ab = ((e_ae_b)/(λ_abKT)). We obtained the excess free energy until the third virial coefficient to the order ξ⁶_ab. The effects of symmetry are taken into account, also the exchange effects were taken in account. We compared our results with others. Upon comparison with Ebeling et al. one observes that the considerable contribution of the third virial coefficient in the region of high temperatures (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Statistical mechanics of linear molecules II: Density expansions,intermolecular potentials and virial coefficients

In the Percus-Yevick and convolution-hypernetted-chain equations obtained in the previous paper, a density expansion for the correlation functions is introduced. To first order in the density, the so-obtained equations are identical and exact. By solving these, the pair correlation functions for linear molecules are obtained explicitly to first order in the density and for arbitrary order in the potential perturbation expansion. From these, the second and third virial coefficients can be extracted for all orders. A generalized charge is defined and used to give generalized multipole expansions for the intermolecular potential. Explicit expressions for this potential model are given up to fourth order. It is shown how the correlation functions, and second and third virial coefficients can be obtained to fourth order for any intermolecular potential with the same perturbation structure.     

Computation of high-order virial coefficients in high-dimensional hard-sphere fluids by Mayer sampling

The Mayer sampling method was used to compute the virial coefficients of high-dimensional hard-sphere fluids. The first 64 virial coefficients for dimensions 12 < D ? 100 were obtained to high precision, and several lower dimensional virial coefficients were computed. The radii of convergence of the virial series in 13, 15, 17 and 19 dimensions agreed well with the analytical results from the Percus–Yevick closure.     

Ab initio calculation of the refractivity and hyperpolarizability second virial coefficients of neon gas

Second virial coefficients for the density dependence of a number of electric properties are calculated for neon gas. Employing an accurate CCSD(T) potential for the Ne² van der Waals dimer and interaction-induced electric dipole polarizabilities and hyperpolarizabilities obtained from CCSD response theory, we evaluated the dielectric, refractivity, Kerr and ESHG second virial coefficients using both a semiclassical and a quantum statistical approach. The results cover a wide range of temperatures and are expected to be more reliable than the available experimental and empirical data. Quantum effects are found to be important only for temperatures below 100 K. The frequency-dependence of the refractivity virial coefficient is found to be small, but not negligible. For frequencies in the visible region it accounts for a few percent of the final results. For the ESHG virial coefficient of neon, frequency dependence is found to be very important, accounting for 20–25% of the second virial coefficient at the typical frequencies employed in experiments.     

Quantum corrections to thermodynamics of polar hard sphere fluids

The quantum corrections to the thermodynamic properties of polar hard sphere fluids and fluid mixtures are estimated taking into account the influence of dipole and quadrupole moments. Expressions are given for the second virial coefficient, free energy and pressure and results are given for different values ofμ* andϑ*. The first order quantum correction arises due to the translational contribution only. The quantum effect increases with density,μ* andϑ*. Numerical results are also estimated for binary mixtures of (i) hard spheres and dipole hard spheres and (ii) hard spheres and quadrupole hard spheres. The ‘excess’ free energy for dipole hard sphere binary mixture is also reported. It is found that the ‘excess’ quantum effect depends on the concentration and the particle diameter ratio and increases with increase ofμ* andϑ*.     

Two Hard Spheres in a Spherical Pore: Exact Analytic Results in Two and Three Dimensions

The partition function and the one- and two-body distribution functions are evaluated for two hard spheres with different sizes constrained into a spherical pore. The equivalent problem for hard disks is addressed too. We establish a relation valid for any dimension between these partition functions, second virial coefficient for inhomogeneous systems in a spherical pore, and third virial coefficients for polydisperse hard spheres mixtures. Using the established relation we were able to evaluate the cluster integral b ₂(V) related with the second virial coefficient for the Hard Disc system into a circular pore. Finally, we analyse the behaviour of the obtained expressions near the maximum density.     

Second virial coefficient in one dimension as a function of asymptotic quantities

A result from Dodd and Gibbs (J. Math. Phys., 15, 41 (1974)) for the second virial coefficient of particles in one dimension, subject to delta-function interactions, has been obtained by direct integration of the wave functions. It is shown that this result can be obtained from a phase shift formalism, if one also includes the contribution of oscillating terms. The result is important in work to follow, for the third virial coefficient, for which a similar formalism is being developed. We examine a number of fine points in the quantum mechanical formalisms.     

Critical properties of non-spherical molecule fluids from the virial expansion

Critical constants of pure fluids (as important reference data in constructing vapour-liquid phase diagrams and basic input of various estimation methods) were determined for systems of non-spherical Kihara molecules; values of the critical temperature, density, compression factor and pressure of fluids composed of prolate and oblate molecules were evaluated from the fourth-order virial expansion. The second and third virial coefficients of the Kihara molecules were determined by applying the recently proposed method in which the effect of molecular core geometry and functional dependence of a pair interaction on the surface-surface distance are factorized and the former contribution determined from a formula for the corresponding hard convex body virial coefficient. The virial expansion for non-spherical Kihara molecules is applied to determine the critical constants of n-alkanes (methane to octane) and cyclic hydrocarbons (cyclopentane, cyclohexane, benzene and naphthalene); a fair agreement with experimental data was found.     

Thermodynamics of the Stockmayer fluid in an applied field

The thermodynamic properties of the Stockmayer fluid in an applied field are studied using theory and computer simulation. Theoretical expressions for the second and third virial coefficients are obtained in terms of the dipolar coupling constant (λ, measuring the strength of dipolar interactions as compared to thermal energy) and dipole–field interaction energy (α, being proportional to the applied field strength). These expressions are tested against numerical results obtained by Mayer sampling calculations. The expression for the second virial coefficient contains terms up to λ⁴, and is found to be accurate over realistic ranges of dipole moment and temperature, and over the entire range of the applied field strength (from zero to infinity). The corresponding expression for the third virial coefficient is truncated at λ³, and is not very accurate: higher order terms are very difficult to calculate. The virial coefficients are incorporated in to a thermodynamic theory based on a logarithmic representation of the Helmholtz free energy. This theory is designed to retain the input virial coefficients, and account for some higher order terms in the sense of a resummation. The compressibility factor is obtained from the theory and compared to results from molecular dynamics simulations with a typical value λ = 1. Despite the mathematical approximations of the virial coefficients, the theory captures the effects of the applied field very well. Finally, the vapour–liquid critical parameters are determined from the theory, and compared to published simulation results; the agreement between the theory and simulations is good.     

The virial coefficients of hard hypersphere binary mixtures

The third, fourth and fifth virial coefficients of hard hypersphere binary mixtures with dimensionality d = 4, 5 have been calculated for size ratios R ≥ 0.1, R = ≡ σ₂₂/σ₁₁, where σ _ii is the diameter of component i. The composition independent partial virial coefficients have been evaluated by Monte Carlo integration of the corresponding Mayer modified star diagrams. The results are compared with the predictions of Santos, S., Yuste, S. B., and Lopez de Haro, M., 1999, Molec. Phys., 96, 1 of the equation of state of a multicomponent mixture of hard hyperspheres, and the good agreement gives strong support to the validity of that recipe.