© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1911.03753 · Physical review. E

Transport, correlations, and chaos in a classical disordered anharmonic chain

Manoj Kumar, Anupam Kundu, Manas Kulkarni, David A. Huse, Abhishek Dhar

Published Aug 20, 2020• 14 citations

Abstract

We explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size N, disorder strength Δ, and temperature T. The conductivity κ_{N}, obtained for finite length (N), saturates to a value κ_{∞}>0 in the large N limit, for all values of disorder strength Δ and temperature T>0. We show evidence that for any Δ>0 the conductivity goes to zero faster than any power of T in the (T/Δ)→0 limit, and find that the form κ_{∞}∼e^{-B|ln(CΔ/T)|^{3}} fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be κ_{N}<κ_{∞} due to boundary resistance when the bulk conductivity is high (the weak disorder case), or κ_{N}>κ_{∞} due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

DOI: 10.1103/physreve.102.022130

Related papers

Energy transport in Sachdev-Ye-Kitaev networks coupled to thermal baths

Cristian Zanoci, Brian G. Swingle

We develop a general framework for studying the equilibrium and non-equilibrium properties of arbitrary networks of Sachdev-Ye-Kitaev clusters coupled to thermal baths. We proceed to apply this technique to the problem of energy transport, which is known to be diffusive due to the strange metal behavior of these models. We use the external baths to impose a temperature gradient in the system and study the emerging non-equilibrium steady state using the Schwinger-Keldysh formalism. We consider two different configurations for the baths, implementing either a boundary or bulk driving, and show that the latter leads to a significantly faster convergence to the steady state. This setup allows us to compute both the temperature and frequency dependence of the diffusion constant. At low temperatures, our results agree perfectly with the previously known values for diffusivity in the conformal limit. We also establish a relationship between energy transport and quantum chaos by showing that the diffusion constant is upper bounded by the chaos propagation rate at all temperatures. Moreover, we find a simple analytical form for the non-equilibrium Green's functions in the linear response regime and use it to derive exact closed-form expressions for the diffusion constant in various limits. We mostly focus on uniform one-dimensional chains, but we also discuss higher-dimensional generalizations.

arXiv:2109.03268•Physical Review Research•Apr 1, 2022•14 citations

Dynamical conductivity of disordered quantum chains

Shintaro Takayoshi, Thierry Giamarchi

Abstract We study the transport properties of a one-dimensional quantum system with disorder. We numerically compute the frequency dependence of the conductivity of a fermionic chain with nearest-neighbor interaction and a random chemical potential by using the Chebyshev matrix product state (CheMPS) method. As a benchmark, we investigate the noninteracting case first. Comparison with exact diagonalization and analytical solutions demonstrates that the results of CheMPS are reliable over a wide range of frequencies. We then calculate the dynamical conductivity spectra of the interacting system for various values of the interaction and disorder strengths. In the high-frequency regime, the conductivity decays as a power law, with an interaction-dependent exponent. This behavior is qualitatively consistent with the bosonized field theory predictions, although the numerical evaluation of the exponent shows deviations from the analytically expected values. We also compute the characteristic pinning frequency at which a peak in the conductivity appears. We confirm that it is directly related to the inverse of the localization length, even in the interacting case. We demonstrate that the localization length follows a power law of the disorder strength with an exponent dependent on the interaction, and find good quantitative agreement with the field theory predictions. In the low-frequency regime, we find a behavior consistent with the one of the noninteracting system $$\omega ^{2}(\ln \omega )^{2}$$ ω 2 ( ln ω ) 2 independently of the interaction. We discuss the consequences of our finding for experiments in cold atomic gases. Graphical abstract

arXiv:2206.00023•The European Physical Journal D•Nov 1, 2022•3 citations

Transport of weakly disordered spin chains at infinite temperature

Ilia Khait, Snir Gazit, Norman Y. Yao +1 more

We study the disordered Heisenberg spin chain, which exhibits many body localization at strong disorder, in the weak to moderate disorder regime. A continued fraction calculation of dynamical correlations is devised, using a variational extrapolation of recurrents. Good convergence for the infinite chain limit is shown. We find that the local spin correlations decay at long times as $C \sim t^{-\beta}$, while the conductivity exhibits a low frequency power law $\sigma \sim \omega^{\alpha}$. The exponents depict sub-diffusive behavior $ \beta 0 $ at all finite disorders, and convergence to the scaling result, $\alpha+2\beta = 1$, at large disorders.

•Mar 21, 2016•0 citations

Spin transport of weakly disordered Heisenberg chain at infinite temperature

Ilia Khait, Snir Gazit, Norman Y. Yao +1 more

We study the disordered Heisenberg spin chain, which exhibits many-body localization at strong disorder, in the weak to moderate disorder regime. A continued fraction calculation of dynamical correlations is devised, using a variational extrapolation of recurrents. Good convergence for the infinite chain limit is shown. We find that the local spin correlations decay at long times as $C\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\beta}}$, whereas the conductivity exhibits a low-frequency power law $\ensuremath{\sigma}\ensuremath{\sim}{\ensuremath{\omega}}^{\ensuremath{\alpha}}$. The exponents depict subdiffusive behavior $\ensuremath{\beta}<1/2,\phantom{\rule{0.28em}{0ex}}\ensuremath{\alpha}>0$ at all finite disorders and convergence to the scaling result $\ensuremath{\alpha}+2\ensuremath{\beta}=1$ at large disorders.

arXiv:1603.06588•Physical review. B./Physical review. B•Jun 17, 2016•81 citations

Steady-state quantum transport through an anharmonic oscillator strongly coupled to two heat reservoirs

Tianqi Chen, Vinitha Balachandran, Chu Guo +1 more

We investigate the transport properties of an anharmonic oscillator, modeled by a single-site Bose-Hubbard model, coupled to two different thermal baths using the numerically exact thermofield based chain-mapping matrix product states (TCMPS) approach. We compare the effectiveness of TCMPS to probe the nonequilibrium dynamics of strongly interacting system irrespective of the system-bath coupling against the global master equation approach in Gorini-Kossakowski-Sudarshan-Lindblad form. We discuss the effect of on-site interactions, temperature bias as well as the system-bath couplings on the steady-state transport properties. Last, we also show evidence of non-Markovian dynamics by studying the nonmonotonicity of the time evolution of the trace distance between two different initial states.

arXiv:2004.05017•Physical review. E•Jul 31, 2020•24 citations

Heat conduction in one-dimensional aperiodic quantum Ising chains

Wenjuan Li, Peiqing Tong

The heat conductivity of nonperiodic quantum Ising chains whose ends are connected with heat baths at different temperatures are studied numerically by solving the Lindblad master equation. The chains are subjected to a uniform transverse field h, while the exchange coupling J{m} between the nearest-neighbor spins takes the two values J{A} and J{B} arranged in Fibonacci, generalized Fibonacci, Thue-Morse, and period-doubling sequences. We calculate the energy-density profile and energy current of the resulting nonequilibrium steady states to study the heat-conducting behavior of finite but large systems. Although these nonperiodic quantum Ising chains are integrable, it is clearly found that energy gradients exist in all chains and the energy currents appear to scale as the system size ~N{α}. By increasing the ratio of couplings, the exponent α can be modulated from α > -1 to α < -1 corresponding to the nontrivial transition from the abnormal heat transport to the heat insulator. The influences of the temperature gradient and the magnetic field to heat conduction have also been discussed.

Physical Review E•Mar 23, 2011•7 citations

We explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size N, disorder strength Δ, and temperature T. The conductivity κ_{N}, obtained for finite length (N), saturates to a value κ_{∞}>0 in the large N limit, for all values of disorder strength Δ and temperature T>0. We show evidence that for any Δ>0 the conductivity goes to zero faster than any power of T in the (T/Δ)→0 limit, and find that the form κ_{∞}∼e^{-B|ln(CΔ/T)|^{3}} fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be κ_{N}<κ_{∞} due to boundary resistance when the bulk conductivity is high (the weak disorder case), or κ_{N}>κ_{∞} due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

We develop a general framework for studying the equilibrium and non-equilibrium properties of arbitrary networks of Sachdev-Ye-Kitaev clusters coupled to thermal baths. We proceed to apply this technique to the problem of energy transport, which is known to be diffusive due to the strange metal behavior of these models. We use the external baths to impose a temperature gradient in the system and study the emerging non-equilibrium steady state using the Schwinger-Keldysh formalism. We consider two different configurations for the baths, implementing either a boundary or bulk driving, and show that the latter leads to a significantly faster convergence to the steady state. This setup allows us to compute both the temperature and frequency dependence of the diffusion constant. At low temperatures, our results agree perfectly with the previously known values for diffusivity in the conformal limit. We also establish a relationship between energy transport and quantum chaos by showing that the diffusion constant is upper bounded by the chaos propagation rate at all temperatures. Moreover, we find a simple analytical form for the non-equilibrium Green's functions in the linear response regime and use it to derive exact closed-form expressions for the diffusion constant in various limits. We mostly focus on uniform one-dimensional chains, but we also discuss higher-dimensional generalizations.

Abstract We study the transport properties of a one-dimensional quantum system with disorder. We numerically compute the frequency dependence of the conductivity of a fermionic chain with nearest-neighbor interaction and a random chemical potential by using the Chebyshev matrix product state (CheMPS) method. As a benchmark, we investigate the noninteracting case first. Comparison with exact diagonalization and analytical solutions demonstrates that the results of CheMPS are reliable over a wide range of frequencies. We then calculate the dynamical conductivity spectra of the interacting system for various values of the interaction and disorder strengths. In the high-frequency regime, the conductivity decays as a power law, with an interaction-dependent exponent. This behavior is qualitatively consistent with the bosonized field theory predictions, although the numerical evaluation of the exponent shows deviations from the analytically expected values. We also compute the characteristic pinning frequency at which a peak in the conductivity appears. We confirm that it is directly related to the inverse of the localization length, even in the interacting case. We demonstrate that the localization length follows a power law of the disorder strength with an exponent dependent on the interaction, and find good quantitative agreement with the field theory predictions. In the low-frequency regime, we find a behavior consistent with the one of the noninteracting system

\omega ^{2}(\ln \omega )^{2}

ω 2 ( ln ω ) 2 independently of the interaction. We discuss the consequences of our finding for experiments in cold atomic gases. Graphical abstract

We study the disordered Heisenberg spin chain, which exhibits many body localization at strong disorder, in the weak to moderate disorder regime. A continued fraction calculation of dynamical correlations is devised, using a variational extrapolation of recurrents. Good convergence for the infinite chain limit is shown. We find that the local spin correlations decay at long times as

C \sim t^{-\beta}

, while the conductivity exhibits a low frequency power law

\sigma \sim \omega^{\alpha}

. The exponents depict sub-diffusive behavior $ \beta 0 $ at all finite disorders, and convergence to the scaling result,

\alpha+2\beta = 1

, at large disorders.

We study the disordered Heisenberg spin chain, which exhibits many-body localization at strong disorder, in the weak to moderate disorder regime. A continued fraction calculation of dynamical correlations is devised, using a variational extrapolation of recurrents. Good convergence for the infinite chain limit is shown. We find that the local spin correlations decay at long times as $C\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\beta}}$, whereas the conductivity exhibits a low-frequency power law $\ensuremath{\sigma}\ensuremath{\sim}{\ensuremath{\omega}}^{\ensuremath{\alpha}}$. The exponents depict subdiffusive behavior $\ensuremath{\beta}<1/2,\phantom{\rule{0.28em}{0ex}}\ensuremath{\alpha}>0$ at all finite disorders and convergence to the scaling result $\ensuremath{\alpha}+2\ensuremath{\beta}=1$ at large disorders.

We investigate the transport properties of an anharmonic oscillator, modeled by a single-site Bose-Hubbard model, coupled to two different thermal baths using the numerically exact thermofield based chain-mapping matrix product states (TCMPS) approach. We compare the effectiveness of TCMPS to probe the nonequilibrium dynamics of strongly interacting system irrespective of the system-bath coupling against the global master equation approach in Gorini-Kossakowski-Sudarshan-Lindblad form. We discuss the effect of on-site interactions, temperature bias as well as the system-bath couplings on the steady-state transport properties. Last, we also show evidence of non-Markovian dynamics by studying the nonmonotonicity of the time evolution of the trace distance between two different initial states.

The heat conductivity of nonperiodic quantum Ising chains whose ends are connected with heat baths at different temperatures are studied numerically by solving the Lindblad master equation. The chains are subjected to a uniform transverse field h, while the exchange coupling J{m} between the nearest-neighbor spins takes the two values J{A} and J{B} arranged in Fibonacci, generalized Fibonacci, Thue-Morse, and period-doubling sequences. We calculate the energy-density profile and energy current of the resulting nonequilibrium steady states to study the heat-conducting behavior of finite but large systems. Although these nonperiodic quantum Ising chains are integrable, it is clearly found that energy gradients exist in all chains and the energy currents appear to scale as the system size ~N{α}. By increasing the ratio of couplings, the exponent α can be modulated from α > -1 to α < -1 corresponding to the nontrivial transition from the abnormal heat transport to the heat insulator. The influences of the temperature gradient and the magnetic field to heat conduction have also been discussed.