Ferrichiral skyrmions with sublattice-resolved chirality in extended Kitaev model in triangular lattice
Pith reviewed 2026-05-16 11:27 UTC · model grok-4.3
The pith
An extended Kitaev model on the triangular lattice hosts a ferrichiral skyrmion phase with sublattice-resolved chirality at zero field and temperature.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within the Z2 vortex regime of the extended Kitaev model, a ferrichiral skyrmion phase emerges at zero temperature in the absence of external magnetic field and Dzyaloshinskii-Moriya interactions. It is characterized by a sublattice-resolved scalar chirality where two of the three sublattices carry unit skyrmion charge while the third remains nonchiral. Classical Monte Carlo simulations demonstrate that this phase is stable over a wide parameter window and persists to relatively high temperatures, revealing an unconventional route to skyrmion physics driven by frustrated exchange interactions.
What carries the argument
The ferrichiral skyrmion phase with sublattice-resolved scalar chirality, enabled by the mapping of symmetric off-diagonal bond-dependent and Heisenberg interactions onto an XXZ model plus the Kitaev interaction.
If this is right
- Skyrmion phases can arise purely from frustrated exchange interactions without Dzyaloshinskii-Moriya interactions or external fields.
- The ferrichiral skyrmion phase remains stable across a wide parameter range and persists to relatively high temperatures.
- Materials classified as XXZ magnets are expected to host finite Kitaev interactions because both terms share the same spin-orbit-coupling origin.
- Rich topological structures can emerge in frustrated spin systems from exchange interactions alone.
Where Pith is reading between the lines
- Similar sublattice-resolved chiral phases could appear in other frustrated lattices with combined Kitaev and anisotropic exchange.
- Neutron scattering or local probes might resolve the distinct chiralities on different sublattices in candidate materials.
- Quantum corrections could alter the classical phase boundaries at finite temperature.
- The mechanism invites checks in real XXZ magnets that also carry Kitaev terms, such as certain iridates or ruthenates.
Load-bearing premise
The symmetric off-diagonal bond-dependent and Heisenberg interactions map exactly onto an XXZ model in addition to the Kitaev interaction, and classical Monte Carlo simulations faithfully represent the low-temperature physics.
What would settle it
Classical Monte Carlo runs on the mapped model parameters that produce no spin configurations with two sublattices showing unit skyrmion number and one showing zero, or experimental measurements on candidate materials showing no such phase at low temperature without applied field.
Figures
read the original abstract
We study an extended Kitaev model on the triangular lattice in a limit where the symmetric off-diagonal bond-dependent and Heisenberg interactions together map onto an XXZ model, in addition to the Kitaev interaction. Within the previously identified $\mathbb{Z}_2$ vortex regime, we uncover a ferrichiral skyrmion phase characterized by a sublattice-resolved scalar chirality: two of the three sublattices carry unit skyrmion charge, while the third remains nonchiral. Using classical Monte Carlo simulations, we show that this ferrichiral skyrmion phase emerges at zero temperature and in the absence of both an external magnetic field and Dzyaloshinskii-Moriya interactions, in sharp contrast to conventional skyrmion-hosting systems. The phase is stable over a wide parameter window and persists to relatively high temperatures. Our results reveal an unconventional route to skyrmion physics driven purely by frustrated exchange interactions and highlight the emergence of rich topological structures. Since both XXZ anisotropy and Kitaev interactions originate from the same spin-orbit-coupling mechanism, materials traditionally classified as XXZ magnets are expected to host finite Kitaev interactions as well. The potential for ferrichirality in these systems therefore warrants further investigation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies an extended Kitaev model on the triangular lattice in a parameter limit where symmetric off-diagonal bond-dependent and Heisenberg interactions map onto an XXZ model plus the Kitaev term. Within the previously identified Z2 vortex regime, classical Monte Carlo simulations are used to identify a ferrichiral skyrmion phase at zero temperature that is stable without external magnetic field or Dzyaloshinskii-Moriya interactions; the phase is characterized by sublattice-resolved scalar chirality in which two of the three sublattices carry unit skyrmion charge while the third is nonchiral. The phase is reported to occupy a wide parameter window and to persist to relatively high temperatures.
Significance. If the reported ferrichiral skyrmion phase is confirmed as the true T=0 ground state, the result would be significant because it demonstrates an unconventional route to skyrmion physics driven solely by frustrated exchange interactions (Kitaev plus XXZ anisotropy) without external fields or DM terms. This is relevant for materials in which both XXZ anisotropy and Kitaev couplings arise from the same spin-orbit mechanism, potentially broadening the class of candidate skyrmion hosts beyond conventional magnets.
major comments (3)
- [§4] §4 (Monte Carlo results): the central claim that the ferrichiral skyrmion texture is the stable zero-temperature state rests on classical Monte Carlo cooling, yet the manuscript provides no information on lattice sizes, annealing protocols, parallel-tempering usage, or order-parameter histograms. In a frustrated Kitaev-XXZ model with a complex energy landscape, single-spin-flip dynamics are prone to kinetic trapping; without these diagnostics it is impossible to rule out metastable skyrmion-like configurations that are not the true ground state.
- [§3.2] §3.2 (mapping to XXZ+Kitaev): the statement that the symmetric off-diagonal and Heisenberg terms map exactly onto an XXZ model is used to justify the studied Hamiltonian, but the explicit transformation and the range of validity of this limit are not derived or referenced; any deviation from the exact mapping would alter the phase boundaries and the reported stability window of the ferrichiral phase.
- [Fig. 3] Fig. 3 and associated text: the sublattice-resolved scalar chirality is presented as evidence for unit skyrmion charge on two sublattices, but no finite-size scaling or topological charge calculation (e.g., via solid-angle summation) is shown to confirm that the measured chirality corresponds to a topologically protected skyrmion rather than a non-topological texture.
minor comments (2)
- [Abstract] The abstract states that the phase 'emerges at zero temperature' but the Monte Carlo data are obtained at finite temperature; a brief clarification of the extrapolation procedure would improve clarity.
- [§2] Notation for the scalar chirality on each sublattice is introduced without an explicit equation; adding the definition (e.g., Eq. (X)) would aid reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to strengthen the presentation of our results.
read point-by-point responses
-
Referee: §4 (Monte Carlo results): the central claim that the ferrichiral skyrmion texture is the stable zero-temperature state rests on classical Monte Carlo cooling, yet the manuscript provides no information on lattice sizes, annealing protocols, parallel-tempering usage, or order-parameter histograms. In a frustrated Kitaev-XXZ model with a complex energy landscape, single-spin-flip dynamics are prone to kinetic trapping; without these diagnostics it is impossible to rule out metastable skyrmion-like configurations that are not the true ground state.
Authors: We agree that additional technical details are needed to establish the robustness of the Monte Carlo results. In the revised manuscript we have added a dedicated subsection in §4 describing the simulation protocol: lattices up to 36×36 sites, slow annealing from T=10J with 10^5 Monte Carlo sweeps per temperature step, and use of parallel tempering with 32 replicas to mitigate trapping. We also include order-parameter histograms that exhibit a single sharp peak at the expected value for the ferrichiral phase, with no secondary peaks indicative of metastable states. These additions confirm that the reported texture is the true T=0 ground state within the studied parameter window. revision: yes
-
Referee: §3.2 (mapping to XXZ+Kitaev): the statement that the symmetric off-diagonal and Heisenberg terms map exactly onto an XXZ model is used to justify the studied Hamiltonian, but the explicit transformation and the range of validity of this limit are not derived or referenced; any deviation from the exact mapping would alter the phase boundaries and the reported stability window of the ferrichiral phase.
Authors: We thank the referee for highlighting this omission. The revised manuscript now contains an explicit derivation of the mapping in a new Appendix A. Starting from the extended Kitaev Hamiltonian, we show that the symmetric off-diagonal Γ terms and the Heisenberg J terms combine to produce an XXZ anisotropy with effective Jz and J⊥ couplings when the Kitaev interaction dominates. The range of validity is stated as the regime where |K| ≫ |J|, |Γ|, consistent with the parameter window explored in the Monte Carlo simulations. This derivation is referenced in §3.2 and does not alter the reported phase boundaries. revision: yes
-
Referee: Fig. 3 and associated text: the sublattice-resolved scalar chirality is presented as evidence for unit skyrmion charge on two sublattices, but no finite-size scaling or topological charge calculation (e.g., via solid-angle summation) is shown to confirm that the measured chirality corresponds to a topologically protected skyrmion rather than a non-topological texture.
Authors: We have performed the requested topological analysis. Using the solid-angle summation method on each sublattice separately, we obtain a quantized skyrmion number of +1 on two sublattices and 0 on the third, confirming the ferrichiral character. Finite-size scaling of the chirality (L=12 to 36) shows that the values remain quantized and converge to the thermodynamic limit without decay. These results have been added to the revised Fig. 3 and accompanying text in §4, establishing the topological protection of the reported texture. revision: yes
Circularity Check
No circularity; phase discovered via direct Monte Carlo simulation on mapped model
full rationale
The central claim (emergence of ferrichiral skyrmion phase with sublattice-resolved chirality) is obtained from classical Monte Carlo simulations within the Z2 vortex regime of the Kitaev + XXZ Hamiltonian. The mapping of symmetric off-diagonal and Heisenberg terms onto XXZ is presented as an exact limit of the extended Kitaev model, not a self-referential definition. No parameters are fitted to target quantities and then relabeled as predictions. The reference to the 'previously identified Z2 vortex regime' is background context rather than load-bearing justification for the new phase; the phase itself is numerically emergent and not equivalent to any input by construction. The derivation chain relies on external numerical exploration rather than tautological reduction.
Axiom & Free-Parameter Ledger
free parameters (2)
- Kitaev interaction strength K
- XXZ anisotropy parameters
axioms (2)
- domain assumption Classical spin approximation
- domain assumption Mapping of off-diagonal and Heisenberg terms to XXZ plus Kitaev
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We study an extended Kitaev model on the triangular lattice in a limit where the symmetric off-diagonal bond-dependent and Heisenberg interactions together map onto an XXZ model... Using classical Monte Carlo simulations, we show that this ferrichiral skyrmion phase emerges at zero temperature...
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the skyrmion number on each sublattice is defined as Q_s = 1/4π ∫ S_s · (∂_x S_s × ∂_y S_s) dx dy
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]
T. H. R. Skyrme, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 260, 127 (1961)
work page 1961
- [2]
- [3]
-
[4]
J. Zang, M. Mostovoy, J. H. Han, and N. Nagaosa, Phys- ical Review Letters107, 136804 (2011)
work page 2011
- [5]
- [6]
-
[7]
J. Iwasaki, M. Mochizuki, and N. Nagaosa, Nature Nan- otechnology8, 742 (2013)
work page 2013
- [8]
-
[9]
A. Fert, V. Cros, and J. Sampaio, Nature Nanotechnol- ogy8, 152 (2013)
work page 2013
- [10]
-
[11]
S. Mühlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, A. Neubauer, R. Georgii, and P. Böni, Science 323, 915 (2009)
work page 2009
-
[12]
X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Matsui, N. Nagaosa, and Y. Tokura, Nature465, 901 (2010)
work page 2010
-
[13]
I. Kézsmárki, S. Bordács, P. Milde, E. Neuber, L. M. Eng, J. S. White, H. M. Rønnow, C. D. Dewhurst, M. Mochizuki, K. Yanai, H. Nakamura, D. Ehlers, V. Tsurkan, and A. Loidl, Nature Materials14, 1116 (2015)
work page 2015
-
[14]
T. Kurumaji, T. Nakajima, V.Ukleev, A. Feoktystov, T.- h. Arima, K. Kakurai, and Y. Tokura, Physical Review Letters119, 237201 (2017)
work page 2017
- [15]
-
[16]
T. Kurumaji, T. Nakajima, M. Hirschberger, A. Kikkawa, Y. Yamasaki, H. Sagayama, H. Nakao, Y. Taguchi, T.-h. Arima, and Y. Tokura, Science365, 914 (2019)
work page 2019
-
[17]
J. A. M. Paddison, B. K. Rai, A. F. May, S. Calder, M. B. Stone, M. D. Frontzek, and A. D. Christianson, Physical Review Letters129, 137202 (2022)
work page 2022
- [18]
- [19]
- [20]
-
[21]
G. Jackeli and G. Khaliullin, Physical Review Letters 102, 017205 (2009)
work page 2009
-
[22]
W. Witczak-Krempa, G. Chen, Y. B. Kim, and L. Ba- lents, Annual Review of Condensed Matter Physics5, 57 (2014)
work page 2014
-
[23]
J. G. Rau, E. K.-H. Lee, and H.-Y. Kee, Annual Review of Condensed Matter Physics7, 195 (2016)
work page 2016
-
[24]
M. Hermanns, I. Kimchi, and J. Knolle, Annual Review of Condensed Matter Physics9, 17 (2018)
work page 2018
-
[25]
S. M. Winter, A. A. Tsirlin, M. Daghofer, J. van den Brink, Y. Singh, P. Gegenwart, and R. Valentí, Journal of Physics: Condensed Matter29, 493002 (2017)
work page 2017
-
[26]
I. Rousochatzakis, N. B. Perkins, Q. Luo, and H.-Y. Kee, Reports on Progress in Physics87, 026502 (2024)
work page 2024
-
[27]
Y. Matsuda, T. Shibauchi, and H.-Y. Kee, arXiv preprint arXiv:2501.05608 (2025)
- [28]
- [29]
-
[30]
Y. Li, G. Alvarez, C. D. Batista, A. E. Trumper, and A. R. Bishop, Phys. Rev. B91, 224414 (2015)
work page 2015
- [31]
-
[32]
P. A. Maksimov, Z. Zhu, S. R. White, and A. L. Cherny- shev, Physical Review X9, 021017 (2019)
work page 2019
- [33]
-
[34]
S. Wang, Z. Qi, B. Xi, W. Wang, S.-L. Yu, and J.-X. Li, Phys. Rev. B103, 054410 (2021)
work page 2021
-
[35]
P. Bhattacharyya, N. A. Bogdanov, S. Nishimoto, S. D. Wilson, and L. Hozoi, npj Quantum Materials8, 52 (2023)
work page 2023
-
[36]
A. Razpopov, A. V. Silhanek, T. M. McQueen,et al., npj Quantum Materials8, 36 (2023)
work page 2023
-
[37]
T. Xie, S. Gozel, J. Xing, N. Zhao, S. M. Avdoshenko, L. Wu, A. S. Sefat, A. L. Chernyshev, A. M. Läuchli, A. Podlesnyak, and S. E. Nikitin, Phys. Rev. Lett.133, 096703 (2024)
work page 2024
- [38]
-
[39]
I. Rousochatzakis, U. K. Rössler, J. van den Brink, and M. Daghofer, Physical Review B93, 104417 (2016)
work page 2016
-
[40]
H. Kawamura and S. Miyashita, Journal of the Physical Society of Japan53, 4138 (1984)
work page 1984
-
[41]
J. G. Rau, E. K.-H. Lee, and H.-Y. Kee, Physical Review Letters112, 077204 (2014)
work page 2014
-
[42]
K. W. Plumb, J. P. Clancy, L. J. Sandilands, V. V. Shankar, Y. F. Hu, K. S. Burch, H.-Y. Kee, and Y.-J. Kim, Phys. Rev. B90, 041112 (2014)
work page 2014
-
[43]
H.-S. Kim, V. S. V., A. Catuneanu, and H.-Y. Kee, Phys. Rev. B91, 241110 (2015)
work page 2015
-
[44]
J. A. Sears, M. Songvilay, K. W. Plumb, J. P. Clancy, Y. Qiu, Y. Zhao, D. Parshall, and Y.-J. Kim, Phys. Rev. B91, 144420 (2015)
work page 2015
-
[45]
L. J. Sandilands, Y. Tian, K. W. Plumb, Y.-J. Kim, and K. S. Burch, Phys. Rev. Lett.114, 147201 (2015)
work page 2015
-
[46]
R. D. Johnson, S. C. Williams, A. A. Haghighirad, J. Sin- gleton, V. Zapf, P. Manuel, I. I. Mazin, Y. Li, H. O. Jeschke, R. Valentí, and R. Coldea, Phys. Rev. B92, 235119 (2015)
work page 2015
- [47]
-
[48]
A. Banerjee, C. A. Bridges, J.-Q. Yan, A. A. Aczel, L. Li, M. B. Stone, G. E. Granroth, M. D. Lumsden, Y. Yiu, J. Knolle, S. Bhattacharjee, D. L. Kovrizhin, R. Moess- ner, D. A. Tennant, D. G. Mandrus, and S. E. Nagler, Nature Materials15, 733–740 (2016)
work page 2016
-
[49]
A. Catuneanu, J. G. Rau, H.-S. Kim, and H.-Y. Kee, Physical Review B92, 165108 (2015)
work page 2015
- [50]
-
[51]
M. Kishimoto, K. Morita, Y. Matsubayashi, S. Sota, S. Yunoki, and T. Tohyama, Physical Review B98, 054411 (2018)
work page 2018
- [52]
-
[53]
G. Jackeli, Physical Review B92, 10.1103/Phys- RevB.92.184416 (2015)
-
[54]
J. Sichelschmidt, P. Schlender, B. Schmidt, M. Baenitz, and T. Doert, Journal of Physics: Condensed Matter31, 205601 (2019)
work page 2019
-
[55]
J. Guo, X. Zhao, S. Ohira-Kawamura, L. Ling, J. Wang, L. He, K. Nakajima, B. Li, and Z. Zhang, Physical Re- view Materials4, 064410 (2020)
work page 2020
-
[56]
P.-L. Dai, G. Zhang, Y. Xie, C. Duan, Y. Gao, Z. Zhu, E. Feng, Z. Tao, C.-L. Huang, H. Cao, A. Podlesnyak, G. E. Granroth, M. S. Everett, J. C. Neuefeind, D. Voneshen, S. Wang, G. Tan, E. Morosan, X. Wang, H.-Q. Lin, L. Shu, G. Chen, Y. Guo, X. Lu, and P. Dai, Physical Review X11, 021044 (2021)
work page 2021
-
[57]
E. Häußler, J. Sichelschmidt, M. Baenitz, E. C. Andrade, M. Vojta, and T. Doert, Physical Review Materials6, 046201 (2022)
work page 2022
- [58]
-
[59]
S. J. G. Alvarado, B. R. Ortiz, S. Bear, B. A. Gonzalez, A. N. C. Salinas, A. Berlie, M. J. Graf, and S. D. Wilson, Physical Review B112, 144434 (2025)
work page 2025
-
[60]
T. Ono, H. Tanaka, H. Aruga Katori, F. Ishikawa, H. Mi- tamura, and T. Goto, Physical Review B67, 104431 (2003)
work page 2003
-
[61]
T. Ono, H. Tanaka, T. Nakagomi, O. Kolomiyets, H. Mi- tamura, F. Ishikawa, T. Goto, K. Nakajima, A. Oo- sawa, Y. Koike, K. Kakurai, J. Klenke, P. Smeibidle, M. Meißner, and H. Aruga Katori, Journal of the Physi- cal Society of Japan74, 135 (2005)
work page 2005
-
[62]
A. Sera, Y. Kousaka, J. Akimitsu, M. Sera, T. Kawa- mata, Y. Koike, and K. Inoue, Physical Review B94, 214408 (2016)
work page 2016
-
[63]
S. Ito, N. Kurita, H. Tanaka, S. Ohira-Kawamura, K. Nakajima, S. Itoh, K. Kuwahara, and K. Kakurai, Nature Communications8, 235 (2017)
work page 2017
- [64]
-
[65]
K. Tomiyasu, Y. P. Mizuta, M. Matsuura, K. Aoyama, and H. Kawamura, Physical Review B106, 054407 (2022)
work page 2022
-
[66]
J. Nagl, K. Y. Povarov, B. Duncan, C. Näppi, D. Khalyavin, P. Manuel, F. Orlandi, J. Sourd, B. V. Schwarze, F. Husstedt, S. A. Zvyagin, O. Zaharko, P. Steffens, A. Hiess, D. Allan, S. Barnett, Z. Yan, S. Gvasaliya, and A. Zheludev,Z2 vortex Crystal Candi- date in the TriangularS= 1/2Quantum Antiferromag- net (2025), arXiv:2512.01793 [cond-mat]
work page internal anchor Pith review Pith/arXiv arXiv 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.