Noise is unavoidable in realistic quantum systems, and resulting states are generally mixed. A key question regarding deep thermalization in these mixed states, whether a universal limiting distribution exists, has remained unanswered. This study answers that question for quantum dynamics containing local quantum channels representing errors. The late-time projected ensemble of mixed states is completely determined by the normalized Choi state of the error channel; it exhibits universality independent of microscopic details within the unitary dynamics. These findings hold analytically for dual-unitary circuits, global Haar-random dynamics, and local Haar-random circuits at large local dimensions.

Error-Driven Equilibration Reveals Final States Depend Only on Error Characteristics

The team at California Institute of Technology and University of Ljubljana have demonstrated exponentially fast convergence towards limiting ensembles, a previously unattainable feat. They observed mixed-state projected ensembles now converge within timescales inversely proportional to error strength, crossing a threshold for practical applications where equilibration was once too slow to measure. This research establishes that the final distribution of states is entirely determined by characteristics inherent in local quantum errors during complex system evolution: specifically, the normalized Choi state, regardless of initial conditions or underlying dynamics.

Analytical calculations applied to three distinct types of systems established this universality: dual-unitary circuits allowing precise analysis, global Haar-random dynamics representing fully chaotic behaviour, and large dimensional local Haar-random circuits increasing realism. Numerical tests using qubit circuits, the basic units of quantum information, and kicked Ising dynamics further confirmed these findings across different physical scenarios. Error strength dictates how quickly equilibrium is reached; exponentially fast convergence towards this limiting ensemble achieves timescales inversely proportional to its magnitude.

Conditional States Reveal Error Channel Effects During Chaotic Quantum Dynamics

By considering states conditioned on measurements of the surroundings, we address which properties of a local quantum error remain observable after chaotic evolution spreads it through a many-body system. Measuring the complement B of retained region A yields conditional states ρNA(α) with Born probabilities Pα.

Their average corresponds to the reduced density matrix used in conventional thermalization, while higher moments resolve fluctuations among outcomes and reveal data about the error channel determining this distribution when noise acts during evolution and every final measurement outcome on B is retained; modern quantum simulators can read out most of a many-body system per experimental run and record all measurement outcomes.

For generic chaotic dynamics, the projected ensemble approaches a limiting distribution dependent only on coarse-grained properties like symmetries or conserved charges. This form of equilibration is known as deep thermalization. Under suitable conditions, Haar statistics emerge while conservation laws, symmetries, and kinetic constraints select charge-resolved ensembles with Gaussian analogues for free fermions and continuous-variable bosons. Constraints on resources such as coherence can even prevent deep thermalization despite an otherwise thermal reduced state.

Emergent randomness has aided benchmarking, randomstate generation, and noise learning, all relying on knowing which ensemble emerges from the dynamics; these have been measured in Rydberg and superconducting simulators. We previously introduced the mixed-state projected ensemble (MSPE) by considering incomplete measurements on subsystem B, equivalently interpreted as completely depolarizing channels acting during or after quantum dynamics. This MSPE approaches a generalised Hilbert, Schmidt ensemble characterised by a single parameter encoding error number.

Initially mixed states evolving unitarily provide another setting whose limiting ensemble retains initial density matrix spectrum. Partial bath measurements generate mixed conditional states through incomplete observation. These results leave open questions about generic quantum errors: partial depolarization, dephasing, and amplitude damping retain different information about input states and cannot be described solely by one parameter like output purity.

In weak noise limits, bit-string statistics from noisy random circuits become nearly independent of the specific error channel; therefore determining which random-density-matrix ensemble emerges after scrambling, or how much microscopic information remains at late times, is unclear. Here we study mixed-state deep thermalization inserting dilute local channels into bulk dynamics while retaining all outcomes on B. Unlike protocols based on incompleteness, this generates mixness entirely through physical error channels.

We find that the MSPE’s defining information is encoded only in the normalized Choi state of the channel after sufficient scrambling because errors are local, this information stays finite and serves as macroscopic data characterising the late-time ensemble.

The resulting limiting distribution is a Bornweighted normalized correlated Wishart ensemble, which we call the channel-induced ensemble (CIE). It’s universal because for fixed error channel it remains independent of microscopic unitary details or initial pure states; different channels select differing ensembles with their respective Choi states serving as relevant retained data post thermalization.

A unitary error recovers Haar while complete depolarization reproduces our earlier generalised Hilbert, Schmidt result, partial depolarization interpolates between these limits; general CPTP maps select other spectra. Approaching the limit contains additional information: in solvable settings, equal spectra yield same limit but transient corrections depend on input/output structure of the Choi state.

Under replica continuation assumptions, exponential convergence is obtained, and its moment-dependent prefactors distinguish relaxation rate from accuracy time needed. We analytically derive CIE convergence within dual-unitary circuits, extended to noisy versions with space/time unitality constraints at cost of added complexity. Numerical simulations using qubit circuits and kicked Ising dynamics further support this conjecture.

Section II defines the ensemble and main result; Sections III & V present analytical derivations while section VI details numerical tests before discussing extensions in VII. We assume complete measurement on B retaining every outcome. The system is initialised to pure state undergoing unitary evolution interspersed with dilute local error channels {N1, N2,.}. As a result, the final state becomes mixed denoted ρ(t).

We say MSPE deeply thermalizes to a reference if, for every fixed k, ∆(k) ξ vanishes at late times; decay of ∆(k) ξ indicates approach rate.

We denote corresponding probability measure by dνM(X) and resulting ensemble EM={PM(X), ρM(X)}. This Bornweighted normalized correlated Wishart ensemble, Proposition 1 states that dilute local quantum errors do not destroy deep thermalization; instead they select the limiting CIE determined solely by error channel’s normalized Choi state independent of microscopic unitary details. The physical picture remains simple: Unitary evolution produces a single many-body state, an ensemble arises after measuring subsystem B, sufficiently scrambling dynamics sample Haar ensembles as in pure-state deep thermalization.

Local Error Analysis Constrains Broader System Implications

While these findings definitively link error characteristics to predictable outcomes within noisy systems, benefiting practical device construction, limitations remain due to focus on *local* errors affecting only individual qubits or small groups thereof. This raises questions about how results change when noise correlates across larger portions of the system; previous research suggests long-range interactions can dramatically alter scrambling dynamics potentially invalidating universality assumptions demanding new theoretical frameworks. However acknowledging that these findings rest upon *local* errors does not diminish their immediate value, understanding even limited imperfections scrambles information effectively.

This research demonstrated that quantum dynamics with local errors converge to a universal limiting distribution determined solely by the normalized Choi state of the error channel, meaning the final mixed-state ensemble is independent of the initial unitary details. The study established this result analytically for several types of circuits and verified it numerically in qubit systems and kicked Ising dynamics, finding exponentially fast convergence in some cases.

These results extend the concept of deep thermalization to include noisy quantum systems, revealing a direct connection between quantum errors and predictable random ensembles. This work clarifies how noise impacts scrambling within these systems using characteristics defined by dimension d2n where n equals number of channels.

👉 More information
🗞 Mixed State Deep Thermalization under Generic Quantum Errors
✍️ Xie-Hang Yu (California Institute of Technology); Pavel Kos (University of Ljubljana)
🧠 ArXiv: https://arxiv.org/abs/2610.31580

More like this