Search

Due to a phenomenon of many-body localization (MBL), the strong disorder may significantly slow down or even completely hinder the thermalization of quantum many-body systems. A sufficiently deep quasiperiodic potential may also inhibit thermalization. In this work, we numerically demonstrate direct differences in the behavior of standard ergodicity-breaking indicators at the MBL crossover in random and quasiperiodic systems. Our key finding is the exponential increase in the sharpness of the MBL crossover with system size for quasiperiodic systems, a trend that is only linear in disordered systems. The strong tendency towards a non-analytic behavior in quasiperiodic systems is consistent with the existence of dynamical regimes with sharply defined boundaries or an MBL phase transition. It highlights the importance of quasiperiodic systems for our understanding of many-body dynamics., Comment: 5pp+supplement, Comments welcome!!

Magic state resources or non-stabilizerness quantify the beyond-Clifford operations necessary for universal quantum computing. How rapidly are magic resources generated by generic many-body dynamics under constraints of locality? We address this problem by exploring magic spreading in brick-wall random unitary circuits. Inspired by the algebraic structure of the Clifford group, we propose a scalable measure of non-stabilizerness, the Calderbank-Shor-Steane entropy, which generalizes the notion of stabilizer entropy and mirrors its qualitative behavior. This metric enables the investigation of non-stabilizerness dynamics for systems of up to N = 1024 qudits. Our main finding is that magic resources equilibrate on timescales logarithmic in system size N, akin to anticoncentration and Hilbert space delocalization measures, but differently from entanglement entropy. We conjecture that our findings describe the phenomenology of non-stabilizerness growth in a broad class of chaotic many-body systems., Comment: 4.25 pages + 6. Comments are welcome!

Many-body localization (MBL) hinders the thermalization of quantum many-body systems in the presence of strong disorder. In this work, we study the MBL regime in bond-disordered spin-1/2 XXZ spin chain, finding the multimodal distribution of entanglement entropy in eigenstates, sub-Poissonian level statistics, and revealing a relation between operators and initial states required for examining the breakdown of thermalization in the time evolution of the system. We employ a real space renormalization group scheme to identify these phenomenological features of the MBL regime that extend beyond the standard picture of local integrals of motion relevant for systems with disorder coupled to on-site operators. Our results pave the way for experimental probing of MBL in bond-disordered spin chains., Comment: comments most welcome!!!

Inverse Participation Ratios (IPRs) and the related Participation Entropies quantify the spread of a quantum state over a selected basis of the Hilbert space, offering insights into the equilibrium and non-equilibrium properties of the system. In this work, we propose three quantum algorithms to estimate IPRs on multi-qubit and multi-qudit quantum devices. The first algorithm allows for the estimation of IPRs in the computational basis by single-qubit measurements, while the second one enables measurement of IPR in the eigenbasis of a selected Hamiltonian, without the knowledge about the eigenstates of the system. Next, we provide an algorithm for IPR in the computational basis for a multi-qudit system. We discuss resources required by the algorithms and benchmark them by investigating the one-axis twisting protocol, the thermalization in a deformed PXP model, and the ground state of a spin-$1$ AKLT chain in a transverse field.

Unitary dynamics of a quantum system initialized in a selected basis state yields, generically, a state that is a superposition of all the basis states. This process, associated with the quantum information scrambling and intimately tied to the resource theory of coherence, may be viewed as a gradual delocalization of the system's state in the Hilbert space. This work analyzes the Hilbert space delocalization under dynamics of random quantum circuits, which serve as a minimal model of chaotic dynamics of quantum many-body systems. We employ analytical methods based on the replica trick and Weingarten calculus to investigate the time evolution of the participation entropies which quantify the Hilbert space delocalization. We demonstrate that the participation entropies approach, up to a fixed accuracy, their long-time saturation value in times that scale logarithmically with the system size. Exact numerical simulations and tensor network techniques corroborate our findings., Comment: 17 pages, 3 figures. Invited Feature Paper for Entropy. Comments are welcome!

Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. Persistent finite size drifts towards ergodicity consistently emerge in spectral properties of disordered many-body systems, excluding naive single-parameter scaling hypothesis and preventing comprehension of the status of the MBL phase. The drifts are related to tendencies towards thermalization and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength suggests the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations., Comment: 72 pages, 14 figures, several references added and discussions clarified. Comments are welcome!

We discuss the dependence of the critical properties of the Anderson model on the dimension $d$ in the language of $\beta$-function and renormalization group recently introduced in Ref.[arXiv:2306.14965] in the context of Anderson transition on random regular graphs. We show how in the delocalized region, including the transition point, the one-parameter scaling part of the $\beta$-function for the fractal dimension $D_{1}$ evolves smoothly from its $d=2$ form, in which $\beta_2\leq 0$, to its $\beta_\infty\geq 0$ form, which is represented by the regular random graph (RRG) result. We show how the $\epsilon=d-2$ expansion and the $1/d$ expansion around the RRG result can be reconciled and how the initial part of a renormalization group trajectory governed by the irrelevant exponent $y$ depends on dimensionality. We also show how the irrelevant exponent emerges out of the high-gradient terms of expansion in the nonlinear sigma-model and put forward a conjecture about a lower bound for the fractal dimension. The framework introduced here may serve as a basis for investigations of disordered many-body systems and of more general non-equilibrium quantum systems., Comment: 13 pages, 12 figures. Comments are welcome!

We study the XXZ model with a random magnetic field in contact with a weakly disordered spin chain, acting as a finite thermal bath. We revise Fermi's golden rule description of the interaction between the thermal bath and the XXZ spin chain, contrasting it with a nonperturbative quantum avalanche scenario for the thermalization of the system. We employ two-point correlation functions to define the extent $\xi_d$ of the thermalized region next to the bath. Unbounded growth of $\xi_d$ proportional to the logarithm of time or faster is a signature of an avalanche. Such behavior signifies the thermalization of the system, as we confirm numerically for a generic initial state in the ergodic and critical regimes of the XXZ spin chain. In the many-body localized regime, a clear termination of avalanches is observed for specifically prepared initial states and, surprisingly, is not visible for generic initial product states. Additionally, we extract the localization length of the local integrals of motion and show that a bath made out of a weakly disordered XXZ chain has a similar effect on the system as a bath modeled by a Hamiltonian from a Gaussian orthogonal ensemble of random matrices. We also comment on the result of the earlier study (Phys. Rev. B 108, L020201 (2023)), arguing that the observed thermalization is due to external driving of the system and does not occur in the autonomous model. Our work reveals experimentally accessible signatures of quantum avalanches and identifies conditions under which termination of the avalanches may be observed., Comment: Authors' accepted manuscript

An important question of quantum information is to characterize genuinely quantum (beyond-Clifford) resources necessary for universal quantum computing. Here, we use the Pauli spectrum to quantify how magic, beyond Clifford, typical many-qubit states are. We first present a phenomenological picture of the Pauli spectrum based on quantum typicality and then confirm it for Haar random states. We then introduce filtered stabilizer entropy, a magic measure that can resolve the difference between typical and atypical states. We proceed with the numerical study of the Pauli spectrum of states created by random circuits as well as for eigenstates of chaotic Hamiltonians. We find that in both cases the Pauli spectrum approaches the one of Haar random states, up to exponentially suppressed tails. Our results underscore differences between typical and atypical states from the point of view of quantum information., Comment: 4.5 pages + 4 SM pages, Comments are welcome!

The Gorini-Kossakowski-Sudarshan-Lindblad master equation (ME) governs the density matrix of open quantum systems (OQSs). When an OQS is subjected to weak continuous measurement, its state evolves as a stochastic quantum trajectory, whose statistical average solves the ME. The ensemble of such trajectories is termed an unraveling of the ME. We propose a method to operationally distinguish unravelings produced by the same ME in different measurement scenarios, using nonlinear averages of observables over trajectories. We apply the method to the paradigmatic quantum nonlinear system of resonance fluorescence in a two-level atom. We compare the Poisson-type unraveling, induced by direct detection of photons scattered from the two-level emitter, and the Wiener-type unraveling, induced by phase-sensitive detection of the emitted field. We show that a quantum-trajectory-averaged variance is able to distinguish these measurement scenarios. We evaluate the performance of the method, which can be readily extended to more complex OQSs, under a range of realistic experimental conditions., Comment: 5 pages, 3 figures, with supplementary material, revised version following the published version

The eigenstate thermalization hypothesis (ETH) is a successful theory that establishes the criteria for ergodicity and thermalization in isolated quantum many-body systems. In this work, we investigate the thermalization properties of spin-$ 1/2 $ XXZ chain with linearly-inhomogeneous interactions. We demonstrate that introduction of the inhomogeneous interactions leads to an onset of quantum chaos and thermalization, which, however, becomes inhibited for sufficiently strong inhomogeneity. To exhibit ETH, and to display its breakdown upon varying the strength of interactions, we probe statistics of energy levels and properties of matrix elements of local observables in eigenstates of the inhomogeneous XXZ spin chain. Moreover, we investigate the dynamics of the entanglement entropy and the survival probability which further evidence the thermalization and its breakdown in the considered model. We outline a way to experimentally realize the XXZ chain with linearly-inhomogeneous interactions in systems of ultracold atoms. Our results highlight a mechanism of emergence of ETH due to insertion of inhomogeneities in an otherwise integrable system and illustrate the arrest of quantum dynamics in presence of strong interactions., Comment: 13 pages, 10 figures

We study the influence of feedback operations on the dynamics of $(d+1)$-dimensional monitored random quantum circuit. Competition between unitary dynamics and measurements leads to an entanglement phase transition, while the feedback steers the dynamics towards an absorbing state, yielding an absorbing state phase transition. Building on previous results in one spatial dimension [Phys. Rev. Lett. 130, 120402 (2023)], we discuss the interplay between the two types of transitions for $d \ge 2$ in the presence of (i) short-range feedback operations or (ii) additional global control operations. In both cases, the absorbing state transition belongs to the $d$-dimensional directed percolation universality class. In contrast, the entanglement transition depends on the feedback operation type and reveals the dynamics' inequivalent features. The entanglement and absorbing state phase transition remain separated for short-range feedback operations. When global control operations are applied, we find the two critical points coinciding; nevertheless, the universality class may still differ, depending on the choice of the control operation., Comment: 10 pages, Comments are welcome! "11th Workshop on Quantum Chaos and Localisation Phenomena" Warsaw, May 2023

Understanding how errors deteriorate the information encoded in a many-body quantum system is a fundamental problem with practical implications for quantum technologies. Here, we investigate a class of encoding-decoding random circuits subject to local coherent and incoherent errors. We analytically demonstrate the existence of a phase transition from an error-protecting phase to an error-vulnerable phase occurring when the error strength is increased. This transition is accompanied by R\'enyi entropy transitions and by onset of multifractal features in the system. Our results provide a new perspective on storing and processing quantum information, while the introduced framework enables an analytic understanding of a dynamical critical phenomenon in a many-body system., Comment: Revised version, incoherent and site-dependent errors included, 4+14 pages, comments welcome!

The quantum sun model is an interacting model that exhibits sharp signatures of ergodicity breaking phase transition. Here, we show that the model exhibits a many-body mobility edge. We provide analytical arguments for its existence, complemented by the state-of-the-art numerical simulations analysing gap ratios, Thouless times as well as entanglement entropy of eigenstates. We also introduce the quantum sun model with particle number conservation, and we argue that it shares many similarities with his unrestricted predecessor., Comment: 4pp+suppl, version close to accepted in PRB Lett

Understanding the nature of entanglement growth in many-body systems is one of the fundamental questions in quantum physics. Here, we study this problem by characterizing the entanglement fluctuations and distribution of $(d+1)$ qubit lattice evolved under a random unitary circuit. Focusing on Clifford gates, we perform extensive numerical simulations of random circuits in $1\le d\le 4$ dimensions. Our findings demonstrate that properties of growth of bipartite entanglement entropy are characterized by the roughening exponents of a $d$-dimensional membrane in a $(d+1)$ elastic medium., Comment: 4 pages

Universal quantum computing requires nonstabilizer (magic) quantum states. Quantifying the nonstabilizerness and relating it to other quantum resources is vital for characterizing the complexity of quantum many-body systems. In this work, we prove that a quantum state is a stabilizer if and only if all states belonging to its Clifford orbit have a flat probability distribution on the computational basis. This implies, in particular, that multifractal states are nonstabilizers. We introduce multifractal flatness, a measure based on the participation entropy that quantifies the wave-function distribution flatness. We demonstrate that this quantity is analytically related to the stabilizer entropy of the state and present several examples elucidating the relationship between multifractality and nonstabilizerness. In particular, we show that the multifractal flatness provides an experimentally and computationally viable nonstabilizerness certification. Our work unravels the direct relation between the nonstabilizerness of a quantum state and its wave-function structure., Comment: 6 pages, Comments are welcome!

Many-body unitary dynamics interspersed with repeated measurements display a rich phenomenology hallmarked by measurement-induced phase transitions. Employing feedback-control operations that steer the dynamics toward an absorbing state, we study the entanglement entropy behavior at the absorbing state phase transition. For short-range control operations, we observe a transition between phases with distinct sub-extensive scalings of entanglement entropy. In contrast, the system undergoes a transition between volume-law and area-law phases for long-range feedback operations. The fluctuations of entanglement entropy and of the order parameter of the absorbing state transition are fully coupled for sufficiently strongly entangling feedback operations. In that case, entanglement entropy inherits the universal dynamics of the absorbing state transition. This is, however, not the case for arbitrary control operations, and the two transitions are generally distinct. We quantitatively support our results by introducing a framework based on stabilizer circuits with classical flag labels. Our results shed new light on the problem of observability of measurement-induced phase transitions., Comment: 4+6pp, comments welcome!

We investigate dynamics of a single mobile impurity immersed in a bath of Anderson localized particles and focus on the regime of relatively strong disorder and interactions. In that regime, the dynamics of the system is particularly slow, suggesting, at short times, an occurrence of many-body localization. Considering longer time scales, we show that the latter is a transient effect and that, eventually, the impurity spreads sub-diffusively and induces a gradual delocalization of the Anderson insulator. The phenomenology of the system in the considered regime of slow dynamics includes a sub-diffusive growth of mean square displacement of the impurity, power-law decay of density correlation functions of the Anderson insulator and a power-law growth of entanglement entropy in the system. We observe a similar regime of slow dynamics also when the disorder in the system is replaced by a sufficiently strong quasi-periodic potential., Comment: AAM version to be published in PRB

Using local density correlation functions for a one-dimensional spin system, we introduce a correlation function difference (CFD) which compares correlations on a given site between a full system of size $L$ and its restriction to $\ell

The interplay between unitary dynamics and local quantum measurements results in unconventional non-unitary dynamical phases and transitions. In this paper we investigate the dynamics of $(d+1)$-dimensional hybrid stabilizer circuits, for $d=1,2,3$. We characterize the measurement-induced phases and their transitions using large-scale numerical simulations focusing on entanglement measures, purification dynamics, and wave-function structure. Our findings demonstrate the measurement-induced transition in $(d+1)$ spatiotemporal dimensions is conformal and close to the percolation transition in $(d+1)$ spatial dimensions., Comment: 19pp, 16 figs, comments welcome!

We perform a thorough and complete analysis of the Anderson localization transition on several models of random graphs with regular and random connectivity. The unprecedented precision and abundance of our exact diagonalization data (both spectra and eigenstates), together with new finite size scaling and statistical analysis of the graph ensembles, unveils a universal behavior which is described by two simple, integer, scaling exponents. A by-product of such analysis is a reconciliation of the tension between the results of perturbation theory coming from strong disorder and earlier numerical works, which seemed to suggest that there should be a non-ergodic region above a given value of disorder $W_{E}$ which is strictly less than the Anderson localization critical disorder $W_C$, and that of other works which suggest that there is no such region. We find that, although no separate $W_{E}$ exists from $W_C$, the length scale at which fully developed ergodicity is found diverges like $|W-W_C|^{-1}$, while the critical length over which delocalization develops is $\sim |W-W_C|^{-1/2}$. The separation of these two scales at the critical point allows for a true non-ergodic, delocalized region. In addition, by looking at eigenstates and studying leading and sub-leading terms in system size-dependence of participation entropies, we show that the former contain information about the non-ergodicity volume which becomes non-trivial already deep in the delocalized regime. We also discuss the quantitative similarities between the Anderson transition on random graphs and many-body localization transition., Comment: v3, expanded discussions, 36 pages + references, 17 figures, comments welcome

We study many-body localization (MBL) transition in disordered Floquet systems using a polynomially filtered exact diagonalization (POLFED) algorithm. We focus on disordered kicked Ising model and quantitatively demonstrate that finite size effects at the MBL transition are less severe than in the random field XXZ spin chains widely studied in the context of MBL. Our conclusions extend also to other disordered Floquet models, indicating smaller finite size effects than those observed in the usually considered disordered autonomous spin chains. We observe consistent signatures of the transition to MBL phase for several indicators of ergodicity breaking in the kicked Ising model. Moreover, we show that an assumption of a power-law divergence of the correlation length at the MBL transition yields a critical exponent $\nu \approx 2$, consistent with the Harris criterion for 1D disordered systems., Comment: version accepted in PRB, clarifications in Eq. 3 and Fig. 5 added, comments welcome

Highly excited eigenstates of quantum many-body systems are typically featureless thermal states. Some systems, however, possess a small number of special, low-entanglement eigenstates known as quantum scars. We introduce a quantum-inspired machine learning platform based on a Quantum Variational Autoencoder (QVAE) that detects families of scar states in spectra of many-body systems. Unlike a classical autoencoder, QVAE performs a parametrized unitary operation, allowing us to compress a single eigenstate into a smaller number of qubits. We demonstrate that the autoencoder trained on a scar state is able to detect the whole family of scar states sharing common features with the input state. We identify families of quantum many-body scars in the PXP model beyond the $\mathbb{Z}_2$ and $\mathbb{Z}_3$ families and find dynamically decoupled subspaces in the Hilbert space of disordered, interacting spin ladder model. The possibility of an automatic detection of subspaces of scar states opens new pathways in studies of models with a weak breakdown of ergodicity and fragmented Hilbert spaces., Comment: Author accepted manuscript to be published in PRB

We study time dynamics of 1D disordered Heisenberg spin-1/2 chain focusing on a regime of large system sizes and a long time evolution. This regime is relevant for observation of many-body localization (MBL), a phenomenon that is expected to freeze the dynamics of the system and prevent it from reaching thermal equilibrium. Performing extensive numerical simulations of the imbalance, a quantity often employed in the experimental studies of MBL, we show that the regime of a slow power-law decay of imbalance persists to disorder strengths exceeding by at least a factor of 2 the current estimates of the critical disorder strength for MBL. Even though we investigate time evolution up to few thousands tunneling times, we observe no signs of the saturation of imbalance that would suggest freezing of system dynamics and provide a smoking gun evidence of MBL. We demonstrate that the situation is qualitatively different when the disorder is replaced by a quasiperiodic potential. In this case, we observe an emergence of a pattern of oscillations of the imbalance that is stable with respect to changes in the system size. This suggests that the dynamics of quasiperiodic systems remain fully local at the longest time scales we reach provided that the quasiperiodic potential is sufficiently strong. Our study identifies challenges in an unequivocal experimental observation of the phenomenon of MBL., Comment: Author accepted version to appear in PRB

We analyze the finite-size scaling of the average gap-ratio and the entanglement entropy across the many-body localization (MBL) transition in one dimensional Heisenberg spin-chain with quasi-periodic (QP) potential. By using the recently introduced cost-function approach, we compare different scenarios for the transition using exact diagonalization of systems up to 22 lattice sites. Our findings suggest that the MBL transition in the QP Heisenberg chain belongs to the class of Berezinskii-Kosterlitz-Thouless (BKT) transition, the same as in the case of uniformly disordered systems as advocated in recent studies. Moreover, we observe that the critical disorder strength shows a clear sub-linear drift with the system-size as compared to the linear drift seen in random disordered models, suggesting that the finite-size effects in the MBL transition for the QP systems are less severe than that in the random disordered scenario. Moreover, deep in the ergodic regime, we find an unexpected double-peak structure of distribution of on-site magnetizations that can be traced back to the strong correlations present in the QP potential., Comment: Version accepted for publication in Phys. Rev. B; 11 pages, 9 figures

We investigate the structure of many-body wave functions of 1D quantum circuits with local measurements employing the participation entropies. The leading term in system size dependence of participation entropy indicates a model dependent multifractal scaling of the wave-functions at any non-zero measurement rate. The sub-leading term contains universal information about measurement-induced phase transitions and plays the role of an order parameter, being constant non-zero in the error correcting phase and vanishing in the quantum Zeno phase. We provide robust numerical evidence investigating a variety of quantum many-body systems, and provide an analytical interpretation of this behavior expressing the participation entropy in terms of partition functions of classical statistical models in 2D., Comment: 4+16 pages, added new results on universality and extended discussions, comments welcome

The level statistics in the transition between delocalized and localized {phases of} many body interacting systems is {considered}. We recall the joint probability distribution for eigenvalues resulting from the statistical mechanics for energy level dynamics as introduced by Pechukas and Yukawa. The resulting single parameter analytic distribution is probed numerically {via Monte Carlo method}. The resulting higher order spacing ratios are compared with data coming from different {quantum many body systems}. It is found that this Pechukas-Yukawa distribution compares favorably with {$\beta$--Gaussian ensemble -- a single parameter model of level statistics proposed recently in the context of disordered many-body systems.} {Moreover, the Pechukas-Yukawa distribution is also} only slightly inferior to the two-parameter $\beta$-h ansatz shown {earlier} to reproduce {level statistics of} physical systems remarkably well., Comment: Version accepted for publication in Fritz Haake memorial volume of Journal of Physics A

Quantum systems evolving unitarily and subject to quantum measurements exhibit various types of non-equilibrium phase transitions, arising from the competition between unitary evolution and measurements. Dissipative phase transitions in steady states of time-independent Liouvillians and measurement induced phase transitions at the level of quantum trajectories are two primary examples of such transitions. Investigating a many-body spin system subject to periodic resetting measurements, we argue that many-body dissipative Floquet dynamics provides a natural framework to analyze both types of transitions. We show that a dissipative phase transition between a ferromagnetic ordered phase and a paramagnetic disordered phase emerges for long-range systems as a function of measurement probabilities. A measurement induced transition of the entanglement entropy between volume law scaling and sub-volume law scaling is also present, and is distinct from the ordering transition. The two phases correspond to an error-correcting and a quantum-Zeno regimes, respectively. The ferromagnetic phase is lost for short range interactions, while the volume law phase of the entanglement is enhanced. An analysis of multifractal properties of wave function in Hilbert space provides a common perspective on both types of transitions in the system. Our findings are immediately relevant to trapped ion experiments, for which we detail a blueprint proposal based on currently available platforms., Comment: accepted version, comments welcome

The transition between ergodic and many-body localized phases is expected to occur via an avalanche mechanism, in which \emph{ergodic bubbles} that arise due to local fluctuations in system properties thermalize their surroundings leading to delocalization of the system, unless the disorder is sufficiently strong to stop this process. We propose an algorithm based on neural networks that allows to detect the ergodic bubbles using experimentally measurable two-site correlation functions. Investigating time evolution of the system, we observe a logarithmic in time growth of the ergodic bubbles in the MBL regime. The distribution of the size of ergodic bubbles converges during time evolution to an exponentially decaying distribution in the MBL regime, and a power-law distribution with a thermal peak in the critical regime, supporting thus the scenario of delocalization through the avalanche mechanism. Our algorithm permits to pin-point quantitative differences in time evolution of systems with random and quasiperiodic potentials, as well as to identify rare (Griffiths) events. Our results open new pathways in studies of the mechanisms of thermalization of disordered many-body systems and beyond., Comment: version accepted in PRB Lett

We study the impact of quenched disorder on the dynamics of locally constrained quantum spin chains, that describe 1D arrays of Rydberg atoms in both frozen (Ising-type) and dressed (XY-type) regime. Performing large-scale numerical experiments, we observe no trace of many-body localization even at large disorder. Analyzing the role of quenched disorder terms in constrained systems we show that they act in two, distinct and competing ways: as an on-site disorder term for the basic excitations of the system, and as an interaction between excitations. The two contributions are of the same order, and as they compete (one towards localization, the other against it), one does never enter a truly strong disorder, weak interaction limit, where many-body localization occurs. Such a mechanism is further clarified in the case of XY-type constrained models: there, a term which would represent a bona-fide local quenched disorder term acting on the excitations of the clean model must be written as a series of non-local terms in the unconstrained variables. Our observations provide a simple picture to interpret the role of quenched disorder that could be immediately extended to other constrained models or quenched gauge theories., Comment: 4 + 4 pages, comments welcome

Many-body localization (MBL) behavior is analyzed {in an extended Bose-Hubbard model with quasiperiodic infinite-range interactions. No additional disorder is present. Examining level statistics and entanglement entropy of eigenstates we show that a significant fraction of eigenstates of the system is localized in the presence of strong interactions. In spite of this, our results suggest that the system becomes ergodic in the standard thermodynamic limit in which the energy of the system is extensive. At the same time, the MBL regime seems to be stable if one allows for a super-extensive scaling of the energy. We show that our findings can be experimentally verified by studies of time dynamics in many-body cavity quantum electrodynamics setups. The "quench spectroscopy" is a particularly effective tool that allows us to systematically study energy dependence of time dynamics and to investigate a mobility edge in our system., Comment: not too many interesting pages

The many-body localization transition for Heisenberg spin chain with a speckle disorder is studied. Such a model is equivalent to a system of spinless fermions in an optical lattice with an additional speckle field. Our numerical results show that the many-body localization transition in speckle disorder falls within the same universality class as the transition in an uncorrelated random disorder, in contrast to the quasiperiodic potential typically studied in experiments. This hints at possibilities of experimental studies of the role of rare Griffiths regions and of the interplay of ergodic and localized grains at the many-body localization transition. Moreover, the speckle potential allows one to study the role of correlations in disorder on the transition. We study both spectral and dynamical properties of the system focusing on observables that are sensitive to the disorder type and its correlations. In particular, distributions of local imbalance at long times provide an experimentally available tool that reveals the presence of small ergodic grains even deep in the many-body localized phase in a correlated speckle disorder., Comment: version close to accepted in Phys. Rev. B

Thermalization of random-field Heisenberg spin chain is probed by time evolution of density correlation functions. Studying the impacts of average energies of initial product states on dynamics of the system, we provide arguments in favor of the existence of a mobility edge in the large system-size limit., Comment: 7+3 pages, 7+4 figures. Close to published version

Polynomially filtered exact diagonalization method (POLFED) for large sparse matrices is introduced. The algorithm finds an optimal basis of a subspace spanned by eigenvectors with eigenvalues close to a specified energy target by a spectral transformation using a high order polynomial of the matrix. The memory requirements scale better with system size than in the state-of-the-art shift-invert approach. The potential of POLFED is demonstrated examining many-body localization transition in 1D interacting quantum spin-1/2 chains. We investigate the disorder strength and system size scaling of Thouless time. System size dependence of bipartite entanglement entropy and of the gap ratio highlights the importance of finite-size effects in the system. We discuss possible scenarios regarding the many-body localization transition obtaining estimates for the critical disorder strength., Comment: 4+5 pages, version accepted in Physical Review Letters, comments welcome

We analyze the localization properties of the disordered Hubbard model in the presence of a synthetic magnetic field. An analysis of level spacing ratio shows a clear transition from ergodic to many-body localized phase. The transition shifts to larger disorder strengths with increasing magnetic flux. Study of dynamics of local correlations and entanglement entropy indicates that charge excitations remain localized whereas spin degree of freedom gets delocalized in the presence of the synthetic flux. This residual ergodicity is enhanced by the presence of the magnetic field with dynamical observables suggesting incomplete localization at large disorder strengths. Furthermore, we examine the effect of quantum statistics on the local correlations and show that the long-time spin oscillations of a hard-core boson system are destroyed as opposed to the fermionic case., Comment: delayed version close to the published one

Spectral statistics of disordered systems encode Thouless and Heisenberg time scales whose ratio determines whether the system is chaotic or localized. Identifying similarities between system size and disorder strength scaling of Thouless time for disordered quantum many-body systems with results for 3D and 5D Anderson models, we argue that the two-parameter scaling breaks down in the vicinity of the transition to the localized phase signalling subdiffusive dynamics., Comment: 2nd version, several minor changes in text and discussions expanded, 4+1 pages, 3+1 figures, comments welcome

We compare accuracy of two prime time evolution algorithms involving Matrix Product States - tDMRG (time-dependent density matrix renormalization group) and TDVP (time-dependent variational principle). The latter is supposed to be superior within a limited and fixed auxiliary space dimension. Surprisingly, we find that the performance of algorithms depends on the model considered. In particular, many-body localized systems as well as the crossover regions between localized and delocalized phases are better described by tDMRG, contrary to the delocalized regime where TDVP indeed outperforms tDMRG in terms of accuracy and reliability. As an example, we study many-body localization transition in a large size Heisenberg chain. We discuss drawbacks of previous estimates [Phys. Rev. B 98, 174202 (2018)] of the critical disorder strength for large systems., Comment: 11 pages, 17 figures, accepted version, comments welcome

We numerically study level statistics of disordered interacting quantum many-body systems. A two-parameter plasma model which controls level repulsion exponent $\beta$ and range $h$ of interactions between eigenvalues is shown to reproduce accurately features of level statistics across the transition from ergodic to many-body localized phase. Analysis of higher order spacing ratios indicates that the considered $\beta$-$h$ model accounts even for long range spectral correlations and allows to obtain a clear picture of the flow of level statistics across the transition. Comparing spectral form factors of $\beta$-$h$ model and of a system in the ergodic-MBL crossover, we show that the range of effective interactions between eigenvalues $h$ is related to the Thouless time which marks the onset of quantum chaotic behavior of the system. Analysis of level statistics of random quantum circuit which hosts chaotic and localized phases supports the claim that $\beta$-$h$ model grasps universal features of level statistics in transition between ergodic and many-body localized phases also for systems breaking time-reversal invariance., Comment: normalization in eq. 4 added, 10pp, 12 figs accepted in PRB

The level dynamics across the many body localization transition is examined for XXZ-spin model with a random magnetic field. We compare different scenaria of parameter dependent motion in the system and consider measures such as level velocities, curvatures as well as their fidelity susceptibilities. Studying the ergodic phase of the model we find that the level dynamics does not reveal the commonly believed universal behavior after rescaling the curvatures by the level velocity variance. At the same time, distributions of level curvatures and fidelity susceptibilities coincide with properly rescaled distributions for Gaussian Orthogonal Ensemble of random matrices. Profound differences exists depending on way the level dynamics is imposed in the many-body localized phase of the model in which the level dynamics can be understood with the help of local integrals of motion., Comment: version close to that accepted in PRB

We show that a one-dimensional Hubbard model with all-to-all coupling may exhibit many-body localization in the presence of local disorder. We numerically identify the parameter space where many-body localization occurs using exact diagonalization and finite-size scaling. The time evolution from a random initial state exhibits features consistent with the localization picture. The dynamics can be observed with quantum gases in optical cavities, localization can be revealed through the time-dependent dynamics of the light emitted by the resonator., Comment: 30pp. significantly modified and expanded version

The fidelity susceptibility measures sensitivity of eigenstates to a change of an external parameter. It has been fruitfully used to pin down quantum phase transitions when applied to ground states (with extensions to thermal states). Here we propose to use the fidelity susceptibility as a useful dimensionless measure for complex quantum systems. We find analytically the fidelity susceptibility distributions for Gaussian orthogonal and unitary universality classes for arbitrary system size. The results are verified by a comparison with numerical data., Comment: 2nd version, 5+5pp, comments welcome

Level statistics of systems that undergo many--body localization transition are studied. An analysis of the gap ratio statistics from the perspective of inter- and intra-sample randomness allows us to pin point differences between transitions in random and quasi-random disorder, showing the effects due to Griffiths rare events for the former case. It is argued that the transition in the case of random disorder exhibits universal features that are identified by constructing an appropriate model of intermediate spectral statistics which is a generalization of the family of short-range plasma models. The considered weighted short-range plasma model yields a very good agreement both for level spacing distribution including its exponential tail and the number variance up to tens of level spacings outperforming previously proposed models. In particular, our model grasps the critical level statistics which arise at disorder strength for which the inter-sample fluctuations are the strongest. Going beyond the paradigmatic examples of many-body localization in spin systems, we show that the considered model also grasps the level statistics of disordered Bose- and Fermi-Hubbard models. The remaining deviations for long-range spectral correlations are discussed and attributed mainly to the intricacies of level unfolding., Comment: 19pp. enlarged by including 1807.06983; version accepted in Phys. Rev. B

Spectral statistics of systems that undergo many--body localization transition are studied. An analysis of the gap ratio statistics from the perspective of inter- and intra-sample randomness allows us to pin point differences between transitions in random and quasi-random disorder, showing the effects due to Griffiths rare events for the former case. It is argued that the transition for a random disorder exhibits universal features that are identified by constructing an appropriate model of intermediate spectral statistics which is a generalization of the family of short-range plasma models. The model incorporates the inter- and intra-sample fluctuations and faithfully reproduces level spacing distributions as well as number variance during the transition from ergodic to many--body localized phase. In particular, it grasps the critical level statistics which arise at disorder strength for which the fluctuations are the strongest., Comment: This contribution has been incorporated in the enlarged version of arXiv:1808.02795 without changes in the reported results

Many-body localization for a system of bosons trapped in a one dimensional lattice is discussed. Two models that may be realized for cold atoms in optical lattices are considered. The model with a random on-site potential is compared with previously introduced random interactions model. While the origin and character of the disorder in both systems is different they show interesting similar properties. In particular, many-body localization appears for a sufficiently large disorder as verified by a time evolution of initial density wave states as well as using statistical properties of energy levels for small system sizes. Starting with different initial states, we observe that the localization properties are energy-dependent which reveals an inverted many-body localization edge in both systems (that finding is also verified by statistical analysis of energy spectrum). Moreover, we consider computationally challenging regime of transition between many body localized and extended phases where we observe a characteristic algebraic decay of density correlations which may be attributed to subdiffusion (and Griffiths-like regions) in the studied systems. Ergodicity breaking in the disordered Bose-Hubbard models is compared with the slowing-down of the time evolution of the clean system at large interactions., Comment: expanded second version, comments welcome

We study many-body localization in a one dimensional optical lattice filled with bosons. The interaction between bosons is assumed to be random, which can be realized for atoms close to a microchip exposed to a spatially fluctuating magnetic field. Close to a Feshbach resonance, such controlled fluctuations can be transfered to the interaction strength. We show that the system reveals an inverted mobility edge, with mobile particles at the lower edge of the spectrum. A statistical analysis of level spacings allows us to characterize the transition between localized and excited states. The existence of the mobility edge is confirmed in large systems, by time dependent numerical simulations using tDMRG. A simple analytical model predicts the long time behavior of the system., Comment: 6pp. invited talk at "8th Workshop on Chaos and Localization Phenomena" Warsaw, May 2017

We address the question whether the super-Heisenberg scaling for quantum estimation is realizable. We unify the results of two approaches. In the first one, the original system is compared with its copy rotated by the parameter dependent dynamics. If the parameter is coupled to the one-body part of the Hamiltonian the precision of its estimation is known to scale at most as $N^{-1}$ (Heisenberg scaling) in terms of the number of elementary subsystems used, $N$. The second approach considers fidelity at criticality often leading to an apparent super-Heisenberg scaling. However, scaling of time needed to ensure adiabaticity of the evolution brings back the the Heisenberg limit. We illustrate the general theory on a ferromagnetic Heisenberg spin chain which exhibits such super-Heisenberg scaling of fidelity around the critical value of the magnetic field. Even an elementary estimator represented by a single-site magnetization already outperforms the Heisenberg behavior providing the $N^{-1.5}$ scaling. In this case Fisher information sets the ultimate scaling as $N^{-1.75}$ which can be saturated by measuring magnetization on all sites simultaneously. We discuss universal scaling predictions of the estimation precision offered by such observables, both at zero and finite temperatures, and support them with numerical simulations in the model. We provide an experimental proposal of realization of the considered model via mapping the system to ultra-cold bosons in periodically shaken optical lattice. We explicitly derive that the Heisenberg limit is recovered when time needed for preparation of quantum states involved is taken into acocunt., Comment: Again significantly modified version, 15pages. 5figures,106 references comments WELCOME

The possibility of observing many body localization of ultracold atoms in a one dimensional optical lattice is discussed for random interactions. In the non-interacting limit, such a system reduces to single-particle physics in the absence of disorder, i.e. to extended states. In effect the observed localization is inherently due to interactions and is thus a genuine many-body effect. In the system studied, many-body localization manifests itself in a lack of thermalization visible in temporal propagation of a specially prepared initial state, in transport properties, in the logarithmic growth of entanglement entropy as well as in statistical properties of energy levels., Comment: 5pp, 4figs. version close to published one

We study the ground-state properties of bosons loaded into the $p$-band of a one dimensional optical lattice. We show that the phase diagram of the system is substantially affected by the anharmonicity of the lattice potential. In particular, for a certain range of tunneling strength, the full many-body ground state of the system becomes degenerate. In this region, an additional symmetry of the system, namely the parity of the occupation number of the chosen orbital, is spontaneously broken. The state with nonvanishing staggered angular momentum, which breaks the time-reversal symmetry, becomes the true ground state of the system., Comment: 5 pages, 2 figures