Recent from talks
Nothing was collected or created yet.
Virial coefficient
View on WikipediaVirial coefficients appear as coefficients in the virial expansion of the pressure of a many-particle system in powers of the density, providing systematic corrections to the ideal gas law. They are characteristic of the interaction potential between the particles and in general depend on the temperature. The second virial coefficient depends only on the pair interaction between the particles, the third () depends on 2- and non-additive 3-body interactions, and so on.
Derivation
[edit]The first step in obtaining a closed expression for virial coefficients is a cluster expansion[1] of the grand canonical partition function
Here is the pressure, is the volume of the vessel containing the particles, is the Boltzmann constant, is the absolute temperature, is the fugacity, with the chemical potential. The quantity is the canonical partition function of a subsystem of particles:
Here is the Hamiltonian (energy operator) of a subsystem of particles. The Hamiltonian is a sum of the kinetic energies of the particles and the total -particle potential energy (interaction energy). The latter includes pair interactions and possibly 3-body and higher-body interactions. The grand partition function can be expanded in a sum of contributions from one-body, two-body, etc. clusters. The virial expansion is obtained from this expansion by observing that equals . In this manner one derives
- .
These are quantum-statistical expressions containing kinetic energies. Note that the one-particle partition function contains only a kinetic energy term. In the classical limit the kinetic energy operators commute with the potential operators and the kinetic energies in numerator and denominator cancel mutually. The trace (tr) becomes an integral over the configuration space. It follows that classical virial coefficients depend on the interactions between the particles only and are given as integrals over the particle coordinates.
The derivation of higher than virial coefficients becomes quickly a complex combinatorial problem. Making the classical approximation and neglecting non-additive interactions (if present), the combinatorics can be handled graphically as first shown by Joseph E. Mayer and Maria Goeppert-Mayer.[2]
They introduced what is now known as the Mayer function:
and wrote the cluster expansion in terms of these functions. Here is the interaction potential between particle 1 and 2 (which are assumed to be identical particles).
Definition in terms of graphs
[edit]The virial coefficients are related to the irreducible Mayer cluster integrals through
The latter are concisely defined in terms of graphs.
The rule for turning these graphs into integrals is as follows:
- Take a graph and label its white vertex by and the remaining black vertices with .
- Associate a labelled coordinate k to each of the vertices, representing the continuous degrees of freedom associated with that particle. The coordinate 0 is reserved for the white vertex
- With each bond linking two vertices associate the Mayer f-function corresponding to the interparticle potential
- Integrate over all coordinates assigned to the black vertices
- Multiply the end result with the symmetry number of the graph, defined as the inverse of the number of permutations of the black labelled vertices that leave the graph topologically invariant.
The first two cluster integrals are
The expression of the second virial coefficient is thus:
where particle 2 was assumed to define the origin (). This classical expression for the second virial coefficient was first derived by Leonard Ornstein in his 1908 Leiden University Ph.D. thesis.
See also
[edit]- Boyle temperature – temperature at which the second virial coefficient vanishes
- Excess property
- Compressibility factor
References
[edit]- ^ Hill, T. L. (1960). Introduction to Statistical Thermodynamics. Addison-Wesley. ISBN 9780201028409.
{{cite book}}: ISBN / Date incompatibility (help) - ^ Mayer, J. E.; Goeppert-Mayer, M. (1940). Statistical Mechanics. New York: Wiley.
Further reading
[edit]- Dymond, J. H.; Smith, E. B. (1980). The Virial Coefficients of Pure Gases and Mixtures: a Critical Compilation. Oxford: Clarendon. ISBN 0198553617.
- Hansen, J. P.; McDonald, I. R. (1986). The Theory of Simple Liquids (2nd ed.). London: Academic Press. ISBN 012323851X.
- http://scitation.aip.org/content/aip/journal/jcp/50/10/10.1063/1.1670902
- http://scitation.aip.org/content/aip/journal/jcp/50/11/10.1063/1.1670994
- Reid, C. R., Prausnitz, J. M., Poling B. E., Properties of gases and liquids, IV edition, Mc Graw-Hill, 1987
Virial coefficient
View on GrokipediaOverview
Definition
Virial coefficients are temperature-dependent parameters that quantify the deviations from ideal gas behavior in dilute gases arising from intermolecular interactions.[2] These coefficients appear in the virial expansion of thermodynamic properties, providing a systematic way to account for non-ideal effects at low densities where higher-order terms become negligible.[4] The general form of the virial equation of state for pressure is where is the number density, is Boltzmann's constant, and is the temperature.[2] By convention, the first virial coefficient is set to , representing the ideal gas contribution, while the higher-order coefficients (for ) incorporate the cumulative effects of interactions involving increasing numbers of particles.[4] The term "virial" originates from the virial theorem in classical mechanics, formulated by Rudolf Clausius in 1870.Virial Expansion
The virial expansion represents the equation of state of a non-ideal gas as a power series in density, enabling the quantification of intermolecular interactions beyond the ideal gas limit. The compressibility factor , which measures deviations from ideal behavior and is given by (where is pressure, is volume, is the number of particles, is Boltzmann's constant, and is temperature), takes the form Here, denotes the number density, and , , , and higher coefficients are the virial coefficients, which depend solely on temperature.[3] This density-based expansion is particularly advantageous for analyzing real gases at low to moderate densities, where it captures the effects of finite particle size and attractive forces that cause to deviate from unity.[5] Alternative formulations of the virial expansion exist to suit different applications. In terms of molar volume , the expansion becomes , where is the gas constant and the primed coefficients relate to the unprimed ones by scaling factors involving Avogadro's number.[6] Another variant expands the pressure in powers of fugacity, useful for phase equilibria and mixture properties at low pressures.[3] These forms, originally proposed by Kamerlingh Onnes in 1901, provide flexibility in fitting experimental data for various gas systems.[7] Each term in the virial series carries a physical interpretation tied to the order of particle interactions: the leading term of 1 reflects non-interacting particles, accounts for pairwise correlations (e.g., exclusions and attractions between two molecules), incorporates three-body effects, and higher terms represent increasingly complex multi-particle correlations.[3] The expansion is valid primarily at low densities, where is small enough that higher-order terms diminish rapidly, ensuring convergence and approximating the behavior of dilute non-ideal gases effectively.[5] This makes it a cornerstone tool for understanding and modeling real gas thermodynamics in regimes inaccessible to the ideal gas law.[6] The coefficients , , and so on correspond to the second, third, and higher virial coefficients introduced in the definition section.Historical Development
Virial Theorem
The virial theorem was introduced by Rudolf Clausius in 1870 as a generalization of the center-of-mass theorem, originally formulated to apply to mechanical systems and their relation to heat processes.[8] Clausius derived it from the equations of motion for a system of particles, emphasizing its utility in understanding equilibrium states where forces and positions balance dynamically.[8] The theorem states that for a stable system of discrete particles bound by internal forces, twice the time-averaged total kinetic energy equals the time average of the sum over all particles of the scalar product of their position vectors and the forces acting on them: where denotes the total kinetic energy, is the position vector of the -th particle, and is the total force on that particle.[9] This relation arises from considering the time derivative of the moment of inertia-like quantity , where is the momentum, and averaging over a sufficiently long period where boundary terms vanish.[9] For bound systems, such as those in periodic orbits or confined within finite volumes, the time-averaged form holds because the system's motion is ergodic or recurrent, allowing the average to capture the equilibrium distribution of configurations.[9] In these cases, the theorem provides insight into the balance between kinetic and potential energies without requiring explicit solutions to the equations of motion. This mechanical principle extends to gases in thermal equilibrium, where the virial relates the average kinetic energy of molecules to the pressure they exert on container walls, forming a basis for kinetic theory derivations.[9]Introduction of Coefficients
The virial theorem, established in the classical mechanics of the 19th century, laid the groundwork for analyzing the pressure-volume relations in gases, which later influenced the development of expansions for real gases beyond the ideal gas law. In the early 20th century, Heike Kamerlingh Onnes introduced the virial expansion empirically in 1901 to describe the equation of state for real gases at low densities, expressing the compressibility factor as a power series in density with coefficients that capture deviations from ideality. This approach generalized the ideal gas law by incorporating higher-order terms, allowing for better fitting of experimental data on gases like oxygen and nitrogen under varying pressures and temperatures.[10] A significant theoretical advancement occurred in 1932 when George E. Uhlenbeck and Leonard Gropper formalized the virial expansion within quantum statistical mechanics, deriving expressions for the coefficients in non-ideal Einstein-Bose and Fermi-Dirac gases and highlighting quantum corrections to classical behavior.[11] Building on this, Joseph E. Mayer and collaborators provided a rigorous justification in the 1930s and 1940s by linking the virial coefficients to intermolecular potentials through statistical mechanical derivations, demonstrating how the coefficients arise from integrals over particle interactions.[12] Following World War II, computations of virial coefficients increasingly incorporated quantum effects, particularly for low-temperature gases like helium, enabling more accurate predictions of thermodynamic properties through semiclassical approximations and direct quantum mechanical evaluations.[13]Theoretical Foundations
Derivation from Partition Function
In classical statistical mechanics, the virial coefficients are derived from the canonical ensemble partition function for a system of indistinguishable particles interacting via a potential energy . The partition function is expressed as where , is Planck's constant, is the particle mass, and the integrals are over momenta and positions in volume .[14][15] The momentum integrals factorize and yield , resulting in with the thermal de Broglie wavelength and the reduced configurational integral . For non-interacting particles, , recovering the ideal gas partition function.[14][15] The pressure follows from the thermodynamic relation Applying Stirling's approximation for large , , so At low densities, interactions cause to deviate from 1, and the cluster expansion of in powers of the density generates the virial series. Specifically, the reduced configurational integral admits an expansion leading to the virial equation of state with virial coefficients for and , where the are cluster coefficients from irreducible integrals over the potential . This expansion, first systematically developed by Joseph E. Mayer and Maria Goeppert Mayer, captures deviations from ideal gas behavior due to pairwise and higher-order interactions.[14][16] The cluster expansion provides a computational framework for evaluating the , but the virial form directly stems from the density-dependent logarithmic structure of .[14]Cluster Expansion Approach
The Ursell-Mayer cluster expansion represents a foundational method in statistical mechanics for deriving the virial coefficients through a perturbative treatment of particle interactions in the partition function. Originally developed by H. D. Ursell in 1927 for general many-body systems and extended by J. E. Mayer in the 1930s to classical gases, the approach decomposes the logarithm of the grand partition function into contributions from irreducible clusters of particles. Specifically, , where is the grand partition function, is the fugacity (with , the chemical potential, and the thermal wavelength), and denotes the reduced -particle cluster integral defined as , or equivalently in Mayer form without the Ursell function . Here, sums all connected products of the Mayer f-functions over irreducible diagrams involving particles, ensuring the expansion captures only topologically connected interaction configurations.[17] The connection to virial coefficients arises by expressing the pressure and density in terms of the fugacity series and inverting to obtain the virial expansion , where is the density. The coefficients relate directly to the cluster integrals via in standard reduced notation, highlighting how each emerges from the -particle irreducible clusters. This algebraic structure allows practical computation of virial coefficients by evaluating finite sums of diagrams, starting from the lowest-order terms.[17][18] The expansion's utility stems from its convergence properties in dilute systems, where the density is low and interactions are weak perturbations around the ideal gas limit. Higher-order cluster terms scale as , becoming negligible as , which justifies truncation at low orders (typically up to third or fourth virial coefficients) for gases near ideal behavior; rigorous bounds on the remainder ensure analytic continuation beyond the radius of convergence under suitable conditions on the potential. For light gases such as helium or hydrogen, quantum mechanical effects necessitate corrections to the classical cluster integrals, incorporating Slater sums or path-integral representations to account for diffraction and exchange statistics at low temperatures, where these modifications can alter virial coefficients by up to 10-20% for the second term.[19][20]Graphical Methods
Mayer Cluster Integrals
The Mayer f-bond formalism provides a graphical representation for the cluster integrals in the virial expansion of classical many-body systems, facilitating the systematic evaluation of higher-order terms. Central to this approach is the Mayer f-function, defined for a pair of particles and as , where is the inverse temperature, is Boltzmann's constant, is the temperature, is the pairwise interaction potential, and is the distance between the particles. This function captures the deviation from ideal non-interacting behavior, being zero for large separations where and negative (or positive depending on the potential) in regions of significant interaction.[21] The l-th order cluster integral is then expressed as where the integral is over the coordinates of particles in volume , and the sum runs over all connected labeled graphs formed by the l particles with bonds corresponding to the functions in each graph. This formulation arises from expanding the configurational integral in the canonical partition function using products of the f-functions, which encode pairwise correlations; these cluster integrals enter the activity expansion and are related to the virial coefficients via series inversion between activity and density. The factor accounts for the indistinguishability of particles, and the ensures volume independence.[21][22] A key topological reduction in this formalism restricts the sum to connected graphs only, excluding disconnected configurations. Disconnected graphs factorize into products of integrals over their separate connected components, which are already incorporated into lower-order cluster integrals in the overall expansion of the partition function; including them would lead to overcounting. This reduction simplifies computations by focusing on irreducible, connected structures that cannot be partitioned into independent subclusters. For direct computation of virial coefficients, diagrams without articulation points (irreducible clusters) are often used.[21][23]Diagram Interpretation
In Mayer diagrams, the topological structure encodes the nature of particle interactions within the virial expansion, where vertices represent particles and edges denote Mayer f-bonds, defined as , capturing deviations from ideal gas behavior due to pairwise potentials. Articulation points, or vertices whose removal disconnects the graph, signify locations where multiple sub-clusters converge, thereby representing multi-body correlations that cannot be decomposed into simpler pairwise interactions without altering the overall connectivity. Similarly, bridges—edges whose removal severs the diagram into disconnected components—highlight critical interaction pathways that enforce higher-order correlations among particles, ensuring the diagram's contribution to the virial coefficient reflects irreducible collective effects rather than separable parts.[24][25] For the third virial coefficient, distinct diagram topologies illustrate varying interaction geometries: chain diagrams, consisting of linear sequences of bonds (e.g., particles 1–2–3 connected sequentially), depict acyclic, tree-like multi-body interactions that emphasize sequential correlations along a path. In contrast, ring diagrams feature closed loops (e.g., particles 1–2–3 with bonds forming a triangle), capturing cyclic interactions that introduce feedback loops in particle correlations, contributing to a more compact representation of three-particle effects in the coefficient. These topological differences determine the symmetry and integrability of the corresponding cluster integrals, with rings often yielding higher symmetry factors.[24][26] The classical Mayer diagram approach assumes pairwise additivity and Boltzmann statistics, imposing limitations when applied to quantum gases, where exchange effects and quantum statistics (e.g., Bose or Fermi) disrupt the validity of the expansion beyond dilute regimes, particularly near condensation or degeneracy. Additionally, the standard formulation neglects explicit three-body potentials, requiring extensions to nonadditive interactions for accurate descriptions of systems like quantum fluids or molecular gases with significant many-body forces.[27][28]Specific Coefficients
Second Virial Coefficient
The second virial coefficient, denoted , quantifies the leading-order deviation from ideal gas behavior due to pairwise molecular interactions in a dilute gas.[29] It arises in the virial expansion of the equation of state as the coefficient of the inverse density term, capturing how intermolecular forces modify the pressure relative to the ideal case.[30] Physically, measures the net effect of attractive and repulsive pair interactions between molecules. For systems dominated by repulsive forces, such as hard spheres of diameter , is positive and equals , reflecting the excluded volume per pair that reduces the effective free space available to the gas.[30] In contrast, when attractive interactions prevail, as in gases with potential wells, becomes negative, indicating enhanced molecular clustering that increases pressure deviations from ideality. The explicit expression for in terms of the isotropic pair potential is derived from the cluster expansion in statistical mechanics: where , is Boltzmann's constant, is temperature, and the integral is over all space.[29] For spherically symmetric potentials, this simplifies to The integrand is the Mayer -function for pairs, which is negative in repulsive regions and positive in attractive regions.[29] The value of depends strongly on temperature, transitioning from negative at low (attractive dominance) to positive at high (repulsive dominance). The Boyle temperature is defined as the unique temperature where , at which the gas exhibits ideal-like behavior over a moderate pressure range due to the balance of attractive and repulsive forces. For typical gases, scales with the depth of the potential well, often exceeding room temperature.Third and Higher Coefficients
The third virial coefficient (also denoted ) accounts for three-body interactions in the virial expansion, extending beyond pairwise additivity to capture correlated effects among triplets of particles. It is given by an integral over configurations involving two Mayer -bonds, where the -function is defined as with and the pair potential; these configurations correspond to the irreducible triangle diagram (all three particles mutually connected) and the chain diagram (three particles in a linear arrangement with bonds between consecutive pairs).[29] The triangle term emphasizes cyclic three-body correlations, while the chain term highlights sequential interactions, both essential for accurate modeling of non-ideal gas behavior at moderate densities.[29] For higher orders, the th virial coefficient follows a general combinatorial structure derived from the Mayer cluster expansion: , where the sum is over all one-particle irreducible cluster integrals with bonds connecting labeled points, and each , with the symmetry factor of the graph.[29] This pattern reflects the summation over all connected graphs with bonds, ensuring the expansion captures multi-particle correlations without overcounting due to the factorial prefactor and bond weighting.[29] Computing third and higher virial coefficients presents significant challenges owing to the exponential growth in the number and complexity of irreducible diagrams; for instance, while involves only two diagrams, requires three, ten, and 56, with the topological diversity complicating analytical evaluation and numerical integration. To address the limited convergence radius of the virial series, approximation methods such as Padé resummation are applied, which construct rational functions from partial sums of coefficients to analytically continue the expansion and estimate thermodynamic properties at higher densities.[31] These techniques have proven effective for model systems like hard spheres, where exact low-order coefficients enable reliable extrapolation.[31]Applications
Real Gas Equation of State
The virial equation of state provides a perturbative correction to the ideal gas law for describing the behavior of real gases, particularly at low to moderate densities where intermolecular interactions become significant. It expresses the compressibility factor , where is pressure, is molar volume, is the gas constant, and is temperature, as a power series in inverse molar volume: Here, , , and higher-order coefficients capture pairwise, three-body, and multi-body interactions, respectively, with the leading term recovering the ideal gas limit () as . This form was first proposed empirically by Heike Kamerlingh Onnes in 1901 to fit experimental data for gases like helium and hydrogen, demonstrating its utility in quantifying deviations from ideality. The equation is derived from statistical mechanics via cluster expansions, ensuring thermodynamic consistency for dilute systems. A practical application involves approximating the van der Waals equation of state using the first two virial coefficients, where the excluded volume and attractive interactions correspond roughly to , with and as van der Waals constants; this truncation accurately reproduces van der Waals behavior at low densities. For critical point analysis, truncating the virial series to the second or third order enables estimation of the critical compressibility factor , aligning with experimental values for simple gases and highlighting the role of virial coefficients in phase transition phenomena near criticality. Such approximations are particularly insightful for spherical molecules, where higher terms become negligible below the Boyle temperature. Despite its precision at low densities, the virial equation exhibits limitations at high densities, where the power series diverges due to the accumulation of higher-order terms, failing to capture liquid-like behavior or phase equilibria accurately. In such regimes, cubic equations of state like the Peng-Robinson model are preferred, as they provide a closed-form expression valid across a broader range of conditions, including supercritical states, while still incorporating virial-like corrections in their low-density limits. This transition underscores the virial form's role as a foundational tool for benchmarking more comprehensive models in thermodynamic simulations.Binary Mixtures
In binary mixtures, the virial expansion generalizes the pure-component form to account for interactions between like and unlike molecules, enabling the description of non-ideal behavior in multicomponent gases. The pressure is expressed as , where are the number densities of components , are the second virial coefficients for pure components, and (with ) are the cross virial coefficients capturing unlike-pair interactions.[32] This form arises from the cluster expansion in statistical mechanics, where the cross terms ensure the equation of state reflects composition-dependent deviations from ideality.[33] The cross second virial coefficient for a binary mixture of components 1 and 2 is given by the integral over the Mayer function for unlike pairs:where , is the pairwise interaction potential between unlike molecules separated by , and the integral is over all space. For central potentials with spherical symmetry, this simplifies to
This expression parallels the pure-component virial but uses the unlike-pair potential, highlighting the role of intermolecular forces in mixture properties.[34] These cross coefficients are essential for modeling phase equilibria in binary gas mixtures, such as those involving nitrogen and oxygen in air, where they contribute to accurate fugacity calculations and vapor-liquid equilibrium predictions at elevated pressures. For instance, virial expansions with cross terms have been applied to hydrocarbon mixtures like methane-ethane to compute thermodynamic properties relevant to natural gas processing. To estimate without direct computation from potentials, combining rules are employed for interaction parameters, such as the Lorentz-Berthelot rules for Lennard-Jones potentials: and , which parameterize from pure-component data.[35][36][37]
Experimental Aspects
Measurement Techniques
The primary experimental approach for determining virial coefficients involves pressure-volume-temperature (P-V-T) measurements, particularly at low densities where higher-order terms in the virial expansion become negligible. In this method, the compressibility factor is measured for a gas under controlled isothermal conditions, and data are fitted to the virial equation , allowing extraction of the second virial coefficient from the linear term and the third virial coefficient from quadratic fits at slightly higher densities. A widely adopted technique is the Burnett method, where the gas undergoes successive isothermal expansions between two fixed volumes, with pressures recorded after each expansion to compute density ratios without direct volume measurement; this yields precise and values with uncertainties often below 0.1% for noble gases like helium and argon.[38][39] Acoustic methods provide an indirect yet highly accurate means to infer virial coefficients through measurements of the speed of sound in dilute gases, leveraging thermodynamic relations that connect acoustic properties to the equation of state. Using spherical or cylindrical acoustic resonators, the speed of sound is determined from resonance frequencies at varying pressures and temperatures, typically in the range 100–500 K and up to several MPa; the data are analyzed via the acoustic virial expansion , where and are acoustic virial coefficients related to the thermodynamic and by factors involving the ideal-gas heat capacity ratio . This approach excels for real gases like hydrogen and nitrogen, offering resolutions better than 0.01% and minimal adsorption errors compared to direct P-V-T setups.[40][41] Modern techniques, such as laser interferometry and molecular beam scattering, enable determination of virial coefficients by probing intermolecular potentials at the molecular level, often yielding data that can be compared to theoretical predictions from statistical mechanics. Laser interferometry measures refractivity virial coefficients by detecting phase shifts in laser light passing through gas samples of varying densities in a differential interferometer setup; the refractive index expansion provides and analogs via the Lorentz-Lorenz relation, with applications to atomic gases like neon and argon achieving precisions of 0.5% or better. Complementarily, molecular beam scattering experiments collide crossed molecular beams to map differential cross-sections, from which pair potentials are inverted and integrated to compute ; this method has been pivotal for systems like argon, validating potential models with experimental scattering data at collision energies up to 100 meV. These approaches are particularly valuable for validating theoretical expressions, as potential-derived virials closely match P-V-T results for simple gases.[42][43]Tabulated Values
One of the primary resources for compiled virial coefficient data is the series of critical compilations by J. H. Dymond and E. B. Smith, beginning with their 1969 work and expanded in the 1980 edition, which systematically tabulates second and third virial coefficients for numerous pure gases based on experimental measurements up to that period.[44] This was followed by a NIST sequel in 2002, "Virial Coefficients of Pure Gases," edited by Dymond, Smith, and others, incorporating additional data and refinements for pure gases through the late 1990s.[45] More recent updates appear in NIST's 2021 compilation for mixtures, building on the pure gas foundations, while the REFPROP database integrates virial data for practical thermophysical calculations up to the 2020s.[46][47] Representative examples from these compilations illustrate typical values for the second virial coefficient . For nitrogen () at 300 K, cm³/mol, reflecting attractive intermolecular forces dominant at this temperature.[48] For noble gases, trends show quantum effects influencing ; for instance, helium exhibits positive deviations at low temperatures due to exchange effects, while argon at 300 K has cm³/mol, with path-integral calculations confirming quantum corrections up to several percent.[49]| Gas | Temperature (K) | (cm³/mol) | Source Compilation |
|---|---|---|---|
| Nitrogen | 300 | -4.2 | Dymond & Smith (1980) |
| Argon | 300 | -21.7 | NIST (2002) |
| Helium | 300 | 11.9 | NIST (2002) |
