Order-by-disorder from Schwinger bosons in a frustrated honeycomb ferromagnet
Pith reviewed 2026-05-21 19:02 UTC · model grok-4.3
The pith
Generalized Schwinger-boson theory shows quantum fluctuations stabilize double-zigzag order via order-by-disorder in the J1-J3 honeycomb model.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The double-zigzag (dZZ) phase emerges in a narrow parameter range of the ferro-antiferromagnetic J1–J3 Heisenberg model on the honeycomb lattice. It is stabilized by quantum fluctuations through an order-by-disorder mechanism within the generalized Schwinger-boson mean-field theory. This finding agrees with recent DMRG calculations and characterizes the magnetic excitations relevant to inelastic neutron scattering on BaCo2(AsO4)2.
What carries the argument
Generalized Schwinger-boson mean-field theory (gSBMFT) that treats ferromagnetic and antiferromagnetic interactions equally to uncover the order-by-disorder selection of the dZZ phase.
Load-bearing premise
The generalized Schwinger-boson mean-field theory accurately captures the quantum fluctuations selecting the dZZ order without bias from decoupling scheme or bosonic representation.
What would settle it
A DMRG calculation on larger clusters that finds no stable dZZ phase in the predicted narrow parameter range would falsify the central claim.
read the original abstract
The cobalt-based honeycomb magnet BaCo$_2$(AsO$_4$)$_2$ (BCAO) has recently emerged as a promising platform for frustrated magnetism beyond conventional paradigms. Neutron-scattering experiments and first-principles calculations have revealed an unexpected double-zigzag (dZZ) magnetically ordered ground state, whose microscopic origin remains under active debate. Here, we revisit this problem within a ferro--antiferromagnetic $J_1$--$J_3$ Heisenberg model on the honeycomb lattice using a generalized Schwinger-boson mean-field theory (gSBMFT) that treats ferromagnetic and antiferromagnetic interactions on equal footing. This approach, combined with exact diagonalization (ED), allows us to demonstrate the emergence of the dZZ phase in a narrow parameter range, stabilized by quantum fluctuations through an order-by-disorder mechanism, in good agreement with recent density-matrix renormalization-group (DMRG) results. We further characterize the associated magnetic excitations and discuss their relevance to recent inelastic neutron-scattering (INS) measurements on BCAO.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the double-zigzag (dZZ) magnetic order observed in BaCo₂(AsO₄)₂ within a J₁–J₃ Heisenberg model on the honeycomb lattice, where J₁ is ferromagnetic and J₃ antiferromagnetic. Using a generalized Schwinger-boson mean-field theory (gSBMFT) that treats both interaction signs on equal footing, together with exact diagonalization, the authors demonstrate that quantum fluctuations select the dZZ phase via an order-by-disorder mechanism in a narrow parameter window, consistent with recent DMRG results. Magnetic excitations are characterized and compared to inelastic neutron scattering data on the material.
Significance. If the central result holds, the work supplies a microscopic mechanism for the unexpected dZZ ground state in BCAO and illustrates how order-by-disorder can resolve degeneracy in frustrated magnets with competing ferro- and antiferromagnetic couplings. The explicit agreement with DMRG benchmarks and the excitation spectra add value for interpreting experiments on honeycomb cobaltates.
major comments (2)
- [gSBMFT formulation and comparison to ED] gSBMFT formulation section: the claim that the saddle-point equations capture unbiased quantum fluctuations selecting dZZ rests on the specific decoupling of ferromagnetic J₁ versus antiferromagnetic J₃ bonds. Without explicit comparison to alternative decouplings or bosonic representations, it remains possible that the mean-field ansatz preferentially stabilizes dZZ, weakening the order-by-disorder interpretation.
- [Phase diagram and ED checks] Phase diagram and ED checks section: the dZZ phase occupies a narrow J₁/J₃ window. The absence of systematic finite-size scaling for the mean-field energies and lack of error bars on the reported energies make it difficult to confirm that dZZ is robustly selected over competing orders in this window.
minor comments (2)
- [gSBMFT formulation] The notation for the bosonic bond parameters and Lagrange multipliers could be unified across equations to improve readability.
- [Magnetic excitations] Figure captions for the excitation spectra should explicitly state the momentum path used in the calculations.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which help us clarify the presentation of our results. We address the two major comments point by point below.
read point-by-point responses
-
Referee: [gSBMFT formulation and comparison to ED] gSBMFT formulation section: the claim that the saddle-point equations capture unbiased quantum fluctuations selecting dZZ rests on the specific decoupling of ferromagnetic J₁ versus antiferromagnetic J₃ bonds. Without explicit comparison to alternative decouplings or bosonic representations, it remains possible that the mean-field ansatz preferentially stabilizes dZZ, weakening the order-by-disorder interpretation.
Authors: We appreciate this observation on the decoupling scheme. In our gSBMFT we adopt a bond-selective decoupling that assigns the ferromagnetic channel to J₁ bonds and the antiferromagnetic channel to J₃ bonds; this choice is dictated by the requirement that each mean-field parameter reflects the physical tendency of the corresponding interaction sign. We have added a short explanatory paragraph in the revised gSBMFT section justifying why a uniform decoupling would be unphysical for a mixed ferro-antiferromagnetic model and would raise the variational energy. The order-by-disorder interpretation is further corroborated by the quantitative agreement between the gSBMFT phase boundaries and both our ED spectra and the independent DMRG results cited in the manuscript. A exhaustive benchmark against every alternative bosonic representation lies outside the scope of the present work. revision: yes
-
Referee: [Phase diagram and ED checks] Phase diagram and ED checks section: the dZZ phase occupies a narrow J₁/J₃ window. The absence of systematic finite-size scaling for the mean-field energies and lack of error bars on the reported energies make it difficult to confirm that dZZ is robustly selected over competing orders in this window.
Authors: The narrow stability window of the dZZ phase is an intrinsic feature of the order-by-disorder mechanism in this model and is consistent with the DMRG phase diagram we reference. Because the gSBMFT saddle-point equations are solved in the thermodynamic limit, the mean-field energies themselves do not carry finite-size corrections. For the ED checks we have compared energies on clusters of 12, 18 and 24 sites; the relative ordering remains stable. We agree that the presentation can be improved by quantifying uncertainties. In the revision we add error bars derived from the variance across the different clusters and boundary conditions, together with a brief discussion of residual finite-size effects. revision: yes
Circularity Check
Standard gSBMFT with external DMRG/ED benchmarks shows no reduction to self-inputs
full rationale
The paper applies generalized Schwinger-boson mean-field theory to a J1-J3 honeycomb model and reports dZZ selection via order-by-disorder in a narrow window, with direct comparison to independent DMRG and ED results. No load-bearing step reduces the reported stability or mechanism to a parameter fitted inside the same equations or to a self-citation chain. The decoupling is presented as treating ferro- and antiferromagnetic bonds on equal footing without the result being forced by construction. This is the typical honest non-finding for a mean-field plus benchmark study.
Axiom & Free-Parameter Ledger
free parameters (1)
- J1/J3 ratio
axioms (2)
- domain assumption Spins can be faithfully represented by Schwinger bosons with the constraint enforced at mean-field level
- domain assumption The honeycomb lattice with only J1 and J3 interactions suffices to describe BCAO magnetism
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
generalized Schwinger-boson mean-field theory (g-SBMFT) that treats ferromagnetic and antiferromagnetic interactions on equal footing... triplet pairing amplitudes... order-by-disorder mechanism
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
-
[1]
−1 2. .2−1. . . .1 .(A1) If we now define the singlet operator ˆOA = 1√ 6 . . . . .1−1. .−1 1. . . . . ,(A2) it is straightforward to verify that ˆA† ˆA=− 1 4 ˆp† ij ˆO† A ˆTp ˆOAˆpij,(A3) with ˆA= 1 2(bi↑bj↓ −b i↓bj↑) the singlet pairing operator. In a similar way, the correspondence with the boson hopping ˆB= 1 2(b† i↑bj↑ +b † i↓bj↓) is g...
-
[2]
1 .(A5) This concludes the connection with the standard SBMFT since we have established the well known relationˆSiˆSj =: ˆB† ˆB:− ˆA† ˆAwhen the projections are used. It is thus clear that depending the definition of the constraint pairing and hopping operators ˆOp,h, any type of mean field can be constructed in a very flexible way. Annexe B: Diagona...
-
[3]
J. Chaloupka, G. Jackeli, and G. Khaliullin, Kitaev- Heisenberg Model on a Honeycomb Lattice : Possible Exotic Phases in Iridium Oxides A 2 IrO 3, Physical Re- view Letters105, 027204 (2010)
work page 2010
-
[4]
S. Sinn, C. H. Kim, B. H. Kim, K. D. Lee, C. J. Won, J. S. Oh, M. Han, Y. J. Chang, N. Hur, H. Sato, B.-G. Park, C. Kim, H.-D. Kim, and T. W. Noh, Electronic Structure of the Kitaev Materialα-RuCl 3 Probed by Photoemis- sion and Inverse Photoemission Spectroscopies, Scientific Reports6, 39544 (2016)
work page 2016
- [5]
-
[6]
Y. Kasahara, T. Ohnishi, Y. Mizukami, O. Tanaka, S. Ma, K. Sugii, N. Kurita, H. Tanaka, J. Nasu, Y. Mo- tome, T. Shibauchi, and Y. Matsuda, Majorana quanti- zation and half-integer thermal quantum Hall effect in a Kitaev spin liquid, Nature559, 227 (2018)
work page 2018
-
[7]
M. Songvilay, J. Robert, S. Petit, J. A. Rodriguez-Rivera, W. D. Ratcliff, F. Damay, V. Bal´ edent, M. Jim´ enez-Ruiz, P. Lejay, E. Pachoud, A. Hadj-Azzem, V. Simonet, and C. Stock, Kitaev interactions in the co honeycomb anti- ferromagnetsN a 3Co2SbO6 andN a 2Co2T eO6, Physical Review B102, 224429 (2020)
work page 2020
-
[8]
P. A. Maksimov, S. Jiang, L. P. Regnault, and A. L. Chernyshev, Strong kitaev interaction in baco 2(aso4)2, Phys. Rev. Lett.135, 066703 (2025)
work page 2025
-
[9]
A. Devillez, J. Robert, E. Lhotel, R. Ballou, C. Cave- nel, F. D. Romero, Q. Faure, H. Jacobsen, J. Lass, D. G. Mazzone, U. B. Hansen, M. Enderle, S. Raymond, S. D. Brion, V. Simonet, and M. Songvilay, Bond-dependent interactions and ill-ordered state in the honeycomb co- baltate baco2(aso4)2, arXiv:2503.22300 (2025)
-
[10]
L. P. Regnault and J. Rossat-Mignod, Phase transitions in quasi two-dimensional planar magnets (Springer Ne- therlands, 1990)
work page 1990
-
[11]
P. A. Maksimov, A. V. Ushakov, Z. V. Pchelkina, Y. Li, S. M. Winter, and S. V. Streltsov, Ab initio guided mini- mal model for the kitaev materialBaCo 2(AsO4)2 : Im- portance of direct hopping, third-neighbor exchange, and quantum fluctuations, Physical Review B106, 165131 (2022)
work page 2022
-
[12]
P. A. Maksimov, Proximity-induced sequence of field transitions in the kitaev candidateBaCo 2(AsO4)2, Phy- sical Review B108, L180405 (2023)
work page 2023
-
[13]
X. Liu and H.-Y. Kee, Non-kitaev versus kitaev honey- comb cobaltates, Physical Review B107, 054420 (2023)
work page 2023
-
[14]
S. Lee, S. Zhang, S. M. Thomas, L. Pressley, C. A. Bridges, E. S. Choi, V. S. Zapf, S. M. Winter, and M. Lee, Quantum order by disorder is a key to understanding the magnetic phases ofBaCo 2AsO4)2, npj Quantum Mate- rials10, 11 (2025)
work page 2025
-
[15]
S. Das, S. Voleti, T. Saha-Dasgupta, and A. Parame- kanti, Xy magnetism, kitaev exchange, and long-range frustration in theJ eff= 1 2 honeycomb cobaltates, Physi- cal Review B104, 134425 (2021)
work page 2021
- [16]
-
[17]
Auerbach, Interacting electrons and quantum magne- tism (Springer New York, NY, 2012)
A. Auerbach, Interacting electrons and quantum magne- tism (Springer New York, NY, 2012)
work page 2012
-
[18]
D. P. Arovas and A. Auerbach, Functional integral theo- ries of low-dimensional quantum heisenberg models, Phy- sical Review B38, 316 (1988)
work page 1988
-
[19]
S. Sachdev, Kagom´ e- and triangular-lattice heisenberg antiferromagnets : Ordering from quantum fluctuations and quantum-disordered ground states with unconfined 15 bosonic spinons, Physical Review B45, 12377 (1992)
work page 1992
-
[20]
Ferromagnetic frustrated spin systems on the square lattice: a Schwinger boson study
H. Feldner, D. C. Cabra, and G. L. Rossini, Ferromagne- tic frustrated spin systems on the square lattice : Schwin- ger boson study, Physical Review B84, 214406 (2011), arXiv : 1106.5018
work page internal anchor Pith review Pith/arXiv arXiv 2011
- [21]
-
[22]
J. C. Halimeh and M. Punk, Spin structure factors of chi- ral quantum spin liquids on the kagome lattice, Physical Review B94, 104413 (2016)
work page 2016
- [23]
-
[24]
B. Schneider, J. C. Halimeh, and M. Punk, Projective symmetry group classification of chiral𭟋 2 spin liquids on the pyrochlore lattice : Application to the spin- 1 2 xxz heisenberg model, Physical Review B105, 125122 (2022)
work page 2022
-
[25]
M. Kargarian, A. Langari, and G. A. Fiete, Unusual magnetic phases in the strong interaction limit of two- dimensional topological band insulators in transition me- tal oxides, Physical Review B86, 205124 (2012)
work page 2012
-
[26]
X. Yang and F. Wang, Schwinger boson spin-liquid states on square lattice, Phys. Rev. B94, 035160 (2016)
work page 2016
-
[27]
F. Wang and A. Vishwanath, Spin-liquid states on the triangular and kagom´ e lattices : A projective-symmetry- group analysis of schwinger boson states, Phys. Rev. B 74, 174423 (2006)
work page 2006
-
[28]
F. Wang, Schwinger boson mean field theories of spin li- quid states on a honeycomb lattice : Projective symmetry group analysis and critical field theory, Phys. Rev. B82, 024419 (2010)
work page 2010
-
[29]
J. Merino and A. Ralko, Role of quantum fluctuations on spin liquids and ordered phases in the heisenberg model on the honeycomb lattice, Physical Review B97, 205112 (2018)
work page 2018
- [30]
-
[31]
C. Wu, B. Chen, X. Dai, Y. Yu, and Z.-B. Su, Schwinger- boson mean-field theory of the heisenberg ferrimagnetic spin chain, Physical Review B60, 1057 (1999)
work page 1999
- [32]
- [33]
-
[34]
A. E. Trumper, L. O. Manuel, C. J. Gazza, and H. A. Ceccatto, Schwinger-Boson Approach to Quantum Spin Systems : Gaussian Fluctuations in the “Natural” Gauge, Physical Review Letters78, 2216 (1997)
work page 1997
- [35]
-
[36]
F. Wang and R. Tao, Schwinger-boson approach to the two-dimensional antiferromagnetic heisenberg model beyond the mean-field approximation, Physical Review B61, 3508 (2000)
work page 2000
-
[37]
L. O. Manuel, C. J. Gazza, A. E. Trumper, and H. A. Ceccatto, Heisenberg model with Dzyaloshinskii-Moriya interaction : A mean-field Schwinger-boson study, Physi- cal Review B54, 12946 (1996)
work page 1996
- [38]
- [39]
-
[40]
K. Mondal and C. Kadolkar, Schwinger boson mean- field theory of the kagome heisenberg antiferromagnet with dzyaloshinskii-moriya interactions, Physical Review B95, 134404 (2017)
work page 2017
-
[41]
A. Ralko and J. Merino, Chiral bosonic quantum spin liquid in the integer-spin heisenberg-kitaev model, Phy- sical Review B110, 134402 (2024)
work page 2024
-
[42]
D. Sasamoto and J. Nasu, Schwinger boson theory fors= 1 kitaev quantum spin liquids (2025), arXiv:2509.25761 [cond-mat.str-el]
-
[43]
L.-P. Regnault, C. Boullier, and J. Lorenzo, Polarized- neutron investigation of magnetic ordering and spin dy- namics in BaCo2(AsO4)2 frustrated honeycomb-lattice magnet, Heliyon4, 10.1016/j.heliyon.2018.e00507 (2018), publisher : Elsevier
-
[44]
J. C. Halimeh and M. Punk, Spin structure factors of chiral quantum spin liquids on the kagome lattice, Phys. Rev. B94, 104413 (2016)
work page 2016
-
[45]
H. Zhang and C. A. Lamas, Exotic disordered phases in the quantumJ 1-J2 model on the honeycomb lattice, Phys. Rev. B87, 024415 (2013)
work page 2013
-
[46]
D.-V. Bauer and J. O. Fjærestad, Schwinger-boson mean- field study of theJ 1−J2 heisenberg quantum antiferro- magnet on the triangular lattice, Phys. Rev. B96, 165141 (2017)
work page 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.