Modern quantum experiments can probe many-body wavefunctions through projective measurements, providing snapshots of individual many-body configurations. A fundamental question is whether the intrinsic structure of their measurement distributions can reveal quantum dynamics beyond predefined observables. Here, we address this question at the many-body ergodicity-breaking transition using two complementary characteristics of nonequilibrium wavefunction snapshots: their binary intrinsic dimension (BID) and the topology of wavefunction networks constructed from Hamming distances. We show that the critical point is characterized by an extensive but submaximal BID, and scale-free network connectivity with highly connected hubs. These signatures distinguish the ergodic, critical, and nonergodic regimes, and provide experimentally accessible probes of ergodicity breaking from projective measurement snapshots.
Modern quantum experiments can probe many-body wavefunctions through projective measurements, providing snapshots of individual many-body configurations. A fundamental question is whether the intrinsic structure of their measurement distributions can reveal quantum dynamics beyond predefined observables. Here, we address this question at the many-body ergodicity-breaking transition using two complementary characteristics of nonequilibrium wavefunction snapshots: their binary intrinsic dimension (BID) and the topology of wavefunction networks constructed from Hamming distances. We show that the critical point is characterized by an extensive but submaximal BID, and scale-free network connectivity with highly connected hubs. These signatures distinguish the ergodic, critical, and nonergodic regimes, and provide experimentally accessible probes of ergodicity breaking from projective measurement snapshots.
The local Hilbert space dimension is a fundamental parameter in quantum many-body systems, but its role in controlling the ergodicity-breaking transition remains less explored. We systematically study this question by introducing toy models of qudits in zero dimensions that exhibit an ergodicity-breaking transition. Using the hybridization condition, we analytically obtain the dependence of the critical value of the interactions on the local Hilbert space dimension. Numerical studies of spectral correlation measures show agreement with the theoretical predictions, hence confirming the existence of the ergodicity-breaking transition in these models. We further introduce a distinct higher spin version of the model and numerically obtain the same critical behavior, demonstrating that the observed critical point is independent of the microscopic structure of interactions considered.
The fate of ergodicity and localization in disordered, interacting quantum systems in two spatial dimensions remains an outstanding open problem in non-equilibrium physics. Theoretical and numerical approaches are severely constrained by exponential Hilbert space growth. Here, we leverage the high data-acquisition rates of superconducting quantum processors to extensively sample Hilbert space configurations, effectively imaging the state as a wavefunction network (WFN). Tracking the bitstring return probability across varied disorder strengths, we find heavy-tailed distributions and power-law typical decay characteristic of glassy relaxation. Furthermore, the distribution of bitstring probabilities transitions from a Porter-Thomas form at low disorder to an eight-decade power law at stronger disorder. In the reconstructed configuration space, WFN connectivity is compatible with a random geometric network at low disorder and becomes scale-free at stronger disorder. Finally, a multifractal analysis based on the binary intrinsic dimension (BID) reveals a non-ergodic regime where the wavefunction support saturates to a finite fraction of the total Hilbert space. Together, these results rule out a direct ergodic-to-localized transition, pointing instead to an extended multifractal regime. By visualizing wavefunction dynamics directly in configuration space, this work complements traditional real-space studies, showing the potential of quantum processors to resolve long-standing questions in non-equilibrium physics.
Quantum many-body systems can exhibit distinct regimes where dynamics is either ergodic, dynamically exploring an extensive region of available state-space, or non-ergodic, where the dynamics may be restricted. An example is the many-body localization (MBL) transition, where disorder induces non-ergodic behaviour. Most measures of ergodicity notably rely on global quantities, such as level spacing statistics. We explore the ability for a subsystem to probe the ergodicity of dynamics via measurement of local observables, and use expected results from random matrix theory (RMT) as a benchmark for the ergodic regime. We exploit two predictions from RMT as ergodicity is broken: the time evolution of the quantum Fisher information, and a fluctuation-dissipation relation. These are investigated in three different ergodicity breaking mechanisms, namely, as a consequence of transition to integrability, MBL, and Quantum Many-Body Scars (QMBS). We show that the predicted behaviour from RMT can be used as a potential witness for transition to non-ergodic behaviour from the measurement of local observables alone.
The interplay between quenched disorder and interaction effects opens the possibility in a closed quantum many-body system of a phase transition at finite energy density between an ergodic phase, which is governed by the laws of statistical physics, and a localized one, in which the degrees of freedom are frozen and ergodicity breaks down. The possible existence of a quantum phase transition at finite energy density is strongly questioning our understanding of the fundamental laws of nature and has generated an active field of research called many-body localization. This thesis consists of three parts and is dedicated to the understanding and characterization of the phenomenon of many-body localization, approaching it from complementary facets. In particular, borrowing methods and tools from different fields, we analyze timely problems. The first part of the thesis is devoted to detecting the many-body localization transition and to characterize both the ergodic and the localized phase it separates. Here we provide a characterization from two different perspectives: the first one is based on the study of local entanglement properties. In the second one, using tools from quantum-chaos theory, we attempt to answer the question of understanding time-irreversibility, and thus probing the breaking of ergodicity. We analyze experimentally viable observables. Moreover, we propose two different quantities to distinguish an Anderson insulating phase from a many-body localized one, which is one of the issues in experiments. The second part focuses on understanding the existence of a putative subdiffusive multifractal phase. Analyzing the quantum dynamics of the system in this region of the phase diagram, we point out the importance of finite-size effects, questioning the existence of this multifractal phase. We speculate with a possible scenario in which the diffusivity and thus ergodicity could be restored in the thermodynamic limit. Furthermore, we find that the propagation is highly non-Gaussian, which could have an important effect on understanding the critical point of the according transition. We tackle this problem also from a different angle. A possible toy-model to understand many-body localization entails the Anderson model on a random-regular graph. Also in the latter model the possible existence of an intermediate multifractal phase has been conjectured. There, studying the survival return probability of a particle with time, we give a new characterization of multifractal phases and give indication of the possible existence of this phase. Nevertheless, we also outline possible caveats. In the last part of this thesis we study the interplay between symmetry and correlated disorder in a non-interacting fermionic system. We show another possible mechanism for breaking localization. In particular, we focus on studying information and particle transport, emphasizing how the two types of propagation can be different.
Analog quantum simulators natively return large collections of single-shot projective measurements, the wave-function snapshots. These are high-dimensional configurations drawn from the Born distribution of the state they prepare. These collections carry far more information than is retained when they are simply compressed into low-order observables fixed in advance. Moreover, this procedure requires knowing beforehand what one is looking for. This thesis chooses to take a different route, asking what structure a dataset of snapshots possesses as a set of points in configuration space, and what kind of physics that structure reflects. The tools used here are drawn from un supervised learning and network theory, namely principal component analysis (PCA), intrinsic-dimension estimators (TWO-NN and BID), and the wave-function networks obtained by promoting snapshots to nodes and activating links according to the met ric structure of configuration space. They are chosen under two constraints, namely interpretability at every step, and the smallest possible amount of physical input: no order parameter, no symmetry, and no knowledge of the Hamiltonian. At equilibrium, wave-function networks are used to build a classification of one dimensional quantum phases based on snapshots alone. The scheme first selects, in the spirit of Occam’s razor, the measurement basis carrying the minimal information content, and then analyses the networks constructed in that basis. Paramagnetic and symmetry-broken phases display low complexity and “probability” networks, symmetry protected topological phases high complexity and Erdős–Rényi networks, and critical points high complexity and scale-free networks. Interestingly, gapped and gapless mat ter are treated on the same footing, and real-space decimation enters the scheme to resolve massive phases rather than flowing towards criticality. Out of equilibrium, the dynamical exponent governing transport of a conserved quan tity is shown to be imprinted on the singular-value spectrum of the snapshot matrix. In particular, it is shown that the PCA entropy, and in fact the leading eigenvalue alone, grows with the same exponent as the underlying spin or energy transport, at modest system sizes, evolution times and sample numbers. Ergodicity breaking is then ad dressed through the binary intrinsic dimension (BID) of the snapshots, in the quantum sun model: the BID grows and saturates at the Thouless time, reaching its maximum in the ergodic phase, remaining system-size independent in the non-ergodic one, and growing sub-maximally at criticality, where the data occupy an extensive but lower dimensional manifold. The associated networks reproduce the equilibrium dichotomy, being random in the ergodic phase and scale-free at the transition, whose location can be identified by the BID within a few percent of spectral estimates. Preliminary experimental observations, not reported here, show analogous features. Taken together, these results indicate that datasets orders of magnitude smaller than the Hilbert-space dimension suffice to identify collective properties of quantum matter through the study of the structure of the data points in configuration space. Moreover, the same structural classes recur across ground states and quenched dynamics, with no order parameter supplied and under realistic experimental imperfections. The final chapter elaborates on another aspect of quantum simulator capabilities, that is, that to realize gauge theory dynamics. There, one-dimensional U(1) lattice gauge theories with a topological θ term are mapped exactly onto kinetically con strained fermions in an emergent Stark potential. The mapping accounts for the slow relaxation reported in recent Rydberg experiments: prethermalization near an inte grable point at θ = 0, and Stark localization with Hilbert-space fragmentation at θ ̸ =0.
Non-ergodic quantum many-body dynamics offers a route to persistent quantum coherence far from equilibrium, beyond the conventional expectations of thermalizing statistical mechanics. Most known examples have been identified through intuition, analogy, or numerical search, rather than by systematic Hamiltonian design. This leaves few general methods for engineering robust non-ergodic dynamics in disorder-free interacting quantum many-body systems, especially in higher dimensions where enhanced connectivity generally disfavors ergodicity breaking. Here, we use an eigenstate-to-Hamiltonian construction approach to systematically engineer two-dimensional quantum spin Hamiltonians with tunable non-ergodic properties. Starting from a common parent model and input eigenstate, we construct two target Hamiltonians on the square lattice that differ in their interaction geometry: axial versus diagonal. We show that the resulting interaction geometry and coupling pattern control the Hilbert-space structure and dynamics. The axial model forms a connected Hilbert-space network and exhibits an anisotropy-driven crossover from ergodic to non-ergodic behavior through the suppression of resonances. In contrast, the diagonal model exhibits Hilbert-space fragmentation and, within its largest irreducible sector, signatures of quantum many-body scar dynamics including towers of highly localized, low-entanglement eigenstates, and long-lived coherent revivals. These results establish a direct link between Hamiltonian design, Hilbert-space network structure, and emergent non-ergodic dynamics.