Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1211.2236 · Physical Review B

Nonequilibrium thermal transport and its relation to linear response

Christoph Karrasch, Roni Ilan, Joel E. Moore

Published Nov 15, 2013• 101 citations
arXivHTMLPDF

Abstract

We study the real-time dynamics of spin chains driven out of thermal equilibrium by an initial temperature gradient ${T}_{L}\ensuremath{\ne}{T}_{R}$ using density matrix renormalization group methods. We demonstrate that the nonequilibrium energy current saturates fast to a finite value if the linear-response thermal conductivity is infinite, i.e., if the Drude weight $D$ is nonzero. Our data suggest that a nonintegrable dimerized chain might support such dissipationless transport ($D>0$). We show that the steady-state value ${J}_{E}$ of the current for arbitrary ${T}_{L}\ensuremath{\ne}{T}_{R}$ is of the functional form ${J}_{E}=f({T}_{L})\ensuremath{-}f({T}_{R})$, i.e., it is completely determined by the linear conductance. We argue for this functional form, which is essentially a Stefan-Boltzmann law in this integrable model; for the XXX ferromagnet, $f$ can be computed via the thermodynamic Bethe ansatz in good agreement with the numerics. Inhomogeneous systems exhibiting different bulk parameters as well as Luttinger liquid boundary physics induced by single impurities are discussed briefly.

DOI
10.1103/physrevb.88.195129

Related papers

Nonequilibrium thermal transport and vacuum expansion in the Hubbard model

Christoph Karrasch

One of the most straightforward ways to study thermal properties beyond linear response is to monitor the relaxation of an arbitrarily large left-right temperature gradient ${T}_{L}\ensuremath{-}{T}_{R}$. In one-dimensional systems which support ballistic thermal transport, the local energy currents $\ensuremath{\langle}j(t)\ensuremath{\rangle}$ acquire a nonzero value at long times, and it was recently investigated whether or not this steady state fulfills a simple additive relation $\ensuremath{\langle}j(t\ensuremath{\rightarrow}\ensuremath{\infty})\ensuremath{\rangle}=f({T}_{L})\ensuremath{-}f({T}_{R})$ in integrable models. In this paper, we probe the nonequilibrium dynamics of the Hubbard chain using density matrix renormalization group (DMRG) numerics. We show that the above form provides an effective description of thermal transport in this model; violations are below the finite-time accuracy of the DMRG. As a second setup, we study how an initially equilibrated system radiates into different nonthermal states (such as the vacuum).

arXiv:1612.04964•Physical review. B./Physical review. B•Mar 27, 2017•11 citations

Numerical method for nonlinear steady-state transport in one-dimensional correlated conductors

Malcolm Einhellinger, Alex Cojuhovschi, Eric Jeckelmann

We present a method for investigating the steady-state transport properties of one-dimensional correlated quantum systems. Using a procedure based on our analysis of finite-size effects in a related classical model ($LC$ line) we show that stationary currents can be obtained from transient currents in finite systems driven out of equilibrium. The nonequilibrium dynamics of correlated quantum systems is calculated using the time-evolving block decimation method. To demonstrate our method we determine the full $I$-$V$ characteristic of the spinless fermion model with nearest-neighbor hopping ${t}_{\text{H}}$ and interaction ${V}_{\text{H}}$ using two different setups to generate currents (turning on/off a potential bias). Our numerical results agree with exact results for noninteracting fermions (${V}_{\text{H}}=0$). For interacting fermions we find that in the linear regime $eV\ensuremath{\ll}4{t}_{\text{H}}$ the current $I$ is independent from the setup and our numerical data agree with the predictions of the Luttinger liquid theory combined with the Bethe Ansatz solution. For larger potentials $V$ the steady-state current depends on the current-generating setup and as $V$ increases we find a negative differential conductance with one setup while the currents saturate at finite values in the other one. Both effects are due to finite renormalized bandwidths.

arXiv:1201.5323•Physical Review B•Jun 21, 2012•24 citations

Zero-frequency transport properties of one-dimensional spin- 12 systems

Fabian Heidrich‐Meisner, A. Honecker, D. C. Cabra +1 more

We report a detailed analysis of the Drude weights for both thermal and spin transport in one-dimensional spin-$\frac{1}{2}$ systems by means of exact diagonalization and analytic approaches at finite temperatures. Transport properties are studied first for the integrable $\mathrm{XXZ}$ model and second for various nonintegrable systems such as the dimerized chain, the frustrated chain, and the spin ladder. We compare our results obtained by exact diagonalization and mean-field theory with those of the Bethe ansatz, bosonization, and other numerical studies in the case of the anisotropic Heisenberg model both in the gapless and gapped regime. In particular, we find indications that the Drude weight for spin transport is finite in the thermodynamic limit for the isotropic chain. For the nonintegrable models, a finite-size analysis of the numerical data for the Drude weights is presented, covering the entire parameter space of the dimerized and frustrated chain. We also discuss which conclusions can be drawn from bosonization regarding the question of whether the Drude weights are finite or not. One of our main results is that the Drude weights vanish in the thermodynamic limit for nonintegrable models.

arXiv:cond-mat/0304595•Physical review. B, Condensed matter•Oct 22, 2003•198 citations

Spin and energy currents in integrable and nonintegrable spin- 1 2 chains: A typicality approach to real-time autocorrelations

Robin Steinigeweg, Jochen Gemmer B, Wolfram Brenig

We use the concept of typicality to study the real-time dynamics of spin and energy currents in spin-$1/2$ models in one dimension and at nonzero temperatures. These chains are the integrable XXZ chain and a nonintegrable modification due to the presence of a staggered magnetic field oriented in the $z$ direction. In the framework of linear response theory, we numerically calculate autocorrelation functions by propagating a single pure state, drawn at random as a typical representative of the full statistical ensemble. By comparing to small-system data from exact diagonalization (ED) and existing short-time data from time-dependent density matrix renormalization group, we show that typicality is satisfied in finite systems over a wide range of temperature and is fulfilled in both integrable and nonintegrable systems. For the integrable case, we calculate the long-time dynamics of the spin current and extract the spin Drude weight for large systems outside the range of ED. We particularly provide strong evidence that the high-temperature Drude weight vanishes at the isotropic point. For the nonintegrable case, we obtain the full relaxation curve of the energy current and determine the heat conductivity as a function of magnetic field, exchange anisotropy, and temperature.

arXiv:1408.6837•Physical Review B•Mar 9, 2015•64 citations

Hubbard-to-Heisenberg crossover (and efficient computation) of Drude weights at low temperatures

Christoph Karrasch

We illustrate how finite-temperature charge and thermal Drude weights of one-dimensional systems can be obtained from the relaxation of initial states featuring global (left–right) gradients in the chemical potential or temperature. The approach is tested for spinless interacting fermions as well as for the Fermi-Hubbard model, and the behavior in the vicinity of special points (such as half filling or isotropic chains) is discussed. We present technical details on how to implement the calculation in practice using the density matrix renormalization group and show that the non-equilibrium dynamics is often less demanding to simulate numerically and features simpler finite-time transients than the corresponding linear response current correlators; thus, new parameter regimes can become accessible. As an application, we determine the thermal Drude weight of the Hubbard model for temperatures T which are an order of magnitude smaller than those reached in the equilibrium approach. This allows us to demonstrate that at low T and half filling, thermal transport is successively governed by spin excitations and described quantitatively by the Bethe ansatz Drude weight of the Heisenberg chain.

arXiv:1611.00573•New Journal of Physics•Feb 27, 2017•29 citations

Luttinger liquid physics from the infinite-system density matrix renormalization group

Christoph Karrasch, Joel E. Moore

We study one-dimensional spinless fermions at zero and finite temperature $T$ using the density-matrix renormalization group. We consider nearest- as well as next-nearest-neighbor interactions; the latter render the system inaccessible by a Bethe ansatz treatment. Using an infinite-system algorithm we demonstrate the emergence of Luttinger liquid physics at low energies for a variety of static correlation functions as well as for thermodynamic properties. The characteristic power-law suppression of the momentum distribution $n(k)$ function at $T=0$ can be directly observed over several orders of magnitude. At finite temperature, we show that $n(k)$ obeys a scaling relation. The Luttinger liquid parameter and the renormalized Fermi velocity can be extracted from the density response function, the specific heat, or the susceptibility without the need to carry out any finite-size analysis. We illustrate that the energy scale below which Luttinger liquid power laws manifest vanishes as the half-filled system is driven into a gapped phase by large interactions.

arXiv:1207.0011•Physical Review B•Oct 31, 2012•33 citations
We study the real-time dynamics of spin chains driven out of thermal equilibrium by an initial temperature gradient TL\ensuremath≠TR{T}_{L}\ensuremath{\ne}{T}_{R}TL​\ensuremath=TR​ using density matrix renormalization group methods. We demonstrate that the nonequilibrium energy current saturates fast to a finite value if the linear-response thermal conductivity is infinite, i.e., if the Drude weight DDD is nonzero. Our data suggest that a nonintegrable dimerized chain might support such dissipationless transport (D>0D>0D>0). We show that the steady-state value JE{J}_{E}JE​ of the current for arbitrary TL\ensuremath≠TR{T}_{L}\ensuremath{\ne}{T}_{R}TL​\ensuremath=TR​ is of the functional form JE=f(TL)\ensuremath−f(TR){J}_{E}=f({T}_{L})\ensuremath{-}f({T}_{R})JE​=f(TL​)\ensuremath−f(TR​), i.e., it is completely determined by the linear conductance. We argue for this functional form, which is essentially a Stefan-Boltzmann law in this integrable model; for the XXX ferromagnet, fff can be computed via the thermodynamic Bethe ansatz in good agreement with the numerics. Inhomogeneous systems exhibiting different bulk parameters as well as Luttinger liquid boundary physics induced by single impurities are discussed briefly.
One of the most straightforward ways to study thermal properties beyond linear response is to monitor the relaxation of an arbitrarily large left-right temperature gradient ${T}_{L}\ensuremath{-}{T}_{R}$. In one-dimensional systems which support ballistic thermal transport, the local energy currents $\ensuremath{\langle}j(t)\ensuremath{\rangle}$ acquire a nonzero value at long times, and it was recently investigated whether or not this steady state fulfills a simple additive relation $\ensuremath{\langle}j(t\ensuremath{\rightarrow}\ensuremath{\infty})\ensuremath{\rangle}=f({T}_{L})\ensuremath{-}f({T}_{R})$ in integrable models. In this paper, we probe the nonequilibrium dynamics of the Hubbard chain using density matrix renormalization group (DMRG) numerics. We show that the above form provides an effective description of thermal transport in this model; violations are below the finite-time accuracy of the DMRG. As a second setup, we study how an initially equilibrated system radiates into different nonthermal states (such as the vacuum).
We present a method for investigating the steady-state transport properties of one-dimensional correlated quantum systems. Using a procedure based on our analysis of finite-size effects in a related classical model (LCLCLC line) we show that stationary currents can be obtained from transient currents in finite systems driven out of equilibrium. The nonequilibrium dynamics of correlated quantum systems is calculated using the time-evolving block decimation method. To demonstrate our method we determine the full III-VVV characteristic of the spinless fermion model with nearest-neighbor hopping tH{t}_{\text{H}}tH​ and interaction VH{V}_{\text{H}}VH​ using two different setups to generate currents (turning on/off a potential bias). Our numerical results agree with exact results for noninteracting fermions (VH=0{V}_{\text{H}}=0VH​=0). For interacting fermions we find that in the linear regime $eV\ensuremath{\ll}4{t}_{\text{H}}$ the current III is independent from the setup and our numerical data agree with the predictions of the Luttinger liquid theory combined with the Bethe Ansatz solution. For larger potentials VVV the steady-state current depends on the current-generating setup and as VVV increases we find a negative differential conductance with one setup while the currents saturate at finite values in the other one. Both effects are due to finite renormalized bandwidths.
We report a detailed analysis of the Drude weights for both thermal and spin transport in one-dimensional spin-12\frac{1}{2}21​ systems by means of exact diagonalization and analytic approaches at finite temperatures. Transport properties are studied first for the integrable XXZ\mathrm{XXZ}XXZ model and second for various nonintegrable systems such as the dimerized chain, the frustrated chain, and the spin ladder. We compare our results obtained by exact diagonalization and mean-field theory with those of the Bethe ansatz, bosonization, and other numerical studies in the case of the anisotropic Heisenberg model both in the gapless and gapped regime. In particular, we find indications that the Drude weight for spin transport is finite in the thermodynamic limit for the isotropic chain. For the nonintegrable models, a finite-size analysis of the numerical data for the Drude weights is presented, covering the entire parameter space of the dimerized and frustrated chain. We also discuss which conclusions can be drawn from bosonization regarding the question of whether the Drude weights are finite or not. One of our main results is that the Drude weights vanish in the thermodynamic limit for nonintegrable models.
We use the concept of typicality to study the real-time dynamics of spin and energy currents in spin-1/21/21/2 models in one dimension and at nonzero temperatures. These chains are the integrable XXZ chain and a nonintegrable modification due to the presence of a staggered magnetic field oriented in the zzz direction. In the framework of linear response theory, we numerically calculate autocorrelation functions by propagating a single pure state, drawn at random as a typical representative of the full statistical ensemble. By comparing to small-system data from exact diagonalization (ED) and existing short-time data from time-dependent density matrix renormalization group, we show that typicality is satisfied in finite systems over a wide range of temperature and is fulfilled in both integrable and nonintegrable systems. For the integrable case, we calculate the long-time dynamics of the spin current and extract the spin Drude weight for large systems outside the range of ED. We particularly provide strong evidence that the high-temperature Drude weight vanishes at the isotropic point. For the nonintegrable case, we obtain the full relaxation curve of the energy current and determine the heat conductivity as a function of magnetic field, exchange anisotropy, and temperature.
We illustrate how finite-temperature charge and thermal Drude weights of one-dimensional systems can be obtained from the relaxation of initial states featuring global (left–right) gradients in the chemical potential or temperature. The approach is tested for spinless interacting fermions as well as for the Fermi-Hubbard model, and the behavior in the vicinity of special points (such as half filling or isotropic chains) is discussed. We present technical details on how to implement the calculation in practice using the density matrix renormalization group and show that the non-equilibrium dynamics is often less demanding to simulate numerically and features simpler finite-time transients than the corresponding linear response current correlators; thus, new parameter regimes can become accessible. As an application, we determine the thermal Drude weight of the Hubbard model for temperatures T which are an order of magnitude smaller than those reached in the equilibrium approach. This allows us to demonstrate that at low T and half filling, thermal transport is successively governed by spin excitations and described quantitatively by the Bethe ansatz Drude weight of the Heisenberg chain.
We study one-dimensional spinless fermions at zero and finite temperature TTT using the density-matrix renormalization group. We consider nearest- as well as next-nearest-neighbor interactions; the latter render the system inaccessible by a Bethe ansatz treatment. Using an infinite-system algorithm we demonstrate the emergence of Luttinger liquid physics at low energies for a variety of static correlation functions as well as for thermodynamic properties. The characteristic power-law suppression of the momentum distribution n(k)n(k)n(k) function at T=0T=0T=0 can be directly observed over several orders of magnitude. At finite temperature, we show that n(k)n(k)n(k) obeys a scaling relation. The Luttinger liquid parameter and the renormalized Fermi velocity can be extracted from the density response function, the specific heat, or the susceptibility without the need to carry out any finite-size analysis. We illustrate that the energy scale below which Luttinger liquid power laws manifest vanishes as the half-filled system is driven into a gapped phase by large interactions.