Generic van der Waals equation of state, modified free volume theory of diffusion, and viscosity of simple liquids

The shear viscosity formula derived by the density fluctuation theory in previous papers is computed for argon, krypton, and methane by using the self-diffusion coefficients derived in the modified free volume theory with the help of the generic van der Waals equation of state. In the temperature regime near or above the critical temperature, the density dependence of the shear viscosity can be accounted for by ab initio calculations with the self-diffusion coefficients provided by the modified free volume theory if the minimum (critical) free volume is set equal to the molecular volume and the volume overlap parameter (alpha) is taken about unity in the expression for the self-diffusion coefficient. In the subcritical temperature regime, if the density fluctuation range parameter is chosen appropriately at a temperature, then the resulting expression for the shear viscosity can well account for its density and temperature dependence over the ranges of density and temperature experimentally studied. In the sense that once the density fluctuation range is fixed at a temperature, the theory can account for the experimental data at other subcritical temperatures on the basis of the intermolecular force only; the theory is predictive even in the subcritical regime of temperature. Theory is successfully tested in comparison with experimental data for self-diffusion coefficients and shear viscosity for argon, krypton, and methane.     

Relations between transport coefficients and their density and temperature dependence

Nonequilibrium statistical mechanics via density fluctuation theory predicts relations between the bulk and shear viscosity, thermal conductivity, and self-diffusion coefficient of a fluid. In this Feature Article, we discuss such relations holding for fluids over wide ranges of density and temperature experimentally studied in the laboratory. It is discussed how such relations can be used to successfully compute the density and temperature dependence on the basis of intermolecular interaction potential models with the help of the modified free volume theory and the generic van der Waals equation of state once the parameters in them are determined at a low density or at a subcritical temperature. Although some approximations have been made to derive them, they represent a reliable molecular theory of transport coefficients over the entire density and temperature ranges of fluids--namely, gases and liquids--a theory hitherto unavailable in the kinetic theory of liquids and dense gases.     

Voids, generic van der Waals equation of state, and transport coefficients of liquids

In this Perspective, we discuss the role of voids in transport processes in liquids and the manner in which the concept of voids enters the generic van der Waals equation of state and the modified free volume theory. The density fluctuation theory is then discussed and we show how the density fluctuation theory can be made a molecular theory with the help of the modified free volume theory and the generic van der Waals equation of state. The confluence of the aforementioned three theories makes it possible to calculate the transport coefficients of liquids by using the information on the equilibrium pair correlation function, which can be calculated either by an integral equation theory or Monte Carlo simulations. A number of relations between transport coefficients are also presented, which are derived on the basis of the density fluctuation theory. Since they can be used to obtain one transport coefficient from another they can be very useful in handling experimental and theoretical data. An application of the modified free volume theory to polymer melts is discussed as an example for a theory of transport properties of complex liquids.     

Statistical-mechanical theory of rheology: Lennard-Jones fluids

The generalized Boltzmann equation for simple dense fluids gives rise to the stress tensor evolution equation as a constitutive equation of generalized hydrodynamics for fluids far removed from equilibrium. It is possible to derive a formula for the non-Newtonian shear viscosity of the simple fluid from the stress tensor evolution equation in a suitable flow configuration. The non-Newtonian viscosity formula derived is applied to calculate the non-Newtonian viscosity as a function of the shear rate by means of statistical mechanics in the case of the Lennard-Jones fluid. For that purpose we have used the density-fluctuation theory for the Newtonian viscosity, the modified free volume theory for the self-diffusion coefficient, and the generic van der Waals equation of state to compute the mean free volume appearing in the modified free volume theory. Monte Carlo simulations are used to calculate the pair-correlation function appearing in the generic van der Waals equation of state and shear viscosity formula. To validate the Newtonian viscosity formula obtained we first have examined the density and temperature dependences of the shear viscosity in both subcritical and supercritical regions and compared them with molecular-dynamic simulation results. With the Newtonian shear viscosity and thermodynamic quantities so computed we then have calculated the shear rate dependence of the non-Newtonian shear viscosity and compared it with molecular-dynamics simulation results. The non-Newtonian viscosity formula is a universal function of the product of reduced shear rate (gamma*) times reduced relaxation time (taue*) that is independent of the material parameters, suggesting a possibility of the existence of rheological corresponding states of reduced density, temperature, and shear rate. When the simulation data are reduced appropriately and plotted against taue*gamma* they are found clustered around the reduced (universal) non-Newtonian viscosity formula. Thus we now have a molecular theory of non-Newtonian shear viscosity for the Lennard-Jones fluid, which can be implemented with a Monte Carlo simulation method for the pair-correlation function.     

Modified free volume theory of self-diffusion and molecular theory of shear viscosity of liquid carbon dioxide

In previous work on the density fluctuation theory of transport coefficients of liquids, it was necessary to use empirical self-diffusion coefficients to calculate the transport coefficients (e.g., shear viscosity of carbon dioxide). In this work, the necessity of empirical input of the self-diffusion coefficients in the calculation of shear viscosity is removed, and the theory is thus made a self-contained molecular theory of transport coefficients of liquids, albeit it contains an empirical parameter in the subcritical regime. The required self-diffusion coefficients of liquid carbon dioxide are calculated by using the modified free volume theory for which the generic van der Waals equation of state and Monte Carlo simulations are combined to accurately compute the mean free volume by means of statistical mechanics. They have been computed as a function of density along four different isotherms and isobars. A Lennard-Jones site-site interaction potential was used to model the molecular carbon dioxide interaction. The density and temperature dependence of the theoretical self-diffusion coefficients are shown to be in excellent agreement with experimental data when the minimum critical free volume is identified with the molecular volume. The self-diffusion coefficients thus computed are then used to compute the density and temperature dependence of the shear viscosity of liquid carbon dioxide by employing the density fluctuation theory formula for shear viscosity as reported in an earlier paper (J. Chem. Phys. 2000, 112, 7118). The theoretical shear viscosity is shown to be robust and yields excellent density and temperature dependence for carbon dioxide. The pair correlation function appearing in the theory has been computed by Monte Carlo simulations.     

A new correlation on predicting self- and mutual-diffusion coefficient of Lennard–Jones chain fluid

Based on the Chapman–Enskog theory of diffusion and molecular dynamics simulation data for Lennard–Jones chain (LJC) fluid, a new semi-empirical correlation for calculating the self-diffusion coefficient of LJC fluid is proposed. The new correlation introduces in two correction functions with six fitting parameter to modify the impact of intermolecular repulsive and attractive potential energy on molecular friction coefficient. The new correlation represents the experimental self-diffusion coefficients with an average absolute deviation (AAD) of 3.46% for 23 polyatomic compounds (1102 experimental data points) over wide ranges of temperature and pressure. On this basis, the van der Waals mixing rule is adopted to calculate the mutual-diffusion coefficient of binary LJC fluid. By comparison of calculated results with the experimental data of 12 binary LJC systems over wide range of temperature and composition, the average absolute deviation is 6.98% which verifies the accuracy and the effectiveness of the new correlation.     

A perturbation method for the Ornstein-Zernike equation and the generic van der Waals equation of state for a square well potential model

We calculate the generic van der Waals parameters A and B for a square well model by means of a perturbation theory. To calculate the pair distribution function or the cavity function necessary for the calculation of A and B, we have used the Percus-Yevick integral equation, which is put into an equivalent form by means of the Wiener-Hopf method. This latter method produces a pair of integral equations, which are solved by a perturbation method treating the Mayer function or the well width or the functions in the square well region exterior to the hard core as the perturbation. In the end, the Mayer function times the well width is identified as the perturbation parameter in the present method. In this sense, the present perturbation method is distinct from the existing thermodynamic perturbation theory, which expands the Helmholtz free energy in a perturbation series with the inverse temperature treated as an expansion parameter. The generic van der Waals parameters are explicitly calculated in analytic form as functions of reduced temperature and density. The van der Waals parameters are recovered from them in the limits of vanishing density and high temperature. The equation of state thus obtained is tested against Monte Carlo simulation results and found reliable, provided that the temperature is in the supercritical regime. By scaling the packing fraction with a temperature-dependent hard core, it is suggested to construct an equation of state for fluids with a temperature-dependent hard core that mimicks a soft core repulsive force on the basis of the equation of state derived for the square well model.     

Molecular theory of thermal conductivity of the Lennard-Jones fluid

In this paper the thermal conductivity of the Lennard-Jones fluid is calculated by applying the combination of the density-fluctuation theory, the modified free volume theory of diffusion, and the generic van der Waals equation of state. A Monte Carlo simulation method is used to compute the equilibrium pair-correlation function necessary for computing the mean free volume and the coefficient in the potential-energy and virial contributions to the thermal conductivity. The theoretical results are compared with our own molecular dynamics simulation results and with those reported in the literature. They agree in good accuracy over wide ranges of density and temperature examined in molecular dynamics simulations. Thus the combined theory represents a molecular theory of thermal conductivity of the Lennard-Jones fluid and by extension simple fluids, which enables us to compute the nonequilibrium quantity by means of the Monte Carlo simulations for the equilibrium pair-correlation function.     

基于扩散势能和自由体积的自扩散系数模型

A new model for self-diffusion coe±cients was proposed based on both the concepts of molecular free volume and activation energy. The unknown parameters of this model were clearly defined and compared with the Chapman-Enskog model. At the same time a new method for calculating activation energy was devised and applied to the new model. In addition, the free volume was defined by implementing the generic van der Waals equation of state, the radial distribution function of which was obtained by using the MorsaliGoharshadi empirical formula. Under the same conditions, the new model was better than the original free volume model.     

The use of the cohen–turnbull model for calculation of the self-diffusion coefficients of liquid dichloroalkanes

The paper presents the self-diffusion coefficients calculated for liquid dichloroalkanes C₆H₁₂Cl₂, C₈H₁₆Cl₂, C₁₀H₂₂Cl₂ and C₁₂H₂₄Cl₂, with the use of the Cohen and Turnbull model. Determination of self-diffusion coefficients permits a separate analysis of intra- and intermolecular motions and provides information on geometrical and dynamical properties of molecules. The self-diffusion coefficients of selected dichloroalkanes have been determined by X-ray diffraction and compared with the corresponding NMR results. The suitability of the Cohen–Turnbull model of the translating motion for prediction of self-diffusion coefficients for molecules whose shape significantly differs from the spherical symmetry is analysed. Angular distributions of X-ray scattered intensity were measured, and differential radial distribution functions of electron density (DRDFs) were calculated. The mean coordination numbers were obtained from the area delimited by the minima of the DRDFs, and their dependence on the length of the methylene chain is also presented subsequently. On the basis of the DRDFs the average free volume of the molecules and total free volume of the liquids were calculated. The activation volume of the diffusion was found to make about 0.6 of the van der Waals volume of the molecule. As expected the diffusion coefficients decrease with increasing molecular weight. The equation relating the self-diffusion coefficient with the volume of the coordination spheres in the liquid has been derived.     

A note on excess Gibbs energy models,equations of state and the local composition concept

A density-dependent local composition expression for the residual energy is derived from a generalized NRTL expression for the excess energy and the van der Waals fluid theory. Integration of this expression yields a volume-dependent expression for the Helmholtz energy from which equations of state utilizing the local composition concept are derived and which in the high-density limit contain the well-known activity coefficient models.The local composition versions of the Carnahan—Starling—van der Waals, the Redlich—Kwong—Soave and the Peng—Robinson equations of state are derived. It is further shown that the group contribution versions of the NRTL, the Wilson and the UNIQUAC excess models may be derived from the generalized NRTL expression for the residual energy when applied to groups instead of molecules.It is thus demonstrated that all current local composition activity-coefficient models can be derived from a local composition version of the van der Waals equation of state using different sets of assumptions. In the same way the van Laar, the Scatchard—Hildebrand and the Flory—Huggins activity coefficient models are obtained from the van der Waals equation of state using the original mixing rules.     

Computational studies on thermodynamic properties, effective diameters, and free volume of argon using an ab initio potential

A quantum mechanical derived ab initio interaction potential for the argon dimer was tested in molecular simulations to reproduce the thermophysical properties of the vapor-liquid phase equilibria using the Gibbs ensemble Monte Carlo simulations as well as the liquid and supercritical equation of state using the NVT Monte Carlo simulations. The ab initio interaction potential was taken from the literature. A recently developed theory [R. Laghaei et al., J. Chem. Phys. 124, 154502 (2006)] was used to compute the effective diameters of argon in fluid phases and the results were subsequently applied in the generic van der Waals theory to compute the free volume of argon. The calculated densities of the coexisting phases, the vapor pressure, and the equation of state show excellent agreement with experimental values. The effective diameters and free volumes of argon are given over a wide range of densities and temperatures. An empirical formula was used to fit the effective diameters as a function of density and temperature. The computed free volume will be used in future investigations to calculate the transport properties of argon.     

On the decisive role of finite chain extensibility and global interactions in networks

The stress-strain dependence of dry networks at unidirectional extension and compression is studied. The phenomenological van der Waals equation of state is compared with different molecular models in order to provide an interpretation of the van der Waals corrections. It is shown that the stress-strain behavior predicted by the phantom, Langevin, and constrained junction fluctuation models are altogether covered by the van der Waals approach. The relationship between the suppression of junction fluctuation parameter introduced by Dossin and Graessley and the van der Waals corrections has been worked out. The effect of junction functionality on the small strain modulus as well as on the second Mooney-Rivlin coefficient is also presented.     

The Effect of van der Waals Forces on the Deposition of Highly Dispersed Aerosol Particles on Ultrafine Fibers

The effect of van der Waals forces on the collection of highly dispersed aerosol particles with ultrafine fiber filters was studied theoretically. The capture coefficient was found from the numerical solution of the equation of convective diffusion with the account of the particle size, the effect of van der Waals forces acting between a particle and a fiber, and the gas slip effect at the surface of ultrafine fibers. It was shown that allowance being made for van der Waals forces markedly affects the capture coefficient within the maximal particle penetration range and that the radius of the most penetrating particles decreases with the rising effect of these forces.     

Corresponding-states laws for protein solutions

The solvent around protein molecules in solutions is structured and this structuring introduces a repulsion in the intermolecular interaction potential at intermediate separations. We use Monte Carlo simulations with isotropic, pair-additive systems interacting with such potentials. We test if the liquid-liquid and liquid-solid phase lines in model protein solutions can be predicted from universal curves and a pair of experimentally determined parameters, as done for atomic and colloid materials using several laws of corresponding states. As predictors, we test three properties at the critical point for liquid-liquid separation: temperature, as in the original van der Waals law, the second virial coefficient, and a modified second virial coefficient, all paired with the critical volume fraction. We find that the van der Waals law is best obeyed and appears more general than its original formulation: A single universal curve describes all tested nonconformal isotropic pair-additive systems. Published experimental data for the liquid-liquid equilibrium for several proteins at various conditions follow a single van der Waals curve. For the solid-liquid equilibrium, we find that no single system property serves as its predictor. We go beyond corresponding-states correlations and put forth semiempirical laws, which allow prediction of the critical temperature and volume fraction solely based on the range of attraction of the intermolecular interaction potential.     

Theory of excluded volume equation of state: higher approximations and new generation of equations of state for entire density range

A novel theory of an equation of state based on excluded volume and formulated in two preceding papers for gases and gaseous mixtures is extended to the entire density range by considering higher (beginning from the third) approximations of the theory. The algorithm of constructing higher approximations is elaborated. Equations of state are deduced using the requirement of maximum simplicity and contain a single free parameter to be chosen by reason of convenience or simplicity or to be used as a fitting parameter with respect to the computer simulation database. In this way, precise equations of state are derived for the hard-sphere fluid in the entire density range. On the side, the theory reproduces most known earlier equations of state for hard spheres and determines their place in the hierarchy of approximations. Equations of state for van der Waals fluids are also presented, and their critical parameters are estimated.     

Monte Carlo simulations of thermodynamic and structural properties of Mie(14,7) fluids

The vapor-liquid phase envelope of Mie(14,7) fluids is determined by the Gibbs ensemble Monte Carlo (MC) simulation technique. The NVT-MC simulation method is then utilized to compute the equation of state and the pair correlation function over a wide range of densities and temperatures. The effective diameters are calculated via the virial minimization method and the results are applied as the repulsion-attraction splitting distance within the generic van der Waals (GvdW) theory to compute the mean free volume. The density and temperature dependence of these parameters are studied and discussed. The results for the effective diameter, and the GvdW parameters are fitted to analytical functions of density and temperature. An examination of the results for the fluid phase equilibria of argon shows excellent agreement with empirical data for the densities of the coexisting phases, the vapor pressure, and the critical point. The computed free volumes are used to compute the diffusion coefficient of argon and the results are compared with experimental data.     

Two-step evaluation of binding energy and potential energy surface of van der Waals complexes

Evaluation of intermolecular distance and binding energy (BE) of van der Waals complex/cluster at ab initio level of theory is computationally demanding when many monomers are involved. Starting from MP2 energy, we reached a two-step evaluation method of BE of van der Waals complex/cluster through reasonable approximations; BE = BE(HF) + sum Mi> Mj{BE (Mi- Mj)(MP2 or MP2.5) - BE(Mi-Mj)(HF)} where HF represents the Hartree-Fock calculation, Mi, Mj, etc. are interacting monomers, and MP2.5 represents the arithmetic mean of MP2 and MP3. The first term is the usual BE of the complex/cluster evaluated at the HF level. The second term is the sum of the difference in two-body BE between the correlated and HF levels of theory. This equation was applied to various van der Waals complexes consisting of up-to-four monomers at MP2 and MP2.5 levels of theory. We found that this method is capable of providing precise estimate of the BE and reproducing well the potential energy surface of van der Waals complexes/clusters; the maximum error of the BE is less than 1 kcal/mol and 1% in most cases except for several limited cases. The origins of error in these cases are discussed in detail.     

Finite extensibility and orientational effects in rubbers. A physical interpretation of the ‘van der Waals equation’ for the elastic force

We propose a physical interpretation of the so-called van der Waals equation of state for rubbers, which gives a relation between the force and the deformation. On a phenomenological basis this equation takes the finite extensibility and a non-defined interaction into account. Here the fininte extensibility is discussed for the dilute case (no entanglements) and the highly entangled limit. The intramolecular interactions are described by orientational effects. The resulting equation of state for the force shows the same features as the van der Waals equation.Dedicated to Prof. H.-G. Kilian on the occasion of his 60th birthday.