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

Temperature-dependent energy diffusion in chaotic spin chains

Cristian Zanoci, Brian G. Swingle

Published Mar 25, 2021• 10 citations
arXivHTMLPDF

Abstract

We study the temperature dependence of energy diffusion in two chaotic gapped quantum spin chains, a tilted-field Ising model and an XZ model, using an open-system approach. We introduce an energy imbalance by coupling the chain to thermal baths at its boundary and study the nonequilibrium steady states of the resulting Lindblad dynamics using a matrix product operator ansatz for the density matrix. We define an effective local temperature profile by comparing local reduced density matrices in the steady state with those of a uniform thermal state. We then measure the energy current for a variety of driving temperatures and extract the temperature dependence of the energy diffusion constant. For the Ising model, we are able to study temperatures well below the energy gap and find a regime of dilute excitations where rare three-body collisions control energy diffusion. A kinetic model correctly predicts the observed exponential increase of the energy diffusion constant at low temperatures. For the XZ model, we are only able to access intermediate to high temperatures relative to the energy gap, and we show that the data are well described by an expansion around the infinite temperature limit. We also discuss the limitations of the particular driving scheme and suggest that lower temperatures can be accessed using larger baths.

DOI
10.1103/physrevb.103.115148

Related papers

Temperature dependence of energy transport in the Z 3 chiral clock model

Yongchan Yoo, Brian G. Swingle

We employ matrix product state simulations to study energy transport within the nonintegrable regime of the one-dimensional ${\mathbb{Z}}_{3}$ chiral clock model. To induce a nonequilibrium steady state throughout the system, we consider open system dynamics with boundary driving featuring jump operators with adjustable temperature and footprint in the system. Given a steady state, we diagnose the effective local temperature by minimizing the trace distance between the true local state and the local state of a uniform thermal ensemble. Via a scaling analysis, we extract the transport coefficients of the model at relatively high temperatures above both its gapless and gapped low-temperature phases. In the medium-to-high temperature regime we consider, diffusive transport is observed regardless of the low-temperature physics. We calculate the temperature dependence of the energy diffusion constant as a function of model parameters, including in the regime where the model is quantum critical at the low temperature. Notably, even within the gapless regime, an analysis based on power series expansion implies that intermediate-temperature transport can be accessed within a relatively confined setup. Although we are not yet able to reach temperatures where quantum critical scaling would be observed, our approach is able to access the transport properties of the model over a broad range of temperatures and parameters. We conclude by discussing the limitations of our method and potential extensions that could expand its scope, for example, to even lower temperatures.

arXiv:2311.00046•Physical review. B./Physical review. B•Jun 3, 2024•1 citation

Coexistence of energy diffusion and local thermalization in nonequilibrium XXZ spin chains with integrability breaking

Juan José Mendoza-Arenas, Stephen R. L. Clark, Dieter Jaksch

In this work we analyze the simultaneous emergence of diffusive energy transport and local thermalization in a nonequilibrium one-dimensional quantum system, as a result of integrability breaking. Specifically, we discuss the local properties of the steady state induced by thermal boundary driving in a XXZ spin chain with staggered magnetic field. By means of efficient large-scale matrix product simulations of the equation of motion of the system, we calculate its steady state in the long-time limit. We start by discussing the energy transport supported by the system, finding it to be ballistic in the integrable limit and diffusive when the staggered field is finite. Subsequently, we examine the reduced density operators of neighboring sites and find that for large systems they are well approximated by local thermal states of the underlying Hamiltonian in the nonintegrable regime, even for weak staggered fields. In the integrable limit, on the other hand, this behavior is lost, and the identification of local temperatures is no longer possible. Our results agree with the intuitive connection between energy diffusion and thermalization.

arXiv:1410.5838•Physical Review E•Apr 23, 2015•35 citations

Real-time and real-space spin and energy dynamics in one-dimensional spin- 1 2 systems induced by local quantum quenches at finite temperatures

Christoph Karrasch, Joel E. Moore, Fabian Heidrich‐Meisner

We study the spin and energy dynamics in one-dimensional spin-$\frac{1}{2}$ systems induced by local quantum quenches at finite temperatures using a time-dependent density matrix renormalization group method. System sizes are chosen large enough to ensure that the time-dependent data for the accessible time scales represent the behavior in the thermodynamic limit. As a main result, we observe a ballistic spreading of perturbations of the energy density in the integrable spin-$\frac{1}{2}$ XXZ chain for all temperatures and exchange anisotropies, related to the divergent thermal conductivity in this model and the exact conservation of the energy current. In contrast, the spin dynamics is ballistic in the massless phase, but shows a diffusive behavior at high temperatures in the easy-axis phase in the case of a vanishing background spin density. We extract a quantitative estimate for the spin-diffusion constant from the time dependence of the spatial variance of the spin density, which agrees well with values obtained from current-current correlation functions using an Einstein relation. Interestingly, the diffusion constant approaches a constant value deep in the easy-axis regime. As an example for nonintegrable models, we consider two-leg ladders, for which we observe indications of diffusive energy and spin dynamics. The relevance of our results for recent experiments with quantum magnets and bosons in optical lattices is discussed.

arXiv:1312.2938•Physical Review B•Feb 28, 2014•138 citations

Density dynamics in translationally invariant spin- 12 chains at high temperatures: A current-autocorrelation approach to finite time and length scales

Robin Steinigeweg, Jochen Gemmer B

We investigate transport in several translationally invariant spin-$\frac{1}{2}$ chains in the limit of high temperatures. We concretely consider spin transport in the anisotropic Heisenberg chain, the pure Heisenberg chain within an alternating field, and energy transport in an Ising chain which is exposed to a tilted field. Our approach is essentially based on a connection between the evolution of the variance of an inhomogeneous nonequilibrium density and the current-autocorrelation function at finite times. Although this relationship is not restricted to the case of diffusive transport, it allows to extract a quantitative value for the diffusion constant in that case. By means of numerically exact diagonalization we indeed observe diffusive behavior in the considered spin chains for a range of model parameters and confirm the diffusion coefficients which were obtained for these systems from nonequilibrium bath scenarios.

arXiv:0906.4457•Physical Review B•Nov 2, 2009•68 citations

Diffusive and Subdiffusive Spin Transport in the Ergodic Phase of a Many-Body Localizable System

Marko Žnidarič, Antonello Scardicchio, Vipin Kerala Varma

We study high temperature spin transport in a disordered Heisenberg chain in the ergodic regime. By employing a density matrix renormalization group technique for the study of the stationary states of the boundary-driven Lindblad equation we are able to study extremely large systems (400 spins). We find both a diffusive and a subdiffusive phase depending on the strength of the disorder and on the anisotropy parameter of the Heisenberg chain. Studying finite-size effects, we show numerically and theoretically that a very large crossover length exists that controls the passage of a clean-system dominated dynamics to one observed in the thermodynamic limit. Such a large length scale, being larger than the sizes studied before, explains previous conflicting results. We also predict spatial profiles of magnetization in steady states of generic nondiffusive systems.

arXiv:1604.08567•Physical Review Letters•Jul 22, 2016•287 citations

Dynamics of nonequilibrium thermal entanglement

Ilya Sinayskiy, Francesco Petruccione, Daniel Klaus Burgarth

The dynamics of a simple spin chain (two spins) coupled to bosonic baths at different temperatures is studied. The analytical solution for the reduced density matrix of the system is found. The dynamics and temperature dependence of spin-spin entanglement is analyzed. It is shown that the system converges to a steady state. If the energy levels of the two spins are different, the steady-state concurrence assumes its maximum at unequal bath temperatures. It is found that a difference in local energy levels can make the steady-state entanglement more stable against high temperatures.

arXiv:0807.0379•Physical Review A•Dec 2, 2008•78 citations
We study the temperature dependence of energy diffusion in two chaotic gapped quantum spin chains, a tilted-field Ising model and an XZ model, using an open-system approach. We introduce an energy imbalance by coupling the chain to thermal baths at its boundary and study the nonequilibrium steady states of the resulting Lindblad dynamics using a matrix product operator ansatz for the density matrix. We define an effective local temperature profile by comparing local reduced density matrices in the steady state with those of a uniform thermal state. We then measure the energy current for a variety of driving temperatures and extract the temperature dependence of the energy diffusion constant. For the Ising model, we are able to study temperatures well below the energy gap and find a regime of dilute excitations where rare three-body collisions control energy diffusion. A kinetic model correctly predicts the observed exponential increase of the energy diffusion constant at low temperatures. For the XZ model, we are only able to access intermediate to high temperatures relative to the energy gap, and we show that the data are well described by an expansion around the infinite temperature limit. We also discuss the limitations of the particular driving scheme and suggest that lower temperatures can be accessed using larger baths.
We employ matrix product state simulations to study energy transport within the nonintegrable regime of the one-dimensional Z3{\mathbb{Z}}_{3}Z3​ chiral clock model. To induce a nonequilibrium steady state throughout the system, we consider open system dynamics with boundary driving featuring jump operators with adjustable temperature and footprint in the system. Given a steady state, we diagnose the effective local temperature by minimizing the trace distance between the true local state and the local state of a uniform thermal ensemble. Via a scaling analysis, we extract the transport coefficients of the model at relatively high temperatures above both its gapless and gapped low-temperature phases. In the medium-to-high temperature regime we consider, diffusive transport is observed regardless of the low-temperature physics. We calculate the temperature dependence of the energy diffusion constant as a function of model parameters, including in the regime where the model is quantum critical at the low temperature. Notably, even within the gapless regime, an analysis based on power series expansion implies that intermediate-temperature transport can be accessed within a relatively confined setup. Although we are not yet able to reach temperatures where quantum critical scaling would be observed, our approach is able to access the transport properties of the model over a broad range of temperatures and parameters. We conclude by discussing the limitations of our method and potential extensions that could expand its scope, for example, to even lower temperatures.
In this work we analyze the simultaneous emergence of diffusive energy transport and local thermalization in a nonequilibrium one-dimensional quantum system, as a result of integrability breaking. Specifically, we discuss the local properties of the steady state induced by thermal boundary driving in a XXZ spin chain with staggered magnetic field. By means of efficient large-scale matrix product simulations of the equation of motion of the system, we calculate its steady state in the long-time limit. We start by discussing the energy transport supported by the system, finding it to be ballistic in the integrable limit and diffusive when the staggered field is finite. Subsequently, we examine the reduced density operators of neighboring sites and find that for large systems they are well approximated by local thermal states of the underlying Hamiltonian in the nonintegrable regime, even for weak staggered fields. In the integrable limit, on the other hand, this behavior is lost, and the identification of local temperatures is no longer possible. Our results agree with the intuitive connection between energy diffusion and thermalization.
We study the spin and energy dynamics in one-dimensional spin-12\frac{1}{2}21​ systems induced by local quantum quenches at finite temperatures using a time-dependent density matrix renormalization group method. System sizes are chosen large enough to ensure that the time-dependent data for the accessible time scales represent the behavior in the thermodynamic limit. As a main result, we observe a ballistic spreading of perturbations of the energy density in the integrable spin-12\frac{1}{2}21​ XXZ chain for all temperatures and exchange anisotropies, related to the divergent thermal conductivity in this model and the exact conservation of the energy current. In contrast, the spin dynamics is ballistic in the massless phase, but shows a diffusive behavior at high temperatures in the easy-axis phase in the case of a vanishing background spin density. We extract a quantitative estimate for the spin-diffusion constant from the time dependence of the spatial variance of the spin density, which agrees well with values obtained from current-current correlation functions using an Einstein relation. Interestingly, the diffusion constant approaches a constant value deep in the easy-axis regime. As an example for nonintegrable models, we consider two-leg ladders, for which we observe indications of diffusive energy and spin dynamics. The relevance of our results for recent experiments with quantum magnets and bosons in optical lattices is discussed.
We investigate transport in several translationally invariant spin-12\frac{1}{2}21​ chains in the limit of high temperatures. We concretely consider spin transport in the anisotropic Heisenberg chain, the pure Heisenberg chain within an alternating field, and energy transport in an Ising chain which is exposed to a tilted field. Our approach is essentially based on a connection between the evolution of the variance of an inhomogeneous nonequilibrium density and the current-autocorrelation function at finite times. Although this relationship is not restricted to the case of diffusive transport, it allows to extract a quantitative value for the diffusion constant in that case. By means of numerically exact diagonalization we indeed observe diffusive behavior in the considered spin chains for a range of model parameters and confirm the diffusion coefficients which were obtained for these systems from nonequilibrium bath scenarios.
We study high temperature spin transport in a disordered Heisenberg chain in the ergodic regime. By employing a density matrix renormalization group technique for the study of the stationary states of the boundary-driven Lindblad equation we are able to study extremely large systems (400 spins). We find both a diffusive and a subdiffusive phase depending on the strength of the disorder and on the anisotropy parameter of the Heisenberg chain. Studying finite-size effects, we show numerically and theoretically that a very large crossover length exists that controls the passage of a clean-system dominated dynamics to one observed in the thermodynamic limit. Such a large length scale, being larger than the sizes studied before, explains previous conflicting results. We also predict spatial profiles of magnetization in steady states of generic nondiffusive systems.
The dynamics of a simple spin chain (two spins) coupled to bosonic baths at different temperatures is studied. The analytical solution for the reduced density matrix of the system is found. The dynamics and temperature dependence of spin-spin entanglement is analyzed. It is shown that the system converges to a steady state. If the energy levels of the two spins are different, the steady-state concurrence assumes its maximum at unequal bath temperatures. It is found that a difference in local energy levels can make the steady-state entanglement more stable against high temperatures.