English

Comparing numerical methods for hydrodynamics in a one-dimensional lattice spin model

Strongly Correlated Electrons 2024-12-06 v2 Quantum Gases

Abstract

In ergodic quantum spin chains, locally conserved quantities such as energy or particle number generically evolve according to hydrodynamic equations as they relax to equilibrium. We investigate the complexity of simulating hydrodynamics at infinite temperature with multiple methods: time evolving block decimation (TEBD), TEBD with density matrix truncation (DMT), the recursion method with a universal operator growth hypothesis (R-UOG), and operator-size truncated (OST) dynamics. Density matrix truncation and the OST dynamics give consistent dynamical correlations to t=60/Jt=60/J and diffusion coefficients agreeing within 1%. TEBD only converges for t20t\lesssim 20, but still produces diffusion coefficients accurate within 1%. The universal operator growth hypothesis fails to converge and only matches other methods on short times. We see no evidence of long-time tails in either DMT or OST calculations of the current-current correlator, although we cannot rule out that they appear at longer times. We do observe power-law corrections to the energy density correlator. At finite wavelength, we observe a crossover from purely diffusive, overdamped decay of the energy density, to underdamped oscillatory behavior similar to that of cold atom experiments. We dub this behavior "hot band second sound" and offer a microscopically-motivated toy model.

Keywords

Cite

@article{arxiv.2310.06886,
  title  = {Comparing numerical methods for hydrodynamics in a one-dimensional lattice spin model},
  author = {Stuart Yi-Thomas and Brayden Ware and Jay D. Sau and Christopher David White},
  journal= {arXiv preprint arXiv:2310.06886},
  year   = {2024}
}

Comments

18 pages, 13 figures

R2 v1 2026-06-28T12:46:20.698Z