We study how perturbations affect dynamics of integrable many-body quantum systems, causing transition from integrability to chaos. Looking at spin transport in the Heisenberg chain with impurities we find that in the thermodynamic limit transport gets diffusive already at an infinitesimal perturbation. Small extensive perturbations therefore cause an immediate transition from integrability to chaos. Nevertheless, there is a remnant of integrability encoded in the dependence of the diffusion constant on the impurity density, namely, at small densities it is proportional to the square root of the inverse density, instead of to the inverse density as would follow from Matthiessen's rule. We show that Matthiessen's rule has to be modified in nonballistic systems. Results also highlight a nontrivial role of interacting scattering on a single impurity, and that there is a regime where adding more impurities can actually increase transport.
We study how perturbations affect dynamics of integrable many-body quantum systems, causing transition from integrability to chaos. Looking at spin transport in the Heisenberg chain with impurities we find that in the thermodynamic limit transport gets diffusive already at an infinitesimal perturbation. Small extensive perturbations therefore cause an immediate transition from integrability to chaos. Nevertheless, there is a remnant of integrability encoded in the dependence of the diffusion constant on the impurity density, namely, at small densities it is proportional to the square root of the inverse density, instead of to the inverse density as would follow from Matthiessen's rule. We show that Matthiessen's rule has to be modified in nonballistic systems. Results also highlight a nontrivial role of interacting scattering on a single impurity, and that there is a regime where adding more impurities can actually increase transport.
We study the breaking of integrability by a finite density of dilute impurities, specifically the emerging diffusive transport. Provided the distance between impurities (localized perturbations) is large, one would expect the scattering rates to be additive, and therefore, the resistivity to be proportional to the number of impurities (the so-called Matthiessen's rule). We show that this is, in general, not the case. If transport is anomalous in the original integrable system without impurities, the diffusion constant in the nonintegrable system at low impurity density has a nontrivial power-law dependence on the impurity density, with the power being determined by the dynamical scaling exponent of anomalous transport. We also find a regime at high impurity density in which, counterintuitively, adding more impurities to an already diffusive system increases transport rather than decreases it.
We show that the onset of quantum chaos at infinite temperature in two many-body one-dimensional lattice models, the perturbed spin-1/2 XXZ and Anderson models, is characterized by universal behavior. Specifically, we show that the onset of quantum chaos is marked by maxima of the typical fidelity susceptibilities that scale with the square of the inverse average level spacing, saturating their upper bound, and that the strength of the integrability- or localization-breaking perturbation at these maxima decreases with increasing system size. We also show that the spectral function below the ``Thouless'' energy (in the quantum-chaotic regime) diverges when approaching those maxima. Our results suggest that, in the thermodynamic limit, arbitrarily small integrability- or localization-breaking perturbations result in quantum chaos in the many-body quantum systems studied here.
Integrable quantum systems of finite size are generically robust against weak enough integrability-breaking perturbations, but become quantum chaotic and thermalizing if the integrability breaking is strong enough. We argue that the onset of quantum chaos can be described as a Fock-space delocalization process, with the eigenstates of the integrable system being taken as the Fock states. The integrability-breaking perturbation introduces hopping in this Fock space, and chaos sets in when this hopping delocalizes the many-body eigenstates in this space. Depending on the range of the dominant Fock-space hopping, delocalization can occur either through a crossover or via a transition that becomes sharp in the appropriate large-system dynamic limit. In either case, the perturbation strength at the onset of chaos scales to zero in the usual thermodynamic limit, with a size dependence that we estimate analytically and compute numerically for a few specific models. We also identify two intermediate finite-size-dependent regimes: There is generally an intermediate nonchaotic regime in which integrability is broken strongly enough to produce some systemwide many-body resonances but not enough to thermalize the system. In spatially extended systems (but not in quantum dots) there is also a crossover or transition between chaotic regimes where the ratio of the system size to the mean free path of the quasiparticles of the integrable system is small versus large compared to unity.
We study transport in spin chains employing the Thouless approach based on the level sensitivity to the boundary conditions, R. Although spin transport in the integrable easy-axis XXZ model is diffusive, the corresponding R is much closer to ballistic chains than to chaotic diffusive systems. In the case of the grand canonical ensemble, this observation can be rigorously justified, while in the case of the canonical ensemble it can be demonstrated by numerical calculations. Integrability-breaking perturbation (IBP) strongly reduces R, which reveals a pronounced minimum at the crossover from anomalous diffusive to normal dissipative transport. This minimum coincides with the onset of the universality of the random matrix theory. Results for various IBPs suggest a discontinuous jump of the spin conductivity in the thermodynamic limit, and moreover that its value is proportional to the strength of the IBP.
Understanding the physics of the integrable spin-1/2 XXZ chain has witnessed substantial progress, due to the development and application of sophisticated analytical and numerical techniques. In particular, infinite-temperature magnetization transport has turned out to range from ballistic, over superdiffusive, to diffusive behavior in different parameter regimes of the anisotropy. Since integrability is rather the exception than the rule, a crucial question is the change of transport under integrability-breaking perturbations. This question includes the stability of superdiffusion at the isotropic point and the change of diffusion constants in the easy-axis regime. In our work, we study this change of diffusion constants by a variety of methods and cover both, linear response theory in the closed system and the Lindblad equation in the open system, where we throughout focus on periodic boundary conditions. In the closed system, we compare results from the recursion method to calculations for finite systems and find evidence for a continuous change of diffusion constants over the full range of perturbation strengths. In the open system weakly coupled to baths, we find diffusion constants in quantitative agreement with the ones in the closed system in a range of nonweak perturbations, but disagreement in the limit of weak perturbations. Using a simple model in this limit, we point out the possibility of a diverging diffusion constant in such an open system.
We study the ballistic-to-diffusive transition induced by the weak breaking of integrability in a boundary-driven XXZ spin chain. Studying the evolution of the spin current density Js as a function of the system size L, we show that, accounting for boundary effects, the transition has a nontrivial universal behavior close to the XX limit. It is controlled by the scattering length ${L}^{*}\ensuremath{\propto}{V}^{\ensuremath{-}2}$, where V is the strength of the integrability-breaking term. In the XXZ model, the interplay of interactions controls the emergence of a transient ``quasiballistic'' regime at length scales much shorter than L∗. This parametrically large regime is characterized by a strong renormalization of the current which forbids a universal scaling, unlike the XX model. Our results are based on matrix product operator numerical simulations and agree with perturbative analytical calculations.