Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Scaling of diffusion constants in the spin-12 XX ladder

R. Steinigeweg, F. Heidrich-Meisner, J. Gemmer, K. Michielsen, and H. De Raedt
Phys. Rev. B 90, 094417 – Published 25 September 2014
PDFHTMLExport Citation
Abstract
Authors
Article Text
  • INTRODUCTION
  • MODEL AND DEFINITIONS
  • NUMERICAL METHOD: DYNAMICAL TYPICALITY
  • RESULTS
  • CONCLUSIONS
  • ACKNOWLEDGMENTS
    • References

    Abstract

    We study the dynamics of spin currents in the spin-12 XX ladder at finite temperature. Within linear response theory, we numerically calculate autocorrelation functions for quantum systems larger than what is accessible with exact diagonalization using the concept of dynamical quantum typicality. While the spin Drude weight vanishes exponentially quickly with increasing system size, we show that this model realizes standard diffusive dynamics. Moreover, we unveil the existence of three qualitatively different dependencies of the spin-diffusion coefficient on the rung-coupling strength, resulting from a crossover from exponential to Gaussian dissipation as the rung coupling increases, in agreement with analytical predictions. We further discuss the implications of our results for experiments with cold atomic gases.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Received 13 June 2014
    • Revised 11 September 2014

    DOI:https://doi.org/10.1103/PhysRevB.90.094417

    ©2014 American Physical Society

    Authors & Affiliations

    R. Steinigeweg1,*, F. Heidrich-Meisner2, J. Gemmer3, K. Michielsen4,5, and H. De Raedt6

    • 1Institute for Theoretical Physics, Technical University Braunschweig, D-38106 Braunschweig, Germany
    • 2Department of Physics and Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, D-80333 Munich, Germany
    • 3Department of Physics, University of Osnabrück, D-49069 Osnabrück, Germany
    • 4Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, D-52425 Jülich, Germany
    • 5RWTH Aachen University, D-52056 Aachen, Germany
    • 6Department of Applied Physics, Zernike Institute for Advanced Materials, University of Groningen, NL-9747AG Groningen, Netherlands

    • *r.steinigeweg@tu-bs.de

    Article Text (Subscription Required)

    Click to Expand

    References (Subscription Required)

    Click to Expand
    Issue

    Vol. 90, Iss. 9 — 1 September 2014

    Reuse & Permissions
    Access Options
    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1

      Spin-current autocorrelation function C(t) for r=J⊥/J∥=1 and β=0: (a) Relaxation curve for large systems N≤36 and times tJ∥≤50. (b) Saturation at a very small Drude weight is only visible in a semilog plot of (a) for very long times tJ∥≫50. (c) Since the Drude weight C¯≲O(1%), D(t)→const for large times.

      Reuse & Permissions
    • Figure 2
      Figure 2

      High-temperature spin Drude weight C¯(t1,t2), extracted at very long times [t1J∥,t2J∥]=[300,400], for different coupling ratios r=J⊥/J∥=0.25, 0.5, 0.75, and 1 (symbols). The finite-size scaling follows an exponential A(r)e−γN [solid line, γ≠γ(r) and A(r)=1/r], suggesting C¯→0 for N→∞.

      Reuse & Permissions
    • Figure 3
      Figure 3

      Spin-diffusion constant D vs r=J⊥/J∥ for β→0. There are apparently three scaling regimes: (i) r≪1: D∝1/r2, (ii) 1≲r≲2: D∝1/r, and (iii) r≫1: D=const. The 1/r curve results from Eq. (7); see text.

      Reuse & Permissions
    • Figure 4
      Figure 4

      Spin-current autocorrelation function C(t) for β→0: qualitatively different regimes, depending on r=J⊥/J∥. (a) Weak-r regime: C(t) decays exponentially, resulting in D∝1/r2, as expected from perturbation theory. (b) Strong-r regime: The decay curve agrees with the Gaussian prediction of Eq. (7), in line with the generic behavior D∝1/r suggested in Refs. [39, 40]. For large r, revivals occur, resulting in D=const.

      Reuse & Permissions
    • Figure 5
      Figure 5

      Spin-current autocorrelation function C(t) in the XX ladder for rung couplings (a) r=0.5, (b) r=1.5, and (c) r=4 and different temperatures βJ∥≤1. In all cases, N=28. Apparently, C(t) is qualitatively the same for temperatures down to βJ∥∼0.5.

      Reuse & Permissions
    • Figure 6
      Figure 6

      Spin-current autocorrelation function C(t) in the XX ladder for rung couplings (a) r=0.25, (b) r=1.5, and (c) r=4 in the high-temperature limit β→0 and 〈Sz〉=0, Sz=0, Sz=1. In all cases, N=34. Apparently, the cases of half filling Sz=0 and almost half filling Sz=1 are already practically identical for finite N. No significant difference from 〈Sz〉=0 is visible for r=0.25, while differences for larger r are finite-size effects, as illustrated in Fig. 7.

      Reuse & Permissions
    • Figure 7
      Figure 7

      Spin-current autocorrelation function C(t) in the XX ladder for (a) 〈Sz〉=0 and (b) Sz=0 for a strong rung coupling r=J⊥/J∥=1.5 in the high-temperature limit β→0. (c) The resulting time-dependent diffusion coefficient, as extracted from C(t) depicted in (a) and (b). Although finite-size results for 〈Sz〉=0 and Sz=0 differ from each other, they seem to converge to the same value in the thermodynamic limit. Apparently, the convergence of the 〈Sz〉=0 data is much faster in time, and there are no finite-size effects up to times tJ∥∼10 when comparing N=28 and N=34. Hence, our choice of t2J∥=8.3 for extracting D from D(t) is reasonable.

      Reuse & Permissions
    ×

    Sign up to receive regular email alerts from Physical Review B

    Log In

    Cancel
    ×

    Search

    Article Lookup

    Paste a citation or DOI
    Enter a citation
    ×