Spin-nematic order induced superconductivity

Naoum Karchev

REVIEW 4 major objections 5 minor 25 references

Spin-nematic order induced superconductivity

T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read Spin-nematic (quadrupolar) order turns the fermion-quadrupolar interaction into a four-fermion pairing interaction that yields p-wave superconductivity; the sign of the coupling fixes equal-spin versus opposite-spin pairing, and the critical coupling vanishes at zero chemical potential.

desk verdict Novel mechanism, but the central κcr(μ) result is an artifact of a wrong single-particle dispersion; the idea deserves peer review, the phase diagram does not. read the letter →

arxiv 1908.02514 v2 pith:YOGT2DPS submitted 2019-08-07 cond-mat.supr-con

Naoum Karchev This is my paper · ORCID

classification cond-mat.supr-con

keywords superconductivitychemicalcouplingfermionsinteractionphasepotentialwhen

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

Most theories of unconventional superconductivity explain pairing through magnetic fluctuations: electrons exchange magnetic ripples and pair up. This paper starts from a different place. It considers a crystal in which electron spins form a spin-nematic pattern: the spins are stretched along preferred axes, like tiny rods, but produce no net magnetization. This state, known as quadrupolar or nematic order, is documented in the iron-based superconductor FeSe.

The author writes a model in which conduction electrons interact with this nematic order through an interaction built from the same mathematical structures used for the spin order. Inside the nematic phase, that interaction becomes a direct pairing force between electrons, acting in the odd-parity 'p-wave' channel. Solving the resulting mean-field equations gives two regimes. For a positive coupling constant, the paired electrons have parallel spins; for a negative one, opposite spins. At zero chemical potential, which corresponds to half filling, any nonzero coupling, however small, creates a superconducting state; this is the familiar Cooper instability. Away from half filling, the model demands a threshold: the coupling must exceed a critical value that grows with the chemical potential, shown in Fig. 2.

The calculation is analytic plus the numerical solution of gap equations on a cubic lattice. The paper does not compare quantitatively with FeSe: no transition temperature, gap size, or coupling estimate is derived from experiment. The nematic order is treated as a rigid background, and its fluctuations, believed to be important in FeSe, are discarded.

Extended reading notes

Core claim

The central assertion: in the nematic phase, the fermion-quadrupolar interaction (16) reduces to the four-fermion pairing interaction (18), and the mean-field gap equations (25)-(26) then yield p-wave superconductivity with a spin structure selected by the sign of the coupling: 'When the coupling constant is positive the superconductivity is p-wave with spin-parallel paired fermions. When it is negative the superconductivity is p-wave and fermions are spin-antiparallel paired. For a system with zero chemical potential, even a very small coupling can bind fermions into bound state that leads to the superconductivity. When the chemical potential is non-zero the system possesses quantum critical transition from normal spin-nematic phase to phase where superconductivity coexists with spin-nematicity' (Abstract; Eqs. 22-27; Figs. 1-2). If correct, the model shows that static quadrupolar order alone can produce odd-parity pairing with a critical coupling that grows with chemical potential.

Load-bearing premise

The nematic order parameter is imported from the classical limit of the spin-only Hamiltonian and is treated as a rigid background: <Sxx>=<Syy>=1 is inserted in place of the operator in (16)-(17), and the fermions never feed back into the nematic order, while all dynamical quadrupolar fluctuations are discarded. If the fermion subsystem substantially renormalizes the nematic order, or if the pairing glue is the quadrupolar fluctuation propagator rather than the static order parameter, the central result collapses or changes channel. Location: Section 'Fermion-Quadrupolar Interaction', Eqs. (9)-(10) and (16)-(18); the experimental motivation (FeSe with strong nematic fluctuations) makes the neglect of fluctuations especially fragile.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Colleague—

Two things to know. The paper's core construction is new and worth a look: it builds a fermion-quadrupolar interaction out of the triplet-pairing d-vector, then shows that in a nematic background this reduces to a four-fermion pairing term with a sign-selected spin structure. The sign rule—κ>0 equal-spin, κ<0 opposite-spin—is concrete and testable. That part is good. But the paper's central phase diagram claim, the finite μ-dependent critical coupling, is wrong. It comes from a dispersion error.

Look at Eqs. (25)-(27). The paper defines ε_k = -2t(coskx+cosky+coskz), so the single-particle energy is ξ_k = ε_k - μ. The gap equations instead use (-2t ε_k - μ). That double-counts t and is dimensionally off. With the correct ξ_k, the linearized gap equation is 1 = |κ| ∫ e_k^2 / |ξ_k|. For any μ inside the band, the Fermi surface has nonzero measure and e_k = sinkx+sinky+sinkz does not vanish identically on it, so the integral diverges logarithmically. Cooper's theorem then gives κcr=0 for every in-band μ, not just μ=0. The finite κcr that grows with μ in Fig. 2 is almost certainly a numerical artifact of an unregularized divergent integral or the wrong dispersion. That removes the advertised quantum critical transition from the normal nematic to the coexisting superconducting-nematic phase.

The reader's other concerns are real but secondary. The no-coexistence claim after Eq. (25) is asserted, not proved. The numerics lack mesh and regulator details. And treating the nematic order as a rigid background ignores fluctuations, which is a strong assumption for FeSe. But those could be fixed or tempered. The dispersion error is load-bearing and needs a full re-derivation.

Give the paper credit for not fitting parameters: κcr(μ) is computed, not tuned. And the reduction from (16) to (18) is genuinely clever. But as written, the main quantitative result does not survive.

Who is this for? Someone interested in nematicity-induced pairing mechanisms could mine the construction and the sign rule. The paper deserves a serious referee, not because the conclusions are right, but because the idea is novel and the error is identifiable and fixable. I'd send it out, with a directive to redo the gap equations with the correct dispersion and either prove or drop the coexistence claim.

Would I cite it? Not as it stands. Would I bring it to reading group? Maybe, as a cautionary tale about dispersion mistakes.

Best, [Name]

Referee Report

4 major / 5 minor

Summary. The manuscript studies a spin-fermion model on a cubic lattice, combining a spin-1 Heisenberg/biquadratic Hamiltonian, a tight-binding fermion band, and a fermion-quadrupolar interaction. In the ferroquadrupolar/nematic phase the quadrupolar operator is replaced by its classical expectation value, reducing the interaction to a four-fermion pairing interaction. A Hartree-Fock treatment with a T1u gap ansatz yields spin-triplet p-wave superconductivity, with the sign of the coupling selecting equal-spin versus opposite-spin pairing. The paper claims that at zero chemical potential an arbitrarily small coupling produces superconductivity, while for nonzero chemical potential there is a quantum critical coupling separating the normal nematic phase from a phase where superconductivity coexists with nematicity, and it connects the low-pressure regime of FeSe to this scenario.

Significance. The proposed mechanism—static quadrupolar order generating odd-parity pairing—is conceptually interesting and the paper makes a specific, falsifiable prediction for the chemical-potential dependence of the critical coupling. The model is simple and the T1u channel selection is clearly identified, and the paper does not fit any parameter to data. However, the central quantitative prediction is not reliable because the gap equations use a single-particle dispersion inconsistent with the Hamiltonian, and the advertised finite critical coupling for in-band chemical potentials appears to be an artifact of that inconsistency.

major comments (4)

[Superconductivity, Eqs. (19)-(27)] The single-particle dispersion in the gap equations is inconsistent with the Hamiltonian. Equation (19) defines ε_k = -2t(coskx+cosky+coskz), so the quasiparticle energy should contain ξ_k = ε_k - μ. Equation (23) omits μ entirely, and Eqs. (25)-(27) use (-2t ε_k - μ), which is neither ξ_k nor any dimensionally consistent expression. All subsequent results, including Figs. 1 and 2, need to be recomputed with the correct dispersion.
[Eq. (27) and Fig. 2] With the correct dispersion, the linearized gap equation is 1 = |κ| ∫ d^3k/(2π)^3 (sin kx + sin ky + sin kz)^2 / |ε_k - μ|. For any μ inside the band (-6t < μ < 6t), the Fermi surface has nonzero measure, and the integrand diverges logarithmically near it, giving κcr = 0 rather than a finite value. The finite κcr(μ) shown in Fig. 2, and the quantum critical transition claimed in the abstract and Conclusion, are therefore not supported by the corrected equation. This is a load-bearing error in the central claim.
[Fermion-Quadrupolar Interaction, Eqs. (16)-(18)] The reduction of the fermion-quadrupolar interaction to the four-fermion interaction (18) is not shown; since this is the basis of the pairing channel, the algebra should be provided. More importantly, the replacement of S^ba by its classical expectation value in Eq. (17) is made once and never updated: the fermions have no feedback on the nematic order, and dynamical quadrupolar fluctuations are discarded. Given that the paper claims coexistence of nematicity and superconductivity and motivates the model with FeSe, where nematic fluctuations are strong, the rigidity of the nematic background is a significant limitation that should be justified or relaxed.
[Figs. 1 and 2] No numerical grid, momentum cutoff, or regularization is specified for the integrals in Eqs. (25)-(27). Because the κcr integral is divergent at the Fermi surface for in-band μ, any finite value obtained from it is regularization-dependent, and the results in Figs. 1 and 2 cannot be reproduced as presented.

minor comments (5)

[After Eq. (24)] The statement that the system (25) has no solution with both Δp and Δap nonzero is asserted without proof; either provide a derivation or a reference.
[References] Reference [5] is misformatted ("Aarts Steglich" should be Steglich et al.), and several other references contain typographical or formatting errors.
[Fig. 1 caption] The caption is ambiguous: "upper graph/lower one" combined with two colors makes it difficult to identify which curve corresponds to which sign of the coupling and which value of μ; please label the curves directly.
[Eq. (20)] The identity <c_{k↑}c_{-k↑}> = <c_{k↓}c_{-k↓}> is used without justification; in a nematic background that breaks spin-rotation symmetry, this equality should be verified explicitly.
[Conclusion] The statement that the low-pressure region of FeSe is "well described by the theory" is qualitative only; no direct comparison with the measured Tc or pressure dependence is provided.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the model's effective pairing interaction follows algebraically from the assumed spin-fermion coupling, and the gap/κcr results are computed, not fitted.

full rationale

The paper makes no empirical fits and contains no author self-citations. Its central chain is: given a spin-fermion model Hfq constructed from fermion bilinears D (Eq. 13) and quadrupolar operators, replacing the spin-quadrupolar operator by its nematic expectation value (10) yields an explicit four-fermion interaction (18); Hartree-Fock decoupling and a T1u gap ansatz (21) then yield gap equations (25)-(26) whose solutions and κcr(μ) are computed rather than tuned. The fact that D in (13) is already a spin-triplet pairing operator means the model is designed to have a pairing channel, but that is a modeling choice, not a circular reduction: the superconducting order is not inserted as an input; it emerges as a nonzero self-consistent solution only for the appropriate sign of κ. The finite-μ critical coupling is obtained from Eq. (27), a parameter-free integral, so no fitted input is renamed as a prediction. Concerns about the rigid-nematic mean-field treatment or the dispersion used in the gap equations are physics/correctness issues, not circularity.

Assumptions & free parameters 3 free parameters · 6 assumptions · 1 invented entities

The central claim rests on the invented fermion-quadrupolar coupling, a rigid classical nematic background, and two hand-chosen ansatze (the T1u form factor and the equal-spin identity). No parameter is fitted to experimental data: the ledger is a list of assumptions and ansatze, not of fitted values.

free parameters (3)

T1u gap form factor (Eq. 21) = Delta_p or Delta_ap times (sin kx + sin ky + sin kz)
The odd-parity gap is assumed to have the T1u configuration, chosen by hand after citing ref. [14]; no derivation rules out other odd-parity channels on the cubic lattice.
Equal-spin gap identity (Eq. 20) = <ck_up c_-k_up> = <ck_down c_-k_down>
Fixes the two equal-spin gap amplitudes equal; justified only by the zero-magnetization symmetry of the nematic state and stated without derivation. It shapes the sign-selection result.
Model couplings kappa and mu (scanned) = kappa/t in [-10,10]; mu/t in [0,5]
Model inputs scanned to produce the phase diagram; they are not fitted to any experimental data, which keeps the circularity burden low but also prevents any quantitative comparison with FeSe.

assumptions (6)

domain assumption Semiclassical (large-S) analysis of the bilinear-biquadratic spin Hamiltonian (1) gives a ferroquadrupolar ground state for 5pi/4 < gamma < 3pi/2 with <Sxx>=<Syy>=1 (Eqs. 9-10).
Borrowed from refs. [12,13]; the paper does not rederive this and relies on it to set the whole nematic background of the calculation.
ad hoc to paper The quadrupolar operator S^ba is replaced by its classical expectation value <S^ba> and is never updated with the fermions present (replacement in Eq. 17).
The fermions see a rigid order parameter; back-reaction and all nematic fluctuations are neglected, even though the paper's own motivation (FeSe) features strong nematic fluctuations.
domain assumption Identity (20): <c_{k up} c_{-k up}> = <c_{k down} c_{-k down}>.
Follows from time-reversal symmetry in the zero-magnetization nematic state but is stated without derivation; it forces the two equal-spin gap amplitudes to be equal.
standard math Hartree-Fock (mean-field) decoupling of the four-fermion interaction (18) into the quasiparticle Hamiltonian (19).
Standard BCS procedure; the paper does not quantify its validity for the large couplings considered (kappa/t up to 10).
domain assumption One-band model: only the dxy electrons of iron are retained.
The crystal-field splitting argument in the text is qualitative; a multi-orbital treatment could change the pairing channel and the critical coupling.
ad hoc to paper The fermion vector D^a (13) couples to the quadrupolar tensor exactly like the spin-1 Schwinger-boson operators (15)-(16).
This is the invented interaction at the heart of the model; no microscopic derivation from a realistic FeSe Hamiltonian is offered, and the signs and coefficients in (18) depend on it.

invented entities (1)

Fermion-quadrupolar interaction Hfq (Eq. 16)
purpose: Couples the triplet-pairing d-vector D^a to the spin-quadrupolar tensor; its condensation in the nematic phase produces the four-fermion pairing interaction (18).
The interaction is postulated by analogy with the Schwinger-boson representation of spin-1 operators; no independent experimental or microscopic handle fixes its form or strength, so it is an assumption of the model rather than a derived or measured quantity.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Spin-nematic order induced superconductivity." pith.science (2026). https://pith.science/paper/YOGT2DPS

@misc{pith2026190802514,
  author       = {Pith},
  title        = {Pith review of: Spin-nematic order induced superconductivity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YOGT2DPS}},
  note         = {Machine review of arXiv:1908.02514}
}

read the original abstract

We explore a spin-fermion model with fermion-spin-quadrupolar interaction. In a nematic phase, this interaction reduces to a four-fermion interaction that is a basis of superconductivity. When the coupling constant is positive the superconductivity is p-wave with spin-parallel paired fermions. When it is negative the superconductivity is p-wave and fermions are spin-antiparallel paired. For a system with zero chemical potential, even a very small coupling can bind fermions into bound state that leads to the superconductivity. When the chemical potential is non-zero the system possesses quantum critical transition from normal spin-nematic phase to phase where superconductivity coexists with spin-nematicity. The value of the quantum critical fermion-spin-nematicity coupling constant depends on the chemical potential.

Figures

Figures reproduced from arXiv: 1908.02514 by the authors.

**Figure 2.** (Color online) The critical constant κcr/t as a function of the chemical potential µ/t Conclusion. – The appearance of unconventional superconductivity in Fe-based systems is commonly thought to arise from a spin fluctuation pairing mechanism [15–18]. However, the unique properties of nonmagnetic F eSe provides a case to test this view and gives an opportunity for theoretical investigation of a new mechanism of su… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages

[1]

V. L. Ginzburg, Sov.Phys.JETP, 4 (1957) 153

work page 1957
[2]

Saxena, P

S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. Haselwimmer, M. Steiner, E. Pugh, I. Walker, S. Julian, P. Monthoux, G. Lonzarich, A. Huxley, I. Sheikin, D. Braithwaite, and J. Flouquet , Nature (London), 406 (2000) 587

work page 2000
[3]

D. Aoki, A. Huxley, E. Ressouche, D. Braithwaite, J. Flouquet, J-P. Brison, E.Lhotel, and C. Paulsen, Nature (London), 413 (1974) 613

work page 1974
[4]

J. G. Bednorz and K. A. M¨uller, Z. Phys, B 64 (1986) 189

work page 1986
[5]

Aarts Steglich, C

J. Aarts Steglich, C. D. Bredl, W. Lieke, D. Meschede, W. Franz, and H. Schafer, Phys. Rev. Lett, 43 (1979) 1892

work page 1979
[6]

Hegger C

H. Hegger C. Petrovic, E. G. Moshopoulou, M. F. Hundley, J. L. Sarrao, Z. Fisk, and J. D. Thompson, Phys. Rev. Lett. , 84 (2000) 4986

work page 2000
[7]

Muramatsu et al., J

T. Muramatsu et al., J. Phys. Soc. Jpn. , 70 (2001) 3362

work page 2001
[8]

Yashima, S

M. Yashima, S. Kawasaki, H. Mukuda, Y. Kitaoka, H. Shishido, R. Settai, and Y. Onuki , Phys. Rev. , B 76 (2007) 020509R

work page 2007

Show all 25 references

[9]

A. F. Andreev, I. A. Grishchuck, Sov. Phys. JETP , 60 (1984) 267

1984
[10]

T. M. McQueen, A. J.Williams, P.W. Stephens, J. Tao, Y. Zhu, V. Ksenofontov, F. Casper, C. Felser, and R. J. Cava, Phys. Rev. Lett. , 103 (2009) 057002

2009
[11]

F. C. Hsu, J. Y. Luo, K. W. Yeh, T. K. Chen, T. W. Huang, et al., PNAS, 105 (2008) 14262

2008
[12]

Papaniculau, Nucl

N. Papaniculau, Nucl. Phys. , B 240 (1984) 281

1984
[13]

Andreas L¨auchli, Fr´ed´eric Mila, and Karlo Penc, Phys. Rev. Lett , 97 (2006) 087205

2006
[14]

Raghu, S

S. Raghu, S. A. Kivelson, and D. J. Scalapino , Phys. Rev., B 81 (2010) 224505

2010
[15]

Johnston, Adv

D.C. Johnston, Adv. Phys, 59 (2010) 803

2010
[16]

Mazin, Nature, 464 (2010) 183

I.I. Mazin, Nature, 464 (2010) 183

2010
[17]

Graser, T

S. Graser, T. A. Maier, P. J. Hirshfeld, D. J. Scalapino,, New J. Phys. , 11 (2009) 025016

2009
[18]

A.V, Chubukov , Annual Review of Condensed Matter Physics, 3 (2012) 57

2012
[19]

Medvedev, T

S. Medvedev, T. M. McQueen, I. A. Troyan, T. Palasyuk T, M. I. Eremets , et al. , Nat. Mater. , 8 (2009) 630

2009
[20]

Terashima, N

T. Terashima, N. Kikugawa, S. Kasahara, T. Watashige, T. Shibauchi, et al. , J. Phys. Soc. Jpn. , 84 (2015) 063701

2015
[21]

Y. Sun, S. Pyon, T. Tamegai,, Phys. Rev., B 93 (2016) 104502

2016
[22]

Kothapalli, A

K. Kothapalli, A. E. B¨ohmer, W. T. Jayasekara, B. G. Ueland, P. Das, et al., Nature. Commun., 7 (2016) 12728

2016
[23]

Bendele, A

M. Bendele, A. Amato, K. Conder, M. Elender, H. Keller, H.-H. Klauss, H. Luetkens, E. Pom- jakushina,A.Raselli, and R.Khasanov , Phys. Rev. Lett., 104 (2010) 087003

2010
[24]

Bendele,A

M. Bendele,A. Ichsanow, Yu. Pashkevich, L. Keller, Th. Str ¨assle, A. Gusev, E. Pomjakushina, K. Conder, R. Khasanov, and H. Keller , Physical Review , B 85 (2012) 064517

2012
[25]

T. Imai, K. Ahilan, F. L. Ning, T. M. McQueen, and R. J. Cava, Phys. Rev. Lett. , 102 (2009) 177005. p-5

2009

Pith tools

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

[1] [1]

V. L. Ginzburg, Sov.Phys.JETP, 4 (1957) 153

work page 1957

[2] [2]

Saxena, P

S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. Haselwimmer, M. Steiner, E. Pugh, I. Walker, S. Julian, P. Monthoux, G. Lonzarich, A. Huxley, I. Sheikin, D. Braithwaite, and J. Flouquet , Nature (London), 406 (2000) 587

work page 2000

[3] [3]

D. Aoki, A. Huxley, E. Ressouche, D. Braithwaite, J. Flouquet, J-P. Brison, E.Lhotel, and C. Paulsen, Nature (London), 413 (1974) 613

work page 1974

[4] [4]

J. G. Bednorz and K. A. M¨uller, Z. Phys, B 64 (1986) 189

work page 1986

[5] [5]

Aarts Steglich, C

J. Aarts Steglich, C. D. Bredl, W. Lieke, D. Meschede, W. Franz, and H. Schafer, Phys. Rev. Lett, 43 (1979) 1892

work page 1979

[6] [6]

Hegger C

H. Hegger C. Petrovic, E. G. Moshopoulou, M. F. Hundley, J. L. Sarrao, Z. Fisk, and J. D. Thompson, Phys. Rev. Lett. , 84 (2000) 4986

work page 2000

[7] [7]

Muramatsu et al., J

T. Muramatsu et al., J. Phys. Soc. Jpn. , 70 (2001) 3362

work page 2001

[8] [8]

Yashima, S

M. Yashima, S. Kawasaki, H. Mukuda, Y. Kitaoka, H. Shishido, R. Settai, and Y. Onuki , Phys. Rev. , B 76 (2007) 020509R

work page 2007