In statistical mechanics, the properties of a system of noninteracting particles are described using the partition function. For N noninteracting particles, the system is described by the Hamiltonian and the partition function can be calculated as From the partition function, one can calculate the Helmholtz free energy and, from that, all the thermodynamic properties of the system, like the entropy, the internal energy, the chemical potential etc. When the particles of the system interact, an exact calculation of the partition function is usually not possible. For low density, the interactions can be approximated with a sum of two-particle potentials: For this interaction potential, the partition function can be written as and the free energy is where Q is the configuration integral:
The configuration integral cannot be calculated analytically for a general pair potential . One way to calculate the potential approximately is to use the Mayer cluster expansion. This expansion is based on the observation that the exponential in the equation for can be written as a product of the form Next, define the Mayer function by. After substitution, the equation for the configuration integral becomes: The calculation of the product in the above equation leads into a series of terms; the first is equal to one, the second term is equal to the sum over i and j of the terms, and the process continues until all the higher order terms are calculated. Each term must appear only once. With this expansion it is possible to find terms of different order, in terms of the number of particles that are involved. The first term is the non-interaction term, the second term corresponds to the two-particle interactions, the third to the two-particle interactions amongst 4 particles, and so on. This physical interpretation is the reason this expansion is called the cluster expansion: the sum can be rearranged so that each term represents the interactions within clusters of a certain number of particles. Substituting the expansion of the product back into the expression for the configuration integral results in a series expansion for : Substituting in the equation for the free energy, it is possible to derive the equation of state for the system of interacting particles. The equation will have the form which is known as the virial equation, and the components are the virial coefficients. Each of the virial coefficients corresponds to one term from the cluster expansion. Keeping only the two-particle interaction term, it can be shown that the cluster expansion, with some approximations, gives the Van der Waals equation. This can be applied further to mixtures of gases and liquid solutions.
