The occupation number is a key observable for diagnosing thermalization, as it connects directly to standard statistical laws such as Fermi-Dirac, Bose-Einstein, and Boltzmann distributions. In the context of spin systems, it represents the population of the sublevels of the magnetization in the z direction. We use this quantity to probe the onset of thermalization in an isolated one-dimensional quantum spin-1 Ising model with transverse and longitudinal fields and in its classical counterpart. Thermalization is achieved when the long-time average of the occupation number converges to the microcanonical prediction as the chain length L increases, consistent with the emergence of ergodicity. However, the finite-size scaling analysis in the quantum model is challenged by the exponential growth of the Hilbert space with L. To overcome this limitation, we turn to the classical model, which enables access to much larger system sizes. By tracking the dynamics of individual spins on their three-dimensional Bloch spheres and employing tools from random matrix theory, we establish a quantitative criterion for classical ergodicity in interacting spin systems. We find that deviations from classical ergodicity decay algebraically with system size. This power-law scaling then provides a quantitative bound on the approach to thermal equilibrium in the quantum model.
The occupation number is a key observable for diagnosing thermalization, as it connects directly to standard statistical laws such as Fermi-Dirac, Bose-Einstein, and Boltzmann distributions. In the context of spin systems, it represents the population of the sublevels of the magnetization in the z direction. We use this quantity to probe the onset of thermalization in an isolated one-dimensional quantum spin-1 Ising model with transverse and longitudinal fields and in its classical counterpart. Thermalization is achieved when the long-time average of the occupation number converges to the microcanonical prediction as the chain length L increases, consistent with the emergence of ergodicity. However, the finite-size scaling analysis in the quantum model is challenged by the exponential growth of the Hilbert space with L. To overcome this limitation, we turn to the classical model, which enables access to much larger system sizes. By tracking the dynamics of individual spins on their three-dimensional Bloch spheres and employing tools from random matrix theory, we establish a quantitative criterion for classical ergodicity in interacting spin systems. We find that deviations from classical ergodicity decay algebraically with system size. This power-law scaling then provides a quantitative bound on the approach to thermal equilibrium in the quantum model.
The occupation number is a key observable for diagnosing thermalization, as it connects directly to standard statistical laws such as Fermi--Dirac, Bose--Einstein, and Boltzmann distributions. In the context of spin systems, it represents the population of the sublevels of the magnetization in the z-direction. We use this quantity to probe the onset of thermalization in the isolated quantum and classical one-dimensional spin-1 Ising model with transverse and longitudinal fields. Thermalization is achieved when the long-time average of the occupation number converges to the microcanonical prediction as the chain length L increases, consistent with the emergence of ergodicity. However, the finite-size scaling analysis in the quantum model is challenged by the exponential growth of the Hilbert space with L. To overcome this limitation, we turn to the corresponding classical model, which enables access to much larger system sizes. By tracking the dynamics of individual spins on their three-dimensional Bloch spheres and employing tools from random matrix theory, we establish a quantitative criterion for classical ergodicity in interacting spin systems. We find that deviations from classical ergodicity decay algebraically with system size. This power-law scaling then provides a quantitative bound on the approach to thermal equilibrium in the quantum model.
Abstract A study of quantum-classical correspondence for the onset of thermalization in an isolated system of four interacting particles is presented. Provided that classical dynamics of the system is chaotic, thermalization is present and fits the expected statistical mechanics results. Quantum classical correspondence is established for the classical analogue of the single particle distributions of occupation numbers. The latter quantity is useful to test thermalization and it allows to obtain the microcanonical Boltzmann factor in the quantum and semi-classical cases. The correspondence between them and the statistical mechanics results are surprisingly good considering that we have only four interacting particles. The resulting Boltzmann factor for the occupation number distribution depends on the eigenfunction width, which is a semi-classical parameter explicitly calculated and consistent with the quantum one.
We investigate the relation between the classical ergodicity and the quantum eigenstate thermalization in the fully connected Ising ferromagnets. In the case of spin-1/2, an expectation value of an observable in a single-energy eigenstate coincides with the long-time average in the underlying classical dynamics, which is a consequence of the Wentzel-Kramers-Brillouin approximation. In the case of spin-1, the underlying classical dynamics is not necessarily ergodic. In that case, it turns out that, in the thermodynamic limit, the statistics of the expectation values of an observable in the energy eigenstates coincides with the statistics of the long-time averages in the underlying classical dynamics starting from random initial states sampled uniformly from the classical phase space. This feature seems to be a general property in semiclassical systems, and the result presented here is crucial in discussing equilibration, thermalization, and dynamical transitions of such systems.
From the analysis of the relaxation process of isolated lattice many-body quantum systems quenched far from equilibrium, we deduce a criterion for predicting when they are certain to thermalize. It is based on the algebraic behavior $\ensuremath{\propto}{t}^{\ensuremath{-}\ensuremath{\gamma}}$ of the survival probability at long times. We show that the value of the power-law exponent $\ensuremath{\gamma}$ depends on the shape and filling of the weighted energy distribution of the initial state. Two scenarios are explored in detail: $\ensuremath{\gamma}\ensuremath{\ge}2$ and $\ensuremath{\gamma}<1$. Exponents $\ensuremath{\gamma}\ensuremath{\ge}2$ imply that the energy distribution of the initial state is ergodically filled and the eigenstates are uncorrelated, so thermalization is guaranteed to happen. In this case, the power-law behavior is caused by bounds in the energy spectrum. Decays with $\ensuremath{\gamma}<1$ emerge when the energy eigenstates are correlated and signal lack of ergodicity. They are typical of systems undergoing localization due to strong onsite disorder and are found also in clean integrable systems.
In this Thesis we discussed several issues related to the dynamics of isolated quantum many-body systems. The brief overview of Chapter 1 highlighted a series of open fundamental theoretical questions which are crucial for the investigation of nonequilibrium phenomena. Among them, we mainly focused on the emergence of thermal behaviors and on the characterization of dynamical phase transitions and of the associated universal dynamics. To this end, we analyzed some paradigmatic models, such as the Fermi-Hubbard model, the quantum Ising chain, and the O(N) model. -Chapter 2: Absence of thermalization in a Fermi liquid. In this Chapter we report two results that question the common expectation of thermalization in a Landau-Fermi liquid after an interaction quench. We first consider the perturbative expansion in the interaction strength of the long-wavelength structure factor S(q) and show that it does not satisfy the hypothesis that steady-state averages correspond to thermal ones. In particular, S(q) has an analytical component ∼ const. + O(q^2) for q → 0, compatible with thermalization, but it retains also a nonanalytic term ∼ |q| characteristic of a Fermi liquid at zero-temperature. In real space, this nonanalyticity turns in a power-law decay of the density-density correlations, in contrast with the exponential decay associated to thermalization. We next analyze the low-density case, where one can obtain results valid at any order in interaction but at leading in the density, and find that in the steady state the momentum distribution jump at the Fermi surface is strictly finite, though smaller than at equilibrium. -Chapter 3: Nonadiabatic stationary behavior in a driven Ising chain. In this Chapter we discuss the emergence of a nonadiabatic behavior in the dynamics of the order parameter in a low-dimensional system subject to a linear ramp of one of its parameters within a gapped phase, which should be the most favorable situation for an adiabatic evolution. We study in details the dynamics of the spontaneous magnetization m_x(t) in a quantum Ising chain after a linear variation in time of the transverse field within the ordered phase. In particular, focusing on the value of the order parameter at the end of the ramp m_x(τ), we show that the smaller the switching rate of the transverse field is the closer m_x(τ) gets to its ground state value. Nonetheless, this apparently small correction to adiabaticity eventually leads to the disruption of the order exponentially fast in the subsequent time evolution, no matter how slowly the ramp is performed. -Chapter 4: Aging and coarsening in isolated quantum systems after a quench. This Chapter analyzes the nonequilibrium dynamics of an isolated quantum system after a sudden quench to the dynamical critical point, where the emergence of scaling and universal exponents is expected, due to the absence of time and length scales. We explore these features for a bosonic interacting scalar field theory with O(N) symmetry in the large-N limit, where the model is exactly solvable and the exponents and scaling forms of the relevant two-time correlation functions can be analytically calculated. Moreover, we provide evidence of the emergence of a dynamical scaling behavior also for quenches below the dynamical critical point, associated with coarsening. We find that the latter case is characterized by the same scaling functions as those describing the critical case, yet with different exponents. -Chapter 5: Dynamical transitions and statistics of excitations. In this Chapter we study the dynamics of the O(N) model (introduced in the previous Chapter) resulting from a different protocol: a linear ramp of one of its parameters. We find that the presence of a dynamical phase transition, as well as its critical properties, are robust against the change of the protocol. We show that a characterization based on the critical dimensions and exponents would suggest that the dynamical phase transition is analogous to the equilibrium thermal one. However, its nonequilibrium nature becomes evident analyzing the statistics of excitations produced in the ramp process. In particular, the critical properties are encoded in the fluctuations in the number of excitations, which display qualitatively different behaviors depending on the ramp being performed above, at, or below the dynamical critical point. These behaviors bear no dependence on the duration of the protocol.
Ergodicity in quantum many-body systems is---despite its fundamental importance---still an open problem. Many-body localization provides a general framework for quantum ergodicity and may therefore offer important insights. However, the characterization of many-body localization through simple observables is a difficult task. In this article, we introduce a measure for distances in Hilbert space for spin-21 systems that can be interpreted as a generalization of the Anderson localization length to many-body Hilbert space. We show that this many-body localization length is equivalent to a simple local observable in real space, which can be measured in experiments of superconducting qubits, polar molecules, Rydberg atoms, and trapped ions. By using the many-body localization length and a necessary criterion for ergodicity that it provides, we study many-body localization and quantum ergodicity in power-law-interacting Ising models subject to disorder in a transverse field. Based on the nonequilibrium dynamical renormalization group, numerically exact diagonalization, and an analysis of the statistics of resonances, we find a many-body localized phase at infinite temperature for small power-law exponents. Within the applicability of these methods, we find no indications of a delocalization transition.