{"id":"d8e8e7ae-463a-43c1-9d03-bc765a7ad580","arxiv_id":"2506.06198","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"A corrected Stratton-based theory and a fast numerical implementation (GETELEC-2.0) for thermal-field emission from semiconductors, with predictions about band-bending-dependent current saturation and electron energy distributions.","lead":"This paper presents corrected equations and new software for calculating how electrons escape from semiconductor surfaces under strong electric fields. The authors use it to explain why valence-band electrons are rarely seen in experiments and to predict non-linear current-voltage behavior of semiconductor field emitters.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The valence-band correction rests on an integration-by-parts identity that is false as written: Eq. (11)'s lFD gives -lFD' = (1/kT)(1-fFD), not fFD, so Eq. (15) predicts zero valence current at low T. This must be fixed before the Stratton correction is claimed.","rationale":"The reader's weakest assumption concerned the imposed band-bending scenarios, which is a legitimate limitation of the applications but is not the most load-bearing issue for the paper's central claim. The paper's headline theoretical result is Eq. (15), a correction to Stratton's valence-band emission equation. That result depends on an integration-by-parts identity in the Appendix. As printed, the identity is false: with lFD defined by Eq. (11), -lFD' is not the Fermi function, and the resulting Eq. (15) would suppress valence-band emission at low temperature, contradicting the paper's own Figs. 5 and 8. This is an internal inconsistency in the derivation, not merely a disagreement with previous literature. It is concrete, checkable, and directly connected to the claimed correction. The concern does not require assuming bad faith; it may be a correctable sign error, but it must be corrected and re-verified before the central claim can be accepted. The reader's conditional verdict remains appropriate: the paper is not acceptable in its current form, but a focused revision could resolve the issue. I therefore recommend keeping the verdict conditional rather than moving to reject or accept, and I disagree with the reader's choice of the weakest assumption because a more fundamental problem sits in the derivation of Eq. (15) itself.","tokens_in":12257,"tokens_out":21158,"duration_ms":217286,"concrete_test":"Take the i-Ge parameters of Section IV and set T -> 0. Evaluate the valence-band current two ways: (a) directly from the transformed double integral jV = L ∫_{-∞}^{EV} fFD(E) ∫_{αbar E - αEV}^{E} D(Ez) dEz dE with fFD(E) = 1/(1+e^{(E-Ef)/kT}) by numerical quadrature; (b) from Eq. (15) using lFD exactly as defined in Eq. (11). If (b) is essentially zero while (a) is large because the valence band is fully occupied, the central identity fails. As an analytic check, verify whether -d/dE ln(1+e^{(E-Ef)/kT}) equals fFD(E); it does not. Then repeat the comparison against Stratton's Eq. (A.61) to see whether the extra ln term is actually required.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (15) is the central theoretical claim: it removes a ln term from Stratton's valence-band formula and underpins the cryogenic valence-emission predictions. The Appendix derivation (Eqs. A.25, A.55) replaces the Fermi occupation fFD(E) by -lFD'(E). With lFD defined in Eq. (11) as -ln(1/(1+e^{(E-Ef)/kT})) = ln(1+e^{(E-Ef)/kT}), one gets -lFD'(E) = -(1/kT) e^{(E-Ef)/kT}/(1+e^{(E-Ef)/kT}) = -(1/kT)(1-fFD(E)), not fFD(E). The identity used in the integration by parts is therefore algebraically false as written. If the intended supply function is the standard N(E)=kT ln(1+e^{(Ef-E)/kT}), whose derivative is -fFD(E), then Eq. (11) has the wrong sign in the exponent; that is not a harmless typo. For the i-Ge parameters in Section IV (EV = -5.1 eV, EF = -4.75 eV), lFD(E) = ln(1+e^{(E-Ef)/kT}) is exponentially small for all E below EV at low temperature, so Eq. (15) gives essentially zero valence-band current, directly contradicting the paper's own conclusion that valence-band emission dominates at cryogenic temperatures. Since the difference between Eq. (A.60) and Stratton's Eq. (A.61) is exactly the boundary term of this integration by parts, the claimed removal of the extra ln term is not established until the sign convention and boundary term are corrected and checked.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reviews Stratton's thermal-field emission theory for semiconductors, claims to correct the valence-band emission equation, and presents GETELEC-2.0, an extension of an existing numerical emission code. The new code evaluates conduction- and valence-band emission integrals, Nottingham heat, and electron energy distributions using a precomputed fourth-order polynomial fit to the JWKB Gamow function. The authors validate this interpolation against full numerical integration, simulate intrinsic Ge under three imposed band-bending scenarios, and discuss field, temperature, and Nottingham-heat dependences. The central theoretical claim is that Eq. (15) removes an erroneous logarithmic boundary term from Stratton's valence-band formula, which they argue explains the scarcity of observed valence-band emission and predicts valence-band-dominated emission at cryogenic temperatures.","tokens_in":12598,"tokens_out":9159,"duration_ms":95340,"significance":"If the corrected valence-band integral is rigorously established, the paper would provide a useful correction to a long-standing formula and a practical computational tool for semiconductor field-emission simulation, including energy distributions and Nottingham heating. The numerical strategy of tabulating polynomial fits to the Gamow function is transparent and is validated over a broad parameter range (RMSE 0.01531 for current density and 0.02987 for Nottingham heat), and the authors explicitly acknowledge the main physical limitations: the deferred self-consistent band-bending calculation, the assumed replacement-electron energy, and the neglect of surface states. However, the headline theoretical derivation currently contains an algebraic sign error in the supply function that, as written, makes the central valence-band result false at cryogenic temperatures, so the significance cannot be assessed until the derivation and all subsequent figures are corrected.","major_comments":[{"comment":"The supply function is defined with the wrong sign. Eq. (11) gives lFD(E) = ln(1 + e^((E-Ef)/kT)), whose derivative is lFD'(E) = (1/kT)(1 - fFD(E)). Therefore the identity used in Eqs. (A.25) and (A.55), namely that the emitted-current integrand can be written as -lFD'(E) g(E), is algebraically false: -lFD'(E) = -(1/kT)(1 - fFD(E)), not fFD(E). The function whose derivative is -fFD(E) is ln(1 + e^((Ef-E)/kT)), i.e., the standard supply function with the opposite sign in the exponent. With the printed definition, Eq. (15) gives essentially zero valence-band current at low temperature for E below EV, directly contradicting the paper's claim that valence-band emission dominates at cryogenic temperatures. In addition, the integration by parts in Eq. (A.27) and Eq. (A.57) is missing the kT prefactor that must appear when the derivative of a dimensionless log-supply function is used. The sign convention and the kT factor must be corrected, and the derivation of Eq. (15) must be redone, before the claimed correction of Stratton's valence-band equation can be accepted.","section":"Section II.B, Eq. (11); Appendix Eqs. (A.25), (A.55)"},{"comment":"The valence-band change of variables is internally inconsistent. Eq. (A.38) does not follow from Eqs. (A.33)-(A.36): substituting Ebar = Ev - epsilon and epsilon = EV - E into DV(Ebar + Ebar_r) gives a different argument than the printed DV(EV - E + pr^2/2m). Eq. (A.45) also has the wrong sign on the pz^2 term: from E = EV - (pz^2 + pr^2)/(2m*) one obtains E = EV - pz^2/(2m*) - (m/m*)Er, not E = EV - (m/m*)Er + pz^2/(2m*). Although the resulting integration limit in Eq. (A.46) may survive this particular sign error, the derivation as a whole is not reproducible. These errors are load-bearing because Eq. (15) is obtained from this chain of substitutions, and the claimed boundary-term difference between Eq. (A.60) and Stratton's Eq. (A.61) cannot be verified until the change of variables, the Jacobian orientation, and the integration-by-parts boundary term are all corrected.","section":"Appendix A.2, Eqs. (A.38) and (A.45)"},{"comment":"The field-dependent predictions are computed from imposed band-bending scenarios, not from a self-consistent solution. Section II.A explicitly defers the Poisson-continuity-heat self-consistent loops to a forthcoming publication, and Section IV states that all results use the zero-current approximation with pre-set band positions. Consequently, the saturation in Fig. 6, the conduction/valence crossover, and the temperature dependence in Fig. 8 are predictive only under an assumed, unvalidated mapping between the applied field and the band displacement (EC-EF). A different, physically correct mapping could change the crossover point and the saturation behavior even if the emission integrals were exact. The authors acknowledge this limitation in the conclusions, but the manuscript's abstract and Section IV present these results as explanations of experimental observations. The claims should be reframed as illustrative calculations for prescribed band geometries, or the model should be coupled to at least a simplified self-consistent band-bending calculation and the mapping validated.","section":"Section IV, Figs. 6 and 8"}],"minor_comments":[{"comment":"The abstract contains the typo 'Key works' instead of 'Key words', and in Section IV.1 the sentence 'it has such a low magnitude in comparison with that of the emission from the valence band' should read 'conduction band'.","section":"Abstract and Section IV.1"},{"comment":"The symbol Q in the barrier expression is never defined; if it denotes the image-potential coefficient, it should be written explicitly (e.g., Q = q^2/(16*pi*epsilon_0)) so that the barrier model is self-contained.","section":"Eq. (5)"},{"comment":"The expression U(z) = F R(gamma-1)z + z^2/(gamma z + R(gamma-1)) is ambiguous as typeset; parentheses should be added to make the denominator unambiguous, and the parameters should be related to Eq. (5).","section":"Eq. (20)"},{"comment":"A figure caption is duplicated: the caption 'Figure 2: Evolution of the Nottingham heat' refers to a later figure, not the schematic in Section II.A; renumber the figures.","section":"Figure numbering"},{"comment":"The effective masses are quoted as 'm*e = 0.98m kg and m*h = 0.59m kg'; this should be written as m*e = 0.98 m_0 and m*h = 0.59 m_0, where m_0 is the free-electron mass, to avoid dimensional confusion.","section":"Section IV.1"},{"comment":"The text refers to Eqs. (21-22) and Eq. (23) for the barrier approximation and asymptotic forms, but these equations are not all displayed or clearly numbered in the manuscript; please check the equation numbering and display all referenced expressions.","section":"Section III"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (11) is elementary but central: it changes the predicted sign and magnitude of the valence-band current at cryogenic temperatures and invalidates the headline correction to Stratton's formula as written. I am not recommending rejection because the error appears fixable, but the authors must correct the supply function, re-derive the integration by parts, and regenerate all temperature- and field-dependent figures before the paper can be considered for publication. The deferred self-consistent band-bending calculation is also a serious limitation for the experimental claims and should be addressed in framing if not in implementation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nYou should know two things about this paper before spending time on it. First, the computational upgrade to GETELEC is real: the polynomial pre-tabulation of the Gamow function reproduces full JWKB currents and Nottingham heats to RMSE ~0.02 with a sixfold speedup, and the validation against GETELEC-1.0 over 8000 random parameter sets is a reasonable check. Second, the theoretical headline — a correction to Stratton's valence-band emission integral — does not survive contact with the paper's own definitions. Eq. (11) defines lFD(En) = -ln(1/(1+e^{(En-Ef)/kT})) = ln(1+e^{(En-Ef)/kT}). Its derivative is (1/kT)(1-fFD), not -fFD. The appendix replaces fFD by -lFD' in the integration by parts (Eqs. A.25 and A.55). With the stated lFD, that substitution gives a negative supply function, and at cryogenic temperatures the valence-band integral (Eq. 15) is exponentially suppressed, contradicting the paper's own conclusion that valence emission dominates at low T. The difference between Eq. (A.60) and Stratton's Eq. (A.61) is exactly the boundary term of this integration by parts, so the claimed correction is not established. Either the sign in Eq. (11) should be reversed (Ef - En) or the derivation needs a different supply function; as written, it is internally inconsistent.\n\nThe rest of the conduction-band derivation is standard and follows Stratton. The paper honestly states its main limitations: no self-consistent band bending, replacement energy set to EF, smooth surface. But the results section leans on imposed band-bending scenarios without validating the field-to-bending mapping, and the abstract's claim of reproducing experimental I-V behavior is not backed by any experimental comparison in the text. The code is referenced but not accessible in the manuscript.\n\nAll of this is fixable. The polynomial interpolation approach is a genuine practical contribution, and a corrected derivation could make the Stratton claim interesting. But the central algebra must be fixed, verified step-by-step (or machine-checked), and ideally accompanied by one experimental dataset before the paper can be believed.\n\nMy recommendation: send it to peer review, but flag the sign issue prominently. A good referee can sort out whether the correction is real or a typo. Don't cite it in its current form.","headline":"The paper has a useful computational framework, but the headline Stratton correction rests on a Fermi-Dirac sign error that, as written, kills the valence-band current at low temperature.","tokens_in":13138,"tokens_out":2624,"would_cite":false,"duration_ms":23068,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["79.70.+q"],"model":"deepseek-v4-flash","headline":"Correcting Stratton's 1964 valence-band emission equation removes a spurious logarithmic term that changes predicted currents at cryogenic temperatures.","keywords":["electron emission","semiconductors","thermal-field emission","valence band","Fowler-Nordheim","JWKB","Nottingham heat","GETELEC"],"falsifier":"Numerically integrate the full momentum-space integral in Eq. (14) (without the analytic integration-by-parts step) for a representative flat-band semiconductor barrier at cryogenic temperature, and compare it against Eq. (15) and Stratton's Eq. (A.61); whichever expression matches the direct numerical integral is the correct field-emission equation.","tokens_in":12062,"feed_emoji":"⚡","tokens_out":15125,"duration_ms":120836,"temperature":0.7,"pith_summary":"This paper corrects the equation for electron field emission from the valence band of a semiconductor, originally derived by Stratton in 1964. The correction removes a logarithmic term that appears in the original expression after an integration by parts; the paper argues this term is spurious and does not follow from the full momentum-space integral. If the correction is right, it matters most at cryogenic temperatures, where the valence band is the only populated source of electrons and the erroneous term would otherwise distort predicted currents. The authors implement the corrected equations in GETELEC-2.0, a computational tool that uses pre-tabulated polynomial fits to the tunneling barrier to reproduce full JWKB results about six times faster, and use it to reproduce non-linear I-V curves and explain the apparent absence of valence-band electrons in experiments.","feed_headline":"Valence-band emission equation corrected after 60 years","feed_subtitle":"At cryogenic temperatures the fix changes predicted currents and explains missing valence-band electron spectra.","key_machinery":"The load-bearing object is the corrected valence-band emission integral (Eq. 15), obtained by applying integration by parts to the full momentum-space integral and discarding a boundary term that Stratton's 1964 derivation retained. The computational speedup comes from pre-tabulating the Gamow function $G(E_n; F, R, \\gamma)$ as a fourth-order polynomial in $E_n$ whose coefficients are stored over a grid of the barrier parameters $F$ (local field), $R$ (curvature radius), and $\\gamma$ (far-field ratio), so that the JWKB transmission coefficient $D(E_n)=1/(1+e^{G(E_n)})$ can be evaluated in closed form at runtime. This makes the full emission integrals, including the corrected valence-band term, cheap enough to use inside 3D device simulations.","core_discovery":"The central result is a corrected expression for the valence-band contribution to thermal-field emission, Eq. (15): $j_V = L \\int_{-\\infty}^{E_V} l_{FD}(E)\\left(D(E) - \\bar{\\alpha} D(\\bar{\\alpha} E - \\alpha E_V)\\right) dE$, where $\\alpha = m^*/m$, $\\bar{\\alpha} = 1 + m^*/m$, $L = qm/(2\\pi^2\\hbar^3)$, and $D$ is the JWKB-Kemble transmission coefficient. This differs from Stratton's Eq. (A1) by the absence of an additional term, $\\ln\\left\\{1 + e^{E_V/kT}\\right\\}$ times an integral of $D$ over the valence-band edge, which the authors show arises from an incorrect partial integration in the original derivation. They also report that a fourth-order polynomial fit to the Gamow function $G(E_n; F, R, \\gamma)$ reproduces the full JWKB current density and Nottingham heat with RMSE $0.01531$ and $0.02987$ respectively, at roughly sixfold speedup. With this tool, they find that the emitted electron energy distribution transitions from metal-like to conduction-band-dominated to valence-band-dominated as the field decreases, and that a measurable crossover between conduction and valence contributions could reveal the semiconductor band gap.","pith_inferences":["If the removal of the logarithmic term is correct, published analyses of cryogenic semiconductor emission that used Stratton's original expression may need to be revisited, and the discrepancy might show up most strongly in Nottingham heating at low temperature.","The paper's imposed band-bending scenarios are a necessary simplification; the natural next step is to couple GETELEC-2.0 with self-consistent Poisson-continuity-heat loops and check whether the predicted saturation and crossover points move, since the field-to-bending mapping is currently assumed.","The polynomial-fitting strategy for $G(E_n)$ could be exported to any emission model where the barrier is a smooth function of a few parameters, making full JWKB accuracy affordable in large 3D emitter simulations.","The dynamic-range explanation for missing valence-band spectra is testable: run a high-dynamic-range electron energy analyzer at the predicted conduction/valence crossover field and look for the two-peak plus band-gap structure."],"forward_implications":["At cryogenic temperatures, the corrected equation changes the predicted magnitude of valence-band emission, which may explain why valence-band electrons have not been observed in experiments and points to the conditions under which they should be.","The conduction and valence contributions cross at a specific applied field; at that crossover, the model predicts the band gap should become visible in the emitted electron energy distribution.","The model predicts that n-type semiconductor emitters also exhibit current saturation at high fields once the surface degenerates, extending a saturation effect that has been observed experimentally in p-type materials.","Because the polynomial-fit approach reproduces full JWKB results across all emission regimes, the same speedup strategy is applicable to photo-emission and to emitters with more complex supply functions."],"supporting_citations":[{"why":"Stratton's 1955 paper that first formulated field emission from semiconductors, providing the starting equations for the present work.","marker":"[9]"},{"why":"Stratton's 1962 theory of field emission from semiconductors, from which the conduction-band emission equation (Eq. 8) is taken.","marker":"[10]"},{"why":"Stratton's 1964 energy-distributions paper whose valence-band equation (Eq. A1) is corrected here by removing a spurious logarithmic term.","marker":"[15]"},{"why":"Kemble's JWKB method that gives the transmission coefficient used in all emission integrals.","marker":"[16]"},{"why":"Kyritsakis and Djurabekova's GETELEC-1.0 framework and the empirical barrier expression (Eq. 20) that the new version extends.","marker":"[7]"},{"why":"Kyritsakis and Xanthakis' generalized Fowler-Nordheim equation for nanoscopic emitters, providing the asymptotic barrier forms used in numerical JWKB integration.","marker":"[17]"},{"why":"Kyritsakis' general form of the tunneling barrier for sharp emitters, which justifies parametrizing the barrier by F, R, and γ.","marker":"[18]"},{"why":"Baskin, Lvov, and Fursey's theory of field emission from semiconductors, the source of the zero-current approximation adopted in the simulations.","marker":"[11]"},{"why":"Choueib et al.'s experimental current saturation in H-passivated Si nanowires, which the model reproduces and explains.","marker":"[1]"}],"fun_headline_variants":["Valence-band term fixed in 60-year-old emission equation","Corrected theory explains missing valence electrons","Faster thermal-field emission simulations for semiconductors","Semiconductor emission equation corrected after decades"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions stand on the assumption that the three imposed band-bending scenarios (flat band, $E_C = E_F$, $E_C \\ll E_F$) correctly represent how a real semiconductor's bands respond to the applied field, and that mapping is never derived or validated.","fun_headline_variants_meta":{"raw":{"variants":["Valence-band term fixed in 60-year-old emission equation","Corrected theory explains missing valence electrons","Faster thermal-field emission simulations for semiconductors","Semiconductor emission equation corrected after decades"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1505,"prompt_tokens":958,"completion_tokens":547,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":490}},"tokens_in":574,"tokens_out":547,"duration_ms":6433,"temperature":1.0,"reasoning_tokens":490,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T06:00:01.501317+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the full momentum-space integral in Eq. (14) (without the analytic integration-by-parts step) for a representative flat-band semiconductor barrier at cryogenic temperature, and compare it against Eq. (15) and Stratton's Eq. (A.61); whichever expression matches the direct numerical integral is the correct field-emission equation.","supporting_citations":[{"cited_title":"Field emission from semiconductors,","cited_arxiv_id":null,"evidence_quote":"Stratton's 1955 paper that first formulated field emission from semiconductors, providing the starting equations for the present work."},{"cited_title":"Energy distributions of field emitted electrons,","cited_arxiv_id":null,"evidence_quote":"Stratton's 1964 energy-distributions paper whose valence-band equation (Eq. A1) is corrected here by removing a spurious logarithmic term."},{"cited_title":"A Contribution to the Theory of the B. W. K. Method,","cited_arxiv_id":null,"evidence_quote":"Kemble's JWKB method that gives the transmission coefficient used in all emission integrals."},{"cited_title":"A general computational method for electron emission and thermal effects in field emitting nanotips,","cited_arxiv_id":null,"evidence_quote":"Kyritsakis and Djurabekova's GETELEC-1.0 framework and the empirical barrier expression (Eq. 20) that the new version extends."},{"cited_title":"Derivation of a generalized Fowler-Nordheim equation for nanoscopic field -emitters,","cited_arxiv_id":null,"evidence_quote":"Kyritsakis and Xanthakis' generalized Fowler-Nordheim equation for nanoscopic emitters, providing the asymptotic barrier forms used in numerical JWKB integration."},{"cited_title":"General form of the tunneling barrier for nanometrically sharp electron emitters,","cited_arxiv_id":null,"evidence_quote":"Kyritsakis' general form of the tunneling barrier for sharp emitters, which justifies parametrizing the barrier by F, R, and γ."}],"review_version":1}