🔬 condensed matter

Fluctuations of topological charges in two-dimensional classical Heisenberg model through high-temperature and low-temperature expansions

This paper extends the analysis of topological charge fluctuations from the 2D XY model to the 2D classical Heisenberg model, demonstrating through high- and low-temperature expansions that these fluctuations scale with the area of a region at high temperatures but follow a perimeter law at low temperatures.

Original authors: Shan-Chang Tang

Published 2026-08-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Shan-Chang Tang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world made of tiny, invisible magnets, each one a spinning top that can point in any direction. In the realm of physics, scientists study how these tops behave when they are packed together in a flat, two-dimensional grid. Sometimes, these grids act like a calm, ordered army where everyone points the same way. Other times, they are a chaotic riot of spinning. The big question is: what happens when you heat them up or cool them down? Does the chaos just get worse, or does something magical happen where the order suddenly snaps into a new kind of structure? This is the story of phase transitions, a concept that explains why ice melts or water boils, but here, it's about the invisible dance of magnetic spins.

To understand the paper's journey, we need two main characters. First, there are "topological charges." Think of these as the number of times the spins twist around a central point, like a whirlpool in a pond. If you have a whirlpool spinning clockwise and another spinning counter-clockwise right next to it, they might cancel each other out. Second, there is the idea of "fluctuation." Imagine you are counting how many whirlpools are inside a specific square patch of the pond. If the water is calm, the number stays steady. If the water is turbulent, the number jumps around wildly. The paper asks: How wildly does this number jump as we change the temperature?

The author of this paper, Shan-Chang Tang, decided to investigate this question using a specific model called the "2D classical Heisenberg model." In this model, the spins are like 3D arrows that can point anywhere in space, not just flat on a piece of paper. While scientists already knew that a similar model (the XY model) has a famous transition driven by these whirlpools, the Heisenberg model is trickier. Some theories suggested that complex, stable "skyrmions" (like tiny, stable tornadoes of spins) might drive a similar transition here. However, previous computer simulations hinted that these tornadoes might not actually exist in this model, leaving a mystery: Is there a transition at all, and if so, what does it look like?

The paper tackles this mystery by using two different mathematical "flashlights" to look at the system from opposite ends of the temperature spectrum. They don't just run a computer simulation this time; they use pure math to calculate exactly how the "whirlpool count" fluctuates when the system is very hot and when it is very cold.

When they shine their light on the high-temperature side (where the spins are jiggling wildly), they found that the fluctuation of these topological charges grows in direct proportion to the area of the region they are looking at. Imagine a large, flat field. If you count the number of random, unconnected whirlpools in a square patch, the more ground you cover (the area), the more the count jumps around. It's like counting raindrops hitting a roof: the bigger the roof, the more the total number varies. This "area law" suggests that at high temperatures, the topological defects are free and unbound, wandering around independently.

Then, they switched their flashlight to the low-temperature side (where the spins are mostly calm and aligned). Here, the story changes completely. The math showed that the fluctuation is no longer proportional to the area, but to the perimeter (the edge length) of the region. Imagine a crowd of people holding hands in a circle. If you look at the edge of the circle, the wiggles and movements are most noticeable there. The "noise" of the topological charges is confined to the boundary. This "perimeter law" suggests that at low temperatures, the defects are tightly bound together, canceling each other out in the middle of the region and only showing their presence at the edges.

The paper concludes that these two distinct behaviors—area law at high heat and perimeter law at low cold—strongly suggest that a transition does indeed happen in the 2D Heisenberg model, much like the famous Kosterlitz-Thouless transition in the simpler XY model. However, the author is careful to note that while the math proves the existence of this different behavior on either side, they haven't calculated the exact temperature where the switch happens yet. They also explicitly rule out the idea that stable, complete skyrmions are the drivers here, noting that previous simulations found no evidence of them. Instead, the transition seems to be driven by the binding and unbinding of simpler, vortex-like defects.

In short, this paper uses advanced math to confirm that the 2D Heisenberg model has a hidden switch. At high temperatures, the topological chaos spreads everywhere (area), but at low temperatures, it gets locked down to the edges (perimeter). It's a bit like a party where, when the music is loud and fast, everyone is dancing wildly across the whole dance floor, but when the music slows down, everyone pairs up and stands still, only the people at the edge of the room are still swaying. The author has mapped out the dance moves for both extremes, proving that a transition exists, even if they haven't yet pinpointed the exact moment the music changes.

Technical Summary: Fluctuations of Topological Charges in the 2D Classical Heisenberg Model

Problem Statement
The paper addresses the nature of phase transitions in the two-dimensional (2D) classical Heisenberg model. While the Kosterlitz-Thouless (KT) transition in the 2D XY model is well-understood as being driven by the binding and unbinding of vortices, the situation in the Heisenberg model (where spins have three components) is more complex. Although skyrmions are often cited as potential topological defects driving a transition, they are not stable in the 2D Heisenberg model, and simulations have shown only vortex-like defects at low temperatures. Despite this, the author investigates whether the fluctuation of topological charges, defined as $\chi = \langle Q^2 \rangle - \langle Q \rangle^2$ , can serve as an indicator of a phase transition in the Heisenberg model, building upon previous Monte Carlo simulation results that suggested an area law at high temperatures and a perimeter law at low temperatures.

Methodology
The author provides a theoretical derivation of the topological charge fluctuations using two distinct expansion methods, avoiding the reliance on numerical simulations used in their prior work:

High-Temperature Expansion:
- The system is treated on a discrete lattice with the Hamiltonian $H = -J \sum_{\langle ij \rangle} S_i \cdot S_j$ .
- The average value of observables is expanded in powers of $\beta = 1/k_B T$ using a diagrammatic technique.
- The trace of spin functions is calculated hierarchically. The author defines $\mu_l = \text{Tr}(H^l)$ and $\nu_l = \text{Tr}(O[S]H^l)$ , where $O[S]$ represents the topological charge density squared.
- Specific diagrams are constructed to calculate the coefficients $\alpha_l$ up to the fourth order in $\beta$ . The calculation involves evaluating traces of products of spin components (e.g., $S_x^2 S_y^2$ ) using rotational invariance and specific integral identities (Appendix A).
- The fluctuation $\chi_L$ for an $L \times L$ region is expanded by summing correlations between plaquette charges $\rho_{ij}$ .
Low-Temperature Expansion:
- The system is treated in the continuous limit, assuming spins fluctuate around a specific direction (chosen as the $z$ -axis).
- The Hamiltonian is approximated as two decoupled massless free scalar fields ( $m_x, m_y$ ) representing small deviations from the ordered state.
- The topological charge $Q$ for a region $\Omega$ is expressed using Green's theorem as a line integral over the boundary $\partial \Omega$ .
- The fluctuation $\chi_L$ is calculated using the correlation functions of the scalar fields, $G_{xx}(|r-r'|)$ , which are derived via functional integration with infrared and ultraviolet cutoffs (Appendix B).

Key Results
The theoretical expansions yield the following scaling behaviors for the fluctuation of topological charges $\chi_L$ within a region of linear size $L$ :

High-Temperature Regime: The expansion up to the fourth order in $\beta$ reveals that the leading term of the fluctuation is proportional to $L^2$ . Since the area of the region scales as $L^2$ , the fluctuation obeys an area law. Specifically, the leading term is derived as $\chi_L \approx L^2 (4/9)$ .
Low-Temperature Regime: In the continuous limit, the calculation of the boundary integral yields a result where the fluctuation is proportional to $L$ . Since the perimeter of the region scales as $L$ , the fluctuation obeys a perimeter law. The derivation shows $\chi_L \propto L / (\beta^2 J^2)$ .

Significance and Claims
The paper claims to theoretically validate the results previously obtained via Monte Carlo simulations. By demonstrating that the fluctuation of topological charges transitions from an area law at high temperatures to a perimeter law at low temperatures, the author argues that this behavior strongly suggests the existence of a phase transition in the 2D classical Heisenberg model, analogous to the KT transition in the XY model.

The author explicitly states that while these results support the existence of a transition, the specific transition temperature has not been calculated in this work and remains a subject for further research. The paper does not propose new experimental applications or future implications beyond the theoretical confirmation of the scaling laws and the suggestion of a transition mechanism similar to the XY model.

More like this