{"id":"39cc0b1f-fa9b-4203-97e6-994ec3f1a2c9","arxiv_id":"2505.21234","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Room-temperature phonon-limited mobilities are predicted as roughly 10^3 to 10^4 cm2/Vs in graphyne and about 10^2 cm2/Vs in graphdiyne using the DFTBephy workflow.","lead":"Using a fast approximate quantum method called DFTB, this paper computes how strongly electrons and holes scatter from atomic vibrations in two carbon sheet materials, graphyne and graphdiyne, and converts those scattering rates into charge-carrier mobilities. The numbers matter for choosing semiconducting 2D carbon materials for nanoelectronics, and the method is cheap enough to screen other large-unit-cell materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper never validates DFTB electron-phonon matrix elements against DFT; given the 2-3 order-of-magnitude graphdiyne mobility discrepancy, the claimed quantitative accuracy of the SERTA mobilities is not established.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing concern: DFTB EPC matrix elements are never benchmarked against DFT, so the accuracy claim is unsupported, and the graphdiyne discrepancy is the concrete symptom of this gap. My stress-test pass confirms this reading. I considered whether the concern is merely 'outside current consensus' or 'internally inconsistent'; it is neither. The paper is internally coherent, and the disagreement with DFT-based literature is acknowledged in Table I, but the authors' dismissive explanation (zeroth-order deformation-potential theory) is not tested. The proposed test is concrete and feasible: DFT EPC calculations for graphynes are available in the literature, and a direct comparison would settle whether the DFTB workflow is quantitatively reliable or only qualitatively indicative. Independent support in the paper—good structural, elastic, and phonon benchmarks, reproducible code availability, and a clear workflow—supports the qualitative conclusions and the methodological demonstration, but it does not establish the quantitative accuracy of the SERTA mobilities for graphdiyne. Therefore, the conditional verdict remains appropriate: the paper should be accepted with the condition that the EPC accuracy be demonstrated or the claims be softened to qualitative predictions.","tokens_in":17032,"tokens_out":2240,"duration_ms":26291,"concrete_test":"Compute the EPC matrix elements and SERTA scattering rates for graphyne and graphdiyne with a DFT-based method (e.g., EPW or Quantum ESPRESSO using PBE and comparable k/q grids), and compare g_mn^lambda and tau_SERTA directly to the DFTBephy values. If the DFT-based tau for graphdiyne remains close to 0.04/0.14 ps, the discrepancy lies in Ref. [34] and the claim survives; if the DFT-based tau is one to two orders of magnitude larger, the mio-1-1 EPC is the culprit and the central accuracy claim fails. As a cheaper first step, run the same DFTB-versus-DFT EPC comparison on graphene, where reliable DFT EPC results already exist, to isolate whether the DFTB Hamiltonian and overlap gradients are quantitatively correct for sp2/sp carbon.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the DFTBephy workflow computes electron-phonon couplings and phonon-limited mobilities efficiently and accurately. The supporting benchmarks cover lattice constants, elastic moduli, phonons, and band edges (Section III, Figs. 3-5), but there is no direct comparison between DFTB and DFT for the EPC matrix elements g_mn^lambda(k,q) or for the resulting SERTA scattering rates. This is the load-bearing gap because EPC depends on the gradients of the Hamiltonian and overlap matrices supplied by the mio-1-1 Slater-Koster set, and on the electronic eigenstates, which differ substantially from DFT: the DFTB band gaps are 1.38 eV (GY) and 1.52 eV (GDY) versus 0.46 eV and 0.48 eV from DFT-PBE, and the DFTB effective masses are systematically larger (Section III, Table II). Without a direct EPC benchmark, the graphyne relaxation-time agreement with Ref. [36] (within 21-54%) does not demonstrate that the EPC matrix elements are correct; it could be fortuitous. The graphdiyne discrepancy in Table I is more serious: DFTBephy gives tau = 0.04 ps (holes) and 0.14 ps (electrons), while the DFT-based values from Ref. [34] are 17.49 ps and 1.91 ps, and the mobilities differ by two to three orders of magnitude. The manuscript attributes this to zeroth-order deformation-potential theory in Ref. [34], but this explanation is not tested. If instead the DFTB EPC is inaccurate—for example, due to the use of mio-1-1 without self-consistent charges or due to incorrect Hamiltonian/overlap gradients for sp/sp2 carbon—then all reported SERTA mobilities are called into question. The paper's own conclusions therefore rest on an unvalidated premise, and the mixed transport-level benchmark does not resolve it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a density-functional tight-binding (DFTB) based workflow, implemented in the DFTBephy package, for computing electron-phonon couplings (EPCs) and phonon-limited charge carrier mobilities in two γ-graphyne polymorphs, graphyne (GY) and graphdiyne (GDY). The authors benchmark DFTB against DFT-PBE for lattice constants, elastic moduli, phonon dispersions, and band structures; they then compute relaxation times within the self-energy relaxation-time approximation (SERTA) and mobilities within SERTA, the constant relaxation-time approximation (CRTA), and analytical parabolic- and Kane-band models. They report SERTA relaxation times of 0.63 ps (holes) and 1.69 ps (electrons) for GY, and 0.04 ps and 0.14 ps for GDY, with mobilities of order 10^3–10^4 cm2/Vs for GY and 10^2 cm2/Vs for GDY. The central claim is that the DFTB approach is both computationally efficient and quantitatively accurate for EPC and transport in these materials.","tokens_in":17468,"tokens_out":2970,"duration_ms":34158,"significance":"If the reported mobilities are correct, the paper would demonstrate that a tight-binding-based workflow can replace computationally expensive DFT-based EPC calculations for large 2D carbon allotropes, substantially lowering the cost of transport predictions for materials with large unit cells. The paper is commendable for using a reproducible, open-source workflow and for providing access to the data via Zenodo. The benchmarks for structural, elastic, and phonon properties are useful in their own right. However, the transport-level claim of 'accurate' is not yet established: the EPC matrix elements are never compared against a DFT reference, and the GDY mobilities differ from prior DFT-based calculations by two to three orders of magnitude. The significance of the paper therefore hinges on whether this validation gap can be closed.","major_comments":[{"comment":"","section":"Section III, Fig. 3–5 and Table I"},{"comment":"","section":"Section II.B and II.E; Table I"},{"comment":"","section":"Section III (Electronic structure) and Table II (SI)"}],"minor_comments":[{"comment":"","section":"Abstract and title"},{"comment":"","section":"Section II.D, Eq. (10)"},{"comment":"","section":"Section II.E, computational details"},{"comment":"","section":"Section III, scattering rates"},{"comment":"","section":"Section III, mobilities"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be of interest to the computational materials community, but the absence of any direct validation of the DFTB EPC matrix elements against a DFT reference is a serious gap, especially given the order-of-magnitude discrepancy for GDY. I would encourage the editor to request a revision that either adds such a benchmark or substantially tempers the 'accurate' claim. The provided data availability statement and open-source workflow are strengths that should be highlighted in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this: the paper is a solid, useful application of the authors' previously published DFTBephy method to gamma-graphyne and graphdiyne, but the central claim that it yields \"accurate\" SERTA mobilities is not actually demonstrated. The specific relaxation times and mobility curves are new, but the method is prior work, and the graphyne numbers land close to earlier Wannier-DFT results while the graphdiyne numbers are off by two to three orders of magnitude.\n\nWhat the paper does well: the structural, elastic, and phonon benchmarks against DFT are thorough and look good. The SI includes a sensible test of several Slater-Koster parametrizations, and the code and data are deposited, so the results are reproducible. The analytical Kane-band comparison is a nice sanity check and gives physical insight into how nonparabolicity and energy-dependent scattering affect the mobility.\n\nThe soft spot is exactly where the stress-test puts it. The entire transport calculation rests on DFTB electron-phonon matrix elements, but the paper never compares those matrix elements or the resulting scattering rates against a DFT-based reference. It validates lattice constants, phonons, band edges, and effective masses, but EPC depends on Hamiltonian gradients and eigenstates, and those are known to differ: DFTB band gaps are about 1 eV larger than PBE, and the DFTB effective masses are 2-4 times larger. The GY relaxation times agree with Ref. [36] within 20-54%, but that agreement could be coincidental, with a too-large mass and an inaccurate EPC partly canceling. The GDY discrepancy is not a minor detail: tau = 0.04 ps instead of 17.49 ps for holes. Blaming zeroth-order deformation-potential theory in Ref. [34] without testing that alternative is hand-waving. I also missed convergence tests for the k- and q-meshes and any error bars on the mobilities.\n\nDespite those reservations, this is a coherent, well-written computational study. The authors are up front about the GDY discrepancy and about the band-gap differences. It deserves a serious referee, but the referee should push for a direct DFT EPC benchmark—even for one of the two materials or a smaller model system—and for a quantitative resolution of the GDY mismatch. With those additions, this would be a dependable standard reference for DFTB-based transport in 2D carbon allotropes.","headline":"A competent application of the authors' own DFTBephy code to graphynes, but the quantitative accuracy claim is under-supported because the electron-phonon matrix elements are never validated against DFT.","tokens_in":786,"tokens_out":759,"would_cite":true,"duration_ms":32278,"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 claims that a density-functional tight-binding (DFTB) workflow, implemented in the DFTBephy package, can compute electron-phonon couplings and phonon-limited charge carrier mobilities for gamma-graphyne and gamma-graphdiyne at a…","keywords":["graphyne","graphdiyne","electron-phonon coupling","charge carrier mobility","density functional tight binding","SERTA","Boltzmann transport equation","Kane band model"],"falsifier":"A DFT-based electron-phonon coupling calculation for the valence and conduction band edges of graphyne, compared head-to-head with the DFTB matrix elements, would confirm or refute this premise; if the resulting SERTA mobilities differ by more than a factor of two, the DFTB assumption fails.","tokens_in":16859,"feed_emoji":"⚡","tokens_out":8494,"duration_ms":75343,"temperature":0.7,"pith_summary":"The paper claims that a density-functional tight-binding (DFTB) workflow, implemented in the DFTBephy package, can compute electron-phonon couplings and phonon-limited charge carrier mobilities for gamma-graphyne and gamma-graphdiyne at a fraction of the cost of full DFT. Using the self-energy relaxation-time approximation (SERTA), it reports hole (electron) relaxation times of 0.63 ps (1.69 ps) for graphyne and 0.04 ps (0.14 ps) for graphdiyne at room temperature, with mobilities reaching roughly $10^4$ cm$^2$/Vs for electrons in graphyne and around $10^2$ cm$^2$/Vs for both carriers in graphdiyne. These values place graphyne between graphene and MoS2 and graphdiyne close to MoS2. The study also compares constant relaxation-time, parabolic-band, and Kane-band analytical models to show where non-parabolicity and energy-dependent scattering matter.","feed_headline":"Graphyne electrons hit 10⁴ cm²/Vs","feed_subtitle":"Fast DFTB workflow maps phonon-limited mobility for large-carbon 2D sheets.","key_machinery":"The central machinery is the DFTBephy workflow, which builds electron-phonon coupling matrix elements from finite-difference gradients of the real-space Hamiltonian and overlap matrices with respect to atomic displacements, Fourier-transforms them, and combines them with phonon eigenvectors from Phonopy. The coupling matrix elements then feed the self-energy relaxation-time approximation (SERTA) scattering rate, Eqs. (6)-(7), and the Boltzmann conductivity formula, Eq. (8). A secondary piece is the Kane-band model, a non-parabolic dispersion correction parameterized by an alpha value per band, which is used for analytical mobility estimates.","core_discovery":"On its own terms, the paper establishes that the DFTBephy workflow, using the mio-1-1 parameter set without self-consistent charges, yields quantitatively reliable electron-phonon couplings for gamma-graphynes. The central quantitative results are SERTA relaxation times of 0.63 ps for holes and 1.69 ps for electrons in graphyne, and 0.04 ps and 0.14 ps in graphdiyne, giving room-temperature mobilities of about $3.5 \\times 10^3$ cm$^2$/Vs (holes) and $9.6 \\times 10^3$ cm$^2$/Vs (electrons) in graphyne, and about $4.9 \\times 10^2$ and $7.4 \\times 10^2$ cm$^2$/Vs in graphdiyne. The paper also finds that for holes, acoustic phonon scattering dominates near the band edge, while for electrons optical phonon scattering is roughly three times stronger in graphyne. It further claims that the Kane-band analytical model with an energy-dependent relaxation time reproduces the SERTA trend better than parabolic or constant-relaxation approximations.","pith_inferences":["If a DFT-based EPC benchmark were run for graphyne, it would settle whether mio-1-1's electron-phonon matrix elements are quantitative; the paper's own benchmarks cover geometry, phonons, and band edges but not EPC itself.","Assuming the DFTB couplings are accurate, the two-to-three-order-of-magnitude gap with the graphdiyne literature (Ref. [34]) indicates those earlier deformation-potential estimates overestimate mobility by omitting full phonon scattering.","The same workflow should transfer to other sp/sp2 carbon allotropes and large-unit-cell covalent organic frameworks, enabling high-throughput screening of phonon-limited transport.","Because DFTB band gaps here are about 1 eV larger than PBE, a hybrid-functional or GW correction could shift carrier densities and scattering phase space; testing this would reveal how sensitive the reported mobilities are to the gap error."],"forward_implications":["Graphyne's phonon-limited electron mobility reaches about $10^4$ cm$^2$/Vs at 300 K, between graphene and MoS2, making the material a plausible 2D transistor channel.","Graphdiyne's mobility, around $10^2$ cm$^2$/Vs, is comparable to MoS2, so its larger pores and lower stiffness do not push transport into a useless range.","In graphyne, optical phonon scattering is the dominant electron-scattering channel, so engineering that stiffens optical modes (e.g., strain or isotope substitution) could raise electron mobility.","The DFTB workflow makes electron-phonon coupling calculations tractable for large-unit-cell 2D carbon allotropes, where conventional DFT-based EPC is computationally prohibitive.","The Kane-band model with energy-dependent relaxation times tracks SERTA more closely than constant relaxation time, so it can serve as a low-cost screening tool for carrier-density-dependent mobility."],"supporting_citations":[{"why":"Supplies the DFTB-based electron-phonon coupling calculation scheme and the workflow used throughout.","marker":"[20]"},{"why":"Provides the mio-1-1 Slater–Koster parameter set that defines the DFTB Hamiltonian and overlap matrix elements.","marker":"[26]"},{"why":"Supplies the small-displacement phonon calculations and force constants used for the dynamical matrix and phonon eigenvectors.","marker":"[22]"},{"why":"Provides the DFTB+ code used to generate the real-space Hamiltonian and overlap matrices.","marker":"[24]"},{"why":"Gives the SERTA and Boltzmann transport formalism for relaxation times and conductivities.","marker":"[15]"},{"why":"Gives the comparison graphyne relaxation times and mobilities (the paper's values are 21-54% higher).","marker":"[36]"},{"why":"Gives the graphdiyne literature mobilities that the paper's SERTA values are two to three orders of magnitude below.","marker":"[34]"}],"fun_headline_variants":["DFTBephy computes phonon-limited mobility in graphynes","Graphyne electrons hit 10⁴ cm²/Vs via DFTBephy","Fast DFTB maps graphyne mobility: electrons 10⁴, holes 10³","Graphyne electron mobility: 10⁴ cm²/Vs, between graphene and MoS2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the mio-1-1 DFTB parameter set, without self-consistent charges, reproduces the electron-phonon matrix elements of graphynes quantitatively, an assumption the paper tests only indirectly through structural, phonon, and band-edge benchmarks and through mobility comparisons that disagree wildly for graphdiyne.","fun_headline_variants_meta":{"raw":{"variants":["DFTBephy computes phonon-limited mobility in graphynes","Graphyne electrons hit 10⁴ cm²/Vs via DFTBephy","Fast DFTB maps graphyne mobility: electrons 10⁴, holes 10³","Graphyne electron mobility: 10⁴ cm²/Vs, between graphene and MoS2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001036,"raw_usage":{"total_tokens":4430,"prompt_tokens":1085,"completion_tokens":3345,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":701,"completion_tokens_details":{"reasoning_tokens":3255}},"tokens_in":701,"tokens_out":3345,"duration_ms":28396,"temperature":1.0,"reasoning_tokens":3255,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:33:40.310610+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A DFT-based electron-phonon coupling calculation for the valence and conduction band edges of graphyne, compared head-to-head with the DFTB matrix elements, would confirm or refute this premise; if the resulting SERTA mobilities differ by more than a factor of two, the DFTB assumption fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the DFTB-based electron-phonon coupling calculation scheme and the workflow used throughout."},{"cited_title":"Ponc´ e, W","cited_arxiv_id":null,"evidence_quote":"Gives the SERTA and Boltzmann transport formalism for relaxation times and conductivities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the comparison graphyne relaxation times and mobilities (the paper's values are 21-54% higher)."}],"review_version":1}