Cavity correlation and bridge functions at high density and near the critical point: a test of second-order Percus–Yevick theory

ABSTRACT

Cavity correlation functions, pair correlation functions, and bridge functions for the Lennard-Jones fluid are calculated from first Percus–Yevick (PY) theory and second-order Percus– Yevick (PY2) theory, molecular dynamics, and grand canonical Monte Carlo techniques. We find that the PY2 theory is significantly more accurate than the PY theory, especially at high density and near the critical point. The pair correlation function near the critical point has the expected slowly decaying long-range behaviour. However, we do not observe any long-range behaviour in the bridge function for the state points near the critical point we have simulated. However, we do note that the bridge function, which is usually negative near r = 0, becomes positive as r → 0. This behaviour is seen for the bridge functions computed from both PY2 and molecular dynamics, but not from PY.     

Non-ergodicity transition and multiple glasses in binary mixtures: on the accuracy of the input static structure in the mode coupling theory

We examine the question of the accuracy of the static correlation functions used as input in the mode coupling theory (MCT) of non-ergodic states in binary mixtures. We first consider hard-sphere mixtures and compute the static pair structure from the Ornstein-Zernike equations with the Percus-Yevick closure and more accurate ones that use bridge functions deduced from Rosenfeld's fundamental measures functional. The corresponding MCT predictions for the non-ergodicity lines and the transitions between multiple glassy states are determined from the long-time limit of the density autocorrelation functions. We find that while the non-ergodicity transition line is not very sensitive to the input static structure, up to diameter ratios D(2)/D(1)?=?10, quantitative differences exist for the transitions between different glasses. The discrepancies with the more accurate closures become even qualitative for sufficiently asymmetric mixtures. They are correlated with the incorrect behavior of the PY structure at high size asymmetry. From the example of ultra-soft potential it is argued that this issue is of general relevance beyond the hard-sphere model.     

Further property of Lennard-Jones fluid: Thermal conductivity

In this paper, we calculate the thermal conductivity of noble gases, methane, and three noble gas mixtures including He+Kr, He+Xe, and Kr+Xe assuming they obey Lennard-Jones (LJ) (12-6) model potential. One of the required quantities to calculate the thermal conductivity of these systems is the pair correlation function. Therefore, we solve numerically the Ornstein-Zernike (OZ) integral equation using the mean spherical approximation (MSA) to obtain the pair correlation functions. We use these functions to obtain the thermal conductivity, then compare our results with the available data. According to the results obtained from the present work for pure and mixtures of LJ fluids reveals that the integral equations method is suitable for predicting the thermal conductivity of this class of fluid.     

Heat capacity of square-well fluids of variable width

We have obtained by Monte Carlo NVT simulations the constant-volume excess heat capacity of square-well fluids for several temperatures, densities and potential widths. Heat capacity is a thermodynamic property much more sensitive to the accuracy of a theory than other thermodynamic quantities, such as the compressibility factor. This is illustrated by comparing the reported simulation data for the heat capacity with the theoretical predictions given by the Barker-Henderson perturbation theory as well as with those given by a non-perturbative theoretical model based on Baxter's solution of the Percus-Yevick integral equation for sticky hard spheres. Both theories give accurate predictions for the equation of state. By contrast, it is found that the Barker-Henderson theory strongly underestimates the excess heat capacity for low to moderate temperatures, whereas a much better agreement between theory and simulation is achieved with the non-perturbative theoretical model, particularly for small well widths, although the accuracy of the latter worsens for high densities and low temperatures, as the well width increases.     

Modified thermodynamic perturbation theory for fused–sphere dimer fluids

Wertheim's thermodynamic perturbation theory of first order (TPT1) is based on the approximation that the monomer–monomer distribution functions can be approximated by the reference fluid distribution functions regardless of the amount of bonding. This is remarkably accurate for chains formed by tangent spheres, but no longer valid for chains of fused spheres. This constitutes the reason for the inadequacy of TPT1 for fused sphere chains. We present a systematic modification of TPT1, the path integral perturbation method, that takes into account the variations of the distribution functions with extent of bonding. We demonstrate the accuracy of the theory for mixtures of hard spheres and diatomics over a range of extent of bonding (pure monomers to pure dimers) and degree of fusion (bond length 0–1). We found that the choice of reference fluid was decisive for the accuracy of the model's predictions. The proposed theory can accurately predict the properties of mixtures of hard spheres and diatomics, and of the pure fused diatomic fluids. The results from the path integral theory are in excellent agreement with simulation results, and compare favourably with the results from the Tildesley–Streett and the Boublík–Nezbeda equations of state.     

Time-dependent properties of the fluid mixture Ne/Kr in the region of the critical mixing point at about 2000 bar

In continuation of our previous molecular dynamics calculations [1], the dynamic structure of the fluid mixture Ne/Kr is investigated in the region of its upper critical mixing point at high pressures. The time-dependent pair-distribution functions have been found to change anomalously near the critical point, relative to the range far above. These changes can be explained by a large increase in the characteristic time of the system close to its critical point, in agreement with experimental findings on other liquids and binary mixtures and theoretical predictions. Our calculations yield a correlation time which is about four times larger than for the non-critical system, in accord with light scattering experiments on similar systems. It is therefore confirmed again that MDC, using effective Lennard-Jones (12–6) potentials, is able to simulate satisfactorily many features of the real Ne/Kr system, including its gas-gas equilibrium of the second kind.     

An experimental investigation of the pressure and concentration dependence of muonic Coulomb capture and cascade in gases

Muonic Coulomb capture ratios and x-ray intensity patterns in four binary mixtures of the gases N₂, Ne, Ar, Kr, and Xe have been measured at three atomic ratios. An influence of the concentration has been established. The Lyman series intensity patterns of pure N₂, Ne, Ar, Kr, and Xe were measured at pressures between 0.4 and 51 bar and found to depend on the pressure. Possible explanations are discussed.     

Excess Gibbs energy for six binary solid solutions of molecularly simple substances

In this paper we apply the method developed in a previous study of Ar + CH₄ to the evaluation of the excess Gibbs energy G^{E S} for solid solutions of two molecularly simple components The method depends on combining information on the excess Gibbs energy G^{E L} for the liquid mixture of the two components with a knowledge of the (T, x) solid-liquid phase diagram Certain thermal properties of the pure substances are also needed G^{E S} has been calculated for binary mixtures of Ar + Kr, Kr + CH₄, CO + N₂, Kr + Xe, Ar + N₂ and Ar + CO. In general, but not always, the solid mixtures are more non-ideal than the liquid mixtures of the same composition at the same temperature Except for the Kr + CH₄ system, the ratio $\text{r =}\text{G}{}^{\mathrm{E\; S}}\text{G}{}^{\mathrm{E\; L}}$ is larger the richer the solution in the component with the smaller molecules     

Electrical and optical characteristics of inert-gas plasmas in spatially nonuniform electrode systems

The characteristics of unipolar and bipolar corona discharges in mixtures He/Ar/Kr, He/Ar/Xe, and He/Kr/Xe at atmosphere pressures in electrode systems of the point-plane and point-grid type distributed along the length. It is shown that bipolar corona discharges can be used in the preionization systems of high-pressure inert-gas-atom lasers. The spectral characteristics of the plasma are investigated, and the dependence of the relative intensity of the radiation in the lines Ar(4s-4p), Kr(5s-5p), and Xe(6s-6p) on the composition and pressure of the working mixtures and on the value of the corona discharge current are determined. Zh. Tekh. Fiz. 67, 15–18 (November 1997)     

TPT2 and SAFTD equations of state for mixtures of hard chain copolymers

We present the second-order thermodynamic perturbation theory (TPT2) and the dimer statistical associating fluid theory (SAFTD) equations of state for mixtures consisting of hetero-nuclear hard chain molecules based on extensions of Wertheim's theory for associating fluids. The second-order perturbation theory, TPT2, is based on the hard sphere mixture reference fluid. SAFTD is an extension of TPT1 (= SAFT) and is based on the non-spherical (hard disphere mixture) reference fluid. The TPT2 equation of state requires only the contact values of the hard sphere mixture site-site correlation functions, while the SAFTD equation of state requires the contact values of site-site correlation functions of both hard sphere and hard disphere mixtures. We test several approximations for site-site correlation functions of hard disphere mixtures and use these in the SAFTD equation of state to predict the compressibility factor of copolymers. Since simulation data are available only for a few pure copolymer systems, theoretical predictions are compared with molecular simulation results for the compressibility factor of pure hard chain copolymer systems. Our comparisons show a very good performance of TPT2, which is found to be more accurate than TPT1 (= SAFT). Using a modified Percus-Yevick site-site correlation function SAFTD is found to represent a significant improvement over SAFT and is slightly more accurate than TPT2. Comparison of SAFTD with generalized Flory dimer (GFD) theory shows that both are equivalent at intermediate to high densities for the compressibility factor of copolymer systems investigated here.     

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.     

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     

Thermodynamic properties of the mixtures of some ionic liquids with alcohols using a simple equation of state

In the present work, we have used a simple equation of state called the GMA EoS to calculate the density of three ionic liquid mixtures including 1-butyl-3-methylimidazolum hexafluorophosphate, [BMIM] [PF₆] + methanol, 1-butyl-3-methylimidazolum tetrafluoroborate, [BMIM] [BF₄] + methanol, and [BMIM] [BF₄] + ethanol at different temperatures, pressures, and compositions. The isothermal compressibility, excess molar volumes, and excess Gibbs molar energy of these mixtures have been computed using this equation of state. The values of statistical parameters show that the GMA EoS can predict these thermodynamic properties very well within the experimental errors. The results show that isothermal compressibility of ionic liquids is lower than alcohols and the effect of temperature and pressure on the isothermal compressibility of ionic liquids is lower than alcohols. The excess molar volumes and excess molar Gibbs energy for these ionic liquid mixtures with alcohols are all negative at various temperatures and pressures over the whole composition range. The results have been interpreted in terms of intermolecular interactions and structural factors of the ionic liquids and alcohols.     

Weighted-density approaches for polymer structure at solid–polymer interfaces

Weighted-density approximations (WDAs), which are based on the weighting function for the second-order direct correlation functions (DCFs) of the uniform polymeric fluids, have been developed to investigate the structural and thermodynamic properties of polymer melts at interfaces. The advantage is the simplicity of calculation of the weighting functions and their accuracies in the applications. They were applied to study the local density distributions and adsorption isotherms of the freely jointed tangent hard-sphere chain, Yukawa chain, and hard-sphere chain mixture in slit pores. The polymer reference interaction model (PRISM) integral equation with the Percus–Yevick (PY) closure has been used to calculate the second-order DCF of the polymeric fluids required as inputs. The mean-field approximation (MFA) has been used to calculate the weighting function for the attractive contribution of a freely jointed tangent Yukawa chain fluid, having attraction among the beads. The calculated results show that (i) for the freely jointed tangent hard-sphere chain, the present theory is in excellent agreement with the computer simulations over a wide range of chain lengths and bulk densities, (ii) the WDA approach for the attraction provides an accurate method for the local density distributions of a freely jointed tangent Yukawa chain fluid, and that (iii) the present theory also yields a reasonably good result for the structural properties of the freely jointed hard-sphere chain mixtures composed of the chain and monomer.     

Freezing and melting of binary mixtures confined in a nanopore

This paper reports on a Grand Canonical Monte Carlo study of the freezing and melting of Lennard–Jones Ar/Kr mixtures confined in a slit pore composed of two strongly attractive structureless walls. For all molar compositions and temperatures, the pore, which has a width of 1.44?nm, accommodates two contact layers and one inner layer. Different wall/fluid interactions are considered, corresponding to pore walls that have a larger affinity for either Ar or Kr. The solid/liquid phase diagram of the confined mixture is determined and results compared with data for the bulk mixture. The structure of the confined mixture is studied using 2D order parameters and both positional g(r) and bond orientational G₆(r) pair correlation functions. It is found that in the confined solid phase, both the contact and inner layers have a hexagonal crystal structure. It is shown that the freezing temperature of the Ar/Kr confined mixture is higher than the bulk freezing point for all molar compositions. Also, it is found that the freezing temperature becomes larger as the ratio α of the wall/fluid to the fluid/fluid interactions increases, in agreement with previous simulation studies on pure substances confined in nanopores. In the case of pore walls having a stronger affinity for Kr atoms (ε _Ar/W<ε _Kr/W), it is observed that both the contact and inner layers of the confined mixture undergo, at the same temperature, a transition from the liquid phase to the crystal phase. The freezing of Ar/Kr mixtures confined between the walls having a stronger affinity for Ar (ε _Ar/W?>?ε _Kr/W) is more complex: for Kr molar concentration lower than 0.35, we observe the presence of an intermediate state between all layers being 2D hexagonal crystals and all the layers being liquid. This intermediate state consists of a crystalline contact layer and a liquid-like inner layer. It is also shown that the qualitative variations of the increase of freezing temperature with the molar composition depend on the affinity of the pore wall for the different components. These results confirm that, in addition to the parameter α the ratio of the wall/fluid interactions for the two species, η=?_Ar/W/?_Kr/W, is a key variable in determining the freezing and melting behaviour of the confined mixture.     

Lyman alpha (1216 Å) generation in krypton

We measured the nonlinear susceptibilities of Kr up to fifth order via nonresonant third-harmonic generation (THG) of a frequency-doubled dye laser at 3648 Å. In agreement with theory the optimum intensity conversion efficiency saturates at 10^-5 in pure Kr. The phase matching of the third harmonic and fundamental wave breaks up owing to intensity-dependent changes of refractive indices. Phase matching in high-density Kr/Ar mixtures raised the 3648 Å to 1216 Å third harmonic conversion obtained in pure Kr by more than two orders of magnitude, as predicted by theory but hitherto never experimentally confirmed. With laser intensities of 5 × 10⁹ W/cm² at 3648 Å we measured conversion efficiencies of 1.4 × 10^-4 limited by gas breakdown. Taking into account the Lyman alpha absorption of Kr₂ and KrAr molecules in the high-density regime, a conversion efficiency of 10^-3 is expected from extrapolating the experimental results to Kr/Ar mixtures of up to 3 × 10³ torr. It is limited by self-focusing, assuming the minimum threshold intensity of gas breakdown is 10¹¹ W/cm², as predicted by theory.     

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.     

矩形均匀电子束在激光介质中的能量沉积的模拟计算

本文主要论述了0.5MeV,单能垂直入射的电子束在三种不同密度的氟化氪准分子激光介质中的能量沉积。计算是采用Monte Carlo方法的MCSED程序,光轴方向采用周期性边界条件,因此能够给出平行于和垂直于电子束入射方向的沉积能量的空间分布。本文给出了在三种不同密度下的出射电子的角分布。用该程序对垂直入射的,初始能量为1MeV的电子在半无限大Al靶中的沉积能量的计算结果与ONETRAN程序的结果及实验结果的比较表明MCSED程序的计算结果是可靠的。     

Longitudinal and bulk viscosities of binary fluid mixtures

Expressions for zeroth, second and fourth sum rules of longitudinal and bulk stress auto correlation functions have been derived for binary fluid mixtures. Longitudinal and bulk viscosities of an Ar–Kr mixture have been calculated using Mori's memory function formalism coupled with the sum rules of longitudinal and bulk stress auto correlation functions. The results obtained are compared with the molecular dynamics simulation. Mass dependence of the longitudinal and bulk viscosities has been studied for different compositions of an isotopic mixture at different densities and temperatures. For very large mass ratio, the longitudinal and bulk viscosities of the isotopic mixture are more dependent on mole fraction than on mass.