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:1612.04964 · Physical review. B./Physical review. B

Nonequilibrium thermal transport and vacuum expansion in the Hubbard model

Christoph Karrasch

Published Mar 27, 2017• 11 citations
arXivHTMLPDF

Abstract

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).

DOI
10.1103/physrevb.95.115148

Related papers

Nonequilibrium thermal transport and its relation to linear response

Christoph Karrasch, Roni Ilan, Joel E. Moore

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.

arXiv:1211.2236•Physical Review B•Nov 15, 2013•101 citations

Thermal Conductivity of the One-Dimensional Fermi-Hubbard Model

Christoph Karrasch, Dante Marvin Kennes, Fabian Heidrich‐Meisner

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.

arXiv:1506.05788•Physical Review Letters•Sep 9, 2016•24 citations

Diffusive high-temperature transport in the one-dimensional Hubbard model

Tomaž Prosen, Marko Žnidarič

We consider charge and spin transport in the one-dimensional Hubbard model at infinite temperature, half-filling, and zero magnetization. Implementing matrix-product-operator simulations of the nonequilibrium steady states of boundary-driven open Hubbard chains for up to 100 sites we find clear numerical evidence of diffusive transport for weak driving and any (nonzero and finite) value of the interaction $U$. In addition, nonequilibrium steady states are characterized by nonzero spin-up spin-down density-density correlations. For maximal driving, on the other hand, we find an insulating behavior and cosine-shaped density profiles.

arXiv:1203.1727•Physical Review B•Sep 14, 2012•71 citations

Time evolution of the Luttinger model with nonuniform temperature profile

Edwin Langmann A, Joel L. Lebowitz, Vieri Mastropietro +1 more

We study the time evolution of a one-dimensional interacting fermion system described by the Luttinger model starting from a nonequilibrium state defined by a smooth temperature profile $T(x)$. As a specific example we consider the case when $T(x)$ is equal to ${T}_{L}$ $({T}_{R})$ far to the left (right). Using a series expansion in $\ensuremath{\epsilon}=2({T}_{R}\ensuremath{-}{T}_{L})/({T}_{L}+{T}_{R})$, we compute the energy density, the heat current density, and the fermion two-point correlation function for all times $t\ensuremath{\ge}0$. For local (delta-function) interactions, the first two are computed to all orders, giving simple exact expressions involving the Schwarzian derivative of the integral of $T(x)$. For nonlocal interactions, breaking scale invariance, we compute the nonequilibrium steady state (NESS) to all orders and the evolution to first order in $\ensuremath{\epsilon}$. The heat current in the NESS is universal even when conformal invariance is broken by the interactions, and its dependence on ${T}_{L,R}$ agrees with numerical results for the $XXZ$ spin chain. Moreover, our analytical formulas predict peaks at short times in the transition region between different temperatures and show dispersion effects that, even if nonuniversal, are qualitatively similar to ones observed in numerical simulations for related models, such as spin chains and interacting lattice fermions.

arXiv:1701.06620•Physical review. B./Physical review. B•Jun 23, 2017•48 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

Quantum Nonequilibrium Dynamics: Transport, Entanglement, and Thermalization

HyungWon Kim

We study three aspects of quantum nonequilibrium dynamics; (1) transport of conserved quantities, (2) entanglement spreading, and (3) construction of local operators which slowly relax to thermal equilibrium. Motivated by recent progresses in ultracold atom experiments, we first analyze transport phenomena of a population imbalanced two-component fermi gas with arbitrary strength of inter-species interaction in three dimension. Using the Boltzmann kinetic equation at dilute regime, we obtain the transport coefficients of linear responses to gradients of temperature and chemical potential imbalance. We identify the magneto-caloric effects, and determine how these effects depend on interaction strength and population imbalance. Then, we propose an experimental protocol to observe these effects in an experiment with ultracold atoms. Next, we study entanglement spreading in a one-dimensional quantum Ising chain with longitudinal and transverse fields, which is diffusive and nonintegrable. Fully diagonalizing the Hamiltonian matrix, we explicitly show that the entanglement spreading is ballistic, thus faster than diffusive transport of conserved quantities. We provide a local spreading picture of entanglement entropy in terms of logarithmic negativity. Then, we discuss a role of energy conservation in entanglement spreading through analyzing a Floquet system. Lastly, we construct local operators that relax slowly to equilibrium with the same one-dimensional Hamiltonian. By showing that the Hamiltonian satisfies the eigenstate thermalization hypothesis, we first conclude that this model relaxes local operators to thermal equilibrium. Then, we systematically construct local operators that relax slower than conventional diffusive modes.

•Jan 1, 2014•1 citation
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 TL\ensuremath−TR{T}_{L}\ensuremath{-}{T}_{R}TL​\ensuremath−TR​. In one-dimensional systems which support ballistic thermal transport, the local energy currents \ensuremath⟨j(t)\ensuremath⟩\ensuremath{\langle}j(t)\ensuremath{\rangle}\ensuremath⟨j(t)\ensuremath⟩ acquire a nonzero value at long times, and it was recently investigated whether or not this steady state fulfills a simple additive relation \ensuremath⟨j(t\ensuremath→\ensuremath∞)\ensuremath⟩=f(TL)\ensuremath−f(TR)\ensuremath{\langle}j(t\ensuremath{\rightarrow}\ensuremath{\infty})\ensuremath{\rangle}=f({T}_{L})\ensuremath{-}f({T}_{R})\ensuremath⟨j(t\ensuremath→\ensuremath∞)\ensuremath⟩=f(TL​)\ensuremath−f(TR​) 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 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 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 ${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, 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.
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.
We consider charge and spin transport in the one-dimensional Hubbard model at infinite temperature, half-filling, and zero magnetization. Implementing matrix-product-operator simulations of the nonequilibrium steady states of boundary-driven open Hubbard chains for up to 100 sites we find clear numerical evidence of diffusive transport for weak driving and any (nonzero and finite) value of the interaction UUU. In addition, nonequilibrium steady states are characterized by nonzero spin-up spin-down density-density correlations. For maximal driving, on the other hand, we find an insulating behavior and cosine-shaped density profiles.
We study the time evolution of a one-dimensional interacting fermion system described by the Luttinger model starting from a nonequilibrium state defined by a smooth temperature profile T(x)T(x)T(x). As a specific example we consider the case when T(x)T(x)T(x) is equal to TL{T}_{L}TL​ (TR)({T}_{R})(TR​) far to the left (right). Using a series expansion in $\ensuremath{\epsilon}=2({T}_{R}\ensuremath{-}{T}_{L})/({T}_{L}+{T}_{R})$, we compute the energy density, the heat current density, and the fermion two-point correlation function for all times $t\ensuremath{\ge}0$. For local (delta-function) interactions, the first two are computed to all orders, giving simple exact expressions involving the Schwarzian derivative of the integral of T(x)T(x)T(x). For nonlocal interactions, breaking scale invariance, we compute the nonequilibrium steady state (NESS) to all orders and the evolution to first order in $\ensuremath{\epsilon}$. The heat current in the NESS is universal even when conformal invariance is broken by the interactions, and its dependence on TL,R{T}_{L,R}TL,R​ agrees with numerical results for the XXZXXZXXZ spin chain. Moreover, our analytical formulas predict peaks at short times in the transition region between different temperatures and show dispersion effects that, even if nonuniversal, are qualitatively similar to ones observed in numerical simulations for related models, such as spin chains and interacting lattice fermions.
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 three aspects of quantum nonequilibrium dynamics; (1) transport of conserved quantities, (2) entanglement spreading, and (3) construction of local operators which slowly relax to thermal equilibrium. Motivated by recent progresses in ultracold atom experiments, we first analyze transport phenomena of a population imbalanced two-component fermi gas with arbitrary strength of inter-species interaction in three dimension. Using the Boltzmann kinetic equation at dilute regime, we obtain the transport coefficients of linear responses to gradients of temperature and chemical potential imbalance. We identify the magneto-caloric effects, and determine how these effects depend on interaction strength and population imbalance. Then, we propose an experimental protocol to observe these effects in an experiment with ultracold atoms. Next, we study entanglement spreading in a one-dimensional quantum Ising chain with longitudinal and transverse fields, which is diffusive and nonintegrable. Fully diagonalizing the Hamiltonian matrix, we explicitly show that the entanglement spreading is ballistic, thus faster than diffusive transport of conserved quantities. We provide a local spreading picture of entanglement entropy in terms of logarithmic negativity. Then, we discuss a role of energy conservation in entanglement spreading through analyzing a Floquet system. Lastly, we construct local operators that relax slowly to equilibrium with the same one-dimensional Hamiltonian. By showing that the Hamiltonian satisfies the eigenstate thermalization hypothesis, we first conclude that this model relaxes local operators to thermal equilibrium. Then, we systematically construct local operators that relax slower than conventional diffusive modes.