In frustrated pyrochlore antiferromagnets such as the chromium spinels ZnCr2O4 and MgCr2O4 and the breathing pyrochlore LiGaCr4O8, inelastic neutron scattering shows nearly dispersionless 'molecular-like' excitations coexisting with conventional magnons, a puzzle now two decades old. Shang Gao (Phys. Rev. B 110, 214420 (2024)) showed that a dynamical site-phonon model, in which classical spins couple to harmonic atomic displacements through distance-dependent exchange, generates effective four-spin interactions that give a unified quantitative account of both excitation types, using a phenomenological spin-lattice coupling b of order 0.1-0.2. What remains open is the microscopic origin of these non-bilinear interactions: whether the spin-lattice mechanism, with a first-principles value of b set by the exchange derivative dJ/dr and the elastic constant c, actually operates in these compounds, or whether purely electronic multi-spin exchange (ring exchange from higher-order hopping, or multipolar exchange) produces equivalent four-spin terms. The problem is to determine this origin quantitatively, compound by compound: extract the spin-lattice parameters from first principles and/or experiment, quantify or bound the competing electronic multi-spin mechanisms admissible for the same compounds, decide which mechanism dominates where, and validate the attribution against spectroscopy. The source also proposes generalizing the dynamical site-phonon description to related magnetoelastic frustrated magnets (NiGa2S4, alpha-CaCr2O4); that generalization is not required.

Background

The pyrochlore lattice is a three-dimensional network of corner-sharing tetrahedra. In the chromium spinel magnets ZnCr2O4 and MgCr2O4, magnetic Cr3+ ions (spin S = 3/2) occupy this lattice, and geometric frustration — the lattice cannot accommodate all antiferromagnetic nearest-neighbor bonds on a tetrahedron simultaneously — suppresses conventional order, producing instead a cooperative magnetostructural transition at low temperature. Inelastic neutron scattering (INS) measures the dynamical spin structure factor S(Q, omega), the intensity of spin fluctuations at momentum transfer Q and energy omega. Since the work of Lee et al. (2000, 2002), INS on these compounds has revealed, alongside dispersive spin-wave (magnon) branches, a set of nearly flat, weakly dispersing modes concentrated near zone boundaries. These were interpreted as collective oscillations of small spin clusters carved out of the lattice hexamers (six spins) and heptamers (seven spins) and are called molecular-like excitations, by analogy with the discrete vibrational modes of a free molecule. The puzzle is that free-molecule models treat these flat modes and the magnons as two uncorrelated kinds of dynamics, while experiment shows them coexisting and interacting in one material; moreover, Heisenberg models with fitted further-neighbor exchange reproduce the paramagnetic-phase scattering of MgCr2O4 (Bai et al. 2019) but generate no molecular-like modes in the ordered phase.

A first family of proposed mechanisms is spin-lattice coupling. In bond-phonon (exchange-striction) models the length of each magnetic bond modulates its exchange, which upon static distortion yields effective biquadratic spin interactions and a spin-Peierls-like order-by-distortion transition (Tchernyshyov, Moessner and Sondhi 2002; Penc, Shannon and Shiba 2004, who explained the half-magnetization plateau). These static theories account for the magnetostructural transition but, as Gao (2024) demonstrated quantitatively, not for the molecular-like dynamics. The site-phonon model (Bergman, Shindou, Fiete and Balents 2006; Wang and Vishwanath 2008; Aoyama and Kawamura 2016, 2019) instead attaches a harmonic displacement u_i to every magnetic site:

H=∑⟨ij⟩J ⁣(∣rij(0)+ui−uj∣)Si⋅Sj+c2∑i∣ui∣2,\mathcal{H}=\sum_{\langle ij\rangle}J\!\left(\left|\mathbf{r}^{(0)}_{ij}+\mathbf{u}_{i}-\mathbf{u}_{j}\right|\right)\mathbf{S}_{i}\cdot\mathbf{S}_{j}+\frac{c}{2}\sum_{i}\left|\mathbf{u}_{i}\right|^{2},

where J(r) is a distance-dependent exchange and c an elastic constant. Integrating out the harmonic displacements generates two effective four-spin (non-bilinear) terms,

HSL=−Jˉb∑⟨i,j⟩(Si⋅Sj)2−Jˉb2∑i∑j≠k∈N(i)e^ij⋅e^ik(Si⋅Sj)(Si⋅Sk),\mathcal{H}_{\mathrm{SL}}=-\bar{J}b\sum_{\langle i,j\rangle}\left(\mathbf{S}_{i}\cdot\mathbf{S}_{j}\right)^{2}-\frac{\bar{J}b}{2}\sum_{i}\sum_{j\neq k\in N(i)}\hat{\mathbf{e}}_{ij}\cdot\hat{\mathbf{e}}_{ik}\left(\mathbf{S}_{i}\cdot\mathbf{S}_{j}\right)\left(\mathbf{S}_{i}\cdot\mathbf{S}_{k}\right),

with b=1cJ[dJdr∣r=∣rij(0)∣]2b=\frac{1}{cJ}\left[\left.\frac{dJ}{dr}\right|_{r=\left|\mathbf{r}^{(0)}_{ij}\right|}\right]^{2} a dimensionless spin-lattice coupling; the second term is geometry-dependent, involving the cosines between bond directions e^ij⋅e^ik\hat{\mathbf{e}}_{ij}\cdot\hat{\mathbf{e}}_{ik}. Gao (2024) evolved this Hamiltonian with classical Landau-Lifshitz spin dynamics, treating the four-spin contribution to the molecular field within the random phase approximation (RPA, here the factorization of a four-spin expectation into products of instantaneous two-spin correlations), and showed that with a phenomenological b≈0.2b\approx0.2 for MgCr2O4 and b=0.1b=0.1 (with exchanges J=10.4J = 10.4 meV and J′=6.2J' = 6.2 meV) for the breathing pyrochlore LiGaCr4O8 (a compound where up- and down-pointing tetrahedra carry different exchanges; Okamoto et al. 2013, He et al. 2021), the model reproduces the coexistence of magnons and molecular-like modes on both sides of the ordering temperature. All of these b values were chosen by hand to fit the spectra; no microscopic value of b, dJ/dr, or c has ever been extracted for any of these compounds.

The alternative microscopic sources of four-spin interactions named in the source paper are documented only in other material families: biquadratic exchange arising from multiple hopping of electrons in spin-1 chains (Mila and Zhang 2000) and in iron pnictides (Wysocki, Belashchenko and Antropov 2011); systematic first-principles schemes for deriving multi-spin terms in correlated insulators (Hoffmann and Blugel 2020); spin-lattice-driven trimerization in the S = 1 kagome antiferromagnet Na2Ti3Cl8 (Paul, Chung, Birol and Changlani 2020); and multipolar (higher-order multipole) exchange between spin-orbit-coupled moments in ordered double perovskites (Chen, Pereira and Balents 2010). None has been quantified for the chromium pyrochlores. The most recent first-principles landmark, Jaubert, Iqbal and Jeschke (2025), mapped density-functional-theory total energies onto precise bilinear Heisenberg exchange constants for the cubic, tetragonal and orthorhombic structures of ZnCr2O4 and MgCr2O4, supporting a three-dimensional spin-Peierls mechanism for their magnetostructural transition, but it contains no four-spin, ring-exchange, or b-type couplings and does not address the molecular-like excitations. Whether the four-spin interactions that produce the molecular-like modes come from spin-lattice coupling or from purely electronic multi-spin exchange, and how large each contribution is in each compound, therefore remains undetermined. The source paper also flags two related triangular-lattice magnets where lattice degrees of freedom are known to matter but spin dynamics remain unexplained: NiGa2S4 (Nakatsuji et al. 2005) and alpha-CaCr2O4 (Toth et al. 2012).

Problem Statement

Establish, quantitatively and compound by compound, the microscopic origin of the effective four-spin (non-bilinear) interactions that generate the molecular-like excitations of frustrated pyrochlore magnets. For the compounds in which these excitations are documented — namely the spinels ZnCr2O4 and MgCr2O4 and the breathing pyrochlore LiGaCr4O8 — this requires three jointly supporting results:

(1) Determine the spin-lattice coupling from first principles and/or experiment: the exchange derivative dJdr\frac{dJ}{dr}, the elastic constant cc, and hence the dimensionless coupling b=1cJ[dJdr]2b=\frac{1}{cJ}\left[\frac{dJ}{dr}\right]^{2} entering the four-spin Hamiltonian HSL\mathcal{H}_{\mathrm{SL}} obtained by integrating out the harmonic site displacements (defined in the Background), with uncertainty estimates, and test whether the microscopic bb is consistent with the phenomenological b≈0.1b\approx0.1-0.20.2 needed to reproduce the measured spectra. The values of dJdr\frac{dJ}{dr} and cc must be obtained independently of the spectroscopic fits that produced the phenomenological bb — they may not be reverse-engineered from the measured spectra or from the fitted bb itself — and the consistency test must be carried out with the propagated uncertainties of both the microscopic and the phenomenological value stated explicitly.

(2) Quantify the competing purely electronic four-spin interactions for the same compounds to the extent that their electronic structure admits them: the ring (cyclic) exchange generated by higher-order electron-hopping processes of the Cr3+ (3d3, S=3/2S=3/2) moments, and multipolar exchange terms insofar as the spin-orbit entanglement of those moments makes such terms admissible, including the symmetry structure of the terms that are admitted (for example, whether it reproduces or lacks the bond-direction-cosine geometry e^ij⋅e^ik\hat{\mathbf{e}}_{ij}\cdot\hat{\mathbf{e}}_{ik} of the spin-lattice four-spin term). The spin-lattice mechanism of (1) together with the electronic mechanisms admitted for these compounds in (2) constitute the candidate set over which the accounting below is taken; a mechanism may be left out of that set only on the basis of a stated microscopic or symmetry argument, not by omission.

(3) On this basis, decide which mechanism dominates in which material, and validate the attribution against the magnetic spectroscopy on which the molecular-like excitations are documented, namely the coexistence of dispersive magnons and flat molecular-like modes in S(Q,ω)S(\mathbf{Q},\omega) across the ordering transition; lattice-side signatures such as phonon anomalies and spin-phonon splittings are supporting discriminants where they are available, but are not required components of the validation. This includes deciding whether the dynamical spin-lattice (site-phonon) interpretation of the molecular-like excitations survives quantitatively, and it may require going beyond the classical-spin, random-phase-approximation treatment used so far.

Beyond this completion boundary, the source proposes generalizing the dynamical site-phonon description to a larger family of magnetoelastic frustrated magnets in which the lattice degree of freedom is known to matter but the spin dynamics remain unexplained, notably the triangular-lattice compounds NiGa2S4 and alpha-CaCr2O4; that generalization is an extension of a validated description, not a required deliverable of this problem.

The problem is complete when, for the three compounds named above, a material-resolved quantitative accounting exists in which the contribution of each candidate mechanism of the set defined in (1) and (2) to the effective four-spin interactions is computed or bounded with uncertainties tight enough to rank the mechanisms against one another, and the dominant origin of the molecular-like excitations is thereby established and shown to be consistent with the measured spectra.

References

  1. Shang Gao, "Dynamic Spin-Lattice Coupling and Statistical Interpretation for the Molecular-Like Excitations in Frustrated Pyrochlores", Phys. Rev. B 110, 214420 (2024). DOI: 10.1103/PhysRevB.110.214420; arXiv:2404.05523; https://doi.org/10.1103/PhysRevB.110.214420
  2. S.-H. Lee, C. Broholm, T. H. Kim, W. Ratcliff, S.-W. Cheong, "Local Spin Resonance and Spin-Peierls-like Phase Transition in a Geometrically Frustrated Antiferromagnet", Phys. Rev. Lett. 84, 3718 (2000). DOI: 10.1103/PhysRevLett.84.3718; https://doi.org/10.1103/PhysRevLett.84.3718
  3. S. H. Lee, C. Broholm, W. Ratcliff, G. Gasparovic, Q. Huang, T. H. Kim, S. W. Cheong, "Emergent excitations in a geometrically frustrated magnet", Nature 418, 856 (2002). DOI: 10.1038/nature00964; https://doi.org/10.1038/nature00964
  4. O. Tchernyshyov, R. Moessner, S. L. Sondhi, "Order by Distortion and String Modes in Pyrochlore Antiferromagnets", Phys. Rev. Lett. 88, 067203 (2002). DOI: 10.1103/PhysRevLett.88.067203; https://doi.org/10.1103/PhysRevLett.88.067203
  5. Karlo Penc, Nic Shannon, Hiroyuki Shiba, "Half-Magnetization Plateau Stabilized by Structural Distortion in the Antiferromagnetic Heisenberg Model on a Pyrochlore Lattice", Phys. Rev. Lett. 93, 197203 (2004). DOI: 10.1103/PhysRevLett.93.197203; https://doi.org/10.1103/PhysRevLett.93.197203
  6. D. Bergman, R. Shindou, G. A. Fiete, L. Balents, "Models of degeneracy breaking in pyrochlore antiferromagnets", Phys. Rev. B 74, 134409 (2006). DOI: 10.1103/PhysRevB.74.134409; https://doi.org/10.1103/PhysRevB.74.134409
  7. F. Wang, A. Vishwanath, "Spin Phonon Induced Collinear Order and Magnetization Plateaus in Triangular and Kagome Antiferromagnets: Applications to CuFeO2", Phys. Rev. Lett. 100, 077201 (2008). DOI: 10.1103/PhysRevLett.100.077201; https://doi.org/10.1103/PhysRevLett.100.077201
  8. Kazushi Aoyama, Hikaru Kawamura, "Spin-Lattice-Coupled Order in Heisenberg Antiferromagnets on the Pyrochlore Lattice", Phys. Rev. Lett. 116, 257201 (2016). DOI: 10.1103/PhysRevLett.116.257201; https://doi.org/10.1103/PhysRevLett.116.257201
  9. Kazushi Aoyama, Hikaru Kawamura, "Spin ordering induced by lattice distortions in classical Heisenberg antiferromagnets on the breathing pyrochlore lattice", Phys. Rev. B 99, 144406 (2019). DOI: 10.1103/PhysRevB.99.144406; https://doi.org/10.1103/PhysRevB.99.144406
  10. Yoshihiko Okamoto, G. J. Nilsen, J. Paul Attfield, Zenji Hiroi, "Breathing Pyrochlore Lattice Realized in A-Site Ordered Spinel Oxides LiGaCr4O8 and LiInCr4O8", Phys. Rev. Lett. 110, 097203 (2013). DOI: 10.1103/PhysRevLett.110.097203; https://doi.org/10.1103/PhysRevLett.110.097203
  11. Zheng He, Yiqing Gu, Hongliang Wo, Yu Feng, Die Hu, Yiqing Hao, Yimeng Gu, Helen C. Walker, Devashibhai T. Adroja, Jun Zhao, "Neutron Scattering Studies of the Breathing Pyrochlore Antiferromagnet LiGaCr4O8", Phys. Rev. Lett. 127, 147205 (2021). DOI: 10.1103/PhysRevLett.127.147205; https://doi.org/10.1103/PhysRevLett.127.147205
  12. X. Bai, J. A. M. Paddison, E. Kapit, S. M. Koohpayeh, J.-J. Wen, S. E. Dutton, A. T. Savici, A. I. Kolesnikov, G. E. Granroth, C. L. Broholm, J. T. Chalker, M. Mourigal, "Magnetic Excitations of the Classical Spin Liquid MgCr2O4", Phys. Rev. Lett. 122, 097201 (2019). DOI: 10.1103/PhysRevLett.122.097201; https://doi.org/10.1103/PhysRevLett.122.097201
  13. F. Mila, F.-C. Zhang, "On the origin of biquadratic exchange in spin 1 chains", Eur. Phys. J. B 16, 7 (2000). DOI: 10.1007/s100510070242; https://doi.org/10.1007/s100510070242
  14. Aleksander L. Wysocki, Kirill D. Belashchenko, Vladimir P. Antropov, "Consistent model of magnetism in ferropnictides", Nat. Phys. 7, 485 (2011). DOI: 10.1038/nphys1933; https://doi.org/10.1038/nphys1933
  15. Arpita Paul, Chia-Min Chung, Turan Birol, Hitesh J. Changlani, "Spin-lattice Coupling and the Emergence of the Trimerized Phase in the S=1 Kagome Antiferromagnet Na2Ti3Cl8", Phys. Rev. Lett. 124, 167203 (2020). DOI: 10.1103/PhysRevLett.124.167203; https://doi.org/10.1103/PhysRevLett.124.167203
  16. M. Hoffmann, S. Blugel, "Systematic derivation of realistic spin models for beyond-Heisenberg solids", Phys. Rev. B 101, 024418 (2020). DOI: 10.1103/PhysRevB.101.024418; https://doi.org/10.1103/PhysRevB.101.024418
  17. G. Chen, R. Pereira, L. Balents, "Exotic phases induced by strong spin-orbit coupling in ordered double perovskites", Phys. Rev. B 82, 174440 (2010). DOI: 10.1103/PhysRevB.82.174440; https://doi.org/10.1103/PhysRevB.82.174440
  18. S. Nakatsuji, Y. Nambu, H. Tonomura, O. Sakai, S. Jonas, C. Broholm, H. Tsunetsugu, Y. Qiu, Y. Maeno, "Spin Disorder on a Triangular Lattice", Science 309, 1697 (2005). DOI: 10.1126/science.1114727; https://doi.org/10.1126/science.1114727
  19. S. Toth, B. Lake, K. Hradil, T. Guidi, K. C. Rule, M. B. Stone, A. T. M. N. Islam, "Magnetic Soft Modes in the Distorted Triangular Antiferromagnet alpha-CaCr2O4", Phys. Rev. Lett. 109, 127203 (2012). DOI: 10.1103/PhysRevLett.109.127203; https://doi.org/10.1103/PhysRevLett.109.127203
  20. Ludovic D. C. Jaubert, Yasir Iqbal, Harald O. Jeschke, "Spin-Peierls Transition in the Frustrated Spinels ZnCr2O4 and MgCr2O4", Phys. Rev. Lett. 134, 086702 (2025). DOI: 10.1103/PhysRevLett.134.086702; arXiv:2405.13908; https://doi.org/10.1103/PhysRevLett.134.086702