Fact-checked by Grok 2 weeks ago

Quantum statistical mechanics

Quantum statistical mechanics is the branch of physics that integrates the principles of quantum mechanics with statistical mechanics to describe the thermodynamic properties and collective behavior of systems composed of many indistinguishable particles, such as atoms or electrons, where quantum effects like wave-particle duality and superposition play a dominant role.^[1] Unlike classical statistical mechanics, it employs quantum statistics to account for particle indistinguishability, leading to distinct distribution functions for bosons and fermions that govern occupation numbers in energy states.^[2] The field originated in the mid-1920s amid the formulation of quantum mechanics, with Satyendra Nath Bose's 1924 paper deriving Planck's law for blackbody radiation using a novel counting method for indistinguishable photons, which Einstein extended in 1925 to massive particles, establishing Bose-Einstein statistics.^[3] Independently, in 1926, Enrico Fermi and Paul Dirac developed Fermi-Dirac statistics for particles with half-integer spin, incorporating the Pauli exclusion principle to explain electron behavior in atoms and metals.^[4] These foundational works resolved inconsistencies in classical theories, such as the ultraviolet catastrophe in radiation and the failure to explain atomic specific heats at low temperatures.^[2] Central to quantum statistical mechanics are the quantum ensembles, including the microcanonical (fixed energy), canonical (fixed temperature), and grand canonical (fixed temperature and chemical potential), formulated using the density operator \rho whose expectation values are computed via traces over Hilbert space, \langle A \rangle = \mathrm{Tr}(\rho A).^[5] The partition function Z = \mathrm{Tr}(e^{-\beta H}), where H is the Hamiltonian and \beta = 1/(k_B T), enables derivation of thermodynamic quantities like free energy and entropy.^[2] Notable applications include the prediction and realization of Bose-Einstein condensation (BEC) in dilute atomic gases, first achieved experimentally in 1995, which exhibits macroscopic quantum coherence and superfluidity.^[4] For fermions, it underpins the theory of degenerate electron gases in white dwarfs and metals, explaining electrical conductivity and magnetism, while in interacting systems, it supports models of superconductivity via BCS theory.^[6] Quantum statistical mechanics also extends to quantum field theory for relativistic particles and non-equilibrium processes, influencing fields from condensed matter physics to cosmology.^[7]

Foundations of Quantum Statistics

Relation to classical statistical mechanics

Quantum statistical mechanics developed as a natural extension of classical statistical mechanics to resolve discrepancies between classical predictions and experimental observations of atomic and subatomic phenomena, such as blackbody radiation and specific heats of solids.^[8] In the classical framework, Ludwig Boltzmann introduced the ergodic hypothesis around 1871–1877, asserting that an isolated system in thermal equilibrium will, over sufficiently long times, traverse all accessible microstates in phase space with equal probability, thereby equating time averages of observables to ensemble averages.^[9] This dynamical assumption underpinned Boltzmann's derivation of the entropy as S = k \ln W, where W is the number of microstates corresponding to a macrostate, linking microscopic irreversibility to macroscopic thermodynamics.^[10] Building on this, Josiah Willard Gibbs formalized the ensemble theory in 1902, representing equilibrium states via probability densities over the entire phase space of many-particle systems to compute thermodynamic potentials systematically.^[11] A fundamental difference lies in the descriptive framework: classical statistical mechanics operates in a continuous phase space spanned by position \mathbf{q} and momentum \mathbf{p} coordinates for each particle, where deterministic trajectories evolve under Hamilton's equations, and fluctuations are treated probabilistically due to ignorance of initial conditions.^[12] In contrast, quantum statistical mechanics employs the Hilbert space of wavefunctions or state vectors, where particles lack simultaneous definite positions and momenta due to the Heisenberg uncertainty principle, and states can superpose, leading to interference and non-classical correlations. Observables become self-adjoint operators, with measurements yielding eigenvalues probabilistically according to Born's rule, replacing classical point particles with delocalized quantum entities whose collective behavior requires accounting for intrinsic quantum fluctuations beyond mere ensemble averaging.^[13] The partition function exemplifies this transition, serving as a generating function for thermodynamic properties in both regimes. Classically, for an N-particle system, it is

Z = \frac{1}{N! h^{3N}} \int e^{-\beta H(\mathbf{q}, \mathbf{p})} \, d^{3N}\mathbf{q} \, d^{3N}\mathbf{p},

where \beta = 1/(k_B T), H is the classical Hamiltonian, h is Planck's constant introduced for dimensional consistency and to handle indistinguishability via the $1/N! factor, and the integral averages the Boltzmann weight over phase space volume element d\Gamma.^[13] Quantum mechanically, the partition function becomes the trace over the Hilbert space,

Z = \mathrm{Tr} \left( e^{-\beta \hat{H}} \right) = \sum_n e^{-\beta E_n},

summing over the discrete energy eigenstates |n\rangle of the quantum Hamiltonian \hat{H}, inherently incorporating quantum level spacing and statistics for identical particles without ad hoc corrections. In the semiclassical limit of high temperatures or large quantum numbers, the quantum Z approaches the classical form via the Weyl correspondence, validating the continuity between the theories.^[13] The evolution of statistical descriptions also parallels closely: the classical Liouville equation, \partial_t \rho + \{\rho, H\}_{\mathrm{PB}} = 0, dictates the incompressible flow of the phase-space probability density \rho(\mathbf{q}, \mathbf{p}, t) under the Poisson bracket \{\cdot, \cdot\}_{\mathrm{PB}}, conserving total probability and phase volume.^[12] Quantum mechanically, John von Neumann introduced the analogous equation in 1929 for the density operator \hat{\rho}(t), i\hbar \partial_t \hat{\rho} = [\hat{H}, \hat{\rho}], where the commutator replaces the Poisson bracket, ensuring unitary evolution that preserves the trace (total probability) and von Neumann entropy, thus extending Liouville's theorem to the quantum domain. The density operator \hat{\rho} generalizes the classical distribution by accommodating pure states, mixed states from ensembles, and decoherence effects.^[14]

Quantum description of many-particle systems

In quantum mechanics, the configuration of a many-particle system is described in a composite Hilbert space formed by the tensor product of the individual single-particle Hilbert spaces. For a system of N particles, the total Hilbert space \mathcal{H} is given by \mathcal{H} = \bigotimes_{i=1}^N \mathcal{H}_i, where \mathcal{H}_i denotes the Hilbert space associated with the i-th particle, capturing its degrees of freedom such as position and spin. This construction allows the state of the entire system to be represented as a vector in a space whose dimension grows exponentially with N, reflecting the vast number of possible configurations.^[15] Observables in this framework are represented by Hermitian operators acting on \mathcal{H}, ensuring that their eigenvalues correspond to measurable real values. Key examples include the position operator \hat{\mathbf{r}}_i and momentum operator \hat{\mathbf{p}}_i for each particle i, which satisfy the canonical commutation relations [\hat{r}_{i,\alpha}, \hat{p}_{i,\beta}] = i\hbar \delta_{i j} \delta_{\alpha \beta}, where \alpha, \beta label spatial components and the indices i, j distinguish particles. The total Hamiltonian operator \hat{H}, which dictates the dynamics via the Schrödinger equation i\hbar \frac{\partial}{\partial t} |\psi\rangle = \hat{H} |\psi\rangle, typically takes the form

\hat{H} = \sum_{i=1}^N \frac{\hat{\mathbf{p}}_i^2}{2m} + V(\{\hat{\mathbf{r}}_i\}),

combining the kinetic energy terms for each particle with a potential energy operator V that accounts for interactions, such as pairwise potentials between particles or external fields. This operator is Hermitian, guaranteeing real energy eigenvalues.^[16]^[17] Pure states of the system are normalized state vectors |\psi\rangle \in \mathcal{H} with \langle \psi | \psi \rangle = 1, which can be expressed as linear superpositions of basis states according to the superposition principle, a cornerstone of quantum mechanics that enables interference effects absent in classical descriptions. For distinguishable particles, such as those in different internal states or trapped separately, the full tensor product structure permits product states like |\psi\rangle = \bigotimes_i |\psi_i\rangle, where each factor describes an individual particle. In contrast, for indistinguishable particles, the superposition principle requires states to be constructed as appropriate linear combinations that account for particle identity, leading to fundamentally different quantum behaviors compared to the distinguishable case, though the specific projection onto symmetric or antisymmetric subspaces is addressed in treatments of quantum statistics. The non-commutativity embodied in the position-momentum relations introduces the Heisenberg uncertainty principle, \Delta r \Delta p \geq \hbar/2, which imposes fundamental limits on simultaneous measurements and underpins the probabilistic nature essential to statistical mechanics in quantum systems.^[15]^[18]

Density Matrix Formalism

Definition and properties of the density operator

In quantum statistical mechanics, the density operator, often denoted as \rho, serves as the central mathematical tool for describing the state of a quantum system in a probabilistic or ensemble sense, extending beyond the pure state vectors of standard quantum mechanics. Introduced by John von Neumann in his foundational work on the probabilistic structure of quantum theory, the density operator encapsulates both pure and mixed states within a unified framework.^[19] For a pure state represented by a normalized state vector |\psi\rangle in the system's Hilbert space, the density operator is defined as the outer product \rho = |\psi\rangle\langle\psi|. This operator projects onto the subspace spanned by |\psi\rangle and satisfies the normalization condition \operatorname{Tr}(\rho) = 1. In the case of a mixed state, which arises when the system is known to be in one of several pure states |\psi_i\rangle with respective probabilities p_i (satisfying \sum_i p_i = 1 and p_i \geq 0), the density operator takes the general form \rho = \sum_i p_i |\psi_i\rangle\langle\psi_i|. This construction interprets the mixed state as an average over an ensemble of pure states, reflecting incomplete knowledge or statistical ignorance about the precise preparation of the system.^[19]^[20] The density operator possesses several intrinsic properties that ensure its utility in statistical descriptions. It is Hermitian (\rho^\dagger = \rho), guaranteeing real eigenvalues, and positive semi-definite, meaning all eigenvalues \lambda_k \geq 0, with the eigenvalues themselves interpretable as probabilities (corresponding to the p_i in a suitable basis). The trace condition \operatorname{Tr}(\rho) = 1 enforces normalization, independent of the choice of basis, underscoring the basis-independent nature of \rho. For equilibrium states, \rho is often diagonal in the energy eigenbasis, simplifying computations for thermal ensembles.^[21]^[22] A key distinction between pure and mixed states lies in the idempotency property: for pure states, \rho^2 = \rho, indicating that \rho acts as a projector onto a single-dimensional subspace, whereas for mixed states, \rho^2 \neq \rho and \operatorname{Tr}(\rho^2) < 1, quantifying the degree of mixture or "purity" of the state. The eigenvalues of \rho directly represent the probabilities associated with the eigenstates, providing a probabilistic interpretation akin to classical statistical mechanics but rooted in quantum superposition. In many-particle systems, \rho operates on the tensor product Hilbert space of all particles, allowing for the treatment of correlations and indistinguishability.^[20]^[21]

Expectation values and probabilities

In quantum statistical mechanics, the density operator provides a systematic way to compute the average values of physical observables for systems that may be in mixed states. For an observable represented by a Hermitian operator A, the expectation value is given by

\langle A \rangle = \mathrm{Tr}(\rho A),

where \rho is the density operator and the trace is taken over the system's Hilbert space. This expression generalizes the pure-state case \langle \psi | A | \psi \rangle to ensembles, as \rho = \sum_i p_i | \psi_i \rangle \langle \psi_i | with probabilities p_i. Equivalently, in an orthonormal basis \{ |\psi_k\rangle \}, it expands to

\langle A \rangle = \sum_k \langle \psi_k | \rho | \psi_k \rangle.

This formulation, introduced by von Neumann, ensures that expectation values are linear and satisfy the properties of quantum averages even for incomplete knowledge of the system state.^[23] The density operator also determines measurement probabilities. For an observable A with spectral decomposition A = \sum_a a P_a, where P_a is the projector onto the eigenspace of eigenvalue a, the probability of obtaining outcome a is

P(a) = \mathrm{Tr}(\rho P_a).

This follows from the Born rule extended to mixed states and yields the full probability distribution over possible measurement results. In the basis expansion, P(a) = \sum_k \langle \psi_k | \rho P_a | \psi_k \rangle, directly linking to the state's statistical description. Von Neumann derived this as part of the probabilistic interpretation of quantum ensembles.^[23] For composite systems, such as those interacting with an environment, the reduced density operator describes the relevant subsystem by tracing out the degrees of freedom of the other part. If the total system AB has density operator \rho_{AB}, the reduced density operator for subsystem A is

\rho_A = \mathrm{Tr}_B (\rho_{AB}),

where the partial trace over B sums over its basis states: \rho_A = \sum_j \langle \phi_j^B | \rho_{AB} | \phi_j^B \rangle for orthonormal basis \{ |\phi_j^B\rangle \} of B. This operation preserves the trace and positivity of \rho_A, enabling the computation of subsystem expectation values \langle A \rangle = \mathrm{Tr}_A (\rho_A A) without full knowledge of the environment, crucial for open quantum systems in statistical mechanics. The concept arises naturally in the density formalism for multipartite states.^[24] A concrete example is the position probability density in a one-dimensional system, given by the diagonal matrix element \langle x | \rho | x \rangle, which represents the probability of finding the particle at position x upon measurement. Integrating over regions yields spatial probabilities, analogous to the classical probability density but incorporating quantum correlations via \rho. This application highlights how the density operator unifies statistical predictions across position and momentum observables.^[23]

Von Neumann entropy

In quantum statistical mechanics, the von Neumann entropy serves as the quantum analog of classical thermodynamic entropy, quantifying the uncertainty or mixedness of a quantum state described by a density operator \rho. It is defined as

S(\rho) = -k \operatorname{Tr}(\rho \ln \rho),

where k is Boltzmann's constant and \operatorname{Tr} denotes the trace over the Hilbert space. Equivalently, in the eigenbasis of \rho with eigenvalues \lambda_i \geq 0 (satisfying \sum_i \lambda_i = 1), it takes the form

S(\rho) = -k \sum_i \lambda_i \ln \lambda_i.

This expression was introduced by John von Neumann as a measure of entropy for quantum ensembles, extending the concept to account for both pure and mixed states.^[25] The von Neumann entropy exhibits several key properties that mirror and generalize those of classical entropy. It is non-negative, S(\rho) \geq 0, with equality holding if and only if \rho represents a pure state (i.e., \rho = |\psi\rangle\langle\psi| for some state vector |\psi\rangle). For product states of independent systems, \rho = \rho_1 \otimes \rho_2, the entropy is additive: S(\rho) = S(\rho_1) + S(\rho_2). These properties arise from the spectral decomposition of \rho and the concavity of the function x \ln x.^[25] Von Neumann's formulation generalizes the classical Shannon entropy, which applies to probability distributions, to the quantum domain via the density operator; the Shannon entropy emerges in the diagonal limit where \rho is classical. This connection highlights the entropy's role in bridging quantum mechanics and information theory. The entropy achieves its maximum value of S(\rho) \leq k \ln \dim(\mathcal{H}) for a d-dimensional Hilbert space \mathcal{H} when \rho is the maximally mixed state \rho = I/d, representing complete ignorance of the quantum state.^[25]

Equilibrium Ensembles

Microcanonical ensemble

In quantum statistical mechanics, the microcanonical ensemble provides the foundational description for an isolated system with a fixed total energy E, volume V, and number of particles N, assuming all accessible quantum states within a narrow energy shell around E are equally probable.^[26] This ensemble corresponds to the quantum analog of the classical microcanonical ensemble, where the system's dynamics are constrained to a hypersurface of constant energy in phase space.^[27] The density operator \rho for the quantum microcanonical ensemble is defined as the uniform projector onto the subspace of energy eigenstates with energies approximately equal to E:

\rho = \frac{1}{\Omega(E)} \sum_{E_i \approx E} |E_i\rangle \langle E_i|,

where \Omega(E) denotes the degeneracy, or number of such eigenstates, ensuring \mathrm{Tr}(\rho) = 1.^[26] In a more precise formulation for a sharp energy constraint, it takes the form \rho(E) = [\delta](/page/Delta)(H - E) / \Omega(E), where H is the Hamiltonian operator.^[26] The role of the partition function is played by the density of states \Omega(E), which approximates the trace \Omega(E) \approx \mathrm{Tr}[\delta(H - E)], counting the effective number of accessible microstates at energy E.^[26] This quantity encodes the phase space volume available to the system and forms the basis for thermodynamic quantities in the ensemble.^[26] The entropy S of the microcanonical ensemble is given by S = [k](/page/K) \ln \Omega(E), where [k](/page/K) is Boltzmann's constant, which coincides with the von Neumann entropy S = -[k](/page/K) \mathrm{Tr}(\rho \ln \rho) for this uniform density operator.^[27] The temperature T emerges thermodynamically from the relation $1/T = \partial S / \partial E, reflecting how the density of states grows with energy.^[26] In quantum systems, ergodicity is established through von Neumann's quantum ergodic theorem, which asserts that for an isolated system in a pure energy eigenstate, the long-time average of an observable equals the microcanonical ensemble average, provided the system satisfies suitable mixing conditions.^[27] This theorem justifies the use of ensemble averages to compute equilibrium properties from the dynamics of individual states.^[27]

Canonical ensemble

The canonical ensemble in quantum statistical mechanics describes an isolated quantum system of fixed particle number that is in thermal contact with a large heat reservoir at temperature T, allowing energy exchange while maintaining thermal equilibrium. This setup extends the microcanonical ensemble by incorporating temperature effects through the reservoir, enabling the calculation of finite-temperature properties such as average energies and heat capacities. The formalism relies on the density operator to represent the statistical state of the system, capturing the probabilistic nature of quantum measurements in thermal environments.^[28] The density operator for the canonical ensemble is given by

\rho = \frac{e^{-\beta H}}{Z},

where H is the Hamiltonian of the system, \beta = 1/(kT) with k the Boltzmann constant, and Z is the partition function defined as the trace

Z = \mathrm{Tr}\left(e^{-\beta H}\right).

This form arises from maximizing the von Neumann entropy subject to constraints on the average energy, ensuring \mathrm{Tr}(\rho) = 1 and \mathrm{Tr}(\rho H) = \langle H \rangle. The partition function Z normalizes the operator and serves as the central quantity for thermodynamic derivations.^[29]^[30] From the partition function, key thermodynamic potentials follow directly. The Helmholtz free energy is

F = -kT \ln Z,

which connects statistical mechanics to classical thermodynamics via Legendre transforms. The average energy is obtained by differentiating the free energy, yielding

\langle H \rangle = -\frac{\partial \ln Z}{\partial \beta}.

These relations allow computation of fluctuations and response functions, such as the heat capacity C = \partial \langle H \rangle / \partial T.^[28] The canonical ensemble can be derived from the microcanonical ensemble by considering the total system as the small system of interest coupled to a large reservoir with fixed total energy E_\mathrm{tot}. The probability of the system having energy E is proportional to the microcanonical density of states of the reservoir at E_\mathrm{tot} - E, approximated for large reservoir size as P(E) \propto e^{S_\mathrm{res}(E_\mathrm{tot} - E)/k} \approx e^{-\beta E} using the reservoir's temperature \beta^{-1}. Energy fluctuations in the system scale as \sqrt{N} for N particles, negligible relative to the reservoir's extensive energy, justifying the fixed-T approximation.^[29]^[30] A representative example is the quantum harmonic oscillator with Hamiltonian H = \hbar \omega (a^\dagger a + 1/2), where the energy eigenvalues are \epsilon_n = \hbar \omega (n + 1/2) for n = 0, 1, 2, \dots. The partition function is the geometric series

Z = \sum_{n=0}^\infty e^{-\beta \hbar \omega (n + 1/2)} = \frac{e^{-\beta \hbar \omega / 2}}{1 - e^{-\beta \hbar \omega}} = \frac{1}{2 \sinh(\beta \hbar \omega / 2)},

leading to average energy \langle H \rangle = \hbar \omega (1/2 + 1/(e^{\beta \hbar \omega} - 1)), which interpolates between zero-point energy at low T and classical equipartition at high T. This illustrates quantum corrections to classical behavior in thermal systems.

Grand canonical ensemble

The grand canonical ensemble provides a statistical description of a quantum mechanical system that can exchange both energy and particles with a large reservoir, maintaining fixed temperature T, volume V, and chemical potential \mu. This setup is ideal for modeling open quantum systems, such as gases in contact with particle reservoirs, where the number of particles N fluctuates around an average value. Unlike the canonical ensemble, which fixes N, the grand canonical formulation accounts for these fluctuations naturally, making it essential for studying phenomena like quantum phase transitions in dilute gases.^[31] The equilibrium state of the system is represented by the density operator

\hat{\rho} = \frac{e^{-\beta (\hat{H} - \mu \hat{N})}}{\Xi},

where \beta = 1/(k_B T), \hat{H} is the system's Hamiltonian, \hat{N} is the particle number operator, k_B is Boltzmann's constant, and \Xi is the grand partition function given by

\Xi = \Tr \left[ e^{-\beta (\hat{H} - \mu \hat{N})} \right].

This density operator ensures that expectation values of observables are computed as \langle \hat{O} \rangle = \Tr (\hat{\rho} \hat{O}), with the trace taken over the Hilbert space of the system. The chemical potential \mu controls the average particle density, analogous to how temperature governs energy distribution.^[32] Key thermodynamic quantities derive from \Xi. The average particle number is

\langle N \rangle = \frac{1}{\beta} \left( \frac{\partial \ln \Xi}{\partial \mu} \right)_{\beta, V} = k_B T \left( \frac{\partial \ln \Xi}{\partial \mu} \right)_{T, V},

which determines the equilibrium particle density. Particle number fluctuations are quantified by the variance

\Delta N^2 = \langle N^2 \rangle - \langle N \rangle^2 = k_B T \left( \frac{\partial \langle N \rangle}{\partial \mu} \right)_{T, V},

reflecting the compressibility of the system and becoming negligible relative to \langle N \rangle^2 for large systems. These relations stem from the logarithmic derivatives of the grand potential \Omega = -k_B T \ln \Xi, which equals the negative of the pressure-volume product, \Omega = -p V.^[33]^[34] In the grand canonical ensemble, the system is treated as coupled to a much larger reservoir, justifying the approximation of fixed \mu and T; for sufficiently large reservoirs, this yields results equivalent to the canonical ensemble in the thermodynamic limit where N fluctuations are suppressed. For non-interacting particles, the grand partition function simplifies to a product over single-particle energy levels \epsilon_k: for bosons, \Xi = \prod_k \frac{1}{1 - e^{-\beta (\epsilon_k - \mu)}}; for fermions, \Xi = \prod_k \left(1 + e^{-\beta (\epsilon_k - \mu)}\right).^[35] This form highlights the role of quantum statistics in determining equilibrium properties without requiring a full many-body diagonalization.

Indistinguishable Particles and Quantum Statistics

Symmetrization principle

In quantum mechanics, identical particles are indistinguishable, meaning that the physical state of a system must remain unchanged under the exchange of any two such particles, up to a phase factor. This requirement is formalized by the exchange operator \hat{P}_{ij}, which swaps the labels of particles i and j in the many-particle wavefunction or state vector. For a valid state |\psi\rangle in the many-body Hilbert space, the action of the exchange operator yields \hat{P}_{ij} |\psi\rangle = \pm |\psi\rangle, where the + sign corresponds to bosons and the - sign to fermions.^[36] The symmetrization principle extends this to the full permutation group of N particles, requiring that the wavefunction or state transform according to irreducible representations of the symmetric group S_N: totally symmetric for bosons and totally antisymmetric for fermions. This distinction arises from the spin-statistics theorem, which connects the symmetry type to the particle's intrinsic spin—integer spin for bosons and half-integer spin for fermions—without relying on relativistic considerations here.^[36] To construct such states, antisymmetric wavefunctions for fermions are built using Slater determinants, which provide an orthonormal basis ensuring the required antisymmetry under permutations. For example, for N fermions in single-particle orbitals \phi_k(\mathbf{r}), the state is given by

\Psi(\mathbf{r}_1, \dots, \mathbf{r}_N) = \frac{1}{\sqrt{N!}} \det \begin{pmatrix} \phi_1(\mathbf{r}_1) & \cdots & \phi_1(\mathbf{r}_N) \\ \vdots & \ddots & \vdots \\ \phi_N(\mathbf{r}_1) & \cdots & \phi_N(\mathbf{r}_N) \end{pmatrix},

originally introduced to satisfy the antisymmetry condition. The bosonic analog employs permanents, replacing the determinant with a sum over permutations without sign factors, yielding totally symmetric states.^[37] These symmetry requirements imply that observables in identical particle systems are unaffected by non-symmetric operators, as only the symmetric or antisymmetric components contribute to expectation values. Arbitrary phase factors under exchange are physically irrelevant, as they do not alter measurable probabilities or matrix elements.^[36]

Bose-Einstein and Fermi-Dirac statistics

In quantum statistical mechanics, Bose-Einstein statistics governs the distribution of indistinguishable bosons, which obey Bose-Einstein commutation relations, while Fermi-Dirac statistics applies to indistinguishable fermions, which follow anticommutation relations. These statistics emerge when treating non-interacting particles in the grand canonical ensemble, where the total grand partition function factors into independent contributions from each single-particle state labeled by momentum or quantum number k with energy \varepsilon_k. The key result is the average occupation number \langle n_k \rangle for each state, which differs fundamentally from the classical Maxwell-Boltzmann limit due to quantum indistinguishability. The derivation begins with the single-particle Hamiltonian for mode k, H_k = \varepsilon_k n_k, where n_k is the number operator. For bosons, the annihilation and creation operators a_k and a_k^\dagger satisfy the commutation relation [a_k, a_k^\dagger] = 1, allowing the number states to have eigenvalues n_k = 0, 1, 2, \dots. The grand partition function for this mode is then

\Xi_k = \sum_{n_k=0}^\infty e^{-\beta n_k (\varepsilon_k - \mu)} = \frac{1}{1 - e^{-\beta (\varepsilon_k - \mu)}},

where \beta = 1/(k_B T) and \mu is the chemical potential, provided the series converges. The average occupation number is obtained as \langle n_k \rangle = \frac{1}{\beta} \frac{\partial \ln \Xi_k}{\partial \mu}, yielding

\langle n_k \rangle = \frac{1}{e^{\beta (\varepsilon_k - \mu)} - 1}.

This form was first derived by extending Bose's counting of indistinguishable quanta to an ideal gas.^[38] For fermions, the operators c_k and c_k^\dagger satisfy the anticommutation relation \{c_k, c_k^\dagger\} = 1, restricting the eigenvalues to n_k = 0 or $1 due to the Pauli exclusion principle. The grand partition function simplifies to

\Xi_k = \sum_{n_k=0}^1 e^{-\beta n_k (\varepsilon_k - \mu)} = 1 + e^{-\beta (\varepsilon_k - \mu)},

and the average occupation number is

\langle n_k \rangle = \frac{1}{e^{\beta (\varepsilon_k - \mu)} + 1}.

This ensures \langle n_k \rangle \leq 1 for all states, directly enforcing the exclusion principle. The statistics were introduced independently by Fermi and Dirac in their quantization of the ideal gas, emphasizing antisymmetric wavefunctions for identical particles.^[39] The chemical potential \mu is bounded by \mu < \min_k \varepsilon_k for bosons to prevent divergence in the occupation number (e.g., \mu < [0](/page/0) if energies are measured from zero). For fermions at T=0, the distribution reduces to a step function \langle n_k \rangle = \theta(\mu - \varepsilon_k), with \mu \leq [0](/page/0) ensuring no occupation above the vacuum level in the non-degenerate limit, though degeneracy allows positive \mu up to the Fermi energy. These distributions highlight the quantum enhancement of low-energy occupations for bosons and suppression for fermions compared to classical statistics.

Ideal quantum gases

Ideal quantum gases consist of non-interacting particles obeying either Bose-Einstein or Fermi-Dirac statistics, leading to distinct thermodynamic behaviors compared to classical gases, particularly at low temperatures or high densities where quantum effects dominate. The occupation numbers \langle n_k \rangle follow the respective distributions derived from symmetrization principles, enabling the computation of macroscopic properties like pressure and energy from sums over single-particle states. These systems serve as foundational models for understanding quantum degeneracy in dilute atomic gases and electron systems in solids. For the ideal Bose gas of non-interacting bosons, the equation of state is expressed through the pressure

P = \frac{k_B T}{\lambda^3} g_{5/2}(z),

where \lambda = h / \sqrt{2 \pi m k_B T} is the thermal de Broglie wavelength, z = e^{\beta \mu} is the fugacity (\beta = 1 / k_B T, \mu the chemical potential), and the Bose-Dirac integral is g_\nu(z) = \sum_{l=1}^\infty z^l / l^\nu. This relation, independent of volume above the condensation point, reflects the saturability of excited states and was originally derived by Einstein in his analysis of quantum statistics for monatomic gases.^[40] In contrast, the ideal Fermi gas of non-interacting fermions at zero temperature fills all states up to the Fermi energy \epsilon_F = \frac{\hbar^2}{2m} (3 \pi^2 n)^{2/3}, where n is the particle number density. The resulting zero-temperature pressure, known as degeneracy pressure, is P = \frac{2}{5} n \epsilon_F, arising from the Pauli exclusion principle's constraint on state occupancy and providing kinetic support against collapse even without interactions. This ground-state configuration was introduced by Fermi in his quantization of identical particles and refined by Sommerfeld for conduction electrons.^[39]^[41] At low temperatures, the specific heat C_V of the ideal Fermi gas is linear in T, C_V = \gamma T with \gamma = \frac{\pi^2 k_B^2}{3} g(\epsilon_F) and g(\epsilon) the density of states at the Fermi level, due to thermal smearing of the Fermi surface over an energy window \sim k_B T. For the ideal Bose gas above the condensation temperature, C_V approaches the classical value \frac{3}{2} N k_B but incorporates an exponential gap in excitation probabilities from the fugacity z < 1, leading to suppressed contributions from higher states.^[42]^[43] Quantum corrections to the classical ideal gas law emerge when the de Broglie wavelength \lambda becomes comparable to the interparticle spacing n^{-1/3}, quantified by the small parameter n \lambda^3. The first-order virial expansion for the pressure is P = n k_B T \left[ 1 \pm \frac{1}{2^{5/2}} n \lambda^3 + \mathcal{O}((n \lambda^3)^2) \right], where the + sign applies to bosons and - to fermions, capturing the initial enhancement (bosons) or suppression (fermions) of pressure due to statistical correlations.^[44]

References

[1]
Quantum Statistical Mechanics
The philosophy of statistical mechanics is that real systems are complicated. At best you have an approximate description of them. For example, you model the ...
[2]
Review of quantum statistical mechanics (Chapter 1)
The purpose of this chapter is to give a mini-review of the basic concepts of quantum statistical mechanics as applied to noninteracting systems of particles.
[3]
As the world looks for quantum solutions, Bose statistics turns 100
Feb 9, 2024 · Interactions between Bose and Einstein resulted in a new statistical tool, the Bose-Einstein statistics, introduced for photons in 1924 by Bose ...Missing: original | Show results with:original
[4]
Quantum Statistical Mechanics - Book - IOPscience
This book establishes the foundations of non-equilibrium quantum statistical mechanics in order to support students and academics in developing and building ...Examples And Applications... · Probability In Quantum... · Evolution Of The Canonical...
[5]
On Quantum Statistical Mechanics: A Study Guide - Majewski - 2017
Nov 23, 2017 · In this paper, we will try to give an overview and road map to the area of quantum statistical mechanics without becoming too diverted by ...
[6]
Quantum Statistical Mechanics - an overview | ScienceDirect Topics
Quantum statistical mechanics is the study of probability potentials of individual particles as well as bulk properties.
[7]
Quantum information and statistical mechanics: an introduction to ...
Jun 28, 2013 · This is a short review on an interdisciplinary field of quantum information science and statistical mechanics.
[8]
On the foundations of statistical mechanics - ScienceDirect.com
Jul 6, 2025 · Statistical mechanics (SM) finds its motivation in the study of the behavior and properties of macroscopic bodies.
[9]
[PDF] Boltzmann's Ergodic Hypothesis, a Conjecture for Centuries?
May 16, 1994 · Boltzmann's more than one hundred year old ergodic hypothesis is given. The existing main results, the majority of which is connected with the ...
[10]
Boltzmann's Ergodic Hypothesis - jstor
28). In a short paper (1868 a) Boltzmann discusses the relation between the atom- istic structure of matter and the description of physical processes ...
[11]
[PDF] Elementray Principles in Statistical Mechanics. - Project Gutenberg
Jan 22, 2016 · This PDF file ... Mechanics, by Josiah Willard Gibbs. *** END OF THIS PROJECT GUTENBERG EBOOK ELEMENTARY PRINCIPLES STATISTICAL MECHANICS ***.
[12]
[PDF] General linear dynamics – quantum, classical or hybrid - arXiv
Mar 18, 2011 · Schrödinger equation can be transferred to the cases of Liouville or von Neumann equations with the help of Liouville space notions. This ...
[13]
On the Relation Between Classical and Quantum Statistical Mechanics
Abstract. Classical and quantum statistical mechanics are compared in the high temperature limit β = ί/kT^O. While this limit is rather trivial for spin.
[14]
[PDF] A New Look at the Quantum Liouville Theorem - PDXScholar
Aug 9, 2020 · It is further shown that such a formulation will provide the exact correspondence between classical and quantum statistical mechanics with the ...
[15]
[PDF] Chapter 4 Introduction to many-body quantum mechanics
The many-body Hilbert space can thus be split into orthogonal subspaces, one in which particles pick up a − sign and are called fermions, and the other where ...
[16]
[PDF] Many–Particle Systems
An example for a single–particle operator is the kinetic energy operator ... In all these systems the Hamiltonian is a sum of one– and two–particle operators.
[17]
[PDF] Notes on Quantum Mechanics - Physics
Evaluate the commutator [x, p] of the position operator, x, and the momentum operator, p = −ih ∂/∂x. Problem 5.4. Consider all the possible commutators ...
[18]
[PDF] LECTURE NOTES ON QUANTUM MECHANICS
These notes are for the two-semester graduate level quantum mechanics class taught at Michi- gan State University. Although they are more terse than a typical ...
[19]
Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik
Neumann, J. von. "Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik." Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen ...Missing: John pdf
[20]
[PDF] The Quantum Density Matrix and Its Many Uses - arXiv
Aug 16, 2023 · The axiomatic formulation of quantum mechanics assumes a description of ... description in terms of a density matrix with specific properties.
[21]
[PDF] 5.74 Introductory Quantum Mechanics II - MIT OpenCourseWare
Properties of the density matrix. 1) ρ is Hermetian: ρ ... The equation of motion for the density matrix follows naturally from the definition of ρ and the.
[22]
[PDF] Chapter 9 Density Matrices
Definition 9.1. The density matrix ρ for the pure state |ψ > is given by ρ := |ψ ><ψ |. This density matrix has the following properties: I) ρ2 = ρ projector.
[23]
[PDF] Von Neumann's 1927 Trilogy on the Foundations of Quantum ... - arXiv
May 20, 2025 · In his first solo paper on quantum mechanics, (von Neumann. 1927a), von Neumann would essentially reject the methodology of the Jordan ...
[24]
[PDF] Chapter 7. Open Quantum Systems
This chapter discusses the effects of weak interaction of quantum systems with their environment. Some parts of this material are on the fine line between ...
[25]
Statistical entropy of quantum systems
### Summary of von Neumann Entropy from arXiv:2412.15316
[26]
[PDF] Statistical Mechanics I - MIT OpenCourseWare
• Microcanonical ensemble: As the internal energy has a fixed value E, a density matrix that includes this constraint is ρ(E) = δ(H −. E) . (VI.84). Ω(E). In ...
[27]
[PDF] Proof of the Ergodic Theorem and the H-Theorem in Quantum ...
The main notions of statistical mechanics are re-interpreted in a quantum-mechanical way, the ergodic theorem and the H-theorem are formulated and proven ( ...
[28]
[PDF] Physics 127a: Class Notes - Lecture 12: Quantum Statistical ...
density matrix in the canonical ensemble. ¯ρ = e−βH . (36). Let us evaluate ¯ρ in the coordinate representation ¯ρ(x,x0;β). A convenient way to do this is to ...
[29]
[PDF] Lecture 7: Ensembles
An important point about the canonical ensemble is that we derived a result about the system only. The partition function is a sum over microstates of the ...<|control11|><|separator|>
[30]
[PDF] Statistical Mechanics at Fixed Temperature (Canonical Ensemble)
Everything we do will be correct. 6-23. Page 24. Partition function of quantum harmonic oscillator. Z1. ≡ Z(one oscillator) = X states,n=0,1,2... e−β n. = e.
[31]
[PDF] Chapter 10 Grand canonical ensemble
The grand canonical ensemble is a generalization of the canonical ensemble where the restriction to a definite number of particles is removed.
[32]
Grand Canonical Ensemble - an overview | ScienceDirect Topics
... of an open, isothermal quantum system. The density operator ρ ˆ for a quantum grand canonical ensemble is given by. (120) ρ ˆ ( N ˆ ) = Ξ − 1 ( V , β , μ ) ...
[33]
[PDF] The Grand Canonical Ensemble
Apr 10, 2001 · T,V,µ fixed, E and N fluctuate. Ω = −kT lnΞ = −kT ln Z, where Ω is the grand potential and. Ξ(= Z), is the grand partition function. Toc.
[34]
[PDF] 12. Grandcanonical Ensemble - DigitalCommons@URI
Dec 16, 2015 · . Fluctuations in particle number: hN2i−hNi2 = N2 V kBTκT . An alternative expression for hN2i−hNi2 is calculated in exercise [tex95]. The ...
[35]
[PDF] Chapter 8 - Identical Particles - MIT OpenCourseWare
May 8, 2022 · IDENTICAL PARTICLES IN QUANTUM MECHANICS. 153. Figure 8.2: Example ... The symmetrization postulate solves the exchange degeneracy problem.
[36]
[2106.14679] Permanent variational wave functions for bosons - arXiv
Jun 25, 2021 · We study the performance of permanent states (the bosonic counterpart of the Slater determinant state) as approximating functions for bosons.Missing: mechanics | Show results with:mechanics
[37]
[PDF] Quantum Theory of a Monoatomic Ideal Gas A translation of ...
Jan 9, 2015 · Page 1. Quantum Theory of a Monoatomic Ideal Gas. A translation of Quantentheorie des einatomigen idealen Gases (Einstein, 1924). A. Einstein.
[38]
[PDF] On Quantizing an Ideal Monatomic Gas - Gilles Montambaux
The aim of the present paper is to present a method of quantization of an ideal gas which, according to our opinion, is as independent of arbitrary.Missing: Sul perfetto monoatomico
[39]
[PDF] Quantentheorie des einatomigen idealen Gases.
Quantentheorie des einatomigen idealen Gases. Zweite Abhandlung. Von A. EINSTEIN. paper_1925_00. Instituut-Lorentz, Leiden University.
[40]
[PDF] Free electron Fermi gas (Sommerfeld, 1928)
Each point can have 2 electrons (because of spin). After filling in. N electrons, the result is a spherical sea of electrons called the. Fermi sphere.
[41]
[PDF] PHYS 445 Lecture 30 - Low temperature Fermi gas 30
Heat capacity for ideal Fermi gas. This lecture contains a somewhat lengthy mathematical examination of non-interacting. Fermi gases at low temperatures.
[42]
Specific heat of an ideal Bose gas above the Bose condensation ...
Sep 1, 2004 · The formula yields the exact values of the specific heat and its first derivative at and has an accuracy of at least 0.3% over the full range ...
[43]
[PDF] 3. Quantum Gases - DAMTP
Recall that in quantum mechanics, particles come in two classes: bosons and fermions. ... in terms of the chemical potential and temperature: N = N(µ, T).