Search

Discrete time crystalline phases have attracted significant theoretical and experimental attention in the last few years. Such systems require a seemingly impossible combination of nonadiabatic driving and a finite-entropy long-time state, which, surprisingly, is possible in nonergodic systems. Previous works have often relied on disorder for the required nonergodicity; here, we describe the construction of a discrete time crystal (DTC) phase in nondisordered, nonintegrable Ising-type systems. After discussing the conditions for interacting and periodically driven systems to display such phases in general, we propose a concrete model and then provide approximate analytical arguments and direct numerical evidence that it satisfies the conditions and displays a DTC phase robust to local periodic perturbations., Comment: 5 pages, 3 figures

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization (Stark-MBL) at least up to time which scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a non-monotonic dependence of the dynamical exponent on the electric field. As the field increases towards a critical value proportional to the Hubbard interaction strength, transport slows down, becoming sub-diffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting super-diffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

In delocalized systems, particle number fluctuations, also known as quantum surface roughness, and the mean-square displacement exhibit a temporal power-law growth followed by a saturation to a system-size-dependent value. We use simple scaling arguments to show that these quantities satisfy the Family-Vicsek scaling law and derive a dynamic scaling relation between the dynamical exponents, assuming that the saturation times of both quantities scale similarly with the system size. This relation clarifies the mechanism behind quantum surface roughness growth and suggests that diffusive quantum many-body systems belong to the Edwards-Wilkinson universality class. Moreover, it provides a convenient way to assess quantum transport in cold-atom experiments. We numerically verify our results by studying two non-interacting models and one interacting model having regimes with distinct dynamical exponents., Comment: Published versoion

Quantum many-body scars are atypical, highly nonthermal eigenstates embedded in a sea of thermal eigenstates that have been observed in, for example, kinetically constrained quantum many-body models. These special eigenstates are characterized by a bipartite entanglement entropy that scales as most logarithmically with the subsystem size. We use numerical optimization techniques to investigate if quantum many-body scars of the experimentally relevant PXP model can be well approximated by Gaussian states. Gaussian states are described by a number of parameters that scales quadratically with system size, thereby having a much lower complexity than generic quantum many-body states, for which this number scales exponentially. We find that while quantum many-body scars can typically be well approximated by (symmetrized) Gaussian states, this is not the case for ergodic (thermal) eigenstates. This observation suggests that the non-ergodic part of the PXP Hamiltonian is related to certain quadratic parent Hamiltonians, thereby hinting on the origin of the quantum many-body scars., Comment: 8 pages, 5 figures

We provide bounds on temporal fluctuations around the infinite-time average of out-of-time-ordered and time-ordered correlators of many-body quantum systems without energy gap degeneracies. For physical initial states, our bounds predict the exponential decay of the temporal fluctuations as a function of the system size. We numerically verify this prediction for chaotic and interacting integrable spin-1/2 chains, which satisfy the assumption of our bounds. On the other hand, we show analytically and numerically that for the XX model, which is a noninteracting system with gap degeneracies, the temporal fluctuations decay polynomially with system size for operators that are local in the fermion representation and decrease exponentially in the system size for non-local operators. Our results demonstrate that the decay of the temporal fluctuations of correlators cannot be used as a reliable metric of chaos or lack thereof., Comment: 14 pages, 4 figures

We study quantum transport in a quasiperiodic Aubry-Andr\'e-Harper (AAH) model induced by the coupling of the system to a Markovian heat bath. We find that coupling the heat bath locally does not affect transport in the delocalized and critical phases, while it induces logarithmic transport in the localized phase. Increasing the number of coupled sites at the central region introduces a transient diffusive regime, which crosses over to logarithmic transport in the localized phase and in the delocalized regime to ballistic transport. On the other hand, when the heat bath is coupled to equally spaced sites of the system, we observe a crossover from ballistic and logarithmic transport to diffusion in the delocalized and localized regimes, respectively. We propose a classical master equation, which captures our numerical observations for both coupling configurations on a qualitative level and for some parameters, even on a quantitative level. Using the classical picture, we show that the crossover to diffusion occurs at a time that increases exponentially with the spacing between the coupled sites, and the resulting diffusion constant decreases exponentially with the spacing., Comment: 11 pages, 7 figs

We prove that spin chains symmetric under a combination of mirror and spin-flip symmetries and with a nondegenerate spectrum show finite spin transport at zero total magnetization and infinite temperature. We demonstrate this numerically using two prominent examples: the Stark many-body localization system and the symmetrized many-body localization system. We provide evidence of delocalization at all energy densities and show that the delocalization mechanism is robust to breaking the symmetry. We use our results to construct two localized systems which, when coupled, delocalize each other., Comment: Main text: 4+$\epsilon$ pages, 4 figures. Supplemental material: 7 pages, 3 figures

We study the chaotic properties of a large-spin XXZ chain with onsite disorder and a small number of excitations above the fully polarized state. We show that while the classical limit, which is reached for large spins, is chaotic, enlarging the spin suppresses quantum chaos features. We attribute this suppression to the occurrence of large and slow clusters of onsite excitations, and propose a way to facilitate their fragmentation by introducing additional decay channels. We numerically verify that the introduction of such relaxation channels restores chaoticity for large spins, so that only three excitations are required to achieve strong level repulsion and ergodic eigenstates., Comment: 7 pages, 5 figures

The Rosenzweig-Porter random matrix ensemble serves as a qualitative phenomenological model for the level statistics and fractality of eigenstates across the many-body localization transition in static systems. We propose a unitary (circular) analogue of this ensemble, which similarly captures the phenomenology of many-body localization in periodically driven (Floquet) systems. We define this ensemble as the outcome of a Dyson Brownian motion process. We show numerical evidence that this ensemble shares some key statistical properties with the Rosenzweig-Porter ensemble for both the eigenvalues and the eigenstates., Comment: 7 pages, 3 figures

Using numerically exact methods we study transport in an interacting spin chain which for sufficiently strong spatially constant electric field is expected to experience Stark many-body localization. We show that starting from a generic initial state, a spin-excitation remains localized only up to a finite delocalization time, which depends exponentially on the size of the system and the strength of the electric field. This suggests that bona fide Stark many-body localization occurs only in the thermodynamic limit. We also demonstrate that the transient localization in a finite system and for electric fields stronger than the interaction strength can be well approximated by a Magnus expansion up-to times which grow with the electric field strength., Comment: 7 pages, 4 figures and supplementary

We study the dynamical behavior of the Anderson insulator in the presence of a local noise. We show that the noise induces logarithmically slow energy and entanglement growth, until the system reaches an infinite-temperature state, where both quantities saturate to extensive values. The saturation value of the entanglement entropy approaches the average entanglement entropy over all possible product states. At infinite temperature, we find that a density excitation spreads logarithmically with time, without any signs of asymptotic diffusive behavior. In addition, we provide a theoretical picture which qualitatively reproduces the phenomenology of particle transport., Comment: 7 pages, 4 figures; added references and discussion of the entanglement entropy asymptotic value

We propose a mechanism to suppress heating in periodically driven many-body quantum systems by employing sufficiently long-range interactions and experimentally relevant initial conditions. The mechanism is robust to local perturbations and does \emph{not} rely on disorder or high driving frequencies. Instead, it makes use of an approximate fragmentation of the many-body spectrum of the non-driven system into bands, with band gaps that grow with the system size. We show that when these systems are driven, there is a regime where \emph{decreasing} the driving frequency \emph{decreases} heating and entanglement build-up. This is demonstrated numerically for a prototypical system of spins in one dimension, but the results can be readily generalized to higher dimensions., Comment: 5 pages, 4 figures

Isolated many-body quantum systems quenched far from equilibrium can eventually equilibrate, but it is not yet clear how long they take to do so. To answer this question, we use exact numerical methods and analyze the entire evolution, from perturbation to equilibration, of a paradigmatic disordered many-body quantum system in the chaotic regime. We investigate how the equilibration time depends on the system size and observables. We show that if dynamical manifestations of spectral correlations in the form of the correlation hole ("ramp") are taken into account, the time for equilibration scales exponentially with system size, while if they are neglected, the scaling is better described by a power law with system size, though with an exponent larger than what is expected for diffusive transport., Comment: Journal version; 10 pages, 4 figures

We numerically investigate the minimum number of interacting particles, which is required for the onset of strong chaos in quantum systems on a one-dimensional lattice with short-range and long-range interactions. We consider multiple system sizes which are at least three times larger than the number of particles and find that robust signatures of quantum chaos emerge for as few as 4 particles in the case of short-range interactions and as few as 3 particles for long-range interactions, and without any apparent dependence on the size of the system., Comment: 8 pages, 4 figures, SciPost submission

We analyze the saturation value of the bipartite entanglement and number entropy starting from a random product state deep in the MBL phase. By studying the probability distributions of these entropies we find that the growth of the saturation value of the entanglement entropy stems from a significant reshuffling of the weight in the probability distributions from the bulk to the exponential tails. In contrast, the probability distributions of the saturation value of the number entropy are converged with system size, and exhibit a sharp cut-off for values of the number entropy which correspond to one particle fluctuating across the boundary between the two halves of the system. Our results therefore rule out slow particle transport deep in the MBL phase and confirm that the slow entanglement entropy production stems uniquely from configurational entanglement., Comment: 4 pages, 7 figures

We analyze and discuss convergence properties of a numerically exact algorithm tailored to study the dynamics of interacting two-dimensional lattice systems. The method is based on the application of the time-dependent variational principle in a manifold of binary and quaternary Tree Tensor Network States. The approach is found to be competitive with existing matrix product state approaches. We discuss issues related to the convergence of the method, which could be relevant to a broader set of numerical techniques used for the study of two-dimensional systems., Comment: 11 pages, 10 figures

Using a numerically exact technique we study spin transport and the evolution of spin-density excitation profiles in a disordered spin-chain with long-range interactions, decaying as a power-law, $r^{-\alpha}$ with distance and $\alpha<2$. Our study confirms the prediction of recent theories that the system is delocalized in this parameters regime. Moreover we find that for $\alpha>3/2$ the underlying transport is diffusive with a transient super-diffusive tail, similarly to the situation in clean long-range systems. We generalize the Griffiths picture to long-range systems and show that it captures the essential properties of the exact dynamics., Comment: 5 pages, 5 figures, revised and added entanglement entropy

The disordered XXZ model is a prototype model of the many-body localization transition (MBL). Despite numerous studies of this model, the available numerical evidence of multifractality of its eigenstates is not very conclusive due severe finite size effects. Moreover it is not clear if similarly to the case of single-particle physics, multifractal properties of the many-body eigenstates are related to anomalous transport, which is observed in this model. In this work, using a state-of-the-art, massively parallel, numerically exact method, we study systems of up to 24 spins and show that a large fraction of the delocalized phase flows towards ergodicity in the thermodynamic limit, while a region immediately preceding the MBL transition appears to be multifractal in this limit. We discuss the implication of our finding on the mechanism of subdiffusive transport., Comment: 13 pages, 8 figures

Locality imposes stringent constraints on the spreading of information in nonrelativistic quantum systems, which is reminiscent of a "light-cone," a casual structure arising in their relativistic counterparts. Long-range interactions can potentially soften such constraints, allowing almost instantaneous long jumps of particles, thus defying causality. Since interactions decaying as a power-law with distance, $r^{-\alpha}$, are ubiquitous in nature, it is pertinent to understand what is the fate of causality and information spreading in such systems. Using a numerically exact technique we address these questions by studying the out-of-time-order correlation function of a representative generic system in one-dimension. We show that while the interactions are long-range, their effect on information spreading is asymptotically negligible as long as $\alpha>1$. In this range we find a complex compound behavior, where after a short transient a fully local behavior emerges, yielding asymptotic "light-cones" virtually indistinguishable from "light-cones" in corresponding local models. The long-range nature of the interaction is only expressed in the power-law leaking of information from the "light-cone," with the same exponent as the exponent of the interaction, $\alpha$. Our results directly imply that all previously obtained rigorous bounds on information spreading in long-range interacting systems are not tight, and thus could be improved., Comment: 9 pages, 9 figures (corrects inconsistency in the Introduction, other minor changes)

We numerically study spin transport and nonequilibrium spin-density profiles in a clean one-dimensional spin-chain with long-range interactions, decaying as a power-law,$r^{-\alpha}$ with distance. We find two distinct regimes of transport: for $\alpha<1/2$, spin excitations relax instantaneously in the thermodynamic limit, and for $\alpha>1/2$, spin transport combines both diffusive and superdiffusive features. We show that while for $\alpha>3/2$ the spin diffusion coefficient is finite, transport in the system is never strictly diffusive, contrary to corresponding classical systems., Comment: 8 pages, 10 figures

Local observables in generic periodically driven closed quantum systems are known to relax to values described by periodic infinite temperature ensembles. At the same time, ergodic static systems exhibit anomalous thermalization of local observables and satisfy a modified version of the eigenstate thermalization hypothesis (ETH), when disorder is present. This raises the question, how does the introduction of disorder affect relaxation in periodically driven systems? In this work, we analyze this problem by numerically studying transport and thermalization in an archetypal example. We find that thermalization is anomalous and is accompanied by subdiffusive transport with a disorder dependent dynamical exponent. Distributions of matrix elements of local operators in the eigenbases of a family of effective time-independent Hamiltonians, which describe the stroboscopic dynamics of such systems, show anomalous departures from predictions of ETH signaling that only a modified version of ETH is satisfied. The dynamical exponent is shown to be related to the scaling of the variance of these distributions with system size., Comment: 7 pages, 8 figures, version published in Phys. Rev. B (Rapid Communication)

Recently it has been shown that interparticle interactions\emph ongenerically\emph default destroy dynamical localization in periodically driven systems, resulting in diffusive transport and heating. In this work we rigorously construct a family of interacting driven systems which are dynamically localized and effectively decoupled from the external driving potential. We show that these systems exhibit tunable periodic or quasiperiodic revivals of the many-body wavefunction and thus\emph onof all\emph default physical observables. By numerically examining spinless fermions on a one dimensional lattice we show that the analytically obtained revivals of such systems remain stable for finite systems with open boundary conditions while having a finite lifetime in the presence of static spatial disorder. We find this lifetime to be inversely proportional to the disorder strength., Comment: 4 pages, 4 figures

We study the applicability of the time-dependent variational principle in matrix product state manifolds for the long time description of quantum interacting systems. By studying integrable and nonintegrable systems for which the long time dynamics are known we demonstrate that convergence of long time observables is subtle and needs to be examined carefully. Remarkably, for the disordered nonintegrable system we consider the long time dynamics are in good agreement with the rigorously obtained short time behavior and with previous obtained numerically exact results, suggesting that at least in this case the apparent convergence of this approach is reliable. Our study indicates that while great care must be exercised in establishing the convergence of the method, it may still be asymptotically accurate for a class of disordered nonintegrable quantum systems., Comment: We trade the discussion of a diffusive integrable system in favor of a discussion of diffusive nonintegrable system, which better highlights the outcome of our work

Using a numerically exact method we study the stability of dynamical localization to the addition of interactions in a periodically driven isolated quantum system which conserves only the total number of particles. We find that while even infinitesimally small interactions destroy dynamical localization, for weak interactions density transport is significantly suppressed and is asymptotically diffusive, with a diffusion coefficient proportional to the interaction strength. For systems tuned away from the dynamical localization point, even slightly, transport is dramatically enhanced and within the largest accessible systems sizes a diffusive regime is only pronounced for sufficiently small detunings., Comment: Scipost resubmission. 14 pages, 4 figures. Changes to the figures. Corrects a few typos

Using a combination of numerically exact and renormalization-group techniques we study the nonequilibrium transport of electrons in an one-dimensional interacting system subject to a quasiperiodic potential. For this purpose we calculate the growth of the mean-square displacement as well as the melting of domain walls. While the system is nonintegrable for all studied parameters, there is no on finite region default of parameters for which we observe diffusive transport. In particular, our model shows a rich dynamical behavior crossing over from superdiffusion to subdiffusion. We discuss the implications of our results for the general problem of many-body localization, with a particular emphasis on the rare region Griffiths picture of subdiffusion., Comment: 6 pages, 5 figures. A more detailed analysis of the dynamical exponents extraction and discussion of the relevant times. Adds a log-derivative for the FRG section

Entanglement growth and out-of-time-order correlators (OTOC) are used to assess the propagation of information in isolated quantum systems. In this work, using large scale exact time-evolution we show that for weakly disordered nonintegrable systems information propagates behind a ballistically moving front, and the entanglement entropy growths linearly in time. For stronger disorder the motion of the information front is algebraic and sub-ballistic and is characterized by an exponent which depends on the strength of the disorder, similarly to the sublinear growth of the entanglement entropy. We show that the dynamical exponent associated with the information front coincides with the exponent of the growth of the entanglement entropy for both weak and strong disorder. We also demonstrate that the temporal dependence of the OTOC is characterized by a fast\emph onnonexponential\emph default growth, followed by a slow saturation after the passage of the information front. Finally,we discuss the implications of this behavioral change on the growth of the entanglement entropy., Comment: 6 pages, 5 figures

Recent studies point towards nontriviality of the ergodic phase in systems exhibiting many-body localization (MBL), which shows subexponential relaxation of local observables, subdiffusive transport and sublinear spreading of the entanglement entropy. Here we review the dynamical properties of this phase and the available numerically exact and approximate methods for its study. We discuss in which sense this phase could be considered ergodic and present possible phenomenological explanations of its dynamical properties. We close by analyzing to which extent the proposed explanations were verified by numerical studies and present the open questions in this field., Comment: Review article (25 pages, 2 figures). Adds missing references to the Appendix on the relation between the mean-square displacement and the current-current correlation functions

Many-body localization (MBL) in a one-dimensional Fermi Hubbard model with random on-site interactions is studied. While for this model all single-particle states are trivially delocalized, it is shown that for sufficiently strong disordered interactions the model is many-body localized. It is therefore argued that MBL does not necessary rely on localization of the single-particle spectrum. This model provides a convenient platform to study pure MBL phenomenology, since Anderson localization in this model does not exist. By examining various forms of the interaction term a dramatic effect of symmetries on charge transport is demonstrated. A possible realization in a cold atom experiments is proposed., Comment: 10 pages, 9 figures

It is commonly believed that quantum isolated systems satisfying the eigenstate thermalization hypothesis (ETH) are diffusive. We show that this assumption is too restrictive, since there are systems that are asymptotically in a thermal state, yet exhibit anomalous, subdiffusive thermalization. We show that such systems satisfy a modified version of the ETH ansatz and derive a general connection between the scaling of the variance of the offdiagonal matrix elements of local operators, written in the eigenbasis of the Hamiltonian, and the dynamical exponent. We find that for subdiffusively thermalizing systems the variance scales more slowly with system size than expected for diffusive systems. We corroborate our findings by numerically studying the distribution of the coefficients of the eigenfunctions and the offdiagonal matrix elements of local operators of the random field Heisenberg chain, which has anomalous transport in its thermal phase. Surprisingly, this system also has non-Gaussian distributions of the eigenfunctions, thus directly violating Berry's conjecture., Comment: 5 pages, 3 figures; generalized derivations and introduced analogy with Thouless time

We study the real-time dynamics of a two-dimensional Anderson--Hubbard model using nonequilibrium self-consistent perturbation theory within the second-Born approximation. When compared with exact diagonalization performed on small clusters, we demonstrate that for strong disorder this technique approaches the exact result on all available timescales, while for intermediate disorder, in the vicinity of the many-body localization transition, it produces quantitatively accurate results up to nontrivial times. Our method allows for the treatment of system sizes inaccessible by any numerically exact method and for the complete elimination of finite size effects for the times considered. We show that for a sufficiently strong disorder the system becomes nonergodic, while for intermediate disorder strengths and for all accessible time scales transport in the system is strictly subdiffusive. We argue that these results are incompatible with a simple percolation picture, but are consistent with the heuristic random resistor network model where subdiffusion may be observed for long times until a crossover to diffusion occurs. The prediction of slow finite-time dynamics in a two-dimensional interacting and disordered system can be directly verified in future cold atoms experiments, Comment: Title change and minor changes in the text

We study the infinite temperature dynamics of a prototypical one-dimensional system expected to exhibit many-body localization. Using numerically exact methods, we establish the dynamical phase diagram of this system based on the statistics of its eigenvalues and its dynamical behavior. We show that the nonergodic phase is reentrant as a function of the interaction strength, illustrating that localization can be reinforced by sufficiently strong interactions even at infinite temperature. Surprisingly, within the accessible time range, the ergodic phase shows subdiffusive behavior, suggesting that the diffusion coefficient vanishes throughout much of the phase diagram in the thermodynamic limit. Our findings strongly suggest that Wigner-Dyson statistics of eigenvalue spacings may appear in a class of ergodic but subdiffusive systems., Comment: 5 pages, 3 figures. (updated references and journal reference)

Following the field theoretic approach of Basko et al., Ann. Phys. 321, 1126 (2006), we study in detail the real-time dynamics of a system expected to exhibit many-body localization. In particular, for time scales inaccessible to exact methods, we demonstrate that within the second-Born approximation that the temporal decay of the density-density correlation function is non-exponential and is consistent with a finite value for $t\to\infty$, as expected in a non-ergodic state. This behavior persists over a wide range of disorder and interaction strengths. We discuss the implications of our findings with respect to dynamical phase boundaries based both on exact diagonalization studies and as well as those established by the methods of Ref. 1., Comment: 5 pages, 3 figures