SciPost Submission Page

The classical two-dimensional Heisenberg model revisited: An $SU(2)$-symmetric tensor network study

by Philipp Schmoll, Augustine Kshetrimayum, Jens Eisert, Roman Orus, Matteo Rizzi

Submission summary

As Contributors: Philipp Schmoll
Arxiv Link: (pdf)
Date submitted: 2021-06-29 08:59
Submitted by: Schmoll, Philipp
Submitted to: SciPost Physics
Academic field: Physics
  • Condensed Matter Physics - Computational
  • Quantum Physics
  • Statistical and Soft Matter Physics
Approaches: Theoretical, Computational


The classical Heisenberg model in two spatial dimensions constitutes one of the most paradigmatic spin models, taking an important role in statistical and condensed matter physics to understand magnetism. Still, despite its paradigmatic character and the widely accepted ban of a (continuous) spontaneous symmetry breaking, controversies remain whether the model exhibits a phase transition at finite temperature. Importantly, the model can be interpreted as a lattice discretization of the $O(3)$ non-linear sigma model in $1+1$ dimensions, one of the simplest quantum field theories encompassing crucial features of celebrated higher-dimensional ones (like quantum chromodynamics in $3+1$ dimensions), namely the phenomenon of asymptotic freedom. This should also exclude finite-temperature transitions, but lattice effects might play a significant role in correcting the mainstream picture. In this work, we make use of state-of-the-art tensor network approaches, representing the classical partition function in the thermodynamic limit over a large range of temperatures, to comprehensively explore the correlation structure for Gibbs states. By implementing an $SU(2)$ symmetry in our two-dimensional tensor network contraction scheme, we are able to handle very large effective bond dimensions of the environment up to $\chi_E^\text{eff} \sim 1500$, a feature that is crucial in detecting phase transitions. With decreasing temperatures, we find a rapidly diverging correlation length, whose behaviour is apparently compatible with the two main contradictory hypotheses known in the literature, namely a finite-$T$ transition and asymptotic freedom, though with a slight preference for the second.

Current status:
Editor-in-charge assigned

Submission & Refereeing History

You are currently on this page

Submission 2106.06310v2 on 29 June 2021

Reports on this Submission

Anonymous Report 1 on 2021-8-12 (Invited Report)


In this manuscript, the authors reported a large-scale tensor network numerical study of the two-dimensional (2D) classical Heisenberg model with O(3) symmetry.

The classical Heisenberg model is a paradigmatic model in statistical physics and condensed matter. While extensive analytical and numerical efforts were devoted to this model in the past, it remains a controversial issue whether the 2D classical O(3) Heisenberg model on the square lattice exhibits a Berezinskii-Kosterlitz-Thouless (BKT) phase transition at finite temperature.

Taking advantage of their tensor network library incorporating SU(2) symmetry, the authors address this important open problem. After casting the partition function into a tensor network, a symmetry-preserving corner transfer matrix renormalization group (CTMRG) approach is employed to (approximately) contract the tensor network in the thermodynamic limit. Within this framework, various physical quantities, including spin-spin correlation functions, thermodynamic entropy, and entanglement measures, have been computed. In particular, a rapidly diverging correlation length is observed at low temperature, which is consistent with the two contradictory hypotheses in the literature: finite-temperature BKT transition and asymptotic freedom (i.e., absence of finite-temperature transition).

With decreasing temperatures, the truncation of angular momentum basis in the partition function causes increasingly large errors. However, this is unavoidable for carrying out CTMRG contractions. It is then extremely difficult to conclusively answer whether the transition temperature is exactly zero or very low. Thus, this leaves room for future investigations.

In my opinion, the results reported in this manuscript surely represent state-of-the-art tensor network calculations and certainly shed new light on an important open problem. The manuscript is very well written. It would also motivate further development of tensor network methods and applying them to other paradigmatic models. Therefore, I would warmly recommend its publication in SciPost.

Below I have two minor comments:

1) In page 8, the claim that “the model is Bethe integrable” is a bit confusing. As far as I know, the computation of the partition function for the one-dimensional classical Heisenberg model relies on (directly) diagonalizing the transfer matrix in the angular momentum basis, so the Bethe ansatz technique is not used here, nor in Refs. [71,72].

2) In page 9, the tensor network representation of the classical O(3) Heisenberg model reminds me of an earlier paper [Liu et al., Phys. Rev. D 88, 056005 (2013)], where a similar tensor network representation was derived in Sec. III B. If they coincide, it might make sense to add a reference.

  • validity: top
  • significance: high
  • originality: high
  • clarity: high
  • formatting: perfect
  • grammar: excellent

Login to report or comment