{"id":"833ce02a-54ba-49f2-b766-693d0c1ac1ba","arxiv_id":"2502.06265","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper recomputes barrier transmission for quadratic-band-crossing semimetals including previously omitted evanescent waves, yielding corrected conductance and Fano factor predictions.","lead":"Electrons in certain semimetals have an unusual quadratic relation between momentum and energy, and this paper redoes the calculation of how they tunnel through a voltage barrier after finding that earlier work dropped the decaying solutions of the wave equation. With those decaying, or evanescent, waves included, the predicted transmission, electrical conductance, and shot-noise Fano factor change for these materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 2d evanescent-mode treatment survives scrutiny, but the 3d calculation uses spinors in Eq. (24) that are singular at normal incidence and lack specified normalizations, so the 3d transmission results are not well-defined as presented.","rationale":"Good-faith reading: the paper's central correction is to include evanescent waves and impose full derivative matching. In 2d this is internally consistent: the dispersion relation yields four kx values; for imaginary kx the spinor Ψ_- is the eigenvector of H with eigenvalue +E, exactly as used; and integrating the second-order ODE twice gives the stated conditions. The reader's identified weakest assumption (boundary conditions) therefore does not land as stated. However, the 3d part is not independently verifiable from the manuscript. Eq. (24) has coordinate singularities at k⊥=0 and omits normalizations, so the matching calculation is underspecified. That is a genuine load-bearing gap for the 3d results. Since the reader already issued a conditional verdict, the appropriate recommendation is UNCHANGED: the verdict remains conditional, with the condition now directed at the 3d eigenvector basis.","tokens_in":13077,"tokens_out":27973,"duration_ms":253316,"concrete_test":"Compute the normalized Ψ_{+,1} and Ψ_{+,2} from Eq. (24) in the limit k⊥→0 along two different azimuthal directions (e.g., kx→0 with ky=0 versus ky→0 with kx=0). If the spinor limits (or the resulting 16x16 matching matrix determinant) differ, the 3d scattering amplitudes at normal incidence are ill-defined. Also check unitarity R+T=1 for representative parameters in Figs. 7-8.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The boundary-condition worry raised by the reader is probably not the critical flaw: for the Hamiltonian (7), a step potential V(x) yields continuity of both spinor components and their derivatives, and the evanescent modes with ε=-E (using the Ψ_- branch) do satisfy Hψ=Eψ, so the 2d correction is plausible. The load-bearing gap is in Sec. III. The eigenvectors (24) contain explicit 1/(kx - i ky) and 1/(kx + i ky) factors, so they are singular at k⊥=0. Normal incidence (θ=0) is included in the polar plots in Figs. 7-8, and θ=0 lies inside the conductance integral (33). The normalization factors N_{±,s} are not given, and no limiting procedure is specified. As a result, the 16 matching equations (31) are not well-defined at exactly k⊥=0, and the plotted T(θ=0) values cannot be independently reproduced; if the limit is direction-dependent, the 3d transmission amplitudes are not trustworthy. This supports the reader's conditional verdict, but for a different reason.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper re-examines transmission through a rectangular electrostatic barrier in semimetals with quadratic band-crossing points (QBCPs). The author argues that Ref. [11] discarded evanescent solutions of the Schrödinger equation, so the boundary conditions were only partially satisfied; keeping the evanescent modes gives a square system of eight matching equations in two dimensions and sixteen in three dimensions. The transport quantities T(E,V0,phi), conductivity sigma, and Fano factor F are recomputed and shown in polar and line plots, with a closed low-energy expression for the 2d case in Eq. (14). The paper is framed as a correction of previously published results.","tokens_in":13262,"tokens_out":19521,"duration_ms":182087,"significance":"If the 2d calculation is correct, it provides a nontrivial correction to a published transport calculation and yields an explicit falsifiable prediction at small E in Eq. (14), including the absence of Klein tunneling. The 2d formalism is standard and the inclusion of evanescent modes is physically reasonable; I agree with the skeptical assessment that the boundary-condition matching itself is not the central weakness. The 3d section, however, is not in a usable state as written: the outgoing evanescent term in Eq. (30) grows exponentially, the eigenvectors in Eq. (24) are singular at normal incidence with unspecified normalizations, and the conductance formula relies on an unproved statement that interband transmission vanishes. The paper also does not provide code, data, or full transmission amplitudes, which would be needed to replace previously published curves. No circularity concern attaches to the central calculation: the transmittance follows from solving a boundary-value problem with no tunable parameters.","major_comments":[{"comment":"The transmitted evanescent term in the right-side wavefunction is written with the factor e^{kappa(z-L)}, which grows without bound for z>L because kappa = sqrt(2mE + k_perp^2) > 0. This makes the 3d scattering state non-normalizable and does not correspond to an outgoing evanescent mode; the analogous 2d term in Eq. (11) correctly uses e^{-kappa(x-L)}. Since this term enters the matching equations Eq. (31), the 3d results in Figs. 7-9 are not well-defined as written.","section":"III.A, Eq. (30)"},{"comment":"The four 3d eigenvectors contain explicit factors 1/(k_x - i k_y), 1/(k_x + i k_y), and their squares, so they are singular at normal incidence k_perp = 0; the normalization factors N_{+-s} are not given. Since theta = 0 is included in the polar plots and is an endpoint of the conductance integral Eq. (33), the transmission amplitudes are not well-defined at normal incidence. The author should specify the orthonormalization convention, the branch of the square root for imaginary k_z, and a k_perp -> 0 limiting procedure, and show that the limit is independent of the azimuthal direction in the transverse plane.","section":"III.A, Eq. (24)"},{"comment":"The assertion that |t_{n,2}| = 0 for an incident Psi_{+,1} state is stated without derivation. The boundary-condition system Eq. (31) couples all four spinor components, and no symmetry is displayed that would obviously decouple the two positive-energy bands at fixed k_perp; a scalar barrier can in general scatter between degenerate bands. Because Eq. (33) uses only |t_{n,1}|^2 for the conductance, this assumption is load-bearing. Please provide the symmetry argument or a numerical check that includes |t_{n,2}|^2 before the 3d conductance can be accepted.","section":"III.B, footnote 2 and Eq. (33)"},{"comment":"The paper's central quantitative claim is that the previous transmission coefficients from Ref. [11] are incorrect, yet the new amplitudes are not shown ('extremely long') and no code or data are supplied. The low-energy formula Eq. (14) provides one check, but it does not validate the intermediate-energy curves in Figs. 2-9. I request the full amplitudes, or at least machine-readable data or scripts that solve Eqs. (12) and (31), so the corrected curves can be independently verified.","section":"II.A and III.A, Eqs. (11)-(13) and (30)-(32)"}],"minor_comments":[{"comment":"The last boundary condition should read d_x phi_M(L) = d_x phi_R(L), not d_x phi_L(x)|_{x=L}; Eq. (31) has the analogous typo (and an x versus z slip in the derivative argument).","section":"II.A, Eq. (12)"},{"comment":"The denominator in the Fano factor should be sum_n T_n, not sum_{n~} T_{n~}; the tilde over n appears to be a typographical artifact.","section":"I, Eq. (2)"},{"comment":"The lower-panel captions of Figs. 3, 5, and 8 label the reflection coefficient as T(E,V0,phi) or T(E,V0,theta); these should be R(E,V0,phi) and R(E,V0,theta), respectively.","section":"Fig. 3, Fig. 5, Fig. 8"},{"comment":"The section headings 'F ormalism' and 'T ransmission coefficients' contain TeX-spacing artifacts and should be corrected.","section":"II and III headings"},{"comment":"The special case E = V0 is not discussed; inside the barrier the four k_x solutions change their character and the independent modes can coalesce, so a brief remark about this non-generic point would help.","section":"II.A and III.A, barrier-region modes"}],"recommendation":"major_revision","confidential_remarks":"The paper is heavily self-referential, but I do not see circularity in the central boundary-value calculation; my concerns are technical. The 3d section needs substantive revision (growing evanescent factor, singular/un-normalized spinors, and the unproved decoupling of the two bands) before the claims can be evaluated. The 2d part is more credible and may survive once the reproducibility gap is addressed. I would not reject at this stage because the issues appear fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plainly: this is a self-correction and the central 2d fix is real. The earlier Annals paper on QBCP barriers missed the evanescent branches of the quadratic dispersion; this version includes them, redoes the matching, and recomputes transmission, conductance, and Fano factor. The method is standard plane-wave matching with no free parameters, and the small-E formula Eq. (14) is a concrete new output. The paper is honest about the earlier miss, and the citations for analogous evanescent-mode treatments (bilayer graphene, multi-Weyl semimetals) are appropriate.\n\nThe 2d part holds up under scrutiny. The boundary conditions—continuity of both spinor components and their x-derivatives—are the natural ones for the Hamiltonian in Eq. (7), and using the Ψ− branch with ε = −E gives legitimate evanescent solutions. The reader's worry about a missing interface condition does not land.\n\nThe soft spot is the 3d section. The eigenvectors in Eq. (24) contain 1/(kx − i ky) and 1/(kx + i ky) factors, and the normalizations N±,s are never given. At k⊥ = 0 (normal incidence), the spinors are singular. Normal incidence is plotted in Figs. 7–8 and lies inside the k⊥ integral for conductance, Eq. (33). Without a specified limiting procedure or normalized basis, the 16 matching equations are not well-defined at exactly θ = 0, and the 3d transmission curves cannot be independently reproduced. That is a load-bearing gap, not a cosmetic one.\n\nMinor issues: the full transmission amplitudes are withheld as “extremely long,” there is no code or data, and there are typos in the boundary-condition subscripts and figure captions (e.g., ∂xφL in the second boundary condition of Eq. (12)).\n\nVerdict: worth a serious referee, with a request to fix the 3d basis or restrict the paper to 2d. People working on QBCP mesoscopics should use the 2d corrected coefficients.","headline":"Self-correction that fixes a real 2d QBCP tunneling error; the 3d section is undone by singular, unnormalized spinor basis.","tokens_in":13801,"tokens_out":2509,"would_cite":true,"duration_ms":20541,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.23.-b","73.40.Gk","72.10.-d"],"model":"deepseek-v4-flash","headline":"For quadratic-band-crossing semimetals, the previous tunneling calculation missed evanescent waves and is now corrected.","keywords":["quadratic band crossing","evanescent waves","rectangular potential barrier","transmission coefficient","Fano factor","conductance","Luttinger semimetal","Klein tunneling"],"falsifier":"Solve the same rectangular-barrier scattering problem on a lattice QBCP model by direct numerical integration of the time-independent Schrödinger equation and compare the transmission at normal incidence with the small-$E$ formula $T=4E\\,\\sec^2\\phi\\,\\csc^2(L\\sqrt{2mV_0})/V_0$; a mismatch would show that the assumed derivative-matching conditions are incomplete or incorrect.","tokens_in":12842,"feed_emoji":"⚛️","tokens_out":8223,"duration_ms":67287,"temperature":0.7,"pith_summary":"This paper corrects the theory of tunneling through a rectangular potential barrier in semimetals whose bands touch at a quadratic band-crossing point (QBCP), in both two and three dimensions. The earlier treatment of this problem left out the evanescent-wave solutions, so the wavefunction could not satisfy all the boundary conditions at the barrier edges. With those solutions included, the paper recomputes the transmission and reflection coefficients, the ballistic conductance, and the Fano factor for representative parameters. If the correction is right, the published curves for these quantities are outdated and the corrected results, including a new small-energy formula, are the valid ones. The conclusion that QBCPs do not show Klein tunneling survives the correction.","feed_headline":"Missing evanescent waves recalculate semimetal barrier tunneling","feed_subtitle":"Corrected transmission, conductance, and Fano-factor curves for 2D and 3D quadratic-band-crossing materials.","key_machinery":"The load-bearing object is the piecewise scattering wavefunction that includes evanescent modes with imaginary wavevectors, $k_x=\\pm i\\sqrt{2mE+q_n^2}$ outside the barrier and similarly inside, in addition to the propagating modes $k_x=\\pm\\sqrt{2mE-q_n^2}$. The paper defines a QBCP as a nodal point where the low-energy bands disperse quadratically in momentum, and the Hamiltonian is quadratic in derivatives. The matching procedure integrates the Schrödinger equation across the barrier edges twice, yielding continuity of the two spinor components and their $x$-derivatives at $x=0$ and $x=L$; in 3d the same conditions are imposed on the four spinor components at $z=0$ and $z=L$. This gives exactly eight equations for the eight coefficients in 2d and sixteen in 3d, closing the system that the earlier work left open.","core_discovery":"The central claim is that a correct scattering solution for a rectangular potential in a QBCP semimetal must include evanescent plane-wave parts alongside the propagating ones, for the incoming, reflected, and transmitted waves and inside the barrier. For the two-dimensional model with Hamiltonian $H^{\\mathrm{kin}}_{2d}(k_x,k_y)=\\frac{1}{2m}[2k_x k_y\\,\\sigma_x+(k_y^2-k_x^2)\\sigma_z]$, the dispersion $\\varepsilon_{2d}(k_x,k_y)=(k_x^2+k_y^2)/(2m)$ gives four wavevector solutions for fixed transverse momentum, two real and two imaginary; dropping the imaginary pair leaves the derivative-matching equations underdetermined. The paper keeps all eight amplitude coefficients in 2d (sixteen in 3d), fixes them by continuity of both spinor components and their $x$- (or $z$-) derivatives at the two barrier edges, and obtains the transmission amplitude $t_n$. It then computes $T=|t_n|^2$, the conductance, and the Fano factor via the Landauer formula, and derives the small-$E$ limit $T(E,V_0,\\phi)=4E\\,\\sec^2\\phi\\,\\csc^2(L\\sqrt{2mV_0})/V_0+O(E^2)$ in 2d. The corrected results are shown in polar and conductivity plots, and compared with normal electron gases and graphene.","pith_inferences":["Editorial inference: the same evanescent-wave omission could affect other QBCP transport calculations, such as tunneling through delta-function potentials or Josephson-junction Andreev spectra; the paper explicitly lists those problems as next steps, so a consistent treatment will need the full 8- or 16-dimensional matching.","Editorial inference: the low-energy scaling $\\sigma\\propto E$ and $F\\to 1$ in 2d is sharp enough to test experimentally in a ballistic QBCP device; a measurement of conductance versus gate voltage at low temperature would distinguish the corrected curves from the old ones and from graphene's $E$-linear conductivity.","Editorial inference: the boundary-condition counting suggests a general recipe for any band structure whose group velocity is not linear: keep all wavevector solutions of the characteristic equation, including complex ones, and match enough derivatives to close the system.","Editorial inference: for tilted or anisotropic quadratic dispersions, the number and form of evanescent modes will change, and the same matching procedure could be run again to produce corrected transmission formulas for those variants."],"forward_implications":["The published $T(E,V_0,\\phi)$, $\\sigma(E,V_0)$, and $F(E,V_0)$ curves for QBCP barriers are superseded by the corrected results in Figs. 2–9.","In two dimensions, the conductivity vanishes and the Fano factor approaches unity as $E\\to 0$, so a QBCP junction becomes sub-Poissonian at low energy while still suppressing transmission below the Sharvin limit.","QBCPs continue to show no Klein tunneling: transmission is not unity at normal incidence, in contrast to Dirac, Weyl, and triple-point fermions.","The appearance of evanescent waves is tied to the quadratic-in-momentum dispersion along the transport direction; the same feature occurs in bilayer graphene, semi-Dirac, and multi-Weyl semimetals, so the method carries over to those problems.","For 3d QBCPs, the calculation shows how to handle the two degenerate conduction bands by assigning separate amplitudes to each, with only the $\\Psi_{+1}$ incident component needed because the cross-amplitude vanishes by symmetry."],"supporting_citations":[{"why":"the earlier treatment of QBCP tunneling whose omission of evanescent waves this paper corrects.","marker":"[11]"},{"why":"introduces the low-energy QBCP Hamiltonian used for the 2d calculation.","marker":"[3]"},{"why":"provides the analogous evanescent-wave solutions for multi-Weyl semimetals, cited to justify the presence of complex wavevectors.","marker":"[25]"},{"why":"gives the semi-Dirac tight-binding model and low-energy behavior, an analogous nonlinear dispersion where evanescent waves appear.","marker":"[27]"},{"why":"shows extra evanescent waves in semi-Dirac Josephson junctions, cited as an analogous situation.","marker":"[28]"},{"why":"supplies the graphene tunneling result used as the no-Klein-tunneling contrast.","marker":"[32]"},{"why":"is the source of the wavefunction-matching procedure used to impose the boundary conditions.","marker":"[33]"},{"why":"provides the Landauer conductance and Fano-factor formulas used to convert transmission into transport observables.","marker":"[34]"}],"fun_headline_variants":["Evanescent waves correct semimetal barrier tunneling","Semimetal tunneling: evanescent waves not optional","Rectangular barrier transmission fixed in semimetals","Quadratic semimetal barrier: evanescent modes matter","Proper evanescent waves reshape semimetal conductance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that requiring both the wavefunction and its derivative to be continuous at the two barrier edges is the complete and correct set of interface conditions for the quadratic Hamiltonian; if the true conditions differ, the solved coefficients and all resulting transmission curves would change.","fun_headline_variants_meta":{"raw":{"variants":["Evanescent waves correct semimetal barrier tunneling","Semimetal tunneling: evanescent waves not optional","Rectangular barrier transmission fixed in semimetals","Quadratic semimetal barrier: evanescent modes matter","Proper evanescent waves reshape semimetal conductance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1341,"prompt_tokens":940,"completion_tokens":401,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":323}},"tokens_in":556,"tokens_out":401,"duration_ms":4148,"temperature":1.0,"reasoning_tokens":323,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T16:12:50.963065+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the same rectangular-barrier scattering problem on a lattice QBCP model by direct numerical integration of the time-independent Schrödinger equation and compare the transmission at normal incidence with the small-$E$ formula $T=4E\\,\\sec^2\\phi\\,\\csc^2(L\\sqrt{2mV_0})/V_0$; a mismatch would show that the assumed derivative-matching conditions are incomplete or incorrect.","supporting_citations":[{"cited_title":"Mandal, Tunneling in Fermi systems with quadratic band crossing points, Annals of Physics 419, 168235 (2020)","cited_arxiv_id":null,"evidence_quote":"the earlier treatment of QBCP tunneling whose omission of evanescent waves this paper corrects."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the low-energy QBCP Hamiltonian used for the 2d calculation."},{"cited_title":"Deng, H.-F","cited_arxiv_id":null,"evidence_quote":"provides the analogous evanescent-wave solutions for multi-Weyl semimetals, cited to justify the presence of complex wavevectors."},{"cited_title":"Banerjee, R","cited_arxiv_id":null,"evidence_quote":"gives the semi-Dirac tight-binding model and low-energy behavior, an analogous nonlinear dispersion where evanescent waves appear."},{"cited_title":"Mandal, Andreev bound states in Josephson junctions of semi-Dirac semimetals, Physica B: Condensed Matter 683, 415918 (2024)","cited_arxiv_id":null,"evidence_quote":"shows extra evanescent waves in semi-Dirac Josephson junctions, cited as an analogous situation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the graphene tunneling result used as the no-Klein-tunneling contrast."},{"cited_title":"Salehi and S","cited_arxiv_id":null,"evidence_quote":"is the source of the wavefunction-matching procedure used to impose the boundary conditions."},{"cited_title":"Tworzyd lo, B","cited_arxiv_id":null,"evidence_quote":"provides the Landauer conductance and Fano-factor formulas used to convert transmission into transport observables."}],"review_version":1}