Mixed-State Deep Thermal Ensemble

• Mixed-State Deep Thermal Ensemble is a framework describing probability distributions over quantum mixed states conditioned by measurement or thermodynamic constraints, extending traditional pure-state thermalization.
• It employs techniques such as maximum-entropy purification, augmented conditioning, and moment analysis to reveal universal behavior in chaotic quantum dynamics.
• The ensemble framework is operationally significant for quantum information tasks, including teleportation protocols and variational Gibbs state preparation on quantum devices.

Searching arXiv for the cited works to ground the article and verify bibliographic details. arXiv search query: (Yu et al., 12 May 2025) Mixed state deep thermalization arXiv search query: (Sherry et al., 18 Jul 2025) Do mixed states exhibit deep thermalisation? arXiv search query: (Ippoliti et al., 2022) Solvable model of deep thermalization with distinct design times arXiv search query: (Nakata et al., 2014) Thermal states of random quantum many-body systems arXiv search query: (Miller, 1 Aug 2025) Statistical Mechanics of Random Mixed State Ensembles with Fixed Energy Mixed-state deep thermal ensemble denotes a family of constructions in which the object of interest is not a single reduced density matrix but a probability distribution over mixed states, typically attached to a local subsystem and generated by conditioning, tracing, or imposing thermodynamic constraints. In the narrow sense used in recent work on deep thermalisation, it is the universal late-time ensemble of mixed subsystem states produced either by incomplete measurements on a complementary region or by unitary evolution from an initially mixed state, after reformulating the maximum-entropy principle on an augmented purification (Yu et al., 12 May 2025, Sherry et al., 18 Jul 2025). In a broader statistical-mechanical sense, closely related ensembles arise from random Hamiltonians at finite temperature, from fixed-energy measures over density matrices, and from variational multiscale ansätze for Gibbs-state preparation (Nakata et al., 2014, Miller, 1 Aug 2025, Sewell et al., 2022). This suggests that the topic sits at the intersection of quantum thermalisation, random-state theory, and nonequilibrium quantum information.

1. Conceptual framework and basic definitions

Deep thermalisation is stronger than ordinary thermalisation. In the projected-ensemble formulation, a bipartite pure state $\ket{\Psi}_{AB}$ defines an ensemble of conditional pure states on $A$ by

$$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$

with

$$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$

The $k$-th moment operator is

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$

and exact Haar-like behaviour means $\rho^{(k)}$ matches the Haar moment operator $\rho_H^{(k)}$ (Ippoliti et al., 2022).

For mixed states, the conditioning procedure no longer produces pure states in general. One formulation defines the projected ensemble from a mixed state $\varrho$ as

$${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$

with

$A$0

A related construction, the mixed state projected ensemble (MSPE), appears when measurements on the complement are incomplete, so that some outcomes are lost and the conditional state on $A$1 is mixed rather than pure (Sherry et al., 18 Jul 2025, Yu et al., 12 May 2025).

The moment language remains central in the mixed-state setting. In the MSPE approach, the $A$2-th moment is

$A$3

while in the mixed-state deep-thermalisation approach the natural reference ensemble is

$A$4

whose moments encode the universal mixed-state distribution associated with the initial state $A$5 (Yu et al., 12 May 2025, Sherry et al., 18 Jul 2025).

2. Failure of the pure-state paradigm and the augmented-system construction

A central 2025 result is that the standard pure-state deep-thermalisation framework fails sharply for mixed initial states, even with infinitesimal initial mixedness. For the Haar-scrambled reference ensemble $A$6, the second moment is

$A$7

which implies the average purity

$A$8

Equality holds only if $A$9 is pure. For any finite mixedness, however small, $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$0, so the ensemble is not an ensemble of pure states. It therefore cannot be the Haar ensemble of pure states, and it cannot be the usual Scrooge ensemble of pure states either (Sherry et al., 18 Jul 2025).

The replacement paradigm introduces a mixed-state deep thermal ensemble by purification. One purifies the initial mixed state as

$$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$1

constructs a maximum-entropy projected ensemble of pure states on the augmented system $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$2, and then traces out $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$3. The resulting moment formula is

$$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$4

This is exactly the form of a Scrooge ensemble on $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$5, followed by partial trace over $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$6 (Sherry et al., 18 Jul 2025).

The dependence on the initial state is explicit and spectral. The general $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$7-moment is

$$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$8

For $$\mathcal{E} = \{ (p(z), \ket{\psi_z}_A):\ z = 1, \dots d_B \},$$9, the dependence reduces to $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$0, so the distance from the Haar ensemble increases monotonically with the second Rényi entropy

$$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$1

for $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$2. In the special case where all nonzero eigenvalues of $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$3 are equal, the ensemble reduces to a generalised Hilbert-Schmidt ensemble,

$$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$4

This establishes a maximum-entropy principle for mixed states that is fundamentally distinct from the pure-state case (Sherry et al., 18 Jul 2025).

3. Mixed state projected ensembles and incomplete measurement

A second major route to mixed-state deep thermal ensembles begins from imperfect measurements. In the MSPE construction, a many-body system is partitioned into subsystem $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$5 and complementary subsystem $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$6, with $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$7 further split into contiguous regions $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$8. Outcomes on $$p(z) = \| \, {}_B\!\braket{z}{\Psi}_{AB} \|^2, \qquad \ket{\psi_z}_A = {}_B\!\braket{z}{\Psi}_{AB} / \sqrt{p(z)}.$$9 and $k$0 are retained, while the outcomes on $k$1 sites in $k$2 are lost. If $k$3 denotes the retained outcome string, then

$k$4

with

$k$5

The MSPE is

$k$6

If no measurement outcomes are lost, $k$7, then $k$8 is pure and the construction reduces to the usual pure-state projected ensemble (Yu et al., 12 May 2025).

The solvable setting studied for this problem is a $k$9D dual-unitary brick-wall circuit with an initial state given by a tensor product of Bell pairs,

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$0

and local measurements in the Bell-like basis

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$1

Lost outcomes act as a totally depolarizing channel or partial trace, which is the mechanism that converts conditional pure states into mixed states (Yu et al., 12 May 2025).

The late-time universality statement is explicit: for any fixed ratio $$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$2, the MSPE deep thermalizes to a generalized Hilbert-Schmidt ensemble,

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$3

When $$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$4, this becomes the usual Hilbert-Schmidt ensemble of random density matrices. The derivation uses the permutation-operator basis

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$5

with expansion

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$6

and asymptotic coefficients

$$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$7

where $$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$8 is the number of cycles of the permutation $$\rho^{(k)} = \sum_z p(z)\, (\ket{\psi_z}\bra{\psi_z})^{\otimes k},$$9 (Yu et al., 12 May 2025).

The convergence is exponentially fast in the circuit depth. For large $\rho^{(k)}$0 and $\rho^{(k)}$1,

$\rho^{(k)}$2

and

$\rho^{(k)}$3

Accordingly, the thermalization time for fixed accuracy $\rho^{(k)}$4 scales as

$\rho^{(k)}$5

This establishes a mixed-state version of design-time separation in a fully solvable setting (Yu et al., 12 May 2025).

4. Design times, GAP/Scrooge ensembles, and solvable chaotic dynamics

An exactly solvable precursor to the mixed-state theory studies deep thermalization in a model consisting of a finite subsystem $\rho^{(k)}$6, a small bottleneck subsystem $\rho^{(k)}$7, and a large bath $\rho^{(k)}$8, with dynamics generated by alternating Haar-random gates on $\rho^{(k)}$9 and $\rho_H^{(k)}$0. After projective measurement of the bath, the projected ensemble on $\rho_H^{(k)}$1 approaches a universal limit. In the thermodynamic limit $\rho_H^{(k)}$2, the limiting ensemble on $\rho_H^{(k)}$3 is the Scrooge or Gaussian Adjusted Projected (GAP) ensemble, with induced density matrix

$\rho_H^{(k)}$4

and projected-ensemble moment

$\rho_H^{(k)}$5

This is the maximally entropic ensemble compatible with $\rho_H^{(k)}$6 (Ippoliti et al., 2022).

The same work makes the distinction between regular and deep thermalization sharp. Regular thermalization concerns only

$\rho_H^{(k)}$7

while deep thermalization requires convergence of all projected-ensemble moments toward Haar moments. The deviation is measured by

$\rho_H^{(k)}$8

and the design time $\rho_H^{(k)}$9 is the earliest time when $\varrho$0. Since $\varrho$1, one has

$\varrho$2

Near infinite temperature,

$\varrho$3

and the exact result for the design times is

$\varrho$4

For $\varrho$5,

$\varrho$6

while

$\varrho$7

These formulas show that deep thermalization occurs later than ordinary thermalization and that higher moments equilibrate more slowly (Ippoliti et al., 2022).

This framework is not itself a mixed-state theory, because its conditional states are pure. However, it supplies the structural template later generalized to mixed-state deep thermalization: projected ensembles, maximum-entropy limiting distributions, and a hierarchy of moment-dependent equilibration times. The mixed-state theories of 2025 can therefore be read as extensions of the GAP/Scrooge paradigm to settings with incomplete measurement or finite initial entropy (Ippoliti et al., 2022, Sherry et al., 18 Jul 2025).

5. Thermal-state ensembles, random density matrices, and fixed-energy constraints

A broader line of work studies thermal or constrained ensembles of density matrices without deriving them from measurement-conditioned subsystem states. For a random Hamiltonian $\varrho$8, the thermal state is

$\varrho$9

and the ensemble of thermal states is $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$0. Randomness is quantified by the mixed-state $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$1-design distance

$${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$2

with

$${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$3

For random global Hamiltonians $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$4, the distance decreases monotonically as $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$5 increases,

$${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$6

interpolating between the maximally mixed state at $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$7 and the Haar-random ground-state ensemble as $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$8. The ensemble becomes an $${\cal E}_{\rm PE}[\varrho] := \{p(z_B), \varrho_A(z_B)\}_{z_B},$$9-approximate $A$00-design at

$A$01

equivalently at temperature

$A$02

For random $A$03-local Hamiltonians, by contrast,

$A$04

so the ensemble is always an exact state $A$05-design, but it does not become a $A$06-design for $A$07. Numerical evidence indicates two temperature regimes separated by a singular point, with values reported roughly around

$A$08

for different interaction ranges (Nakata et al., 2014).

A different generalization treats density matrices themselves as microstates. In that approach, the space of states is $A$09, equipped with a measure $A$10, and one imposes the fixed-average-energy constraint

$A$11

with density of states

$A$12

The ensemble average of an observable is

$A$13

and the average density matrix

$A$14

obeys

$A$15

The associated thermodynamic quantities are

$A$16

and

$A$17

A useful response identity is

$A$18

for the perturbation $A$19 (Miller, 1 Aug 2025).

For a qubit with

$A$20

the average density matrix is

$A$21

which can also be written in canonical form only in this special two-level case. For $A$22 non-interacting spins in a magnetic field,

$A$23

the Bures-Hall ensemble exhibits a phase-transition-like phenomenon without interactions, and the relative energy fluctuations do not vanish generically in the thermodynamic limit. This suggests that mixed-state ensemble geometry alone can generate nonstandard thermodynamic behaviour (Miller, 1 Aug 2025).

6. Operational significance: teleportation, imperfect measurements, and multiscale Gibbs ansätze

Mixed-state deep thermal ensembles are operationally relevant because they arise naturally when measurements are lossy or when the initial state has finite entropy. In the MSPE setting, the states composing the ensemble can mediate teleportation of quantum information from the right end of the system to the left subsystem $A$24. The relevant diagnostic is the annealed Rényi-$A$25 conditional entropy

$A$26

Its values distinguish three regimes: $A$27

$A$28

and

$A$29

The sharp transition occurs at

$A$30

For $A$31, one has the teleportation phase; for $A$32, the decoupled phase. In the large-$A$33 limit, the exact $A$34 result is

$A$35

This connects mixed-state deep thermalization directly to conditional entropy and information transfer (Yu et al., 12 May 2025).

A distinct operational manifestation is variational Gibbs-state preparation. The thermal multi-scale entanglement renormalization ansatz (TMERA) modifies DMERA by injecting each new qubit not in the pure state $A$36 but in a single-qubit mixed state

$A$37

and evolving via

$A$38

Because the circuit is unitary, the final state is a product-spectrum thermal ansatz. The mode energies are taken as

$A$39

and the single-qubit entropy is

$A$40

The free-energy objective is

$A$41

In the paper’s interpretation, this realizes a mixed-state deep thermal ensemble over wavepacket modes localized to different scales (Sewell et al., 2022).

The benchmark system is the 1D critical transverse-field Ising model

$A$42

equivalently the free-Majorana Hamiltonian

$A$43

For $A$44, TMERA produces global fidelities $A$45 for $A$46-site systems across all temperatures; even $A$47 gives fidelity above $A$48 across all temperatures. The per-site infidelity behaves approximately as

$A$49

At low temperatures, the energy above the ground state scales roughly as $A$50, and the entropy per qubit scales roughly as $A$51. The largest reported entropy mismatch for $A$52 is about $A$53 bits per qubit near $A$54, and the largest energy mismatch is about $A$55 per qubit near $A$56 (Sewell et al., 2022).

Taken together, these results show that mixed-state deep thermal ensembles are not merely abstract limiting distributions. They encode experimentally relevant imperfect-measurement statistics, furnish solvable universality classes for chaotic dynamics, connect to random density-matrix measures such as the generalized Hilbert-Schmidt, Hilbert-Schmidt, Bures-Hall, and GAP/Scrooge ensembles, and provide constructive ansätze for thermal state preparation on quantum devices (Yu et al., 12 May 2025, Sherry et al., 18 Jul 2025, Sewell et al., 2022).