---
title: Hilbert–Schmidt Ensemble Overview
url: https://www.emergentmind.com/topics/hilbert-schmidt-ensemble
type: topic
---

# Hilbert–Schmidt Ensemble Overview

Searching arXiv for relevant papers on Hilbert–Schmidt ensembles, induced measures, and related random density matrix results.
Using arXiv search to retrieve recent and foundational papers relevant to the Hilbert–Schmidt ensemble.
The Hilbert–Schmidt ensemble is the unitarily invariant, flat Euclidean measure on the convex set of density matrices, defined by the Hilbert–Schmidt metric
\[
d_{HS}(\rho,\sigma)=\sqrt{\operatorname{Tr}\big[(\rho-\sigma)^2\big]}.
\]
In random-matrix and quantum-information formulations, it is the square case of the induced ensemble \(\mu_{N,K}\), obtained either by partial tracing a Haar-random pure state on \(\mathbb{C}^N\otimes \mathbb{C}^K\) or, equivalently, by normalizing a Wishart matrix \(W=GG^\dagger\) built from a Ginibre matrix \(G\). The Hilbert–Schmidt specialization corresponds to \(K=N\), so that
\[
\rho=\frac{GG^\dagger}{\operatorname{Tr}(GG^\dagger)}
\]
is distributed according to the fixed-trace Wishart–Laguerre measure and provides the canonical “structureless” ensemble of random mixed states [1010.3570, 2008.05153].

## 1. Definition as an induced and fixed-trace Wishart ensemble

For a bipartite Hilbert space \(\mathcal H_A\otimes\mathcal H_B\) of dimensions \(N\times K\), one may draw a Haar-random pure state \(|\psi\rangle\in\mathcal H_A\otimes\mathcal H_B\) and define the reduced state
\[
\rho_A=\operatorname{Tr}_B\big(|\psi\rangle\langle\psi|\big).
\]
This construction yields the induced measures \(\mu_{N,K}\) on the set \(\mathcal M_N\) of \(N\times N\) density matrices. The Hilbert–Schmidt ensemble is the square case \(K=N\). An equivalent formulation expands \(|\psi\rangle\) in a product basis with coefficients \(X_{ij}\), arranges them into a matrix \(G=(X_{ij})\in\mathbb C^{N\times K}\), and sets \(W=GG^\dagger\). The reduced state then has normalized Wishart form
\[
\rho=\frac{W}{\operatorname{Tr}W}=\frac{GG^\dagger}{\operatorname{Tr}(GG^\dagger)}.
\]
This equivalence is central both conceptually and computationally [1010.3570].

In the formulation used for exact distance calculations, the induced fixed-trace measure on density matrices is
\[
P(\rho)=C(\det \rho)^\alpha \delta(\operatorname{Tr}\rho-1),
\qquad
\alpha=\frac{\beta}{2}(K-N+1)-1,
\]
where \(N\) is the system dimension, \(K\) is the environment or ancilla dimension, and \(\beta\in\{1,2\}\) is the Dyson index specifying real-symmetric (\(\beta=1\)) or complex-Hermitian (\(\beta=2\)) ensembles. The associated Wishart construction uses
\[
W=GG^\dagger,
\]
with \(G\) an \(N\times K\) real or complex Ginibre matrix distributed as
\[
P_G(G)=\left(\frac{\beta}{2\pi}\right)^{\beta NK/2}\exp\!\left(-\frac{\beta}{2}\operatorname{Tr}(GG^\dagger)\right).
\]
With this convention, the real variance of each independent entry is \(1\) for \(\beta=1\), and the complex variance is \(1\) per complex component for \(\beta=2\) [2008.05153].

Two equivalent sampling procedures therefore generate Hilbert–Schmidt random density matrices. The first is the Ginibre/Wishart recipe: sample \(G\in\mathbb C^{N\times N}\) with i.i.d. complex normal entries, set \(W=GG^\dagger\), and output \(\rho=W/\operatorname{Tr}W\). The second is the Haar-pure-state construction: sample \(U\in U(NK)\), set \(|\psi\rangle=U|0,0\rangle\), compute \(\rho_A=\operatorname{Tr}_B(|\psi\rangle\langle\psi|)\), and output \(\rho_A\). Proposition 1 in the structured-ensemble framework further shows stability under partial traces: tracing out a factor from an induced ensemble \(\mu_{M,K}\) yields \(\mu_{N,KL}\) when \(M=NL\) [1010.3570].

## 2. Spectral form, unitary invariance, and basic moments

The induced and Hilbert–Schmidt ensembles are unitarily invariant. In eigenvalue–eigenvector coordinates, the measure factorizes into Haar measure on eigenvectors and a joint density on eigenvalues. In the complex case \((\beta=2)\), the induced ensemble has joint eigenvalue density
\[
P(\{\lambda_i\})\propto
\delta\!\left(\sum_{i=1}^N \lambda_i-1\right)
\prod_{i=1}^N \lambda_i^{K-N}
\prod_{1\le i<j\le N}(\lambda_i-\lambda_j)^2.
\]
For the Hilbert–Schmidt ensemble, \(K=N\), so the factor \(\prod_i \lambda_i^{K-N}\) becomes unity. Equivalently, in the matrix formulation, the dependence on a fixed comparator reduces to spectral invariants because of unitary invariance [1010.3570, 2008.05153].

The ensemble averages that organize many exact calculations are especially simple. One has
\[
E[\rho]=\frac{I_N}{N},
\]
and the average purity is
\[
E[\operatorname{Tr}(\rho^2)]
=
\frac{\beta(N+K+2/\beta-1)}{\beta NK+2}.
\]
For \(\beta=2\),
\[
E[\operatorname{Tr}(\rho^2)]=\frac{N+K}{NK+1},
\]
and in the Hilbert–Schmidt specialization \(K=N\),
\[
E[\operatorname{Tr}(\rho^2)]_{\mathrm{HS}}=\frac{2N}{N^2+1}.
\]
If \(\sigma\) is fixed with \(\operatorname{Tr}\sigma=1\), then
\[
E[\operatorname{Tr}(\rho\sigma)]=\frac{1}{N},
\]
and for two independent induced density matrices \(\rho_1,\rho_2\),
\[
E[\operatorname{Tr}(\rho_1\rho_2)]=\frac{1}{N}.
\]
These identities are the basis of the compact Hilbert–Schmidt distance formulas [2008.05153].

In the large-dimension regime \(N\to\infty\) with fixed aspect ratio \(c=K/N\), the empirical eigenvalue density of the rescaled eigenvalues \(x=N\lambda\) converges to the Marchenko–Pastur law
\[
\pi^{(1)}_c(x)
=
\max(1-c,0)\,\delta(0)
+
\frac{\sqrt{4c-(x-c-1)^2}}{2\pi x},
\qquad
x\in[x_-,x_+],
\quad
x_\pm=1+c\pm 2\sqrt c.
\]
For the Hilbert–Schmidt case \(c=1\), this reduces to
\[
p_{\mathrm{HS}}(\mu)=\frac{1}{2\pi}\sqrt{\frac{4}{\mu}-1},
\qquad
\mu\in[0,4].
\]
The average von Neumann entropy of the reduced state for a Haar-random pure state on \(N\times K\) is
\[
\mathbb E[S(\rho_A)]
=
\sum_{j=K+1}^{NK}\frac{1}{j}
-
\frac{N-1}{2K},
\]
so that in the Hilbert–Schmidt specialization
\[
\mathbb E[S(\rho)]_{\mathrm{HS}}
=
\ln N-\frac{1}{2}+O\!\left(\frac{\ln N}{N}\right).
\]
This suggests that Hilbert–Schmidt random states are typically highly mixed, with entropy close to \(\ln N\) but purity of order \(1/N\) [1010.3570].

## 3. Exact mean-square Hilbert–Schmidt distances

The Hilbert–Schmidt distance between two density matrices is
\[
D_{HS}(\rho,\sigma)=\|\rho-\sigma\|_2=\sqrt{\operatorname{Tr}\big[(\rho-\sigma)^2\big]},
\]
and its squared form is
\[
D_{HS}^2(\rho,\sigma)=\operatorname{Tr}(\rho^2)+\operatorname{Tr}(\sigma^2)-2\operatorname{Tr}(\rho\sigma).
\]
The exact results available for the Hilbert–Schmidt ensemble concern the mean of the squared distance, \(E[D_{HS}^2]\), not the mean of the unsquared distance \(E[D_{HS}]\). Since \(E[\sqrt{D_{HS}^2}]\neq \sqrt{E[D_{HS}^2]}\) in general, this distinction is essential [2008.05153].

For a random density matrix \(\rho\) from the induced ensemble and a fixed density matrix \(\sigma\), the exact mean-square Hilbert–Schmidt distance is
\[
E[D_{HS}^2(\rho,\sigma)]
=
\operatorname{Tr}(\sigma^2)
+
\frac{\beta(N+K+2/\beta-1)}{\beta NK+2}
-
\frac{2}{N}.
\]
Equivalently,
\[
E[D_{HS}^2(\rho,\sigma)]
=
E[\operatorname{Tr}(\rho^2)]
+\operatorname{Tr}(\sigma^2)-\frac{2}{N}.
\]
For \(\beta=2\),
\[
E[D_{HS}^2(\rho,\sigma)]
=
\operatorname{Tr}(\sigma^2)+\frac{N+K}{NK+1}-\frac{2}{N},
\]
and in the Hilbert–Schmidt case \(K=N\),
\[
E[D_{HS}^2(\rho,\sigma)]_{\mathrm{HS}}
=
\operatorname{Tr}(\sigma^2)+\frac{2N}{N^2+1}-\frac{2}{N}.
\]
If \(\sigma\) is pure, \(\operatorname{Tr}(\sigma^2)=1\); if \(\sigma=I/N\) is maximally mixed, \(\operatorname{Tr}(\sigma^2)=1/N\). A notable structural feature is that no spectral detail of \(\sigma\) beyond \(\operatorname{Tr}(\sigma^2)\) enters the final formula [2008.05153].

For two independent induced density matrices \(\rho_1,\rho_2\), possibly with different ancilla dimensions \(K_1,K_2\),
\[
E[D_{HS}^2(\rho_1,\rho_2)]
=
\frac{\beta(N+K_1+2/\beta-1)}{\beta NK_1+2}
+
\frac{\beta(N+K_2+2/\beta-1)}{\beta NK_2+2}
-
\frac{2}{N}.
\]
In the complex Hilbert–Schmidt specialization,
\[
E[D_{HS}^2(\rho_1,\rho_2)]_{\mathrm{HS}}
=
\frac{4N}{N^2+1}-\frac{2}{N}.
\]
For large \(N\), this is asymptotically \(2/N+O(1/N^2)\), matching known large-\(N\) behavior [2008.05153].

The same paper derives corresponding exact formulas for non-normalized Wishart matrices. If \(W\) is Wishart with \(K\) degrees of freedom and \(H\) is a fixed Hermitian matrix,
\[
E[\|W-H\|_2^2]
=
NK(N+K+2/\beta-1)+\operatorname{Tr}(H^2)-2K\operatorname{Tr}(H),
\]
using
\[
E[\operatorname{Tr}(W^2)]=NK(N+K+2/\beta-1),
\qquad
E[\operatorname{Tr}(WH)]=K\operatorname{Tr}(H).
\]
For two Wishart matrices \(W_1,W_2\),
\[
E[\|W_1-W_2\|_2^2]
=
N\Big[(K_1+K_2)(N+2/\beta-1)+(K_1-K_2)^2\Big].
\]
These Wishart results are mapped to the fixed-trace formulas through an exact Laplace-transform procedure: one introduces an auxiliary variable in \(\delta(\operatorname{Tr}\rho-t)\), takes the Laplace transform, rescales \(\rho=(\beta/2s)W\), and then inverts the transform at \(t=1\). No Weingarten calculus is required; unitary invariance and Gaussian/Wishart moments suffice [2008.05153].

## 4. Asymptotic regimes and relation to other invariant ensembles

In the Hilbert–Schmidt ensemble \((\beta=2, K=N)\), the average purity satisfies
\[
E[\operatorname{Tr}(\rho^2)]
=
\frac{2N}{N^2+1}
=
\frac{2}{N}+O\!\left(\frac{1}{N^3}\right).
\]
Consequently, for a fixed comparator \(\sigma\),
\[
E[D_{HS}^2(\rho,\sigma)]
=
\operatorname{Tr}(\sigma^2)
+
\frac{2N}{N^2+1}
-
\frac{2}{N}
=
\operatorname{Tr}(\sigma^2)+O\!\left(\frac{1}{N^2}\right).
\]
If \(\sigma\) is pure, the mean-square distance approaches \(1\); if \(\sigma=I/N\), it behaves as \(1/N+O(1/N^2)\) and therefore vanishes as \(N\to\infty\). For two independent Hilbert–Schmidt states,
\[
E[D_{HS}^2(\rho_1,\rho_2)]
=
\frac{2}{N}+O\!\left(\frac{1}{N^2}\right).
\]
A plausible implication is that high-dimensional Hilbert–Schmidt states are typically close to maximally mixed and only weakly separated from each other in Hilbert–Schmidt geometry [2008.05153].

The Hilbert–Schmidt ensemble is also the baseline from which several structured invariant ensembles are defined. In the broader family
\[
W_{k,s}=\big(p_1U_1+\cdots+p_kU_k\big)\,G_1\cdots G_s,
\qquad
\rho=\frac{W_{k,s}W_{k,s}^\dagger}{\operatorname{Tr}(W_{k,s}W_{k,s}^\dagger)},
\]
the choice \(k=1\), \(s=1\) produces the Hilbert–Schmidt ensemble, whereas \(k=2\), \(s=1\), \(p_1=p_2=\tfrac12\) yields the Bures ensemble. This places Hilbert–Schmidt measure at the simplest point in a hierarchy of unitarily invariant random-state models [1010.3570].

The distinction from the Bures ensemble is explicit at the level of eigenvalue statistics. The Bures eigenvalue joint density is
\[
P_B(\lambda_1,\dots,\lambda_N)
=
C_N^B\,
\delta\!\Big(\sum_{i=1}^N\lambda_i-1\Big)\,
\prod_{i=1}^N \lambda_i^{-1/2}\,
\prod_{i<j}\frac{(\lambda_i-\lambda_j)^2}{\lambda_i+\lambda_j},
\]
whereas the Hilbert–Schmidt density lacks both the \(\lambda_i^{-1/2}\) factor and the \((\lambda_i+\lambda_j)^{-1}\) denominator. The constructive Bures algorithm correspondingly uses both a Ginibre matrix and a Haar unitary,
\[
\rho_{\mathrm B}
=
\frac{(\mathbb 1+U)GG^\dagger(\mathbb 1+U^\dagger)}
{\operatorname{Tr}\big[(\mathbb 1+U)GG^\dagger(\mathbb 1+U^\dagger)\big]}.
\]
This suggests that the Hilbert–Schmidt ensemble should be regarded not merely as another unitary-invariant measure, but as the Euclidean benchmark against which more strongly purity-biased monotone-metric ensembles are compared [1010.3570].

Beyond Bures, the structured-ensemble paper describes superpositions of \(k\) random maximally entangled states and product-Ginibre constructions obtained from selective measurements in maximally entangled bases. Their limiting spectral laws include the arcsine law, the family \(\nu_k\), and the Fuss–Catalan distributions \(\pi^{(s)}\), while the Hilbert–Schmidt ensemble reappears as the Marchenko–Pastur baseline at \(s=1\) [1010.3570].

## 5. Low-dimensional determinantal moments and separability

For \(2\times 2\) bipartite systems, the Hilbert–Schmidt ensemble has a particularly explicit determinantal theory. The two-rebit state space is \(9\)-dimensional, the two-qubit state space is \(15\)-dimensional, and the Hilbert–Schmidt measure is the unitarily invariant flat measure on these convex bodies. In eigenvalue variables \(\lambda_i\ge 0\), \(\sum_i\lambda_i=1\), the Hilbert–Schmidt density is proportional to
\[
P^{HS}_\alpha(\lambda)\propto \prod_{i<j}|\lambda_i-\lambda_j|^{2\alpha},
\]
with \(\alpha=\tfrac12\) for two-rebits and \(\alpha=1\) for two-qubits. The normalization constants of the eigenvalue-simplex densities are \(80640\) for \(\alpha=\tfrac12\) and \(378378000\) for \(\alpha=1\) [1207.1297].

In \(2\times 2\), the Peres–Horodecki criterion is necessary and sufficient: a state is separable if and only if its partial transpose is positive semidefinite, equivalently if and only if \(|\rho^{PT}|\ge 0\). The determinant range is
\[
|\rho^{PT}|\in\left[-\frac{1}{16},\frac{1}{256}\right],
\]
while the “balanced” variable satisfies
\[
|\rho^{PT}|\cdot |\rho|
\in
\left[-\frac{1}{110592},\frac{1}{256^2}\right].
\]
These bounds permit moment-based reconstruction of the \(|\rho^{PT}|\) distribution under Hilbert–Schmidt measure [1207.1297].

Exact determinantal moments are available in closed form. For two-qubits \((\alpha=1)\),
\[
\big\langle |\rho|^k\big\rangle^{HS}_{\alpha=1}
=
\frac{108972864000\;
\Gamma(k+1)\Gamma(k+2)\Gamma(k+3)\Gamma(k+4)}
{\Gamma\!\big(4(k+4)\big)}.
\]
More generally, the mixed ratios
\[
R_n(k;\alpha)
=
\frac{\langle |\rho^{PT}|^n |\rho|^k\rangle^{HS}_\alpha}
{\langle |\rho|^k\rangle^{HS}_\alpha}
\]
are rational functions of \(k\) that can be written as ratios of degree-\(3n\) polynomials. For \(n=1\), the two-qubit and two-rebit cases are
\[
R_1(k;1)=
\frac{k^3+6k^2-k-42}{256k^3+3456k^2+15536k+23256}
=
\frac{k\big(k(k+6)-1\big)-42}{8(2k+9)(4k+17)(4k+19)},
\]
and
\[
R_1\!\big(k;\tfrac12\big)
=
\frac{(k-1)\big(k(2k+11)+16\big)}
{32(k+3)(4k+11)(4k+13)}.
\]
The first Hilbert–Schmidt moments of \(|\rho^{PT}|\) are
\[
\langle |\rho^{PT}|\rangle^{HS}_{\alpha=1}=-\frac{7}{3876},
\qquad
\langle |\rho^{PT}|\rangle^{HS}_{\alpha=1/2}=-\frac{1}{858}.
\]
For the balanced variable,
\[
\langle |\rho^{PT}|\cdot |\rho| \rangle^{HS}_{\alpha=1/2}=0,
\qquad
\langle |\rho^{PT}|\cdot |\rho| \rangle^{HS}_{\alpha=1}
=
-\frac{1}{4576264}.
\]
A central contrast is that the analogous Bures ratios appear as degree-\(5n\) polynomials in \(k\), not degree-\(3n\) [1207.1297].

These exact Hilbert–Schmidt moments support inverse moment reconstructions of the \(|\rho^{PT}|\) density, for example via Legendre-polynomial reconstruction of Provost. In this way one obtains exact Hilbert–Schmidt separability probabilities:
\[
P^{HS}_{sep}=\frac{29}{64}\quad\text{for two-rebits},
\qquad
P^{HS}_{sep}=\frac{8}{33}\quad\text{for two-qubits},
\qquad
P^{HS}_{sep}=\frac{26}{323}\quad\text{for two-quaterbits}.
\]
This makes the Hilbert–Schmidt ensemble a quantitatively tractable reference measure for entanglement geometry in low-dimensional bipartite state spaces [1207.1297].

## 6. Qutrit Hilbert–Schmidt geometry and classicality indicators

For qutrits, the Hilbert–Schmidt ensemble can be developed directly from the metric-induced measure on
\[
\mathfrak P_3=\{\varrho\in M_3(\mathbb C)\mid \varrho=\varrho^\dagger,\ \varrho\ge 0,\ \operatorname{tr}\varrho=1\},
\]
with invariance under the adjoint action of \(SU(3)\). The Hilbert–Schmidt metric is
\[
g_{HS}\propto \operatorname{Tr}(d\varrho\otimes d\varrho),
\qquad
ds_{HS}^2=\operatorname{Tr}(d\varrho^2).
\]
For the regular stratum of nondegenerate spectra,
\[
\varrho=U\,\operatorname{diag}(r_1,r_2,r_3)\,U^\dagger,
\qquad
U\in U(3)/U(1)^3,
\]
and the Hilbert–Schmidt eigenvalue joint density on the ordered simplex is
\[
P^{HS}(r_1,r_2,r_3)\propto
\delta\!\left(1-\sum_{j=1}^3 r_j\right)\prod_{j<k}(r_j-r_k)^2.
\]
For degenerate strata of partition type \((k_1,\dots,k_s)\), the density becomes
\[
P^{HS}_{k_1,\dots,k_s}(r_1,\dots,r_s)
\propto
\delta\!\left(1-\sum_{i=1}^s k_i r_i\right)
\prod_{i<j}(r_i-r_j)^{2k_i k_j},
\]
with Haar-distributed angular variables on the corresponding coset \(U(3)/(U(k_1)\times\cdots\times U(k_s))\) [2208.13908].

The qutrit state space is stratified as
\[
\mathfrak P_3=
\mathfrak P_{[T^3]}
\cup
\mathfrak P_{[S(U(2)\times U(1))]}
\cup
\mathfrak P_{[SU(3)]}.
\]
The regular stratum \(\mathfrak P_{[T^3]}\) has nondegenerate spectra and orbit space equal to the interior of an ordered simplex face. The degenerate stratum \(\mathfrak P_{[S(U(2)\times U(1))]}\) corresponds to spectral types \((2,1)\) and \((1,2)\). The maximally mixed state is the single-point stratum \(\mathfrak P_{[SU(3)]}\) [2208.13908].

Within this geometry, “classical states” are those with everywhere positive qutrit Wigner function:
\[
\mathfrak P_3^{(+)}
=
\{\varrho\in\mathfrak P_3\mid W_\varrho(z)\ge 0,\ \forall z\in\Omega_3\}.
\]
The associated geometric-probability indicator is
\[
\mathcal Q_3
=
\frac{\mathrm{Volume}(\text{Classical States})}
{\mathrm{Volume}(\text{All States})},
\]
with a stratum-wise version \(\mathcal Q_3[H_\alpha]\). Because of \(SU(3)\) invariance, the ratios reduce to eigenvalue integrals over orbit-space regions determined by Wigner-function positivity. For the Hilbert–Schmidt ensemble, the paper derives closed forms on both the regular and degenerate strata. On the regular stratum,
\[
\mathcal Q^{HS}_{[T^3]}(\zeta)
=
\frac{20\cos^2\!\big(\zeta-\frac{\pi}{6}\big)+1}
{128\left(4\cos^2\!\big(\zeta-\frac{\pi}{6}\big)-1\right)^5},
\]
with representative values
\[
\mathcal Q^{HS}_{[T^3]}(0)
=
\mathcal Q^{HS}_{[T^3]}(\pi/3)
=
\frac{1}{256},
\qquad
\mathcal Q^{HS}_{[T^3]}(\pi/6)
=
\frac{21}{31104}.
\]
On the degenerate stratum,
\[
\mathcal Q^{HS}_{[S(U(2)\times U(1))]}(\zeta)
=
\frac{1}{1056}
\left[
\csc^5\!\Big(\zeta+\frac{\pi}{6}\Big)+\sec^5(\zeta)
\right],
\]
with
\[
\mathcal Q^{HS}_{[S(U(2)\times U(1))]}(0)
=
\mathcal Q^{HS}_{[S(U(2)\times U(1))]}(\pi/3)
=
\frac{1}{32},
\]
and
\[
\mathcal Q^{HS}_{[S(U(2)\times U(1))]}(\pi/6)\approx 0.003887.
\]
The global Hilbert–Schmidt indicator has minimum
\[
\min \mathcal Q_3(\zeta)=0.0006751
\quad\text{at}\quad
\zeta=\pi/6,
\]
and satisfies
\[
\mathcal Q_3(0)-\mathcal Q_3(\pi/3)=0.
\]
Thus the Hilbert–Schmidt qutrit indicator is exactly symmetric in the moduli parameter \(\zeta\), while the corresponding Bures and Bogoliubov–Kubo–Mori indicators display mild symmetry breaking [2208.13908].

The qutrit comparison also shows a consistent ordering across regular and degenerate strata:
\[
\text{HS} > \text{Bures} > \text{BKM}
\]
for the probability of Wigner-function positivity. Degenerate strata are more classical than the regular stratum for all three ensembles. This suggests that, within the geometric-probability framework used there, the Hilbert–Schmidt ensemble places comparatively more weight on regions of state space compatible with positive Wigner functions, while higher-symmetry strata amplify that tendency [2208.13908].

Source: https://www.emergentmind.com/topics/hilbert-schmidt-ensemble