{"id":"ebe06e76-126b-4583-9fe2-bfe2ea4a5a36","arxiv_id":"1909.00609","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors derive a dielectric-function-method gap equation for a Dirac cone and compute plasmon-mediated superconducting critical temperatures in graphene that rise to the millikelvin range with density.","lead":"This paper applies the dielectric function method, a non-adiabatic theory of superconductivity, to single-layer graphene with a linear Dirac dispersion. It derives the gap equation and predicts critical temperatures in the microkelvin-to-millikelvin range, which differ from the BCS result at low doping.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (11) contains an internally inconsistent cutoff: Λ is called a wavevector but set to 8 eV, so the high-momentum part of the RPA polarization—and therefore K(0,0) and every Tc in Fig. 2—rests on an unexplained unit conversion.","rationale":"The reader identifies the same weakest assumption: the RPA density response in Eq. (11) and the ambiguous cutoff. I agree that this is the single most load-bearing input. Every numerical result in Fig. 2 and the abstract's quantitative claim depend on the screening function, and the paper neither derives Eq. (11) nor resolves the wavevector-versus-energy meaning of Λ. The BCS comparison is also uncontrolled, but it affects only the qualitative 'better suited to low densities' conclusion; the cutoff affects the DFM curve itself. Because the DFM framework is established and the defect is a fixable clarification rather than a demonstrated conceptual error, the conditional verdict remains appropriate: the numerical claims should not be taken at face value until Eq. (11) and the cutoff convention are verified with an explicit unit-consistent derivation and recomputed Tc values.","tokens_in":7682,"tokens_out":12060,"duration_ms":172813,"concrete_test":"Independent check: re-derive the zero-temperature RPA response of graphene from Eq. (10) (or directly from Wunsch et al.) keeping units explicit, and compute K(0,0) and Tc for n=10^10 cm^-2, κ=1 using (i) Λ as an energy cutoff 8 eV converted by k_c=8 eV/(ħv_F) and (ii) a lattice momentum cutoff k_c=π/a_CC≈2.2 Å^-1. If the two Tc values agree within a factor 2, the concern is minor; if they differ by orders of magnitude, Eq. (11)'s cutoff ambiguity undermines the reported Fig. 2 and the comparison with BCS.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All computed critical temperatures flow from the kernel K(ω,ω′) in Eq. (4), whose only material input is the RPA dielectric function (13), built from the density response χ(q,iΩ) in Eq. (11). The paper does not show the reduction from Eq. (10) to Eq. (11), and the printed expression is dimensionally inconsistent as it stands: q, Ω, k and y are said to be dimensionless in units of k_F and ε_F, while the upper limit of the interband k-integral is written Λ/k_F with Λ described as a cutoff wavevector and then set to Λ≈8 eV [25]. An energy of 8 eV is not a wavevector unless a factor ħv_F is supplied; in graphene ħv_F≈6.6 eV·Å, so the implied momentum cutoff is about 1.2 Å^-1, whereas a lattice cutoff such as π/a_CC≈2.2 Å^-1 is a different number. At the lowest densities in Fig. 2 (k_F≈1.8×10^7 m^-1), either assignment puts the upper limit near 10^3 k_F, so the second integral in Eq. (11) samples a very large momentum region. The choice changes the high-momentum part of χ(q,iΩ), hence K(0,0), and because Tc is exponential in -1/λ through Eq. (6), even a factor-of-two change in the cutoff-dependent screening can shift Tc by orders of magnitude. This is the most load-bearing step for the numerical claim in Fig. 2 and for the asserted contrast with BCS. The DFM formalism itself is established, and Wunsch et al. is a sound source, so the issue is a missing, checkable derivation and an inconsistent cutoff convention, not a conceptual failure of the method.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript applies the dielectric function method (DFM) to calculate the critical temperature for plasmon-mediated Cooper pairing in single-layer graphene with a linear Dirac dispersion. The authors write down a DFM gap equation for the Dirac cone (Eqs. (3)-(4)), combine it with an RPA dielectric function that includes a lattice cutoff and a dielectric constant kappa (Eqs. (8)-(13)), and extract Tc versus carrier density from Eq. (6). Their central result is that Tc lies in the microkelvin-to-millikelvin range for n between 10^10 and 10^12 cm^-2 and has a weaker density dependence than the BCS result of Kopnin and Sonin, which they take as evidence that DFM is better suited than BCS at low doping. The paper also argues that the graphene kernel resembles that of a 3D parabolic band rather than a 2D parabolic band.","tokens_in":8076,"tokens_out":13967,"duration_ms":123081,"significance":"The question addressed is relevant: low-carrier-density superconductors are precisely the regime where DFM is designed to improve on BCS, and graphene provides a clean linear-dispersion model. The DFM/RPA formalism is established, the manuscript is largely self-contained, and I found no circularity in the argument: the input polarization function comes from the independent RPA literature. However, the numerical claim is not yet verifiable because the reduction from Eq. (10) to Eq. (11) is not shown, the cutoff in Eq. (11) is dimensionally inconsistent, and the BCS comparison in Fig. 2 is not controlled. If these points are fixed, the paper can be a useful contribution.","major_comments":[{"comment":"The definition and use of the cutoff Lambda are inconsistent. In Eq. (11) the variables q, Omega, k, and y are dimensionless, so the upper limit of the second integral over k must be a pure number. The text calls Lambda 'a cutoff wavevector' and then sets Lambda about 8 eV [25]. If Lambda is an energy, the dimensionless upper limit is Lambda/(hbar v_F k_F), which equals about 680 at n=10^10 cm^-2 (k_F about 1.8 x 10^7 m^-1, epsilon_F about 0.012 eV) and about 68 at n=10^12 cm^-2. If Lambda is intended as a wavevector, the value 8 eV is meaningless without a conversion via hbar v_F. This order-of-magnitude change in the upper limit across Fig. 2 affects the interband part of chi(q,iOmega), hence epsilon(q,iOmega), the kernel K(0,0), and every Tc through Eq. (6). The authors should state the units of Lambda, show the conversion used in the numerics, and check the sensitivity of Fig. 2 to the cutoff model; they should also justify evaluating the Dirac-cone and RPA expressions up to momenta of several hundred k_F at the lowest densities.","section":"II.B, Eq. (11)"},{"comment":"The passage from the general DFM gap equation (2) to the Dirac-cone equations (3)-(4), and then to the normalized-gap equation (5) and the coupling parameter lambda in Eq. (7), is not shown. The text states that after using the linear dispersion and converting the summation to integration the gap equation 'becomes' Eq. (3), and the Fredholm equation is simply written down. Because the prefactor (omega' + epsilon_F) in Eq. (4) and the Heaviside term in Eq. (7) are not self-evident, and because K and lambda directly determine Tc, a reader cannot reproduce the main numerical results without these derivations or explicit references that contain them in the same notation. I ask the authors to provide the missing steps or exact citations.","section":"II.A, Eqs. (3)-(7)"},{"comment":"The comparison with BCS is not controlled. The caption of Fig. 2 states that the BCS curve is computed for a coupling constant of 0.01 and a Debye window of 5 meV, while the DFM curves are computed for kappa = 1.00-1.08 with no equivalent adjustable parameters. The abstract's claim of a 'significantly different behaviour' of Tc as a function of carrier density is therefore based on a comparison between two models that are not matched in band structure, screening, or interaction parameters. Please either recompute the BCS curve with the same Dirac dispersion and RPA screening as the DFM calculation, or explicitly present the BCS curve as an illustrative curve with its parameters clearly stated.","section":"III, Fig. 2"}],"minor_comments":[{"comment":"The notation 'T <~ Tc' is unclear; the manuscript should use a precise inequality such as T much less than Tc or T <= Tc.","section":"II.A, Eq. (2)"},{"comment":"The symbol lambda is used both for the band index in Eq. (1) and for the coupling parameter in Eq. (6); please use different symbols to avoid ambiguity.","section":"II.A, Eqs. (1) and (6)"},{"comment":"The text says the density range is 10^10-10^12 cm^-2, but the x-axis in Fig. 2 appears to start at 10^11 cm^-2; please make the text and figure consistent.","section":"III and Fig. 2"},{"comment":"Reference [3] contains a typo: 'quantim Hall effect' should be 'quantum Hall effect'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first DFM gap equation for a Dirac cone, and it gives a genuine analytic observation—the kernel for graphene resembles the 3D parabolic case rather than the 2D parabolic one. That's worth remembering.\n\nThe paper is a straightforward application of a known formalism. The authors carry over the DFM machinery from Kirzhnits/Takada, replace the parabolic dispersion with the linear Dirac cone, and plug in the RPA polarization from Wunsch et al. The resulting Tc curves are microkelvin to millikelvin, rising with density and falling with dielectric constant. The physics is plausible, and the authors are careful to mention the BKT caveat and the neglect of flexural phonons. There is no sign of curve-fitting or self-citation abuse; the method itself is established.\n\nThe soft spots are real but not fatal. The biggest is the cutoff Λ in Eq. (11). The text calls it a cutoff wavevector and then sets it to 8 eV. An energy is not a wavevector without a conversion factor. If the intent is Λ = 8 eV / ħ v_F ≈ 1.2 Å^-1, the paper should say so. The second integral in Eq. (11) runs to Λ/k_F, which at the lowest densities is hundreds of k_F. The magnitude of the polarization at those momenta will change with the cutoff, and since Tc depends exponentially on 1/λ, a factor-of-two uncertainty in the cutoff can change Tc by orders of magnitude. This is not a minor typo; it is a missing piece of the calculation. The stress-test note is right on this.\n\nThe second soft spot is the BCS comparison. They plot the Kopnin-Sonin result with fixed coupling (0.01) and Debye window (5 meV), which are not derived from the same electron-plasmon interaction. So the 'significantly different behaviour' is not a controlled comparison. It might well be true that DFM is better suited at low densities, but this figure doesn't demonstrate it. A BCS calculation using the same plasmon kernel with a small Debye window would make the point.\n\nEqs. (5) and (7) are asserted without derivation. They follow the Zubarev method, but the reader has to take them on faith. A short appendix or a few intermediate steps would help.\n\nAll in all, this is a solid piece of work with one load-bearing ambiguity. It deserves peer review, but the cutoff must be clarified and the comparison redone before the central claim is convincing.","headline":"First DFM gap equation for a Dirac cone, with a real analytic observation, but a dimensionally sloppy cutoff and an uncontrolled BCS comparison weaken the numerical claim.","tokens_in":8583,"tokens_out":6421,"would_cite":false,"duration_ms":319200,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper derives a DFM gap equation for graphene's Dirac cone and predicts plasmon-mediated Cooper pairing with critical temperatures in the microkelvin-to-millikelvin range, with density dependence different from BCS.","keywords":["graphene superconductivity","plasmon-mediated pairing","dielectric function method","Dirac cone gap equation","random phase approximation","critical temperature","low carrier density","two-dimensional superconductivity"],"falsifier":"Recompute the kernel (4) using an independent finite-temperature random-phase-approximation polarization function and a cutoff converted consistently as a wavevector, then compare the resulting $T_c$ versus density curve to Fig. 2; a material change in magnitude or density dependence would show that the predicted temperatures are artifacts of the screening model or the cutoff.","tokens_in":7487,"feed_emoji":"❄️","tokens_out":10103,"duration_ms":86035,"temperature":0.7,"pith_summary":"The paper applies the dielectric function method, a non-adiabatic weak-coupling formalism that does not assume a small Debye window around the Fermi level, to single-layer graphene. It derives the gap equation for graphene's linear Dirac dispersion and computes the critical temperature of plasmon-mediated Cooper pairing as a function of carrier density, using an RPA-screened Coulomb interaction that includes plasmons. For densities between $10^{10}$ and $10^{12}\\,\\mathrm{cm}^{-2}$, the predicted critical temperatures lie in the microkelvin-to-millikelvin range, increase as the carrier density rises, and decrease when the environmental dielectric constant grows. The predicted density dependence is significantly gentler than the BCS result, which the paper takes as evidence that the dielectric function method is better suited than BCS or standard electron-phonon theory at low carrier densities.","feed_headline":"Plasmon pairing can make graphene superconducting at microkelvins","feed_subtitle":"A non-BCS dielectric-function calculation puts the transition in the microkelvin-to-millikelvin range and ties it to gate-tunable density.","key_machinery":"The carrying object is the DFM kernel $K(\\omega,\\omega')$ in Eq. (4), which enters the gap equation (3) through the frequency-dependent interaction $V(q,i\\Omega)=2\\pi e^2/(|q|\\epsilon(q,i\\Omega))$. The RPA dielectric function (13) uses the zero-temperature graphene density-density response function (11) to supply screening and plasmon effects, while the linear Dirac dispersion manifests as a prefactor that makes the kernel vanish at the band edge, giving it a 3D-parabolic-like shape rather than the usual 2D shape. The normalized gap function follows from the integral equation (5), and $T_c$ is obtained from Eq. (6) with the coupling parameter $\\lambda$ of Eq. (7).","core_discovery":"The central claim is that the dielectric function method, applied to a Dirac cone and combined with the random-phase-approximation dielectric function, yields the correct weak-coupling description of superconductivity in low-doped single-layer graphene. The derived gap equation (3) with kernel (4) replaces the BCS Debye-window interaction by a frequency-dependent screened interaction, and solving the normalized gap equation gives critical temperatures in the millikelvin-to-microkelvin range. Unlike the BCS curve for Dirac electrons, the DFM critical temperature rises smoothly with carrier density and is suppressed by an increased dielectric constant, so the paper concludes that DFM is more appropriate than standard BCS and Migdal-Eliashberg approaches in this regime.","pith_inferences":["An experimental distinction could come from gate-tuning: measuring the slope of $T_c$ versus carrier density in a clean suspended sample would separate DFM from BCS without needing to resolve the absolute temperature scale.","The same Dirac-cone kernel construction should transfer to other two-dimensional Dirac and Weyl materials and to gapped Dirac systems; the vanishing of the kernel at the band edge is a structural prediction that any alternative theory should reproduce.","Because the RPA dielectric function is taken at zero temperature while the gap equation is used at small but finite $T_c$, a finite-temperature polarization function could shift the low-density results; this is a natural next calculation.","The BCS comparison uses a fixed Debye window and coupling constant; matching those parameters to graphene's actual phonon spectrum would test whether the reported discrepancy is structural or partly parametric."],"forward_implications":["If the central claim is right, clean single-layer graphene at densities near $10^{12}\\,\\mathrm{cm}^{-2}$ should show a superconducting transition around the millikelvin scale, with the transition temperature dropping as the sheet is placed on higher-dielectric-constant substrates.","The density dependence of $T_c$ is a fingerprint: DFM predicts a gentler, monotonic rise with carrier density, whereas the BCS curve for Dirac electrons rises more steeply for the same parameters.","The formalism extends BCS to interaction regions comparable to the Fermi energy, so low-density graphene becomes a test bed for non-adiabatic pairing rather than a violation of the small-window assumption.","Because the calculation addresses pair breaking, the actual 2D phase transition is expected to be BKT-limited; at weak coupling the BKT temperature should lie close to the computed $T_c$."],"supporting_citations":[{"why":"supplies the zero-temperature RPA density-density response function for graphene used to build the dielectric function in Eq. (13)","marker":"[24]"},{"why":"provides the weak-coupling DFM gap equation that the Dirac-cone derivation extends","marker":"[15]"},{"why":"introduces the dielectric function method as a way to describe superconductivity beyond a small Debye window","marker":"[14]"},{"why":"gives the BCS critical temperature for Dirac electrons used as the comparison curve in Fig. 2","marker":"[11]"},{"why":"supports the DFM's applicability to two-dimensional electron systems and the numerical procedure used here","marker":"[16]"},{"why":"supplies the lattice cutoff value that enters the interband polarization integral","marker":"[25]"},{"why":"argues the electron-plasmon mechanism is favorable at low doping, motivating the interaction potential","marker":"[12]"},{"why":"sets the lower bound of the carrier-density range through electron-hole puddle disorder","marker":"[30]"}],"fun_headline_variants":["Graphene superconductivity from plasmonic pairing, not BCS","Dielectric function method puts graphene's Tc in microkelvins","Plasmonic pairing drives graphene's superconductivity at low density","Graphene's Tc rises smoothly with density via plasmonic pairing","Non-BCS plasmonic approach predicts graphene superconductivity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation rests on the zero-temperature random-phase-approximation screening model for graphene and on the conversion of the 8 eV lattice cutoff into a wavevector for an integral upper limit; if either input is wrong, the kernel and every computed critical temperature change.","fun_headline_variants_meta":{"raw":{"variants":["Graphene superconductivity from plasmonic pairing, not BCS","Dielectric function method puts graphene's Tc in microkelvins","Plasmonic pairing drives graphene's superconductivity at low density","Graphene's Tc rises smoothly with density via plasmonic pairing","Non-BCS plasmonic approach predicts graphene superconductivity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000763,"raw_usage":{"total_tokens":3335,"prompt_tokens":845,"completion_tokens":2490,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":461,"completion_tokens_details":{"reasoning_tokens":2402}},"tokens_in":461,"tokens_out":2490,"duration_ms":15775,"temperature":1.0,"reasoning_tokens":2402,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:42:35.431198+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the kernel (4) using an independent finite-temperature random-phase-approximation polarization function and a cutoff converted consistently as a wavevector, then compare the resulting $T_c$ versus density curve to Fig. 2; a material change in magnitude or density dependence would show that the predicted temperatures are artifacts of the screening model or the cutoff.","supporting_citations":[{"cited_title":"Dynamical polarization of graphene at finite doping,","cited_arxiv_id":null,"evidence_quote":"supplies the zero-temperature RPA density-density response function for graphene used to build the dielectric function in Eq. (13)"},{"cited_title":"Plasmon mechanism of superconductivity in two- and three-dimensional electron systems,","cited_arxiv_id":null,"evidence_quote":"provides the weak-coupling DFM gap equation that the Dirac-cone derivation extends"},{"cited_title":"The description of superconductivity in terms of dielectric response function,","cited_arxiv_id":null,"evidence_quote":"introduces the dielectric function method as a way to describe superconductivity beyond a small Debye window"},{"cited_title":"BCS superconductivity of Dirac electrons in graphene layers,","cited_arxiv_id":null,"evidence_quote":"gives the BCS critical temperature for Dirac electrons used as the comparison curve in Fig. 2"},{"cited_title":"Superconductivity in the two- dimensional electron gas induced by high-energy optical phonon mode and large polarization of the SrTiO 3 substrate,","cited_arxiv_id":null,"evidence_quote":"supports the DFM's applicability to two-dimensional electron systems and the numerical procedure used here"},{"cited_title":"Electronic properties of disordered two-dimensional carbon,","cited_arxiv_id":null,"evidence_quote":"supplies the lattice cutoff value that enters the interband polarization integral"},{"cited_title":"Superconducting state of pure and doped graphene,","cited_arxiv_id":null,"evidence_quote":"argues the electron-plasmon mechanism is favorable at low doping, motivating the interaction potential"},{"cited_title":"Observa- tion of electron-hole puddles in graphene using a scanning single-electron transistor,","cited_arxiv_id":null,"evidence_quote":"sets the lower bound of the carrier-density range through electron-hole puddle disorder"}],"review_version":1}