Recognition: 2 theorem links
· Lean TheoremPhase-slip residual-order spin state in FeSe
Pith reviewed 2026-05-13 02:42 UTC · model grok-4.3
The pith
FeSe is governed by a nearly degenerate manifold of phase-slip defects embedded in a stripe background rather than a single magnetic order.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
FeSe realizes a residual-order spin state in which local stripe antiferromagnetic correlations survive but long-range phase coherence is lost through a manifold of phase-slip defects. Mixed PBE plus r2SCAN calculations show that multiple slip configurations lie within a narrow energy window. Non-local magneto-elastic coupling redistributes domain-wall energy across the lattice while competing magnetic interactions truncate the real-space coherence length to roughly ten moments. Spectrally weighted superpositions of these states recover the experimental line shapes of both stripe- and Neel-type spin fluctuations.
What carries the argument
The residual-order spin state (ROSS) generated by phase-slip defects inside a stripe background, obtained from PBE+r2SCAN energy landscapes and weighted simulations of the static spin structure factor S(q).
If this is right
- The lack of static long-range order in FeSe follows directly from the near-degeneracy of the phase-slip states.
- Superpositions of the slip configurations simultaneously explain both stripe- and Neel-type spin fluctuations seen in neutron scattering.
- Non-local lattice coupling and competing exchange interactions fix the magnetic coherence length at an optimal scale of about ten moments.
- The residual-order spin state supplies a concrete microscopic basis for spin-fluctuation pairing theories of superconductivity.
Where Pith is reading between the lines
- The same phase-slip mechanism may operate in other iron chalcogenides or pnictides that display mixed spin-fluctuation spectra.
- Local probes with spatial resolution below ten lattice spacings could directly image the predicted short coherence length.
- External tuning of the narrow energy window by pressure or doping could systematically alter the superconducting transition temperature.
Load-bearing premise
The PBE and r2SCAN mixed calculations together with spectrally weighted S(q) simulations correctly capture the true energy ordering and momentum-space signatures of the phase-slip manifold without large functional-dependent errors.
What would settle it
A measurement that directly resolves either the real-space coherence length of the magnetic correlations or the small energy splittings between distinct phase-slip configurations in FeSe would confirm or refute the predicted manifold.
Figures
read the original abstract
Clarifying the magnetic ground state is essential for analysing unconventional superconductivity, because microscopic magnetic order provides one of the basic starting assumptions for spin-fluctuation pairing theories. FeSe exhibits pronounced stripe- and Neel-type spin fluctuations yet lacks static long-range order, posing a long-standing puzzle. By combining PBE and r2SCAN mixed exchange-correlation calculations with spectrally weighted simulations of the static spin structure factor S(q), we show that FeSe is not governed by a single magnetic configuration but by a nearly degenerate manifold of phase-slip defects embedded in a stripe background. We term this state a residual-order spin state (ROSS): a spin state that retains local stripe-like antiferromagnetic correlations but loses long-range phase coherence because of phase slips. Multiple slip configurations are compressed into an exceptionally narrow energy window. Non-local magneto-elastic coupling redistributes domain-wall formation energy through the lattice, whereas competing magnetic interactions truncate the real-space coherence length to an optimal scale of about ten moments. Weighted superpositions of these slip states reproduce the momentum-space line shapes of both stripe- and Neel-type spin fluctuations observed by inelastic neutron scattering, providing a microscopic basis for superconductivity models built on spin fluctuations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes that FeSe lacks a single magnetic ground state and is instead described by a 'residual-order spin state' (ROSS): a nearly degenerate manifold of phase-slip defects embedded in a stripe antiferromagnetic background. Mixed PBE+r2SCAN DFT calculations are used to identify multiple slip configurations lying within a narrow energy window above the stripe state; non-local magneto-elastic coupling and competing interactions are invoked to set an optimal real-space coherence length of ~10 moments. Weighted superpositions of the spin structure factors S(q) from these configurations are shown to reproduce the momentum-space line shapes of both stripe- and Néel-type spin fluctuations measured by inelastic neutron scattering.
Significance. If substantiated, the ROSS picture supplies a concrete real-space mechanism for the coexistence of stripe and Néel fluctuations without long-range order in FeSe, directly relevant to spin-fluctuation theories of its superconductivity. The approach of combining mixed exchange-correlation functionals with spectrally weighted S(q) simulations is technically novel and could be applied to other iron-based systems where magnetic degeneracy is suspected.
major comments (4)
- [Abstract / §3] Abstract and §3 (computational details): the reproduction of experimental INS line shapes is achieved via 'spectrally weighted' superpositions, yet the manuscript provides neither the explicit procedure for determining the weights nor a quantitative comparison (e.g., raw vs. fitted spectra with error bars or χ² values). Because the weights are not derived from first principles, the central claim that the phase-slip manifold explains the observed S(q) rests on an adjustable fitting step whose robustness is untested.
- [§2 / §4] §2 and §4 (energy landscape): the assertion of an 'exceptionally narrow energy window' (~few meV) for the phase-slip manifold is obtained exclusively with the PBE+r2SCAN mixed functional. Given the well-documented sensitivity of FeSe stripe vs. Néel ordering energies to exchange-correlation choice, the manuscript should demonstrate that the near-degeneracy survives under alternative functionals or fully self-consistent treatments; otherwise the microscopic basis for the ROSS state may be functional-dependent.
- [§4] §4 (coherence length): the real-space coherence length is stated to be 'truncated to an optimal scale of about ten moments' by competing interactions, yet this scale appears to function as a free parameter in the spectral simulation. Explicit justification of its value, together with a sensitivity analysis showing how S(q) line shapes change when the length is varied by ±2–3 moments, is required to establish that the reproduction is not tuned.
- [§5] §5 (comparison with experiment): while the manuscript states that weighted superpositions 'reproduce' both stripe and Néel peaks, no quantitative metrics (peak positions, widths, relative intensities with uncertainties) or direct overlay of calculated vs. measured S(q) are supplied. This omission prevents assessment of whether the agreement is within experimental resolution or merely qualitative.
minor comments (2)
- [Introduction] The acronym ROSS is introduced in the abstract but its full expansion ('residual-order spin state') should be restated on first use in the main text for clarity.
- [Figure captions] Figure captions for the S(q) plots should include the precise coherence length, weight values, and any broadening parameters used in the simulation so that the results are reproducible from the text alone.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable suggestions. We have carefully considered each comment and revised the manuscript to improve clarity, provide additional details, and strengthen the presentation of our results. Below we respond point by point.
read point-by-point responses
-
Referee: [Abstract / §3] Abstract and §3 (computational details): the reproduction of experimental INS line shapes is achieved via 'spectrally weighted' superpositions, yet the manuscript provides neither the explicit procedure for determining the weights nor a quantitative comparison (e.g., raw vs. fitted spectra with error bars or χ² values). Because the weights are not derived from first principles, the central claim that the phase-slip manifold explains the observed S(q) rests on an adjustable fitting step whose robustness is untested.
Authors: We agree that the explicit procedure for determining the weights and a quantitative comparison of the fits should be included. In the revised manuscript, we detail the least-squares minimization procedure used to obtain the weights from the experimental INS line shapes. We also provide direct comparisons of the simulated and measured S(q) spectra, including χ² values and error bars to assess the quality of the reproduction. revision: yes
-
Referee: [§2 / §4] §2 and §4 (energy landscape): the assertion of an 'exceptionally narrow energy window' (~few meV) for the phase-slip manifold is obtained exclusively with the PBE+r2SCAN mixed functional. Given the well-documented sensitivity of FeSe stripe vs. Néel ordering energies to exchange-correlation choice, the manuscript should demonstrate that the near-degeneracy survives under alternative functionals or fully self-consistent treatments; otherwise the microscopic basis for the ROSS state may be functional-dependent.
Authors: We recognize the importance of checking functional dependence given the known sensitivities in FeSe. The mixed PBE+r2SCAN approach was employed as it provides improved accuracy for the magneto-elastic effects central to our model. We have added calculations using the standard PBE functional in the revised §4, which show that the phase-slip configurations remain closely spaced in energy, although the precise window width varies. This supports the robustness of the ROSS picture. revision: partial
-
Referee: [§4] §4 (coherence length): the real-space coherence length is stated to be 'truncated to an optimal scale of about ten moments' by competing interactions, yet this scale appears to function as a free parameter in the spectral simulation. Explicit justification of its value, together with a sensitivity analysis showing how S(q) line shapes change when the length is varied by ±2–3 moments, is required to establish that the reproduction is not tuned.
Authors: The coherence length is physically motivated by the balance of non-local magneto-elastic coupling and competing magnetic interactions, as explained in §4. To address the concern, we have included in the revised manuscript a sensitivity analysis where the coherence length is varied by ±3 moments. The resulting S(q) line shapes remain consistent with the experimental data within the reported resolution, confirming that the choice is not finely tuned. revision: yes
-
Referee: [§5] §5 (comparison with experiment): while the manuscript states that weighted superpositions 'reproduce' both stripe and Néel peaks, no quantitative metrics (peak positions, widths, relative intensities with uncertainties) or direct overlay of calculated vs. measured S(q) are supplied. This omission prevents assessment of whether the agreement is within experimental resolution or merely qualitative.
Authors: We agree that quantitative metrics enhance the comparison. The revised manuscript now includes direct overlays of the calculated and experimental S(q) along the relevant momentum directions, as well as tabulated values for peak positions, widths, and relative intensities with associated uncertainties. These additions allow for a direct evaluation of the agreement quality. revision: yes
Circularity Check
Reproduction of INS line shapes depends on data-fitted weights for phase-slip superpositions
specific steps
-
fitted input called prediction
[Abstract]
"Weighted superpositions of these slip states reproduce the momentum-space line shapes of both stripe- and Neel-type spin fluctuations observed by inelastic neutron scattering, providing a microscopic basis for superconductivity models built on spin fluctuations."
The weights in the superpositions are not computed from the PBE+r2SCAN energies or any first-principles procedure but are selected to match the observed INS line shapes; the reproduction is therefore achieved by construction through adjustment of the weights rather than constituting an independent prediction.
full rationale
The DFT-derived near-degeneracy of phase-slip states is an independent computational result. However, the key claim that weighted superpositions explain both stripe and Néel fluctuations reduces to choosing weights to match experimental S(q) data rather than deriving them from the energies or first principles. This is a moderate instance of fitted input called prediction. No self-definitional, self-citation load-bearing, or other patterns are present in the provided text.
Axiom & Free-Parameter Ledger
free parameters (1)
- coherence length scale
axioms (1)
- domain assumption PBE and r2SCAN mixed functionals accurately rank the energies of stripe and phase-slip configurations in FeSe
invented entities (1)
-
residual-order spin state (ROSS)
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
nearly degenerate manifold of phase-slip defects... weighted superpositions... reproduce... stripe- and Néel-type spin fluctuations
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
PBE and r2SCAN mixed exchange-correlation calculations
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. Dahm et al., Strength of the spin-fluctuation- mediatedpairinginteractioninahigh-temperaturesu- perconductor, Nature Physics 5, 217–221 (2009)
work page 2009
-
[2]
P. Hlobil et al., Tracing the electronic pairing glue in unconventionalsuperconductorsviainelasticscanning tunneling spectroscopy,Phys. Rev. Lett. 118, 167001 (2017)
work page 2017
-
[3]
D. J. Scalapino, A common thread: The pairing inter- action for unconventional superconductors,Rev. Mod. Phys. 84, 1383 (2012)
work page 2012
-
[4]
P.Monthoux,A.V.Balatsky,D.Pines,Self-consistent pairing in a two-dimensional Hubbard model,Phys. Rev. Lett. 72, 1874 (1994)
work page 1994
-
[5]
K. Kuroki et al., Unconventional pairing originating from the disconnected Fermi surfaces of supercon- ducting 𝐿𝑎𝐹𝑒 𝐴𝑠𝑂 1− 𝑥 𝐹𝑥, Phys. Rev. Lett.101,087004 (2008)
work page 2008
-
[6]
I. I. Mazin et al., Unconventional superconductiv- ity with a sign reversal in the order parameter of 𝐿𝑎𝐹𝑒 𝐴𝑠𝑂 1− 𝑥 𝐹𝑥, Phys. Rev. Lett.101,057003(2008)
work page 2008
-
[7]
K.Kurokietal.,Pnictogenheightasapossibleswitch between high-𝑇𝑐 nodeless and low-𝑇𝑐 nodal pairings in the iron-based superconductors,Phys. Rev. B 79, 224511 (2009)
work page 2009
-
[8]
R.Liuetal.,Spincorrelationsinthenematicquantum disordered state of FeSe,Nature Communications 16, 5212 (2025)
work page 2025
-
[9]
W. Li et al., Phase separation and magnetic order in K-dopedironselenidesuperconductor, Nature Physics 8, 126–130 (2012)
work page 2012
-
[10]
S. H. Pan et al., Microscopic electronic inhomogene- ity in the high-𝑇𝑐 superconductor 𝐵𝑖2𝑆𝑟2𝐶𝑎𝐶𝑢 2𝑂8+𝑥, Nature 413, 282–285 (2001)
work page 2001
-
[11]
Q.Wangetal.,MagneticgroundstateofFeSe, Nature Communications 7, 12182 (2016)
work page 2016
-
[12]
Q. Wang et al., Strong interplay between stripe spin fluctuations, nematicity and superconductivity in FeSe, Nature Materials 15, 159–163 (2016)
work page 2016
-
[13]
Chen et al., Anisotropic spin fluctuations in de- twinned FeSe,Nature Materials 18, 709–716 (2019)
T. Chen et al., Anisotropic spin fluctuations in de- twinned FeSe,Nature Materials 18, 709–716 (2019)
work page 2019
-
[14]
Baek et al., Orbital-driven nematicity in FeSe, Nature Materials 14, 210–214 (2015)
S.-H. Baek et al., Orbital-driven nematicity in FeSe, Nature Materials 14, 210–214 (2015)
work page 2015
-
[15]
J. K. Glasbrenner et al., Effect of magnetic frustration on nematicity and superconductivity in iron chalco- genides, Nature Physics 11, 953–958 (2015)
work page 2015
-
[16]
K. Kothapalli et al., Strong cooperative coupling of pressure-induced magnetic order and nematicity in FeSe, Nature Communications 7, 12728 (2016)
work page 2016
-
[17]
J.-P.Sun etal., Dome-shaped magneticorder compet- ing with high-temperature superconductivity at high pressures in FeSe,Nature Communications 7, 12146 (2016)
work page 2016
-
[18]
P.S.Wangetal.,Pressure-inducedstripe-orderantifer- romagnetism and first-order phase transition in FeSe, Phys. Rev. Lett. 117, 237001 (2016)
work page 2016
-
[19]
K. Matsuura etal., Maximizing𝑇𝑐 by tuning nematic- ity and magnetism in 𝐹𝑒𝑆𝑒 1− 𝑥 𝑆𝑥 superconductors, 7 Nature Communications 8, 1143 (2017)
work page 2017
-
[20]
Bendele et al., Pressure induced static magnetic order in superconducting 𝐹𝑒𝑆𝑒 1− 𝑥, Phys
M. Bendele et al., Pressure induced static magnetic order in superconducting 𝐹𝑒𝑆𝑒 1− 𝑥, Phys. Rev. Lett. 104, 087003 (2010)
work page 2010
-
[21]
Bendele et al., Coexistence of superconductivity andmagnetismin 𝐹𝑒𝑆𝑒 1− 𝑥 underpressure, Phys
M. Bendele et al., Coexistence of superconductivity andmagnetismin 𝐹𝑒𝑆𝑒 1− 𝑥 underpressure, Phys. Rev. B 85, 064517 (2012)
work page 2012
-
[22]
R. Zhou et al., Singular magnetic anisotropy in the nematic phase of FeSe,npj Quantum Materials 5, 43 (2020)
work page 2020
-
[23]
Cao et al., Antiferromagnetic ground state with pair-checkerboard order in FeSe, Phys
H.-Y. Cao et al., Antiferromagnetic ground state with pair-checkerboard order in FeSe, Phys. Rev. B 91, 020504(R) (2015)
work page 2015
-
[24]
K.Liu,Z.-Y.Lu,T.Xiang,Nematicantiferromagnetic states in bulk FeSe,Phys. Rev. B 93, 205154 (2016)
work page 2016
-
[25]
K.Liu,B.-J.Zhang,Z.-Y.Lu,First-principlesstudyof magneticfrustrationinFeSeepitaxialfilmsonSrTiO 3, Phys. Rev. B 91, 045107 (2015)
work page 2015
-
[26]
A. Subedi et al., Density functional study of FeS, FeSe, and FeTe: electronic structure, magnetism, phonons, and superconductivity, Phys. Rev. B 78, 134514 (2008)
work page 2008
-
[27]
L. A. Myers et al., Stripe antiferromagnetic ground- state configuration of FeSe revealed by density func- tional theory,Phys. Rev. B 111, L121109 (2025)
work page 2025
-
[28]
P.O.Sprauetal.,Discoveryoforbital-selectiveCooper pairing in FeSe,Science 357, 75–80 (2017)
work page 2017
-
[29]
T. Hashimoto et al., Superconducting gap anisotropy sensitivetonematicdomainsinFeSe, Nature Commu- nications 9, 282 (2018)
work page 2018
-
[30]
J. Kang and R. M. Fernandes, Superconductivity in FeSe: the role of nematic order,Phys. Rev. Lett. 120, 267001 (2018)
work page 2018
-
[31]
R. M. Fernandes and A. J. Millis, Suppression of su- perconductivitybyNéel-typemagneticfluctuationsin theironpnictides, Phys. Rev. Lett.111,127001(2013)
work page 2013
-
[32]
H. C. Xu et al., Highly anisotropic and twofold sym- metric superconducting gap in nematically ordered 𝐹𝑒𝑆𝑒 0.93𝑆0.07, Phys. Rev. Lett. 117, 157003 (2016)
work page 2016
-
[33]
G. Kresse and J. Furthm"uller, Efficient iterative schemes for ab initio total-energy calculations using aplane-wavebasisset, Phys. Rev. B 54,11169(1996)
work page 1996
-
[34]
G. Kresse and J. Furthm"uller, Efficiency of ab-initio total energy calculations for metals and semiconduc- tors using a plane-wave basis set,Comput. Mater. Sci. 6, 15–50 (1996)
work page 1996
-
[35]
P. E. Bl"ochl, Projector augmented-wave method, Phys. Rev. B 50, 17953 (1994)
work page 1994
-
[36]
G.KresseandD.Joubert,Fromultrasoftpseudopoten- tials to the projector augmented-wave method,Phys. Rev. B 59, 1758 (1999)
work page 1999
-
[37]
J.W.Furnessetal.,Accurateandnumericallyefficient 𝑟2SCAN meta-generalized gradient approximation,J. Phys. Chem. Lett. 11, 8208 (2020)
work page 2020
-
[38]
J. P. Perdew, K. Burke, M. Ernzerhof, Generalized gradient approximation made simple,Phys. Rev. Lett. 77, 3865 (1996)
work page 1996
-
[39]
J. Sun, A. Ruzsinszky, J. P. Perdew, Strongly con- strained and appropriately normed semilocal density functional, Phys. Rev. Lett. 115, 036402 (2015)
work page 2015
-
[40]
M. A. L. Marques et al., Libxc: a library of exchange and correlation functionals for density functional the- ory, Comput. Phys. Commun. 183, 2272 (2012)
work page 2012
-
[41]
S. Lehtola et al., Recent developments in Libxc—A comprehensivelibraryoffunctionalsfordensityfunc- tional theory,SoftwareX 7, 1–5 (2018)
work page 2018
-
[42]
J. Pelliciari et al., Evolution of spin excitations from bulk to monolayer FeSe,Nature Communications 12, 3122 (2021)
work page 2021
-
[43]
Sun et al., Structural-transition-induced quasi two- dimensional Fermi surface in FeSe,Phys
Y. Sun et al., Structural-transition-induced quasi two- dimensional Fermi surface in FeSe,Phys. Rev. B 94, 134505 (2016)
work page 2016
-
[44]
Wang et al., Quasi-one-dimensional spin excita- tions in the iron pnictide 𝑁𝑎𝐹𝑒 0.53𝐶𝑢 0.47 𝐴𝑠, Phys
Y. Wang et al., Quasi-one-dimensional spin excita- tions in the iron pnictide 𝑁𝑎𝐹𝑒 0.53𝐶𝑢 0.47 𝐴𝑠, Phys. Rev. Lett. 136, 066503 (2026)
work page 2026
-
[45]
K. L. Seyler et al., High-speed antiferromagnetic do- main walls driven by coherent spin waves, Nature Communications 16, 9836 (2025)
work page 2025
-
[46]
A. Vindigni et al., Stripe width and non-local domain walls in the two-dimensional dipolar frustrated Ising ferromagnet, Phys. Rev. B 77, 092414 (2008)
work page 2008
-
[47]
V.Jamei,S.A.Kivelson,B.Spivak,Universalaspects of Coulomb-frustrated phase separation, Phys. Rev. Lett. 94, 056805 (2005)
work page 2005
-
[48]
Dai, Antiferromagnetic order and spin dynamics in iron-based superconductors,Rev
P. Dai, Antiferromagnetic order and spin dynamics in iron-based superconductors,Rev. Mod. Phys. 87, 855 (2015)
work page 2015
-
[49]
R.Coldeaetal.,Spinwavesandelectronicinteractions in 𝐿𝑎 2𝐶𝑢𝑂 4, Phys. Rev. Lett. 86, 5377 (2001)
work page 2001
-
[50]
Kreisel et al., Spin excitations in a model of FeSe withorbitalordering, Phys
A. Kreisel et al., Spin excitations in a model of FeSe withorbitalordering, Phys. Rev. B 92,224515(2015)
work page 2015
-
[51]
A. Kreisel et al., Itinerant approach to magnetic neu- tron scattering of FeSe: effect of orbital selectivity, Phys. Rev. B 98, 214518 (2018). 8
work page 2018
-
[52]
M. Ma et al., Prominent role of spin-orbit coupling in FeSe revealed by inelastic neutron scattering,Phys. Rev. X 7, 021025 (2017)
work page 2017
-
[53]
J. A. M. Paddison and A. L. Goodwin, Empirical magneticstructuresolutionoffrustratedspinsystems, Phys. Rev. Lett. 108, 017204 (2012)
work page 2012
-
[54]
E. Granado, J. W. Lynn, R. F. Jardim, M. S. Torikachvili, Two-dimensional magnetic correlations andpartiallong-rangeorderingeometricallyfrustrated 𝑆𝑟2𝑌 𝑅𝑢𝑂6, Phys. Rev. Lett. 110, 017202 (2013)
work page 2013
-
[55]
B. A. Frandsen et al., Verification of Anderson su- perexchange in MnO via magnetic pair distribution function analysis and ab initio theory,Phys. Rev. Lett. 116, 197204 (2016). Methods. First-principles calculations were performed using the Vi- enna ab initio Simulation Package (VASP 6.5.0) within the collinear spin-polarized framework. A plane-wave cuto...
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.