{"id":"81141810-80d7-4377-8898-8fde1e39188c","arxiv_id":"2502.08538","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"First-principles DFT/EPW calculations extract acoustic deformation potentials for Si, diamond, and cBN; Si matches literature, diamond exceeds previous values, and cBN supports a prior estimate.","lead":"This paper computes acoustic deformation potentials for silicon, diamond, and cubic boron nitride using density functional theory and electron-phonon coupling calculations, then fits them to standard angular formulas. The values agree with known silicon data, diamond comes out higher than most earlier estimates, and the cBN value supports a prior estimate, which matters for simulating electron transport in new power-electronics materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central extraction presumes g/|q| is in the long-wavelength linear regime, but the paper never reports the q magnitudes or a convergence test; this is the load-bearing assumption behind every reported Xi_ADP.","rationale":"I read the paper as a method paper: extract acoustic deformation potentials from EPW electron-phonon matrix elements by evaluating Eq. 8 and fitting the angular forms in Eq. 10. The Si results (Fit 1 and Fit 2) and effective masses are a genuine validation, and the phonon dispersions match experiment. The central mathematical relationship, Eq. 9, is the load-bearing point. The reader's weakest assumption correctly identifies that the paper never reports the q magnitudes or demonstrates the small-q linear regime. I agree this is the most serious issue. If g/|q| is not constant over the q range used, then the fitted Xi_d, Xi_u, and Xi_ADP are artifacts of the chosen q sampling; because the q values are not stated, the results cannot be reproduced or assessed. The proposed q-plateau scan would settle this. I also note two secondary concerns that do not change the verdict: (i) the theoretical motivation in Eqs. 3-5 is written for the LA dilatation coupling, while the TA values are extracted from the same expression despite the fact that div(u) = 0 for a transverse wave; this deserves a clearer justification, although the numerical g/q extraction may still be valid if the shear coupling is linear in q. (ii) The cBN long-range subtraction is described in one sentence and the azimuthal scatter is large; a convergence test over azimuthal angles and q would improve confidence. Since the reader already made the verdict conditional on resolving the small-q issue, my stress-test leaves the verdict unchanged.","tokens_in":15168,"tokens_out":8779,"duration_ms":93545,"concrete_test":"For a single valley and fixed direction in Si (e.g., theta = pi/4, phi = 0), compute Xi(|q|) = sqrt(2 m0 omega_q / hbar) |g_nu(k,q)| / |q| for q magnitudes 0.005, 0.01, 0.02, 0.05, and 0.1 in units of 2*pi/a, using the same EPW settings. If the values do not form a plateau as |q| -> 0, or if a finite-difference derivative of sqrt(2 m0 omega_q / hbar) |g| with respect to |q| at the working q differs from the ratio, then Eq. 9 is not satisfied and all fit parameters change. Repeat this test for at least one polar angle per material; if the Si plateau holds but the diamond or cBN plateaus do not, the method only works for Si.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The method's central output is Xi_ADP obtained from Eq. 8 and Eq. 9: Xi = sqrt(2 m0 omega / hbar) |g| / |q| in the limit as |q| -> 0. All reported values for Si, diamond, and cBN depend on evaluating this ratio at finite q points selected by the EPW Wannier interpolation. The manuscript never states the actual |q| values, how many q points per polar angle were used, or whether g was verified to be linear in q in that range. If the q values are not small enough, the assumed identity g proportional to |q| fails: the fitted Xi_d and Xi_u in Eqs. 10a/10b, and hence Xi_LA, Xi_TA, and Xi_ADP via Eqs. 11-12, are not well-defined. This is not merely a missing error bar; it questions whether the extraction procedure produces a unique number at all without a specified q. The Si agreement with literature is reassuring, but it cannot substitute for a convergence test, because a single q value can accidentally match. The same issue affects the diamond and cBN results and the claimed predictive power for materials with sparse data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript develops a first-principles procedure for extracting acoustic deformation potentials (ADPs) for the conduction band minima of silicon, diamond, and cubic boron nitride. Using DFT and DFPT via Quantum Espresso and EPW, the authors compute electron-phonon coupling matrix elements at fixed polar and azimuthal angles and convert them to angle-dependent deformation potentials through Eq. 8 (with the definition in Eq. 9), fit the results to the Herring–Vogt angular forms in Eqs. 10a and 10b, average over polar angle via Eqs. 11a and 11b, and combine longitudinal and transverse contributions via Eq. 12. For Si, the fitted Ξ_ADP values (9.52 and 8.56 eV) agree with literature values; for diamond they obtain larger values (21.42 and 20.04 eV) than most empirical reports; for cBN they obtain 12.21 and 10.23 eV, which bracket the estimate used by Siddiqua et al. The central claim is that this approach provides ADP scattering parameters without requiring empirical mobility or time-of-flight fits.","tokens_in":15456,"tokens_out":3142,"duration_ms":33578,"significance":"If the method is sound, it offers a practical route to transport parameters for materials with sparse experimental data, supporting Monte Carlo and Boltzmann transport simulations. The work has several strengths: it does not use experimental mobilities to fix the deformation potentials, so the central extraction is not circular; it benchmarks against Si, a well-characterized material; the fits shown in Figs. 2, 4, and 6 capture the dominant angular trends; and the authors have deposited the data underlying the figures in a public repository. The main significance, however, is conditional on the method's operational validity, particularly whether the finite-q evaluation used in Eq. 8 is demonstrably in the linear regime and whether the polar correction for cBN is adequately specified. The large spread in literature values for diamond and the absence of experimental data for cBN make the predictive claim rest heavily on these technical details.","major_comments":[{"comment":"The manuscript never reports the magnitude of the phonon wavevectors |q| used in evaluating the ratio g/|q| in Eq. 8, nor does it present a convergence study showing that the computed values are in the small-q linear regime where Eq. 9 applies. Since every fitted value of Ξ_d and Ξ_u and hence every reported Ξ_ADP depends on this finite-q ratio, the extraction procedure is operationally underdefined without such a test. A single q point per polar angle could accidentally match literature values (as in Si), but this does not establish that the method yields a unique material parameter. I request a convergence test over a range of |q| values, or a plot showing g behaves linearly in |q| for the q values used, for at least one material.","section":"§II.B, Eq. 9, and §III.A–C"},{"comment":"For cBN, the only statement about the polar correction is that \"the long range effects were subtracted before the deformation potential was calculated.\" No details are given about the subtraction scheme, the Fröhlich parameters, how the q→0 limit was handled, or the sensitivity of the resulting acoustic matrix elements to the subtraction. This is load-bearing because the longitudinal acoustic deformation potential is extracted from the same matrix elements that require the polar correction, and the cBN result is used as evidence of predictive power. Please specify the polar-correction method used in the EPW calculations and quantify its effect on the reported Ξ_ADP.","section":"§III.C"},{"comment":"No uncertainties are reported for the fitted deformation potentials, despite the visible scatter in the angular data and the two fitting procedures yielding substantially different results (Fit 1 vs Fit 2: 21.42 vs 20.04 eV for diamond; 12.21 vs 10.23 eV for cBN). The azimuthal bars in Figs. 4 and 6 and the differences between fitting protocols should be propagated into the final Ξ_ADP values. Without error estimates, the claimed agreement for cBN with the Siddiqua et al. estimate (which the authors themselves note may be coincidental) cannot be quantitatively assessed, and the diamond discrepancy with literature remains only qualitatively interpreted.","section":"§III.B, §III.C, Tables II and III, Eqs. 11–12"},{"comment":"The angular averaging leading to Eqs. 11a and 11b assumes a uniform distribution of polar angles and an isotropic combination via Eq. 12, but the manuscript does not justify this averaging scheme as the physically appropriate one for scattering in these valleys. This is an assumption of the method rather than a derived result; showing that Eq. 12 reproduces the full first-principles scattering rate in a model test, or at least discussing the sensitivity of Ξ_ADP to the averaging choice, would strengthen the central claim.","section":"§II.A, Eqs. 10–12"}],"minor_comments":[{"comment":"The text states that the diamond conduction band minimum is located at \"roughly 73%\" of the distance between Γ and X, but the caption of Fig. 3(a) states the minimum is at \"82% along the Δ direction.\" Please reconcile this inconsistency.","section":"§III.B, text before Fig. 3"},{"comment":"The phrase \"The black dashed line is then the plot of equation 10b for Ξ_u = 16.4 eV\" is unexplained: neither Fit 1 (19.12 eV) nor Fit 2 (15.89 eV) gives 16.4 eV. Please correct this value or clarify which fit it corresponds to.","section":"§III.C, §II.A, Fig. 6"},{"comment":"The sentence describing the q-point grid says it was used \"in conjuncture with\" a 10×10×10 grid; this should be \"in conjunction with.\" There are also several typographical errors in the introduction (e.g., \"ab inito,\" \"catagorising\") that should be corrected.","section":"§II.B"},{"comment":"The table headings \"Fit 1\" and \"Fit 2\" are not self-explanatory without reading the text; please add a footnote defining the two fitting procedures in each table caption.","section":"§III, Table I and Table II"},{"comment":"Equation 9 defines the deformation potential as a derivative with respect to |q|, but the absolute value notation is applied only to the matrix-element prefactor and not to the matrix element itself. For clarity, specify whether the derivative is taken on the magnitude of the complex matrix element or on its real part.","section":"§II.A, Eq. 9"}],"recommendation":"major_revision","confidential_remarks":"The manuscript addresses a useful and timely topic, and the non-circular extraction from first-principles electron-phonon matrix elements is a genuine strength. However, the missing q-convergence analysis and the inadequate description of the polar correction for cBN are load-bearing gaps that prevent the method from being adopted as presented. The authors should be encouraged to add a convergence study and error estimates; with those additions, the paper would make a solid contribution. The paper's scope is appropriate for this journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does one thing and does it cleanly: it extracts acoustic deformation potentials for Si, diamond, and cBN from EPW electron-phonon matrix elements, fitting the full angular dependence rather than picking special points as Li et al. and Yang et al. did. The Si benchmark lands close to literature values, which lends the workflow initial credibility. The diamond and cBN numbers are new and directly useful for anyone running Monte Carlo transport simulations on ultra-wide-bandgap materials. Credit where due: the data is deposited, the figures are readable, and the text is refreshingly plain about what is assumed and what is uncertain. No circularity—no experimental mobility is fed into the fits.\n\nThe soft spots are real but not fatal. The most important is the one the stress-test highlights: equation 9 defines the ADP as a |q|→0 limit, but the paper never reports the actual q magnitudes, the number of q points, or any test that g/|q| is flat in the range used. That is not a nitpick. The entire output rests on that ratio being in the linear regime, and a single arbitrary q could shift every reported Ξ_d, Ξ_u, and Ξ_ADP. The Si agreement is reassuring, but it is not a substitute for a convergence test. There are also no error bars on any fitted parameter, and the cBN polar correction is described in one sentence. The cBN comparison to Siddiqua et al. is against a \"reasonable estimate,\" so saying \"excellent agreement\" overstates the validation. The diamond value sits well above all prior reports; the authors' explanation about free parameters in empirical fits is plausible, but they do not demonstrate it quantitatively.\n\nEven with these gaps, the central method is sensible and the paper does not overclaim. The data release and the Si benchmark show good faith, and the cBN result is explicitly flagged as an estimate support rather than a measurement. The likely fix is a short convergence study plus error bars and a bit more detail on the polar subtraction—none of that is conceptually hard.\n\nWho is this for? Anyone doing transport simulations of diamond, cBN, or other sparse-data UWBG materials who needs a defensible ADP input. A serious referee should see this; it is a legitimate incremental contribution, not a desk reject. My recommendation is to send it to peer review with the expectation of a revision that addresses the q-convergence issue and reports uncertainties.","headline":"A useful, honest incremental method paper with a real reproducibility gap: the finite-q extraction is never checked for convergence in the long-wavelength limit.","tokens_in":15947,"tokens_out":2201,"would_cite":false,"duration_ms":24662,"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":"Acoustic deformation potentials can now be computed purely from first principles, giving scattering parameters for semiconductors without experimental data.","keywords":["acoustic deformation potential","electron-phonon coupling","density functional perturbation theory","first principles","silicon","diamond","cubic boron nitride","charge transport simulation"],"falsifier":"Recompute the silicon, diamond, or cBN deformation potentials using phonon wavevectors whose magnitudes are repeatedly halved while keeping all other settings fixed; if the fitted $\\Xi_{\\mathrm{ADP}}$ from Eqs. (11) and (12) changes by more than a few percent as $|q|\\to 0$, the reported values are not converged and the comparison with literature is not yet settled.","tokens_in":14954,"feed_emoji":"💎","tokens_out":8086,"duration_ms":76742,"temperature":0.7,"pith_summary":"This paper tries to establish that the acoustic deformation potential—the material parameter controlling how electrons scatter off lattice vibrations—can be obtained from first-principles calculations instead of by fitting to experimental transport data. The authors compute the electron-phonon coupling matrix with density-functional theory and density-functional perturbation theory, evaluate the deformation potential at twenty polar angles between electron and phonon wavevectors, and fit the standard angular expressions to get angle-averaged values. For silicon the result (8.56–9.52 eV) matches the accepted range; for diamond it is larger than prior empirical values (20.0–21.4 eV versus 8.7–17.7 eV); for cubic boron nitride it falls near the estimate used in earlier simulations (10.2–12.2 eV). If the claim is right, transport simulations no longer need experimental scattering parameters for every new material.","feed_headline":"Acoustic deformation potentials now come from first principles","feed_subtitle":"A DFT/DFPT method yields scattering parameters for materials with no experimental data, matching silicon and supporting cBN estimates.","key_machinery":"The central object is the ratio $\\sqrt{2m_0\\omega/\\hbar}\\, g^\\nu_{m,n}(k,q)/|q|$, whose $|q|\\to 0$ limit is the acoustic deformation potential. That ratio is computed from DFT/DFPT electron-phonon matrix elements via Wannier-Fourier interpolation. The angle-dependent fits of $\\Xi_{\\mathrm{LA}}(\\theta)$ and $\\Xi_{\\mathrm{TA}}(\\theta)$ to the computed points separate $\\Xi_d$ and $\\Xi_u$, and Eqs. (11)–(12) average them into the scalar $\\Xi_{\\mathrm{ADP}}$ that enters the acoustic scattering rate.","core_discovery":"In the paper's own terms, the acoustic deformation potential at a band edge is the small-$|q|$ gradient of the electron-phonon coupling matrix element, so it can be read directly out of DFT/DFPT outputs. The authors calculate $\\Xi_{\\mathrm{adp}}$ from Eq. (8) at twenty polar angles with azimuthal averaging, fit the longitudinal and transverse angular dependences $\\Xi_{\\mathrm{LA}}(\\theta)=\\Xi_d+\\Xi_u\\cos^2\\theta$ and $\\Xi_{\\mathrm{TA}}(\\theta)=\\Xi_u\\sin\\theta\\cos\\theta$, and combine the fitted $\\Xi_d$ and $\\Xi_u$ via Eqs. (11) and (12) into a single $\\Xi_{\\mathrm{ADP}}$. The silicon numbers reproduce the literature, the cBN numbers support the estimate used by earlier workers, and the diamond numbers come out higher than the empirical range, which the paper attributes to the use of deformation potentials as free parameters in prior fits.","pith_inferences":["The paper leaves the q-convergence question open: checking $g/|q|$ at successively smaller phonon wavevectors would confirm or overturn the reported values.","The large azimuthal variation seen for cBN hints that the angle-averaged $\\Xi_{\\mathrm{ADP}}$ may depend on how many azimuthal directions are sampled, so a finer azimuthal grid is a natural robustness test.","If the diamond value near 20–21 eV is used in a mobility simulation, the predicted mobility will be lower than simulations using the empirical 8.7–17.7 eV range, making the discrepancy directly testable.","The same extraction could be applied to hole bands and to additional conduction valleys, where the deformation potential may differ, to build complete first-principles scattering tables."],"forward_implications":["Monte Carlo and Boltzmann-transport simulations can obtain acoustic scattering parameters from DFT/DFPT alone, removing the need to fit to time-of-flight or mobility measurements.","The method gives concrete first-principles $\\Xi_{\\mathrm{ADP}}$ values for diamond ($\\approx$20–21 eV) and cBN ($\\approx$10–12 eV) that can be tested in transport simulations.","Because the same pipeline works for silicon, diamond, and cBN, the approach should transfer to other diamond-structure and zinc-blende semiconductors with sparse experimental data.","The silicon and cBN agreement supports the paper's conclusion that the wide spread of diamond deformation potentials in the literature comes largely from empirical free-parameter fitting."],"supporting_citations":[{"why":"Defines the deformation-potential concept that ties band-edge shifts to strain.","marker":"[9]"},{"why":"Supplies the angular dependences of the longitudinal and transverse deformation potentials used to separate $\\Xi_d$ and $\\Xi_u$.","marker":"[10]"},{"why":"Supplies the Wannier-Fourier interpolation that produces the electron-phonon matrix elements from DFT/DFPT.","marker":"[25]"},{"why":"Prior first-principles extraction of silicon deformation potentials used as a comparison baseline.","marker":"[38]"},{"why":"Prior first-principles silicon calculation emphasizing transverse acoustic phonon contributions, also used as a comparison.","marker":"[39]"},{"why":"Provides the cBN acoustic deformation potential estimate that the paper's cBN result supports.","marker":"[61]"},{"why":"Provides empirical diamond uniaxial and dilatation deformation potentials that the paper compares against and finds larger than.","marker":"[17]"}],"fun_headline_variants":["First-principles acoustic deformation potentials for Si, diamond, cBN","DFT computes deformation potentials without empirical fits","Deformation potentials from first principles, no experiments needed","Acoustic deformation potentials: DFT matches Si, estimates cBN","First principles yield deformation potentials for novel materials"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the finite phonon wavevectors used in the calculation are already small enough that $g/|q|$ is in its linear regime; the paper does not report the q magnitudes or a convergence test for this regime.","fun_headline_variants_meta":{"raw":{"variants":["First-principles acoustic deformation potentials for Si, diamond, cBN","DFT computes deformation potentials without empirical fits","Deformation potentials from first principles, no experiments needed","Acoustic deformation potentials: DFT matches Si, estimates cBN","First principles yield deformation potentials for novel materials"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000155,"raw_usage":{"total_tokens":1223,"prompt_tokens":964,"completion_tokens":259,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":182}},"tokens_in":580,"tokens_out":259,"duration_ms":3225,"temperature":1.0,"reasoning_tokens":182,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T04:42:15.191197+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the silicon, diamond, or cBN deformation potentials using phonon wavevectors whose magnitudes are repeatedly halved while keeping all other settings fixed; if the fitted $\\Xi_{\\mathrm{ADP}}$ from Eqs. (11) and (12) changes by more than a few percent as $|q|\\to 0$, the reported values are not converged and the comparison with literature is not yet settled.","supporting_citations":[{"cited_title":"Takahashi, A","cited_arxiv_id":null,"evidence_quote":"Defines the deformation-potential concept that ties band-edge shifts to strain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the angular dependences of the longitudinal and transverse deformation potentials used to separate $\\Xi_d$ and $\\Xi_u$."},{"cited_title":"Tsukioka, Scattering mechanisms of carriers in natu ral diamond (2001)","cited_arxiv_id":null,"evidence_quote":"Supplies the Wannier-Fourier interpolation that produces the electron-phonon matrix elements from DFT/DFPT."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior first-principles extraction of silicon deformation potentials used as a comparison baseline."},{"cited_title":"Kittel, Introduction to Solid State Physics , 8th ed","cited_arxiv_id":null,"evidence_quote":"Prior first-principles silicon calculation emphasizing transverse acoustic phonon contributions, also used as a comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the cBN acoustic deformation potential estimate that the paper's cBN result supports."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides empirical diamond uniaxial and dilatation deformation potentials that the paper compares against and finds larger than."}],"review_version":1}