{"id":"351936db-6c50-47a9-b617-4e849314fa7d","arxiv_id":"2505.20180","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For fermions on an Abrikosov-Nielsen-Olesen vortex, the vacuum fermion number equals [sgn(m+e sqrt(2) v)+sgn(m-e sqrt(2) v)] n/4, and disk edge states carry charge e/2.","lead":"This paper calculates the quantum fermion number of a vortex with charged fermions, finding a step-like dependence on the fermion mass. The calculation tests a recently proposed heat-kernel method and predicts fractionally charged edge states, which may matter for topological materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Subtraction identity (19) is the load-bearing step; continuum phase shifts at the bag boundary could contaminate the eta invariant and alter Eq (52).","rationale":"The reader's weakest_assumption identifies precisely the step I consider most load-bearing: the unproven subtraction formula (19). The paper's logic is coherent once (19) is accepted: the variational equations are consistent, the edge-state subtraction is carefully done, and the integration constants are anchored to known limits. However, none of the internal consistency checks distinguishes edge-state contributions from continuum boundary contributions, so the central formula (52) remains conditional on (19). I do not see a new fatal flaw beyond what the reader flagged; the same concern is sharpened by noting that in the no-edge-state regime the method implicitly assumes a zero continuum boundary contribution, which is exactly the kind of claim that scattering-phase analysis can test. Therefore the reader's CONDITIONAL verdict should stand, and my stress-test does not move it.","tokens_in":9766,"tokens_out":18484,"duration_ms":222844,"concrete_test":"Consider the near-boundary Hamiltonian (27) with constant m, ṽ = e√2 v, and ã, on a half-line with bag condition (10). Compute the full η invariant of the half-line problem exactly via scattering phase shifts, including both bound and continuum states, and compare the boundary-induced quantity η(half-line) − η(full line) with the η of the effective edge Hamiltonian G(g) given in (43)–(44), over all parameter regimes, especially |m| > ṽ with ε = +sgn(m − ṽ) where no edge states exist. If continuum phase shifts produce any nonvanishing contribution, Eq. (19) fails already in the local approximation and Eq. (52) is not established; if the two quantities agree exactly, the subtraction formula is validated in the regime used.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result (52) rests on the subtraction identity (19), N_R(H) = N(H_DR) − N(H_b), and its variational counterpart (21). The authors state explicitly that this is 'natural to assume (though hard to prove rigorously)'. If the boundary-induced spectral asymmetry is not entirely carried by edge states, then subtracting only N(H_b) leaves a residual contribution and Eq. (52) does not follow from the preceding variational calculation. The danger is concrete: a bag boundary condition also modifies the continuum scattering states near the boundary, producing a boundary layer in the density of states even in regimes with no normalizable edge mode, e.g. |m| > √2 e v with ε = +sgn(m − √2 e v), where the paper sets N(H_b) = 0. A continuum phase-shift contribution to η(0,H_DR) − η(0,H_plane) would violate (19) already in the local near-boundary model. The consistency checks in Sec. 4 (ε-independence, the v → 0 limit, and the m = 0 symmetry) are all compatible with (52), but they do not isolate the boundary contribution and therefore do not test (19) directly.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the vacuum fermion number for charged fermions coupled to an Abrikosov-Nielsen-Olesen vortex in 2+1 dimensions, using the method of Fresneda-de Souza-Vassilevich [arXiv:2305.13606]. The method places the soliton on a disk with bag boundary conditions, computes the variation of the eta invariant through a heat-kernel boundary coefficient, and subtracts the contribution of edge states. On the vortex background the authors obtain Eq. (52): N(H) = [sgn(m+√2 e v)+sgn(m−√2 e v)] n/4, which is independent of the bag boundary parameter ε, reduces to the Niemi-Semenoff result when v→0, and is consistent with the m=0 index theorem. The paper also predicts n fermionic zero modes in the intermediate mass range |m|<√2 e v and claims that the edge states carry fractional charges e/2 or continuously varying charges.","tokens_in":9984,"tokens_out":11968,"duration_ms":134437,"significance":"If Eq. (52) is correct, the paper provides a compact closed-form expression for the vortex-induced fermion number that goes beyond the topological Jackiw-Rebbi/Niemi-Semenoff results and predicts half-integer jumps at |m|=√2 e v. The work also serves as a test of the method of [1], and the demonstration of boundary-condition independence is a nontrivial check. The derivation is explicit, reproduces two independent known limits (v→0 and m=0), and may be useful for zero-mode counting and fractional-charge applications. The main weakness is that the central subtraction identity (19) is explicitly unproven, and the boundary-condition independence checks do not isolate the boundary-layer contribution; in addition, the integration constants leading to (52) are not fully fixed by the limits presented in the paper. The result is therefore plausible but not yet established at the level the central claim requires.","major_comments":[{"comment":"The identity N_R(H)=N(H_DR)−N(H_b) is the load-bearing step of the calculation, and the manuscript explicitly states that it is \"natural to assume (though hard to prove rigorously)\". The consistency checks in Sec. 4 compare only the final total fermion number, so they cannot isolate the boundary contribution: independence of ε would hold even if both N(H_DR) and N(H_b) contained a common spurious boundary-layer contribution from continuum modes. The previously checked magnetic-field example in [1] does not cover the scalar and mass couplings present here. A concrete test would be to compute the continuum phase shifts at the bag boundary and verify directly that η(0,H_DR)−η(0,H_plane) equals −2N(H_b), at least in a regime with no normalizable edge states, e.g. |m|>√2 e v with ε=sgn(m−√2 e v), where the paper sets N(H_b)=0.","section":"Section 2.3, Eq. (19)"},{"comment":"The claim that the three integration constants are fixed by m→±∞ (equivalently, v→0) and m=0 is not supported by Sec. 4. The v→0 limit collapses the intervals m>√2 e v and 0<m<√2 e v into a single region for fixed m>0 (and similarly for negative m), so it cannot separately fix the constants on the two sides of each threshold. The m→∞ limit is not evaluated; the discussion only shows absence of zero modes for large |m|, which does not determine η(0,H). Thus the constants leading to Eq. (52) are not actually derived from the stated limits, and Eq. (52) is not uniquely determined by the variational equation (50) alone.","section":"Section 3.3, text after Eq. (50)"}],"minor_comments":[{"comment":"The notation A∥ and dx∥ is introduced only implicitly; please define A∥ explicitly with the chosen orientation and state the relation between ∫_{S_R} dθ√h δA_j ϵ_{nj} and ∫_{S_R} dx∥ δA∥, to avoid a factor-of-R ambiguity.","section":"Equations (24) and (46)"},{"comment":"The sentence \"We have neglected the corrections to (34) which vanish exponentially fast at R→∞\" should specify the order of the neglected terms and state that the bound is uniform in the tangential momentum κ that contributes to the variation.","section":"Section 3.2, near Eq. (34)"},{"comment":"The index formula (55) is quoted without stating the domain assumptions needed for the index theorem on R²; a sentence clarifying the functional setting (e.g., decay of zero modes) would make the argument more precise.","section":"Section 4, m=0 paragraph"},{"comment":"There are several presentation issues: \"V ortex\" in the Section 3 heading, \"m → ±∞(equivalently\" missing a space, the acknowledgments contain an encoding artifact \"S˜ ao\", and the components of w± in Eq. (36) are split across lines in a way that makes the column vectors hard to read.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The method being tested is the authors' own from arXiv:2305.13606, and the subtraction identity (19) is exactly the core assumption of that method. The independent evidence for (19) in the present paper is limited to global consistency checks; a direct phase-shift or large-mass computation would substantially strengthen the manuscript. The integration-constant gap after Eq. (50) is the most fixable but also the most immediate technical objection. The paper fits the journal's scope and the central claim is defensible, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know up front. The main result is new: for a charged Dirac fermion coupled to the ANO vortex, the vacuum fermion number is N = [sgn(m+e√2v)+sgn(m−e√2v)] n/4 — two mass thresholds, a plateau in the middle (N=0, with n zero modes predicted), and correct reduction to the known limits: Niemi-Semenoff when v→0 and the index theorem at m=0. That is a clean answer to a question the literature only settles piecewise.\n\nNow the soft spot, and it is the load-bearing one. The derivation rests on Eq. (19), which identifies the entire near-boundary fermion density with the edge-state contribution. The authors say it themselves: natural to assume, hard to prove rigorously. The worry is concrete. A bag boundary also modifies the continuum scattering states near the boundary, and scattering phase shifts produce spectral asymmetry even when no normalizable edge state exists. In the |m|>√2 e v regime with ε=+sgn(m−√2 e v), the paper sets the edge contribution to zero and attributes nothing to the continuum. If that continuum piece is nonzero, Eq. (52) does not follow from the variational calculation.\n\nTheir consistency checks are real but do not close that gap. The ε-independence is a genuine check — the edge-state spectra differ between the two boundary conditions and the subtraction compensates exactly. But an ε-independent continuum contamination would survive both cases unseen. The v→0 limit is the strongest external anchor, and the m=0 index argument is standard. The intermediate plateau, which is the genuinely new claim, has no independent verification.\n\nThe rest is in better shape. The edge-state spectrum computation is explicit; the fractional charges (e/2 in the Higgs phase, continuously varying above the threshold) are a solid by-product with a plausible condensed-matter hook. The method being tested is the authors' own from arXiv:2305.13606, so the test is partly self-referential, though on a new system with independent anchors. The integration constants in (51) are fixed by external limits rather than derived — legitimate, but worth saying.\n\nWho should read it: people working on fractional fermion number, vortex zero modes, or heat-kernel methods. It deserves a serious referee; the referee's job is to press on (19) or demand an independent check of the plateau. Send it to review, with the outcome conditional on that.","headline":"A genuinely new closed-form fermion number for the ANO vortex with a novel mass-threshold plateau, but the derivation's load-bearing subtraction (Eq. 19) is explicitly unproven, so the result is conditional rather than established.","tokens_in":10520,"tokens_out":6849,"would_cite":true,"duration_ms":70668,"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":"The paper establishes that the vacuum fermion number on an Abrikosov–Nielsen–Olesen vortex takes the boundary-independent value $N(H) = [\\operatorname{sgn}(m+e\\sqrt{2}v)+\\operatorname{sgn}(m-e\\sqrt{2}v)]\\,n/4$, and uses this to test a…","keywords":["Abrikosov-Nielsen-Olesen vortex","eta invariant","heat kernel expansion","edge states","fermion number fractionization","zero modes","bag boundary conditions","fractional charge"],"falsifier":"Compute the spectral asymmetry directly, for instance by numerically diagonalizing the Dirac Hamiltonian on large disks over a range of masses, and examine whether the plane-limit $\\eta(0,H)$ is independent of the boundary parameter $\\varepsilon$ and jumps by $-n$ exactly at $m=\\pm e\\sqrt{2}v$. Any boundary-condition dependence in the large-$R$ limit, or a jump at a different mass, would refute Eq. (52).","tokens_in":9580,"feed_emoji":"🌀","tokens_out":7029,"duration_ms":67068,"temperature":0.7,"pith_summary":"The paper aims to compute the vacuum expectation value of the fermion number on an Abrikosov–Nielsen–Olesen vortex, for charged fermions coupled to the Higgs and gauge fields. It establishes the exact formula $N(H) = [\\operatorname{sgn}(m+e\\sqrt{2}v)+\\operatorname{sgn}(m-e\\sqrt{2}v)]\\,n/4$, where $n$ is the vortex winding number. This matters because it is a nontrivial test of a recent method that computes fermion number from a heat kernel coefficient plus edge-state contributions, and because the formula predicts jumps of $\\pm n/2$ in $N$ as the fermion mass crosses $\\pm e\\sqrt{2}v$, with $n$ zero modes in the intermediate mass range.","feed_headline":"Vortex fermion number follows a simple two-sign formula","feed_subtitle":"A boundary-independent result that predicts n zero modes in the intermediate mass window.","key_machinery":"The load-bearing object is the spectral asymmetry, or $\\eta$ invariant, $\\eta(s,H)=\\sum_{\\lambda>0}\\lambda^{-s}-\\sum_{\\lambda<0}(-\\lambda)^{-s}$, whose value at $s=0$ is related to fermion number by $N=-\\tfrac{1}{2}\\eta(0,H)$. Under a local variation of the background fields, the variation of $\\eta$ is given by a heat kernel coefficient on the boundary; on the disk this reduces to an integral of $\\delta A_\\parallel$ over the circle $S^1_R$. The edge states, found as exponentially decaying solutions of the near-boundary Dirac equation, combine into an effective one-dimensional boundary Hamiltonian $G(g)=i(\\partial_\\parallel-ig(A_\\parallel+n/R))$ whose $\\eta$ variation is computed from the $d=1$ heat kernel coefficient $a_0$. The difference between the disk result and the edge-state contribution cancels the boundary dependence and yields the mass-dependent sign structure.","core_discovery":"The central claim is that the spectral asymmetry of the Dirac Hamiltonian on the ANO vortex background is exactly captured by a mass-dependent combination of two sign functions. Placing the vortex on a large disk, the paper computes the variation of the $\\eta$ invariant through a boundary heat kernel coefficient, subtracts the contribution of the edge states, and takes the disk radius to infinity. The resulting $\\eta$ invariant is independent of the bag boundary parameter $\\varepsilon$, and the vacuum fermion number follows as $N(H) = [\\operatorname{sgn}(m+e\\sqrt{2}v)+\\operatorname{sgn}(m-e\\sqrt{2}v)]\\,n/4$. Consequently the fermion number takes the values $-n/2$, $0$, or $n/2$ depending on the mass interval, and the discontinuity at $|m|=e\\sqrt{2}v$ indicates $n$ fermionic zero modes on the plane for $-e\\sqrt{2}v<m<e\\sqrt{2}v$.","pith_inferences":["Beyond the paper: if boundary independence holds for a wider class of local boundary conditions, the disk-plus-edge-state scheme becomes a general algorithmic shortcut, computing $N$ from boundary heat kernel coefficients and one-dimensional edge spectra without a full spectral analysis on the plane.","Beyond the paper: the edge-state charges that vary continuously with $m$ suggest tunable fractionally charged edge channels; a lattice or tight-binding simulation of this Dirac–Higgs system could test the predicted edge spectrum and charges directly.","Beyond the paper: the zero-mode prediction for $-e\\sqrt{2}v<m<e\\sqrt{2}v$ could be checked by an explicit index-theoretic count or numerical spectral flow on the plane, providing a test of the method that does not rely on the subtraction formula."],"forward_implications":["The $\\eta$ invariant, and hence the vacuum fermion number, is independent of the bag boundary parameter $\\varepsilon$, providing a consistency check for the heat-kernel/edge-state method.","For $|m|<e\\sqrt{2}v$ the vacuum fermion number vanishes, and the spectral asymmetry jump implies $n$ fermionic zero modes on the plane.","For the vortex on a disk, the edge states carry charge $e/2$ in the intermediate mass range, and charges varying continuously between $0$ and $e$ outside that range.","In the $v\\to0$ limit the formula reduces to the known result $N=n\\,\\operatorname{sgn}(m)/2$ for a planar fermion in an external magnetic field."],"supporting_citations":[{"why":"Supplies the heat-kernel-plus-edge-state subtraction method that the paper tests on the ANO vortex.","marker":"[1]"},{"why":"Establishes the relation $N=-\\tfrac{1}{2}\\eta(0,H)$ between vacuum fermion number and spectral asymmetry used throughout.","marker":"[7]"},{"why":"Defines the Abrikosov vortex background whose quantum properties are studied.","marker":"[8]"},{"why":"Defines the Nielsen–Olesen vortex configuration and its winding number $n$.","marker":"[9]"},{"why":"Provides the variational formula for $\\eta$ in terms of a heat kernel coefficient, used to compute variations on the disk.","marker":"[24]"},{"why":"Supplies the heat kernel expansion coefficients $a_1$ and $a_2$ used for boundary integrals and the zero-mode index.","marker":"[28]"},{"why":"Provides the index theory argument used to show the existence of $n$ zero modes at $m=0$.","marker":"[30]"}],"fun_headline_variants":["Vortex fermion number: a two-sign formula","Two signs decide vortex fermion number","Vortex edge states carry fractional charge","Fermion number on vortex: -n/2, 0, or n/2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The subtraction formula (19) assumes that the fermion density near the disk boundary comes entirely from edge states, so that the plane's fermion number equals the disk's fermion number minus the edge-state contribution; the authors note this is natural to assume though hard to prove rigorously.","fun_headline_variants_meta":{"raw":{"variants":["Vortex fermion number: a two-sign formula","Two signs decide vortex fermion number","Vortex edge states carry fractional charge","Fermion number on vortex: -n/2, 0, or n/2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000556,"raw_usage":{"total_tokens":2613,"prompt_tokens":878,"completion_tokens":1735,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":1669}},"tokens_in":494,"tokens_out":1735,"duration_ms":14491,"temperature":1.0,"reasoning_tokens":1669,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:58:12.832961+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the spectral asymmetry directly, for instance by numerically diagonalizing the Dirac Hamiltonian on large disks over a range of masses, and examine whether the plane-limit $\\eta(0,H)$ is independent of the boundary parameter $\\varepsilon$ and jumps by $-n$ exactly at $m=\\pm e\\sqrt{2}v$. Any boundary-condition dependence in the large-$R$ limit, or a jump at a different mass, would refute Eq. (52).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Nielsen–Olesen vortex configuration and its winding number $n$."},{"cited_title":"Edge states and the $\\eta$ invariant","cited_arxiv_id":"2305.13606","evidence_quote":"Supplies the heat-kernel-plus-edge-state subtraction method that the paper tests on the ANO vortex."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Abrikosov vortex background whose quantum properties are studied."},{"cited_title":"Alvarez-Gaume, S","cited_arxiv_id":null,"evidence_quote":"Provides the variational formula for $\\eta$ in terms of a heat kernel coefficient, used to compute variations on the disk."},{"cited_title":"Fursaev and D","cited_arxiv_id":null,"evidence_quote":"Provides the index theory argument used to show the existence of $n$ zero modes at $m=0$."}],"review_version":1}