Grand canonical ensemble
Statistical mechanics |
---|
In
The thermodynamic variables of the grand canonical ensemble are
Basics
In simple terms, the grand canonical ensemble assigns a probability P to each distinct microstate given by the following exponential:
where N is the number of particles in the microstate and E is the total energy of the microstate. k is
The number Ω is known as the grand potential and is constant for the ensemble. However, the probabilities and Ω will vary if different µ, V, T are selected. The grand potential Ω serves two roles: to provide a normalization factor for the probability distribution (the probabilities, over the complete set of microstates, must add up to one); and, many important ensemble averages can be directly calculated from the function Ω(µ, V, T).
In the case where more than one kind of particle is allowed to vary in number, the probability expression generalizes to
where µ1 is the chemical potential for the first kind of particles, N1 is the number of that kind of particle in the microstate, µ2 is the chemical potential for the second kind of particles and so on (s is the number of distinct kinds of particles). However, these particle numbers should be defined carefully (see the note on particle number conservation below).
The distribution of the grand canonical ensemble is called generalized Boltzmann distribution by some authors.[2]
Grand ensembles are apt for use when describing systems such as the
- Note on formulation
- An alternative formulation for the same concept writes the probability as , using the grand partition function rather than the grand potential. The equations in this article (in terms of grand potential) may be restated in terms of the grand partition function by simple mathematical manipulations.
Applicability
The grand canonical ensemble is the ensemble that describes the possible states of an isolated system that is in thermal and chemical equilibrium with a reservoir (the derivation proceeds along lines analogous to the heat bath derivation of the normal
The condition that the system is isolated is necessary in order to ensure it has well-defined thermodynamic quantities and evolution.[1] In practice, however, it is desirable to apply the grand canonical ensemble to describe systems that are in direct contact with the reservoir, since it is that contact that ensures the equilibrium. The use of the grand canonical ensemble in these cases is usually justified either 1) by assuming that the contact is weak, or 2) by incorporating a part of the reservoir connection into the system under analysis, so that the connection's influence on the region of interest is correctly modeled. Alternatively, theoretical approaches can be used to model the influence of the connection, yielding an open statistical ensemble.
Another case in which the grand canonical ensemble appears is when considering a system that is large and thermodynamic (a system that is "in equilibrium with itself"). Even if the exact conditions of the system do not actually allow for variations in energy or particle number, the grand canonical ensemble can be used to simplify calculations of some thermodynamic properties. The reason for this is that various thermodynamic ensembles (microcanonical, canonical) become equivalent in some aspects to the grand canonical ensemble, once the system is very large.[note 1] Of course, for small systems, the different ensembles are no longer equivalent even in the mean. As a result, the grand canonical ensemble can be highly inaccurate when applied to small systems of fixed particle number, such as atomic nuclei.[4]
Properties
- Uniqueness: The grand canonical ensemble is uniquely determined for a given system at given temperature and given chemical potentials, and does not depend on arbitrary choices such as choice of coordinate system (classical mechanics) or basis (quantum mechanics).[1] The grand canonical ensemble is the only ensemble with constant , V, and T that reproduces the fundamental thermodynamic relation.[5]
- Statistical equilibrium (steady state): A grand canonical ensemble does not evolve over time, despite the fact that the underlying system is in constant motion. Indeed, the ensemble is only a function of the conserved quantities of the system (energy and particle numbers).[1]
- Thermal and chemical equilibrium with other systems: Two systems, each described by a grand canonical ensemble of equal temperature and chemical potentials, brought into thermal and chemical contact[note 2] will remain unchanged, and the resulting combined system will be described by a combined grand canonical ensemble of the same temperature and chemical potentials.[1]
- Maximum entropy: For given mechanical parameters (fixed V), the grand canonical ensemble average of the log-probability (also called the "entropy") is the maximum possible for any ensemble (i.e. probability distribution P) with the same , , etc.[1]
- Minimum grand potential: For given mechanical parameters (fixed V) and given values of T, µ1, …, µs, the ensemble average is the lowest possible of any ensemble.[1]
Grand potential, ensemble averages, and exact differentials
The partial derivatives of the function Ω(µ1, …, µs, V, T) give important grand canonical ensemble average quantities:[1][6]
- the averages of numbers of particles
- the average pressure
- the Gibbs entropy
- and the average energy
Exact differential: From the above expressions, it can be seen that the function Ω has the exact differential
First law of thermodynamics: Substituting the above relationship for ⟨E⟩ into the exact differential of Ω, an equation similar to the first law of thermodynamics is found, except with average signs on some of the quantities:[1]
Thermodynamic fluctuations: The variances in energy and particle numbers are[7][8]
Correlations in fluctuations: The covariances of particle numbers and energy are[1]
Example ensembles
The usefulness of the grand canonical ensemble is illustrated in the examples below. In each case the grand potential is calculated on the basis of the relationship
which is required for the microstates' probabilities to add up to 1.
Statistics of noninteracting particles
Bosons and fermions (quantum)
In the special case of a quantum system of many non-interacting particles, the thermodynamics are simple to compute.[9] Since the particles are non-interacting, one can compute a series of single-particle stationary states, each of which represent a separable part that can be included into the total quantum state of the system. For now let us refer to these single-particle stationary states as orbitals (to avoid confusing these "states" with the total many-body state), with the provision that each possible internal particle property (
Since the particles are non-interacting, we may take the viewpoint that each orbital forms a separate thermodynamic system. Thus each orbital is a grand canonical ensemble unto itself, one so simple that its statistics can be immediately derived here. Focusing on just one orbital labelled i, the total energy for a microstate of N particles in this orbital will be Nϵi, where ϵi is the characteristic energy level of that orbital. The grand potential for the orbital is given by one of two forms, depending on whether the orbital is bosonic or fermionic:
- For fermions, the Pauli exclusion principle allows only two microstates for the orbital (occupation of 0 or 1), giving a two-term series
- For indistinguishability of particles, leading to a geometric series:
In each case the value gives the thermodynamic average number of particles on the orbital: the
Indistinguishable classical particles
In classical mechanics it is also possible to consider indistinguishable particles (in fact, indistinguishability is a prerequisite for defining a chemical potential in a consistent manner; all particles of a given kind must be interchangeable[1]). We again consider placing multiple particles of the same kind into the same microstate of single-particle phase space, which we again call an "orbital". However, compared to quantum mechanics, the classical case is complicated by the fact that a microstate in classical mechanics does not refer to a single point in phase space but rather to an extended region in phase space: one microstate contains an infinite number of states, all distinct but of similar character. As a result, when multiple particles are placed into the same orbital, the overall collection of the particles (in the system phase space) does not count as one whole microstate but rather only a fraction of a microstate, because identical states (formed by permutation of identical particles) should not be overcounted. The overcounting correction factor is the factorial of the number of particles.
The statistics in this case take the form of an exponential power series
the value corresponding to Maxwell–Boltzmann statistics.
Ionization of an isolated atom
The grand canonical ensemble can be used to predict whether an atom prefers to be in a neutral state or ionized state. An atom is able to exist in ionized states with more or fewer electrons compared to neutral. As shown below, ionized states may be thermodynamically preferred depending on the environment. Consider a simplified model where the atom can be in a neutral state or in one of two ionized states (a detailed calculation also includes the degeneracy factors of the states[10]):
- charge neutral state, with N0 electrons and energy E0.
- an oxidizedstate (N0 − 1 electrons) with energy E0 + ΔEI + qϕ
- a reduced state (N0 + 1 electrons) with energy E0 − ΔEA − qϕ
Here ΔEI and ΔEA are the atom's
The grand potential in this case is thus determined by
The quantity −qϕ − µ is critical in this case, for determining the balance between the various states. This value is determined by the environment around the atom.
If one of these atoms is placed into a vacuum box, then −qϕ − µ = W, the
At room temperature, this example finds application in
Meaning of chemical potential, generalized "particle number"
In order for a particle number to have an associated chemical potential, it must be conserved during the internal dynamics of the system, and only able to change when the system exchanges particles with an external reservoir.
If the particles can be created out of energy during the dynamics of the system, then an associated µN term must not appear in the probability expression for the grand canonical ensemble. In effect, this is the same as requiring that µ = 0 for that kind of particle. Such is the case for photons in a black cavity, whose number regularly change due to absorption and emission on the cavity walls. (On the other hand, photons in a highly reflective cavity can be conserved and caused to have a nonzero µ.[13])
In some cases the number of particles is not conserved and the N represents a more abstract conserved quantity:
- Chemical reactions: Chemical reactions can convert one type of molecule to another; if reactions occur then the Ni must be defined such that they do not change during the chemical reaction.
- High energy particle physics: Ordinary particles can be spawned out of pure energy, if a corresponding antiparticle is created. If this sort of process is allowed, then neither the number of particles nor antiparticles are conserved. Instead, N = (particle number - antiparticle number) is conserved.[14][note 3] As particle energies increase, there are more possibilities to convert between particle types, and so there are fewer numbers that are truly conserved. At the very highest energies the only conserved numbers are electric charge, weak isospin, and baryon number − lepton number.
On the other hand, in some cases a single kind of particle may have multiple conserved numbers:
- Closed compartments: In a system composed of multiple compartments that share energy but do not share particles, it is possible to set the chemical potentials separately for each compartment. For example, a electron chemical potential.
- Slow equilibration: In some quasi-equilibrium situations it is possible to have two distinct populations of the same kind of particle in the same location, which are each equilibrated internally but not with each other. Though not strictly in equilibrium, it may be useful to name quasi-equilibrium chemical potentials which can differ among the different populations. Examples: (orthohydrogenchemical potentials.
Precise expressions for the ensemble
The precise mathematical expression for statistical ensembles has a distinct form depending on the type of mechanics under consideration (quantum or classical), as the notion of a "microstate" is considerably different. In quantum mechanics, the grand canonical ensemble affords a simple description since diagonalization provides a set of distinct microstates of a system, each with well-defined energy and particle number. The classical mechanical case is more complex as it involves not stationary states but instead an integral over canonical phase space.
Quantum mechanical
A statistical ensemble in quantum mechanics is represented by a density matrix, denoted by . The grand canonical ensemble is the density matrix[citation needed]
where Ĥ is the system's total energy operator (Hamiltonian), N̂1 is the system's total particle number operator for particles of type 1, N̂2 is the total particle number operator for particles of type 2, and so on. exp is the matrix exponential operator. The grand potential Ω is determined by the probability normalization condition that the density matrix has a trace of one, :
Note that for the grand ensemble, the basis states of the operators Ĥ, N̂1, etc. are all states with multiple particles in Fock space, and the density matrix is defined on the same basis. Since the energy and particle numbers are all separately conserved, these operators are mutually commuting.
The grand canonical ensemble can alternatively be written in a simple form using bra–ket notation, since it is possible (given the mutually commuting nature of the energy and particle number operators) to find a complete basis of simultaneous eigenstates |ψi⟩, indexed by i, where Ĥ|ψi⟩ = Ei|ψi⟩, N̂1|ψi⟩ = N1,i|ψi⟩, and so on. Given such an eigenbasis, the grand canonical ensemble is simply
where the sum is over the complete set of states with state i having Ei total energy, N1,i particles of type 1, N2,i particles of type 2, and so on.
Classical mechanical
In classical mechanics, a grand ensemble is instead represented by a
- The number of particles and thus the number of coordinates n varies between the different phase spaces, and,
- it is vital to consider whether permuting similar particles counts as a distinct state or not.
In a system of particles, the number of degrees of freedom n depends on the number of particles in a way that depends on the physical situation. For example, in a three-dimensional gas of monoatoms n = 3N, however in molecular gases there will also be rotational and vibrational degrees of freedom.
The probability density function for the grand canonical ensemble is:
where
- E is the energy of the system, a function of the phase (N1, … Ns, p1, … pn, q1, … qn),
- h is an arbitrary but predetermined constant with the units of energy×time, setting the extent of one microstate and providing correct dimensions to ρ.[note 4]
- C is an overcounting correction factor (see below), a function of N1, … Ns.
Again, the value of Ω is determined by demanding that ρ is a normalized probability density function:
This integral is taken over the entire available phase space for the given numbers of particles.
Overcounting correction
A well-known problem in the statistical mechanics of fluids (gases, liquids, plasmas) is how to treat particles that are similar or identical in nature: should they be regarded as distinguishable or not? In the system's equation of motion each particle is forever tracked as a distinguishable entity, and yet there are also valid states of the system where the positions of each particle have simply been swapped: these states are represented at different places in phase space, yet would seem to be equivalent.
If the permutations of similar particles are regarded to count as distinct states, then the factor C above is simply C = 1. From this point of view, ensembles include every permuted state as a separate microstate. Although appearing benign at first, this leads to a problem of severely non-extensive entropy in the canonical ensemble, known today as the Gibbs paradox. In the grand canonical ensemble a further logical inconsistency occurs: the number of distinguishable permutations depends not only on how many particles are in the system, but also on how many particles are in the reservoir (since the system may exchange particles with a reservoir). In this case the entropy and chemical potential are non-extensive but also badly defined, depending on a parameter (reservoir size) that should be irrelevant.
To solve these issues it is necessary that the exchange of two similar particles (within the system, or between the system and reservoir) must not be regarded as giving a distinct state of the system.[1][note 5] In order to incorporate this fact, integrals are still carried over full phase space but the result is divided by
which is the number of different permutations possible. The division by C neatly corrects the overcounting that occurs in the integral over all phase space.
It is of course possible to include distinguishable types of particles in the grand canonical ensemble—each distinguishable type is tracked by a separate particle counter and chemical potential . As a result, the only consistent way to include "fully distinguishable" particles in the grand canonical ensemble is to consider every possible distinguishable type of those particles, and to track each and every possible type with a separate particle counter and separate chemical potential.
See also
Notes
- ^ To quote Reif, "For purposes of calculating mean values of physical quantities, it makes no noticeable difference whether a macroscopic system is isolated, or in contact with a reservoir with which it can only exchange energy, or in contact with a reservoir with which it can exchange both energy and particles. [...] In some problems where the constraint of a fixed number of particles is cumbersome, one can thus readily circumvent the complication by approximating the actual situation with [...] the grand canonical distribution."
- ^ Thermal and chemical contact means that the systems are able to exchange energy and particles through a connection. The connection must be weak as to not significantly disturb the systems' microstates.
- ^ Of course, very high temperatures are required for significant thermal generation of particle-antiparticle pairs, e.g., of order 109 K for electron-positron creation, and so this process is not a concern for everyday thermodynamics.
- Planck's constantin order to obtain a semiclassical correspondence with quantum mechanics.
- ^ This can be compared to the canonical ensemble where it is optional to consider particles as distinguishable; this only gives N-dependent error in entropy, which is unobservable as long as N is kept constant. In general, however, there is no such freedom: "when the number of particles in a system is to be treated as variable, the average index of probability for phases generically defined corresponds to entropy." (Gibbs).
References
- ^ a b c d e f g h i j k l m Gibbs, Josiah Willard (1902). Elementary Principles in Statistical Mechanics. New York: Charles Scribner's Sons.
- S2CID 118981017.
- ISBN 9780070518001.
- S2CID 119152931.
- S2CID 221978379.
- ^ "5.1 the Gibbs Distribution".
- ^ "Archived copy" (PDF). Archived from the original (PDF) on 2013-10-19. Retrieved 2013-05-02.
{{cite web}}
: CS1 maint: archived copy as title (link) - ^ "Handout 9. NPT and Grand Canonical Ensembles" (PDF). micro.stanford.edu. 26 January 2011. Retrieved 4 November 2023.
- ISBN 9788120327825.
- ^ ISBN 0198517408.
- .
- ^ "2. Semiconductor Doping Technology".
- .
- S2CID 18182534.