Phase coexistence and interface structure of a two-component Lennard-Jones fluid in porous media: application of Born—Green—Yvon equation

The Born-Green-Yvon equation with smoothed density approximation is used to calculate the liquid-liquid density profiles of a symmetric Lennard-Jones fluid in a hard sphere disordered matrix. The phase diagrams are evaluated for model systems characterized by different matrix densities and compared with the results of theoretical predictions and the Monte Carlo simulations of Gordon, P. A., and Glandt, E. D., 1996, J. chem. Phys., 105, 4257. It was found that increasing the matrix packing fraction reduces the magnitude of the miscibility gap and smooths the density profiles between two coexisting phases.     

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.     

Perturbation and variational approach for the equation of state for hard-sphere and Lennardben Jones fluids

The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennard-Jones (LJ) fluids. A suitable flexible functional form for the radial distribution function G(R) is assumed for the mixture, with R as a variable. The function G(R) has an arbitrary parameter m and a different equation of state can be obtained with a suitable choice of m. For m= 0.75 and m= 0.83 results are close to molecular dynamics (MD) result for pure HS and LJ fluid respectively.     

Adsorption on nanostructured chiral surfaces studied by the Monte Carlo method

A Monte Carlo (MC) lattice gas model of adsorption of a racemic mixture of enantiomers of 1,2-dimethylcyclopropane on a chiral surface with different spatial distribution of active sites was proposed. The calculations were performed on a square lattice for both stepped chiral surfaces and smooth surfaces with chiral patterns of active sites. The adsorbing molecules were assumed to be rigid structures of two types being mirror images one of another. Regardless of the enantiomer type, each molecule was composed of four segments occupying four lattice sites. The chiral surfaces were exposed to equimolar mixture of enantiomers whose individual equilibrium adsorption isotherms were calculated using standard Grand Canonical MC technique. The major purpose of the simulation was to examine how the structure of the surface affects separation of enantiomers, that is, to determine enantioselectivity defined as the ratio of their adsorbed amounts. Additionally, comparison of the enantioselectivities corresponding to the stepped and smooth surfaces was made.     

New accurate binary hard sphere mixture radial distribution functions at contact and a new equation of state

New and very accurate formulae for additive binary hard sphere (HS) mixture radial distribution functions (RDFs) at contact are proposed in a simple analytical form. Using the virial theorem, the formulae also provide a new HS mixture equation of state (EOS). The new RDF formulae are the most accurate currently available. The new EOS is of comparable accuracy with that of Malijevsky, A., and Veverka, J. (1999, Phys. Chem. chem. Phys., 1, 4267), which is the most accurate HS mixture EOS currently available. However, the new EOS proposed here is of much simpler analytical form.     

Far-infra-red absorption of non-polar liquids

A new conformal solution theory using a single pure fluid as a reference substance for the calculation of thermodynamic properties of fluid mixtures is developed. The perturbation theory developed by Weeks, Chandler and Andersen (WCA) and by Verlet and Weis (VW) is used to calculate the reference properties. The mean density approximation and corresponding state principle are used to eliminate the higher order terms in the mixture system and to derive the pseudo-parameters for the reference system. The mixture properties are obtained from the reference properties and their corresponding hard sphere excess functions defined as the properties of the mixture less the value of the properties for the hard sphere mixture.

The excess functions of mixing for several liquid mixtures of Lennard-Jones fluids, obeying the Lorentz-Berthelet rule, are calculated by the new method (VW-HSE). Comparison with the results of other theories and Monte Carlo data shows definite improvement. Since only the properties of a pure reference fluid are directly calculated, the method can be applied to more complicated multicomponent systems without additional computational effort as required by other theories.     

Equation of state of mixtures of simple fluids at high pressures

Recently Lee and Levesque have extended the Weeks-Chandler-Andersen pure fluid perturbation theory to mixtures and have compared the results to the Leonard-Henderson-Barker mixture theory. The results seem to favour the Leonard-Henderson-Barker theory. This and other previous comparisons of mixture theories have been mostly confined to the study of ‘zero pressure’ thermodynamic properties. In this paper we compare the Lee-Levesque, Leonard-Henderson-Barker and the van der Waals one-fluid theories to high pressure equation of state data for helium-xenon mixtures. This system is modelled by a binary mixture of Lennard-Jones fluids and the hard sphere reference system is characterized by the Grundke-Henderson hard sphere mixture radial distribution function. The Lee-Levesque theory compares favourably with experimental equation of state data up to pressures of 2000 atmospheres. The Leonard-Henderson-Barker and van der Waals one theories are satisfactory. Although the van der Waals one theory yields the poorest results, it does offer the advantage of having the greatest ease of computation.

Previous theoretical and Monte Carlo calculations of mixture properties have assumed the geometric mean rule for the mixed interaction energy parameter, ε ₁₂, with consequent disagreement with experimental results. We point out that ε ₁₂ can be determined from mixed second virial coefficient data and use such improved determinations of ε ₁₂. We show that this method yields significantly improved theoretical predictions.     

Comments on the analytical solutions of the sticky hard rod model

The result that follows by taking the Percus-Yevick (PY) approximation from the exact analysis of the sticky hard rod model is shown to differ from that previously obtained by Tago, Y., and Katsura, S., 1975, Can. J. Phys. 53, 2587, and the correct PY compressibility and virial equations of state are given. It is also shown that Born-Green theory yields the exact solution for the sticky hard rod fluid.     

New universal scaling laws of diffusion and Kolmogorov-Sinai entropy in simple liquids

A new universal scaling law relating the self-diffusivities of the components of a binary fluid mixture to their excess entropies is derived using mode coupling theory. These scaling laws yield numerical results, for a hard sphere as well as Lennard-Jones fluid mixtures, in excellent agreement with simulation results even at a low density region, where the empirical scaling laws of Dzugutov [Nature (London) 381, 137 (1996)]] and Hoyt, Asta, and Sadigh [Phys. Rev. Lett. 85, 594 (2001)]] fail completely. A new scaling law relating the Kolmogorov-Sinai entropy to the excess entropy is also obtained.     

2D Ising Model for Enantiomer Adsorption on Achiral Surfaces: L- and D-Aspartic Acid on Cu(111)

The 2D Ising model is well-formulated to address problems in adsorption thermodynamics. It is particularly well-suited to describing the adsorption isotherms predicting the surface enantiomeric excess, ees, observed during competitive co-adsorption of enantiomers onto achiral surfaces. Herein, we make the direct one-to-one correspondence between the 2D Ising model Hamiltonian and the Hamiltonian used to describe competitive enantiomer adsorption on achiral surfaces. We then demonstrate that adsorption from racemic mixtures of enantiomers and adsorption of prochiral molecules are directly analogous to the Ising model with no applied magnetic field, i.e., the enantiomeric excess on chiral surfaces can be predicted using Onsager’s solution to the 2D Ising model. The implication is that enantiomeric purity on the surface can be achieved during equilibrium exposure of prochiral compounds or racemic mixtures of enantiomers to achiral surfaces.     

Virial coefficients and equations of state for mixtures of hard discs,hard spheres and hard hyperspheres

The composition-independent virial coefficients of a d-dimensional binary mixture of (additive) hard hyperspheres following from a recent proposal for the equation of state of the mixture (Santos, A., Yuste, S. B., and López de Haro, M., 1999, Molec. Phys., 96, 1) are examined. Good agreement between theoretical estimates and available exact or numerical results is found for d = 2, 3, 4 and 5, except for mixtures whose components are very disparate in size. A slight modification that remedies this deficiency is introduced and the resummation of the associated virial series is carried out, leading to a new proposal for the equation of state. The case of binary hard sphere mixtures (d = 3) is analysed in some detail.     

Chemical reactions at surfaces: Application of the reaction ensamble monte carlo method

The Monte Carlo simulation method introduced by Smith and Triska [J. Chem. Phys.100 (1994) 3019] is extended to the case of a reacting fluid in contact with a hard wall. The fluid structure for both spherical and nonspherical reaction products is discussed for simple models of reacting hard spheres near a hard wall and near a wall interacting via Lennard-Jones (9,3) potential. In the latter case the investigated model assumes that the probability of a chemical reaction changes with a distance from the surface. It is shown that the applied technique is suitable for the study of reacting nonuniform fluids. This work is supported by KBN under the Grant No. 3 T09A 062 10.     

Is it possible to infer the equation of state of a mixture of hard discs from that of the one-component system?

Based on exact asymptotic properties of the composition-independent virial coefficients of a binary mixture of hard discs in the limits α = σ₂/σ₁ → 0, α → 1 and α → ∞, R. J. Wheatley (1998, Molec. Phys., 93, 965) has recently proposed an approximate interpolation equation for these coefficients. In this note, the equation of state equivalent to this interpolation is obtained, expressing the compressibility factor of the mixture in terms of that of the pure system. An extension to an arbitrary number of components is also given. The equation of state derived here is compared with another one recently proposed by following a different route (Santos, A., Yuste, S. B., and López de Haro, M., 1999, Molec. Phys., 96, 1) and with Monte Carlo simulation results. It is shown that the latter equation is more accurate than the former one, at least for not too disparate mixtures (0.7 < α < 1).     

Charged hard dumbbell in the binding mean-spherical approximation

A simple model for charged hard dumbbell is proposed in the binding mean-spherical approximation (BIMSA). The thermodynamic properties are analytical solutions of the unique screening parameter Γ^B with full association. Critical point and vapour-liquid coexistence curve are identical to those of Kalyuzhnyi, Yu. V., 1998, Molec. Phys., 94, 735, where a site-site integral equation has to be solved. Substituting Γ without association for Γ^B, the BIMSA reduces to the simple interpolation scheme (SIS). A simple interpolation between the SIS and the BIMSA is proposed: this gives the critical point (T_c? = 0.0525, p?_c = 0.0640) which, for the time being, is the closest to the computer simulation results. Similarity between the charged hard dumbbell and the restricted primitive model of electrolyte is also addressed.     

Diffusion within a near-critical nanopore fluid

Results are presented for grand canonical Monte Carlo (GCMC) and both equilibrium and non-equilibrium molecular dynamics simulations (EMD and NEMD) conducted over a range of densities and temperatures that span the two-phase coexistence and supercritical regions for a pure fluid adsorbed within a model crystalline nanopore. The GCMC simulations provided the low temperature coexistence points for the open pore fluid and were used to locate the capillary critical temperature for the system. The equilibrium configurational states obtained from these simulations were then used as input data for the EMD simulations in which the self-diffusion coefficients were computed using the Einstein equation. NEMD colour diffusion simulations were also conducted to validate the use of a system averaged Einstein analysis for this inhomogeneous fluid. In all cases excellent agreement was observed between the equilibrium (linear response theory) predictions for the diffusivities and non-equilibrium colour diffusivities. The simulation results are also compared with a recently published quasi-hydrodynamic theory of Pozhar and Gubbins (Pozhar, L. A., and Gubbins, K. E., 1993, J. Chem. Phys., 99, 8970; 1997, Phys. Rev. E, 56, 5367.). The model fluid and the nature of the fluid wall interactions employed conform to the decomposition of the particle–particle interaction potential explicitly used by Pozhar and Gubbins. The local self-diffusivity was calculated from the local fluid–fluid and fluid wall hard core collision frequencies. While this theory provides reasonable results at moderate pore fluid densities, poor agreement is observed in the low density limit.     

Dielectric properties of a mixture of dipolar hard spheres

A theory for the dielectric constant, ε, of a fluid mixture of dipolar hard spheres is formulated by generalizing the methods developed by Ramshaw and Wertheim for the pure fluid case. The resulting expression for ε depends on the pair distribution functions, g _αβ(r ₁, θ₁, r ₂, θ₂) for a dipolar mixture. Due to the unavailability of exact representations for these dipolar pair distribution functions, the results of the mean spherical approximation are employed in the formalism developed. Numerical results are given for ε as calculated from the pair distribution functions for a spherical volume of macroscopic dimensions. The compositional dependence of the ε obtained in this way for a specific mixture is compared with the corresponding properties of the well established theories of Clausius-Mossotti-Debye and Onsager. In addition, the relative importance of the dipole moment and size of the hard sphere parameters in determining ε for a dipolar mixture (the correlative behaviour of which is described by the mean spherical approximation) is evaluated. It is found that the differences in hard core diameters can be largely ignored, in that ε for an ‘effective’ single component fluid can be given to within 2–5 per cent relative error (at worst) of the mean spherical approximation's result. Such an ‘effective pure fluid’ is described as having the same polarization content as the actual mixture being considered. Thereby, the properties of the effective fluid are determined by the quantity y = 4πβ(m ₁ ² ρ₁ + m ₂ ² ρ₂)/9 where m_i and ρ _i are the dipole moment and number density of component i in the binary mixture, with β = (kT)^-1.     

The stability limit of the fluid phase of polydisperse sticky spheres

It has been shown by Stell (1991, J. statist. Phys., 63, 1203) that at low temperature mono-disperse sticky spheres collapse to form coexisting close-packed solid and infinitely dilute gases. We show that polydisperse sticky spheres also collapse and calculate the collapse temperature. The polydisperse spheres separate into fractions with narrower polydispersities which can then solidify. This is perhaps the first example of a single-peaked polydisperse mixture phase solidifying and separating. It implies that a mixture of polydisperse large hard spheres with much smaller hard spheres does not show fluid—fluid coexistence.     

Statistical fluid theory for associating fluids containing alternating heteronuclear chain molecules

Statistical associating fluid theory of homonuclear dimerized chain fluids and homonuclear monomer-dimer mixture chain fluids are extended to fluids containing alternating heteronuclear chain molecules separately. The proposed models account for the appropriate site-site correlation functions at contact. The modified equations of state show a good agreement with generalized Flory dimer theory and MD simulation data for small and medium size ratio of hard sphere diameters.     

Transport coefficients of hard sphere fluids

New calculations have been made of the self-diffusion coefficient D, the shear viscosity η_s, the bulk viscosity η_b and thermal conductivity λ of the hard sphere fluid, using molecular dynamics (MD) computer simulation. A newly developed hard sphere MD scheme was used to model the hard sphere fluid over a wide range up to the glass transition (～0.57 packing fraction). System sizes of up to 32 000 hard spheres were considered. This set of transport coefficient data was combined with others taken from the literature to test a number of previously proposed analytical formulae for these quantities together with some new ones given here. Only the self-diffusion coefficient showed any substantial N dependence for N < 500 at equilibrium fluid densities (ε 0.494). D increased with N, especially at intermediate densities in the range ε ～ 0.3–0.35. The expression for the packing fraction dependence of D proposed by Speedy, R. J., 1987, Molec. Phys., 62, 509 was shown to fit these data well for N ～ 500 particle systems. We found that the packing fraction ε dependence of the two viscosities and thermal conductivity, generically denoted by X, were represented well by the simple formula X/X ₀ = 1/[1 ? (ε/ε₁)]^m within the equilibrium fluid range 0 > ε > 0.493. This formula has two disposable parameters, ε and m, and X ₀ is the value of the property X in the limit of zero density. This expression has the same form as the Krieger-Dougherty formula (Kreiger, I. M., 1972, Adv. Colloid. Interface Sci., 3, 111) which is used widely in the colloid literature to represent the packing fraction dependence of the Newtonian shear viscosity of monodisperse colloidal near-hard spheres. Of course, in the present case, X _o was the dilute gas transport coefficient of the pure liquid rather than the solvent viscosity. It was not possible to fit the transport coefficient normalized by their Enskog values with such a simple expression because these ratios are typically of order unity until quite high packing fractions and then diverge rapidly at higher values over a relatively narrow density range. At the maximum equilibrium fluid packing fraction ε = 0.494 for both the hard sphere fluid and the corresponding colloidal case a very similar value was found for η_s/η_o ?30–40, suggesting that the ‘crowding’ effects and their consequences for the dynamics in this region of the phase diagram in the two types of liquid have much in common. For the hard sphere by MD, D_o/D ～ 11 at the same packing fraction, possibly indicating the contribution from ‘hydrodynamic enhancement’ of this transport coefficient, which is largely absent for the shear viscosity. Interestingly the comparable ratio for hard sphere colloids is the same.     

Three-phase osmotic equilibria using the Gibbs ensemble simulation method

The Gibbs Ensemble Monte Carlo method has been used to simulate osmotic equilibria for Lennard-Jones mixtures. When the simulations are performed with two independent boxes, one containing solvent and the other a mixture of solute and solvent significantly negative osmotic pressures (Π) develop. Following a sugestion of Powles et al. (1997, Molec. Phys., 90, 665), we have extended these simulations to include a third box and the possibility of modelling three coexisting phases. The simulations show that the two phase equilibria with negative values of 77 are metastable and that the system spontaneously separates into three phases: pure solvent, dilute solute-solvent and dense solute-solvent with a resulting osmotic pressure that is normally small and positive.