Pith. sign in

REVIEW 3 major objections 4 minor 46 references

Anharmonic Collective Oscillations in Isotropic Spin Systems and their Spectroscopic Signatures

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Quartic spin oscillations arise generically in isotropic spiral magnets, and thermal fluctuations turn them into a spin-wave gap that scales as √T in the thermodynamic limit.

desk verdict A genuinely new twist on pseudo-Goldstone physics, but the thermodynamic-limit √T law rests on a fitted entropic coefficient the authors never derive. read the letter →

arxiv 2508.21211 v1 pith:HKGA3FBH submitted 2025-08-28 cond-mat.str-el

Anna Fancelli , Matías G. Gonzalez , Subhankar Khatua , Bella Lake , Michel J. P. Gingras , Jeffrey G. Rau , Johannes Reuther This is my paper · ORCID
classification cond-mat.str-el PACS 75.30.Ds75.10.Hk
keywords spinwavesquarticpotentialspiralsystemsfluctuation-inducedgapfinite-sizeeffectsorderbydisordermoleculardynamicsneutronscatteringsignatures
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Spin waves are usually pictured as small oscillations in a quadratic (harmonic) potential. This paper argues that isotropic spin systems with spiral ground states generically host collective oscillations of a different kind: amplitude fluctuations whose energy cost grows as the fourth power of the amplitude, with no fine-tuning of exchange couplings. Thermal fluctuations convert these quartic modes into a temperature-dependent spin-wave gap at the wave vector of the unselected spiral, the anti-Bragg point. In finite systems the gap grows as T^1/4 at low temperature and then crosses over to √T; because the quartic coefficient shrinks as 1/N, the crossover temperature vanishes in the thermodynamic limit, leaving √T as the observable prediction. If correct, this gives neutron scattering a sharp, parameter-free signature of anharmonicity in spiral magnets and extends order-by-disorder physics to systems without exact accidental zero modes.

What carries the argument

The central object is the quartic oscillation: a spin fluctuation perpendicular to the spiral plane at the anti-Bragg wave vector QaB of the unselected spiral, whose energy starts at quartic order because the δ² term vanishes by lattice symmetries. Around it sits the effective Hamiltonian H = p²/2m + λx⁴ + αT x², in which the αT x² term models entropy from coupling the O(1) quartic modes to the O(N) harmonic modes. Mean-field decoupling, x⁴ ≈ 6⟨x²⟩x², gives a self-consistent gap ∆ = √(αT/m)[1+√(1+2T*/T)]^{1/2} with crossover T* = 6λ/α². The load-bearing step is λ ~ 1/N, which follows from the 1/√N Fourier prefactor of the canonical spin representation and forces T* ~ 1/N, leaving √T scaling

What would settle it

Measure, by inelastic neutron scattering on a bulk spiral magnet such as ZnCr₂Se₄, the spin-wave gap at the anti-Bragg wave vector as a function of temperature below T_c: if the gap follows T^1/4 rather than √T in the bulk, the thermodynamic-limit claim fails. On the simulation side, extract λ(N) from the ground-state energy expansion and α(N) from fits for several system sizes: if α varies with N or T instead of being constant, the formula T* = 6λ/α² and the T* → 0 conclusion would need revision.

Watch

Extended reading notes

Core claim

Two arguments carry the paper. First, symmetry: in a Heisenberg model with a planar spiral ground state, tilting spins out of the spiral plane at the anti-Bragg wave vector QaB (a symmetry-related, unselected spiral) costs no quadratic energy — for the square-lattice model, e = eGS + 0.039δ⁴ + O(δ⁶), and the 3D cubic model behaves the same. Second, size scaling: thermal fluctuations turn this quartic mode into a gapped oscillator described by H = p²/2m + λx⁴ + αT x², whose mean-field solution gives ∆ ~ T^1/4 below T* = 6λ/α² and ∆ ~ √T above. Because λ ~ 1/N, T* ~ 1/N → 0, so in the thermodynamic limit only ∆ ~ √T survives.

Load-bearing premise

The load-bearing premise is that the entropic contribution is exactly a temperature-linear term αT x² with a coefficient α that does not depend on temperature or system size; the paper does not derive this term microscopically (α is fitted to simulation data), and the thermodynamic-limit prediction ∆ ~ √T follows directly from that linear-in-T form together with λ ~ 1/N.

Editorial extensions

If this is right

  • In a bulk spiral magnet the spin-wave gap at the anti-Bragg wave vector should grow as √T across the ordered phase, with the T^1/4 regime confined to small finite systems — a crossover accessible by varying system size.
  • The gap is largest just below the ordering transition, where the order parameter is already weak, so the spectroscopic signature should persist over a broad temperature window.
  • The √T gap does not require exact accidental ground-state zero modes; generic quartic potentials suffice, expanding the class of materials expected to show pseudo-Goldstone-type dynamics.
  • Inelastic neutron scattering on known spiral helimagnets such as ZnCr₂Se₄ and GdPtBi should reveal a temperature-dependent gap at the unselected spiral wave vector, providing direct evidence for anharmonic thermal spin dynamics.

Reading between the lines

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

  • The vanishing of the quadratic term is shown to be Hamiltonian-independent, so the same gap mechanism should appear in any isotropic system with symmetry-related degenerate spiral wave vectors, including incommensurate spirals beyond the commensurate five-site cases simulated here.
  • The thermodynamic-limit conclusion depends on having O(1) quartic modes coupled to O(N) harmonic modes; lattice geometries with a macroscopic number of quartic modes per unit cell could evade the 1/N suppression and retain a T^1/4 regime in the bulk — a testable extension.
  • The entropic coefficient α is fixed by fitting, not derived; a microscopic calculation that explicitly integrates out the harmonic modes would either confirm the linear-in-T form or reveal T- or N-dependence that would modify the predicted √T scaling.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies anharmonic, quartic spin oscillations in isotropic classical Heisenberg models with spiral ground states, in 2D (square lattice, J1-J2-J3) and 3D (cubic lattice with J1-J2-J4). It shows analytically that a perpendicular spin modulation at the anti-Bragg wave vector costs only quartic energy, and combines spin molecular dynamics with Monte Carlo simulations to extract a temperature-dependent gap at that wave vector. The numerical gap is fitted to a phenomenological effective oscillator Heff = p²/2m + λx⁴ + αT x², which yields Δ ∼ T^{1/4} at low T and Δ ∼ √T at higher T, with a crossover scale T* = 6λ/α². The paper argues that λ ∼ 1/N, so T* → 0 in the thermodynamic limit and only the √T gap survives; this is connected to pseudo-Goldstone physics and possible neutron-scattering observations.

Significance. If established, the result is significant: it provides a general mechanism—quartic spin oscillations coupled to a macroscopic number of harmonic modes—by which a spin-wave gap grows as √T in the thermodynamic limit without fine-tuned accidental zero modes. The paper contains several genuine strengths: a clean symmetry-based demonstration that the quadratic term vanishes for the anti-Bragg perturbation (SM S1), a separate numerical check that the gap scales linearly with the quartic perturbation strength (End Matter A, Fig. E1), and a clear finite-size scaling T* ∼ 1/N supported by the inset of Fig. 2(b). The use of a reduced-χ² criterion for the fits (SM S4) is also careful. However, the central thermodynamic-limit claim relies on a phenomenological entropic coefficient α whose microscopic origin and N-independence are not derived; this is the main load-bearing gap.

major comments (3)
  1. [Entropic effects, Eq. (5), and End Matter B] The thermodynamic-limit conclusion Δ ∼ √T depends on the entropic term αT x² having α independent of T and N, but α is never derived. End Matter B derives only the 1/N scaling of λ; SM S6 merely solves the assumed model self-consistently; α is extracted by fitting Eq. (6) to the simulation data. If α ∼ N^{-a}, then T* = 6λ/α² ∼ N^{2a-1}, which need not vanish for a > 1/2, changing the asymptotic law. The paper should either derive α microscopically (e.g., by explicitly integrating out the harmonic modes in the spin-wave expansion) or provide a direct numerical test of its N-independence and T-linearity, for example by measuring the effective x² coefficient for several N.
  2. [2D model and Mermin-Wagner, main text near Eq. (5)] In two dimensions, the ordered spiral state used to define the expansion behind Eq. (5) does not exist in the thermodynamic limit at any T > 0. The paper acknowledges this (footnote [20]) but still claims that only Δ ∼ √T survives in the thermodynamic limit for both 2D and 3D models. The finite-size simulations probe an effectively symmetry-broken regime, and the extrapolation N → ∞ for d = 2 requires justification. The conclusion as stated is rigorous only for the 3D model; the 2D case needs qualification or a separate argument.
  3. [Fig. 2 and SM S4: fitting of the scaling regimes] The apparent crossover from T^{1/4} to T^{1/2} is partly self-fulfilling: Eq. (6) is fitted to the data with two free parameters, and the same fitted function then defines T*. The low-T T^{1/4} behavior is independently supported by the raw log-log slopes and the δ-scaling argument, but the T^{1/2} regime is not directly observed; in Fig. 2(b) the red line is a plain √T function drawn by hand, and in the 3D data (L = 10, 15) the intermediate scaling is even less convincing. A direct analysis of local log-log slopes, or a collapse plot of Δ²/T versus T/T*, would make the claim of a surviving √T regime in the thermodynamic limit much stronger.
minor comments (4)
  1. [Fig. 2] Error bars are not visible in the log-log plots; given that the gap is obtained from Lorentzian fits, it would help to show representative error bars for at least one system size. The inset T* versus 1/N should state whether T* values come from the individual fits and include their uncertainties.
  2. [Reproducibility] No data availability statement or code is provided. Since the central conclusions rely on fits to simulation data, making the underlying gap values and fitting scripts available would substantially aid verification.
  3. [Notation, Eq. (1)] The notation ⟨i<j⟩n for nth-nearest-neighbor sums is unconventional; a word definition or a subscript like (n.n.) would improve readability. Also, the normalization factor N_i in Eq. (3) is introduced only in the text; consider defining it in the equation itself.
  4. [Reference [22]] The Supplemental Material URL is left as a placeholder; this should be resolved in the final version.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the quartic energy expansion and λ∼1/N are microscopically derived; the αT entropic term is a phenomenological input, not a circular reduction.

full rationale

The paper's derivation chain is largely self-contained. The quartic form of the anti-Bragg mode energy is established by an explicit microscopic expansion (SM S1: e(δ)=eGS+0.039δ^4+O(δ^6)), and the key size dependence λ∼1/N is derived from the canonical spin representation and the 1/√N Fourier normalization (End Matter B, Eq. (E9) and following). The T^{1/4} law follows from dimensional analysis of H4 and is independently supported by the sMD ∆∼δ result combined with equipartition (End Matter A). The √T law follows from the phenomenological Heff with an added αT x² term; the coefficient α is not derived from the microscopic Hamiltonian and is obtained by fitting Eq. (6) to the simulations. This means the thermodynamic-limit claim 'only √T survives' rests on an unverified assumption about the T- and N-independence of α—a correctness/extrapolation risk, not a circular reduction. The numerical simulations are independent data, the fitted functional form could in principle have been rejected, and the λ∼1/N result is derived rather than assumed. Self-citations [12,14,15] are contextual and non-load-bearing. Therefore no significant circularity is present.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central derivation introduces no new microscopic entities. The load-bearing extra content is the fitted entropic coefficient α and the mean-field decoupling within the phenomenological oscillator model.

free parameters (2)
  • α (entropic coupling) = Not reported; fitted to gap vs T data via Eq. (6)
    Coefficient of the thermal x² term in the effective Hamiltonian (5). Sets T* = 6λ/α² and controls the crossover between T^{1/4} and T^{1/2}. Determined by fitting Eq. (6) to the numerical gap data in Fig. 2.
  • m (effective mass) = Not reported; enters via prefactor sqrt(αT/m)
    Effective mass of the quartic oscillator in Eq. (5). Appears with α as the overall prefactor in the gap formula (6); fitted to the same data.
assumptions (4)
  • domain assumption Spins are classical unit vectors (S_i ∈ S²)
    The paper states 'the spins are treated classically as three-dimensional vectors with unitary norm'. All simulations and derivations use this classical description; quantum effects are deferred.
  • domain assumption Mean-field decoupling x⁴ ≈ 6⟨x²⟩x² in the effective quartic oscillator
    SM S6 uses this decoupling to close the self-consistency equation for ⟨x²⟩ and derive the gap formula (6). It is an uncontrolled approximation for a single-mode quartic oscillator.
  • ad hoc to paper The entropic contribution to the effective potential is αT x² with α independent of T and N
    Eq. (5) adds this term, described as arising from integrating out other spin-wave modes. No microscopic derivation of α is provided; α is fitted. The thermodynamic-limit ∆ ~ sqrt(T) conclusion follows directly from this linear-in-T form and the 1/N scaling of λ.
  • domain assumption The quartic mode is fully described by retaining only the Fourier components at QaB and -QaB
    End Matter Section B restricts the expansion to k = ±QaB to isolate the quartic mode. This assumes the other modes contribute only through the effective entropic term and do not hybridize with the quartic mode.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Anharmonic Collective Oscillations in Isotropic Spin Systems and their Spectroscopic Signatures." pith.science (2026). https://pith.science/paper/HKGA3FBH

@misc{pith2026250821211,
  author       = {Pith},
  title        = {Pith review of: Anharmonic Collective Oscillations in Isotropic Spin Systems and their Spectroscopic Signatures},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HKGA3FBH}},
  note         = {Machine review of arXiv:2508.21211}
}
read the original abstract

Spin waves are the fundamental excitations in magnetically ordered spin systems and are ubiquitously observed in magnetic materials. However, the standard understanding of spin waves as collective spin oscillations in an effective harmonic potential does not consider the possibility of soft modes, such as those due to an effective quartic potential. In this work, we show that such quartic potentials arise under very general conditions in a broad class of isotropic spin systems without a fine-tuning of the interaction parameters. Considering models with spin spiral ground states in two and three spatial dimensions, we numerically demonstrate that quartic amplitude spin oscillations produce a fluctuation-induced spin-wave gap which grows with temperature according to a characteristic power-law. In conjunction with a phenomenological theory, the present work provides a general theoretical framework for describing soft spin modes, extending the previously discussed spin dynamics in the presence of order-by-disorder, and highlighting the important role of finite-size effects. Our predictions of a temperature-dependent gap in spiral spin systems could be tested in inelastic neutron scattering experiments, providing direct spectroscopic evidence for thermal effects arising from soft spin modes in magnetic materials.

Figures

Figures reproduced from arXiv: 2508.21211 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Zero temperature phase diagram for the [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a) Normalized dynamical structure factor [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

46 extracted references · 41 canonical work pages

  1. [20]

    D. S. Inosov, Y . O. Onykiienko, Y . V . Tymoshenko, A. Akopyan, D. Shukla, N. Prasai, M. Doerr, D. Gorbunov, S. Zherlitsyn, D. J. V oneshen, M. Boehm, V . Tsurkan, V . Fe- lea, A. Loidl, and J. L. Cohn, Magnetic field dependence of low-energy magnons, anisotropic heat conduction, and sponta- neous relaxation of magnetic domains in the cubic helimagnet Zn...

    work page 2020

  2. [1]

    2(b) and 2(c), respectively

    At low temperatures, the gap is well described by the scaling ∆ ∼ T 1/4 for both the 2D and 3D models, as can be seen in Fig. 2(b) and 2(c), respectively

    work page

  3. [2]

    Remarkably, and maybe counter-intuitively, the largest gap is found right below the phase transition, although the order nearly vanishes in this regime

    At higher temperatures, the gap increases faster than ∆ ∼ T 1/4. Remarkably, and maybe counter-intuitively, the largest gap is found right below the phase transition, although the order nearly vanishes in this regime. This indicates that the gap at QaB is a robust feature of the ordered phase

    work page

  4. [3]

    2(b) and 2(c)

    Finite-size effects are more prominent at low T , where ∆ ∼ T 1/4, as seen in Fig. 2(b) and 2(c). The T 1/4 behavior .— To explain the origin of the exponent 1/4 observed in Fig. 2 at low temperatures, we model the sys- tem’s dynamics by a simple classical oscillation in a quartic potential and in thermal equilibrium, described by the Hamil- tonian H4 = p...

    work page

  5. [4]

    W ¨olfle, Quasiparticles in condensed matter systems, Reports on Progress in Physics 81, 032501 (2018)

    P. W ¨olfle, Quasiparticles in condensed matter systems, Reports on Progress in Physics 81, 032501 (2018)

    work page 2018

  6. [5]

    Kittel, Quantum Theory of Solids , 2nd ed

    C. Kittel, Quantum Theory of Solids , 2nd ed. (Wiley, New York, NY , 1991)

    work page 1991

  7. [6]

    A. L. Fetter, Rotating vortex lattice in a Bose-Einstein conden- sate trapped in combined quadratic and quartic radial potentials, Phys. Rev. A 64, 063608 (2001)

    work page 2001

  8. [7]

    O. Gygi, H. G. Katzgraber, M. Troyer, S. Wessel, and G. G. Ba- trouni, Simulations of ultracold bosonic atoms in optical lattices with anharmonic traps, Phys. Rev. A 73, 063606 (2006)

    work page 2006

Show all 46 references
  1. [8]

    Lakshmanan and R

    M. Lakshmanan and R. Sahadevan, Painlev ´e analysis, Lie sym- metries, and integrability of coupled nonlinear oscillators of polynomial type, Physics Reports 224, 1 (1993)

    1993

  2. [9]

    V . M. Bannur, P. K. Kaw, and J. C. Parikh, Statistical mechanics of quartic oscillators, Phys. Rev. E 55, 2525 (1997)

    1997

  3. [10]

    T. Lan, C. W. Li, O. Hellman, D. S. Kim, J. A. Mu˜noz, H. Smith, D. L. Abernathy, and B. Fultz, Phonon quarticity induced by changes in phonon-tracked hybridization during lattice expan- sion and its stabilization of rutileTiO2, Phys. Rev. B92, 054304 (2015)

    2015

  4. [11]

    Wehinger, A

    B. Wehinger, A. Bosak, and P. T. Jochym, Soft phonon modes in rutile TiO2, Phys. Rev. B 93, 014303 (2016)

    2016

  5. [12]

    Cs ´aki, C.-S

    C. Cs ´aki, C.-S. Guan, T. Ma, and J. Shu, Generating a Higgs Potential Quartic Term, Phys. Rev. Lett.124, 251801 (2020)

    2020

  6. [13]

    Arkani-Hamed, A

    N. Arkani-Hamed, A. G. Cohen, E. Katz, and A. E. Nelson, The Littlest Higgs, Journal of High Energy Physics2002, 034 (2002)

    2002

  7. [14]

    Arkani-Hamed, A

    N. Arkani-Hamed, A. G. Cohen, E. Katz, A. E. Nelson, T. Gre- goire, and J. G. Wacker, The Minimal Moose for a Little Higgs, Journal of High Energy Physics 2002, 021 (2002)

    2002

  8. [15]

    J. G. Rau, P. A. McClarty, and R. Moessner, Pseudo-Goldstone Gaps and Order-by-Quantum Disorder in Frustrated Magnets, Phys. Rev. Lett. 121, 237201 (2018)

    2018

  9. [16]

    Gohlke, L

    M. Gohlke, L. E. Chern, H.-Y . Kee, and Y . B. Kim, Emer- gence of nematic paramagnet via quantum order-by-disorder and pseudo-Goldstone modes in Kitaev magnets, Phys. Rev. Res. 2, 043023 (2020)

    2020

  10. [17]

    Khatua, M

    S. Khatua, M. J. P. Gingras, and J. G. Rau, Pseudo-Goldstone Modes and Dynamical Gap Generation from Order by Thermal 6 Disorder, Phys. Rev. Lett. 130, 266702 (2023)

    2023

  11. [18]

    Hickey, J

    A. Hickey, J. G. Rau, S. Khatua, and M. J. P. Gingras, Uni- versal temperature-dependent power law excitation gaps in frus- trated quantum spin systems harboring order-by-disorder (2025), arXiv:2505.18253 [cond-mat.str-el]

    2025 arXiv

  12. [19]

    Y . V . Tymoshenko, Y . A. Onykiienko, T. M¨uller, R. Thomale, S. Rachel, A. S. Cameron, P. Y . Portnichenko, D. V . Efremov, V . Tsurkan, D. L. Abernathy, J. Ollivier, A. Schneidewind, A. Pi- ovano, V . Felea, A. Loidl, and D. S. Inosov, Pseudo-Goldstone Magnons in the Frustr...

    2017

  13. [21]

    A. S. Sukhanov, Y . A. Onykiienko, R. Bewley, C. Shekhar, C. Felser, and D. S. Inosov, Magnon spectrum of the Weyl semimetal half-Heusler compound GdPtBi, Phys. Rev. B 101, 014417 (2020)

    2020

  14. [22]

    Rastelli, A

    E. Rastelli, A. Tassi, and L. Reatto, Non-simple magnetic order for simple Hamiltonians, Physica B+C 97, 1 (1979)

    1979

  15. [23]

    Below a certain temperature, finite systems appear magnetically ordered when the correlation length becomes comparable to the system size

    The absence of a phase transition in 2D follows from the Mermin-Wagner-Hohenberg theorem [28, 29]. Below a certain temperature, finite systems appear magnetically ordered when the correlation length becomes comparable to the system size. However, true long-range order is absen...

  16. [24]

    Seabra, P

    L. Seabra, P. Sindzingre, T. Momoi, and N. Shannon, Novel phases in a square-lattice frustrated ferromagnet : 1 3 - magnetization plateau, helicoidal spin liquid, and vortex crystal, Phys. Rev. B 93, 085132 (2016)

    2016

  17. [25]

    It also includes Refs

    See Supplemental Material at URL-will-be-inserted-by- publisher for the analytical expansion of the energy of the quartic perturbation, details on the sMD and MC simulations, fitting procedures, additional data for the cubic lattice, and the derivation of the quartic oscillato...

  18. [26]

    By considering systems with spin spiral ground states, we avoid such accidental zero modes

    Accidental continuous ground-state degeneracies are known to occur in systems with fine-tuned interactions that display sim- pler ground-state wave vectors such as Q = (0, ±π) and Q = (±π, 0) [33]. By considering systems with spin spiral ground states, we avoid such accidental...

  19. [27]

    J. D. Alzate-Cardona, D. Sabogal-Su ´arez, R. F. L. Evans, and E. Restrepo-Parra, Optimal phase space sampling for Monte Carlo simulations of Heisenberg spin systems, Journal of Physics: Condensed Matter 31, 095802 (2019)

    2019

  20. [28]

    B. J. Alder and T. E. Wainwright, Phase Transition for a Hard Sphere System, The Journal of Chemical Physics 27, 1208 (1957)

    1957

  21. [29]

    Savary, K

    L. Savary, K. A. Ross, B. D. Gaulin, J. P. C. Ruff, and L. Ba- lents, Order by quantum disorder in Er2Ti2O7, Phys. Rev. Lett. 109, 167201 (2012)

    2012

  22. [30]

    S. M. Rezende, Fundamentals of magnonics , V ol. 969 (Springer, 2020)

    2020

  23. [31]

    N. D. Mermin and H. Wagner, Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic heisen- berg models, Phys. Rev. Lett. 17, 1133 (1966)

    1966

  24. [32]

    P. C. Hohenberg, Existence of long-range order in one and two dimensions, Phys. Rev. 158, 383 (1967)

    1967

  25. [33]

    L. D. Landau and E. Lifshitz, On the theory of the dispersion of magnetic permeability in ferromagnetic bodies, Phys. Z. Sowje- tunion 8, 101 (1935)

    1935

  26. [34]

    Ahnert and M

    K. Ahnert and M. Mulansky, Boost C++ Library: Odeint (2012)

    2012

  27. [35]

    Zhang, H

    S. Zhang, H. J. Changlani, K. W. Plumb, O. Tchernyshyov, and R. Moessner, Dynamical Structure Factor of the Three- Dimensional Quantum Spin Liquid Candidate NaCaNi2F7, Phys. Rev. Lett. 122, 167203 (2019)

    2019

  28. [36]

    C. L. Henley, Ordering due to disorder in a frustrated vector antiferromagnet, Phys. Rev. Lett. 62, 2056 (1989)

  29. [37]

    Holstein and H

    T. Holstein and H. Primakoff, Field Dependence of the Intrinsic Domain Magnetization of a Ferromagnet, Phys. Rev. 58, 1098 (1940)

    1940

  30. [38]

    Anharmonic Collective Oscillations in Isotropic Spin Systems and their Spectroscopic Signatures

    This representation is a classical form of the Holstein-Primakoff representation [34] with unit spin length and ai = ( xi + ipi)/ √ 2. End Matter Section A: Numerical simulation of the T 1/4 behavior .— The quartic perturbation in Eq. (3) is energetically softer than all other...

  31. [39]

    First, it is useful to rewrite the perturbed state given in Eq

    Square lattice To illustrate how quartic modes emerge in a spiral ground state perturbed by a perpendicular perturbation with a wave vector corresponding to the non-selected spiral, QaB, we consider a generic spiral with period L. First, it is useful to rewrite the perturbed s...

  32. [40]

    (S1), generalized to the equal superposition of the two quartic perturbations Sz i ∼ δ[cos(QaB1 · ri) + cos(QaB2 · ri)]

    Cubic lattice To compute the energy expansion, we consider a state similar to the one in Eq. (S1), generalized to the equal superposition of the two quartic perturbations Sz i ∼ δ[cos(QaB1 · ri) + cos(QaB2 · ri)]. We rewrite the state as: Squartic i = 1q 1 + δ2(cos2 ( 2π L x) ...

  33. [41]

    These equations correspond to the Landau-Lifshitz equations without the damping term [S2]

    Overview of the method and numerical details The sMD simulations consist of numerically integrating the classical equations of motion for the spins: dSi dt = Si × hi, (S1) where hi = − ∂H ∂Si is the local effective field acting on Si, arising from interactions with neighboring...

  34. [42]

    To compute S(QaB, ω), we follow the approach of Ref

    Computing the dynamical structure factor The central quantity of our study is the dynamical structure factor S(q, ω) = 1 2πN NX i,j=1 Z ∞ −∞ eiωte−iq(ri−rj )⟨Si(0) · Sj(t)⟩dt (S2) evaluated at q = QaB. To compute S(QaB, ω), we follow the approach of Ref. [S4]. Specifically, th...

  35. [43]

    (S1) The value of the gap corresponds to the position of the maximum, given byx0

    Gap extraction from the dynamical structure factor The value of the gap at q = QaB is obtained by fitting S(QaB, ω) with a Lorentzian function multiplied by a scaling factor A: f (x) = A γ π[γ2 + (x − x0)2] . (S1) The value of the gap corresponds to the position of the maximum...

  36. [44]

    (6) of the main text, which describes the gap behavior derived from the phenomenological model

    Fitting the simulated gaps The gap values are fitted using the function in Eq. (6) of the main text, which describes the gap behavior derived from the phenomenological model. At low temperatures, this function predicts ∆ ∼ T 1/4. This scaling is confirmed by the linear fit of ...

  37. [45]

    Dynamical gap as a function of the perturbation strength Similarly to the square lattice discussed in the End Matter, we computeS(QaB, ω), starting from a configuration in which the quartic modes are excited directly, using the expression in Eq. (S7). The sMD simulations show ...

  38. [46]

    As shown in Fig

    Equivalence of the gap at the anti-Bragg points As explained in the main text, the cubic lattice exhibits quartic modes at the two not-selected spiral wave vectors denoted as QaB1 and QaB2. As shown in Fig. S4, the gap at QaB1 and QaB2 is the same, which agrees with the equiva...

    2019

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.