Finite-size anomalies of the Drude weight: Role of symmetries and ensembles

摘要

We revisit the numerical problem of computing the high temperature spin stiffness, or Drude weight, D of the spin-1/2 XXZ chain using exact diagonalization to systematically analyze its dependence on system symmetries and ensemble. Within the canonical ensemble and for states with zero total magnetization, we find D vanishes exacdy due to spin-inversion symmetry for all but the anisotropics ∆_(MN) = cos(πM/N) with N,M є Z~+ coprimes and N > M, provided system sizes L ≥ 2N, for which states with different spinTinversion signature become degenerate due to the underlying sl_2 loop algebra symmetry. All these loop-algebra degenerate states carry finite currents which we conjecture [based on data from the system sizes and anisotropics ≥_(MN) (with N < L/2) available to us] to dominate the grand-canonical ensemble evaluation of D in the thermodynamic limit. Including a magnetic flux not only breaks spin-inversion in the zero magnetization sector but also lifts the loop-algebra degeneracies in all symmetry sectors—this effect is more pertinent at smaller A due to the larger contributions to D coming from the low-magnetization sectors which are more sensitive to the system's symmetries. Thus we genetically find a finite D for fluxed ring's and arbitrary 0 < ≥ < 1 in both ensembles. In contrast, at the isotropic point and in the gapped phase (≥ ≥ 1) D is found to vanish in the thermodynamic limit, independent of symmetry or ensemble. Our analysis demonstrates how convergence to the thermodynamic limit within the gapless phase (∆ < 1) may be accelerated and the finite-size anomalies overcome: D extrapolates nicely in me thermodynamic limit to either the recently computed lower bound or the thermodynamic Bethe ansatz result provided both spin inversion is broken and the additional degeneracies at the ∆_(MN),v anisotropics are lifted.