Diffusion and ballistic transport in one-dimensional quantum systems

Coherent Spin-Phonon Scattering in Facilitated Rydberg Lattices

We investigate the dynamics of a one-dimensional spin system with facilitation constraint that can be studied using Rydberg atoms in arrays of optical tweezer traps. The elementary degrees of freedom of the system are domains of Rydberg excitations that expand ballistically through the lattice. Because of mechanical forces, Rydberg excited atoms are coupled to vibrations within their traps. At zero temperature and large trap depth, it is known that virtually excited lattice vibrations only renormalize the timescale of the ballistic propagation. However, when vibrational excitations are initially present-i.e., when the external motion of the atoms is prepared in an excited Fock state, coherent state or thermal state-resonant scattering between spin domain walls and phonons takes place. This coherent and deterministic process, which is free from disorder, leads to a reduction of the power-law exponent characterizing the expansion of spin domains. Furthermore, the spin domain dynamics is sensitive to the coherence properties of the atoms' vibrational state, such as the relative phase of coherently superimposed Fock states. Even for a translationally invariant initial state the latter manifests macroscopically in a phase-sensitive asymmetric expansion.

Integrability Breaking from Backscattering

We analyze the onset of diffusive hydrodynamics in the one-dimensional hard-rod gas subject to stochastic backscattering. While this perturbation breaks integrability and leads to a crossover from ballistic to diffusive transport, it preserves infinitely many conserved quantities corresponding to even moments of the velocity distribution of the gas. In the limit of small noise, we derive the exact expressions for the diffusion and structure factor matrices, and show that they generically have off diagonal components. We find that the particle density structure factor is non-Gaussian and singular near the origin, with a return probability showing logarithmic deviations from diffusion.

Nonlinear Transport by Bethe Bound States

We consider nonlinear ballistic spin transport in the XXZ spin chain and derive an analytical result for the nonlinear Drude weight D^{(3)} at infinite temperatures. In contrast to the linear Drude weight D^{(1)}, we find that the result not only depends on anisotropy but also on the string length of the quasiparticles transporting the spin current. Our result provides further insights into transport by quasiparticles and raises questions about Luttinger liquid universality.

Stability of Superdiffusion in Nearly Integrable Spin Chains

Superdiffusive finite-temperature transport has been recently observed in a variety of integrable systems with non-Abelian global symmetries. Superdiffusion is caused by giant Goldstone-like quasiparticles stabilized by integrability. Here, we argue that these giant quasiparticles remain long-lived and give divergent contributions to the low-frequency conductivity σ(ω), even in systems that are not perfectly integrable. We find, perturbatively, that σ(ω)∼ω^{-1/3} for translation-invariant static perturbations that conserve energy and σ(ω)∼|logω| for noisy perturbations. The (presumable) crossover to regular diffusion appears to lie beyond low-order perturbation theory. By contrast, integrability-breaking perturbations that break the non-Abelian symmetry yield conventional diffusion. Numerical evidence supports the distinction between these two classes of perturbations.

Exponential damping induced by random and realistic perturbations

Given a quantum many-body system and the expectation-value dynamics of some operator, we study how this reference dynamics is altered due to a perturbation of the system's Hamiltonian. Based on projection operator techniques, we unveil that if the perturbation exhibits a random-matrix structure in the eigenbasis of the unperturbed Hamiltonian, then this perturbation effectively leads to an exponential damping of the original dynamics. Employing a combination of dynamical quantum typicality and numerical linked cluster expansions, we demonstrate that our theoretical findings for random matrices can, in some cases, be relevant for the dynamics of realistic quantum many-body models as well. Specifically, we study the decay of current autocorrelation functions in spin-1/2 ladder systems, where the rungs of the ladder are treated as a perturbation to the otherwise uncoupled legs. We find a convincing agreement between the exact dynamics and the lowest-order prediction over a wide range of interchain couplings.

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.

Slow Relaxation and Diffusion in Holographic Quantum Critical Phases

The dissipative dynamics of strongly interacting systems are often characterized by the timescale set by the inverse temperature τ_{P}∼ℏ/(k_{B}T). We show that near a class of strongly interacting quantum critical points that arise in the infrared limit of translationally invariant holographic theories, there is a collective excitation (a quasinormal mode of the dual black hole spacetime) whose lifetime τ_{eq} is parametrically longer than τ_{P}: τ_{eq}≫T^{-1}. The lifetime is enhanced due to its dependence on a dangerously irrelevant coupling that breaks the particle-hole symmetry and the invariance under Lorentz boosts of the quantum critical point. The thermal diffusivity (in units of the butterfly velocity) is anomalously large near the quantum critical point and is governed by τ_{eq} rather than τ_{P}. We conjecture that there exists a long-lived, propagating collective mode with velocity v_{s}, and in this case the relation D=v_{s}^{2}τ_{eq} holds exactly in the limit Tτ_{eq}≫1. While scale invariance is broken, a generalized scaling theory still holds provided that the dependence of observables on the dangerously irrelevant coupling is incorporated. Our work further underlines the connection between dangerously irrelevant deformations and slow equilibration.

Anomalous relaxation and the high-temperature structure factor of XXZ spin chains

We compute the spin-structure factor of XXZ spin chains in the Heisenberg and gapped (Ising) regimes in the high-temperature limit for nonzero magnetization, within the framework of generalized hydrodynamics, including diffusive corrections. The structure factor shows a hierarchy of timescales in the gapped phase, owing to s-spin magnon bound states ("strings") of various sizes. Although short strings move ballistically, long strings move primarily diffusively as a result of their collisions with short strings. The interplay between these effects gives rise to anomalous power-law decay of the spin-structure factor, with continuously varying exponents, at any fixed separation in the late-time limit. We elucidate the cross-over to diffusion (in the gapped phase) and to superdiffusion (at the isotropic point) in the half-filling limit. We verify our results via extensive matrix product operator calculations.

Ballistic Spin Transport in a Periodically Driven Integrable Quantum System

We demonstrate ballistic spin transport of an integrable unitary quantum circuit, which can be understood either as a paradigm of an integrable periodically driven (Floquet) spin chain, or as a Trotterized anisotropic (XXZ) Heisenberg spin-1/2 model. We construct an analytic family of quasilocal conservation laws that break the spin-reversal symmetry and compute a lower bound on the spin Drude weight, which is found to be a fractal function of the anisotropy parameter. Extensive numerical simulations of spin transport suggest that this fractal lower bound is in fact tight.

Kinetic Theory of Spin Diffusion and Superdiffusion in XXZ Spin Chains

We address the nature of spin transport in the integrable XXZ spin chain, focusing on the isotropic Heisenberg limit. We calculate the diffusion constant using a kinetic picture based on generalized hydrodynamics combined with Gaussian fluctuations: we find that it diverges, and show that a self-consistent treatment of this divergence gives superdiffusion, with an effective time-dependent diffusion constant that scales as D(t)∼t^{1/3}. This exponent had previously been observed in large-scale numerical simulations, but had not been theoretically explained. We briefly discuss XXZ models with easy-axis anisotropy Δ>1. Our method gives closed-form expressions for the diffusion constant D in the infinite-temperature limit for all Δ>1. We find that D saturates at large anisotropy, and diverges as the Heisenberg limit is approached, as D∼(Δ-1)^{-1/2}.

Superdiffusion in One-Dimensional Quantum Lattice Models

We identify a class of one-dimensional spin and fermionic lattice models that display diverging spin and charge diffusion constants, including several paradigmatic models of exactly solvable, strongly correlated many-body dynamics such as the isotropic Heisenberg spin chains, the Fermi-Hubbard model, and the t-J model at the integrable point. Using the hydrodynamic transport theory, we derive an analytic lower bound on the spin and charge diffusion constants by calculating the curvature of the corresponding Drude weights at half-filling, and demonstrate that for certain lattice models with isotropic interactions some of the Noether charges exhibit superdiffusive transport at finite temperature and half-filling.

Hydrodynamic Diffusion in Integrable Systems

We show that hydrodynamic diffusion is generically present in many-body, one-dimensional interacting quantum and classical integrable models. We extend the recently developed generalized hydrodynamic (GHD) to include terms of Navier-Stokes type, which leads to positive entropy production and diffusive relaxation mechanisms. These terms provide the subleading diffusive corrections to Euler-scale GHD for the large-scale nonequilibrium dynamics of integrable systems, and arise due to two-body scatterings among quasiparticles. We give exact expressions for the diffusion coefficients. Our results apply to a large class of integrable models, including quantum and classical, Galilean and relativistic field theories, chains, and gases in one dimension, such as the Lieb-Liniger model describing cold atom gases and the Heisenberg quantum spin chain. We provide numerical evaluations in the Heisenberg XXZ spin chain, both for the spin diffusion constant, and for the diffusive effects during the melting of a small domain wall of spins, finding excellent agreement with time-dependent density matrix renormalization group numerical simulations.

Clinical significance of optical coherence tomography-guided angioplasty on treatment selection

The present study aimed to observe whether optical coherence tomography (OCT)-guided angioplasty is able to provide useful clinical information beyond that obtained by angiography as well as provide recommendations for physicians that may improve treatment selection. This prospective study included 83 patients with coronary artery disease (>18 years) undergoing coronary angiography (CAG) for ST-elevation myocardial infarction (n=13), non-ST-elevation myocardial infarction (n=19), stable angina (n=22), unstable angina (n=10), silent ischemia (n=11), or elective percutaneous coronary intervention (n=8). Following the initial CAG (CAG-pre), the patients underwent OCT before angioplasty (OCT-pre, 24 patients), after angioplasty (OCT-post, 22 patients), or both (37 patients). The thrombus burden, calcification and plaque dissection or rupture were compared between the OCT-pre and CAG-pre recordings. Following angioplasty, stent malapposition, suboptimal stent deployment, suboptimal stent lesion coverage, and edge dissection were compared between OCT-post and CAG-post alone. Among the 83 patients, 45.7% had single-vessel and 54.3% had multiple-vessel disease. OCT pre- and post-angioplasty revealed significantly more information on the procedure than CAG alone. This clinical information changed the clinical strategies in 41/83 (49.4%) patients, including 58 modifications of therapeutic strategy (69.9%, 58/83): Thrombus aspiration in 2 cases (2.4%), administration of glycoprotein IIb/IIIa inhibitors in 8 cases (9.6%), additional balloon inflation in 23 cases (27.7%), additional stent implantation in 17 cases (20.5%), avoiding stent implantation in 4 cases (4.8%), collateral intervention in 2 cases (2.4%), and guidewire reposition in 2 cases (2.4%). In conclusion, OCT-pre and OCT-post provided additional clinical information beyond that obtained by angiography alone, which resulted in modification of the treatment strategies in half of the included patients.

Dynamics of coherence, localization and excitation transfer in disordered nanorings

Self-assembled supramolecular aggregates are excellent candidates for the design of efficient excitation transport devices. Both artificially prepared and natural photosynthetic aggregates in plants and bacteria present an important degree of disorder that is supposed to hinder excitation transport. Besides, molecular excitations couple to nuclear motion affecting excitation transport in a variety of ways. We present an exhaustive study of exciton dynamics in disordered nanorings with long-range interactions under the influence of a phonon bath taking the LH2 system of purple bacteria as a model. Nuclear motion is explicitly taken into account by employing the Davydov ansatz description of the polaron and quantum dynamics are obtained using a time-dependent variational method. We reveal an optimal exciton-phonon coupling that suppresses disorder-induced localization and facilitate excitation de-trapping. This excitation transfer enhancement, mediated by environmental phonons, is attributed to energy relaxation toward extended, low-energy excitons provided by the precise LH2 geometry with anti-parallel dipoles and long-range interactions. An analysis of localization and spectral statistics is followed by dynamic measures of coherence and localization, transfer efficiency and superradiance. Linear absorption, 2D photon-echo spectra and diffusion measures of the exciton are examined to monitor the diffusive behavior as a function of the strengths of disorder and exciton-phonon coupling.

A perspective on quantum integrability in many-body-localized and Yang-Baxter systems

Two of the most active areas in quantum many-particle dynamics involve systems with an unusually large number of conservation laws. Many-body-localized systems generalize ideas of Anderson localization by disorder to interacting systems. While localization still exists with interactions and inhibits thermalization, the interactions between conserved quantities lead to some dramatic differences from the Anderson case. Quantum integrable models such as the XXZ spin chain or Bose gas with delta-function interactions also have infinite sets of conservation laws, again leading to modifications of conventional thermalization. A practical way to treat the hydrodynamic evolution from local equilibrium to global equilibrium in such models is discussed. This paper expands upon a presentation at a discussion meeting of the Royal Society on 7 February 2017. The work described was carried out with a number of collaborators, including Jens Bardarson, Vir Bulchandani, Roni Ilan, Christoph Karrasch, Siddharth Parameswaran, Frank Pollmann and Romain Vasseur.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Solvable Hydrodynamics of Quantum Integrable Systems

The conventional theory of hydrodynamics describes the evolution in time of chaotic many-particle systems from local to global equilibrium. In a quantum integrable system, local equilibrium is characterized by a local generalized Gibbs ensemble or equivalently a local distribution of pseudomomenta. We study time evolution from local equilibria in such models by solving a certain kinetic equation, the "Bethe-Boltzmann" equation satisfied by the local pseudomomentum density. Explicit comparison with density matrix renormalization group time evolution of a thermal expansion in the XXZ model shows that hydrodynamical predictions from smooth initial conditions can be remarkably accurate, even for small system sizes. Solutions are also obtained in the Lieb-Liniger model for free expansion into vacuum and collisions between clouds of particles, which model experiments on ultracold one-dimensional Bose gases.

Microscopic Origin of Ideal Conductivity in Integrable Quantum Models

Nonergodic dynamical systems display anomalous transport properties. Prominent examples are integrable quantum systems, whose exceptional properties are diverging dc conductivities. In this Letter, we explain the microscopic origin of ideal conductivity by resorting to the thermodynamic particle content of a system. Using group-theoretic arguments we rigorously resolve the long-standing controversy regarding the nature of spin and charge Drude weights in the absence of chemical potentials. In addition, by employing a hydrodynamic description, we devise an efficient computational method to calculate exact Drude weights from the stationary currents generated in an inhomogeneous quench from bipartitioned initial states. We exemplify the method on the anisotropic Heisenberg model at finite temperatures for the entire range of anisotropies, accessing regimes that are out of reach with other approaches. Quite remarkably, spin Drude weight and asymptotic spin current rates reveal a completely discontinuous (fractal) dependence on the anisotropy parameter.

Solution of the Lindblad equation for spin helix states

Using Lindblad dynamics we study quantum spin systems with dissipative boundary dynamics that generate a stationary nonequilibrium state with a nonvanishing spin current that is locally conserved except at the boundaries. We demonstrate that with suitably chosen boundary target states one can solve the many-body Lindblad equation exactly in any dimension. As solution we obtain pure states at any finite value of the dissipation strength and any system size. They are characterized by a helical stationary magnetization profile and a ballistic spin current which is independent of system size, even when the quantum spin system is not integrable. These results are derived in explicit form for the one-dimensional spin-1/2 Heisenberg chain and its higher-spin generalizations, which include the integrable spin-1 Zamolodchikov-Fateev model and the biquadratic Heisenberg chain.

Thermal Conductivity of the One-Dimensional Fermi-Hubbard Model

We study the thermal conductivity of the one-dimensional Fermi-Hubbard model at a finite temperature using a density matrix renormalization group approach. The integrability of this model gives rise to ballistic thermal transport. We calculate the temperature dependence of the thermal Drude weight at half filling for various interaction strengths. The finite-frequency contributions originating from the fact that the energy current is not a conserved quantity are investigated as well. We report evidence that breaking the integrability through a nearest-neighbor interaction leads to vanishing Drude weights and diffusive energy transport. Moreover, we demonstrate that energy spreads ballistically in local quenches with initially inhomogeneous energy density profiles in the integrable case. We discuss the relevance of our results for thermalization in ultracold quantum-gas experiments and for transport measurements with quasi-one-dimensional materials.

Heat Conductivity of the Heisenberg Spin-1/2 Ladder: From Weak to Strong Breaking of Integrability

We investigate the heat conductivity κ of the Heisenberg spin-1/2 ladder at finite temperature covering the entire range of interchain coupling J(⊥), by using several numerical methods and perturbation theory within the framework of linear response. We unveil that a perturbative prediction κ∝J(⊥)(-2), based on simple golden-rule arguments and valid in the strict limit J(⊥)→0, applies to a remarkably wide range of J(⊥), qualitatively and quantitatively. In the large J(⊥) limit, we show power-law scaling of opposite nature, namely, κ∝J(⊥)(2). Moreover, we demonstrate the weak and strong coupling regimes to be connected by a broad minimum, slightly below the isotropic point at J(⊥)=J(∥). Reducing temperature T, starting from T=∞, this minimum scales as κ∝T(-2) down to T on the order of the exchange coupling constant. These results provide for a comprehensive picture of κ(J(⊥),T) of spin ladders.

Heat Transport in Spin Chains with Weak Spin-Phonon Coupling

The heat transport in a system of S=1/2 large-J Heisenberg spin chains, describing closely Sr(2)CuO(3) and SrCuO(2) cuprates, is studied theoretically at T≪J by considering interactions of the bosonized spin excitations with optical phonons and defects. Treating rigorously the multiboson processes, we derive a microscopic spin-phonon scattering rate that adheres to an intuitive picture of phonons acting as thermally populated defects for the fast spin excitations. The mean-free path of the latter exhibits a distinctive T dependence reflecting a critical nature of spin chains and gives a close description of experiments. By the naturalness criterion of realistically small spin-phonon interaction, our approach stands out from previous considerations that require large coupling constants to explain the data and thus imply a spin-Peierls transition, absent in real materials.

Probing ultrafast excitation energy transfer of the chlorosome with exciton–phonon variational dynamics

Excitation energy transfer of the chlorosome is investigated using exciton–phonon variational dynamics revealing ultrafast energy relaxation and exciton delocalization on a 100 fs scale.

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.

Far-from-equilibrium spin transport in Heisenberg quantum magnets

We study experimentally the far-from-equilibrium dynamics in ferromagnetic Heisenberg quantum magnets realized with ultracold atoms in an optical lattice. After controlled imprinting of a spin spiral pattern with an adjustable wave vector, we measure the decay of the initial spin correlations through single-site resolved detection. On the experimentally accessible time scale of several exchange times, we find a profound dependence of the decay rate on the wave vector. In one-dimensional systems, we observe diffusionlike spin transport with a dimensionless diffusion coefficient of 0.22(1). We show how this behavior emerges from the microscopic properties of the closed quantum system. In contrast to the one-dimensional case, our transport measurements for two-dimensional Heisenberg systems indicate anomalous superdiffusion.

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.

Spin-current autocorrelations from single pure-state propagation

We demonstrate that the concept of quantum typicality allows for significant progress in the study of real-time spin dynamics and transport in quantum magnets. To this end, we present a numerical analysis of the spin-current autocorrelation function of the antiferromagnetic and anisotropic spin-1/2 Heisenberg chain as inferred from propagating only a single pure state, randomly chosen as a "typical" representative of the statistical ensemble. Comparing with existing time-dependent density-matrix renormalization group data, we show that typicality is fulfilled extremely well, consistent with an error of our approach, which is perfectly under control and vanishes in the thermodynamic limit. In the long-time limit, our results provide for a new benchmark for the enigmatic spin Drude weight, which we obtain from chains as long as L=33 sites, i.e., from Hilbert spaces of dimensions almost O(104) larger than in existing exact-diagonalization studies.

Lower bounds on high-temperature diffusion constants from quadratically extensive almost-conserved operators

We prove a general theorem which provides a strict lower bound on high-temperature Green-Kubo diffusion constants in locally interacting quantum lattice systems, under the assumption of existence of a quadratically extensive almost conserved quantity: an operator whose commutator with the lattice Hamiltonian is localized on the boundary sites only. We explicitly demonstrate and compute such a bound in two important models in one dimension: in the (isotropic) Heisenberg spin 1/2 chain and in the fermionic Hubbard chain.

NMR signature of one-dimensional behavior of 3He in nanopores

We have performed thermodynamic and NMR relaxation time measurements of 3He adsorbed in the pores of the mesoporous molecular sieve MCM-41 at temperatures down to 1.7 K and at a range of frequencies up to 240 kHz. The MCM-41 substrate comprises a uniform array of quasi-1D straight pores with a diameter of 2.3 nm. We preplated the pores with a monolayer of 4He to achieve an effective diameter of 1.6 nm at low temperatures. We made NMR measurements as a function of line density and frequency to investigate the spin dynamics and the effect of dimensionality. We observed T(1) is proportional to ω1/2, which is characteristic of one-dimensional diffusion. At these temperatures this arises from a classical size effect in the narrow pores. Our results demonstrate the possibility to study the spin dynamics of a 1D Tomonaga-Luttinger liquid at lower temperatures, where the 3He liquid will constitute a quantum 1D system.

Ballistic spreading of entanglement in a diffusive nonintegrable system

We study the time evolution of the entanglement entropy of a one-dimensional nonintegrable spin chain, starting from random nonentangled initial pure states. We use exact diagonalization of a nonintegrable quantum Ising chain with transverse and longitudinal fields to obtain the exact quantum dynamics. We show that the entanglement entropy increases linearly with time before finite-size saturation begins, demonstrating a ballistic spreading of the entanglement, while the energy transport in the same system is diffusive. Thus, we explicitly demonstrate that the spreading of entanglement is much faster than the energy diffusion in this nonintegrable system.

Low temperature ballistic spin transport in the S = 1/2 antiferromagnetic Heisenberg chain compound SrCuO2

We report zero and longitudinal magnetic field muon spin relaxation (μSR) measurements of the spin S = 1/2 antiferromagnetic Heisenberg chain material SrCuO2. We find that in a weak applied magnetic field B0 the spin-lattice relaxation rate λ follows a power law λ is proportional to B(0)(-n) with n = 0.9(3). This result is temperature independent for 5 K ≤ T ≤ 300 K. Within conformal field theory and using the Müller ansatz we conclude ballistic spin transport in SrCuO2.

Families of quasilocal conservation laws and quantum spin transport

For fundamental integrable quantum chains with deformed symmetries we outline a general procedure for defining a continuous family of quasilocal operators whose time derivative is supported near the two boundary sites only. The program is implemented for a spin 1/2 XXZ chain, resulting in improved rigorous estimates for the high temperature spin Drude weight.

Matrix-valued Boltzmann equation for the nonintegrable Hubbard chain

The standard Fermi-Hubbard chain becomes nonintegrable by adding to the nearest neighbor hopping additional longer range hopping amplitudes. We assume that the quartic interaction is weak and investigate numerically the dynamics of the chain on the level of the Boltzmann type kinetic equation. Only the spatially homogeneous case is considered. We observe that the huge degeneracy of stationary states in the case of nearest neighbor hopping is lost and the convergence to the thermal Fermi-Dirac distribution is restored. The convergence to equilibrium is exponentially fast. However for small next-nearest neighbor hopping amplitudes one has a rapid relaxation towards the manifold of quasistationary states and slow relaxation to the final equilibrium state.

Eigenvalue statistics of reduced density matrix during driving and relaxation

We study a subsystem of an isolated one-dimensional correlated metal when it is driven by a steady electric field or when it relaxes after driving. We obtain numerically exact reduced density matrix ρ for subsystems which are sufficiently large to give significant eigenvalue statistics and spectra of log(ρ). We show that both for generic as well as for the integrable model, the statistics follows the universality of Gaussian unitary and orthogonal ensembles for driven and equilibrium systems, respectively. Moreover, the spectra of modestly driven subsystems are well described by the Gibbs thermal distribution with the entropy determined by the time-dependent energy only.

Burnett coefficients in quantum many-body systems

The Burnett coefficient B is investigated for transport in one-dimensional quantum many-body systems. Extensive numerical computations in spin-1/2 chains suggest a linear growth with time, B(t)~t, for nonintegrable chains exhibiting diffusive transport. For integrable spin chains in the metallic regime, on the other hand, we find a cubic growth with time, B(t)~-D(m)(2)t(3), with the proportionality constant being simply a square of the Drude weight D(m). The results are corroborated with additional studies in noninteracting quantum chains and in the classical limit of large-spin chains.

Conductivity in the Heisenberg chain with next-to-nearest-neighbor interaction

We consider a spin chain given by the XXZ model with a weak next-to-nearest-neighbor perturbation that breaks its exact integrability. We prove that such a system has an ideal metallic behavior (infinite conductivity), by rigorously establishing strict lower bounds on the zero-temperature Drude weight, which are strictly positive. The proof is based on exact renormalization group methods allowing us to prove the convergence of the expansions and to fully take into account the irrelevant terms, which play an essential role in ensuring the correct lattice symmetries. We also prove that the Drude weight verifies the same parameter-free relations as in the absence of the integrability-breaking perturbation.

Optimal transport and von Neumann entropy in a Heisenberg XXZ chain out of equilibrium

In this paper we investigate the spin currents and the von Neumann entropy (vNE) of a Heisenberg XXZ chain in contact with twisted XY-boundary magnetic reservoirs by means of the Lindblad master equation. Exact solutions for the stationary reduced density matrix are explicitly constructed for chains of small sizes by using a quantum symmetry operation of the system. These solutions are then used to investigate the optimal transport in the chain in terms of the vNE. As a result we show that the maximal spin current always occurs in the proximity of minima of the vNE and for particular choices of parameters (coupling with reservoirs and anisotropy) it can exactly coincide with them. As the coupling is increased, current reversals may occur and in the limit of strong coupling we show that minima of the vNE tend to zero, meaning that the maximal transport is achieved in this case with states that are very close to pure states.

Finite-temperature dynamical density matrix renormalization group and the Drude weight of spin-1/2 chains

We propose an easily implemented approach to study time-dependent correlation functions of one-dimensional systems at finite-temperature T using the density matrix renormalization group. The entanglement growth inherent to any time-dependent calculation is significantly reduced if the auxiliary degrees of freedom which purify the statistical operator are time evolved with the physical Hamiltonian but reversed time. We exploit this to investigate the long-time behavior of current correlation functions of the XXZ spin-1/2 Heisenberg chain. This allows a direct extraction of the Drude weight D at intermediate to large T. We find that D is nonzero--and thus transport is dissipationless--everywhere in the gapless phase. At low temperatures we establish an upper bound to D by comparing with bosonization.

Boltzmann-type approach to transport in weakly interacting one-dimensional fermionic systems

We investigate transport properties of one-dimensional fermionic tight binding models featuring nearest and next-nearest neighbor hopping, where the fermions are additionally subject to a weak short range mutual interaction. To this end we employ a pertinent approach which allows for a mapping of the underlying Schrödinger dynamics onto an adequate linear quantum Boltzmann equation. This approach is based on a suitable projection operator method. From this Boltzmann equation we are able to numerically obtain diffusion coefficients in the case of nonvanishing next-nearest neighbor hopping, i.e., the nonintegrable case, whereas the diffusion coefficient diverges without next-nearest neighbor hopping. For the latter case we analytically investigate the decay behavior of the current with the result that arbitrarily small parts of the current relax arbitrarily slowly which suggests anomalous diffusive transport behavior within the scope of our approach.

Spin transport in the XXZ chain at finite temperature and momentum

We investigate the role of momentum for the transport of magnetization in the spin-1/2 Heisenberg chain above the isotropic point at finite temperature and momentum. Using numerical and analytical approaches, we analyze the autocorrelations of density and current and observe a finite region of the Brillouin zone with diffusive dynamics below a cutoff momentum, and a diffusion constant independent of momentum and time, which scales inversely with anisotropy. Lowering the temperature over a wide range, starting from infinity, the diffusion constant is found to increase strongly while the cutoff momentum for diffusion decreases. Above the cutoff momentum diffusion breaks down completely.

Exact nonequilibrium steady state of a strongly driven open XXZ chain

An exact and explicit ladder-tensor-network ansatz is presented for the nonequilibrium steady state of an anisotropic Heisenberg XXZ spin-1/2 chain which is driven far from equilibrium with a pair of Lindblad operators acting on the edges of the chain only. We show that the steady-state density operator of a finite system of size n is-apart from a normalization constant-a polynomial of degree 2n - 2 in the coupling constant. Efficient computation of physical observables is facilitated in terms of a transfer operator reminiscent of a classical Markov process. In the isotropic case we find cosine spin profiles, 1/n(2) scaling of the spin current, and long-range correlations in the steady state. This is a fully nonperturbative extension of a recent result [Phys. Rev. Lett. 106, 217206 (2011)].

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.

Decay of currents for strong interactions

The decay of current autocorrelation functions is investigated for quantum systems featuring strong interactions. Here the term "interaction" refers to that part of the Hamiltonian causing the (major) decay of the current. On the time scale before the (first) zero crossing of the current, its relaxation is shown to be well described by a suitable perturbation theory in the lowest orders of the interaction strength, particularly if interactions are strong. In this description the relaxation is found to be rather close to a Gaussian decay and the resulting diffusion coefficient approximately scales with the inverse interaction strength. These findings are confirmed by numerical results from exact diagonalization for several one-dimensional transport models including spin transport in the Heisenberg chain with respect to different spin quantum numbers, anisotropy, next-nearest-neighbor interaction, or alternating magnetic field; energy transport in the Ising chain with tilted magnetic field; and transport of excitations in a randomly coupled modular quantum system. The impact of these results for weak interactions is discussed.

Spin transport in a one-dimensional anisotropic Heisenberg model

We analytically and numerically study spin transport in a one-dimensional Heisenberg model in linear-response regime at infinite temperature. It is shown that as the anisotropy parameter Δ is varied spin transport changes from ballistic for Δ<1 to anomalous at the isotropic point Δ=1, to diffusive for finite Δ>1, ending up as a perfect isolator in the Ising limit of infinite Δ. Using perturbation theory for large Δ a quantitative prediction is made for the dependence of diffusion constant on Δ.

Open XXZ spin chain: nonequilibrium steady state and a strict bound on ballistic transport

An explicit matrix product ansatz is presented, in the first two orders in the (weak) coupling parameter, for the nonequilibrium steady state of the homogeneous, nearest neighbor Heisenberg XXZ spin 1/2 chain driven by Lindblad operators which act only at the edges of the chain. The first order of the density operator becomes, in the thermodynamic limit, an exact pseudolocal conservation law and yields-via the Mazur inequality-a rigorous lower bound on the high-temperature spin Drude weight. Such a Mazur bound is a nonvanishing fractal function of the anisotropy parameter Δ for |Δ|<1.

Coherent spin-current oscillations in transverse magnetic fields

We address the coherence of the dynamics of spin-currents with components transverse to an external magnetic field for the spin-1/2 Heisenberg chain. We study current autocorrelations at finite temperatures and the real-time dynamics of currents at zero temperature. Besides a coherent Larmor oscillation, we find an additional collective oscillation at higher frequencies, emerging as a coherent many-magnon effect at low temperatures. Using numerical and analytical methods, we analyze the oscillation frequency and decay time of this coherent current-mode versus temperature and magnetic field.

Distribution of NMR relaxations in a random Heisenberg chain

NMR measurements of the (29)Si spin-lattice relaxation time T(1) were used to probe the spin-1/2 random Heisenberg chain compound BaCu(2)(Si(1-x)Ge(x))(2)O(7). Remarkable differences between the pure (x=0) and the fully random (x=0.5) cases are observed, indicating that randomness generates a distribution of local magnetic relaxations. This distribution, which is reflected in a stretched exponential NMR relaxation, exhibits a progressive broadening with decreasing temperature, caused by a growing inequivalence of magnetic sites. Compelling independent evidence for the influence of randomness is also obtained from magnetization data and Monte Carlo calculations. These results suggest the formation of random-singlet states in this class of materials, as previously predicted by theory.

Correlation-induced memory effects in transport properties of low-dimensional systems

We demonstrate the remnant presence of initial correlations in the steady-state electrical current jS flowing between low-dimensional interacting leads. The leads are described as Luttinger liquids and electrons can tunnel via a quantum point contact. We derive an analytic result for the time-dependent current and show that ground-state correlations have a large impact on the relaxation and long-time behavior. In particular, the I-V characteristic is not reproduced by quenching the interaction in time. We further present a universal formula of jS for an arbitrary sequence of interaction quenches and it is established that jS is history dependent for nonsmooth switching process.

Projection operator approach to spin diffusion in the anisotropic Heisenberg chain at high temperatures

Spin transport in the anisotropic Heisenberg chain is typically investigated theoretically with respect to the finiteness of transport coefficients only. Assuming their finiteness at high temperatures, we develop a concrete quantitative picture of the diffusion constant/(dc-)conductivity as a function of both the anisotropy parameter Δ and the spin quantum number s, going beyond the most commonly considered case s=1/2. Using this picture, we enable the comparison of finite transport coefficients from complementary theoretical methods on a quantitative level, having more significance than the finiteness alone. Our method is essentially based on an application of the time-convolutionless projection operator technique to current autocorrelations. This technique, although being a perturbation theory in Δ, is found to be applicable, even if Δ is not small. This finding supports the applicability to a wider class of strongly interacting many-particle quantum systems.

Physics. STM ready for the time domain

The relaxation of a spin on a surface was imaged with a scanning tunneling microscope at 50-ns intervals.