{"id":"71c66fad-36a0-4991-b5b7-a3d7cdb16d0b","arxiv_id":"2505.24150","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A fully analytic semiclassical derivation gives the Zener pair-production rate and tunnelling current in biased bilayer graphene, including absolute normalization.","lead":"This paper derives a fully analytic formula for the rate of electron-hole pair production (Zener tunnelling) in biased bilayer graphene under an electric field, using a semiclassical analytic-continuation method. The value for readers is a closed-form expression for the tunnelling current with an absolute normalization that earlier hybrid numerical treatments could not fix.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The absolute normalization in Eq. (4.8) rests on the un-derived finite-ky replacement rule (3.9); without an independent check of T(ky), the claimed prefactor is not established.","rationale":"The reader's weakest_assumption identifies exactly the finite-ky replacement rule (3.9), and my reading of the derivation confirms that this rule is the least secure link in the chain leading to Eq. (4.8). The ky=0 matching calculation is internally coherent: the analytic continuation phases, the matching coefficients, and the solution of Eqs. (3.31) all reproduce Eq. (3.32) as written. The problematic step is the assertion that finite ky is recovered by the replacement rule, which is stated without a derivation. The final rate (4.8) depends on the Gaussian width β1 and on the oscillation phase in the bracket, both of which come from that rule. An independent numerical solution of the original ODE at finite ky would settle whether Eq. (3.33) is correct; this is a concrete and feasible test. The paper otherwise has genuine strengths: the Zwaan analytic-continuation framework is applied carefully, the matching conditions are explicit, and the final expression is fully analytic and falsifiable. I therefore do not see a reason to move the verdict away from CONDITIONAL; the reader's requested condition is precisely the check needed. No ad hominem or theatrical framing is intended: the issue is a missing derivation for a load-bearing technical step, not a claim about the authors' process.","tokens_in":31709,"tokens_out":31219,"duration_ms":309082,"concrete_test":"Solve the two-component bilayer Schrödinger equation (3.1) at ε=0 by shooting from x≫x+ with outgoing/decaying boundary conditions, for a fixed S0≈10 and several ky values around the Gaussian width implied by Eq. (4.8). Extract |T(ky)|² and the transmission phase, and compare with Eq. (3.33). If log|T(ky)|² deviates from -√2(S0+Sy) at the Gaussian width by more than ~10%, or if the transmission zero is not at (S0-Sy)=π√2, then rule (3.9) is invalid and the prefactor in Eq. (4.8) must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central new result is the absolute pair-production rate in Eq. (4.8). The prefactor of that rate is controlled entirely by the Gaussian integral over transverse momentum ky, which uses T(ky) from Eq. (3.33). That T(ky) is obtained from the ky=0 matching solution by the replacement rule (3.9), e^{±iπ/4 S0} → e^{±iπ/4(S0 ∓ i Sy)}, stated in Sec. III.A.1 without derivation and applied in Sec. III.C. This rule does two things: it puts the factor e^{-Sy/√2} into the amplitude, fixing the Gaussian width β1, and it puts the phase Sy/√2 into the sine, fixing the interference bracket in Eq. (4.8). A different finite-ky extension—for instance replacing S0 by the full real action S0+Sy in the matching exponentials—would change the interference phase and could change the prefactor; an error in the coefficient of Sy would directly rescale the absolute rate. The paper offers no independent derivation of Eq. (3.9) and no numerical check of Eq. (3.33) against direct solution of the bilayer Schrödinger equation or against the numerical results of Ref. [13]. Since the claimed absolute normalization is the paper's main advance, this unsecured step is the most load-bearing concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a fully analytic calculation of the Zener electron-hole pair-production rate in biased bilayer graphene, using the Zwaan analytic-continuation variant of semiclassical theory. Starting from the two-band Hamiltonian (3.1), the authors construct semiclassical wavefunctions in the three spatial regions, match them by analytic continuation around the turning points, and obtain the ky-dependent transmission amplitude T(ky) in Eq. (3.33). Integrating |T(ky)|^2 over transverse momentum gives the central result, Eq. (4.8), for the pair-production rate per unit area, including an oscillatory interference term. The paper emphasizes that this is an absolute normalization, in contrast to the earlier hybrid numerical-WKB treatment of Ref. [13], and it also discusses the device-length dependence of the wavefunction normalization.","tokens_in":32034,"tokens_out":15772,"duration_ms":139287,"significance":"If the result is correct, Eq. (4.8) fixes the absolute prefactor of the Zener tunnelling rate in biased bilayer graphene, including the oscillatory field dependence coming from quantum interference. This is a nontrivial advance over prior work, which relied on fitting WKB coefficients to numerical solutions. The method itself, analytic continuation in the complex plane, is an elegant and parameter-free technique that avoids turning-point matching to special functions; the constants β0 and β1 are computed from definite integrals, with no fitted parameters. The paper also clarifies the non-standard structure of the semiclassical wavefunctions, in particular the necessity of retaining decaying components in classically allowed regions. However, the central prefactor rests on a finite-ky replacement rule, Eq. (3.9), that is stated without derivation, so the absolute normalization claim is not yet established to the standard claimed.","major_comments":[{"comment":"The finite-ky replacement rule e^{±iπ/4 S0} → e^{±iπ/4(S0 ∓ i Sy)} is stated without derivation. This rule is used in §III.C to obtain T(ky), and it controls both the Gaussian width β1 in the ky integral and the interference phase S0−Sy in Eq. (4.8). Because the absolute normalization of the pair-production rate is the paper's central claim, this step needs a derivation from the ky-dependent semiclassical momentum p_F(x)^2−ky² and the associated analytic continuation around the shifted turning points, or an independent numerical check. As written, the rule is an ad hoc ansatz; a different natural extension, such as replacing S0 by S0−Sy in the matching exponentials, would change the prefactor and the interference phase.","section":"§III.A.1, Eq. (3.9)"},{"comment":"The transmission amplitude and the resulting pair-production rate are not benchmarked against any numerical solution of the bilayer Schrödinger equation (3.1), nor against the numerical WKB results of Ref. [13]. The numerical check in Appendix A is for the Dirac Hamiltonian, not for biased bilayer graphene. Since the paper claims that leading-order WKB provides the absolute prefactor and not just the exponent, an independent check is needed to rule out O(1) corrections to the prefactor in Eq. (4.8).","section":"§III.C, Eq. (3.33) and §IV.2, Eq. (4.8)"},{"comment":"Equation (3.4) as printed, p(x)^2 = ±2m√(F²x²−Δ²−k_y²), is dimensionally inconsistent and cannot be reconciled with the definitions of p_F, p_A, S0, and Sy in Eqs. (3.6)–(3.8). The correct relation should be p(x)^2 = ±2m√(F²x²−Δ²) − k_y², with p denoting the x-component of momentum; this yields the turning points in Eq. (3.5) and the expansion used for Sy. Please correct this equation and ensure all subsequent ky expansions are consistent with it.","section":"Eq. (3.4)"}],"minor_comments":[{"comment":"The displayed text contains repeated copies of the abstract and of Section II, and an earlier draft with unresolved 'Levitov [cite]' placeholders; the published version should be cleaned of these artifacts.","section":"General presentation"},{"comment":"There is a typographical error in 'the distance from the turning point is march larger than the wavelength'; this should read 'much larger'.","section":"§III.A.2"},{"comment":"The phrase 'multiply the wave fanction by √p−∞' should read 'multiply the wave function by √p−∞'.","section":"§IV.3"},{"comment":"The notation gsτ = 4 for the spin-valley degeneracy is not defined; please clarify, e.g. as g_{sv} = 4.","section":"§IV.2, Eq. (4.8)"},{"comment":"In the clockwise continuation line, 'cl : (−ipF(x))^{−1/2} → −i(pA(x))^{1/2}' appears to have the wrong power; it should likely be (pA(x))^{−1/2} to be dimensionally consistent.","section":"Appendix A.3, Eq. (A13)"},{"comment":"The axis label 'J/v exp(πΔ²_k/vF) vs −πΔ²_k/vF' is garbled; please rewrite it so that the exponential rescaling and the horizontal variable are immediately clear.","section":"Fig. 6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to be a poorly cleaned LaTeX file with duplicated text and placeholder citations. The main scientific concern is the un-derived finite-ky replacement rule (3.9), which controls the absolute prefactor claimed as the paper's key advance. If the authors can derive that rule from the ky-dependent semiclassical action and add a numerical validation of T(ky) against the bilayer Schrödinger equation or the earlier results of Ref. [13], the paper could be a useful contribution. Without that, the absolute normalization claim is not established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a real step forward for Zener tunneling in biased bilayer graphene—the first fully analytic, parameter-free expression for the absolute pair-production rate—but the main new number rests on an unproven finite-ky rule. I would send it to review, with a clear request to fix or verify that rule.\n\nWhat's genuinely new: applying Zwaan analytic continuation to the bilayer problem avoids the hybrid numerical-WKB fitting of Nandkishore and Levitov. The matching at ky=0 is worked out in detail, and the observation that decaying components must be kept in allowed regions is correct and important. The normalization section distinguishes long- and short-device regimes and gets the same rate for both, which is a good consistency check. Nothing in the derivation is fitted; beta0 and beta1 are computed integrals. That is the right way to present such a calculation.\n\nThe soft spots. First and most load-bearing: Eq. (3.9), the rule e^{±iπ/4 S0} -> e^{±iπ/4(S0 ∓ i Sy)}, is simply stated as 'required' with no derivation. It controls both the Gaussian width beta1 and the interference phase in the sine, so it fixes the absolute prefactor in Eq. (4.8). A different plausible extension—replacing S0 by the full real action S0+Sy in the matching exponentials—would shift the interference bracket and the prefactor. The paper gives no independent check of T(ky) against direct numerical integration of the bilayer Schrödinger equation, nor against the numerical coefficients in Ref. [13]. The appendix has a numerical check only for the linear Dirac case. This is exactly the gap the reader flagged, and I agree it prevents 'accept' as is.\n\nSecond, the introduction promises a small logarithmic correction to the gap dependence, and the promise is repeated in the text, but I cannot find it in the final rate (4.8), which is a pure exponential in S0 ~ Δ^{3/2}/F. Either the log correction is hiding in beta0 or it was dropped; as written, the paper is internally inconsistent on this point. That should be fixed.\n\nAlso the arXiv text contains repeated/duplicated passages from an earlier draft; it needs a careful edit, but that is cosmetic.\n\nWho is this for: people working on TFETs or field-induced tunneling in 2D materials, and anyone applying the Zwaan method in condensed matter. It deserves a serious referee. My recommendation: major revision, with the finite-ky rule derived (or at least numerically verified against the bilayer ODE) and the log-correction discrepancy resolved.","headline":"A genuinely new fully analytic Zener rate for biased bilayer graphene, but the absolute prefactor rests on an unproven finite-ky replacement rule; worth refereeing after major revision.","tokens_in":32487,"tokens_out":2819,"would_cite":true,"duration_ms":29271,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By analytic continuation of semiclassical wavefunctions around the turning points, this paper derives the absolute pair-production rate for Zener tunnelling in biased bilayer graphene and fixes the normalization of the tunnelling current.","keywords":["Zener tunnelling","biased bilayer graphene","semiclassical analytic continuation","pair production rate","absolute normalisation","tunnelling current","band gap","evanescent wavefunction components"],"falsifier":"Solve the two-component Schrödinger equation (3.1) at fixed $k_y$ with absorbing boundary conditions, extract $|T(k_y)|^2$ from the asymptotic running wave, integrate over $k_y$, and compare the resulting pair-production rate with Eq. (4.8) across a range of $F$ and $\\Delta$; a discrepancy in the prefactor, or in the amplitude or phase of the $1-\\cos(\\sqrt{2}S_0-\\pi/8)/2^{1/4}$ factor, would falsify the absolute normalization.","tokens_in":31507,"feed_emoji":"⚡","tokens_out":7921,"duration_ms":72443,"temperature":0.7,"pith_summary":"The paper claims that Zener tunnelling in biased bilayer graphene—a two-layer carbon sheet with a voltage-tunable band gap—can be computed fully analytically by continuing semiclassical wavefunctions around the classical turning points in the complex plane. The central output is a closed-form pair-production rate per unit area, $\\frac{dn}{dt}=\\frac{m\\Delta^2 e^{-\\sqrt{2}S_0}}{2\\pi(\\pi\\beta_1)^{1/2}}\\left(\\frac{F^2}{m\\Delta^3}\\right)^{3/4}\\left[1-\\frac{\\cos(\\sqrt{2}S_0-\\pi/8)}{2^{1/4}}\\right]$, which fixes the absolute normalization of the tunnelling current rather than just its exponential suppression. Along the way the paper finds that wavefunctions in classically allowed regions must keep exponentially decaying components, otherwise the matching conditions fail. If correct, the result turns earlier exponential estimates into a definite current formula applicable to devices such as tunnel field-effect transistors.","feed_headline":"Analytic continuation fixes Zener tunnelling rate in bilayer graphene","feed_subtitle":"Fully analytic result yields absolute pair-production rate and oscillatory tunnelling current across the gap.","key_machinery":"The machinery is the analytic continuation of semiclassical wavefunctions around the two turning points $x_\\pm$, performed by rotating $x-x_\\pm=\\rho e^{i\\phi}$ through $\\pm\\pi$ so that matching across the barrier never touches the singular points. The exponents are controlled by the actions $S_0$ and $S_y$ of Eq. (3.8); the phase factors $e^{\\pm i\\pi S_0/4}$ carry the tunnel suppression, and the replacement rule $e^{\\pm i\\pi S_0/4}\\to e^{\\pm i\\pi(S_0\\mp iS_y)/4}$ extends the $k_y=0$ solution to finite transverse momentum, turning the integral over $k_y$ into the prefactor $(F^2/m\\Delta^3)^{3/4}$. Matching at $x_-$ and $x_+$ determines both the scattering phase $\\phi=-\\pi/4$ and the transmission amplitude $T(k_y)$, with unphysical exponentially growing components cancelled exactly.","core_discovery":"Starting from the effective two-component Schrödinger equation for biased bilayer graphene, the paper derives the transmission amplitude $T(k_y)=-2e^{i\\pi/4}e^{-(S_0+S_y)/\\sqrt{2}}\\sin\\left(\\frac{S_0-S_y}{\\sqrt{2}}\\right)$, with $S_0=\\beta_0\\sqrt{2m}\\,\\Delta^{3/2}/F$ and $S_y=\\beta_1\\frac{k_y^2}{F}\\sqrt{\\frac{\\Delta}{2m}}$, where $\\beta_0=1.748$ and $\\beta_1=1.198$. Integrating $|T(k_y)|^2$ over transverse momentum yields Eq. (4.8), an absolutely normalized pair-production rate. The paper's structural claim is that the matching problem requires two independent decaying solutions in the forbidden region and that decaying components must be retained in classically allowed regions; omitting them gives inconsistent matching. The matching conditions at both turning points fix the scattering phase and the amplitudes that earlier treatments left to numerical fitting, thereby establishing the absolute normalization and a small correction to the gap dependence of the tunnelling-rate coefficient.","pith_inferences":["The unproved replacement rule (3.9) is the main thing to test: a direct numerical solution of the two-component Schrödinger equation at finite $k_y$ would confirm or correct the prefactor $(F^2/m\\Delta^3)^{3/4}$ and the oscillation amplitude.","If the absolute normalization is correct, existing exponent-only estimates of Zener currents in bilayer-graphene tunnel transistors can be upgraded to full current-voltage curves including interference oscillations.","The same analytic-continuation strategy should transfer to other gapped systems with quartic or non-quadratic dispersion where Airy-function matching is unavailable; the Dirac limit treated in the appendix is one boundary case.","The evanescent components in allowed regions suggest a connection to boundary-localised modes of non-Hermitian type, an analogy the paper itself raises; transport through p-n junctions could search for the predicted transmission zeros."],"forward_implications":["The pair-production rate per unit area is fixed in absolute units by Eq. (4.8), and after multiplying by the spin-valley degeneracy $g_{s\\tau}=4$ the tunnelling current for a device of length $L$ and width $W$ is $I=LW\\,e\\,dn/dt$.","The tunnelling rate oscillates with gap and field through the factor $1-\\cos(\\sqrt{2}S_0-\\pi/8)/2^{1/4}$, so the current is not a monotonic exponential of $1/F$.","Two linearly independent decaying solutions are required inside the barrier, and exponentially decaying components must be retained in classically allowed regions for the matching to be consistent.","At transverse momenta where $T(k_y)=0$ the coefficient of the decaying component remains nonzero, implying localised modes with zero transmission at those $k_y$.","The same normalized pair-production rate results in both the long-device (linear-dispersion) and short-device (quadratic-dispersion) regimes, so the absolute normalization is independent of device length across the two limits."],"supporting_citations":[{"why":"provides the earlier hybrid WKB-plus-numerics treatment of Zener tunnelling in bilayer graphene whose absolute normalization this paper establishes and refines.","marker":"[13]"},{"why":"supplies the analytic-continuation method for computing transmission coefficients via phase integrals around turning points, which is the core technique used here.","marker":"[18]"},{"why":"gives the standard complex-plane continuation rules used to match semiclassical solutions across turning points.","marker":"[19]"},{"why":"provides the bilayer graphene band model and parameters entering the dispersion, effective mass, and device-regime estimates.","marker":"[20]"},{"why":"supports the claim that common-path interference oscillations in Zener tunnelling are universal for higher-than-quadratic dispersions.","marker":"[21]"},{"why":"underpins the treatment of exponentially decaying components in complex WKB that the paper finds essential in classically allowed regions.","marker":"[22]"}],"fun_headline_variants":["Absolute Zener tunnelling rate from analytic continuation","Analytic continuation nails bilayer graphene tunnelling rate","Tunnelling rate fixed by semiclassical analytic continuation","Exact pair production rate in biased bilayer graphene","Bilayer graphene tunnelling: absolute rate via Zwaan method"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the unproved replacement rule (3.9) that converts the zero-transverse-momentum matching solution into finite transverse momentum, together with the assumption that leading-order WKB wavefunctions give the order-one prefactor and not just the exponential suppression.","fun_headline_variants_meta":{"raw":{"variants":["Absolute Zener tunnelling rate from analytic continuation","Analytic continuation nails bilayer graphene tunnelling rate","Tunnelling rate fixed by semiclassical analytic continuation","Exact pair production rate in biased bilayer graphene","Bilayer graphene tunnelling: absolute rate via Zwaan method"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000218,"raw_usage":{"total_tokens":1410,"prompt_tokens":890,"completion_tokens":520,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":442}},"tokens_in":506,"tokens_out":520,"duration_ms":4942,"temperature":1.0,"reasoning_tokens":442,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:32:55.764648+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the two-component Schrödinger equation (3.1) at fixed $k_y$ with absorbing boundary conditions, extract $|T(k_y)|^2$ from the asymptotic running wave, integrate over $k_y$, and compare the resulting pair-production rate with Eq. (4.8) across a range of $F$ and $\\Delta$; a discrepancy in the prefactor, or in the amplitude or phase of the $1-\\cos(\\sqrt{2}S_0-\\pi/8)/2^{1/4}$ factor, would falsify the absolute normalization.","supporting_citations":[{"cited_title":"Spontaneous momentum polarization and diodicity in Bernal bilayer graphene","cited_arxiv_id":"2302.04261","evidence_quote":"provides the earlier hybrid WKB-plus-numerics treatment of Zener tunnelling in bilayer graphene whose absolute normalization this paper establishes and refines."},{"cited_title":"A theory of the electri- cal breakdown of solid dielectrics,","cited_arxiv_id":null,"evidence_quote":"supplies the analytic-continuation method for computing transmission coefficients via phase integrals around turning points, which is the core technique used here."},{"cited_title":"¨Uber das verhalten eines elektrons im homo- genen elektrischen feld nach der relativistischen theorie diracs,","cited_arxiv_id":null,"evidence_quote":"gives the standard complex-plane continuation rules used to match semiclassical solutions across turning points."},{"cited_title":"Consequences of dirac’s theory of positrons,","cited_arxiv_id":null,"evidence_quote":"provides the bilayer graphene band model and parameters entering the dispersion, effective mass, and device-regime estimates."},{"cited_title":"Common-path in- terference and oscillatory zener tunneling in bilayer graphene p-n junctions,","cited_arxiv_id":null,"evidence_quote":"supports the claim that common-path interference oscillations in Zener tunnelling are universal for higher-than-quadratic dispersions."},{"cited_title":"Ground-state decay rate for the zener breakdown in band and mott insulators,","cited_arxiv_id":null,"evidence_quote":"underpins the treatment of exponentially decaying components in complex WKB that the paper finds essential in classically allowed regions."}],"review_version":1}