{"id":"3a7a23f2-d8b7-4b20-958c-abde1342c554","arxiv_id":"2412.14255","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Zero modes with Gaussian, exponential, or polynomial decay are exact solutions when the mass diverges, stays finite, or vanishes at infinity.","lead":"By inverse construction, this paper writes down exact zero-energy wavefunctions and the mass and velocity profiles that produce them, for Gaussian, exponential, and polynomial decay in modified Jackiw-Rebbi and Schrödinger equations, including complex nonhermitian profiles. The result maps each decay class to how the mass behaves at infinity, which is a useful input for distinguishing topological Majorana modes from trivial Andreev states.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim is overstated: Gaussian faster-than-exponential modes occur with constant nonzero mass and linearly diverging velocity (Eqs. 14 and 18), so the Sec. IX statement that faster-than-exponential decay implies divergent mass is false without a bounded-velocity qualification.","rationale":"Good faith reading: the paper is an inverse-method catalogue, not a proof of a no-go theorem. The central construction via Eqs. (5)-(7) is correct, and all three main examples are genuine exact solutions of Eq. (3) or Eq. (4) for the displayed fields. Credit is due for explicit formulas and for the honest limitations in Secs. V and VIII about special coefficient choices. However, the paper advertises a general classification in the abstract and Sec. IX. The single most load-bearing weak point is that the advertised mapping from mass asymptotics to decay class is not universally valid; the authors' own constant-mass Gaussian example is a counterexample unless the statement is restricted to bounded or uniform velocity. This is not a failure of the inverse method but an overreach in the central claim. The reader's weakest assumption concerning ansatz-completeness and the periodic-mode domain issue is related; in particular, the periodic modes in Sec. VIII F are normalizable only on [-pi,pi] and cannot be extended to normalizable modes on R, so that example should be presented as a finite-interval or finite-ring construction rather than an infinite-line array. The appropriate verdict remains CONDITIONAL: the construction is sound, but the central claim and the periodic-mode interpretation require qualification before acceptance.","tokens_in":20904,"tokens_out":9099,"duration_ms":89141,"concrete_test":"Re-evaluate the Gaussian mode of Sec. VIII A with the constant-mass choice m=m_*=-beta and the velocity sv(x)=(beta x+alpha)/2 from Eq. (18). Substitute phi(x)=exp(-beta x^2/2-alpha x) into Eq. (3) and verify that phi''+2s v phi'-m phi=0 identically while m is constant and nonzero and phi decays Gaussian. If the identity holds, the Sec. IX sentence is false as written and must be qualified by the companion-field asymptotics or the assumption of bounded velocity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IX states: 'faster-than-exponential modes correspond to the divergence of the mass term, exponential modes correspond to finite and nonzero asymptotic masses, while slower-than-exponential modes correspond to the mass vanishing at large distances.' This unqualified mapping is contradicted by the paper's own Gaussian example. Taking phi(x)=exp(-beta x^2/2 - alpha x) (Eq. 14), m(x)=m_*=-beta, and sv(x)=(beta x+alpha)/2 (Eq. 18), Eq. (3) is satisfied exactly: phi''+2sv phi'-m phi=0. Here phi decays faster than exponentially, yet m is constant and nonzero, and v diverges linearly. Thus the asymptotic behavior of the mass term alone does not determine the localization class; the velocity profile is equally decisive. The earlier summary in Sec. II contains the necessary qualification 'For uniform Dirac velocity,...', but the abstract and Sec. IX state the classification without it. Additionally, the three ansatz families are not shown to exhaust all possible decay classes, and the paper itself notes (after Eqs. 21-23 and 33-35) that these ansaetze solve only special coefficient choices. The inverse construction is mathematically sound as a solution-generating method; what is not established is the general classification claimed in the conclusions. The central claim therefore needs to be reformulated as a statement about the inverse-constructed families with specified companion-field asymptotics, not about zero modes in general.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an inverse-method approach to construct exact zero-energy solutions of the modified Jackiw-Rebbi equation (Eq. 3) and eigenstates of the associated Schrödinger equation (Eq. 4). For three families of wavefunctions with Gaussian (faster-than-exponential), exponential, and polynomial (slower-than-exponential) decay, the authors derive the corresponding mass m(x), Dirac velocity v(x), and potential V(x) that make those wavefunctions exact solutions. They also construct second linearly independent modes with the same pseudospin via reduction of order (Eq. 10), modes with opposite pseudospin when v(x) is constant (Eq. 11), and fields supporting two modes at finite separation (Eqs. 12-13). Periodic analogues are treated in Sec. VIII F, and the constructions are extended to non-Hermitian fields. The central classification claim, stated in the abstract and Sec. IX, is that faster-than-exponential decay corresponds to a divergent mass term, exponential decay to finite nonzero asymptotic mass, and slower-than-exponential decay to a mass vanishing at infinity.","tokens_in":21288,"tokens_out":3158,"duration_ms":29637,"significance":"If the central claims hold, the paper provides a useful catalogue of exact analytical zero-mode solutions with controlled asymptotic behavior, extending the authors' earlier smooth-domain-wall results to Gaussian and power-law localization. The explicit formulas for companion fields and for second linearly independent modes are valuable and are, as far as the displayed algebra shows, correct. Several of the constructed pairs are consistent with bulk-boundary correspondence expectations, and the non-Hermitian extensions are a distinct contribution. The main weakness is that the general classification in Section IX is stated more broadly than what the inverse construction actually proves; this is fixable but currently affects the paper's central message.","major_comments":[{"comment":"The unqualified statement that 'faster-than-exponential modes correspond to the divergence of the mass term' is contradicted by the paper's own Gaussian example. Taking phi(x) = exp(-beta x^2/2 - alpha x) with m(x) = m_* = -beta and s v(x) = (beta x + alpha)/2, Eq. (18), satisfies Eq. (3) exactly, yet m is constant and nonzero while the velocity diverges linearly. The qualification 'for uniform Dirac velocity' appears in Section II but is missing from the abstract and Section IX. The classification should be restated as a joint statement about the asymptotic behavior of both m(x) and v(x), or explicitly restricted to the uniform-velocity case.","section":"Section IX and abstract; Eqs. (14)-(18)"},{"comment":"The conclusions present the decay-class/mass-asymptotics correspondence as a general classification of zero modes, but the paper analyzes only three ansatz families and explicitly notes that these wavefunctions solve Eq. (3) only for special coefficient choices, not for general polynomial fields. No theorem is given that every normalizable zero mode of Eq. (3) has one of the Gaussian, exponential, or polynomial asymptotics, and mixed decays such as stretched exponentials are not discussed. The classification should be scoped to the inverse-constructed families with stated companion-field asymptotics, not to all zero modes of the modified Jackiw-Rebbi equation.","section":"Section IX; also Eqs. (21)-(23) and (33)-(35)"},{"comment":"The periodic modes in Eqs. (60)-(62) are described as modes localized on an infinite equally spaced array of points. These wavefunctions are periodic and bounded but not normalizable on the real line; they are normalizable only on a finite interval such as -pi <= x <= pi or on a circle. On the infinite line they do not decay, so the assertion that they are localized zero modes on an infinite array is unsupported. The domain of integration and the boundary conditions for these solutions should be stated explicitly, and the corresponding claim in Section IX should be qualified.","section":"Sec. VIII F and Sec. IX (periodic modes)"}],"minor_comments":[{"comment":"There are several typographical errors, including 'Schödinger' for 'Schrödinger' in Section II and other headings, 'satysfing' in Section VII, 'wronksian' in Section V, and 'Schrôdinger' in Section VIII C.","section":"Throughout"},{"comment":"The argument of the erfi function in Eq. (20) appears inconsistent with the preceding expression and with the statement that this mode decays as approximately 1/x; the normalization and argument should be checked.","section":"Eq. (20)"},{"comment":"The phrase 'at large distances |x| > infinity' should read 'at large distances |x| -> infinity'.","section":"Section II"},{"comment":"The inverse method is definitional: Eq. (5) constructs m(x) so that the chosen phi solves Eq. (3) by construction. Stating this explicitly would help readers distinguish the exact solution-generating nature of the method from a predictive derivation of mode asymptotics from general field profiles.","section":"Section IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a natural extension of the authors' Ref. [40], with new content centered on decay classes and non-Hermitian cases. The overstatement in the abstract and Section IX and the periodic-mode domain issue are the main obstacles; both are local fixes that do not require new calculations of a different type. If the authors scope the classification appropriately, the paper should be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWorth your time, but the headline classification is overstated. The inverse method here is transparent and correct: for any chosen wavefunction, Eqs. (5)-(7) produce fields that make it an exact solution of the modified Jackiw-Rebbi or Schrödinger equation. The explicit Gaussian, exponential, polynomial, and periodic families—plus the nonhermitian extensions—are genuinely new relative to the authors' earlier exponential-domain-wall paper. I checked the reduction-of-order and opposite-pseudospin relations; they work. If you work on Majorana nanowires or Jackiw-Rebbi-type models, this is a clean catalog connecting decay shape to asymptotic mass behavior, and the STM suggestion is reasonable.\n\nThe soft spot is the Sec IX conclusion. It states without qualification that faster-than-exponential modes correspond to divergent mass, exponential to finite nonzero mass, slower-than-exponential to vanishing mass. The paper itself contradicts this in the Gaussian example: with m = -beta constant and v(x) = (beta x + alpha)/2, phi ~ exp(-beta x^2/2) solves Eq. (3) exactly, with constant nonzero mass and a linearly diverging velocity. So the decay class is not determined by the mass term alone; the velocity behavior is equally decisive. The Sec II summary has the necessary 'For uniform Dirac velocity' qualifier, but the abstract and conclusions drop it. That is a fixable but real flaw.\n\nThe periodic modes are a second, minor issue: they are normalizable only on -pi <= x <= pi, while the operators are treated on the infinite line. Calling them localized on an infinite array is unsupported without defining the periodic extension.\n\nThe paper does not prove the classification is exhaustive; the ansatz families are special cases, which the authors acknowledge. That limits the framing to a catalog rather than a theorem, but the catalog itself is sound.\n\nI'd send it to review. A serious referee should require the classification to be reformulated with the velocity qualification, and the periodic-mode domain to be clarified. Conditional acceptance after moderate revision.","headline":"Useful exact-solution catalog, but the Sec IX classification overreaches: constant nonzero mass plus linearly diverging velocity already gives a Gaussian mode.","tokens_in":21787,"tokens_out":2716,"would_cite":true,"duration_ms":23873,"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":"The paper claims that faster-than-exponential, exponential, and slower-than-exponential zero modes are classified by whether the mass diverges, stays finite, or vanishes at large distances, with exact solutions constructed by inverse…","keywords":["zero-energy modes","Jackiw-Rebbi equation","Majorana zero modes","topological insulators","topological superconductors","inverse methods","non-Hermitian physics","decay localization"],"falsifier":"Find a normalizable zero mode of $\\varphi''(x)+2s v(x)\\varphi'(x)-m(x)\\varphi(x)=0$ with $m(x)$ tending to a finite nonzero constant as $|x|\\to\\infty$ but with slower-than-exponential polynomial decay; the central claim predicts that such a mode cannot exist.","tokens_in":1907,"feed_emoji":"⚛️","tokens_out":7705,"duration_ms":119403,"temperature":0.7,"pith_summary":"Zero-energy boundary modes in topological insulators, superconductors, and superfluids are usually pictured as exponentially localized, but the same modified Jackiw-Rebbi equation also admits Gaussian and polynomially decaying modes. This paper tries to establish that the decay class is controlled by the asymptotic behavior of the mass term: a mass that diverges at large distances gives faster-than-exponential decay, a finite nonzero mass gives exponential decay, and a mass that vanishes gives polynomial decay. Using inverse methods, the authors take each candidate wavefunction and derive the exact position-dependent mass, Dirac velocity, and Schrödinger potential that admit it, so every displayed mode is a genuine zero-energy solution. The scheme extends to non-Hermitian equations and gives conditions under which real modes can coexist with complex fields. If correct, the classification provides a practical way to distinguish topologically protected exponential Majorana modes from disorder-induced near-zero modes with weaker localization.","feed_headline":"Mass at infinity dictates the decay of zero-energy modes","feed_subtitle":"Exact Gaussian, exponential, and polynomial zero modes classified by the far-field mass.","key_machinery":"The machinery is the inverse method for the modified Jackiw-Rebbi equation. Instead of solving the second-order equation $\\varphi''+2sv\\varphi'-m\\varphi=0$ for $\\varphi$ given $m$ and $v$, one fixes a normalizable $\\varphi$ and algebraically recovers the fields from $m(x)=\\varphi''/\\varphi+2sv\\,\\varphi'/\\varphi$ and $V(x)-E=\\varphi''/\\varphi$; reversing the roles gives $sv(x)=-(\\varphi''-m\\varphi)/(2\\varphi')$. A reduction-of-order formula $\\varphi^s_2(x)=\\varphi^s_1(x)\\int^x e^{-2s\\int^y v\\,dz}/(\\varphi^s_1(y))^2\\,dy$ produces the second same-pseudospin solution, and for constant $v$ the duality $\\varphi^{-s}=e^{2svx}\\varphi^{s}$ gives the opposite-pseudospin partner. A two-mode version solves for $m(x)$ and $v(x)$ from a pair of prescribed modes. This inversion is what lets the paper assert that each displayed Gaussian, exponential, or polynomial wavefunction is an exact solution rather than an approximation, and it is what makes the asymptotic mass the quantity that selects the decay class.","core_discovery":"The central discovery is a three-way correspondence between the asymptotic behavior of the mass term $m(x)$ and the localization class of zero modes of $\\varphi''(x)+2s v(x)\\varphi'(x)-m(x)\\varphi(x)=0$. Gaussian modes $\\varphi\\sim e^{-\\beta x^2/2}$ are exact zero modes when the mass diverges quadratically at infinity, or when the Dirac velocity diverges linearly; exponentially decaying modes $\\varphi\\sim e^{-\\alpha x}/\\cosh^\\gamma(\\beta x)$ correspond to finite nonzero asymptotic masses; polynomially decaying modes $\\varphi\\sim(\\beta^2 x^2+1)^{-n/2}$ correspond to a mass that dips and vanishes at infinity, with the localization length diverging as the topological gap closes. The same inverse construction yields the Schrödinger potential $V(x)-E=\\varphi''/\\varphi$, produces a second linearly independent mode by reduction of order, and extends to non-Hermitian fields where real eigenmodes can coexist with complex masses, velocities, or potentials. The authors also construct pairs of modes centered at finite distance and periodic modes peaked on an evenly spaced array, and they argue that the presence and number of zero modes still follow the bulk-boundary correspondence even when the localization law is set by the far-field data.","pith_inferences":["A testable extension is to use the same inverse construction with stretched-exponential trial wavefunctions $e^{-|x|^\\nu}$; the paper's logic predicts the mass would have to decay algebraically or logarithmically, and finding such families would interpolate between the polynomial and exponential classes.","Because the periodic modes are normalized on the finite interval $-\\pi\\le x\\le\\pi$ while the operators are treated on the infinite line, a reader should regard the infinite-array interpretation as an open question; a fully periodic boundary treatment would settle whether these are true localized zero modes on an infinite lattice.","The classification suggests a diagnostic for trivial versus topological bound states: measure the spatial decay exponent of near-zero-energy peaks, with exponential decay and finite asymptotic mass as the protected signature and Gaussian or polynomial profiles indicating a confining trap or a vanishing gap rather than a genuine topological transition.","Extending the inverse method to lattice models or to dispersions with higher-order momentum terms could reveal whether the mass-asymptotic correspondence persists when the quadratic-momentum approximation is abandoned."],"forward_implications":["If the mass profile of a nanowire or superconductor is measured or engineered, the decay profile of its zero-energy modes is fixed, so the classification can be checked directly in local-density-of-states experiments.","Gaussian zero modes arise naturally in harmonic traps, so cold-atom superfluids and harmonic-shell models of nuclei are settings where the predicted divergent mass should produce superexponential Majorana-like modes.","The bulk-boundary correspondence still dictates the number of zero modes even when their localization is faster or slower than exponential; only the decay law, not the mode count, depends on the far-field behavior of the mass and velocity.","In the non-Hermitian regime, real zero modes can survive even when the mass, velocity, or potential is complex, provided the appropriate reality condition such as $V(x)-E\\in\\mathbb{R}$ holds.","Periodic zero modes localized on an equally spaced array of points are exact solutions, offering a mechanism for zero-energy modes at multiple sites without the energy splitting normally expected from overlapping modes."],"supporting_citations":[{"why":"Supplies the original Jackiw-Rebbi equation that the paper modifies with a quadratic momentum term.","marker":"[34]"},{"why":"Establishes the modified Jackiw-Rebbi equation, the duality to Schrödinger eigenmodes, and the smooth-domain-wall zero modes this paper extends.","marker":"[40]"},{"why":"Defines the $\\mathbb{Z}_2$ topological invariant used to read off the expected number of zero modes from asymptotic fields.","marker":"[13]"},{"why":"Provides the bulk-boundary correspondence used to interpret mode count and localization at interfaces.","marker":"[35]"},{"why":"Supports the topological origin of zero-energy edge states in particle-hole symmetric systems, used alongside the correspondence.","marker":"[36]"},{"why":"Generalizes the bulk-boundary correspondence to topological defects, the framing used for domain walls with spatially varying mass and velocity.","marker":"[37]"},{"why":"Supplies the harmonic-trap cold-atom setup where Gaussian zero modes are physically realizable.","marker":"[41]"},{"why":"Gives the harmonic-trap Majorana corner-mode context that motivates superexponentially localized modes.","marker":"[42]"},{"why":"Documents near-zero-energy bound states of disorder origin that lack exponential localization, the class the classification distinguishes from protected modes.","marker":"[48–63]"},{"why":"Proposed experimental probe: STM imaging of Majorana wavefunctions in atomic chains, which could measure the decay profile of zero modes.","marker":"[69]"}],"fun_headline_variants":["Zero-mode decay dictated by mass at infinity","Exact zero modes: Gaussian, exponential, or polynomial decay","Mass behavior picks zero-mode localization law","Far-field mass sets zero-mode decay: exponential, Gaussian, polynomial"],"cache_read_input_tokens":23808,"weakest_assumption_plain":"The classification rests on the three chosen ansatz families, and the paper does not prove that every normalizable zero mode of Eq. (3) belongs to one of those families.","fun_headline_variants_meta":{"raw":{"variants":["Zero-mode decay dictated by mass at infinity","Exact zero modes: Gaussian, exponential, or polynomial decay","Mass behavior picks zero-mode localization law","Far-field mass sets zero-mode decay: exponential, Gaussian, polynomial"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0004,"raw_usage":{"total_tokens":2117,"prompt_tokens":999,"completion_tokens":1118,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":1055}},"tokens_in":615,"tokens_out":1118,"duration_ms":8887,"temperature":1.0,"reasoning_tokens":1055,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:24:03.872388+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find a normalizable zero mode of $\\varphi''(x)+2s v(x)\\varphi'(x)-m(x)\\varphi(x)=0$ with $m(x)$ tending to a finite nonzero constant as $|x|\\to\\infty$ but with slower-than-exponential polynomial decay; the central claim predicts that such a mode cannot exist.","supporting_citations":[{"cited_title":"Marra and A","cited_arxiv_id":null,"evidence_quote":"Establishes the modified Jackiw-Rebbi equation, the duality to Schrödinger eigenmodes, and the smooth-domain-wall zero modes this paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the $\\mathbb{Z}_2$ topological invariant used to read off the expected number of zero modes from asymptotic fields."},{"cited_title":"Teo and C","cited_arxiv_id":null,"evidence_quote":"Generalizes the bulk-boundary correspondence to topological defects, the framing used for domain walls with spatially varying mass and velocity."},{"cited_title":"Zhou and Z","cited_arxiv_id":null,"evidence_quote":"Supplies the harmonic-trap cold-atom setup where Gaussian zero modes are physically realizable."},{"cited_title":"Kheirkhah, Y","cited_arxiv_id":null,"evidence_quote":"Gives the harmonic-trap Majorana corner-mode context that motivates superexponentially localized modes."},{"cited_title":"Pawlak, M","cited_arxiv_id":null,"evidence_quote":"Proposed experimental probe: STM imaging of Majorana wavefunctions in atomic chains, which could measure the decay profile of zero modes."}],"review_version":1}