Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2601.03253 · Journal of Statistical Physics

Grand-Canonical Typicality

Cedric Igelspacher, Roderich Tumulka, Cornelia Vogel

quant-phSubmitted Jan 6, 2026• Announced Jan 7, 2026• 0 citations
v1 · Jan 6, 2026v2 · Jan 24, 2026v3 · May 21, 2026
arXivHTMLPDF

Abstract

We study how the grand-canonical density matrix arises in macroscopic quantum systems. ``Canonical typicality'' is the known statement that for a typical wave function $Ψ$ from a micro-canonical energy shell of a quantum system $S$ weakly coupled to a large but finite quantum system $B$, the reduced density matrix $\hatρ^S_Ψ=\mathrm{tr}^B |Ψ\rangle\langle Ψ|$ is approximately equal to the canonical density matrix $\hatρ_\mathrm{can}=Z^{-1}_\mathrm{can} \exp(-β\hat{H}^S)$. Here, we discuss the analogous statement and related questions for the \emph{grand-canonical} density matrix $\hatρ_\mathrm{gc}=Z^{-1}_\mathrm{gc} \exp(-β(\hat{H}^S-μ_1 \hat{N}_{1}^S-\ldots-μ_r\hat{N}_{r}^S))$ with $\hat{N}_{i}^S$ the number operator for molecules of type $i$ in the system $S$. This includes (i) the case of chemical reactions and (ii) that of systems $S$ defined by a spatial region which particles may enter or leave. It includes the statements (a) that the density matrix of the appropriate (generalized micro-canonical) Hilbert subspace $H_\mathrm{gmc} \subset H^S \otimes H^B$ (defined by a micro-canonical interval of total energy and suitable particle number sectors), after tracing out $B$, yields $\hatρ_\mathrm{gc}$; (b) that typical $Ψ$ from $H_\mathrm{gmc}$ have reduced density matrix $\hatρ^S_Ψ$ close to $\hatρ_\mathrm{gc}$; and (c) that the conditional wave function $ψ^S$ of $S$ has probability distribution $\mathrm{GAP}_{\hatρ_\mathrm{gc}}$ if a typical orthonormal basis of $H^B$ is used. That is, we discuss the foundation and justification of both the density matrix and the distribution of the wave function in the grand-canonical case. We also extend these considerations to the so-called generalized Gibbs ensembles, which apply to systems for which some macroscopic observables are conserved.

Journal reference
Journal of Statistical Physics 193: 120 (2026)
DOI
10.1007/s10955-026-03656-5
Comments
39 pages LaTeX, no figures

Related papers

Canonical Typicality

Sheldon Goldstein, Joel L. Lebowitz, Roderich Tumulka +1 more

It is well known that a system weakly coupled to a heat bath is described by the canonical ensemble when the composite S + B is described by the microcanonical ensemble corresponding to a suitable energy shell. This is true for both classical distributions on the phase space and quantum density matrices. Here we show that a much stronger statement holds for quantum systems. Even if the state of the composite corresponds to a single wave function rather than a mixture, the reduced density matrix of the system is canonical, for the overwhelming majority of wave functions in the subspace corresponding to the energy interval encompassed by the microcanonical ensemble. This clarifies, expands, and justifies remarks made by Schrödinger in 1952.

arXiv:cond-mat/0511091•Physical Review Letters•Feb 8, 2006•711 citations

Canonical Typicality for Other Ensembles than Micro-canonical

Stefan Teufel, Roderich Tumulka, Cornelia Vogel

Abstract We generalize Lévy’s lemma, a concentration-of-measure result for the uniform probability distribution on high-dimensional spheres, to a much more general class of measures, so-called GAP measures. For any given density matrix $$\rho $$ ρ on a separable Hilbert space $${\mathcal {H}}$$ H , $${\textrm{GAP}}(\rho )$$ GAP ( ρ ) is the most spread-out probability measure on the unit sphere of $${\mathcal {H}}$$ H that has density matrix $$\rho $$ ρ and thus forms the natural generalization of the uniform distribution. We prove concentration-of-measure whenever the largest eigenvalue $$\Vert \rho \Vert $$ ‖ ρ ‖ of $$\rho $$ ρ is small. We use this fact to generalize and improve well-known and important typicality results of quantum statistical mechanics to GAP measures, namely canonical typicality and dynamical typicality. Canonical typicality is the statement that for “most” pure states $$\psi $$ ψ of a given ensemble, the reduced density matrix of a sufficiently small subsystem is very close to a $$\psi $$ ψ -independent matrix. Dynamical typicality is the statement that for any observable and any unitary time evolution, for “most” pure states $$\psi $$ ψ from a given ensemble the (coarse-grained) Born distribution of that observable in the time-evolved state $$\psi _t$$ ψ t is very close to a $$\psi $$ ψ -independent distribution. So far, canonical typicality and dynamical typicality were known for the uniform distribution on finite-dimensional spheres, corresponding to the micro-canonical ensemble, and for rather special mean-value ensembles. Our result shows that these typicality results hold also for $${\textrm{GAP}}(\rho )$$ GAP ( ρ ) , provided the density matrix $$\rho $$ ρ has small eigenvalues. Since certain GAP measures are quantum analogs of the canonical ensemble of classical mechanics, our results can also be regarded as a version of equivalence of ensembles.

Annales Henri Poincaré•Jul 2, 2024•7 citations

Universal fluctuations around typicality for quantum ergodic systems

Michel Bauer, Denis Bernard, Tony Jin

For a quantum system in a macroscopically large volume V, prepared in a pure state and subject to maximally noisy or ergodic unitary dynamics, the reduced density matrix of any sub-system v≪V is almost surely totally mixed. We show that the fluctuations around this limiting value, evaluated according to the invariant measure of these unitary flows, are captured by the Gaussian unitary ensemble (GUE) of random matrix theory. An extension of this statement, applicable when the unitary transformations conserve the energy but are maximally noisy or ergodic on any energy shell, allows to decipher the fluctuations around canonical typicality. According to typicality, if the large system is prepared in a generic pure state in a given energy shell, the reduced density matrix of the sub-system is almost surely the canonical Gibbs state of that sub-system. We show that the fluctuations around the Gibbs state are encoded in a deformation of the GUE whose covariance is specified by the Gibbs state. Contact with the eigenstate thermalization hypothesis is discussed.

arXiv:1907.08081•Physical review. E•Jan 10, 2020•18 citations

Universal fluctuations of reduced density matrices in maximally noisy or ergodic quantum systems and typicality

Michel Bauer, Denis Bernard, Tony Jin

For a quantum system in a macroscopically large volume $V$, prepared in a pure state and subject to maximally noisy or ergodic unitary dynamics, the reduced density matrix of any sub-system $v\ll V$ is almost surely totally mixing. We show that the fluctuations around this limiting value, evaluated according to the invariant measure of these unitary flows, are captured by the Gaussian unitary ensemble (GUE) of random matrix theory. An extension of this statement, applicable when the unitary transformations conserve the energy but are maximally noisy or ergodic on any energy shell, allows to decipher the fluctuations around canonical typicality. According to typicality, if the large system is prepared in a generic pure state in a given energy shell, the reduced density matrix of the sub-system is almost surely the canonical Gibbs state of that sub-system. We show that the fluctuations around the Gibbs state are encoded in a deformation of the GUE whose covariance is specified by the Gibbs state. Contact with the eigenstate thermal hypothesis (ETH) is discussed.

•Jul 18, 2019•0 citations

Universal Probability Distribution for the Wave Function of an Open Quantum System

Sheldon Goldstein, Joel L. Lebowitz, Christian Mastrodonato +2 more

An open quantum system (i.e., one that interacts with its environment) is almost always entangled with its environment; it is therefore usually not attributed a wave function but only a reduced density matrix ρ. Nevertheless, there is a precise way of attributing to it a wave function ψ1, called its conditional wave function, which is a random wave function of the system whose probability distribution µ1 depends on the entangled wave function ψ ∈ H1 ⊗ H2 in the Hilbert space of system and environment together. We prove several universality (or typicality) results about µ1; they show that if the environment is sufficiently large then µ1 does not depend much on the details of ψ and is approximately given by one of the so-called GAP measures. Specifically, for most entangled states ψ with given reduced density matrix ρ1, µ1 is close to GAP (ρ1). We also show that, if the coupling between the system and the environment is weak, then for most entangled states ψ from a microcanonical subspace corresponding to energies in a narrow interval [E, E + δE] (and most bases of H2), µ1 is close to GAP (ρβ) with ρβ the canonical density matrix on H1 at inverse temperature β = β(E). This provides the mathematical justification of the claim that GAP (ρβ) is the thermal equilibrium distribution of ψ1.

•Apr 28, 2011•2 citations

Open quantum systems and the grand canonical ensemble

Benedikt M. Reible, Luigi Delle Site

The celebrated Lindblad equation governs the nonunitary time evolution of density operators used in the description of open quantum systems. It is usually derived from the von Neumann equation for a large system, at given physical conditions, when a small subsystem is explicitly singled out and the rest of the system acts as an environment whose degrees of freedom are traced out. In the specific case of a subsystem with variable particle number, the equilibrium density operator is given by the well-known grand canonical Gibbs state. Consequently, solving the Lindblad equation in this case should automatically yield, without any additional assumptions, the corresponding density operator in the limiting case of statistical equilibrium. Current studies of the Lindblad equation with varying particle number assume, however, the grand canonical Gibbs state a priori: the chemical potential is externally imposed rather than derived from first principles, and hence the corresponding density operator is not obtained as a natural solution of the equation. In this work, we investigate the compatibility of grand canonical statistical mechanics with the derivation of the Lindblad equation. We propose an alternative and complementary approach to the current literature that consists in using a generalized system Hamiltonian which includes a term μN. In a previous paper, this empirically well-known term has been formally derived from the von Neumann equation for the specific case of equilibrium. Including μN in the system Hamiltonian leads to a modified Lindblad equation which yields the grand canonical state as a natural solution, meaning that all the quantities involved are obtained from the physics of the system without any external assumptions.

arXiv:2508.16985•Physical review. E•Aug 7, 2025•3 citations
We study how the grand-canonical density matrix arises in macroscopic quantum systems. ``Canonical typicality'' is the known statement that for a typical wave function ΨΨΨ from a micro-canonical energy shell of a quantum system SSS weakly coupled to a large but finite quantum system BBB, the reduced density matrix ρ^ΨS=trB∣Ψ⟩⟨Ψ∣\hatρ^S_Ψ=\mathrm{tr}^B |Ψ\rangle\langle Ψ|ρ^​ΨS​=trB∣Ψ⟩⟨Ψ∣ is approximately equal to the canonical density matrix ρ^can=Zcan−1exp⁡(−βH^S)\hatρ_\mathrm{can}=Z^{-1}_\mathrm{can} \exp(-β\hat{H}^S)ρ^​can​=Zcan−1​exp(−βH^S). Here, we discuss the analogous statement and related questions for the \emph{grand-canonical} density matrix ρ^gc=Zgc−1exp⁡(−β(H^S−μ1N^1S−…−μrN^rS))\hatρ_\mathrm{gc}=Z^{-1}_\mathrm{gc} \exp(-β(\hat{H}^S-μ_1 \hat{N}_{1}^S-\ldots-μ_r\hat{N}_{r}^S))ρ^​gc​=Zgc−1​exp(−β(H^S−μ1​N^1S​−…−μr​N^rS​)) with N^iS\hat{N}_{i}^SN^iS​ the number operator for molecules of type iii in the system SSS. This includes (i) the case of chemical reactions and (ii) that of systems SSS defined by a spatial region which particles may enter or leave. It includes the statements (a) that the density matrix of the appropriate (generalized micro-canonical) Hilbert subspace Hgmc⊂HS⊗HBH_\mathrm{gmc} \subset H^S \otimes H^BHgmc​⊂HS⊗HB (defined by a micro-canonical interval of total energy and suitable particle number sectors), after tracing out BBB, yields ρ^gc\hatρ_\mathrm{gc}ρ^​gc​; (b) that typical ΨΨΨ from HgmcH_\mathrm{gmc}Hgmc​ have reduced density matrix ρ^ΨS\hatρ^S_Ψρ^​ΨS​ close to ρ^gc\hatρ_\mathrm{gc}ρ^​gc​; and (c) that the conditional wave function ψSψ^SψS of SSS has probability distribution GAPρ^gc\mathrm{GAP}_{\hatρ_\mathrm{gc}}GAPρ^​gc​​ if a typical orthonormal basis of HBH^BHB is used. That is, we discuss the foundation and justification of both the density matrix and the distribution of the wave function in the grand-canonical case. We also extend these considerations to the so-called generalized Gibbs ensembles, which apply to systems for which some macroscopic observables are conserved.
It is well known that a system weakly coupled to a heat bath is described by the canonical ensemble when the composite S + B is described by the microcanonical ensemble corresponding to a suitable energy shell. This is true for both classical distributions on the phase space and quantum density matrices. Here we show that a much stronger statement holds for quantum systems. Even if the state of the composite corresponds to a single wave function rather than a mixture, the reduced density matrix of the system is canonical, for the overwhelming majority of wave functions in the subspace corresponding to the energy interval encompassed by the microcanonical ensemble. This clarifies, expands, and justifies remarks made by Schrödinger in 1952.
Abstract We generalize Lévy’s lemma, a concentration-of-measure result for the uniform probability distribution on high-dimensional spheres, to a much more general class of measures, so-called GAP measures. For any given density matrix ρ\rho ρ ρ on a separable Hilbert space H{\mathcal {H}}H H , GAP(ρ){\textrm{GAP}}(\rho )GAP(ρ) GAP ( ρ ) is the most spread-out probability measure on the unit sphere of H{\mathcal {H}}H H that has density matrix ρ\rho ρ ρ and thus forms the natural generalization of the uniform distribution. We prove concentration-of-measure whenever the largest eigenvalue ∥ρ∥\Vert \rho \Vert ∥ρ∥ ‖ ρ ‖ of ρ\rho ρ ρ is small. We use this fact to generalize and improve well-known and important typicality results of quantum statistical mechanics to GAP measures, namely canonical typicality and dynamical typicality. Canonical typicality is the statement that for “most” pure states ψ\psi ψ ψ of a given ensemble, the reduced density matrix of a sufficiently small subsystem is very close to a ψ\psi ψ ψ -independent matrix. Dynamical typicality is the statement that for any observable and any unitary time evolution, for “most” pure states ψ\psi ψ ψ from a given ensemble the (coarse-grained) Born distribution of that observable in the time-evolved state ψt\psi _tψt​ ψ t is very close to a ψ\psi ψ ψ -independent distribution. So far, canonical typicality and dynamical typicality were known for the uniform distribution on finite-dimensional spheres, corresponding to the micro-canonical ensemble, and for rather special mean-value ensembles. Our result shows that these typicality results hold also for GAP(ρ){\textrm{GAP}}(\rho )GAP(ρ) GAP ( ρ ) , provided the density matrix ρ\rho ρ ρ has small eigenvalues. Since certain GAP measures are quantum analogs of the canonical ensemble of classical mechanics, our results can also be regarded as a version of equivalence of ensembles.
For a quantum system in a macroscopically large volume V, prepared in a pure state and subject to maximally noisy or ergodic unitary dynamics, the reduced density matrix of any sub-system v≪V is almost surely totally mixed. We show that the fluctuations around this limiting value, evaluated according to the invariant measure of these unitary flows, are captured by the Gaussian unitary ensemble (GUE) of random matrix theory. An extension of this statement, applicable when the unitary transformations conserve the energy but are maximally noisy or ergodic on any energy shell, allows to decipher the fluctuations around canonical typicality. According to typicality, if the large system is prepared in a generic pure state in a given energy shell, the reduced density matrix of the sub-system is almost surely the canonical Gibbs state of that sub-system. We show that the fluctuations around the Gibbs state are encoded in a deformation of the GUE whose covariance is specified by the Gibbs state. Contact with the eigenstate thermalization hypothesis is discussed.
For a quantum system in a macroscopically large volume VVV, prepared in a pure state and subject to maximally noisy or ergodic unitary dynamics, the reduced density matrix of any sub-system v≪Vv\ll Vv≪V is almost surely totally mixing. We show that the fluctuations around this limiting value, evaluated according to the invariant measure of these unitary flows, are captured by the Gaussian unitary ensemble (GUE) of random matrix theory. An extension of this statement, applicable when the unitary transformations conserve the energy but are maximally noisy or ergodic on any energy shell, allows to decipher the fluctuations around canonical typicality. According to typicality, if the large system is prepared in a generic pure state in a given energy shell, the reduced density matrix of the sub-system is almost surely the canonical Gibbs state of that sub-system. We show that the fluctuations around the Gibbs state are encoded in a deformation of the GUE whose covariance is specified by the Gibbs state. Contact with the eigenstate thermal hypothesis (ETH) is discussed.
An open quantum system (i.e., one that interacts with its environment) is almost always entangled with its environment; it is therefore usually not attributed a wave function but only a reduced density matrix ρ. Nevertheless, there is a precise way of attributing to it a wave function ψ1, called its conditional wave function, which is a random wave function of the system whose probability distribution µ1 depends on the entangled wave function ψ ∈ H1 ⊗ H2 in the Hilbert space of system and environment together. We prove several universality (or typicality) results about µ1; they show that if the environment is sufficiently large then µ1 does not depend much on the details of ψ and is approximately given by one of the so-called GAP measures. Specifically, for most entangled states ψ with given reduced density matrix ρ1, µ1 is close to GAP (ρ1). We also show that, if the coupling between the system and the environment is weak, then for most entangled states ψ from a microcanonical subspace corresponding to energies in a narrow interval [E, E + δE] (and most bases of H2), µ1 is close to GAP (ρβ) with ρβ the canonical density matrix on H1 at inverse temperature β = β(E). This provides the mathematical justification of the claim that GAP (ρβ) is the thermal equilibrium distribution of ψ1.
The celebrated Lindblad equation governs the nonunitary time evolution of density operators used in the description of open quantum systems. It is usually derived from the von Neumann equation for a large system, at given physical conditions, when a small subsystem is explicitly singled out and the rest of the system acts as an environment whose degrees of freedom are traced out. In the specific case of a subsystem with variable particle number, the equilibrium density operator is given by the well-known grand canonical Gibbs state. Consequently, solving the Lindblad equation in this case should automatically yield, without any additional assumptions, the corresponding density operator in the limiting case of statistical equilibrium. Current studies of the Lindblad equation with varying particle number assume, however, the grand canonical Gibbs state a priori: the chemical potential is externally imposed rather than derived from first principles, and hence the corresponding density operator is not obtained as a natural solution of the equation. In this work, we investigate the compatibility of grand canonical statistical mechanics with the derivation of the Lindblad equation. We propose an alternative and complementary approach to the current literature that consists in using a generalized system Hamiltonian which includes a term μN. In a previous paper, this empirically well-known term has been formally derived from the von Neumann equation for the specific case of equilibrium. Including μN in the system Hamiltonian leads to a modified Lindblad equation which yields the grand canonical state as a natural solution, meaning that all the quantities involved are obtained from the physics of the system without any external assumptions.