{"id":"9dfe68ec-9031-4b4e-8fe8-b55de5a83e10","arxiv_id":"2411.18315","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Measurements and parameter-free calculations show that Pb0.64Bi0.36 is a strong-coupling superconductor with an anisotropic, overlapped three-gap state, and identify the hexagonal structure as the main driver of its unusually large electron-phonon coupling.","lead":"This paper combines new measurements and first-principles calculations to explain why the Pb0.64Bi0.36 alloy superconducts so strongly, with an electron-phonon coupling near 2 and a critical temperature of about 8.6 K. It shows that the hexagonal crystal structure, electron doping, and spin-orbit effects cooperate to boost the coupling, and it predicts an anisotropic, overlapped three-gap superconducting state.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"All theoretical superconductivity quantities rest on ordered-VCA phonons and the pairing kernel; an explicit-disorder supercell test is needed to confirm λ≈2.05 and the overlapped three-gap-like picture.","rationale":"The reader's weakest-assumption analysis identified the same load-bearing point: all phonon, EPC, and SCDFT superconductivity results use an ordered VCA medium, with disorder validated only indirectly through electronic structure and residual resistivity. I found no additional concern more central than this. The SCDFT Tc overestimate and the fitted μ* in the Eliashberg analysis are real but secondary: they are explicitly discussed by the authors, and they do not directly threaten the structural-transition-driven λ enhancement or the gap-structure claim as much as an unverified disorder effect on phonons/EPC would. The paper deserves credit for the KKR-CPA validation, the residual-resistivity match, and the tunneling α2F comparison, all of which make the VCA concern less severe than in a typical VCA study. However, the central claim in the abstract—that the alloy has an overlapped three-gap-like structure with strong anisotropy—is derived from SCDFT on sharp VCA Fermi-surface sheets in a regime where the paper itself argues disorder scattering is strong (ℏ/τ ≫ ℏωD). A direct supercell test would settle whether the ordered-medium assumption is quantitatively safe. Since the reader already marked the verdict CONDITIONAL and the concern is shared rather than newly decisive, the appropriate recommendation is UNCHANGED.","tokens_in":34213,"tokens_out":3285,"duration_ms":36375,"concrete_test":"Construct ordered hexagonal supercells of 32–64 atoms at x=0.36 with explicit Pb/Bi site occupations (several random or special-quasirandom realizations), relax them, and compute DFPT phonons and EPC with the same ultrasoft pseudopotentials, SOC, and convergence parameters as Sec. II. Average λ and α2F over realizations and compare with the VCA values λ=2.05 and α2F in Fig. 13; if feasible, also run a supercell SCDFT calculation to see whether the three-maximum gap histogram survives. A practical threshold: if the ensemble-averaged λ deviates by more than 10% or the gap distribution loses its three-peak structure, the central claim needs revision; if λ and α2F agree within that margin, the VCA assumption is exonerated for the phonon-driven part of the claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that λ≈2.05 is driven primarily by the fcc-to-hexagonal transition, with the SCDFT gap forming an overlapped three-gap-like anisotropic structure—is computed entirely in the ordered virtual-crystal medium of Sec. II. VCA replaces Pb and Bi by an average potential and weighted atomic mass, so phonon frequencies, linewidths, α2F, and the SCDFT pairing kernel/gap distribution contain no chemical disorder. The paper validates VCA only at the electronic-structure level (Sec. IV C): DOS near EF and band positions match KKR-CPA, and disorder-induced τ≈3 fs reproduces the residual resistivity. It does not validate disorder effects on lattice dynamics or on electron-phonon matrix elements. Because the Pb and Bi pseudopotentials differ appreciably (6p2 vs 6p3 occupations), local potential fluctuations can renormalize mode-specific EPC matrix elements and phonon lifetimes beyond the averaged medium. The paper's own disorder analysis reports ℏ/τ = 240 meV ≫ ℏωD = 9 meV, placing the alloy in the strong-scattering regime; this is precisely where Anderson's theorem need not protect the sharp VCA Fermi-surface sheets used for the anisotropic SCDFT gap calculation, so the three-gap-like distribution could be washed out or reshaped. The agreement of the calculated α2F with the tunneling α2F of Pb0.65Bi0.35 in Fig. 13(h) is genuine independent support, as is the reproduction of elemental Pb, but it does not settle whether explicit disorder shifts λ=2.05 or modifies the gap structure. Thus the load-bearing assumption is that the ordered VCA medium captures the phonons and pairing kernel; this is plausible but unchecked.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper combines experimental measurements (magnetic susceptibility, resistivity, specific heat) with first-principles calculations (DFT, DFPT, isotropic Eliashberg, anisotropic SCDFT) to explain the strong-coupling superconductivity of Pb0.64Bi0.36. The authors report Tc = 8.6 K, λ ≈ 2.05, and propose that the fcc-to-hexagonal structural transition is the primary driver of the enhanced electron-phonon coupling, with additional contributions from Bi doping and spin-orbit coupling. They further argue that, unlike Pb which exhibits two well-separated superconducting gaps, the alloy shows an overlapped three-gap-like anisotropic structure, and that disorder-induced electron scattering (τ ≈ 3 fs) reduces Tc, partially explaining the difference between the SCDFT value (9.6 K) and experiment.","tokens_in":34519,"tokens_out":8321,"duration_ms":74782,"significance":"If the conclusions hold, this work provides a compelling microscopic explanation for the record-strong electron-phonon coupling in a bulk ambient-pressure superconductor and a rare theoretical account of multigap anisotropic superconductivity in a disordered alloy. The paper has clear strengths: parameter-free SCDFT calculations for both Pb and the alloy, careful validation of the VCA electronic structure against all-electron KKR-CPA, reproduction of the residual resistivity, and comparison of the computed α2F with tunneling data for a closely related alloy. The clean decomposition of the λ enhancement into structural, doping, and SOC contributions is a valuable conceptual result. The main weakness is that the central gap-structure claim rests on calculations performed in the ordered VCA medium, while the paper's own disorder analysis places the alloy in a strong-scattering regime where the robustness of the anisotropic gap distribution is not self-evident.","major_comments":[{"comment":"","section":"Sec. IV.E.3, Figs. 19-20, Table V"},{"comment":"","section":"Sec. IV.E.5, Fig. 22"}],"minor_comments":[{"comment":"Equation numbering is inconsistent: the Eliashberg equations are numbered (22) in the main text but the text refers to \"Equations (38),\" and the Supplemental Material also numbers them (38). Please unify the numbering and cross-references.","section":"Sec. IV.E.1 and SM IX"},{"comment":"The table header says \"I (THz) defined with Eq. 32,\" but Eq. (32) appears only in the Supplemental Material; in the main text I is defined in Eq. (21). Please correct the cross-reference.","section":"Table IV and Eq. (21)"},{"comment":"The axis labels in Fig. 15 contain garbled text (e.g., \"BP/BF/BA/BC/BJ...\"), apparently a font or encoding artifact. The horizontal axis should be T/Tc and the vertical axis should be ΔCe/(γTc).","section":"Fig. 15"},{"comment":"Reference [20] is an empty placeholder. The pseudopotential details and the PSlibrary reference should be properly cited.","section":"Ref. [20]"},{"comment":"The phrase \"quite clerly overestimated\" contains a typo; it should read \"quite clearly overestimated.\"","section":"Sec. IV.E.2"},{"comment":"The term \"Heviside function\" should be \"Heaviside function.\"","section":"Eq. (22) and SM Eq. (38)"},{"comment":"The statement in Sec. IV.C that \"electron scattering is not a dominating factor for the electronic structure\" appears to be in tension with the later finding ℏ/τ = 240 meV ≫ ℏωD. Please clarify the distinction between the validity of the VCA band centers and the strong scattering regime relevant for superconductivity.","section":"Sec. IV.C and Sec. IV.E.5"},{"comment":"The agreement between the computed α2F and the tunneling-derived α2F of Pb0.65Bi0.35 is stated qualitatively. A quantitative measure (e.g., a logarithmic-average frequency comparison or a weighted residual) would strengthen the validation claim.","section":"Fig. 13(h)"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely of interest to the readership of Physical Review B, and the combination of experiment with parameter-free SCDFT is valuable. The main concern is the VCA-based gap structure in a strongly disordered alloy; an explicit disorder test or a quantitative broadening analysis is needed before the three-gap-like claim can be considered established. The circularity concern about µ* in the isotropic Eliashberg section is real but not fatal, because the central structural conclusions come from SCDFT and DFPT. I recommend major revision rather than rejection, as the issues are identifiable and potentially addressable within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Rough take: this is the first ab initio study of the hexagonal epsilon-phase of Pb-Bi, and the central argument—that the fcc-to-hexagonal transition, not just Bi doping, pushes lambda up to ~2—looks right. The SCDFT gap structure (three overlapping gap maxima on three Fermi-surface sheets) is plausible but less certain because it comes from an ordered virtual-crystal medium.\n\nWhat the paper does well: it combines new experiments (Tc, specific heat, critical fields) with a careful decomposition of lambda into structural, doping, and SOC contributions. The VCA electronic structure is validated against KKR-CPA, and the residual resistivity is reproduced. The calculated alpha2F matches the old tunneling data for Pb0.65Bi0.35, which is independent evidence that the phonon and EPC approximations are on track. The authors are also upfront about the mu* fitting in the isotropic Eliashberg section and about the SCDFT Tc overestimate.\n\nThe main soft spot is exactly the one the stress-test note flags. The phonons, EPC, and the SCDFT pairing kernel are computed in an ordered VCA medium. The paper argues mass disorder is negligible, but chemical disorder in this alloy is strong by their own numbers: hbar/tau = 240 meV, an order of magnitude above the Debye energy. That puts you in a regime where the sharp Fermi-surface sheets used for the anisotropic gap calculation could be washed out, and the three-gap picture might be smeared into a single, broadened gap. The agreement with tunneling alpha2F helps, but it doesn't test the gap distribution. I'd want an explicit-disorder calculation (supercell with real Pb/Bi, or CPA phonons) before betting on the three-gap claim. That is a substantial revision request, not a rejection.\n\nThe isotropic Eliashberg results with fitted mu* are less concerning than the reader's report suggests, because the authors clearly label that section as constrained, and the parameter-free SCDFT is the primary tool. The overestimated SCDFT Tc is a known trend, and their disorder correction points in the right direction even if it doesn't fully close the gap.\n\nBottom line: this paper deserves peer review. It is a serious, honest, well-documented piece of work on a classic system. I'd suggest the referee ask for a disorder test for the phonons and pairing, or at least a frank paragraph on why the ordered medium is justified in that strong-scattering regime. I'd bring it to a reading group focused on phonon engineering or alloy superconductors.","headline":"Solid first-principles account of Pb-Bi's strong coupling; the structural-transition story is convincing, but the three-gap picture needs a disorder test.","tokens_in":35168,"tokens_out":2839,"would_cite":true,"duration_ms":28043,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Pb$_{0.64}$Bi$_{0.36}$ reaches $\\lambda\\approx 2$ because a cubic-to-hexagonal transition softens its phonons, and, unlike lead, it forms an overlapped three-gap superconducting state.","keywords":["strong-coupling superconductor","electron-phonon coupling","multiband superconductivity","Pb-Bi alloy","virtual crystal approximation","spin-orbit coupling","anisotropic superconducting gap","SCDFT"],"falsifier":"Measure the phonon dispersion of hexagonal Pb$_{0.64}$Bi$_{0.36}$ by inelastic neutron or x-ray scattering and compare with the virtual-crystal calculation: if the low-frequency optical modes near the $\\Gamma$ point are significantly harder than predicted (so that $\\hat{\\omega}^2$ does not drop below the fcc Pb value), the claimed structural-transition-driven enhancement of $\\lambda$ would be refuted. Alternatively, high-resolution tunneling or specific-heat measurements that resolve only one gap, or two well-separated gaps, rather than the predicted overlapped three-gap structure, would also disprove the central claim.","tokens_in":33978,"feed_emoji":"⚛️","tokens_out":10626,"duration_ms":85965,"temperature":0.7,"pith_summary":"This paper claims that the hexagonal alloy Pb$_{0.64}$Bi$_{0.36}$ owes its exceptionally strong superconductivity (electron-phonon coupling $\\lambda \\approx 2.05$ and $T_c \\approx 8.6$ K, among the strongest under ambient pressure) primarily to the cubic-to-hexagonal structural transition, which softens the lattice vibrations, with additional contributions from Bi electron doping and spin-orbit coupling. It further claims that lead is a two-gap superconductor with well-separated gaps, whereas the alloy forms an overlapped, strongly anisotropic three-gap-like structure. The evidence combines heat-capacity, magnetization, and resistivity measurements on a polycrystalline sample with first-principles electronic, phonon, and electron-phonon calculations that treat the alloy as an ordered virtual crystal and are validated by KKR-CPA and by comparison with tunneling-derived Eliashberg functions. The work also shows that the chemical disorder inherent to the alloy reduces $T_c$ from the clean limit, a result supported by both the measured residual resistivity and the calculated scattering time.","feed_headline":"Hexagonal switch gives Pb-Bi its lambda approx 2 and a three-gap state","feed_subtitle":"First-principles and heat-capacity data trace the strong coupling to softened phonons, with disorder pulling Tc down.","key_machinery":"The machinery at work is the decomposition of the electron-phonon coupling constant as $\\lambda = 2I/\\hat{\\omega}^2$, where $I = \\int \\omega \\alpha^2F(\\omega)\\,d\\omega$ measures the electronic contribution through phonon linewidths and $\\hat{\\omega}^2$ is the average squared phonon frequency. Applied to a sequence of systems (fcc Pb, hexagonal Pb, and Pb$_{0.64}$Bi$_{0.36}$, each with and without spin-orbit coupling), this ratio isolates how the structural transition, electron doping, and relativistic effects each contribute to the enhancement of $\\lambda$. The phonons and electron-phonon matrix elements are computed with density functional perturbation theory using a virtual crystal approximation for the alloy, validated against KKR-CPA electronic structure and the measured residual resistivity; the anisotropic superconducting gap is solved with density functional theory for superconductors (SCDFT). The same machinery also yields the key qualitative contrast: two well-separated Fermi-surface sheets in Pb give two distinct gap maxima, while three sheets with mixed sp character in the alloy give overlapping gap distributions.","core_discovery":"The central discovery is that the extraordinary electron-phonon coupling of the Pb-Bi alloy comes mainly from the structural transition from fcc to hexagonal (nearly hcp) coordination, and that the superconducting state is a strongly anisotropic, overlapped three-gap-like structure. By writing $\\lambda = 2I/\\hat{\\omega}^2$, the authors separate the electronic factor $I$ (proportional to phonon linewidths) from the phononic factor $\\hat{\\omega}^2$ (average squared phonon frequency). Going from fcc Pb to a hypothetical hexagonal Pb, $\\hat{\\omega}^2$ drops from 1.79 to 1.35 THz$^2$ while $I$ rises slightly, producing a jump in $\\lambda$ from 1.47 to 2.03; adding Bi electrons and spin-orbit coupling changes these factors further ($I$ to 1.45 THz$^2$, $\\hat{\\omega}^2$ to 1.43 THz$^2$), bringing the final $\\lambda$ to 2.05 in the isotropic picture and 2.08 as the Fermi-surface average in SCDFT. On the gap structure, SCDFT solutions give bimodal and clearly separated gap values on the two Fermi surface sheets of Pb (averages 1.30 and 1.44 meV), while the three sheets of the alloy carry overlapping gap distributions (averages 1.84, 1.88, and 1.93 meV) whose spread narrows with temperature, consistent with the deviations from single-gap s-wave behavior seen in specific heat and critical fields. Finally, the short electron scattering time $\\tau \\approx 2.75$ fs extracted from the calculation is shown to lower $T_c$, moving the SCDFT value of 9.6 K toward the experimental 8.6 K.","pith_inferences":["An editor's inference: the structural-transition mechanism identified here could serve as a search strategy for new strong-coupling superconductors, namely alloys that undergo symmetry-lowering transitions with an accompanying softening of the phonon spectrum.","If disorder truly suppresses $T_c$, then reducing chemical disorder, through ordering, thin-film growth, or eliminating defects, might push Pb$_{0.64}$Bi$_{0.36}$ closer to its clean-limit $T_c$, a testable material-engineering route.","The overlapped three-gap structure implies a strongly momentum-dependent order parameter; this could manifest in field-angle-dependent thermal conductivity or in the temperature dependence of the magnetic penetration depth, providing independent experimental checks.","The two-step validation strategy (KKR-CPA for the electronic structure plus measured resistivity for the scattering time) may be transferable to predicting superconductivity in other disordered alloys where virtual-crystal phonons are used."],"forward_implications":["The fcc-to-hexagonal transition is the dominant factor in raising $\\lambda$ from 1.47 in fcc Pb to 2.03 in the hexagonal structure, implying that phonon softening from a structural transition can outweigh electron-count effects in strong-coupling superconductors.","The overlapped three-gap structure explains the experimentally observed deviations from single-gap BCS behavior in the temperature dependence of specific heat and of the lower and upper critical fields.","Chemical disorder suppresses $T_c$: the clean-limit SCDFT value of 9.6 K is reduced toward the measured 8.6 K when scattering at $\\tau \\approx 2.75$ fs is included, suggesting that more ordered samples of the same composition could superconduct at higher temperatures.","Spin-orbit coupling contributes significantly to $\\lambda$ by lowering phonon frequencies (a 31% increase in the alloy relative to the scalar-relativistic case), so relativistic effects are essential for quantitative predictions in heavy-element superconductors.","The computed Eliashberg function for Pb$_{0.64}$Bi$_{0.36}$ agrees with tunneling-derived data for Pb$_{0.65}$Bi$_{0.35}$, validating the virtual-crystal approach for electron-phonon properties of this alloy class."],"supporting_citations":[{"why":"supplies the tunneling-measured Eliashberg function and $\\lambda=2.13$ for Pb$_{0.65}$Bi$_{0.35}$ against which the computed $\\alpha^2F$ and $\\lambda$ are validated","marker":"[5]"},{"why":"first-principles study of cubic Pb-Bi/Tl-Bi alloys that this work extends to the hexagonal phase","marker":"[17]"},{"why":"provides the KKR-CPA method used to validate the virtual-crystal electronic structure and to compute disorder scattering and residual resistivity","marker":"[22]"},{"why":"provides the SCDFT implementation used to compute $T_c$ and the anisotropic $\\mathbf{k}$-resolved gap distribution","marker":"[35]"},{"why":"supplies the Allen-Dynes strong-coupling formula used to extract $\\lambda$ from experimental $T_c$ and to compare with the computed $T_c$","marker":"[49]"},{"why":"provides experimental phonon dispersion for Pb used to validate the computed phonons and Eliashberg functions","marker":"[68]"},{"why":"supplies the ab initio two-band superconducting gap calculations for Pb that this work contrasts with the three-gap-like structure of the alloy","marker":"[78]"}],"fun_headline_variants":["Hexagonal switch gives Pb-Bi strong coupling and three-gap state","Pb-Bi alloy's lambda~2 traced to hexagonal coordination, three-gap superconductivity","Anisotropic three-gap superconductivity in Pb-Bi alloy from hexagonal structure","Disorder lowers Tc in strong-coupling Pb-Bi alloy with three overlapping gaps","Why Pb-Bi has lambda~2: hexagonal structure and anisotropic gaps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the alloy can be represented as an ordered 'virtual crystal': atoms are replaced by an average pseudopotential and the average atomic mass, and the lattice vibrations and pairing are computed in that averaged medium; if real chemical disorder substantially changes the phonons or the pairing beyond what the averaged medium captures, the computed $\\lambda \\approx 2.05$ and the three-gap distribution would shift.","fun_headline_variants_meta":{"raw":{"variants":["Hexagonal switch gives Pb-Bi strong coupling and three-gap state","Pb-Bi alloy's lambda~2 traced to hexagonal coordination, three-gap superconductivity","Anisotropic three-gap superconductivity in Pb-Bi alloy from hexagonal structure","Disorder lowers Tc in strong-coupling Pb-Bi alloy with three overlapping gaps","Why Pb-Bi has lambda~2: hexagonal structure and anisotropic gaps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000952,"raw_usage":{"total_tokens":4165,"prompt_tokens":1151,"completion_tokens":3014,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":767,"completion_tokens_details":{"reasoning_tokens":2910}},"tokens_in":767,"tokens_out":3014,"duration_ms":19950,"temperature":1.0,"reasoning_tokens":2910,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:18:48.239050+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the phonon dispersion of hexagonal Pb$_{0.64}$Bi$_{0.36}$ by inelastic neutron or x-ray scattering and compare with the virtual-crystal calculation: if the low-frequency optical modes near the $\\Gamma$ point are significantly harder than predicted (so that $\\hat{\\omega}^2$ does not drop below the fcc Pb value), the claimed structural-transition-driven enhancement of $\\lambda$ would be refuted. Alternatively, high-resolution tunneling or specific-heat measurements that resolve only one gap, or two well-separated gaps, rather than the predicted overlapped three-gap structure, would also disprove the central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the SCDFT implementation used to compute $T_c$ and the anisotropic $\\mathbf{k}$-resolved gap distribution"},{"cited_title":"Helfand and N","cited_arxiv_id":null,"evidence_quote":"supplies the Allen-Dynes strong-coupling formula used to extract $\\lambda$ from experimental $T_c$ and to compare with the computed $T_c$"},{"cited_title":"Gutowska, A","cited_arxiv_id":null,"evidence_quote":"provides experimental phonon dispersion for Pb used to validate the computed phonons and Eliashberg functions"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the ab initio two-band superconducting gap calculations for Pb that this work contrasts with the three-gap-like structure of the alloy"}],"review_version":1}