Virial derivation

In 1901, H. Kamerlingh Onnes introduced a fundamentally new and improved description of real gas PVT properties in terms of the virial equation of state. [The word virial, deriving from the Latin word viris ( force ) was introduced into physics by R. Clausius, whom we shall meet later.] This equation expresses the compressibility factor Z(Vm, T) in terms of a general power series expansion in inverse molar volume Vm. The starting point for the virial expansion is the ideal limiting behavior (2.12), which can also be expressed as... 

Enthalpies are referred to the ideal vapor. The enthalpy of the real vapor is found from zero-pressure heat capacities and from the virial equation of state for non-associated species or, for vapors containing highly dimerized vapors (e.g. organic acids), from the chemical theory of vapor imperfections, as discussed in Chapter 3. For pure components, liquid-phase enthalpies (relative to the ideal vapor) are found from differentiation of the zero-pressure standard-state fugacities these, in turn, are determined from vapor-pressure data, from vapor-phase corrections and liquid-phase densities. If good experimental data are used to determine the standard-state fugacity, the derivative gives enthalpies of liquids to nearly the same precision as that obtained with calorimetric data, and provides reliable heats of vaporization. 

DERIVATIVE OF NETASTABLE, BOUND AND CHEMICAL CONTRIBUTIONS TO VIRIAL... 

Key = 1 represents an initial calculation for a new system Key 2-5 are subsequent calculations not differing significantly in time requirements Key = 6,7 require temperature derivatives of virial coefficients. 

A direct and transparent derivation of the second virial coefficient follows from the canonical ensemble. To make the notation and argument simpler, we first assume pairwise additivity of the total potential with no angular contribution. The extension to angularly-mdependent non-pairwise additive potentials is straightforward. The total potential... 

The virial pressure equation for hard spheres has a simple fomr detemiined by the density p, the hard sphere diameter a and the distribution fimction at contact g(c+). The derivative of the hard sphere potential is discontinuous at r = o, and... 

This has the fomi of a second virial coefficient in which the Debye screened potential has replaced the Coulomb potential. Expressions for the other excess themiodynamic properties are easily derived. 

Statistical mechanics may be used to derive practical microscopic fomuilae for themiodynamic quantities. A well-known example is tire virial expression for the pressure, easily derived by scaling the atomic coordinates in the canonical ensemble partition fiinction... 

Thus one must rely on macroscopic theories and empirical adjustments for the determination of potentials of mean force. Such empirical adjustments use free energy data as solubilities, partition coefficients, virial coefficients, phase diagrams, etc., while the frictional terms are derived from diffusion coefficients and macroscopic theories for hydrodynamic interactions. In this whole field of enquiry progress is slow and much work (and thought ) will be needed in the future. 

The pressure is usually calculated in a computer simulation via the virial theorem ol Clausius. The virial is defined as the expectation value of the sum of the products of the coordinates of the particles and the forces acting on them. This is usually written iV = X] Pxi where x, is a coordinate (e.g. the x ox y coordinate of an atom) and p. is the first derivative of the momentum along that coordinate pi is the force, by Newton s second law). The virial theorem states that the virial is equal to —3Nk T. 

Although the virial equation itself is easily rationalized on empirical grounds, the mixing rules of Eqs. (4-183) and (4-184) follow rigorously from the methods of statistical mechanics. The temperature derivatives of B and C are given exactly by... 

Although developed for pure materials, this correlation can be extended to gas or vapor mixtures. Basic to this extension is the mixing rule for second virial coefficients and its temperature derivative ... 

When i = J, all equations reduce to the appropriate values for a pure species. When i j, these equations define a set of interaction parameters having no physical significance. For a mixture, values of By and dBjj/dT from Eqs. (4-212) and (4-213) are substituted into Eqs. (4-183) and (4-185) to provide values of the mixture second virial coefficient B and its temperature derivative. Values of and for the mixture are then given by Eqs. (4-193) and (4-194), and values of In i for the component fugacity coefficients are given by Eq. (4-196). 

Gamma/Phi Approach For many XT E systems of interest the pressure is low enough that a relatively simple equation of state, such as the two-term virial equation, is satisfactoiy for the vapor phase. Liquid-phase behavior, on the other hand, may be conveniently described by an equation for the excess Gibbs energy, from which activity coefficients are derived. The fugacity of species i in the liquid phase is then given by Eq. (4-102), written... 

However, if one focuses on the adsorption of a fluid in heterogenous matrices [32,33] and/or on the fluctuations in an adsorbed fluid, it is inevitable to perform developments similar to those above in the grand canonical ensemble. Moreover, this derivation is of importance for the formulation of the virial route to thermodynamics of partially quenched systems. For this purpose, we include only some basic relations of this approach. 

In spite of its simplicity and the visual similarity of this equation to Eq. (7), we would like to note that Eq. (11) leads to a nontrivial thermodynamics of a partially quenched system in terms of correlation functions, see, e.g.. Ref. 25 for detailed discussion. Evidently, the principal route for and to the virial theorem is to exploit the thermodynamics of the replicated system. However, special care must be taken then, because the V and s derivatives do not commute. Moreover, the presence of two different temperatures, Pq and P, requires attention in taking temperature derivatives, setting those temperatures equal, if appropriate, only at the end of the calculations. 

The equilibrium between a pure solid and a gaseous mixture is one of very few classes of solution for which an exact treatment can be made by the methods of statistical mechanics. The earliest work on the theory of such solutions was based on empirical equations, such as those of van der Waals,45 of Keyes,44 and of Beattie and Bridgemann.3 However, the only equation of state of a gas mixture that can be derived rigorously is the virial expansion,46 66... 

Equation 8 may be fitted to those results just described for which the vapor pressure of the pure solid is known. We show graphically the second virial coefficients derived from such fitting and those derived from conventional p-V-T measurements. 

Figure 15 shows the second virial coefficients derived independently by us and by Reuss and Beenakker66 from the measurements of solubility by Dokoupil, van Soest, and Swenker,18 and the coefficients at room temperature from the conventional measurements of Verschoyle,86 Michels and Wassenaar,49 and Michels and Boer-boom.47 These results are sufficient to give unambiguously the parameters of a 12-6 potential... 

Figure 17 shows two sets of virial coefficients—derived by us and by Reuss and Beenakker56 from the measurements of Dokoupil,... 

In the derivation above, we have included the kinetic energy of the nuclei in the Hamiltonian and considered a stationary state. In Eq. II.3, this term has been neglected, and we have instead assumed that the nuclei have given fixed positions. It has been pointed out by Slater34 that, if the nuclei are not situated in the proper equilibrium positions, the virial theorem will appear in a slightly different form. (A variational derivation has been given by Hirschfelder and Kincaid.11)... 

The crucial question is at what value of <)> is the attraction high enough to induce phase separation De Hek and Vrij (6) assume that the critical flocculation concentration is equivalent to the phase separation condition defined by the spinodal point. From the pair potential between two hard spheres in a polymer solution they calculate the second virial coefficient B2 for the particles, and derive from the spinodal condition that if B2 = 1/2 (where is the volume fraction of particles in the dispersion) phase separation occurs. For a system in thermodynamic equilibrium, two phases coexist if the chemical potential of the hard spheres is the same in the dispersion and in the floe phase (i.e., the binodal condition). 

In previous work, we have mainly used the DPM model to investigate the effects of the coefficient of normal restitution and the drag force on the formation of bubbles in fluidized beds (Hoomans et al., 1996 Li and Kuipers, 2003, 2005 Bokkers et al., 2004 Van der Floef et al., 2004), and not so much to obtain information on the constitutive relations that are used in the TFMs. In this section, however, we want to present some recent results from the DPM model on the excess compressibility of the solids phase, which is a key quantity in the constitutive equations as derived from the KTGF (see Section IV.D.). The excess compressibility y can be obtained from the simulation by use of the virial theorem (Allen and Tildesley, 1990). 

This equation may be used to derive the quantum mechanical virial theorem. For this purpose it is necessary to define the kinetic operator... 

The relationship between the exchange potential of DFT and the corresponding energy functional is established through the virial theorem. The two are related via the following relationship derived by Levy and Perdew [23]... 

Exchange identities utilizing the principle of adiabatic connection and coordinate scaling and a generalized Koopmans theorem were derived and the excited-state effective potential was constructed [65]. The differential virial theorem was also derived for a single excited state [66]. 

The theoretical foundations of these rules are, however, rather weak the first one is supposed to result from a formula derived by London for dispersion forces between unlike molecules, the validity of which is actually restricted to distances much larger than r the second one would only be true for molecules acting as rigid spheres. Many authors tried to check the validity of the combination rules by measuring the second virial coefficients of mixtures. It seems that within the experimental accuracy (unfortunately not very high) both rules are roughly verified.24... 

Table 6. Relationship between the fractal dimension dp the exponent for the molar mass dependence of the second virial coefficient and the expected exponent m for the osmotic modulus when the scaling assumptions of Eqs. (93)-(96) are made. The experimental data were derived from the exponents for the second virial coefficient...

This value of kn is actually low by an order of magnitude for dilute suspensions of charged spheres of radius Rg. This is due to the neglect of interchain correlations for c < c in the structure factor used in the derivation of Eqs. (295)-(298). If the repulsive interaction between polyelectrolyte chains dominates, as expected in salt-free solutions, the virial expansion for viscosity may be valid over considerable range of concentrations where the average distance between chains scales as. This virial series may be approxi-... 

When the adsorbent molecides are not independent, we can no longer use the relation (D.2) for the GPF of the system. In this case, we must start from the GPF of the macroscopic system from which we can derive the general form of the BI for any concentration of the adsorbent molecule. The derivation is possible through the McMillan-Mayer theory of solution, but it is long and tedious, even for first-order deviations from an ideal solution. The reason is that, in the general case, the first-order deviations would depend on many second-virial coefficients [the analogue of the quantity B2(T) in Eq. (D.9)]. For each pair of occupancy states, say i and j, there will be a pair potential [/pp(R, i,j), and the corresponding second-virial coefficient... 

---