{"id":"b694359f-5836-4d77-9675-c50bfb962d7c","arxiv_id":"2509.07307","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An ab initio study finds Nb3Sn's superconducting gap is anisotropic yet fully open with 3D orbital pairing, and attributes the Hc2 drop at the martensitic transition to longer coherence lengths from Fermi-velocity redistribution.","lead":"This paper reports a first-principles calculation of the anharmonic vibrations, superconducting gap, and upper critical field of Nb3Sn, the superconductor used in high-field magnets. It argues that the drop in critical field at the structural transition comes from a redistribution of Fermi-surface coherence lengths, and suggests doping rules.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tetragonal phase is an imposed experimental structure, not an ab initio result; Hc2 mechanism depends on it.","rationale":"The reader's weakest assumption exactly captures the most load-bearing risk: the tetragonal phase is not an ab initio output but an experimental input. The paper's own text admits that SSCHA finds no tetragonal minimum, so all tetragonal superconducting properties are computed at a constrained, stressed structure. Since the Hc2 mechanism is based on the difference between cubic and tetragonal Fermi-surface properties, any error in the tetragonal structure propagates directly into the central claim. The PBE0 footnote indicates the energy surface is functional-sensitive, making it plausible that a relaxed tetragonal structure could differ from the experimental reference. The quantitative failure for tetragonal Hc2 (14 T vs 21 T experimental) is consistent with this concern, though it is not itself the root cause. I agree with the reader that the verdict should remain conditional: the work is a substantial advance, but the Hc2 explanation needs either a stabilized tetragonal phase or a clear reframing as a calculation for the experimental geometry rather than a prediction. The proposed test—repeating SSCHA with a functional that stabilizes the tetragonal phase—would settle whether the coherence-length redistribution survives in a self-consistent tetragonal structure.","tokens_in":19270,"tokens_out":5608,"duration_ms":64827,"concrete_test":"Retrain the MLIP on PBE0 (or SCAN) forces and repeat the SSCHA relaxation with the cell and internal coordinates free to distort into a tetragonal stationary point (e.g., enforce P42/mmc symmetry and allow ε and δ to relax under zero external stress). If a tetragonal minimum is found, compare its ε,δ to the experimental values (-0.006,-0.003) and recompute λ, Δk, vF, and ξ0(k) at the relaxed geometry. If the new average ξ0 and Hc2(0) differ by more than 20% from 6.6 nm and 14 T, the original Hc2 mechanism is not robust to the structural assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the Hc2 drop across the martensitic transition is caused by a redistribution of coherence lengths—rests entirely on the superconducting properties of the tetragonal phase. But the tetragonal phase is not a stable solution of the calculation: Sec. II states that 'within SSCHA all tetragonal configurations... relax back to the cubic structure' and that the experimental structure is adopted 'as a reference.' The tetragonal phonons, λ, Δk, vF(k), and ξ0(k) are therefore computed at a non-stationary point of the anharmonic free energy, under a residual internal stress of ~0.2 GPa. The SSCHA phonon spectra for this constrained structure are not guaranteed to describe the physical tetragonal phase; a small change in the lattice parameters or dimerization could substantially alter the Fermi-surface pockets near Γ that the paper identifies as the source of the long-ξ0 regions. The PBE0 test in footnote [43], which finds a ~5 meV/atom tetragonal stabilization, shows that the energy surface is functional-dependent and that the true relaxed tetragonal geometry is not established. If the relaxed tetragonal structure differs from the experimental one, the quantitative Hc2 estimate (14 T vs 21 T experimental) and even the coherence-length mechanism could change. Thus the Hc2 explanation is conditional on an external structural input that the ab initio method does not reproduce.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper combines SSCHA anharmonic lattice dynamics (accelerated by machine-learned potentials) with full-bandwidth anisotropic Migdal-Eliashberg calculations (EPW) to describe cubic and tetragonal Nb3Sn. It claims to resolve three long-standing issues: (i) the martensitic transition is weakly first-order and anharmonicity stabilizes both phases; (ii) the superconducting gap is strongly anisotropic but fully open, with pairing contributions from both longitudinal and transverse Nb d-orbitals; and (iii) the experimentally observed Hc2 drop across the transition is explained by a combination of weaker electron-phonon coupling and a redistribution of Fermi velocities that increases the average Pippard coherence length from 3.9 nm (cubic) to 6.6 nm (tetragonal). The paper additionally proposes that Sn-site doping could enhance transverse-state coupling and raise Tc and Hc2, while Nb-site doping reinforces Hc2 at the cost of Tc.","tokens_in":19597,"tokens_out":5809,"duration_ms":67913,"significance":"If confirmed, this would be the first fully ab initio microscopic account of the superconducting properties of the two phases of Nb3Sn, a technologically central superconductor. The state-of-the-art methodology is a genuine strength: SSCHA-MLIP reproduces the neutron-scattering phonon spectra, including the temperature renormalization of the Gamma_12^+ mode, and the full-bandwidth solution of the Migdal-Eliashberg equations avoids the numerical smearing artifacts that plague earlier DFT studies. The paper also provides a concrete, falsifiable microscopic mechanism for the Hc2 reduction across the martensitic transition, which is a long-standing puzzle in the A15 community. However, the validity of the central Hc2 claim rests on the tetragonal-phase treatment, which the authors themselves report is not a stable minimum of the SSCHA free energy.","major_comments":[{"comment":"The central Hc2 mechanism is computed for a tetragonal structure that the authors' own SSCHA calculation does not stabilize. Sec. II states that 'within SSCHA all tetragonal configurations... relax back to the cubic structure' and that the experimental structure is adopted 'as a reference.' The tetragonal Fermi-surface redistribution, the average coherence length 6.6 nm, and the resulting Hc2 estimate of 14 T in Table II and Eq. (2) are all derived from this non-stationary point. The PBE0 test in footnote [43] indicates that the energy surface is functional-dependent, so the correct relaxed tetragonal geometry is not established. Since a small change in the tetragonal distortion could alter the Gamma-centered FS pockets that are identified as the origin of the long-ξ0 regions, the central claim (iii) is conditional. The authors should either compute the tetragonal superconducting propert","section":"Sec. II and V.B"},{"comment":"The paper claims that the martensitic transition is 'weakly first-order between two nearly-degenerate minima, both stabilized by anharmonic effects' (Sec. IV and Conclusions). This is internally inconsistent with the SSCHA free-energy landscape of Fig. 2(b), which shows a single cubic minimum and no tetragonal stationary point; Sec. II explicitly states that all tetragonal configurations relax back to cubic. A transition between a minimum and a non-stationary point is not a first-order transition. To support conclusion (i), the authors need either a second minimum in the anharmonic free energy (e.g., under finite stress or with a different functional), or they should clearly state that their calculation does not reproduce the martensitic transition and lower the strength of the claim accordingly.","section":"Sec. IV and Conclusions, Fig. 2(b)"},{"comment":"The tetragonal Coulomb pseudopotential μ* is obtained by rescaling the cubic μ* by the DOS ratio (Table S2) rather than from a separate first-principles RPA/KO calculation, unlike the cubic μ*. This choice directly affects the tetragonal gap Δ_k, and through Eq. (1) the coherence length ξ0(k) = ħ v_F / (π Z_k Δ_k), and therefore the predicted Hc2. The qualitative conclusion may survive, but the quantitative estimate (14 T vs 21 T experimental) and even the extent of the coherence-length redistribution depend on this ad hoc rescaling. The authors should report a sensitivity analysis, e.g., using the same μ* for both phases or computing a fully first-principles tetragonal μ*, to demonstrate that the Hc2 mechanism is robust.","section":"Sec. V.A and Table S2"}],"minor_comments":[{"comment":"The main text (Sec. VII) states that PBEsol was used, while the Supplemental Material Section I refers to the 'Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional.' Please clarify which functional was used for the final calculations and, if different, whether results change.","section":"SM Section I vs main text"},{"comment":"The sentence 'The structures were pre-relaxed using and then relaxed within SSCHA' is incomplete; please fix.","section":"Sec. VII"},{"comment":"The notation for the Ginzburg-Landau coherence length is inconsistent: ξ^c_GL(0) in the table header versus ξ_c^GL(0) in the text and the footnote. Please unify. Also define all symbols in the caption consistently.","section":"Table II"},{"comment":"The figure caption describes panels (a-d), (b-e), and (c-f), but the text refers to '(a) and (d)', '(b) and (e)', '(c) and (f)'. The panel labeling appears inconsistent; please check.","section":"Fig. 4 caption"},{"comment":"The sentence 'may justify why the experimentally observed Tc reduction...' is a fragment and reads awkwardly. Integrate it into the main text.","section":"Footnote [61]"},{"comment":"Typo: 'coherence lenght' should be 'coherence length'.","section":"SM Fig. 7 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is transparent about the tetragonal-phase limitation, but that limitation is load-bearing for the paper's most important claim. The SSCHA result that tetragonal Nb3Sn is never a minimum directly contradicts the 'weakly first-order' language, and the ad hoc tetragonal μ* adds further uncertainty to the quantitative Hc2 numbers. If the authors can address the structural issue, the paper would be a strong contribution to the field. I recommend major revision rather than rejection, because the phonon and gap computations are otherwise well executed and the Hc2 mechanism is physically plausible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first SSCHA-MLIP plus full-bandwidth Migdal-Eliashberg study of Nb3Sn, and the vibrational part is very good. The anharmonic phonon dispersions match neutron data convincingly, including the temperature renormalization of the Gamma12 mode, and the comparison with smearing-stabilized DFPT shows why that old trick fails. The gap distribution -- strongly anisotropic but fully open, with both d-parallel and d-perp Nb states pairing -- is a solid, new result that fits recent specific heat and tunneling data. That section alone is worth the paper.\n\nThe soft spot is the one the authors flag themselves but don't fully own. SSCHA relaxes every tetragonal configuration back to cubic; there is no tetragonal minimum at zero pressure. Then the paper adopts the experimental tetragonal cell as a reference and computes all tetragonal phonons, gaps, and coherence lengths on that constrained, stressed structure. That makes the 'weakly first-order transition' claim unsupported -- if the cubic is the only minimum in the calculation, the transition order is not determined. And the Hc2 mechanism, the most advertised result, depends entirely on the tetragonal Fermi surface and coherence-length redistribution from that imposed geometry. The quantitative result also misses: 14 T versus 21 T experimental, a 50% error, while the cubic estimate is within range. The authors note the discrepancy but do not test sensitivity to the assumed lattice parameters or dimerization. Footnote 43, where PBE0 gives roughly 5 meV/atom stabilization for the tetragonal structure, makes clear the energy surface is functional-dependent and the relaxed tetragonal geometry is not established. So the Hc2 story is conditional on an external input, and the 'ab initio answer' framing overstates it.\n\nI also note the tetragonal mu* is rescaled from the cubic value by the DOS ratio (Table S2). That is not a free fit to the target Hc2, but it does inject a chosen input into the tetragonal Tc.\n\nWhere does this leave the paper? It is a state-of-the-art calculation with one genuinely strong section (anharmonic phonons) and one interesting but not yet load-bearing claim (Hc2 mechanism). The authors are transparent about the structural instability, which I credit. A serious referee should engage; the paper deserves peer review, but the transition-order claim needs reframing and the Hc2 mechanism needs either a stabilized tetragonal structure or a sensitivity analysis across plausible structures. I would send it for review and ask for those changes.","headline":"The anharmonic phonons and gap anisotropy are real progress; the Hc2 and first-order claims are built on a tetragonal phase that the calculation itself refuses to stabilize.","tokens_in":20090,"tokens_out":2455,"would_cite":true,"duration_ms":30161,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.25.Dw","74.20.Pq","74.70.Ad"],"model":"deepseek-v4-flash","headline":"This paper presents the first fully ab initio microscopic description of superconducting Nb3Sn and shows that the measured drop of the upper critical field across the martensitic transition is caused by a redistribution of Fermi velocities","keywords":["Nb3Sn","A15 superconductors","upper critical field","anharmonicity","Migdal-Eliashberg","martensitic transition","coherence length","first-principles calculation"],"falsifier":"Compute the tetragonal electronic structure at a genuinely relaxed tetragonal minimum (e.g., using the PBE0 functional that the paper's footnote 43 says gives a larger tetragonal energy gain) and recalculate ξ0 and Hc2; or measure the Fermi-velocity redistribution in tetragonal Nb3Sn directly via quantum oscillations. If the measured clean-limit Hc2 of tetragonal Nb3Sn does not fall near the predicted ~14 T lower bound, the coherence-length mechanism is wrong or the fixed experimental structure is the wrong reference.","tokens_in":19129,"feed_emoji":"🧲","tokens_out":8427,"duration_ms":80737,"temperature":0.7,"pith_summary":"This paper provides the first fully ab initio description of the superconducting state of both the cubic and tetragonal phases of Nb3Sn, the workhorse material for high-field superconducting magnets. It establishes that anharmonic lattice effects, treated with the stochastic self-consistent harmonic approximation, stabilize both structures and reproduce measured phonon spectra. The central discovery is that the experimentally observed drop of the upper critical field Hc2 from 29 T to 21 T across the martensitic transition is not a generic disorder effect; it follows from a calculated redistribution of Fermi velocities that makes the average Pippard coherence length grow from about 3.9 nm to 6.6 nm. The paper also resolves earlier debates on the transition order (weakly first-order) and on the gap structure (strongly anisotropic but fully open, with a three-dimensional pairing mechanism involving both longitudinal and transverse Nb d orbitals). If correct, these results give a material-specific microscopic foundation for optimizing Nb3Sn via targeted doping.","feed_headline":"Nb3Sn's Hc2 drop comes from Fermi-surface reshaping, not disorder","feed_subtitle":"A full ab initio calculation shows the martensitic transition stretches coherence lengths from 3.9 to 6.6 nm, lowering Hc2 from 32 to 14 T.","key_machinery":"The load-bearing machinery is the combination of (i) the SSCHA treatment of anharmonicity, which removes the imaginary phonons that make harmonic theory unstable and reproduces neutron-scattering spectra, and (ii) the momentum-resolved Pippard coherence length ξ0(k) = ħvF(k)/(πZ_kΔ_k), with Z_k = 1 + λ_k, evaluated from the full-bandwidth anisotropic Migdal–Eliashberg solution. The coherence length distribution—controlled mainly by the renormalized Fermi velocity v*_F = v_F/Z rather than by the gap—is what converts the Fermi-surface reshaping into a numerically specific Hc2 value via Hc2(0) = φ0/(2π ξ_GL(0)^2).","core_discovery":"The paper's core claim is that the suppression of Hc2 in tetragonal Nb3Sn is a Fermi-surface effect. Solving the full-bandwidth anisotropic Migdal–Eliashberg equations on anharmonic phonons, the authors compute the momentum-resolved Pippard coherence length ξ0(k) = ħvF(k)/(π Z_k Δ_k) and find that the tetragonal distortion reshapes the Fermi surface so that the Γ-centered electron pockets, short-coherence regions in the cubic phase, become long-coherence regions. The average coherence length increases from 3.9 to 6.6 nm, lowering the clean-limit Hc2(0) from 32 T to 14 T, bracketing the experimental change from 29 T to 21 T. The same calculation yields a strongly anisotropic but nodeless supe","pith_inferences":["The same coherence-length-redistribution mechanism could be tested in other A15 compounds with martensitic transitions (e.g., V3Si), where a similar Hc2 drop is debated; the paper's machinery is transferable.","The claim that the tetragonal phase is stabilized by internal stress suggests that controlled stress engineering (e.g., strain from substrates) could be a cleaner alternative to chemical doping for preserving cubic-phase Hc2.","The paper's clean-limit Hc2 formula neglects Pauli limiting and orbital pair-breaking in the dirty/intermediate regime; a full Hc2(T) curve from the same Fermi-surface data could sharpen the comparison and expose how much of the 14 T underestimate is due to the fixed experimental tetragonal structure.","The predicted doping trade-off between Tc and Hc2 for Nb-site vs Sn-site doping can be tested by measuring the gap anisotropy (e.g., point-contact spectroscopy) in Al-doped vs Ti-doped samples."],"forward_implications":["If correct, the martensitic transition in Nb3Sn is a weakly first-order transition between two anharmonically stabilized structures, not a soft-mode Peierls transition; the tetragonal phase may be stabilized in real samples by internal stress.","The superconducting gap is single-valued and fully open, so two-gap scenarios from early specific-heat and point-contact data are not supported; pairing involves transverse Nb d orbitals and is three-dimensional.","The Hc2 drop across the martensitic transition is intrinsic and Fermi-surface driven: suppressing the transition (e.g., by epitaxial growth or stress-free synthesis) would preserve the higher Hc2.","Sn-site doping (e.g., Al) should raise both Tc and Hc2 by strengthening transverse-state coupling and isotropizing the gap; Nb-site doping (Ti, Ta, Hf) reinforces Hc2 through scattering but lowers Tc by weakening the strongly coupled chain states.","The ab initio clean-limit Hc2(0) of 32 T for cubic Nb3Sn is a lower bound; real samples in the intermediate-to-dirty limit should have slightly higher values, consistent with experiment."],"supporting_citations":[{"why":"Supplies the SSCHA method used to compute anharmonic free energy surfaces and phonons.","marker":"[25]"},{"why":"Supplies the full-bandwidth Migdal–Eliashberg approach used to solve the gap equations without smearing artifacts.","marker":"[26]"},{"why":"Demonstrates the SSCHA+MLIP workflow on NbTi, the methodological template this work transfers to Nb3Sn.","marker":"[15]"},{"why":"The EPW implementation used for the anisotropic Migdal–Eliashberg calculations.","marker":"[55]"},{"why":"Source of experimental Tc, Hc2, and martensitic transition data for Nb3Sn.","marker":"[2]"},{"why":"Documents the Hc2 drop and the elastic-constant anomalies associated with the martensitic transition.","marker":"[9]"},{"why":"Provides experimental λ, mean free paths, and coherence lengths used for comparison.","marker":"[54]"},{"why":"Gives the strong-coupling renormalization of the coherence length and the Hc2 formulas used in the clean limit.","marker":"[59]"},{"why":"Supplies the first-principles RPA and Kukkonen–Overhauser Coulomb pseudopotentials µ* used to set Tc.","marker":"[57]"},{"why":"Documents the harmonic lattice instability that motivates the anharmonic treatment.","marker":"[16]"}],"fun_headline_variants":["First ab initio Nb3Sn model ties Hc2 drop to Fermi-surface shape","Nb3Sn's Hc2 collapse explained by Fermi-surface reshaping","Anharmonic phonons stabilize Nb3Sn, reshape Fermi surface, cut Hc2","Full anisotropic calculation reveals why Nb3Sn's Hc2 plummets","Coherence-length jump from 3.9 to 6.6 nm drives Nb3Sn's Hc2 drop"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"All tetragonal results are computed for the experimental tetragonal geometry even though the calculation's own anharmonic relaxation finds no tetragonal minimum and relaxes every tetragonal starting point back to cubic; if the physically relaxed tetragonal structure differs, the computed Fermi velocities, coherence lengths, and the Hc2 reduction could change substantially.","fun_headline_variants_meta":{"raw":{"variants":["First ab initio Nb3Sn model ties Hc2 drop to Fermi-surface shape","Nb3Sn's Hc2 collapse explained by Fermi-surface reshaping","Anharmonic phonons stabilize Nb3Sn, reshape Fermi surface, cut Hc2","Full anisotropic calculation reveals why Nb3Sn's Hc2 plummets","Coherence-length jump from 3.9 to 6.6 nm drives Nb3Sn's Hc2 drop"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00065,"raw_usage":{"total_tokens":2854,"prompt_tokens":817,"completion_tokens":2037,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":1935}},"tokens_in":561,"tokens_out":2037,"duration_ms":17802,"temperature":1.0,"reasoning_tokens":1935,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T22:27:02.648507+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the tetragonal electronic structure at a genuinely relaxed tetragonal minimum (e.g., using the PBE0 functional that the paper's footnote 43 says gives a larger tetragonal energy gain) and recalculate ξ0 and Hc2; or measure the Fermi-velocity redistribution in tetragonal Nb3Sn directly via quantum oscillations. If the measured clean-limit Hc2 of tetragonal Nb3Sn does not fall near the predicted ~14 T lower bound, the coherence-length mechanism is wrong or the fixed experimental structure is the wrong reference.","supporting_citations":[{"cited_title":"Wu, S.-T","cited_arxiv_id":null,"evidence_quote":"Supplies the SSCHA method used to compute anharmonic free energy surfaces and phonons."},{"cited_title":"Tarantini, F","cited_arxiv_id":null,"evidence_quote":"Demonstrates the SSCHA+MLIP workflow on NbTi, the methodological template this work transfers to Nb3Sn."},{"cited_title":"Mentink, M","cited_arxiv_id":null,"evidence_quote":"Provides experimental λ, mean free paths, and coherence lengths used for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the first-principles RPA and Kukkonen–Overhauser Coulomb pseudopotentials µ* used to set Tc."},{"cited_title":"Sadigh and V","cited_arxiv_id":null,"evidence_quote":"Documents the harmonic lattice instability that motivates the anharmonic treatment."}],"review_version":1}