{"id":"52fd6e1a-dc98-443e-85a5-77b810858931","arxiv_id":"1909.00129","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Using the discrete Wiener-Hopf method, the authors obtain an exact scattered wavefunction for a monoatomic step on square and triangular lattice half-planes, with a far-field approximation validated numerically.","lead":"This paper derives an exact solution for how an electron wave scatters off a single atomic step on the surface of square and triangular lattice crystals, using a nearest-neighbor tight-binding model. The result gives far-field scattering patterns and could inform calculations of surface roughness effects on thin-film resistance.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Far-field asymptotics (4.2) rest on an omitted pole-residue derivation and a contradictory validity interval for T••; the claim that (4.2) matches numerics is not currently verifiable.","rationale":"The paper is a serious and internally coherent attempt to extend the discrete Wiener-Hopf method to step geometries. The exact integral representation (3.1)-(3.2) is supported by a substantial derivation sketch in Appendices D.1 and D.2, and the numerical comparisons are suggestive. The reader's conditional verdict is appropriate because the manuscript omits the derivation of the pole-residue contribution, provides no code or data for the numerical experiments, and relies on several prior-paper results. My stress-test finds the same verification gap to be the most load-bearing: the far-field formula (4.2) is not fully derived, and the manuscript itself flags the missing manipulation in footnote 3. The additional inconsistency between the T•• validity interval in (4.1) and the range stated in Appendix E.2 strengthens the concern, because it leaves unclear whether (4.2) is actually supported for energies such as β−1Eκ=2.52 used in the figures. I am not claiming the result is false; the concrete check of the residue computation and the stationary-phase range would settle the matter. The physical Dirichlet boundary assumption is acknowledged by the authors as a limitation and does not undermine the mathematical claim within the stated model. Since the reader already assigned a conditional verdict, my concern does not move the verdict to a different category, so the verdict remains unchanged.","tokens_in":26081,"tokens_out":27335,"duration_ms":289760,"concrete_test":"Independently compute the residue contributions at zP of the diffraction integral (E.2) for τ=0 (S••) and τ=1 (T••) by expanding the integrand about zP and applying the residue theorem, then compare term-by-term with (4.2c), including the θr branch condition. Separately, rerun the stationary-phase derivation for β−1Eκ∈(7/3,3) using the branch-cut structure of Appendix A.2, checking whether the pole at z=±i alters the saddle-point contour. If the residue formula or the θr criterion differs, recompute the Fig. 13 and Fig. 15 comparison at β−1Eκ=2.52 to determine whether (4.2) still describes the numerical scattered wave.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The far-field approximation (4.2) is a central advertised result, but its reflected-wave components ψ^s_{x,y}|P_s in (4.2c) are flagged in footnote 3 as 'obtained after several manipulations which are omitted.' Because (4.2c) determines the Heaviside shadow-boundary terms, and the numerical comparisons in Figs. 13 and 15 depend on the complete formula (4.2), the claim that (4.2) matches the numerical solution cannot be checked from the text. This is compounded by an internal inconsistency in the stated validity range for T••: (4.1) gives β−1Eκ∈(7/3,3), while Appendix E.2 states the asymptotic analysis follows for β−1Eκ∈(16/3,6). One of these is wrong; if the appendix range is the one actually supported by the cited analysis in [112], then (4.2) is unsupported precisely in the upper-band interval used in Fig. 7(iv), where β−1Eκ=2.52 lies in (7/3,3). The exact integral representation (3.1)-(3.2) is not itself impugned by this concern, but the advertised far-field asymptotics and its graphical validation are.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies scattering of an electronic plane wave by a monoatomic step in a square lattice half-plane and a triangular lattice half-plane, using the nearest-neighbour tight-binding model with Dirichlet boundary conditions on the missing sites. The author formulates the problem as a discrete Wiener-Hopf equation, obtains an exact representation of the scattered wavefunction in (3.1)-(3.2), derives a far-field asymptotic approximation in (4.2), and presents graphical comparisons with a PML-based numerical solution in Figs. 13 and 15. The paper also gives explicit formulas for the step-site wavefunction ψ0,1 for the square lattice and ψ1,1 for the triangular lattice and discusses applications to surface scattering and STM interference patterns.","tokens_in":1593,"tokens_out":1619,"duration_ms":57260,"significance":"If the derivations are correct, the paper provides a useful exact integral-representation solution and a far-field formula for a textbook tight-binding scattering geometry, extending the author's prior Wiener-Hopf analyses of semi-infinite cracks and rigid constraints to the monoatomic step. The main strengths are that no parameter is fitted to the numerics, the numerical solution is an independent check rather than an input, and the unified treatment of square and triangular lattices through the τ parameter is economical. The exact step-site expression (3.3) is an explicit, easily testable result. However, two load-bearing parts of the validation are currently not fully checkable from the manuscript: the derivation of the pole-residue terms in the far-field formula is omitted, and the stated validity interval for that formula is internally contradictory. These issues do not obviously affect the exact integral representation, but they do affect the advertised far-field asymptotics and its numerical verification.","major_comments":[{"comment":"The validity interval for the far-field approximation in the triangular case is stated inconsistently. Equation (4.1) restricts the analysis to β⁻¹Eκ ∈ (−3, 7/3) ∪ (7/3, 3), while Appendix E.2 says the asymptotic analysis follows for β⁻¹Eκ ∈ (−3, 7/3) ∪ (16/3, 6). Since β⁻¹Eκ for the triangular band lies in (−3, 3), the second interval cannot be correct as written. If (16/3, 6) is a typo for (7/3, 3), the error should be corrected; if it is the interval claimed by the supporting analysis in [112], then (4.2) is unsupported in the upper-band interval used in Fig. 7(iv) and Fig. 15(iv), where β⁻¹Eκ = 2.52. Please reconcile the two statements and explicitly confirm that (4.2) is applied only where it is valid.","section":"§4, Eq. (4.1) and Appendix E.2"},{"comment":"The residue contribution ψ^s_{x,y}|_{Ps} is stated in footnote 3 to be 'obtained after several manipulations which are omitted,' and the functions ψrA_{x,y}, ψrB_{x,y}, and the angle θr in (4.2c) are not defined explicitly in the main text (Appendix C identifies ψ^s_{x,y}|_{Ps} with the geometric field ψg only indirectly). Because these terms are part of the total formula (4.2) used in the numerical comparisons in Figs. 13 and 15, the claim that the asymptotic formula matches the numerics cannot be verified from the manuscript. Please provide the derivation of the pole-residue terms, explicit expressions for ψrA, ψrB and θr, and a clear statement of how (4.2c) relates to the geometric solution (C.3).","section":"§4, footnote 3 and Eq. (4.2c)"},{"comment":"The agreement between the numerical solution and the far-field asymptotic approximation is only documented graphically. No error norm, convergence with increasing contour radius R∞, or sensitivity to the PML parameters is reported. To substantiate the central claim of a validated far-field approximation, please include a quantitative error measure (for example, the relative L² or L∞ difference along the discrete contours) and show its behaviour with R∞ for at least the four parameter sets displayed in Figs. 13 and 15.","section":"§5, Figs. 13 and 15"},{"comment":"The multiplicative factors L+(z) and L−(z) are defined only through contour integrals, so (3.1)-(3.2) is an exact integral representation whose practical evaluation still requires further computation. This is not a correctness defect, but the abstract's phrase 'exact solution' should be qualified, and the paper should either outline a numerical strategy for evaluating L± and the outer contour integral in (3.2) or state explicitly that the evaluation follows the procedures of [107,112] with the concrete details needed for reproducibility.","section":"§3 and Appendix B, Eqs. (3.1), (B.1b)"}],"minor_comments":[{"comment":"The one-sided transform notation ψs_{1;−}, ψs_{0;+} is used in (1.8) before the definitions in (A.1) are recalled; consider introducing the notation earlier or adding a pointer.","section":"§1, Eq. (1.8)"},{"comment":"The expressions 'e−e' and 'e+e' are confusing; they should be written as e^{-ε} and e^{+ε} throughout.","section":"§1 and Appendix D, Eqs. (1.10b), (1.11a), (D.1a)"},{"comment":"The Heaviside function H(x) in the boundary-condition discussion is first used without a definition; please specify the convention for H(0), e.g., H(0)=1 if the step site is included.","section":"§1, Eq. (1.6a)"},{"comment":"The sentence 'Because T•• and T••R are \"uncoupled\"' is ambiguous, since the replicated-lattice construction places the two sublattices on a common rectangular grid; please clarify in what sense the two triangular lattices are decoupled in the hopping model.","section":"§2, after Eq. (2.1)"},{"comment":"The strip S is written in terms of ξ1 and ξ2 but these components are not defined; adding ξ = ξ1 + iξ2 with ξ1, ξ2 ∈ R would remove ambiguity.","section":"Appendix E, Eqs. (E.2)-(E.5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is part of a well-established series by the same author on Wiener-Hopf methods for lattice diffraction, and the novelty over [107,108,112,116] is incremental but identifiable (the step geometry and the unified τ parameter). The main concerns are the contradictory validity interval in Section 4 versus Appendix E.2 and the omitted derivation of the pole-residue terms; both are fixable within the paper's scope, so major revision rather than rejection is appropriate. No concerns about citation ethics are raised."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What's actually new: this is the first Wiener-Hopf treatment of a monoatomic step with staggered Dirichlet data on both square and triangular half-planes. The exact integral representation (3.1)-(3.2) is derived from the model, not fitted, and the numerical comparisons in Figs. 13 and 15 are an independent check, not an input. The boundary condition is openly acknowledged in Section 6 as an infinite-barrier idealization, so the physical limitation is not hidden. If you work on discrete lattice diffraction, the square-lattice part is a useful reference solution.\n\nThe soft spots are real and one is load-bearing. First, there is an internal contradiction: (4.1) states the triangular far-field analysis holds for β⁻¹Eκ ∈ (-3,7/3)∪(7/3,3), while Appendix E.2 gives (-3,7/3)∪(16/3,6). These cannot both be right. More seriously, the nearest-neighbor triangular band is [-3,7/3], so both (7/3,3) and (16/3,6) lie outside the bulk pass band. A propagating incident plane wave should not exist there. Fig. 7(iv) uses β⁻¹Eκ = 2.52, which sits in that unphysical interval, so the triangular far-field comparison at that energy cannot be interpreted as bulk scattering. This is not a cosmetic typo; it directly affects the advertised numerical validation.\n\nSecond, the pole-residue contributions in (4.2c) are, by the author's own footnote, obtained after omitted manipulations. Since those terms carry the reflected and shadow-boundary pieces, the claimed match with numerics is not verifiable from the text. Third, no code or data are shipped, the factorization L+ is given only as a contour integral, and the graphical comparisons include no quantitative error estimates. These are smaller issues individually, but together they mean a referee would need to redo substantial parts of the triangular analysis to trust it.\n\nThe square-lattice exact solution stands on firmer ground; the derivation is standard and the omitted steps are more routine there. The central argument is not circular and the author's self-citations are legitimate references to prior method papers.\n\nWho this is for: people working on discrete Wiener-Hopf methods, lattice surface scattering, or thin-film resistivity from an exact-solution angle. The square part deserves citation; the triangular part needs major revision.\n\nRecommendation: send to peer review, but the referee should be told to check the energy-interval contradiction first and require the triangular numerics to be redone inside the true pass band. This is worth referee time, but it is not ready in its current form.","headline":"The square-lattice step solution is a credible, checkable extension of the author's Wiener-Hopf framework, but the triangular-lattice part has a load-bearing internal inconsistency in the energy ranges that undermines the advertised far-field comparison.","tokens_in":26810,"tokens_out":3337,"would_cite":false,"duration_ms":31918,"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":"A bulk electron wave scattered by a single monoatomic step on square and triangular lattice half-planes is solved exactly by a closed contour integral, with a far-field approximation that matches numerics.","keywords":["tight-binding model","Wiener-Hopf method","electron scattering","monoatomic step","lattice half-plane","far-field asymptotics","surface scattering","diffraction"],"falsifier":"An independent, high-resolution numerical solution of the Dirichlet lattice problem (1.6a)–(1.6b) with a different absorbing-boundary scheme can check the Wiener–Hopf formulas (3.1)–(3.2); a separate finite-barrier simulation, comparing $|\\psi_{0,1}|$ versus $\\Theta$ with Figs. 9–11, would test whether the infinite-barrier premise describes real surfaces.","tokens_in":25842,"feed_emoji":"⚛️","tokens_out":9419,"duration_ms":78365,"temperature":0.7,"pith_summary":"This paper establishes that scattering of a bulk electron plane wave by a single monoatomic step on a square or triangular lattice half-plane has an exact solution in the nearest-neighbour tight-binding model. The model treats the surface as an infinitely high barrier, so the wavefunction is exactly zero at the step-adjacent sites. A standard complex-variable technique for half-plane difference equations, the discrete Wiener–Hopf method, reduces the scattered wave to a contour integral, and its far-field approximation to a saddle-point sum plus reflected waves; both are checked against a numerical simulation on a finite grid. If correct, the paper supplies an analytic foundation for computing how conduction electrons scatter off atomic steps, the microscopic input that phenomenological transport models compress into a single specularity parameter.","feed_headline":"Exact solution for electron scattering at a monatomic step","feed_subtitle":"Square and triangular lattice half-planes now have a closed-form scattered wave, verified numerically.","key_machinery":"The load-bearing object is the discrete Wiener–Hopf equation for the one-sided Fourier transform of the scattered field, with kernel $L(z)=\\frac{1}{2}(1+rh/Q)$ for the square lattice and its analogue for the triangular lattice, where $Q$ carries the band structure, $h=\\sqrt{Q-2}$, and $r=\\sqrt{Q+2}$. The kernel is multiplicatively factorized as $L=L_+L_-$, splitting the equation into functions analytic inside and outside a common annulus; Liouville's theorem then fixes the unknown half-range transforms. The scattered field is recovered from the general modal solution $\\psi_y^F=c_1\\lambda^y+c_2\\lambda^{-y}$ with $\\lambda=(r-h)/(r+h)$, and the far field follows by stationary-phase analysis of the resulting diffraction integral. The same machinery handles both lattices, with the triangular case obtained by embedding two copies of the lattice into a rectangular grid.","core_discovery":"The paper's central claim is that the scattered wavefunction $\\psi^s$ at any lattice site is given exactly by the contour integral (3.2), whose integrand is assembled from the factorized Wiener–Hopf kernel and two pole terms; with $\\tau=0$ for the square lattice and $\\tau=1$ for the triangular lattice, the same closed-form expression (3.1) covers both geometries. The far-field approximation (4.2) represents the result as the sum of saddle-point contributions, the diffracted cylindrical wave, plus two reflected plane-wave terms that switch on and off at a critical angle $\\theta_r$, and this asymptotic formula is shown to match the numerical solution on a finite grid with absorbing boundaries. The paper maintains that the solution is exact within the tight-binding model, and that the only essential physical idealisation is the Dirichlet boundary condition at the step.","pith_inferences":["The factorization of the solution into a lattice-dependent factor $K(z)$ and two pole terms suggests that a step-scattering amplitude could be defined as a measurable quantity, giving an angle- and energy-resolved specularity parameter that interpolates between $p=1$ and $p=0$.","A finite-barrier version of the problem would differ from this solution near the step; since the paper identifies the infinite barrier as essential, a perturbative expansion around the Dirichlet solution is the natural next test of how much the hard-wall idealisation matters.","For a periodic array of identical steps, the single-step scattering phase could be combined with Bloch conditions to predict coherent interference features, such as minigaps, in the surface conductance of vicinal surfaces."],"forward_implications":["The exact contour-integral formula (3.2) gives the scattered wavefunction at every lattice site, so no uncontrolled truncation is needed for a single step under the Dirichlet condition.","The far-field formula (4.2) provides explicit angle- and energy-dependent scattering amplitudes, which can be inserted into convolution schemes for rough surfaces to replace the single specularity parameter $p$.","For low incident energies the kernel $L$ approaches 1, and the step effectively acts as a flat geometric mirror; the numerically observed deviations are attributed to imperfect absorbing boundaries rather than model error.","The same Wiener–Hopf construction extends to honeycomb lattices and to lattice waveguides with steps, as the paper states for future work."],"supporting_citations":[{"why":"Provides the Wiener–Hopf method and multiplicative factorization of the kernel that convert the lattice difference equation into an exactly solvable half-plane problem.","marker":"[84]"},{"why":"Supplies the discrete Fourier transform formulation, radiation conditions, and stationary-phase analysis for square-lattice diffraction that this paper adapts.","marker":"[107]"},{"why":"Introduces the replicated triangular lattice construction, the kernel Q(z), and the triangular-lattice saddle-point analysis used for the T•• case.","marker":"[112]"},{"why":"Gives the Wiener–Hopf treatment of a square-lattice half-plane with mixed boundary conditions, providing the structural template for the step problem.","marker":"[116]"},{"why":"Furnishes the physical justification for the vanishing wavefunction at the metal boundary, i.e., the infinite-barrier approximation of the surface.","marker":"[2]"},{"why":"Presents the analogous exact solution for the square lattice with a semi-infinite rigid constraint, whose asymptotic machinery is reused for the scattered far field.","marker":"[108]"},{"why":"Documents the transform manipulations and general solution constructions on which the scattered-wave representation (1.7) is based.","marker":"[121]"}],"fun_headline_variants":["Exact electron scattering at a step, two lattices","Closed-form solution for step scattering, square & triangular","One exact wavefunction for step scattering on two lattices","Step scattering solved exactly for square and triangular"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the electron wavefunction vanishes exactly on the step-adjacent sites, modelling the surface as an infinitely high potential barrier; if a real surface has a finite barrier, the computed scattering amplitudes do not describe it.","fun_headline_variants_meta":{"raw":{"variants":["Exact electron scattering at a step, two lattices","Closed-form solution for step scattering, square & triangular","One exact wavefunction for step scattering on two lattices","Step scattering solved exactly for square and triangular"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001276,"raw_usage":{"total_tokens":5159,"prompt_tokens":824,"completion_tokens":4335,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":440,"completion_tokens_details":{"reasoning_tokens":4272}},"tokens_in":440,"tokens_out":4335,"duration_ms":538051,"temperature":1.0,"reasoning_tokens":4272,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T06:00:52.090440+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent, high-resolution numerical solution of the Dirichlet lattice problem (1.6a)–(1.6b) with a different absorbing-boundary scheme can check the Wiener–Hopf formulas (3.1)–(3.2); a separate finite-barrier simulation, comparing $|\\psi_{0,1}|$ versus $\\Theta$ with Figs. 9–11, would test whether the infinite-barrier premise describes real surfaces.","supporting_citations":[{"cited_title":"On scattering of waves on square lattice half-plane with mixed boundary condition","cited_arxiv_id":null,"evidence_quote":"Gives the Wiener–Hopf treatment of a square-lattice half-plane with mixed boundary conditions, providing the structural template for the step problem."},{"cited_title":"Diffraction of waves on square lattice by semi-infinite rigid constraint","cited_arxiv_id":null,"evidence_quote":"Presents the analogous exact solution for the square lattice with a semi-infinite rigid constraint, whose asymptotic machinery is reused for the scattered far field."}],"review_version":1}