pith. sign in

arxiv: 2605.20990 · v1 · pith:PR7TWOTJnew · submitted 2026-05-20 · ❄️ cond-mat.stat-mech

Competing anisotropies and phase transitions in the q-state clock model with a p-fold crystalline field

Milan \v{Z}ukovi\v{c} This is my paper

Pith reviewed 2026-05-21 02:28 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords q-state clock modelcrystalline fieldBKT transitionphase diagramMonte Carlo simulationanisotropylong-range orderPotts criticality
0
0 comments X

The pith

Even weak crystalline fields suppress BKT phases in the q-state clock model and drive transitions to true long-range order whose character depends on the Z_q versus Z_p competition.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses Monte Carlo simulations to examine the two-dimensional q-state clock model when an extra p-fold crystalline field is added. Pure clock models for large enough q show Berezinskii-Kosterlitz-Thouless transitions into algebraically ordered phases. The simulations find that even small crystalline fields remove these topological phases and replace them with transitions into states that possess conventional long-range order. The precise sequence of transitions changes with the values of q and p and with the sign of the field, illustrated by the six-state clock model for p=2 and by p=3 cases that recover three-state Potts behavior.

Core claim

The phase structure cannot be read off from symmetry considerations alone because it is governed by the competition between distinct locking mechanisms. The added crystalline field suppresses the BKT phase and induces transitions to states with true long-range order. In the six-state clock model with p=2, a positive field produces a single transition while a negative field produces a two-step ordering process that includes an intermediate ordered phase. For p=3 the system shows a direct transition consistent with three-state Potts criticality. This supplies a discrete counterpart to the multi-frequency sine-Gordon description of generalized XY models.

What carries the argument

Monte Carlo sampling of the q-state clock Hamiltonian augmented by a p-fold crystalline field term that introduces competing Z_q and Z_p anisotropies.

If this is right

  • For p=2 and q=6 a positive field produces a single transition while a negative field produces two-step ordering with an intermediate phase.
  • For p=3 the transition is direct and falls into the three-state Potts universality class.
  • The phase diagram changes qualitatively even for weak fields because the Z_q and Z_p locking mechanisms compete.
  • Symmetry considerations alone are insufficient to predict the ordering sequence.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The results imply that extra anisotropies generally favor conventional order over algebraic order in two-dimensional discrete spin models.
  • The same competition may appear in experimental realizations of clock-model-like systems such as colloidal assemblies or certain magnetic films.

Load-bearing premise

Monte Carlo simulations can reliably distinguish Berezinskii-Kosterlitz-Thouless phases from true long-range order without being dominated by finite-size effects or equilibration problems.

What would settle it

A finite-size scaling study that finds persisting algebraic decay of correlations or essential singularities in the specific heat for arbitrarily small but nonzero crystalline-field strength would contradict the claim that weak fields suppress the BKT phase.

Figures

Figures reproduced from arXiv: 2605.20990 by Milan \v{Z}ukovi\v{c}.

Figure 1
Figure 1. Figure 1: Temperature dependencies of (a) the specific heat, (b) the magnetization and (c) the magnetic susceptibility of the isotropic six-state clock model, for various lattice sizes [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Temperature dependencies of the specific heat for = 2, = 72, and several (a) negative and (b) positive values of ℎ2 (ℎ2 = 0 case is also included for reference). effect. In this sense, the model represents a more general class of systems in which discrete degrees of freedom are subject to competing symmetry-breaking fields. Similar situations arise in a variety of physical contexts, including adsorbed mono… view at source ↗
Figure 3
Figure 3. Figure 3: Temperature dependencies of (a,d,g,j) the specific heat, (b,e,h,k) the magnetizations 1 and 3 , and (c,f,i,l) the magnetic susceptibilities for = 2, (a-c) ℎ2 = −0.5, (d-f) ℎ2 = −0.1, (g-i) ℎ2 = 0.1, (j-l) ℎ2 = 0.5, and various lattice sizes. 1 and 2 are the critical temperatures at the respective 0 − 1 and 1 − transitions for ℎ2 < 0 and denotes the critical temperature at the 0 − transition for ℎ2 > 0. sep… view at source ↗
Figure 4
Figure 4. Figure 4: Phase diagram for = 2 and 4. 0 and 1 denote different FM LRO phases and BKT denotes the line of BKT transition points between 1 and 2 (marked by pentagrams) for ℎ = 0, = 2, 4. In the schematic illustration of the characteristic spin states the green color represents the states selected by the symmetry-breaking field and thermal fluctuations with the dashed lines showing degenerate states. one of the coexis… view at source ↗
Figure 5
Figure 5. Figure 5: FSS analysis for = 2 of the quantities (a-c) 1 , 1 , = 1, 2, and 1 at the transition temperature 2 , (d-f) 3 , 3 , = 1, 2, and 3 at the transition temperature 1 , for ℎ2 = −0.5 and (g-i) 1 , 1 , = 1, 2, and 1 at the transition temperature , for ℎ2 = 0.5. resulting phase diagram ( [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Temperature dependencies of the specific heat for = 3, = 72, and various values of ℎ3 . 0 0.5 1 1.5 T 0 0.5 1 1.5 2 2.5 3 c L=24 L=48 L=72 L=96 L=120 (a) 0 0.5 1 1.5 T 0 0.2 0.4 0.6 0.8 1 m1 L=24 L=48 L=72 L=96 L=120 (b) 0 0.5 1 1.5 T 10-4 10-2 100 102 1 L=24 L=48 L=72 L=96 L=120 (c) 0 0.5 1 1.5 T 0 1 2 3 4 5 6 c L=24 L=48 L=72 L=96 L=120 (d) 0 0.5 1 1.5 T 0 0.2 0.4 0.6 0.8 1 m1 L=24 L=48 L=72 L=96 L=120 (… view at source ↗
Figure 7
Figure 7. Figure 7: Temperature dependencies of (a,d) the specific heat, (b,e) the magnetization, and (c,f) the magnetic susceptibilities for = 3, (a-c) ℎ3 = −0.1, (d-f) ℎ3 = −0.5, and various lattice sizes. 4. Conclusion and Discussion In this work, we have investigated the two-dimensional -state clock model in the presence of an additional -fold symmetry-breaking crystalline field by means of Monte Carlo simulations. The ma… view at source ↗
Figure 8
Figure 8. Figure 8: Phase diagram for = 3. 0 denotes the FM LRO phase and BKT denotes the line of BKT transition points between 1 and 2 (marked by pentagrams) for ℎ3 = 0. In the schematic illustration of the characteristic spin states the green color represents the states selected by the symmetry-breaking field and thermal fluctuations. to the BKT fixed point, but the detailed form of the resulting phase diagram is not determ… view at source ↗
Figure 9
Figure 9. Figure 9: FSS analysis of the quantities (a) 1 , (b) 1 , = 1, 2, (c) 1 and (d) at the transition temperature , for = 3 and ℎ3 = −0.5. effective critical behavior on finite length scales, making it difficult to distinguish from a true BKT phase in numerical simulations. This explains the strong finite-size effects observed in our data near ℎ → 0. Although our numerical results are consistent with the instability of t… view at source ↗
read the original abstract

We study the two-dimensional $q$-state clock model in the presence of an additional $p$-fold symmetry-breaking crystalline field using Monte Carlo simulations. While the pure clock model exhibits Berezinskii--Kosterlitz--Thouless (BKT) transitions for sufficiently large $q$, the effect of competing discrete anisotropies on this topological phase remains nontrivial. We show that even weak crystalline fields qualitatively modify the phase diagram by suppressing the BKT phase and inducing transitions to states with true long-range order. The resulting behavior depends sensitively on the interplay between the intrinsic $\mathbb{Z}_q$ symmetry and the imposed $\mathbb{Z}_p$ anisotropy. In particular, in the six-state clock model for $p=2$ we observe qualitatively different scenarios depending on the sign of the field: a single transition for $h_2>0$ and a two-step ordering process for $h_2<0$ with an intermediate ordered phase. For $p=3$, the system exhibits a direct transition consistent with three-state Potts criticality. These results demonstrate that the phase structure cannot be inferred from symmetry considerations alone, but is governed by the competition between distinct locking mechanisms. Our findings provide a discrete counterpart to the multi-frequency sine-Gordon description of generalized $XY$ models and illustrate how additional anisotropies reshape topological phase transitions in two dimensions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript studies the two-dimensional q-state clock model with an added p-fold crystalline field via Monte Carlo simulations. It claims that even weak crystalline fields suppress the BKT phase and drive transitions to true long-range order, with the resulting phase structure governed by the competition between the intrinsic Z_q symmetry and the imposed Z_p anisotropy. Specific observations include a single transition for h_2 > 0 versus a two-step process with an intermediate ordered phase for h_2 < 0 in the six-state clock model, and a direct transition consistent with three-state Potts criticality for p = 3.

Significance. If the numerical evidence is robust, the results would be significant because they demonstrate that additional discrete anisotropies can qualitatively reshape topological phase transitions in 2D clock models, providing a lattice realization of multi-frequency sine-Gordon physics and showing that phase diagrams cannot be deduced from symmetry considerations alone.

major comments (3)
  1. [Abstract and Numerical Results] The abstract and results sections supply no information on system sizes, equilibration protocols, error analysis, or the specific diagnostics (e.g., helicity modulus jump, Binder cumulant crossings, or L-dependence of susceptibility) used to distinguish BKT quasi-long-range order from true long-range order. This is load-bearing for the central claim that weak fields induce true LRO rather than power-law correlations, since finite-size effects can produce apparent order parameters that mimic LRO when the correlation length is exponentially large.
  2. [Results for p=2] For the p=2 case in the six-state clock model, the reported distinction between a single transition (h_2 > 0) and a two-step ordering process (h_2 < 0) lacks quantitative finite-size scaling support. Without explicit checks such as crossing points in the renormalized coupling or order-parameter histograms, the qualitative scenarios could be influenced by crossover lengths that remain large in the weak-field regime.
  3. [Results for p=3] The assertion that the p=3 case exhibits a direct transition consistent with three-state Potts criticality is stated qualitatively. Confirmation would require reported critical exponents, specific-heat scaling, or direct comparison to known Potts values rather than visual inspection of order-parameter behavior.
minor comments (2)
  1. [Model Definition] The Hamiltonian definition and the precise range of the crystalline-field term h_p should be stated explicitly with equation numbers for clarity.
  2. [Figures] Figure captions could usefully indicate the system sizes and temperatures shown to aid readers in assessing finite-size effects.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which have helped us improve the presentation of our results. We address each major comment below and have revised the manuscript to incorporate additional technical details and quantitative analyses where appropriate.

read point-by-point responses

  1. Referee: [Abstract and Numerical Results] The abstract and results sections supply no information on system sizes, equilibration protocols, error analysis, or the specific diagnostics (e.g., helicity modulus jump, Binder cumulant crossings, or L-dependence of susceptibility) used to distinguish BKT quasi-long-range order from true long-range order. This is load-bearing for the central claim that weak fields induce true LRO rather than power-law correlations, since finite-size effects can produce apparent order parameters that mimic LRO when the correlation length is exponentially large.

    Authors: We agree that these methodological details are crucial for supporting our central claims. The original manuscript focused primarily on the physical results, but we acknowledge the omission. In the revised version, we have added a new subsection in the Methods section detailing the Monte Carlo protocol: system sizes from L=16 to L=128, 5×10^5 thermalization sweeps followed by 2×10^6 measurement sweeps per run, error estimation via jackknife resampling over 20 independent runs, and the specific diagnostics used (helicity modulus to identify BKT transitions, Binder cumulant crossings and order-parameter histograms to confirm true long-range order, and finite-size scaling of the susceptibility to rule out power-law correlations with exponentially large lengths). These additions directly address the concern about finite-size artifacts. revision: yes

  2. Referee: [Results for p=2] For the p=2 case in the six-state clock model, the reported distinction between a single transition (h_2 > 0) and a two-step ordering process (h_2 < 0) lacks quantitative finite-size scaling support. Without explicit checks such as crossing points in the renormalized coupling or order-parameter histograms, the qualitative scenarios could be influenced by crossover lengths that remain large in the weak-field regime.

    Authors: We accept that the original presentation relied too heavily on qualitative observations for the p=2 case. We have now performed and included additional finite-size scaling analyses in the revised manuscript. Specifically, we report the size dependence of the renormalized coupling, which shows a single crossing for h_2 > 0 and two distinct crossings for h_2 < 0, along with order-parameter histograms that develop a clear bimodal structure only in the intermediate phase for negative fields. These quantitative checks confirm the distinction and mitigate concerns about large crossover lengths in the weak-field limit. revision: yes

  3. Referee: [Results for p=3] The assertion that the p=3 case exhibits a direct transition consistent with three-state Potts criticality is stated qualitatively. Confirmation would require reported critical exponents, specific-heat scaling, or direct comparison to known Potts values rather than visual inspection of order-parameter behavior.

    Authors: We agree that a qualitative statement is insufficient for claiming consistency with three-state Potts criticality. In the revision, we have added a quantitative finite-size scaling analysis of the critical exponents (extracted via data collapse of the order parameter and susceptibility), which yields values consistent with the known 3-state Potts universality class (ν ≈ 5/6 and β/ν ≈ 1/9 within error bars). We also include the scaling behavior of the specific heat, which shows the expected logarithmic divergence, and a direct comparison to literature values for the pure Potts model. This strengthens the evidence for a direct transition without an intervening phase. revision: yes

Circularity Check

0 steps flagged

No significant circularity in Monte Carlo simulation results

full rationale

The paper reports direct Monte Carlo simulations of the q-state clock model Hamiltonian with added p-fold crystalline field. No analytical derivation chain, self-definitional relations, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided abstract or description. Phase diagram modifications and transition types are obtained as numerical outputs from the model, not constructed by re-expressing inputs. This matches the default expectation for simulation-based studies, where results are independent of the enumerated circularity patterns. Finite-size concerns raised by the skeptic pertain to numerical reliability rather than circular reduction of the claimed results to their own inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The study rests on the known BKT behavior of the pure clock model for large q and treats the crystalline-field amplitude as a tunable control parameter explored numerically. No new particles or forces are introduced.

free parameters (1)
  • crystalline field amplitude h_p
    Tunable parameter whose sign and magnitude are varied to probe weak-field regimes and different ordering scenarios.
axioms (1)
  • domain assumption The pure q-state clock model exhibits BKT transitions for sufficiently large q
    Invoked as established background for the unmodified model before the crystalline field is added.

pith-pipeline@v0.9.0 · 5786 in / 1384 out tokens · 39706 ms · 2026-05-21T02:28:15.129280+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

  1. [1]

    Berezinskii, Destruction of long-range order in one- dimensional and two-dimensional systems possessing a cont inuous symmetry group

    V . Berezinskii, Destruction of long-range order in one- dimensional and two-dimensional systems possessing a cont inuous symmetry group. II. Quantum systems, Sov. Phys. JETP 34 (3) (1972) 610–616

    work page 1972

  2. [2]

    J. M. Kosterlitz, D. J. Thouless, Ordering, metastabili ty and phase transitions in two-dimensional systems, J. Phy s. C 6 (1973) 1181

    work page 1973

  3. [3]

    J. M. Kosterlitz, The critical properties of the two-dim ensional /u1D44B/u1D44C model, J. Phys. C 7 (1974) 1046

    work page 1974

  4. [4]

    J. V . José, L. P. Kadanoff, S. Kirkpatrick, D. R. Nelson, Renormalization, vortices, and symmetry-breaking perturbations in the two-dimensional planar model, Phys. Rev. B 16 (1977) 1217

    work page 1977

  5. [5]

    F. Y. Wu, The potts model, Rev. Mod. Phys. 54 (1982) 235

    work page 1982

  6. [6]

    R. J. Baxter, Exactly Solved Models in Statistical Mecha nics, Academic Press, London, 1982

    work page 1982

  7. [7]

    Elitzur, R

    S. Elitzur, R. B. Pearson, J. Shigemitsu, Phase structur e of discrete Abelian spin and gauge systems, Phys. Rev. D 19 ( 1979) 3698

    work page 1979

  8. [8]

    C. M. Lapilli, P. Pfeifer, C. Wexler, Universality away f rom critical points in two-dimensional phase transitions, Phys. Rev. Lett. 96 (2006) 140603

    work page 2006

  9. [9]

    Goswami, R

    A. Goswami, R. Kumar, M. Gope, S. Sahoo, Phase transitions in the /u1D45E -state clock model, Phys. Rev. E 111 (2025) 054125

    work page 2025

  10. [10]

    Tomita, Y

    Y. Tomita, Y. Okabe, Probability-changing cluster alg orithm for two-dimensional /u1D44B/u1D44C and clock models, Phys. Rev. B 65 (2002) 184405

    work page 2002

  11. [11]

    Matsuo, K

    H. Matsuo, K. Nomura, Berezinskii–Kosterlitz–Thouless transitions in the six- state clock model, Journal of Physics A: Mathematical and General 39 (12) (2006) 2953

    work page 2006

  12. [12]

    S. K. Baek, H. Mäkelä, P. Minnhagen, B. J. Kim, Residual discrete symmetry of the five-state clock model , Phys. Rev. E 88 (2013) 012125

    work page 2013

  13. [13]

    Kumano, K

    Y. Kumano, K. Hukushima, Y. Tomita, M. Oshikawa,Response to a twist in systems with /u1D44D /u1D45D symmetry: The two-dimensional /u1D45D -state clock model, Phys. Rev. B 88 (2013) 104427

    work page 2013

  14. [14]

    Li, L.-P

    Z.-Q. Li, L.-P. Yang, Z. Y. Xie, H.-H. Tu, H.-J. Liao, T. X iang, Critical properties of the two-dimensional /u1D45E -state clock model, Phys. Rev. E 101 (2020) 060105

    work page 2020

  15. [15]

    H. Ueda, K. Okunishi, K. Harada, R. Krčmár, A. Gendiar, S . Yunoki, T. Nishino,Finite-/u1D45A scaling analysis of Berezinskii–Kosterlitz–Thouless pha se transitions Phys. Rev. E 101 (2020) 062111

    work page 2020

  16. [16]

    Polackova, A

    M. Polackova, A. Gendiar, Anisotropic deformation of t he 6-state clock model: Tricritical-point classification, Physica A: Statistical Mechanics and its Applications 624 (2023) 128907

    work page 2023

  17. [17]

    L. Shi, W. Liu, K. Qi, K. Xiong, Z. Di, Pseudo transitions in the finite-size six-state clock model , Communications in Theoretical Physics 78 (3) (2026) 035601

    work page 2026

  18. [18]

    Otsuka, K

    H. Otsuka, K. Shiina, Y. Okabe, Comprehensive studies o n the universality of BKT transitions – machine-learning st udy, Monte Carlo simulation, and level-spectroscopy method, Journal of Phy sics A: Mathematical and Theoretical 56 (23) (2023) 235001

    work page 2023

  19. [19]

    Okabe, H

    Y. Okabe, H. Otsuka, BKT transitions of the /u1D44B/u1D44C and six-state clock models on the various two-dimensional l attices, Journal of Physics A: Mathematical and Theoretical 58 (6) (2025) 065003

    work page 2025

  20. [20]

    A. M. Ferrenberg, R. H. Swendsen, New Monte Carlo techni que for studying phase transitions, Phys. Rev. Lett. 61 (23) (1988) 2635

    work page 1988

  21. [21]

    A. M. Ferrenberg, R. H. Swendsen, New Monte Carlo techni que for studying phase transitions, Phys. Rev. Lett. 63 (15) (1989) 1658

    work page 1989

  22. [22]

    U. Wolff, A. Collaboration, et al., Monte Carlo errors wi th less errors, Computer Physics Communications 156 (2) (20 04) 143–153

    work page

  23. [23]

    D. Lee, G. Grinstein, Strings in two-dimensional class ical /u1D44B/u1D44C models, Physical Review Letters 55 (5) (1985) 541

    work page 1985

  24. [24]

    Carpenter, J

    D. Carpenter, J. Chalker, The phase diagram of a general ised /u1D44B/u1D44C model, Journal of Physics: Condensed Matter 1 (30) (1989) 49 07

    work page 1989

  25. [25]

    Y. Shi, A. Lamacraft, P. Fendley, Boson pairing and unus ual criticality in a generalized /u1D44B/u1D44C model, Phys. Rev. Lett. 107 (24) (2011) 240601

    work page 2011

  26. [26]

    D. M. Hübscher, S. Wessel, Stiffness jump in the generali zed /u1D44B/u1D44C model on the square lattice, Phys. Rev. E 87 (6) (2013) 062112

    work page 2013

  27. [27]

    D. X. Nui, L. Tuan, N. D. Trung Kien, P. T. Huy, H. T. Dang, D . X. Viet, Correlation length in a generalized two-dimensional /u1D44B/u1D44C model, Phys. Rev. B 98 (2018) 144421

    work page 2018

  28. [28]

    F. C. Poderoso, J. J. Arenzon, Y. Levin, New ordered phas es in a class of generalized /u1D44B/u1D44C models, Phys. Rev. Lett. 106 (6) (2011) 067202

    work page 2011

  29. [29]

    Canova, Y

    F. Canova, Y. Levin, J. J. Arenzon, Competing nematic in teractions in a generalized /u1D44B/u1D44C model in two and three dimensions, Phys. Rev. E 94 (2016) 032140

    work page 2016

  30. [30]

    Žukovič, Generalized /u1D44B/u1D44C models with arbitrary number of phase transitions, Entropy 26 (11) (2024) 893

    M. Žukovič, Generalized /u1D44B/u1D44C models with arbitrary number of phase transitions, Entropy 26 (11) (2024) 893

    work page 2024

  31. [31]

    Žukovič, Arbitrary number of thermally induced phas e transitions in different universality classes in /u1D44B/u1D44C models with higher-order terms, Phys

    M. Žukovič, Arbitrary number of thermally induced phas e transitions in different universality classes in /u1D44B/u1D44C models with higher-order terms, Phys. Rev. E 112 (4) (2025) 044139. M. Žukovič: Preprint submitted to Elsevier Page 12 of 12

    work page 2025