A central theme in the field of quantum many-body physics is the characterization of thermal properties, which can be accessed, for instance, through the analysis of the Hamiltonian spectrum. This thesis approaches the problem from both numerical and analytical viewpoints, treating the spectrum as a statistical object whose correlations encode predictive information about the physics of the underlying system. In the first part of this thesis, we review the relation between thermalization, quantum integrability and chaos in one-dimensional quantum many-body systems. We focus on the statistics of energy levels and establish whether quantum integrability can be assessed from a statistical viewpoint. We then consider the problem of intermediate statistics between purely integrable and non-integrable quantum many-body systems. The first contribution of this thesis is an answer to the problem of discriminating between the spectrum of a genuine integrable system and a statistical mixture of, possibly exponentially many, non-ergodic sectors. We show that the spectral decimation algorithm can capture both qualitative and quantitative features of a `characteristic symmetry sector' from an unbiased analysis of numerically available spectral data. When applied to physically relevant systems, the algorithm is able to identify traces of Hilbert space fragmentation and disorder-induced many-body localization. The second part of this thesis is dedicated to a class of Hamiltonians of particular interest to statistical spectroscopy. Such Hamiltonians involve only the permutation of sites without the presence of external interactions. Nonetheless, they generate a highly tunable and numerically implementable class of models. For instance, it is easy to engineer integrability and its breaking while preserving essential structures such as the block-diagonal decomposition. These Hamiltonians are naturally suited to generating unitary random walks on the permutation group, which can, for instance, be compared with their classical counterparts. After a detailed discussion of the numerical implementation of these Hamiltonians and their spectral properties, we depart from the analysis of the system's static properties to answer a genuine out-of-equilibrium question: do quantum coherences accelerate mixing in the computational basis? We consider an exactly solvable class of quantum stochastic walks which, in the absence of quantum coherences, admit a limit coinciding with a classical continuous-time Markov chain. While the trace distance of the associated density matrices from the maximally mixed state remains identical to that in the classical limit, we show that the probability distance in the presence of quantum coherences, which defines mixing in the computational basis, is upper bounded by the classical limit. An analysis of the slowest mode identifies a necessary condition under which quantum coherences and classical mixing become of the same order, thereby defining a slow-mode time scale. Finally, a numerical analysis of quantum stochastic walks based on random transpositions corroborates the analytical findings, in particular the accelerated mixing. We also identify a closed-form expression for the ratio of quantum to classical mixing times, which yields a collapse over numerically accessible system sizes. Overall, this thesis presents a modern approach to the theory of statistical spectroscopy by proposing analytically grounded and numerically implementable protocols to ascertain the thermal properties of a one-dimensional quantum many-body system from the analysis of the spectrum within numerically accessible volumes. It also shows that the analysis of the statistical properties of Hamiltonians can be of practical aid in the study of out-of-equilibrium properties of quantum many-body systems, such as mixing in quantum stochastic walks.
A central theme in the field of quantum many-body physics is the characterization of thermal properties, which can be accessed, for instance, through the analysis of the Hamiltonian spectrum. This thesis approaches the problem from both numerical and analytical viewpoints, treating the spectrum as a statistical object whose correlations encode predictive information about the physics of the underlying system. In the first part of this thesis, we review the relation between thermalization, quantum integrability and chaos in one-dimensional quantum many-body systems. We focus on the statistics of energy levels and establish whether quantum integrability can be assessed from a statistical viewpoint. We then consider the problem of intermediate statistics between purely integrable and non-integrable quantum many-body systems. The first contribution of this thesis is an answer to the problem of discriminating between the spectrum of a genuine integrable system and a statistical mixture of, possibly exponentially many, non-ergodic sectors. We show that the spectral decimation algorithm can capture both qualitative and quantitative features of a `characteristic symmetry sector' from an unbiased analysis of numerically available spectral data. When applied to physically relevant systems, the algorithm is able to identify traces of Hilbert space fragmentation and disorder-induced many-body localization. The second part of this thesis is dedicated to a class of Hamiltonians of particular interest to statistical spectroscopy. Such Hamiltonians involve only the permutation of sites without the presence of external interactions. Nonetheless, they generate a highly tunable and numerically implementable class of models. For instance, it is easy to engineer integrability and its breaking while preserving essential structures such as the block-diagonal decomposition. These Hamiltonians are naturally suited to generating unitary random walks on the permutation group, which can, for instance, be compared with their classical counterparts. After a detailed discussion of the numerical implementation of these Hamiltonians and their spectral properties, we depart from the analysis of the system's static properties to answer a genuine out-of-equilibrium question: do quantum coherences accelerate mixing in the computational basis? We consider an exactly solvable class of quantum stochastic walks which, in the absence of quantum coherences, admit a limit coinciding with a classical continuous-time Markov chain. While the trace distance of the associated density matrices from the maximally mixed state remains identical to that in the classical limit, we show that the probability distance in the presence of quantum coherences, which defines mixing in the computational basis, is upper bounded by the classical limit. An analysis of the slowest mode identifies a necessary condition under which quantum coherences and classical mixing become of the same order, thereby defining a slow-mode time scale. Finally, a numerical analysis of quantum stochastic walks based on random transpositions corroborates the analytical findings, in particular the accelerated mixing. We also identify a closed-form expression for the ratio of quantum to classical mixing times, which yields a collapse over numerically accessible system sizes. Overall, this thesis presents a modern approach to the theory of statistical spectroscopy by proposing analytically grounded and numerically implementable protocols to ascertain the thermal properties of a one-dimensional quantum many-body system from the analysis of the spectrum within numerically accessible volumes. It also shows that the analysis of the statistical properties of Hamiltonians can be of practical aid in the study of out-of-equilibrium properties of quantum many-body systems, such as mixing in quantum stochastic walks.
Abstract Integrability has long served as a cornerstone concept in classical mechanics, where it possesses a precise and unambiguous definition. Extending this notion to the quantum domain, however, remains a far more subtle and elusive problem. In particular, deciding whether a given quantum Hamiltonian—viewed simply as a matrix of its elements—does or does not define an integrable system is far from obvious. Yet this question is crucial: it bears directly on non-equilibrium dynamics, spectral correlations, the behaviour of correlation functions, and other fundamental properties of many-body quantum systems. In this work, we develop a statistical framework for addressing quantum integrability from a purely probabilistic standpoint. Our approach begins with the observation that a necessary signature of integrability is the finite probability of encountering vanishing energy gaps in the spectrum. On this basis, we formulate a twofold protocol capable of distinguishing between integrable and non-integrable Hamiltonians. The first step consists of a systematic Monte Carlo decimation of the spectrum, designed to reveal the emergence (or absence) of Poissonian level spacing statistics. The iterative decimation compresses the Hilbert space exponentially, and its termination point determines whether the spectrum is governed by a mixed distribution or approaches the Poisson limit. In the second step, the distinction is instead obtained by analysing the distributions of k -step gaps which can help in discriminating between Poisson and mixed statistics. This procedure applies to Hamiltonians of arbitrary finite size, regardless of whether their structure involves a finite number of blocks or an exponentially fragmented Hilbert space. As a concrete benchmark, we implement the protocol on a class of quantum Hamiltonians constructed from the permutation group S N , thereby demonstrating both its effectiveness and its broad applicability.
We develop a systematic theory of spectral decimation for quantum many-body Hamiltonians and show that it provides a quantitative probe of emergent symmetries in statistically mixed spectra. Building on an analytical description of statistical mixtures, we derive an explicit expression for the size of a characteristic symmetry sector (CSS), defined as the largest subsequence of levels exhibiting non-Poissonian correlations. The CSS dimension is shown to be the size-biased average of the underlying symmetry sectors, establishing a direct link between spectral statistics and Hilbert-space structure. We apply this framework to two paradigmatic settings: Hilbert-space fragmentation and disorder-induced many-body localization (MBL). In fragmented systems, the CSS reproduces the mixture prediction and isolates correlated subsectors even when the full spectrum appears nearly Poissonian. In the disordered Heisenberg chain, spectral decimation reveals the gradual emergence of integrability through a shrinking CSS, whose statistics exhibit signatures consistent with local integrals of motion. We introduce a characteristic symmetry entropy (CSE) as a finite-size scaling observable and extract, within accessible system sizes, the crossover exponents. Our results establish spectral decimation as a controlled, unbiased and computationally inexpensive diagnostic of hidden structure in many-body spectra, capable of distinguishing between chaotic dynamics, statistical mixtures, and emergent integrability.
A key goal of quantum chaos is to establish a relationship between widely observed universal spectral fluctuations of clean quantum systems and random matrix theory (RMT).Most prominent features of such RMT behavior with respect to a random spectrum, both encompassed in spectral pair correlation function, are statistical suppression of small level spacings (correlation hole) and enhanced stiffness of the spectrum at large spectral ranges.For single particle systems with fully chaotic classical counterparts, the problem has been partly solved by Berry [Proc.R. Soc.London, A400, 229 (1985)] within the so-called diagonal approximation of semiclassical periodic-orbit sums, while the derivation of the full RMT spectral form factor K(t) (Fourier transform of spectral pair correlation function), from semiclassics has been completed by Müller et al. [Phys.Rev. Lett.93, 014103 (2004)].In recent years, the questions of long-time dynamics at high energies, for which the full many-body energy spectrum becomes relevant, are coming at the forefront even for simple many-body quantum systems, such as locally interacting spin chains.Such systems display two universal types of behaviour which are termed as the 'many-body localized phase' and 'ergodic phase'.In the ergodic phase, the spectral fluctuations are excellently described by RMT, even for very simple interactions and in the absence of any external source of disorder.Here we provide a clear theoretical explanation for these observations.We compute K(t) in the leading two orders in t and show its agreement with RMT for non-integrable, time-reversal invariant many-body systems without classical counterparts, a generic example of which are Ising spin 1/2 models in a periodically kicking transverse field.In particular, we relate K(t) to partition functions of a class of twisted classical Ising models on a ring of size t, hence the leading order RMT behavior K(t)2t is a consequence of translation and reflection symmetry of the Ising partition function.
This dissertation investigates the theory of quantum stochastic processes and its applications in quantum many-body physics. The main goal is to analyse complexity-theoretic aspects of both static and dynamic properties of physical systems modelled by quantum stochastic processes. The thesis consists of two parts: the first one addresses the computational complexity of certain quantum and classical divisibility questions, whereas the second one addresses the topic of Hamiltonian complexity theory. In the divisibility part, we discuss the question whether one can efficiently sub-divide a map describing the evolution of a system in a noisy environment, i.e. a CPTP- or stochastic map for quantum and classical processes, respectively, and we prove that taking the nth root of a CPTP or stochastic map is an NP-complete problem. Furthermore, we show that answering the question whether one can divide up a random variable X into a sum of n iid random variables Yi, i.e. X=∑i=1nYi, is poly-time computable; relaxing the iid condition renders the problem NP-hard. In the local Hamiltonian part, we study computation embedded into the ground state of a many-body quantum system, going beyond "history state" constructions with a linear clock. We first develop a series of mathematical techniques which allow us to study the energy spectrum of the resulting Hamiltonian, and extend classical string rewriting to the quantum setting. This allows us to construct the most physically-realistic QMAEXP-complete instances for the LOCAL HAMILTONIAN problem (i.e. the question of estimating the ground state energy of a quantum many-body system) known to date, both in one- and three dimensions. Furthermore, we study weighted versions of linear history state constructions, allowing us to obtain tight lower and upper bounds on the promise gap of the LOCAL HAMILTONIAN problem in various cases. We finally study a classical embedding of a Busy Beaver Turing Machine into a low-dimensional lattice spin model, which allows us to dictate a transition from a purely classical phase to a Toric Code phase at arbitrarily large and potentially even uncomputable system sizes.
Phenomena in quantum chaotic systems such as spectral fluctuations are known to be described remarkably well by the random matrix theory (RMT) in small energy intervals and over large time scales. However, in exchange for analytical tractability, RMT has been constructed to leave out an important element of realistic physical systems { the notion of locality. Since quantum information propagates at a finite speed in a spatially-extended quantum chaotic system, there exists time scales before which RMT is effective. How do we describe the non-equilibrium quantum dynamics at intermediate timescales which are much longer than the local relaxation time but too short for RMT to apply? What are the signatures of many-body quantum chaotic systems beyond the description of RMT? In this thesis, we address these questions by analysing minimal models that incorporate a notion of locality and at the same time retain some analytical tractability of RMT. More specifically, we study Floquet quantum circuits that are random in space but periodic in time. These circuits consist of local and random quantum gates that act on a q-state `spin' at each site. After a brief review of classical chaos, quantum chaos, and quantum circuits in Chapter 1, we present three main results. In Chapter 2, we formulate a diagrammatical approach for computing ensemble averages of observables over random unitary circuits, and derive exact results on spectral statistics and dynamical observables for a Floquet circuit in the large-q limit. In Chapter 3, we further investigate the spectral statistics in spatiallyextended many-body quantum chaotic systems by computing the spectral form factor K(t) analytically and numerically. We show that K(t) follows the RMT description at times longer than the Thouless time tTh and identify the scaling of tTh at the large-q limit of a Floquet circuit. In Chapter 4, we discuss eigenstate correlations for spatially-extended many-body quantum chaotic systems, in terms of the statistical properties of matrix elements of local observables. While the eigenstate thermalization hypothesis (ETH) is known to give an excellent description of these quantities, the phenomenon of scrambling and the butterfly effect imply structure beyond ETH. We determine analytically and verify numerically the universal form of this structure at long distances and small eigenvalue separations for Floquet systems.
This thesis investigates the dynamics of quantum many-body systems in isolation, in three different setups. We start by considering a paradigmatic quantum chaotic model of magnetism, the Heisenberg antiferromagnet in 2 and 3 spatial dimensions. The focus is on obtaining approximations for the full counting statistics (FCS) of a subsystem order parameter, the staggered magnetization, both in and out of equilibrium (following quantum quenches). We employ mappings of spin variables to bosonic ones, together with self-consistent time-dependent mean-field schemes, to compute the FCS via Gaussian techniques. In particular, we introduce a formula for the FCS of any subsystem observable that is at most quadratic in the bosons. We obtain results both in the presence and absence of long-range magnetic order. We then turn to the first of two ergodicity-breaking scenarios, by considering Hilbert space fragmentation (HSF) in charge- and dipole-conserving quantum chains. We determine the critical charge density at which a "freezing" phase transition from weak to strong fragmentation takes place, for any on-site Hilbert space dimension. To characterize the degree of dynamical connectivity in the system, we introduce a novel framework that enables us to obtain analytic results for the strongly fragmented phase, and numerically explore the weakly fragmented one. We derive some critical exponents, and show that the entanglement entropy of typical eigenstates in the strongly fragmented phase obeys an area law. We also uncover the existence of rare eigenstates whose entropy appears to scale beyond any area law. Finally, we investigate dynamical correlators in the Lieb-Liniger gas, one of the simplest interacting integrable theories solvable by Bethe ansatz. We employ Markov chain Monte Carlo methods to sample matrix elements (form factors) that appear in the Lehmann representation of the finite-temperature single-particle bosonic Green’s function. In the physical regimes we consider, sampling is necessary because a full summation of the relevant form factors is computationally intractable. We benchmark our results against exact solutions available in specific limits, and find them to be very accurate.