Pith. Sign in

REVIEW

BKT transitions of the XY and six-state clock models on the various two-dimensional lattices

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Yutaka Okabe, Hiromi Otsuka

Statistical Mechanics · Submitted 13 Jan 2025 · 53 references

Pith's verdict: Useful new kagome/diced BKT temperatures, but the non-universality claim lacks the error bars needed to support it.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 10, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2501.07388 (version 1)
Pith IDT7J7AQDS, Pith's permanent record of this paper and its review
CitationHiromi Otsuka, Yutaka Okabe. "BKT transitions of the XY and six-state clock models on the various two-dimensional lattices." arXiv:2501.07388 (2025). https://pith.science/paper/T7J7AQDS
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

BKT transition temperatures of the XY and six-state clock models on honeycomb, kagome, and diced lattices track the Ising transition at a nearly constant, but not universal, ratio.

The reading

The paper aims to establish numerical BKT transition temperatures for two closely related two-dimensional spin models—the continuous XY model and its six-state discrete clock version—on honeycomb, kagome, and diced lattices, and to compare these with existing square- and triangular-lattice results. Using the ratio of correlation functions at two distances, the universal jump in the helicity modulus, and a machine-learning phase classifier, it reports $T_{\mathrm{BKT}}=0.573, 0.825, 0.965$ for the XY model on the three lattices, and $(T_1,T_2)=(0.455,0.579)$, $(0.645,0.835)$, and $(0.755,0.975)$ for the six-state clock model. Its central claim is that the ratios of these BKT temperatures to the exactly known Ising transition temperature on the same lattice are close across all five lattices—roughly 0.38 for the XY transition and 0.30–0.31 for the lower clock transition—but are not universal. A sympathetic reader would care because this gives a simple, lattice-insensitive rule of thumb for locating BKT transitions and sharpens the empirical relation between BKT physics and Ising criticality.

What carries the argument

The carrier of the numerical argument is the correlation ratio $R(T)=\langle g(L/2)\rangle/\langle g(L/4)\rangle$, where $g(r)$ is the two-spin correlation function; in the BKT phase this ratio obeys the single-variable finite-size scaling $R(T)=\tilde{f}(L/\xi)$ because the correlation function decays as a power law. At the transition the correlation length diverges as $\xi\propto \exp(c/\sqrt{|T-T_c|})$, so the paper searches for the temperature and constant $c$ that collapse data for system sizes $L=48,72,96,144,192$ onto one curve. Two independent cross-checks carry the same conclusion: the helicity modulus $\Upsilon(T)$ crosses the universal line $\Upsilon=2T/\pi$ at the upper transition, and a fully connected neural network trained on long-range correlation configurations classifies the ordered, BKT, and disordered phases with transition estimates consistent with the Monte Carlo results.

Extended reading notes

Core claim

On the paper's own terms, the discovery is a set of numerically estimated BKT temperatures together with a regularity: $T_{\mathrm{BKT}}/T_c^{\mathrm{Ising}}$ for the XY model lies in the narrow band 0.377–0.401 across square, triangular, honeycomb, kagome, and diced lattices, and $T_2/T_c^{\mathrm{Ising}}$ for the six-state clock model lies in 0.381–0.405, while $T_1/T_c^{\mathrm{Ising}}$ sits near 0.300–0.314. The values for the newly studied lattices are $T_{\mathrm{BKT}}=0.573$, $0.825$, and $0.965$ for the XY model on honeycomb, kagome, and diced lattices, and $(T_1,T_2)=(0.455,0.579)$, $(0.645,0.835)$, and $(0.755,0.975)$ for the six-state clock model on the same lattices. The paper explicitly concludes that the normalized values are not universal, even though they are nearly equal.

Load-bearing premise

The whole comparison rests on assuming that the correlation length diverges exactly with the BKT exponential form, with an exponent that goes like one over the square root of the distance from the transition, and that the curve-collapse procedure (carried out without quoted error bars) identifies the transition correctly; if either assumption is wrong, the reported temperatures and ratios could shift.

Editorial extensions

If this is right

  • For the honeycomb, kagome, and diced lattices, the paper supplies numerical BKT temperatures where previously only square- and triangular-lattice data were common: $T_{\mathrm{BKT}}=0.573$, $0.825$, and $0.965$ for the XY model.
  • The six-state clock model has two BKT transitions on each studied lattice, with $T_1/T_c^{\mathrm{Ising}}$ around 0.30–0.31 and $T_2/T_c^{\mathrm{Ising}}$ around 0.38–0.405; the upper transition consistently sits slightly above the XY $T_{\mathrm{BKT}}$ on the same lattice.
  • The normalized ratios are close but not exactly equal, so the paper rules out a simple universal constant while supporting the weaker statement that the ratio is lattice-insensitive to within a few percent.
  • The machine-learning classifier, fed with long-range correlation configurations rather than raw spin configurations, identifies the ordered, BKT, and disordered phases and gives transition estimates consistent with the Monte Carlo results.
  • The honeycomb XY estimate $0.573$ is consistent with two recent independent estimates, $0.576(1)$ and $0.560(9)$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Going beyond the paper: if the near-constant ratio holds on other lattices, it offers a cheap transfer rule—estimate a lattice's XY BKT temperature as roughly $0.39$ times its Ising $T_c$, and the lower six-state clock transition as roughly $0.30$ times—which is testable on lattices not studied here, such as the Lieb or square-octagon lattices.
  • Going beyond the paper: since the machine-learning classifier sees a slightly broader BKT phase than the correlation-ratio analysis, training it on larger configurations or on improved estimators could sharpen the ML boundary and resolve whether the broadening is purely a finite-size artifact.
  • Going beyond the paper: the near-universal ratio invites a theoretical explanation—for instance, via the spin-wave stiffness or a duality argument—that the paper does not attempt; a derivation that produces the ratio from lattice geometry would turn an empirical regularity into a prediction.
Original abstract

In a two-dimensional (2D) spin system, the XY model, characterized by planar rotational symmetry, exhibits a unique phenomenon known as the Berezinskii-Kosterlitz-Thouless (BKT) transition. In contrast, the clock model, which introduces discrete rotational symmetry, exhibits the BKT transition at two different temperatures due to this discreteness. In this study, we numerically investigate the BKT transition for XY and six-state clock models over various two-dimensional lattices. We employ two primary methods: the Monte Carlo method, which analyzes the size dependence of the ratio of the correlation functions for two different distances, and a machine-learning approach to classify the different phases -- namely, the low-temperature ordered phase, the intermediate BKT phase, and the high-temperature disordered phase. We identify the BKT transition temperatures for the XY and six-state models on honeycomb, kagome, and diced lattices. Combined with the previously calculated data for the triangular lattice, we then compare these values with the second-order phase transition temperatures of the 2D Ising model, for which exact solutions are known. Our results indicate that the ratio of the BKT transition temperatures for each lattice relative to the Ising model transition temperatures are close, although the values are not universal.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,734 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.73 · human participants 0.01 · controlled comparison 0.39 · randomized assignment 0.01 · states quantitative result 0.16 · compares with prior results 0.95 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

In a two-dimensional (2D) spin system, the XY model, characterized by planar rotational symmetry, exhibits a unique phenomenon known as the Berezinskii-Kosterlitz-Thouless (BKT) transition. In contrast, the clock model, which introduces discrete rotational symmetry, exhibits the BKT transition at two different temperatures due to this discreteness. In this study, we numerically investigate the BKT transition for XY and six-state clock models over various two-dimensional lattices. We employ two primary methods: the Monte Carlo method, which analyzes the size dependence of the ratio of the correlation functions for two different distances, and a machine-learning approach to classify the different phases -- namely, the low-temperature ordered phase, the intermediate BKT phase, and the high-temperature disordered phase. We identify the BKT transition temperatures for the XY and six-state models on honeycomb, kagome, and diced lattices. Combined with the previously calculated data for the triangular lattice, we then compare these values with the second-order phase transition temperatures of the 2D Ising model, for which exact solutions are known. Our results indicate that the ratio of the BKT transition temperatures for each lattice relative to the Ising model transition temperatures are close, although the values are not universal.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 In a two-dimensional (2D) spin system, the XY model, characterized by planar rotational symmetry, exhibits a unique phenomenon known as the Berezinskii-Kosterlitz-Thouless (BKT) transition. background 1.00 0.04 0.02 asserted 1.00
s1 In contrast, the clock model, which introduces discrete rotational symmetry, exhibits the BKT transition at two different temperatures due to this discreteness. background 1.00 0.05 0.04 asserted 1.00
s2 In this study, we numerically investigate the BKT transition for XY and six-state clock models over various two-dimensional lattices. contribution 1.00 0.94 0.03 asserted 0.90
s3 We employ two primary methods: the Monte Carlo method, which analyzes the size dependence of the ratio of the correlation functions for two different distances, and a machine-learning approach to classify the different phases -- namely, the low-temperature ordered phase, the intermediate BKT phase, and the high-temperature disordered phase. method 1.00 0.94 0.03 asserted 0.99
s4 We identify the BKT transition temperatures for the XY and six-state models on honeycomb, kagome, and diced lattices. result 0.72 0.96 0.08 asserted 1.00
s5 Combined with the previously calculated data for the triangular lattice, we then compare these values with the second-order phase transition temperatures of the 2D Ising model, for which exact solutions are known. contribution 0.58 0.70 0.08 asserted 0.78
s6 Our results indicate that the ratio of the BKT transition temperatures for each lattice relative to the Ising model transition temperatures are close, although the values are not universal. result 0.81 0.95 0.83 hedged possibility 0.92

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. The illustration of (a) triangular, (b) honeycomb, (c) kagome, (d) diced lattices. in
Figure 2. (left panel) The plot of the correlation ratio R(T ) of the XY model for the honeycomb, kagome, and diced lattice. The system sizes are 48, 72, 96, 144, and 192. In the inset, the FSS plots are given, where X(c, T ) = L/ exp(cBKT/( p |T − TBKT|)). (right panel) Helicity modulus for each lattice. We give the straight line (2/π)T . temperature for the honeycomb lattice is TBKT = 0.573. A rough estimates of the BKT temperature for the kagome and diced lattices are 0.825 and 0.965, respectively. These estimates of the BKT temperatures on the honeycomb, kagome and diced lattices are…
Figure 3. (left panel) The plot of the correlation ratio R(T ) of the six-state clock model for the honeycomb, kagome, and diced lattice. The system sizes are 48, 72, 96, 144, and 192. In the inset, the FSS plots are given, where X(c, T ) = L/ exp(c1,2/( p |T − T1,2|)). (right panel) Helicity modulus for each lattice. We give the straight line (2/π)T . a consequence of the finite size effect. The size used by the machine learning method is smaller than that of the Monte Carlo method. Additionally, the logarithmic correction might be less significant for the Monte Carlo method, which deals with
Figure 4. The machine-learning study of the six-state clock models on the honeycomb, kagome, and diced lattices. The output layer averaged over a test set as a function of T are plotted. The system sizes are L = 24, 48, and 72. For the honeycomb lattice, the samples of T in the ranges 0.32 ≤ T ≤ 0.41, 0.46 ≤ T ≤ 0.55 and 0.62 ≤ T ≤ 0.71 are used for training data. For the kagome lattice, the samples of T in the ranges 0.51 ≤ T ≤ 0.60, 0.69 ≤ T ≤ 0.78 and 0.88 ≤ T ≤ 0.97 are used for training data. For the diced lattice, the samples of T in the ranges 0.61 ≤ T ≤ 0.70, 0.79 ≤ T ≤ 0.88 and 0.98…

Figures reproduced from arXiv:2501.07388 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.