The emergence of randomness from unitary quantum dynamics is a central problem across diverse disciplines, ranging from the foundations of statistical mechanics to quantum algorithms and quantum computation. Physical systems are invariably subject to constraints -- from simple scalar symmetries to more complex non-Abelian ones -- that restrict the accessible regions of Hilbert space and obstruct the generation of pure random states. In this work, we show that for systems with noncommuting symmetries such as SU(2), the degree of Haar-like randomization achievable under unitary dynamics is strongly constrained by experimental limitations on state initialization, in particular low-entanglement initial states, rather than by the symmetry-constrained dynamics themselves. Specifically, we show that time-evolved states can, in principle, reproduce Haar-like behavior at the level of finite statistical moments (i.e., those accessible under realistic experimental conditions with a finite number of state copies) provided that the initial state matches the corresponding moments of the conserved operators in the Haar ensemble. However, for the unentangled initial states commonly used in programmable quantum systems, this condition cannot be satisfied. Consequently, even at asymptotically long times in strongly quantum-chaotic regimes, late-time states remain distinguishable from Haar-random states in probes such as entanglement entropy, with deviations from Haar behavior that remain finite with increasing system size. We quantify the maximal entanglement entropy achievable and identify the unentangled initial conditions that yield the most entropic late-time states. Our results show that the combination of non-Abelian symmetry structure and experimental constraints on state preparation can strongly limit the degree of Haar-like randomization achievable at late times.
The emergence of randomness from unitary quantum dynamics is a central problem across diverse disciplines, ranging from the foundations of statistical mechanics to quantum algorithms and quantum computation. Physical systems are invariably subject to constraints -- from simple scalar symmetries to more complex non-Abelian ones -- that restrict the accessible regions of Hilbert space and obstruct the generation of pure random states. In this work, we show that for systems with noncommuting symmetries such as SU(2), the degree of Haar-like randomization achievable under unitary dynamics is strongly constrained by experimental limitations on state initialization, in particular low-entanglement initial states, rather than by the symmetry-constrained dynamics themselves. Specifically, we show that time-evolved states can, in principle, reproduce Haar-like behavior at the level of finite statistical moments (i.e., those accessible under realistic experimental conditions with a finite number of state copies) provided that the initial state matches the corresponding moments of the conserved operators in the Haar ensemble. However, for the unentangled initial states commonly used in programmable quantum systems, this condition cannot be satisfied. Consequently, even at asymptotically long times in strongly quantum-chaotic regimes, late-time states remain distinguishable from Haar-random states in probes such as entanglement entropy, with deviations from Haar behavior that remain finite with increasing system size. We quantify the maximal entanglement entropy achievable and identify the unentangled initial conditions that yield the most entropic late-time states. Our results show that the combination of non-Abelian symmetry structure and experimental constraints on state preparation can strongly limit the degree of Haar-like randomization achievable at late times.
Randomness generation through quantum-chaotic evolution underpins foundational questions in statistical mechanics and applications across quantum information science, including benchmarking, tomography, metrology, and demonstrations of quantum computational advantage. While statistical mechanics successfully captures the temporal averages of local observables, understanding randomness at the level of higher statistical moments remains a daunting challenge, with analytic progress largely confined to random quantum circuit models or fine-tuned systems exhibiting space-time duality. Here we study how much randomness can be dynamically generated by generic quantum-chaotic evolution under physical, non-random Hamiltonians. Combining theoretical insights with numerical simulations, we show that for broad classes of initially unentangled states, the dynamics become effectively Haar-random well before the system can ergodically explore the physically accessible Hilbert space. Both local and highly nonlocal observables, including entanglement measures, equilibrate to their Haar expectation values and fluctuations on polynomial timescales with remarkably high numerical precision, and with the fastest randomization occurring in regions of parameter space previously identified as maximally chaotic. Interestingly, this effective randomization can occur on timescales linear in system size, suggesting that the sub-ballistic growth of Renyi entropies typically observed in systems with conservation laws can be bypassed in non-random Hamiltonians with an appropriate choice of initial conditions.
Initially unentangled states undergoing quantum chaotic dynamics are expected to evolve into featureless states at late times. While this expectation holds true on an averaged, coarse-grained level, it is unclear if this expectation applies to higher statistical moments, as symmetries present in physical systems constrain the exploration of Hilbert space. Here, we study the universal structure of late-time ensembles obtained from unitary dynamics in quantum chaotic systems with symmetries, such as charge or energy conservation. We find that while quantum states do not ergodically explore the entire Hilbert space at late times, the late-time ensemble typically becomes indistinguishable from Haar random states in the thermodynamic limit and at the level of finite statistical moments. Importantly, our results apply to initial states that are easy to prepare in ongoing experiments---specifically, product states---that lie in the middle of the spectrum of quantum chaotic systems. We show that typically these states not only exhibit the same late-time ensemble average as Haar random states, but also the same state-to-state fluctuations and higher statistical moments. In other words, there is no measurement---whether local or nonlocal---that one can perform at the level of finite statistical moments to tell that states are not exploring the entire Hilbert space. Interestingly, within the class of low-entanglement initial states, we also find atypical initial conditions in the middle of the spectrum of Hamiltonians known to be ``maximally chaotic''. Such atypical initial states have smaller variance of the symmetry operator relative to that of Haar random states, and evolve into nonuniversal ensembles that can be distinguished from the Haar ensemble by simple measurements or subsystem properties. In the limiting case of initial states with negligible variance of the symmetry operator (e.g., states with fixed particle number or energy eigenstates), the late-time ensemble has universal behavior captured by constrained random state ensembles. On the one hand, our results reveal that an extremely high level of quantum state randomness can typically be achieved even when dynamics are constrained by symmetries; on the other hand, our results show that one can still achieve strongly nonergodic dynamics at ``infinite temperature'' within quantum chaotic regimes.
Quantum many-body systems display an extraordinary degree of complexity, yet many of their features are universal: they depend not on microscopic details, but on a few fundamental physical aspects such as symmetries. A central challenge is to distill these universal characteristics from model-specific ones. Random quantum states sampled from a uniform distribution, the Haar measure, provide a powerful framework for capturing this typicality. Here, we experimentally study the entanglement and symmetries of random many-body quantum states generated by evolving simple product states under ergodic Floquet models. We find excellent agreement with the predictions from the Haar-random state ensemble. First, we measure the Rényi-2 entanglement entropy as a function of the subsystem size, observing the Page curve. Second, we probe the subsystem symmetries using entanglement asymmetry. Finally, we measure the moments of partially transposed reduced density matrices obtained by tracing out part of the system in the generated ensembles, thereby revealing distinct entanglement phases. Our results offer an experimental perspective on the typical entanglement and symmetries of many-body quantum systems.
Rydberg atom arrays are powerful platforms for studying quantum many-body systems. We consider the Rydberg-Ising Hamiltonian in a ring geometry and numerically study ensembles of states generated by random global pulse sequences subject to hardware constraints and fixed evolution times. We compare the statistical properties of such states with those of Haar-random states within the relevant lattice symmetry sector. In the strong-interaction regime (short interatomic distance), the dynamics is governed by an effective blockade that restricts Hilbert-space exploration and limits entanglement growth. In this regime, level-spacing statistics of the entanglement spectrum are close to random-matrix predictions, while the distribution of measurement probabilities deviates from Porter-Thomas behavior. For weaker interactions (larger interatomic distance), the system approaches Haar-like statistics at long times, as reflected in entanglement entropy, entanglement-spectrum statistics, and the distribution of measurement probabilities. At intermediate interactions, this behavior is observed on experimentally relevant timescales. Motivated by this observation, we investigate whether individual target states can be efficiently prepared using quantum optimal control in this regime. Employing target states drawn from an ensemble with a broad entropy distribution, we observe high fidelities (infidelities between 10^-5 and 3 x 10^-2 for nine spins). Within this ensemble, the fidelity decreases with the entanglement entropy of the target state, indicating that highly entangled targets are more difficult to prepare under the considered finite-time control constraints.
We study the average bipartite entanglement entropy of Haar-random pure states in quantum many-body systems with global SU(2) symmetry, constrained to fixed total spin J and magnetization Jz=0. Focusing on spin-21 lattices and subsystem fractions f<21, we derive a asymptotic expression for the average entanglement entropy up to constant order in the system volume V. In addition to the expected leading volume law term, we prove the existence of a 21logV finite-size correction resulting from the scaling of the Clebsch-Gordon coefficients and compute explicitly the O(1) contribution reflecting angular-momentum coupling within magnetization blocks. Our analysis uses features of random matrix ensembles and provides a fully analytical treatment for arbitrary spin densities, thereby extending Page type results to non-Abelian sectors and clarifying how SU(2) symmetry shapes average entanglement.
Abstract Scrambling is a process by which the state of a quantum system is effectively randomized due to the global entanglement that “hides” initially localized quantum information. Closely related notions include quantum chaos and thermalization. Such phenomena play key roles in the study of quantum gravity, many-body physics, quantum statistical mechanics, quantum information etc. Scrambling can exhibit different complexities depending on the degree of randomness it produces. For example, notice that the complete randomization implies scrambling, but the converse does not hold; in fact, there is a significant complexity gap between them. In this work, we lay the mathematical foundations of studying randomness complexities beyond scrambling by entanglement properties. We do so by analyzing the generalized (in particular Rényi) entanglement entropies of designs, i.e. ensembles of unitary channels or pure states that mimic the uniformly random distribution (given by the Haar measure) up to certain moments. A main collective conclusion is that the Rényi entanglement entropies averaged over designs of the same order are almost maximal. This links the orders of entropy and design, and therefore suggests Rényi entanglement entropies as diagnostics of the randomness complexity of corresponding designs. Such complexities form a hierarchy between information scrambling and Haar randomness. As a strong separation result, we prove the existence of (state) 2-designs such that the Rényi entanglement entropies of higher orders can be bounded away from the maximum. However, we also show that the min entanglement entropy is maximized by designs of order only logarithmic in the dimension of the system. In other words, logarithmic-designs already achieve the complexity of Haar in terms of entanglement, which we also call max-scrambling. This result leads to a generalization of the fast scrambling conjecture, that max-scrambling can be achieved by physical dynamics in time roughly linear in the number of degrees of freedom. This paper is an extended version of Phys. Rev. Lett. 120 (2018) 130502 [1].