{"id":"a58c85bf-f65f-4a2e-a152-29e4ebd83d49","arxiv_id":"1908.02433","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Electron doping can reduce the lattice thermal conductivity of suspended graphene by about 50% at 200 K through phonon-electron scattering, according to first-principles calculations.","lead":"This paper uses computer simulations to show that adding electrons to a sheet of graphene can cut its heat conduction by about half at a temperature of 200 K. The extra electrons make atomic vibrations scatter more, which could help turn heat into electricity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted 50% reduction relies on a rigid-band Fermi shift that places E_F at the LDA Van Hove singularity at 1.8 eV; self-consistent doping could weaken this peak and with it the central claim.","rationale":"The paper is a competent application of established methods (ShengBTE + EPW) to a timely question, and the low-doping validation against Ref. [6] (Fig. 3) gives some confidence in the e-ph matrix elements. The central claim, however, is not the low-doping behavior; it is that at the high carrier densities accessible by gating (n about 4e14 cm^-2; EF about 1.8 eV), phonon-electron scattering cuts the lattice thermal conductivity by roughly half at 200 K. The mechanism is a peak in the electron-phonon coupling at EF=1.8 eV, coinciding with the DOS maximum (Van Hove singularity) in the rigidly shifted LDA band structure (Fig. 7). The entire prediction therefore depends on the rigid-band approximation being quantitatively accurate on the high-energy side of the Dirac cone, well above the region validated in Fig. 3. A self-consistent doped calculation would test this directly. Our concern does not contradict the reader's verdict; it sharpens it: the conditional acceptance should require this test and also an explicit pristine thermal conductivity baseline, which is currently missing. We therefore keep the verdict unchanged.","tokens_in":5113,"tokens_out":13340,"duration_ms":131421,"concrete_test":"Perform self-consistent DFT calculations for electron doping at n=2.86e13 cm^-2 and n=4e14 cm^-2 by changing the electron number in the unit cell with a compensating uniform positive background (or using a suitable supercell), and recompute the band structure, DOS at the Fermi level, and e-ph matrix elements. Then recompute the phonon-electron scattering rates and thermal conductivity at 200 K for the same Fermi energies as in the paper. If the self-consistent DOS at high doping is noticeably lower than the rigid-shift DOS, or if the resulting kappa reduction falls below 50%, the central claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central 50% reduction at 200 K is obtained by comparing EF=0.2 eV with EF=1.8 eV (Sec. III). At EF=1.8 eV, the LDA density of states reaches its maximum (Fig. 7b), i.e., the Fermi level sits at the M-point Van Hove singularity, and the phonon-electron scattering rates become comparable to the anharmonic p-p rates. All rates are computed with Allen's model using a rigid shift of the undoped LDA band structure (Secs. II and III); no self-consistent treatment of the doped electron system is performed. The rigid-band approximation is plausible at low carrier density (validated against Ref. [6] in Fig. 3 at n=2.86e13 cm^-2), but at concentrations corresponding to EF=1.8 eV (n of a few times 10^14 cm^-2) the added carriers strongly modify screening and can renormalize the band structure; the position and strength of the Van Hove singularity may shift or broaden. The calculation's Gaussian smearing (0.025 Ry) also smooths the singularity. If the true doped DOS at E_F is lower than the rigid-shift value, the p-e scattering rates fall and the 50% reduction is not realized. Since the paper provides no convergence test or experimental validation at high doping, this assumption is the least secure pillar of the claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles calculations of the lattice thermal conductivity of suspended graphene including phonon-electron (p-e) scattering, with the Fermi level rigidly shifted to simulate electron doping. Phonon-phonon (p-p) scattering rates are computed with a modified ShengBTE package, and p-e scattering rates with Allen's model using EPW electron-phonon matrix elements. The main claim is that at large electron doping the lattice thermal conductivity can be reduced by about 50% at 200 K, specifically when comparing EF = 0.2 eV with EF = 1.8 eV, where the electron density of states and electron-phonon coupling are maximal. The authors also compare the calculated phonon dispersion with inelastic x-ray data and the electrical resistivity with a prior calculation at n = 2.86e13 cm^-2.","tokens_in":5379,"tokens_out":2959,"duration_ms":33357,"significance":"If the central claim holds, the paper identifies a practical, gate-controlled route to reducing graphene's lattice thermal conductivity, which would be relevant for thermoelectric applications. The work has notable strengths: the calculations involve no post-hoc fitting; the resistivity benchmark in Fig. 3 uses an independent earlier calculation; the phonon dispersion in Fig. 1 is validated against experiment; and the p-p scattering rates in Fig. 2 reproduce the expected asymptotic behavior for acoustic modes. However, the quantitative 50% reduction rests on a rigid-band Fermi shift into the LDA Van Hove singularity, and the manuscript does not report the undoped baseline or any convergence tests. These gaps make the central quantitative claim less secure than the presentation suggests.","major_comments":[{"comment":"The 50% reduction is defined as the difference between EF = 0.2 eV and EF = 1.8 eV, not as a reduction relative to undoped graphene. The abstract states a '~50% reduction of the lattice thermal conductivity at 200 K' without specifying this baseline, which is misleading because the practically relevant quantity for thermoelectric figures of merit is the reduction relative to the pristine, undoped value. The authors should report the undoped lattice thermal conductivity at 200 K and state the baseline explicitly in the abstract and conclusions.","section":"Abstract and Sec. III (Fig. 6)"},{"comment":"The central result relies on the rigid-band approximation at EF = 1.8 eV, where the LDA density of states reaches its maximum at the M-point Van Hove singularity. The manuscript provides no self-consistent treatment of the doped electron system, no assessment of screening-induced band-structure renormalization, and no sensitivity test for the Gaussian smearing of 0.025 Ry, which smooths the singularity. If the true doped DOS at EF is lower than the rigid-shift LDA value, the p-e scattering rates and hence the predicted 50% reduction would be reduced. The authors should test this by performing self-consistent doped calculations (e.g., with a Fermi-level smearing approach) or at least by varying the smearing and showing that the predicted reduction is robust.","section":"Secs. II and III, especially Fig. 7(b)"},{"comment":"No convergence tests are reported for the 140x140x1 q-grid, the 200x200x1 k-grid, the 0.025 Ry Gaussian smearing, or the 10-nearest-neighbor cutoff for third-order force constants. Since the p-e scattering rates depend on electronic states near EF and the DOS peak, numerical broadening may affect the quantitative result. The authors should provide convergence tests with respect to these parameters, particularly for the p-e rates at EF = 1.8 eV.","section":"Sec. II (computational parameters)"},{"comment":"The undoped lattice thermal conductivity is never quoted or compared with the well-known experimental range of 2600-5300 W/mK (Refs. [1-3]). The phonon dispersion comparison (Fig. 1) and the asymptotic behavior of p-p rates (Fig. 2) are encouraging but do not quantitatively validate the BTE setup. Without this baseline, the reader cannot judge whether the computed absolute values, and hence the 50% reduction, are trustworthy.","section":"Sec. III (Figs. 1, 2, and text)"}],"minor_comments":[{"comment":"The phrasing 'Declining the lattice thermal conductivity' is unconventional; 'Reducing' or 'Suppressing' would be clearer.","section":"Abstract"},{"comment":"There is a typo: 'symmetery' should be 'symmetry'.","section":"Fig. 1 caption"},{"comment":"The comparison with Ref. [6] is described only as 'matches well'; more quantitative detail (e.g., a range of T over which agreement holds) would strengthen the validation.","section":"Sec. III (Fig. 3)"},{"comment":"The discrete scattering rates for EF = 0.2 eV are attributed to intervalley scattering at the K point, but the figure is not fully described; a few sentences explaining the origin and numerical representation of these discrete features would improve readability.","section":"Sec. III (Fig. 4)"},{"comment":"The phrase 'first-principle calculations' should be 'first-principles calculations' for consistency.","section":"Sec. II and III"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a straightforward computational study applying the standard EPW + ShengBTE pipeline to add electron-phonon scattering to graphene's lattice thermal conductivity. It is not a breakthrough, but it does produce a concrete new prediction: at 200 K, moving the Fermi level from 0.2 eV to 1.8 eV reduces the lattice kappa by roughly half. That gives an experimentally accessible external knob—gate voltage—for tuning heat transport, so the result is worth taking seriously.\n\nWhat the paper does well: no post-hoc fitting anywhere. The phonon dispersion matches inelastic X-ray data, the anharmonic scattering rates have the expected q-dependence, and the calculated electrical resistivity at n = 2.86e13 cm^-2 reproduces the earlier literature. The trend of electron-phonon coupling with Fermi level follows the density of states in a sensible way. The pipeline is standard and the execution looks careful.\n\nThe soft spots are real but not fatal. First, the headline 50% reduction is not against pristine graphene; it is the difference between EF = 0.2 eV and EF = 1.8 eV. The paper never shows an undoped kappa curve. That leaves the abstract's claim unmoored and makes it impossible to say how much of the reduction is actually \"doping\" versus just a reference-point artifact. Easy fix, but it should have been in the first version.\n\nSecond, there are no convergence tests. The prediction depends on phonon-electron rates becoming comparable to anharmonic rates at low frequencies, yet the paper reports no checks on q-grid, k-grid, smearing, or supercell size. For a computed prediction, this is an omission that a referee should ask about.\n\nThird, and most substantive, is the rigid-band approximation at high doping. The maximum reduction occurs at EF = 1.8 eV, which is exactly the LDA Van Hove singularity. At that carrier density—a few times 10^14 cm^-2—screening and self-consistent band renormalization could shift or broaden the DOS peak. The paper does not discuss this risk, does not test it with a self-consistent doped calculation, and does not report whether the 50% survives if the singularity is weakened. The stress-test note is on target: the central claim is largest exactly where the approximation is least secure. That does not mean the claim is wrong, but it is less solid than the paper lets on.\n\nMinor wording aside (\"externally induced\" in the abstract is a bit generous for a rigid shift), the citation pattern is appropriate: self-citation of ShengBTE is justified, and the prior-work list covers the relevant defect and isotope studies.\n\nThis paper deserves a serious referee. A good referee will ask for the pristine kappa curve, convergence checks, and an honest treatment of the rigid-band limit at high doping. If those come back clean, the prediction is a useful, citable result for thermoelectric applications and 2D thermal transport. I would send it out.","headline":"A clean, plausible first-principles prediction that charge doping can cut graphene's lattice thermal conductivity by ~50% at 200 K, but the missing pristine baseline and unexamined rigid-band assumption at high doping are the real caveats.","tokens_in":5875,"tokens_out":2171,"would_cite":true,"duration_ms":26723,"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":"Charge doping can cut suspended graphene's lattice thermal conductivity by about half at 200 K.","keywords":["graphene","lattice thermal conductivity","charge doping","electron-phonon scattering","phonon Boltzmann transport equation","relaxation-time approximation","first-principles calculation","thermoelectric application"],"falsifier":"Measure the lattice thermal conductivity of a suspended graphene device at 200 K while sweeping a back gate through the density range corresponding to Fermi energies between roughly 0.2 and 2 eV; if the conductivity does not fall by roughly half and pass through a minimum near $E_F \\approx 1.8$ eV, the central claim is contradicted.","tokens_in":4908,"feed_emoji":"⚡","tokens_out":6380,"duration_ms":63812,"temperature":0.7,"pith_summary":"This paper aims to show that charge doping alone can substantially reduce the lattice thermal conductivity of suspended graphene. In undoped graphene, phonon–phonon anharmonic scattering dominates heat-flow resistance, but the authors argue that once the Fermi energy is pushed into the conduction band by an external gate, phonon–electron scattering becomes comparable to phonon–phonon scattering. Combining first-principles force constants with electron–phonon scattering rates in the phonon Boltzmann equation, they find a roughly 50 percent drop in lattice thermal conductivity at 200 K as the Fermi energy moves from 0.2 eV to 1.8 eV above the Dirac point. Because doping is reversible and contact-free, this would give a practical electrical knob for thermoelectric devices made of graphene.","feed_headline":"Charge doping can cut graphene's heat conduction by half","feed_subtitle":"At 200 K, gate-controlled doping can shrink graphene's lattice thermal conductivity by about 50 percent.","key_machinery":"The load-bearing object is the phonon–electron scattering rate computed with the standard relaxation-time model for electron–phonon coupling and added to the phonon Boltzmann transport equation, alongside three-phonon anharmonic rates obtained from second- and third-order force constants. Shifting the Fermi energy rigidly into the unoccupied states changes only the electron occupation factors and the phase space for scattering, so the same first-principles electron–phonon matrix elements produce larger scattering rates as the density of states at the Fermi level grows. The peak coupling near $E_F = 1.8$ eV is what produces the maximum reduction in phonon lifetimes and hence in $\\kappa_{\\mathrm{ph}}$.","core_discovery":"The central claim is that externally induced electron doping can cut the lattice thermal conductivity of suspended graphene by about half at 200 K. The reduction is not monotonic: the conductivity falls as the Fermi energy rises, reaches a minimum near $E_F = 1.8$ eV, and then recovers, because the electron–phonon coupling strength tracks the density of states, which peaks near that energy. At high doping the phonon–electron scattering rates become comparable to the intrinsic anharmonic phonon–phonon rates, especially for the low-frequency acoustic phonons that carry most of the heat. The authors verify their electron–phonon scattering rates by comparing the computed electrical resistivity with published data, and they find the doping effect is strongest at lower temperatures and weakens as temperature increases.","pith_inferences":["A natural testable extension is to measure the thermal conductivity of a suspended graphene device versus back-gate voltage at fixed temperature; the predicted non-monotonic minimum near $E_F \\approx 1.8$ eV would distinguish this mechanism from ordinary impurity scattering.","If the rigid-band shift is replaced by a full self-consistent band structure at high doping, the optimal Fermi energy and the 50 percent figure could shift; comparing the two would show how much of the result depends on the rigid-band assumption.","The paper's emphasis on electron doping leaves hole doping as a nearby asymmetry: because the coupling strength is predicted to be smaller for holes at high Fermi energies, a smaller thermal-conductivity reduction would support the electron–phonon mechanism rather than a generic doping effect."],"forward_implications":["Gating becomes a reversible, contact-free way to lower graphene's lattice thermal conductivity, complementing defect- or isotope-based reduction strategies.","The same phonon Boltzmann framework predicts that the relative reduction is larger at lower temperatures, so 200 K is not the optimum but a representative operating point.","Because the effect peaks near $E_F = 1.8$ eV, experiments should target carrier densities near that Fermi-level shift rather than monotonically higher doping.","The mechanism should transfer to other two-dimensional materials whose Fermi level can be tuned electrostatically, provided their electron–phonon coupling also grows with doping."],"supporting_citations":[{"why":"supplies the experimentally demonstrated gate-tunable electron density that justifies shifting $E_F$ up to about 2.7 eV","marker":"[5]"},{"why":"provides the published electrical-resistivity data used to validate the computed electron–phonon scattering rates","marker":"[6]"},{"why":"establishes that graphene's electron–phonon coupling strength increases with Fermi energy and is larger for electron doping at high energies","marker":"[9]"},{"why":"supplies the phonon Boltzmann transport equation solver used to compute lattice thermal conductivity from interatomic force constants and scattering rates","marker":"[17]"},{"why":"provides the electron–phonon Wannier interpolation code used to obtain the scattering rates on fine grids","marker":"[18]"},{"why":"gives the relaxation-time formula used to turn electron–phonon matrix elements into phonon–electron scattering rates","marker":"[20]"},{"why":"supplies experimental phonon-dispersion data used to validate the harmonic force constants","marker":"[22]"}],"fun_headline_variants":["Charge doping halves graphene lattice heat flow","Doping cuts graphene's thermal conductance by half","Electron doping reduces graphene heat conduction 50%","Graphene heat conduction cut in half by doping"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that rigidly shifting the Fermi energy into the conduction band, leaving the band structure and electron–phonon matrix elements unchanged, accurately gives the phonon–electron scattering rates up to about 1.8 eV above the Dirac point; if that rigid-band picture or the relaxation-time model's coupling to the low-energy flexural phonons is inaccurate, the 50 percent reduction would shrink or move.","fun_headline_variants_meta":{"raw":{"variants":["Charge doping halves graphene lattice heat flow","Doping cuts graphene's thermal conductance by half","Electron doping reduces graphene heat conduction 50%","Graphene heat conduction cut in half by doping"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000283,"raw_usage":{"total_tokens":1586,"prompt_tokens":777,"completion_tokens":809,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":393,"completion_tokens_details":{"reasoning_tokens":750}},"tokens_in":393,"tokens_out":809,"duration_ms":8204,"temperature":1.0,"reasoning_tokens":750,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:44:30.107454+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the lattice thermal conductivity of a suspended graphene device at 200 K while sweeping a back gate through the density range corresponding to Fermi energies between roughly 0.2 and 2 eV; if the conductivity does not fall by roughly half and pass through a minimum near $E_F \\approx 1.8$ eV, the central claim is contradicted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"establishes that graphene's electron–phonon coupling strength increases with Fermi energy and is larger for electron doping at high energies"},{"cited_title":"Significant reduction of lattice thermal conductivity in suspended graphene by charge doping","cited_arxiv_id":"1908.02433","evidence_quote":"supplies the phonon Boltzmann transport equation solver used to compute lattice thermal conductivity from interatomic force constants and scattering rates"},{"cited_title":"Li, Phys","cited_arxiv_id":null,"evidence_quote":"gives the relaxation-time formula used to turn electron–phonon matrix elements into phonon–electron scattering rates"}],"review_version":1}