{"id":"a6aa816e-289b-4e88-9e0a-c6e06cef3e9e","arxiv_id":"2601.10903","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Calculations show Slater-mode phonon pairing explains the orientation dependence of superconductivity at KTaO3 interfaces, but the coupling is too weak to explain the observed Tc by itself.","lead":"Superconductivity in KTaO3 interfaces is studied by computing pairing from a soft Slater phonon mode. The theory reproduces the observed orientation dependence of Tc but shows this phonon alone is too weak, so other phonons must contribute.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The pairing kernel's suppression of back-scattering relies on evaluating the dynamic Rashba vertex at k0=(k+k')/2; if the physical vertex instead depends on the incoming momentum k, the λ values and 111/001 hierarchy may change.","rationale":"The most sensitive step in the argument is the momentum argument of VR in Eq. 3. The paper states that the dynamic Rashba terms are evaluated at k0=(k+k')/2, which produces the forward-scattering peak and complete back-scattering suppression that shape the gap function and λ. But this is an assertion, not a derivation; the microscopic electron-phonon coupling for a polar mode that modulates a Rashba term has a vertex that depends on the electron momentum k (or possibly on k and q separately), not on the pair center-of-mass. If the true vertex is g(k) (times the phonon displacement), then the product g(k,q)g(-k',-q) does not vanish at back-scattering; calculations show it can be comparable to or even larger than the forward-scattering value. This would alter the magnitude of λ, potentially making the Slater mode alone sufficient to explain Tc, and it could also change the relative strength of 111 vs 001 because the Fermi surface topology and the degree of back-scattering differ. The reader's concern about interface transferability of bulk coupling constants is real but only affects the overall scale; the k0 assumption affects the momentum structure and thus the central conclusion about the mechanism. We therefore propose a concrete numerical/analytic test: recompute λ with the vertex evaluated at k (or derive the correct vertex from the tight-binding model of Ref. [12]) and compare. If the results are unchanged, the paper's conclusion is robust; if not, the claim as stated needs revision. Because the current acceptance is based on an unverified load-bearing assumption, we recommend CONDITIONAL acceptance pending this check.","tokens_in":13658,"tokens_out":10482,"duration_ms":110953,"concrete_test":"Recompute the pairing kernel in Eq. 3 with the dynamic Rashba vertex taken at the incoming momentum k (or, better, with the full orbital matrix elements from the tight-binding Hamiltonian of Ref. [12] derived without the k0 symmetrization), keeping everything else fixed. If the resulting λ for 111 and 001 changes substantially, or if back-scattering contributions become significant, the forward-scattering interpretation and the orientation hierarchy are not robust. An even more direct check: re-derive the electron-phonon vertex g_{k,q} for the polar TO1 mode from the model Hamiltonian in Ref. [11], and determine whether the matrix element is of the form f(k) or f(k+q/2).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative results (λ111≈0.054, λ001≈0.003, and the conclusion that the Slater mode alone gives too-small λ) are controlled by the momentum dependence of the dynamic Rashba matrix elements in Eq. 3. The paper states without derivation that these elements are evaluated at k0=(k+k')/2, which makes the interaction vanish for back-scattering (k'=-k) and peak for forward scattering. This is asserted to follow from Ref. [11], but the microscopic electron-phonon vertex for a polar mode coupled to a Rashba term is generally proportional to the electron momentum at one vertex (g_{k,q} ∝ k), not to the pair center-of-mass k0. Under the standard vertex g_{k,q} ∝ k, the product g_{k,q} g_{-k',-q} is nonzero at back-scattering, so the qualitative picture of a forward-scattering-dominated kernel is not guaranteed. Since λ is an integral over the Fermi surface, a significant back-scattering contribution could raise λ (possibly above the ~0.26 needed for Tc≈2 K) and alter the 111/001 ratio, undermining both the absolute-magnitude claim and the orientation-dependence argument. The paper's 'evaluated at k0' step in Section II is thus a load-bearing assumption that is neither derived nor checked.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies superconductivity in the 2DEG at KTaO3 interfaces driven by the soft transverse optic (TO1 or Slater) phonon mode. Starting from a microscopic theory of the dynamic Rashba electron-phonon coupling (Gastiasoro et al., Phys. Rev. B 105, 224503) and ab-initio couplings for bulk KTO (Venditti et al., arXiv:2510.25655), the author constructs a linearized gap equation for a bilayer tight-binding model of the (111) and (001) interfaces. The secular matrix includes the four dynamic Rashba terms and the q-dependent TO1 phonon dispersion. The numerical solutions yield the BCS coupling constant λ and the angular- and band-dependent gap function. The main results are: (1) λ for (111) exceeds that for (001) by one to two orders of magnitude for all parameter choices, consistent with the observed orientation dependence of Tc; (2) the gap function is strongly anisotropic and band-dependent; and (3) λ for (111) is at most ~0.13, well below the ~0.26 needed for Tc≈2 K, indicating that the Slater mode alone cannot account for the absolute Tc and other phonons must contribute. The paper is transparent about the approximations and lists possible improvements.","tokens_in":14010,"tokens_out":15733,"duration_ms":159377,"significance":"This is a valuable contribution that puts the phenomenological Slater-mode pairing scenario on a concrete microscopic footing for KTO interfaces. The calculation uses published ab-initio couplings and two tight-binding models, and no parameter is fitted to Tc; the orientation hierarchy is robust across these choices. The explicit gap-function anisotropy is a new prediction that could be tested by tunneling. The honest reporting that the Slater mode alone gives λ too small is an important negative result, narrowing the search for the pairing glue. The main limitations—interface transferability of the bulk couplings, the static-phonon approximation, and the neglect of other phonon modes—are clearly acknowledged. If the results hold, they will likely steer future work toward including LO phonons and full strong-coupling calculations.","major_comments":[],"minor_comments":[{"comment":"The statement that the dynamic Rashba matrix elements are evaluated at k0=(k+k')/2 is important because it suppresses back-scattering, but it is not derived. The stress-test concern that the vertex might depend on the incoming momentum is addressed by the standard symmetrization of a local coupling (the coupling α(x) p·σ yields a vertex proportional to (k+k')/2), but this should be stated explicitly, with a reference to the specific equation in Ref. [11]. Please add a short derivation or an explicit citation.","section":"Section II, Eq. (3)"},{"comment":"The 1/√3 scaling of t_i for the (111) polarization is stated without derivation. Since this factor directly affects λ_111 and hence the central quantitative claim, please show the projection of the phonon eigenvector onto the coupling tensor. Also, consider a brief comment on whether the same projection applies to all four t_i terms.","section":"Section III (111 case)"},{"comment":"The use of the static phonon propagator with a q-dependent energy is a hybrid approximation; a one-sentence justification in terms of the adiabatic parameter ω(2k_F)/E_F would be helpful.","section":"Eq. (3)"},{"comment":"The labels 'ω(q)' and 'ω(0)' in the table and figure headers are clear in context, but a fuller caption or a note in the text defining the abbreviation (e.g., 'ω(q) means the full Vaks dispersion is used') would improve readability.","section":"Table III and Fig. 5"},{"comment":"Minor typos: 'can by thought of' should be 'can be thought of' in Section II; 'Rasbha' should be 'Rashba' in the text before Table II; spacing in Eq. (3) is irregular.","section":"Throughout"}],"recommendation":"minor_revision","confidential_remarks":"To the editor: The paper is sound and within the scope of the journal. The stress-test concern about k0 is, in my view, resolved by the standard symmetrized vertex; the authors should make this explicit. I recommend minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a real but incremental paper. It takes a phenomenological Slater-mode pairing theory for KTO interfaces, plugs in ab-initio dynamic Rashba couplings from another group, solves the linearized gap equation, and finds λ_111 >> λ_001, consistent with the observed orientation hierarchy. It also honestly reports that λ_111 alone is too small (≈0.05 vs 0.26) to explain Tc ≈ 2 K. That negative result is itself useful.\n\nThe genuinely new content is the explicit microscopic calculation: λ for both interfaces, the band- and angle-resolved gap function, and the demonstration that the orientation ordering survives with ab-initio couplings. The paper is transparent about its approximations and lists the obvious improvements, which is good practice.\n\nWhere the paper is softest: the pairing kernel's strong forward-scattering character, which drives everything, comes from evaluating the Rashba vertex at k0=(k+k')/2. That step is asserted and attributed to Ref. [11], not derived here. If the physical vertex instead depends on the incoming momentum k, back-scattering need not vanish and both the absolute λ and the 111/001 ratio could shift. I think the orientation hierarchy is probably safe because it follows from orbital degeneracy and the quadratic scaling of λ with t0, but this is the place I'd ask a referee to check carefully.\n\nOther caveats are minor and acknowledged: bulk couplings transferred to the interface with only a geometric factor; static phonon approximation; bilayer model. These do not undermine the main conclusion.\n\nOverall, a solid, modest advance. It deserves serious refereeing and would be a good reading-group paper for anyone working on oxide interfaces. I'd cite it in any paper touching KTO superconductivity.\n\nBottom line: send it to peer review, and expect the referee to scrutinize the k0 vertex, but the paper is transparent about its inputs and outputs.","headline":"A credible microscopic application of Slater-mode pairing to KTO interfaces, with an honest negative result on λ; the orientation hierarchy is likely robust, but the forward-scattering kernel rests on a vertex assumption worth checking.","tokens_in":104,"tokens_out":5422,"would_cite":true,"duration_ms":188147,"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":"Pairing through the soft Slater phonon reproduces the orientation dependence of superconductivity in KTaO3 interfaces, but cannot alone produce the observed Tc.","keywords":["KTaO3","superconductivity","2D electron gas","Slater mode","electron-phonon coupling","Rashba coupling","quantum paraelectric","oxide heterostructure"],"falsifier":"Measure the superconducting gap anisotropy on a (111) KTaO3 interface by tunneling: the calculation predicts a gap Δnk that varies strongly with band index and Fermi-surface angle, peaking in specific directions tied to the star-shaped Fermi surface. A gap that is essentially isotropic, or a Tc(111)/Tc(001) ratio that disagrees with the computed λ ratio while other phonon contributions are held fixed, would rule out the Slater mode as the dominant pairing agent.","tokens_in":13567,"feed_emoji":"🔬","tokens_out":5285,"duration_ms":50737,"temperature":0.7,"pith_summary":"This paper tests whether the soft transverse-optic (Slater) phonon—the mode that drives ferroelectricity in quantum paraelectrics—can explain the superconductivity observed in the two-dimensional electron gas at KTaO3 heterostructures. Pairing is mediated by phonon-induced Rashba couplings computed from first principles, and the linearized gap equation is solved for (111) and (001) interfaces using two tight-binding models. The computed BCS coupling constant λ is much larger for (111) than for (001), matching the measured ordering of transition temperatures, and the gap function is strongly anisotropic in band index and Fermi-surface angle. However, for realistic parameters λ is about 0.05, far below the roughly 0.26 needed for a Tc near 2 K, so the Slater mode alone cannot account for the absolute magnitude of Tc. The paper concludes that other phonons must add to the pairing.","feed_headline":"Slater phonon explains KTaO3 interface pairing, falls short on Tc","feed_subtitle":"The same microscopic calculation gives a pairing strength too small for the absolute 2 K transition, pointing to extra phonons.","key_machinery":"The load-bearing object is the dynamic Rashba electron-phonon coupling generated by the TO1 (Slater) mode: phonon-induced hoppings among t2g orbitals via oxygen 2p states, with one spin-independent term (t0) and three spin-dependent terms (tA, tB, tC) taken from first-principles calculations. In the Cooper channel these couplings act at the average momentum k0 = (k+k')/2, which makes the pairing interaction peak at forward scattering and vanish at back scattering. The gap equation is built as a secular matrix whose kernel combines these matrix elements with the Vaks-parameterized phonon dispersion ω(q); solving the linearized gap equation yields the BCS constant λ and the order parameter Δnk","core_discovery":"The central claim is that the soft TO1 (Slater) phonon, acting through dynamic Rashba couplings, is a major contributor to the superconductivity of the KTaO3 two-dimensional electron gas and explains why the (111) interface has a much higher Tc than the (001) interface. Solving the linearized gap equation in a bilayer tight-binding model, the paper finds λ111 ≈ 0.04–0.13 and λ001 ≈ 0.001–0.01 (depending on the tight-binding model and phonon dispersion), reproducing the orientation hierarchy observed in experiment. The gap function Δnk varies significantly with band index and Fermi-surface angle. The same calculation gives λ values well below the ~0.26 needed for the observed Tc ≈ 2 K at n2D","pith_inferences":["If the orientation hierarchy is as robust as argued, engineering the confinement potential—via gate voltage or capping layers—to change orbital occupancy could provide a tunable knob for Tc across interfaces, a testable extension beyond the paper's numbers.","The reliance on bulk-derived dynamic Rashba couplings at the interface is the natural place for a more complete heterostructure calculation (coupled Schrödinger–Poisson plus interface screening) to improve the quantitative estimate of λ.","The forward-scattering nature of the interaction raises the possibility that in higher-order channels or under stronger coupling the gap could develop nodes or even pair-density-wave character; tunneling and specific-heat measurements could look for this.","If the Slater-mode contribution is supplemented by other phonons, combined isotope-effect experiments or terahertz/inelastic-neutron studies of the interfacial phonon spectrum would help decide which modes matter."],"forward_implications":["The orientation hierarchy T111 > T110 > T001 in KTaO3 interfaces follows from the orbital degeneracy of the t2g manifold and the projection of the Slater-mode polarization, giving a concrete orbital-based design rule for enhancing Tc.","The predicted strong band-index and in-plane angular dependence of the superconducting gap is observable by tunneling spectroscopy, providing a direct fingerprint of Slater-mode pairing.","Because λ scales roughly as 1/ω(0) when the phonon dispersion is included, experiments that soften the TO1 mode—such as doping toward the ferroelectric quantum critical point—should raise Tc substantially.","The finding that the Slater mode alone gives λ ≈ 0.05 implies that other phonons, especially orientation-dependent high-energy modes, must contribute to reach the observed Tc values; theories focusing only on the Slater mode will underestimate Tc.","The strong forward-scattering peak suggests that pairing at these interfaces is dominated by small-momentum-transfer processes, which may favor unconventional gap structures."],"fun_headline_variants":["Slater phonon sets KTaO3 gap orientation, but Tc needs more","KTaO3 pairing from Slater mode, too weak for 2 K alone","Orientation-dependent KTaO3 superconductivity, coupling shortfall","Slater mode shapes KTaO3 gap, but other phonons required","KTaO3 interface Tc: Slater phonon explains anisotropy, not magnitude"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The computation assumes the dynamic Rashba couplings derived for bulk KTaO3 apply unchanged at the interface, apart from a geometric 1/√3 factor for the (111) polarization; altered screening, confinement, or structural relaxation at the interface would change the computed λ and the absolute conclusion about Tc.","fun_headline_variants_meta":{"raw":{"variants":["Slater phonon sets KTaO3 gap orientation, but Tc needs more","KTaO3 pairing from Slater mode, too weak for 2 K alone","Orientation-dependent KTaO3 superconductivity, coupling shortfall","Slater mode shapes KTaO3 gap, but other phonons required","KTaO3 interface Tc: Slater phonon explains anisotropy, not magnitude"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00014,"raw_usage":{"total_tokens":971,"prompt_tokens":693,"completion_tokens":278,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":437,"completion_tokens_details":{"reasoning_tokens":175}},"tokens_in":437,"tokens_out":278,"duration_ms":3571,"temperature":1.0,"reasoning_tokens":175,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T10:10:03.649730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the superconducting gap anisotropy on a (111) KTaO3 interface by tunneling: the calculation predicts a gap Δnk that varies strongly with band index and Fermi-surface angle, peaking in specific directions tied to the star-shaped Fermi surface. A gap that is essentially isotropic, or a Tc(111)/Tc(001) ratio that disagrees with the computed λ ratio while other phonon contributions are held fixed, would rule out the Slater mode as the dominant pairing agent.","supporting_citations":[],"review_version":1}