{"id":"a77fa379-de88-4454-a73a-4ffd2e8d7d7c","arxiv_id":"1909.01313","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Spin-orbit coupling leaves the standard s+- pairing intact for typical iron pnictides, but can promote helical triplet pairing in hole-doped systems and s-wave over d-wave pairing in electron-doped systems.","lead":"This paper uses computer models of iron-based superconductors to test whether spin-orbit coupling changes the pairing symmetry of the superconducting state. It finds the usual pairing pattern is largely unchanged for typical pnictides, but new triplet pairing can appear in heavily hole-doped systems and the preferred symmetry switches in electron-doped systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The doping-regime ordering claims rest on eigenvalues the authors themselves call quantitatively unreliable; a full-LGE check is needed before accepting the triplet window or the d/s swap.","rationale":"The reader's weakest assumption identified the RPA framework generally; my concern is narrower and stronger: the authors themselves state that the LGE eigenvalues are quantitatively unreliable, yet the new doping-regime conclusions are specifically about eigenvalue ordering. The undoped s+− robustness claim is less affected, since it is a large-gap statement supported by kernel and gap-function plots. The hole-doped triplet window and the electron-doped swap are subtle near-degeneracy statements, and no numerical eigenvalues are provided for either. Because the manuscript is transparent about these limitations, I do not recommend rejection; the CONDITIONAL verdict already reflects the need for verification. The proposed full-LGE computation would settle whether the ordering is physical or an artifact of the Fermi-surface-projected static approximation. My agreement with the reader is partial: the concerns share the same methodological family, but the load-bearing point here is the internal admission of eigenvalue unreliability plus the absence of eigenvalue data, rather than the general RPA premise alone.","tokens_in":37011,"tokens_out":8670,"duration_ms":95000,"concrete_test":"Solve the full linearized gap equation of Sec. S2 (Eq. S18) without the Fermi-surface projection and intraband-only replacements of Eqs. (S21)-(S24), for the heavily hole-doped case (µ0 = -150 meV, U = 0.80 eV, λSOC = 75 meV) and the heavily electron-doped cases (µ0 = 20 meV and 25 meV, U = 0.80 eV, λSOC = 0 and 75 meV), including interband pairing, finite Matsubara sums, and all Fermi-surface sheets; tabulate the top five eigenvalues. If the helical triplet is not within about 10% of the leading singlet eigenvalue over a finite U range, or if the d/s swap disappears, the claimed SOC-induced ordering changes are artifacts of the projected static RPA.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. IV A 2 the authors state that \"our analysis is quantitatively not reliable with respect to the LGE eigenvalues λ\", and Sec. S2 admits losing \"quantitative control over the precise criterion that determines the onset of a Cooper instability.\" The central doping-regime claims are, however, eigenvalue-ordering statements: the helical pseudo-spin triplet \"can become the leading instability\" in a narrow window near the SDW instability (Sec. IV A 2), and SOC swaps d-wave and s+− singlet solutions in the electron-doped case (Sec. IV A 3). No eigenvalues are tabulated for either ordering change; a single \"factor of two\" gap is quoted for one hole-doped point. Near-degenerate eigenvalues are precisely where the neglected finite-energy, interband, and omitted corner-pocket contributions (all acknowledged in Sec. IV A 2) can reorder the solutions. The authors even note the triplet state might be stabilized by the neglected terms, so the claimed narrow window is not secured. This is a load-bearing gap between the method and the conclusions drawn from it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies how spin-orbit coupling modifies spin-fluctuation mediated pairing in iron-based superconductors. The authors use a ten-band Hubbard-Hund model with DFT-derived hopping parameters for LaFeAsO and FeSe, include an onsite SOC term, and solve the Fermi-surface-projected linearized gap equation (LGE) with an RPA pairing kernel in which SOC enters both the bandstructure and the two-particle vertex. For the generic LaFeAsO-derived band they find that SOC up to 100 meV leaves the s+− pseudospin-singlet solution as the leading instability and does not generate SOC-specific gap oscillations; in the strongly hole-doped case the leading solution is d-wave pseudospin singlet, but an SOC-promoted helical pseudospin-triplet solution can become leading in a narrow parameter window near the SDW instability; in the strongly electron-doped case SOC swaps the leading d-wave singlet for an s-wave singlet. A FeSe band is analyzed separately, with the acknowledged caveat that the RPA susceptibility peaks at (π,π) rather than at the experimental (π,0)/(0,π) wavevectors. The supplementary material provides a detailed derivation of the BdG setup, the pseudospin gauge, and the RPA pairing vertex.","tokens_in":37197,"tokens_out":5062,"duration_ms":48290,"significance":"If the reported eigenvalue orderings are robust, this is a useful systematic study of a question that has received comparatively little attention in multiband iron-based superconductors: the role of SOC in selecting not only the superconducting symmetry class but also the momentum-space gap structure. The paper's strengths include the realistic ten-band modeling, an explicit and carefully documented pseudospin construction, extensive parameter sweeps over U and λSOC, and an honest discussion of the limitations of the Fermi-surface-projected RPA approach. The predicted narrow helical-triplet window in the hole-doped system and the SOC-driven d/s switching in the electron-doped system are concrete, falsifiable statements that could guide more complete numerical work or experiments. However, the central quantitative claims are weakened by the authors' own admission that the LGE eigenvalues are not quantitatively reliable, which is exactly the observable on which the doping-regime ordering conclusions rest. With additional eigenvalue-level evidence or appropriately weakened conclusions, the paper would be a valuable contribution.","major_comments":[{"comment":"The two flagship doping-regime results — that the helical pseudospin-triplet solution can become the leading instability near the SDW boundary in the strongly hole-doped system and that SOC swaps the leading d-wave and s+− pseudospin-singlet solutions under strong electron doping — are both eigenvalue-ordering statements derived from the Fermi-surface-projected LGE, Eq. (3). In Sec. IV A 2 the authors state that \"our analysis is quantitatively not reliable with respect to the LGE eigenvalues λ,\" and Sec. S2 explains that after the Fermi-surface projection and the redefinition of λ, \"we strictly speaking lose quantitative control over the precise criterion that determines the onset of a Cooper instability.\" No eigenvalues are tabulated for either ranking change; the only quantitative datum given is that the subleading triplet starts out with λ smaller than the leading value by a factor of two. Because near-degenerate eigenvalues are exactly where the neglected finite-energy, interband, and omitted corner-pocket contributions (all acknowledged in Sec. IV A 2) can reorder the hierarchy, the central claims are not yet secured. I request a quantitative presentation of λ for the leading singlet and triplet solutions as a function of λSOC and U in both doping regimes, and preferably a check of the ordering against the non-projected or interband-extended LGE; if this is not feasible, the conclusions should be correspondingly weakened.","section":"Sec. IV A 2 and Sec. S2"},{"comment":"The FeSe analysis is based on a model whose static RPA susceptibility is explicitly acknowledged to peak at (π,π) rather than at the experimental (π,0)/(0,π), so the magnetic-fluctuation input that drives the pairing kernel is not the one realized in FeSe. Since the main conclusion drawn from this calculation — that SOC does not induce qualitatively new gap features despite the tiny pockets — is produced by a kernel built from the wrong magnetic spectrum, it is not empirically supported by the calculation as presented. The authors are transparent about this limitation, but the conclusion still overreaches: a weak-coupling baseline with a quantitatively incorrect susceptibility can at most motivate future work. I suggest either removing the FeSe-specific conclusions from the abstract and summary, or supplementing the calculation with a pairing kernel built from a corrected or experimentally constrained susceptibility (or explicitly labeling the FeSe section as a purely illustrative methodological exercise with no FeSe-specific claims).","section":"Sec. IV B and Sec. S6"},{"comment":"The claim that the helical triplet \"may even become the leading instability\" is further weakened by the authors' own observations that the triplet solution \"becomes leading by accident\" in a narrow parameter window and that it \"is possible, however, that the pseudospin triplet solution can be stabilized by further decrease of temperature or finite-energy and inter-band contributions\" that are neglected in the Fermi-surface-projected LGE. Because this statement appears as a central result in the abstract, it needs either stronger numerical evidence — for example, a demonstration that the accidental window persists under a controlled extension of the approximation — or a clear downgrade to a qualitative tendency. As written, the abstract's wording overstates what the calculation establishes.","section":"Sec. IV A 2"}],"minor_comments":[{"comment":"There is a typo in the abstract: \"thes+−\" should read \"the s+−\".","section":"Abstract"},{"comment":"The captions of Figs. 11 and 16 give chemical potentials as \"µ0 = −150 eV\" and \"µ0 = 20 eV\", respectively, whereas the running text correctly uses meV; please make the units consistent.","section":"Fig. 11 and Fig. 16 captions"},{"comment":"The caption of Fig. 21 refers to panels (a)-(c) and (d)-(f), but the figure and the text use only four panels (a)-(d); update the caption to match the actual panel layout.","section":"Fig. 21 caption"},{"comment":"There are small typos in the text: \"LGE eingevalueλ\" in Sec. IV A 1 should be \"LGE eigenvalue λ\", and in Sec. IV A 3 the sentence \"dominated by intra- dxz and intra dxz components\" should presumably read \"intra-dxz and intra-dyz components\".","section":"Secs. IV A 1 and IV A 3"},{"comment":"Reference 46 lists the journal as \"Phys. Phys. Lett.\"; this should be \"Phys. Rev. Lett.\".","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a competent and useful RPA study, and the main methodological machinery is standard. My principal concern is the mismatch between the authors' explicit caveats about the quantitative reliability of the LGE eigenvalues and the strength of the abstract-level claims about triplet promotion and d/s switching. The FeSe section also draws a conclusion from a kernel built on a susceptibility that is known to disagree with experiment. Both issues are fixable with additional analysis or appropriately softened claims, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the systematic reference calculation for SOC in spin-fluctuation pairing in iron-based superconductors that the field was missing. Its strongest result is the negative one: for generic LaFeAsO-like bands, up to 100 meV SOC leaves the leading s± pseudospin singlet essentially unchanged and produces no SOC-specific gap oscillations. The genuinely new findings are the helical triplet tendency in the heavily hole-doped case and the d-wave-to-s-wave swap in the heavily electron-doped case. The paper uses standard RPA with a Fermi-surface-projected linearized gap equation, but the supplementary material is unusually complete—BdG setup, pseudospin gauge construction, the pairing vertex with SOC, susceptibility plots, and orbital-space form factors. The citation pattern is fair: the authors build on Saito et al. and on Vafek-Chubukov and Eugenio-Vafek, and they cite their own earlier susceptibility papers where those results are actually used. They also earn credit for reporting that the weak-coupling RPA susceptibility for FeSe peaks at (π,π) rather than the experimental (π,0)/(0,π), which properly limits what they claim for that material.\n\nThe soft spots are localized. No code or data are shipped, so the parameter sweep is not independently reproducible; that is a real cost for a purely numerical paper. More substantively, the hole-doped triplet claim rests on eigenvalue ordering in a narrow parameter window near an SDW instability. The authors themselves state that the quantitative eigenvalue criterion is not reliable, and they report only a factor-of-two separation at one point. Their abstract hedges with 'may even become the leading instability,' which is honest, but the claim still inherits more fragility than the rest of the paper. A serious referee should ask for the actual eigenvalues and Stoner factors in that window, and ideally a full, non-projected LGE check for one representative point. The electron-doped d/s swap is visually supported by the gap-function plots, though eigenvalues are not tabulated there either; I see that as a minor omission, not a flaw in the result.\n\nWho is this for? People working on multiband pairing theory and FeSC gap structure. The central robustness result is solid, and the doping-regime tendencies are worth knowing even if the narrow triplet window needs more work. I would accept this for peer review and push for a revision that tabulates eigenvalue data or releases the code.","headline":"A systematic and honest RPA reference calculation; the central s± robustness claim is solid, while the hole-doped triplet window is softer than the abstract suggests and deserves more numerical transparency.","tokens_in":37743,"tokens_out":3856,"would_cite":true,"duration_ms":38380,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.20.-z","71.70.Ej","74.70.Xa","74.25.Ha"],"model":"deepseek-v4-flash","headline":"Spin-orbit coupling does not alter the leading $s_{\\pm}$ state in the generic iron-pnictide band, but in strongly hole-doped systems it promotes a helical pseudospin-triplet solution that can become the leading instability near the…","keywords":["iron-based superconductors","spin-orbit coupling","spin-fluctuation pairing","linearized gap equation","random-phase approximation","pseudospin singlet and triplet pairing","s± pairing","helical triplet pairing"],"falsifier":"Re-solve the linearized gap equation for the same ten-band models without the Fermi-surface projection and with the full Matsubara-frequency-dependent pairing vertex, including interband pairing; if the s-wave and d-wave eigenvalues in the strongly electron-doped case, or the singlet and helical-triplet eigenvalues near the spin-density-wave instability in the strongly hole-doped case, switch order, the paper's SOC-driven switching claims would be overturned. A cleaner experimental check is a phase-sensitive measurement of the gap on heavily electron-doped FeSe films with SOC tuned by strain or gating.","tokens_in":36787,"feed_emoji":"⚛️","tokens_out":11206,"duration_ms":104591,"temperature":0.7,"pith_summary":"The paper asks whether spin-orbit coupling (SOC), which is sizable in iron-based superconductors, can change which superconducting-pairing instability the spin-fluctuation mechanism produces. Using a ten-band Hubbard-Hund model with realistic DFT-derived hoppings and an RPA pairing vertex, the authors solve the Fermi-surface-projected linearized gap equation and compare the leading gap solutions with and without SOC. For the generic LaFeAsO-derived band they find that even 100 meV SOC leaves the leading $s_{\\pm}$ pseudospin-singlet state intact and does not create SOC-specific gap oscillations. The interesting changes appear at extreme doping: strong hole doping gives a leading d-wave singlet with a helical pseudospin-triplet solution promoted by SOC that can win in a narrow window near the spin-density-wave instability, while strong electron doping makes SOC swap the leading solution from d-wave to s-wave. These results matter because gap symmetry is what phase-sensitive experiments can test, and they isolate which SOC-related effects actually shape the pairing state.","feed_headline":"Spin-orbit coupling flips the leading pairing symmetry at high doping","feed_subtitle":"In generic bands the s± state survives; heavy hole doping favors a helical triplet, heavy electron doping favors s-wave.","key_machinery":"The Fermi-surface projected linearized gap equation (LGE) in the intraband approximation, whose eigenvalues $\\lambda$ rank the pairing instabilities, with the static spin-fluctuation pairing vertex obtained in the random-phase approximation from a ten-band Hubbard-Hund model. The load-bearing object is the pseudospin degree of freedom: a Kramers-index basis constructed so that, at vanishing SOC, pseudospin reduces to physical spin, while at finite SOC the problem stays block-diagonal in pseudospin singlet-triplet space and the SOC-induced entanglement is carried by the orbital-to-band matrix elements. The pseudospin triplet-$x$ and triplet-$y$ sectors become coupled, lifting the degeneracy of helical solutions, and the paper shows that this coupling, rather than the SOC-induced anisotropy of the spin susceptibility, is what reorganizes the pairing kernel when SOC is present.","core_discovery":"Working with a ten-orbital model in the 2-Fe Brillouin zone, the authors find that the SOC-induced magnetic anisotropy of the paramagnons is not what controls the pairing kernel; instead the orbital-to-band matrix elements, which entangle spin, orbital, and momentum, determine the structure of the kernel. In this setting the leading instability of the undoped and lightly doped generic band remains the $s_{\\pm}$ (pseudo-)spin singlet for SOC up to $\\lambda_{\\mathrm{SOC}} = 100\\,\\mathrm{meV}$, with no SOC-characteristic oscillations in the gap along the Fermi surfaces. In the strongly hole-doped case, where only three $\\Gamma$-centered hole pockets remain, the leading solution is a d-wave pseudospin singlet, but SOC lifts the degeneracy among triplet solutions and promotes a helical pseudospin-triplet solution that becomes the largest eigenvalue in a narrow parameter window near the spin-density-wave instability. In the strongly electron-doped case with only $M$-centered electron pockets, SOC changes the leading solution from d-wave pseudospin singlet to s-wave pseudospin singlet. For the FeSe band the same weak-coupling framework reproduces the known failure of RPA—the spin susceptibility peaks near $(\\pi,\\pi)$ rather than $(\\pi,0)/(0,\\pi)$—and SOC merely shifts the accidental nodes of the $s_{\\pm}$ solution without generating new gap features.","pith_inferences":["Beyond the paper: because the paper traces the SOC effect to orbital-to-band matrix elements, a model with stronger interorbital hybridization might show SOC-induced gap changes even at modest $\\lambda_{\\mathrm{SOC}}$; this could be tested by varying the $d_{xz}/d_{yz}$--$d_{xy}$ hybridization in the same calculation.","Beyond the paper: the narrow triplet-leading window near the spin-density-wave instability suggests that the competition between singlet d-wave and helical triplet should be strongest where pressure or doping tunes the system closest to magnetism, so tuning experiments on hole-doped pnictides could be used to search for the triplet state.","Beyond the paper: the electron-doped d-to-s swap may partly be a Fermi-surface-topology effect, because SOC opens small hole pockets, rather than a pure spin-fluctuation-polarization effect; a calculation that switches SOC on while keeping the Fermi surface fixed would separate the two contributions.","Beyond the paper: the FeSe caveat implies that conclusions from weak-coupling RPA for FeSe should be viewed as a baseline; a similar analysis in a nematic or strong-coupling setting could decide whether SOC affects the measured gap in the opposite direction."],"forward_implications":["For the generic LaFeAsO-type band, experimental observations of strong gap anisotropy along the Fermi surface should not be attributed to SOC within spin-fluctuation theory; the mechanism must lie elsewhere, such as orbital-selective effects or nematicity.","In strongly hole-doped compounds with only $\\Gamma$-centered hole pockets (KFe$_2$As$_2$-like topology), SOC can push a helical pseudospin-triplet state to the top of the eigenvalue hierarchy near a spin-density-wave instability, making phase-sensitive and spin-resolved experiments the decisive test.","In strongly electron-doped systems (monolayer or intercalated FeSe), including SOC favors $s_{\\pm}$ over d-wave pairing, so the d-wave solution found at vanishing SOC is a consequence of neglecting SOC.","The relative insensitivity of the generic-band gap to SOC means that SOC-induced changes in the pairing kernel are controlled by Fermi-surface and orbital-matrix-element effects rather than by paramagnon spin anisotropy alone.","For FeSe, the weak-coupling RPA susceptibility mismatch with experiment persists after adding SOC, so the SOC-induced node shifts in the FeSe gap are not a quantitative material prediction."],"supporting_citations":[{"why":"Supplies the LaFeAsO DFT-derived hopping parameters that define the generic ten-band bandstructure used for most of the calculation.","marker":"Ref. 60"},{"why":"Supplies the FeSe hopping parameters and the self-consistent renormalization procedure used to build the FeSe Fermi surface.","marker":"Ref. 39"},{"why":"Earlier work on SOC-induced magnetic anisotropy in this model; the current paper builds on it to construct the SOC-aware pairing kernel.","marker":"Ref. 49"},{"why":"Theory of triplet pairing with renormalized interactions in the heavily hole-doped case; the paper reproduces and compares that scenario.","marker":"Ref. 46"},{"why":"Symmetry classification of allowed pairing states with SOC for two electron pockets; defines the s-wave, d-wave, and helical p-wave options tested in the electron-doped section.","marker":"Ref. 47"},{"why":"Provides the procedure for reconstructing full Brillouin-zone gap functions from the Fermi-surface projected linearized gap equation.","marker":"Ref. 61"},{"why":"Earlier inclusion of SOC in the superconducting-state spin response; motivates the treatment of SOC within the same framework.","marker":"Ref. 54"},{"why":"Identified that weak-coupling RPA for FeSe puts the spin susceptibility peak at $\\pi,\\pi$ rather than $\\pi,0$, which the FeSe section confirms and uses as a caveat.","marker":"Ref. 65"},{"why":"Example of SOC changing the leading instability from singlet to helical triplet in a related multiband superconductor; motivates the possibility explored for hole-doped FeSCs.","marker":"Ref. 55"}],"fun_headline_variants":["SOC flips pairing symmetry in doped iron superconductors","Heavy doping lets spin-orbit coupling pick the gap symmetry","Spin-orbit coupling redraws pairing symmetry at extreme doping","Spin-orbit coupling reshapes superconducting gaps only at high doping"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a static, Fermi-surface-projected RPA pairing kernel with only intraband pairing correctly orders the superconducting instabilities even in the presence of strong spin-orbit coupling.","fun_headline_variants_meta":{"raw":{"variants":["SOC flips pairing symmetry in doped iron superconductors","Heavy doping lets spin-orbit coupling pick the gap symmetry","Spin-orbit coupling redraws pairing symmetry at extreme doping","Spin-orbit coupling reshapes superconducting gaps only at high doping"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00066,"raw_usage":{"total_tokens":3126,"prompt_tokens":1162,"completion_tokens":1964,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":778,"completion_tokens_details":{"reasoning_tokens":1907}},"tokens_in":778,"tokens_out":1964,"duration_ms":14827,"temperature":1.0,"reasoning_tokens":1907,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:21:32.730060+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-solve the linearized gap equation for the same ten-band models without the Fermi-surface projection and with the full Matsubara-frequency-dependent pairing vertex, including interband pairing; if the s-wave and d-wave eigenvalues in the strongly electron-doped case, or the singlet and helical-triplet eigenvalues near the spin-density-wave instability in the strongly hole-doped case, switch order, the paper's SOC-driven switching claims would be overturned. A cleaner experimental check is a phase-sensitive measurement of the gap on heavily electron-doped FeSe films with SOC tuned by strain or gating.","supporting_citations":[],"review_version":1}