Thermal Density Operators

For a quantum system with Hamiltonian $H$ in equilibrium at absolute temperature $T$, the canonical thermal state is the Gibbs density operator

$$\rho_\beta = \frac{e^{-\beta H}}{Z(\beta)}, \qquad Z(\beta) = \operatorname{Tr} e^{-\beta H},$$

where

$$\beta = \frac{1}{k_{\mathrm B}T}.$$

The partition function $Z$ normalizes the operator:

$$\operatorname{Tr}\rho_\beta=1.$$

This page is the canonical home for the general operator structure of Gibbs states. Density operators as quantum states belong to Core Formalism. The later canonical-ensemble page develops fixed-temperature thermodynamics, and the grand-canonical page replaces $H$ by $H-\mu N$ when particle number can fluctuate.

Assumptions Behind the Formula

In a finite-dimensional Hilbert space, $e^{-\beta H}$ is well defined for every self-adjoint $H$ and every finite real $\beta$. For the usual positive-temperature state, $\beta>0$.

In an infinite-dimensional Hilbert space, more is required:

1. $H$ must be self-adjoint.
2. The exponential $e^{-\beta H}$ must be trace class.
3. The partition function must satisfy

$$0<Z(\beta)<\infty.$$

A Hamiltonian bounded below with a sufficiently fast-growing discrete spectrum often satisfies these conditions for $\beta>0$. They can fail for an unconfined continuum system in infinite volume or for spectra with excessive degeneracy.

When $Z$ diverges, the formal ratio is not a density operator. One must regulate the system, work in finite volume, use densities per volume, or adopt an infinite-system equilibrium formalism.

Operator Exponential

The notation $e^{-\beta H}$ is defined by functional calculus. In finite dimension it can be represented by the convergent series

$$e^{-\beta H} = \sum_{m=0}^{\infty} \frac{(-\beta H)^m}{m!}.$$

The spectral theorem gives the more useful form. If

$$H = \sum_a E_a P_a,$$

where $P_a$ projects onto the eigenspace with energy $E_a$, then

$$e^{-\beta H} = \sum_a e^{-\beta E_a}P_a.$$

Therefore

$$\rho_\beta = \frac{1}{Z} \sum_a e^{-\beta E_a}P_a.$$

The formula is basis independent. The energy eigenbasis merely makes its probabilities visible.

Energy-Basis Representation

For a nondegenerate discrete spectrum,

$$H|E_n\rangle = E_n|E_n\rangle,$$

the thermal state is

$$\rho_\beta = \sum_n p_n |E_n\rangle\langle E_n|,$$

with

$$p_n = \frac{e^{-\beta E_n}}{Z}, \qquad Z = \sum_n e^{-\beta E_n}.$$

Thus an energy measurement returns $E_n$ with Boltzmann probability $p_n$.

The state is generally mixed:

$$\operatorname{Tr}\rho_\beta^2 < 1$$

whenever more than one energy eigenstate has nonzero weight. At finite positive temperature in finite dimension, every energy eigenstate has nonzero weight and $\rho_\beta$ has full rank.

Degeneracy

Let the energy $E_a$ have degeneracy

$$g_a = \operatorname{Tr}P_a.$$

Every normalized vector in that eigenspace receives the same weight per basis state:

$$\rho_\beta P_a = \frac{e^{-\beta E_a}}{Z} P_a.$$

The probability of measuring the energy $E_a$ is

$$\Pr(E_a) = \frac{ g_a e^{-\beta E_a} }{ Z }.$$

The factor $g_a$ belongs to the probability of the energy value, not to the weight of each orthonormal microstate.

Because $\rho_\beta$ is proportional to the identity within a degenerate eigenspace, it does not choose a preferred basis there. A basis rotation within that eigenspace cannot create physical thermal coherence.

Basic Operator Properties

Positivity

For every $|\psi\rangle$,

$$\langle\psi| e^{-\beta H} |\psi\rangle \geq 0.$$

After normalization,

$$\rho_\beta\geq0.$$

Hermiticity

Since $H=H^\dagger$ and $\beta$ is real,

$$\rho_\beta^\dagger = \rho_\beta.$$

Commutation with the Hamiltonian

The Gibbs state is a function of $H$, so

$$[H,\rho_\beta]=0.$$

Stationarity

Under the same time-independent Hamiltonian,

$$\rho_\beta(t) = e^{-iHt/\hbar} \rho_\beta e^{iHt/\hbar} = \rho_\beta.$$

Stationarity is necessary for equilibrium but not sufficient. Any density operator diagonal in the energy basis is stationary, and many stationary states are not Gibbs states.

Invariance Under an Energy-Zero Shift

Replace the Hamiltonian by

$$H' = H+C\mathbb I.$$

Then

$$\begin{aligned} e^{-\beta H'} &= e^{-\beta C} e^{-\beta H}, \\ Z' &= e^{-\beta C}Z. \end{aligned}$$

The scalar factor cancels:

$$\rho_\beta' = \rho_\beta.$$

Thermal probabilities do not depend on the arbitrary zero of energy. Thermodynamic potentials such as $F=-k_{\mathrm B}T\ln Z$ do shift by the corresponding extensive energy constant, so one must distinguish state invariance from potential conventions.

Thermal Expectation Values

For any observable $A$,

$$\langle A\rangle_\beta = \operatorname{Tr} \left( \rho_\beta A \right).$$

In a discrete energy basis,

$$\langle A\rangle_\beta = \frac{1}{Z} \sum_n e^{-\beta E_n} \langle E_n|A|E_n\rangle.$$

Only the energy-basis diagonal matrix elements of $A$ enter this one-time expectation, even when $[A,H]\neq0$. Off-diagonal matrix elements remain essential for time-dependent correlations and response functions.

For a function $f(H)$,

$$\langle f(H)\rangle_\beta = \frac{1}{Z} \sum_n f(E_n)e^{-\beta E_n}.$$

The quantum trace has reduced the calculation to a weighted spectral sum.

Energy Moments from the Partition Function

Differentiating $Z$ gives

$$\frac{\partial Z}{\partial\beta} = -\operatorname{Tr} \left( He^{-\beta H} \right),$$

so the mean energy is

$$U \equiv \langle H\rangle_\beta = -\frac{\partial}{\partial\beta} \ln Z.$$

A second derivative gives the energy variance:

$$\frac{\partial^2}{\partial\beta^2} \ln Z = \langle H^2\rangle_\beta - \langle H\rangle_\beta^2.$$

Thus

$$\frac{\partial U}{\partial\beta} = -\operatorname{Var}_\beta(H) \leq 0.$$

For a temperature-independent Hamiltonian,

$$C = \frac{\partial U}{\partial T} = k_{\mathrm B}\beta^2 \operatorname{Var}_\beta(H) \geq 0.$$

These identities are canonical-ensemble statements. They appear here because they reveal what information is encoded by the operator normalization; the canonical-ensemble page develops their thermodynamic interpretation and qualifications.

Entropy of a Gibbs State

The von Neumann entropy with thermodynamic units is

$$S(\rho) = -k_{\mathrm B} \operatorname{Tr} \left( \rho\ln\rho \right).$$

For the Gibbs state,

$$\ln\rho_\beta = -\beta H-\ln Z,$$

so

$$S(\rho_\beta) = k_{\mathrm B} \left( \beta U+\ln Z \right).$$

Using $T=(k_{\mathrm B}\beta)^{-1}$ gives

$$F = U-TS = -k_{\mathrm B}T\ln Z.$$

This entropy is the entropy of the full thermal mixed state. It should not be confused with the entanglement entropy of a subsystem in a pure state, though reduced thermal states can also carry entanglement and classical correlations.

Relation to Maximum Entropy

Among density operators with fixed normalization and fixed mean energy,

$$\operatorname{Tr}\rho=1, \qquad \operatorname{Tr}(\rho H)=U,$$

the Gibbs state uniquely maximizes the von Neumann entropy under standard finite-dimensional full-rank assumptions.

Introduce multipliers $\alpha$ and $\beta$ and vary

$$\begin{aligned} \mathcal L[\rho] ={}& -\operatorname{Tr}(\rho\ln\rho) \\ &- \alpha \left( \operatorname{Tr}\rho-1 \right) \\ &- \beta \left[ \operatorname{Tr}(\rho H)-U \right]. \end{aligned}$$

Stationarity gives

$$\ln\rho = -(1+\alpha)\mathbb I - \beta H.$$

Exponentiating and imposing unit trace yields

$$\rho = \frac{e^{-\beta H}}{ \operatorname{Tr}e^{-\beta H} }.$$

The multiplier $\beta$ is identified thermodynamically as $1/(k_{\mathrm B}T)$. The maximum-entropy statement does not claim that the system dynamically reaches the Gibbs state. It identifies the least-biased state consistent with the stated constraints. Maximum Entropy Principle gives the general relative-entropy proof, multiplier duality, and boundary cases.

Relative-Entropy Characterization

The Gibbs state also minimizes the nonequilibrium free-energy functional

$$\mathcal F_T(\rho) = \operatorname{Tr}(\rho H) - T S(\rho)$$

at fixed $H$ and $T$. Using quantum relative entropy,

$$D(\rho\Vert\rho_\beta) = \operatorname{Tr} \left[ \rho \left( \ln\rho-\ln\rho_\beta \right) \right],$$

one finds

$$\mathcal F_T(\rho) - \mathcal F_T(\rho_\beta) = k_{\mathrm B}T D(\rho\Vert\rho_\beta) \geq 0.$$

This identity concerns an equilibrium variational principle. Claims about extractable work, thermal operations, or irreversible entropy production require an operational or dynamical framework and belong in quantum thermodynamics.

High-Temperature Limit

For a finite $D$-dimensional Hilbert space,

$$\lim_{\beta\to0} \rho_\beta = \frac{\mathbb I_D}{D}.$$

The state becomes maximally mixed. To first order,

$$\begin{aligned} \rho_\beta ={}& \frac{1}{D} \left[ \mathbb I_D - \beta \left( H-\overline E\,\mathbb I_D \right) \right] \\ &+ O(\beta^2), \end{aligned}$$

where

$$\overline E = \frac{\operatorname{Tr}H}{D}.$$

The subtraction of $\overline E$ preserves normalization at first order.

In an infinite-dimensional space, the formal maximally mixed state $\mathbb I/D$ does not exist. The limit $\beta\to0^+$ can make $Z$ diverge, so high-temperature statements must be made for regulated systems or appropriate observables.

Zero-Temperature Limit

Let the ground energy be $E_0$ and the ground-space projector be $P_0$ with degeneracy $g_0=\operatorname{Tr}P_0$.

As $\beta\to\infty$,

$$\rho_\beta \longrightarrow \frac{P_0}{g_0}.$$

For a nondegenerate ground state,

$$\rho_\beta \longrightarrow |E_0\rangle\langle E_0|.$$

For a degenerate ground space, the Gibbs limit without additional fields or boundary conditions is the equal mixture on that space, not a selected pure ground state.

If the first excited energy is separated by a gap

$$\Delta = E_1-E_0>0,$$

low-temperature corrections are suppressed by factors of order

$$e^{-\beta\Delta}.$$

In a thermodynamic limit the gap can close, degeneracies can grow, and the order of $\beta\to\infty$, volume, and source limits can matter.

Negative Temperatures

The formula permits $\beta<0$ only when

$$Z(\beta) = \operatorname{Tr} e^{-\beta H}$$

remains finite. This typically requires a spectrum bounded above as well as below.

For $\beta<0$, higher-energy states receive larger Boltzmann weights. Such a state is hotter than every positive-temperature Gibbs state in the thermodynamic sense associated with population inversion; it is not colder than zero.

An unbounded-above Hamiltonian such as the harmonic oscillator does not admit a normalizable negative-temperature Gibbs state. Population inversion in a subsystem is not by itself enough; the equilibrium constraints and bounded spectrum must be specified.

Imaginary-Time Interpretation

Real-time evolution is generated by

$$U(t) = e^{-iHt/\hbar}.$$

After the formal substitution

$$t=-i\tau,$$

the propagator becomes

$$e^{-H\tau/\hbar}.$$

Setting

$$\tau = \beta\hbar$$

produces the unnormalized Gibbs operator $e^{-\beta H}$. The partition function is its trace:

$$Z(\beta) = \operatorname{Tr} e^{-H\tau/\hbar}.$$

In path-integral language, the trace identifies the endpoints and produces an imaginary-time circle of circumference $\beta\hbar$. Bosonic and fermionic fields obey different imaginary-time boundary conditions.

This relation is structural, not a claim that a physical system literally evolves for imaginary time. Imaginary Time owns the operator semigroup, trace-class conditions, ensemble generator, and open-versus-closed boundary dictionary. The full path-integral construction lives in Euclidean and Imaginary-Time Path Integrals.

Thermal State Versus Thermalization

The Gibbs operator specifies an equilibrium state. It does not specify how that state is prepared or reached.

A closed isolated system evolves unitarily:

$$\rho(t) = e^{-iHt/\hbar} \rho(0) e^{iHt/\hbar}.$$

Its full spectrum as a density operator is conserved, so generic unitary evolution cannot turn an arbitrary pure global state into the mixed state $\rho_\beta$. Thermalization of an isolated many-body system concerns local observables, reduced states, dephasing, and ensembles, not literal convergence of the full pure state in trace norm.

An open system can relax toward a Gibbs state under appropriate bath, weak-coupling, Markovian, and detailed-balance assumptions. A thermal fixed point is a dynamical conclusion, not part of the definition of $\rho_\beta$.

Quantum Annealing may use a Gibbs distribution only as a declared equilibrium comparator or as a consequence of a licensed dynamical model; it owns the schedule-dependent endpoint and freeze-out record. This page retains the definition and thermodynamic interpretation of equilibrium density operators.

Reduced Equilibrium at Strong Coupling

Let a system $S$ interact with a bath $B$:

$$H_{\mathrm{tot}} = H_S+H_B+H_I.$$

The joint equilibrium state is

$$\rho_{SB}^{\mathrm{eq}} = \frac{ e^{-\beta H_{\mathrm{tot}}} }{ Z_{\mathrm{tot}} }.$$

The reduced system state is

$$\rho_S^{\mathrm{eq}} = \operatorname{Tr}_B \rho_{SB}^{\mathrm{eq}}.$$

When $H_I$ is not negligible, this need not equal

$$\frac{e^{-\beta H_S}}{ \operatorname{Tr}_S e^{-\beta H_S} }.$$

It can be written using a Hamiltonian of mean force, which depends on temperature and coupling. The bare Gibbs state of $H_S$ is an approximation whose regime must be stated.

The canonical open-system caveats live in Strong-Coupling Open-System Effects.

Generalized Exponential States

If several commuting conserved quantities $Q_j$ are constrained, maximum entropy gives

$$\rho \propto \exp \left( -\sum_j\lambda_j Q_j \right).$$

The grand-canonical state is the important example

$$\rho_{\beta,\mu} \propto e^{-\beta(H-\mu N)}.$$

Integrable systems may require a generalized Gibbs ensemble involving many conserved charges. These states are not ordinary canonical Gibbs states unless the additional multipliers vanish or the charges are redundant.

For noncommuting constraints, the variational problem is still meaningful but its operational interpretation and parameterization require care. One should not import classical joint-distribution intuition without checking the operator structure. The general construction and Kubo–Mori response matrix are developed in Maximum Entropy Principle.

Example: Two-Level System

Choose the ground energy to be zero and let

$$H = \Delta |1\rangle\langle1|, \qquad \Delta>0.$$

Then

$$Z = 1+e^{-\beta\Delta}$$

and

$$\rho_\beta = \frac{ |0\rangle\langle0| + e^{-\beta\Delta} |1\rangle\langle1| }{ 1+e^{-\beta\Delta} }.$$

The excited-state population is

$$p_1 = \frac{1}{ e^{\beta\Delta}+1 }.$$

The mean energy is

$$U = \frac{\Delta}{ e^{\beta\Delta}+1 }.$$

As $T\to0^+$, the state approaches $|0\rangle\langle0|$. As $T\to\infty$, both levels have probability $1/2$.

The plus sign in this two-level expression does not make it a Fermi–Dirac gas. It follows from a partition function with exactly two available levels for one finite system.

Example: Harmonic Oscillator

For

$$H = \hbar\omega \left( a^\dagger a+\frac{1}{2} \right),$$

the partition function is

$$\begin{aligned} Z &= \sum_{n=0}^{\infty} e^{-\beta\hbar\omega(n+1/2)} \\ &= \frac{ e^{-\beta\hbar\omega/2} }{ 1-e^{-\beta\hbar\omega} }. \end{aligned}$$

The thermal state is diagonal in number states:

$$\rho_\beta = \left( 1-e^{-\beta\hbar\omega} \right) \sum_{n=0}^{\infty} e^{-n\beta\hbar\omega} |n\rangle\langle n|.$$

The zero-point factor cancels from the normalized state. The mean occupation is

$$\langle a^\dagger a\rangle_\beta = \frac{1}{ e^{\beta\hbar\omega}-1 }.$$

For $\beta>0$, the geometric sum converges. For $\beta<0$, it diverges, illustrating why an unbounded-above spectrum cannot support a negative-temperature Gibbs state.

Example: Free Particle and Volume Regularization

For a free particle on all of $\mathbb R^d$,

$$H = \frac{\mathbf p^2}{2m},$$

the position volume is infinite and the canonical trace diverges. Place the particle in a finite box of volume $V$ first. In the large-box continuum approximation,

$$Z_1 \sim \frac{V}{ \lambda_T^d },$$

where

$$\lambda_T = \sqrt{ \frac{ 2\pi\hbar^2 }{ mk_{\mathrm B}T } }$$

is the thermal de Broglie wavelength.

The divergence proportional to $V$ is physical state counting. It is handled by finite-volume normalization and the thermodynamic limit, not by pretending that $e^{-\beta H}$ is trace class on infinite-volume one-particle space.

Canonical Boundaries

This page owns:

• the definition and spectral form of $\rho_\beta$;
• normalization, positivity, stationarity, and energy-shift invariance;
• general thermal expectation values;
• the fixed-mean-energy Gibbs variation and relative-entropy free-energy characterization;
• high-, low-, and negative-temperature state limits;
• the imaginary-time interpretation;
• trace-class and strong-coupling caveats.

Other pages own:

• density operators as general quantum states: Core Formalism;
• trace, spectral-sum, factorization, and cumulant properties of $Z$: Partition Functions;
• canonical free energies, heat capacity, and fixed-temperature thermodynamics: Canonical Ensemble;
• variable particle number and chemical potential: Grand-Canonical Ensemble;
• suppression of Bose and Fermi exchange corrections in the dilute regime: Classical Limit of Quantum Statistics;
• general constrained-entropy inference and exponential families: Maximum Entropy Principle;
• the thermodynamic-potential network and Legendre transforms: Thermodynamic Potentials;
• equilibrium entropy identification and coarse-graining caveats: Entropy in Quantum Statistical Mechanics;
• closed-system dynamical approach to equilibrium: Relaxation and Thermalization;
• work, heat, resource theories, and fluctuation relations: Quantum Thermodynamics;
• the shared imaginary-time, KMS, spectral, and path-integral roadmap: Finite-Temperature QM Overview;
• the compact-time transform and discrete frequency-sum workflow: Matsubara Formalism Preview;
• full imaginary-time path-integral construction: dynamics and finite-temperature field-theory pages.

Common Mistakes

• Forgetting the normalization $Z$.
• Treating $e^{-\beta H}$ as trace class without checking the spectrum or volume.
• Calling every energy-diagonal stationary state thermal.
• Assuming a Gibbs state is a pure energy eigenstate at finite temperature.
• Multiplying each degenerate microstate by the degeneracy factor $g_a$.
• Interpreting basis rotations within a degenerate eigenspace as thermal coherence.
• Thinking that shifting the energy zero changes $\rho_\beta$.
• Replacing thermal expectation values by unweighted averages over eigenstates.
• Confusing the maximally mixed high-temperature limit with an infinite-dimensional identity state.
• Claiming a selected pure ground state when the zero-temperature Gibbs limit has a degenerate ground space.
• Using negative temperature for a spectrum unbounded above.
• Treating imaginary time as literal laboratory time.
• Assuming the reduced state of a strongly coupled subsystem is the bare Gibbs state of $H_S$.
• Confusing the definition of a thermal state with a proof of thermalization.
• Confusing thermal entropy with entanglement entropy.

Exercises

Energy-shift invariance

Show directly that replacing $H$ by $H+C\mathbb I$ leaves every thermal expectation value unchanged.

Solution

The shifted partition function is

$$Z' = \operatorname{Tr} e^{-\beta(H+C\mathbb I)} = e^{-\beta C}Z.$$

Therefore

$$\rho_\beta' = \frac{ e^{-\beta C}e^{-\beta H} }{ e^{-\beta C}Z } = \rho_\beta.$$

For every observable $A$,

$$\langle A\rangle_\beta' = \operatorname{Tr}(\rho_\beta' A) = \operatorname{Tr}(\rho_\beta A).$$

Degenerate ground space

A finite Hamiltonian has a threefold-degenerate ground energy $E_0$ and a gap to all excited states. Find the positive-temperature Gibbs limit as $T\to0^+$.

Solution

Let $P_0$ project onto the three-dimensional ground space. Factoring $e^{-\beta E_0}$ from numerator and denominator shows that excited-state terms vanish exponentially relative to the ground terms. Hence

$$\lim_{\beta\to\infty} \rho_\beta = \frac{P_0}{3}.$$

The limit is an equal mixture on the ground space. A pure symmetry-broken or boundary-selected ground state requires an additional selection procedure and a specified order of limits.

Two-level purity

For the two-level Hamiltonian

$$H = \Delta|1\rangle\langle1|,$$

compute $\operatorname{Tr}\rho_\beta^2$ and its limits as $T\to0^+$ and $T\to\infty$.

Solution

The probabilities are

$$p_0 = \frac{1}{1+e^{-\beta\Delta}}, \qquad p_1 = \frac{e^{-\beta\Delta}}{ 1+e^{-\beta\Delta} }.$$

Therefore

$$\operatorname{Tr}\rho_\beta^2 = \frac{ 1+e^{-2\beta\Delta} }{ \left( 1+e^{-\beta\Delta} \right)^2 }.$$

As $\beta\to\infty$, the purity tends to $1$. As $\beta\to0$, it tends to $1/2$, the purity of the maximally mixed qubit state.

Entropy identity

Starting from $\ln\rho_\beta=-\beta H-\ln Z$, derive the Gibbs-state entropy and Helmholtz free energy.

Solution

Substitute into the von Neumann entropy:

$$\begin{aligned} S(\rho_\beta) &= -k_{\mathrm B} \operatorname{Tr} \left( \rho_\beta\ln\rho_\beta \right) \\ &= k_{\mathrm B}\beta \operatorname{Tr}(\rho_\beta H) + k_{\mathrm B}\ln Z \operatorname{Tr}\rho_\beta \\ &= k_{\mathrm B} \left( \beta U+\ln Z \right). \end{aligned}$$

Since $k_{\mathrm B}T\beta=1$,

$$\begin{aligned} F &= U-TS \\ &= -k_{\mathrm B}T\ln Z. \end{aligned}$$

Bare Gibbs state at strong coupling

Explain why

$$\rho_S^{\mathrm{eq}} = \operatorname{Tr}_B \frac{ e^{-\beta(H_S+H_B+H_I)} }{ Z_{\mathrm{tot}} }$$

need not equal $e^{-\beta H_S}/Z_S$.

Solution

When $H_I$ is non-negligible, the exponential generally does not factor:

$$e^{-\beta(H_S+H_B+H_I)} \neq e^{-\beta H_S} e^{-\beta H_B}.$$

The joint equilibrium contains system–bath correlations and interaction-energy contributions. Tracing over the bath therefore produces a reduced state governed by a Hamiltonian of mean force rather than necessarily by the bare $H_S$.

The bare Gibbs approximation can become accurate in a suitable weak-coupling regime, but it is an approximation with assumptions, not an operator identity.

Cross-Links

• Volume Laws

• Statistical Ensembles Overview

• Microcanonical Ensemble

• Canonical Ensemble

• Grand-Canonical Ensemble

• Partition Functions

• Thermodynamic Potentials

• Entropy in Quantum Statistical Mechanics

• Thermal Entropy vs Entanglement Entropy

• Maximum Entropy Principle

• Chemical Potential

• Ensemble Formula Sheet

• Many-Body and Quantum Statistical Mechanics

• Thermodynamic Limit

• Density Operators

• Trace Rule for Expectation Values

• Entropy Overview

• Relative Entropy

• Euclidean and Imaginary-Time Path Integrals

• Matsubara Formalism Preview

• Quantum Thermodynamics

• Thermal Master Equations

• Strong-Coupling Open-System Effects

• Statistical Mechanics Checklist

References

• R. K. Pathria and P. D. Beale, Statistical Mechanics, 3rd ed., Elsevier (2011).
• K. Huang, Statistical Mechanics, 2nd ed., Wiley (1987).
• M. Kardar, Statistical Physics of Particles, Cambridge University Press (2007).
• J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press (2020), chapter on identical particles and statistical mechanics.
• J. von Neumann, Mathematical Foundations of Quantum Mechanics, Princeton University Press (1955).
• H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).
• H. B. Callen and T. A. Welton, “Irreversibility and generalized noise,” Physical Review 83, 34–40 (1951).