Local Integrals of Motion in Dipole-Conserving Models with Hilbert Space Fragmentation.

Hilbert space fragmentation is an ergodicity-breaking phenomenon, in which the Hamiltonian shatters into exponentially many dynamically disconnected sectors. In many fragmented systems, these sectors can be labeled by statistically localized integrals of motion, which are nonlocal operators. We study the paradigmatic nearest-neighbor pair hopping model exhibiting the so-called strong fragmentation. We show that this model hosts local integrals of motion (LIOMs), which correspond to frozen density modes with long wavelengths. The latter modes become subdiffusive when longer-range pair hoppings are allowed. Finally, we make a connection with a tilted (Stark) chain. Contrary to the dipole-conserving effective models, the tilted chain is shown to support either a Hamiltonian or dipole moment as an LIOM. Numerical results are obtained from a numerical algorithm, in which finding LIOMs is reduced to a data compression problem.

Generalized Thermalization in Quantum-Chaotic Quadratic Hamiltonians.

Thermalization (generalized thermalization) in nonintegrable (integrable) quantum systems requires two ingredients: equilibration and agreement with the predictions of the Gibbs (generalized Gibbs) ensemble. We prove that observables that exhibit eigenstate thermalization in single-particle sector equilibrate in many-body sectors of quantum-chaotic quadratic models. Remarkably, the same observables do not exhibit eigenstate thermalization in many-body sectors (we establish that there are exponentially many outliers). Hence, the generalized Gibbs ensemble is generally needed to describe their expectation values after equilibration, and it is characterized by Lagrange multipliers that are smooth functions of single-particle energies.

Phenomenology of Spectral Functions in Disordered Spin Chains at Infinite Temperature.

Studies of disordered spin chains have recently experienced a renewed interest, inspired by the question to which extent the exact numerical calculations comply with the existence of a many-body localization phase transition. For the paradigmatic random field Heisenberg spin chains, many intriguing features were observed when the disorder is considerable compared to the spin interaction strength. Here, we introduce a phenomenological theory that may explain some of those features. The theory is based on the proximity to the noninteracting limit, in which the system is an Anderson insulator. Taking the spin imbalance as an exemplary observable, we demonstrate that the proximity to the local integrals of motion of the Anderson insulator determines the dynamics of the observable at infinite temperature. In finite interacting systems our theory quantitatively describes its integrated spectral function for a wide range of disorders.

Quantitative Impact of Integrals of Motion on the Eigenstate Thermalization Hypothesis.

Even though the eigenstate thermalization hypothesis (ETH) may be introduced as an extension of the random matrix theory, physical Hamiltonians and observables differ from random operators. One of the challenges is to embed local integrals of motion (LIOMs) within the ETH. Here we make steps towards a unified treatment of the ETH in integrable and nonintegrable models with translational invariance. Specifically, we focus on the impact of LIOMs on the fluctuations and structure of the diagonal matrix elements of local observables. We first show that nonvanishing fluctuations entail the presence of LIOMs. Then we introduce a generic protocol to construct observables, subtracted by their projections on LIOMs as well as products of LIOMs. The protocol systematically reduces fluctuations and/or the structure of the diagonal matrix elements. We verify our arguments by numerical results for integrable and nonintegrable models.

Spin Subdiffusion in the Disordered Hubbard Chain.

We derive and study an effective spin model that explains the anomalous spin dynamics in the one-dimensional Hubbard model with strong potential disorder. Assuming that charges are localized, we show that spins are delocalized and their subdiffusive transport originates from a singular random distribution of spin exchange interactions. The exponent relevant for the subdiffusion is determined by the Anderson localization length and the density of the electrons. Although the analytical derivations are valid for low particle density, numerical results for the full model reveal a qualitative agreement up to half filling.

Identification of Majorana Modes in Interacting Systems by Local Integrals of Motion.

Recently, there has been substantial progress in methods of identifying local integrals of motion in interacting integrable models or in systems with many-body localization. We show that one of these approaches can be utilized for constructing local, conserved, Majorana fermions in systems with an arbitrary many-body interaction. As a test case, we first investigate a noninteracting Kitaev model and demonstrate that this approach perfectly reproduces the standard results. Then, we discuss how the many-body interactions influence the spatial structure and the lifetime of the Majorana modes. Finally, we determine the regime for which the information stored in the Majorana correlators is also retained for arbitrarily long times at high temperatures. We show that it is included in the regime with topologically protected soft Majorana modes, but in some cases is significantly smaller.

Complete Many-Body Localization in the t-J Model Caused by a Random Magnetic Field.

The many body localization (MBL) of spin-1/2 fermions poses a challenging problem. It is known that the disorder in the charge sector may be insufficient to cause full MBL. Here, we study dynamics of a single hole in one dimensional t-J model subject to a random magnetic field. We show that strong disorder that couples only to the spin sector localizes both spin and charge degrees of freedom. Charge localization is confirmed also for a finite concentration of holes. While we cannot precisely pinpoint the threshold disorder, we conjecture that there are two distinct transitions. Weaker disorder first causes localization in the spin sector. Carriers become localized for somewhat stronger disorder, when the spin localization length is of the order of a single lattice spacing.

Nature of Bosonic Excitations Revealed by High-Energy Charge Carriers.

We address a long-standing problem concerning the origin of bosonic excitations that strongly interact with charge carriers. We show that the time-resolved pump-probe experiments are capable of distinguishing between regular bosonic degrees of freedom, e.g., phonons, and the hard-core bosons, e.g., magnons. The ability of phonon degrees of freedom to absorb essentially an unlimited amount of energy renders relaxation dynamics nearly independent of the absorbed energy or fluence. In contrast, the hard core effects pose limits on the density of energy stored in the bosonic subsystems resulting in a substantial dependence of the relaxation time on the fluence and/or excitation energy. Very similar effects can be observed also in a different setup when the system is driven by multiple pulses.

Identifying local and quasilocal conserved quantities in integrable systems.

We outline a procedure for counting and identifying a complete set of local and quasilocal conserved operators in integrable lattice systems. The method yields a systematic generation of all independent, conserved quasilocal operators related to the time average of local operators with a support on up to M consecutive sites. As an example, we study the anisotropic Heisenberg spin-1/2 chain and show that the number of independent conserved operators grows linearly with M. In addition to the known local operators, there exist novel quasilocal conserved quantities in all the parity sectors. The existence of quasilocal conserved operators is shown also for the isotropic Heisenberg model. Implications for the anomalous relaxation of quenched systems are discussed as well.

Breakdown of the generalized Gibbs ensemble for current-generating quenches.

We establish a relation between two hallmarks of integrable systems: the relaxation towards the generalized Gibbs ensemble (GGE) and the dissipationless charge transport. We show that the former one is possible only if the so-called Mazur bound on the charge stiffness is saturated by local conserved quantities. As an example we show how a non-GGE steady state with a current can be generated in the one-dimensional model of interacting spinless fermions with a flux quench. Moreover, an extended GGE involving the quasilocal conserved quantities can be formulated for this case.

Integrable Mott insulators driven by a finite electric field.

We develop a method for extracting the steady nonequilibrium current from studies of driven isolated systems, applying it to the model of a one-dimensional Mott insulator at high temperatures. While in the nonintegrable model the nonequilibrium conditions can be accounted for by internal heating, the integrability leads to a strongly nonlinear dc response with a vanishingly small dc conductivity in the linear-response regime. The finding is consistent with equilibrium results for the dc limit of the optical conductivity determined in the presence of a weak and decreasing perturbation.

Nonlinear current response of an isolated system of interacting fermions.

Nonlinear real-time response of interacting particles is studied on the example of a one-dimensional tight-binding model of spinless fermions driven by electric field. Using equations of motion and numerical methods we show that for a nonintegrable case at finite temperatures the major effect of nonlinearity can be taken into account within the linear response formalism extended by a renormalization of the kinetic energy due to the Joule heating. On the other hand, integrable systems show on constant driving a different universality with a damped oscillating current whereby the frequency is related but not equal to the Bloch oscillations.

The Fulde-Ferrell-Larkin-Ovchinnikov phase in the presence of pair hopping interaction.

The recent experimental support for the presence of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase in CeCoIn(5) directed attention towards the mechanisms responsible for this type of superconductivity. We investigate the FFLO state in a model where on-site/inter-site pairing coexists with the repulsive pair hopping interaction. The latter interaction is interesting in that it leads to pairing with non-zero momentum of the Cooper pairs even in the absence of the external magnetic field (the so-called η pairing). It turns out that, depending on the strength of the pair hopping interaction, the magnetic field can induce one of two types of the FFLO phase with different spatial modulations of the order parameter. It is argued that the properties of the FFLO phase may give information about the magnitude of the pair hopping interaction. We also show that η pairing and d-wave superconductivity may coexist in the FFLO state. It holds true also for superconductors which, in the absence of magnetic field, are of pure d-wave type.

Inhomogeneity-induced enhancement of the pairing interaction in cuprate superconductors.

Scanning tunneling spectroscopy has recently discovered a positive correlation between the magnitude of the superconducting gap and positions of dopant oxygen atoms in Bi-based cuprates. We propose a microscopic mechanism that could be responsible for this effect. In particular, we demonstrate that the dopant-induced spatial variation of the atomic levels always enhances the superexchange interaction.