{"id":"65e9b4f5-9560-4524-bfe3-b1de0a2f3742","arxiv_id":"2608.09749","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"At Rényi index 1/2, the boundary stabilizer-Rényi response across BDI transitions equals |Δω| ln 2 for the studied free-fermion families.","lead":"A free-fermion chain with open boundaries shows a quantized jump, |Δω| ln 2, in the bulk-subtracted stabilizer Rényi entropy when the topological winding number changes. The result gives a new, wavefunction-based way to count Majorana boundary channels across SPT transitions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Primitive-family quantized endpoint rests on signed-Pfaffian fits; unresolved sign-violating minors are controlled only at the 1e-3–1e-2 level, so the claimed 1e-5 agreement with ln 2 is not yet established for the exact absolute-minor observable.","rationale":"The reader's weakest assumption is exactly the load-bearing concern: for the primitive range-two family, the signed Pfaffian is not an exact absolute-minor sum (Eq. S2.39), and the independent positive-weight control of the exact observable is at the 1e-3–1e-2 level, not the 1e-5 level claimed for the Pfaffian-assisted extrapolation. This is the single point on which the central universal claim depends: the exact multiplication law (Eq. 12) covers only the shifted-decimated family, and the primitive deformation is what demonstrates robustness beyond decimation, self-duality, and special termination. If hidden sign-violating sectors contribute an O(1) boundary difference, the primitive endpoint could shift away from ln 2. The proposed test—extending the positive doubled-Slater bridge to larger L and smaller m—directly checks whether the unresolved sectors grow or remain negligible. The paper's exact identities, ScS.2.2 and Eq. (12), are solid, and the numerical analysis is extensive and honestly bounded; no fraud or inconsistency is alleged. The conditional verdict is therefore appropriate, and no verdict change is needed.","tokens_in":30709,"tokens_out":4536,"duration_ms":44450,"concrete_test":"Evaluate the exact absolute-minor response Δabs_L(m) = ln S_top - ln S_triv for the primitive family λ=0.2 using the positive doubled-Slater bridge at larger sizes and smaller masses than Table S6.1: specifically L = 240, 480, 960 at m = 0.1, 0.05, and 0.025, with enough Monte Carlo samples to reduce between-run uncertainty below 5e-4. Compare Δabs_L(m) with the signed-Pfaffian ΔPf_L(m). If the difference remains below ~1e-3 and shows no growth with L, the hidden-sign-sector concern is mitigated. If the difference grows with L or exceeds ~1e-2, the fitted primitive endpoint is unreliable and the quantized-value claim for this family must be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The exact shifted-decimated family (Eq. 12) is proved, but the universal one-channel constant ln 2 is supported only numerically for the primitive range-two deformation fλ(z;h) = h - z^{-1} + λ z^{-2}. For this family the fixed-selector Pfaffian is not an exact absolute-minor sum: Eq. S2.39 gives an explicit sign-violating minor at L=7. The boundary-pattern decomposition in Sec. S6.2 resolves detected negative sectors only at the 1e-23 level and is explicitly not an upper bound; the positive doubled-Slater bridge in Sec. S6.3 controls the complete absolute-minor difference at representative points (m=0.1, 0.05; L=40–120) only at the 1e-3 to 1e-2 level (Table S6.1). The global Pfaffian fit in Table S6.2 quotes endpoints within 7.4e-6 of ln 2, but the exact observable it approximates is not shown to agree to that precision. If unresolved sign-violating sectors contributed an O(1) boundary term that does not vanish in the massive limit, the true RSPT_{1/2} for the primitive family could differ from ln 2 by order 1e-3 or more, undermining the strong quantization claim. The paper is transparent about this gap, and the qualitative conclusion (boundary response counts |Δω| channels) is plausible, but the quantitative equality to ln 2 for the non-decimated family is conditional on an unproven assumption about hidden sign sectors.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the stabilizer Rényi entropy at Rényi index 1/2 for finite open BDI chains, aiming to isolate a bulk-subtracted boundary response across a mass inversion. The central claim is that this response approaches |Δω| ln 2, where Δω is the change in winding number, so that the boundary term counts the Majorana channels that change the fermionic SPT index. The paper proves an exact finite-size multiplication law for the shifted-decimated family f_{r,d}(z;h)=z^r(h+z^d), showing that the full crossover equals d times the TFI crossover. For a primitive range-two deformation f_λ(z;h)=h−z^{-1}+λz^{-2}, the paper uses Pfaffian evaluations, boundary-pattern decompositions, and a positive doubled-Slater bridge to argue that the one-channel endpoint is ln 2 within numerical accuracy. Extensive supplemental material documents the technical definitions, the limitations of the signed-Pfaffian representation, and the systematic error estimates.","tokens_in":31038,"tokens_out":3308,"duration_ms":26747,"significance":"If the central claim is correct, the paper identifies a genuinely new phenomenon: the O(1) boundary term of the stabilizer-Rényi entropy, not the volume law, carries the change in the fermionic SPT index. The exact finite-size identity for the shifted-decimated family, Eq. (12), is a strong, proved structural result that cleanly separates the background boundary index from the topology-changing channels. The paper is also unusually transparent about its numerical uncertainties: it explicitly states that the boundary-pattern decomposition is not an all-sector upper bound, that the signed-Pfaffian free energy is not the exact absolute-minor sum for the primitive family, and that the independent positive-weight control of the exact observable is at the 10^{-3}–10^{-2} level. These honest limitations are valuable and make the manuscript's claims easy to assess, even where the evidence is incomplete.","major_comments":[{"comment":"The quantitative claim that the one-channel endpoint equals ln 2 to within 7.4×10^{-6} (Table S6.2) applies to the signed-Pfaffian free energy Φ_L, not to the exact absolute-minor sum S_L. Equations (S2.43)–(S2.44) show that the difference is δ_L = ln(1 + 2W_-(L)/A_L), and the paper's own boundary-pattern decomposition gives only detected lower bounds on W_-/A_L, with no all-sector upper bound. The independent positive bridge in Table S6.1 gives offsets |Δ^{abs}_L − Δ^{Pf}_L| up to 3.6×10^{-3} at the representative points, and the paper conservatively assigns a 10^{-2} mixing scale. Thus the exact absolute-minor observable is presently consistent with ln 2 only at the 10^{-3}–10^{-2} level, not at the 10^{-5} level claimed for the Pfaffian-assisted fits. The main text should either provide a rigorous upper bound on the sign-violating sectors, or explicitly state that the primitive-family quantization to ln 2 is established only at the 10^{-3}–10^{-2} level while the 10^{-5} endpoint refers to the signed-Pfaffian surrogate.","section":"S6.2–S6.3, Table S6.1"},{"comment":"The exact multiplication law reduces the universal one-channel constant to the TFI limit, but it does not prove that the TFI limit is ln 2. Equation (S4.15) is stated as a numerically supported asymptotic statement, and the supporting data in Sec. S5.3 show convergence only to numerical accuracy. Since the central claim is the quantized value |Δω| ln 2, the manuscript should more carefully separate the proved structural part (multiplication by |Δω|) from the numerically inferred part (the one-channel constant). In particular, the phrase 'exactly in its multiplication by |Δω| and numerically in the primitive one-channel constant' in the main text is honest, but the abstract and introduction state the full equality as the headline result. A revision should make the logical status of ln 2 explicit in the abstract and introduction, or provide an independent derivation of the one-channel constant.","section":"S4.3, Eqs. (S4.14)–(S4.16)"}],"minor_comments":[{"comment":"The notation M_{1/2} = 2 ln S_L − L ln 2 is clear, but S_L is first defined only implicitly as the absolute-minor sum; a one-sentence reminder in the main text that S_L ≡ D_1(G) would help readers who do not consult the supplement.","section":"Main text, Eq. (2)"},{"comment":"The symbol R^{Pf}_{1/2} appears in the supplemental discussion (Eq. (S6.32)) without being defined in the main text; please define it where the fitted endpoint is first mentioned.","section":"Main text, near Eq. (10) and Sec. S6.6"},{"comment":"The global fit uses 23 parameters for 45 raw values; the near-critical ansatz terms m ln m, m, m^2 ln m, m^2 are introduced ad hoc. The holdout tests are useful, but a brief justification of why these are the only relevant near-critical terms would strengthen the extrapolation.","section":"Sec. S6.6, Eq. (S6.31)"},{"comment":"In Table S6.1, the column header Δ^{abs}_L − Δ^{Pf}_L is clear, but the text says the offsets 'should not be interpreted as pointwise compatible with zero' while also noting they do not establish growth with L; please make the statistical interpretation more precise, including whether the between-run dispersions are one sigma or two sigma.","section":"Sec. S6.3"},{"comment":"The discussion of the Rényi-index dependence is appropriately cautious, but the sentence 'no value is assigned here at α = 4, and no interpolation or locked-side branch is assumed' could be moved earlier in the section to avoid any impression that the paper claims a full Rényi-index interpolation.","section":"Sec. S7.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is transparent and technically careful, and the exact multiplication result is a solid contribution. The main gap is the mismatch between the 10^{-5} precision claimed for the signed-Pfaffian fits and the 10^{-3}–10^{-2} precision of the exact absolute-minor check for the primitive family. This gap directly affects the headline quantitative claim, so a major revision is appropriate. The paper might be acceptable after either (a) strengthening the sign-sector control to a genuine upper bound, or (b) explicitly reframing the universal constant as numerically supported at the 10^{-3} level while keeping the exact multiplication law as the proved statement. I do not see any evidence of circularity or fabrication; the limitations are stated openly in the supplement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The exact result here is real and worth knowing: for the shifted-decimated family f(z)=z^r(h+z^d), the full boundary-response crossover is exactly d times the TFI crossover, at every finite compatible size. That is Eq. (12), and it is proved cleanly. It gives an exact sense in which the response counts the number of Majorana channels that cross the unit circle, independent of the background winding. That part is solid and new.\n\nThe paper also does something rare: it is unusually honest about the limitations of its own numerics. The signed-Pfaffian reduction is not exact for the primitive range-two family; the authors exhibit an explicit sign-violating minor, resolve boundary sectors, and then check the exact absolute-minor difference with an independent positive-weight bridge. That is the right way to handle a broken selector, and the transparency earns real credit.\n\nThe soft spot is exactly the one the stress test identifies. The universal constant |Δω| ln 2 for the primitive non-decimated family rests on a global fit with 23 parameters whose signed-Pfaffian endpoints are within 7e-6 of ln 2. But the exact absolute-minor observable is only shown to agree with the Pfaffian at the 1e-3 to 1e-2 level, at a few representative points, with no proven upper bound on the unresolved sign-violating sectors. So the strong claim—that the primitive one-channel endpoint is ln 2 to five or six digits—is not established. What is established is a robust qualitative conclusion: the boundary response tracks |Δω|, and the endpoint is consistent with ln 2 within a few times 10^-3. That is still a meaningful result, but the paper should not present the 7e-6 agreement as the precision of the physical observable.\n\nNo code or data are shipped, which is a minor complaint given the detail in the supplemental material, but it would help future verification.\n\nWho is this for? People working on stabilizer Rényi entropies, participation entropies, and fermionic SPT diagnostics. They will get the exact multiplication law and a careful numerical study of a plausible universality. The paper deserves a serious referee: the main claim is interesting, the exact part is correct, and the conditional part is clearly flagged. I would not desk-reject it. I would ask the author to either close the gap between the Pfaffian extrapolation and the exact absolute-minor control, or reframe the primitive-family claim at the 1e-3 level of precision actually demonstrated.","headline":"A genuinely new exact finite-size identity plus a plausible but numerically conditional quantization claim; the paper deserves review but the ln 2 for the primitive family is not yet proven at the claimed precision.","tokens_in":31630,"tokens_out":1070,"would_cite":true,"duration_ms":10421,"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 bulk-subtracted boundary term of the stabilizer-Rényi entropy of open free-fermion BDI chains converges to |Δω| ln 2 at Rényi index 1/2, meaning the O(1) boundary response, not the volume law, carries the fermionic SPT index change.","keywords":["stabilizer Rényi entropy","fermionic SPT phases","BDI chain","Majorana zero modes","winding number","Pfaffian reduction","Shannon-Rényi entropy","boundary response"],"falsifier":"Evaluate the exact absolute-minor difference for the primitive family at λ = 0.2 using the positive doubled-Slater bridge at larger sizes than the reported L = 120 and at smaller masses, and check whether Δabs_L − ΔPf_L grows with L beyond the ~3.6×10^−3 level already seen; a growth law in L, or a boundary refinement at window w ≥ 10 revealing negative sectors contributing at order 10^−10 or larger, would break the claimed quantization. Alternatively, build any finite-range BDI deformation with |Δω| = 1 that lies outside the coherent Pfaffian class and fit the same nine-mass global endpoint: a departure from ln 2 exceeding the stated systematic envelope would falsify the universality of the response.","tokens_in":30430,"feed_emoji":"⚛️","tokens_out":10006,"duration_ms":72681,"temperature":0.7,"pith_summary":"This paper claims that a specific boundary term of the stabilizer Rényi entropy — a measure of non-stabilizerness built from the probabilities of Pauli measurements — is quantized across transitions between fermionic symmetry-protected topological (SPT) phases. For open free-fermion BDI chains in one dimension, the paper compares two gapped phases on opposite sides of a mass inversion, subtracts the nonuniversal bulk contribution independently on each side, and finds that the remaining boundary response approaches |Δω| ln 2 at Rényi index 1/2, where |Δω| is the number of Majorana channels transferred (the change in the winding number). If correct, the O(1) boundary term, rather than the dominant volume law, carries the change in the fermionic SPT index, offering a way to read a topological invariant from the Pauli spectrum of the wave function. Exact finite-chain relations for a shifted-decimated family prove the multiplication of the full crossover by d = |Δω|, and a non-solvable deformation with broken bulk duality and local boundary perturbations is fitted to the same one-channel value ln 2 to within 7.4×10^−6. Through an exact correspondence, the same quantization implies a 2|Δω| ln 2 boundary difference in the computational-basis Shannon–Rényi entropy of the associated SSH-type chain.","feed_headline":"At one Rényi index, an entropy boundary term is quantized: |Δω| ln 2","feed_subtitle":"A bulk-subtracted stabilizer-Rényi difference approaches ln 2 per transferred Majorana channel.","key_machinery":"The central object is the α = 1/2 stabilizer-Rényi entropy of a real Gaussian (free-fermion) state, which equals the absolute-minor sum S_L(G) = Σ |det G[I,J]| over the polar factor G of the finite open-chain matrix, with M_{1/2} = 2 ln S_L − L ln 2. In a coherent-sign chamber, this exponentially large sum collapses to a single Pfaffian, S_L(G) = |Pf[R(G) + J]|, allowing polynomial evaluation for large L — the main computational engine. The argument is carried by two complementary families: the shifted-decimated family z^r(h + z^d), whose finite open chain reduces exactly to d independent transverse-field Ising blocks plus r zero-mode blocks, giving the exact identity Δ_{r,d}(x) = d·Δ_TFI(x); and the primitive deformation f_λ(z;h) = h − $z^{{−1}}$ + λ $z^{{−2}}$, chosen because its signed Pfaffian demonstrably fails to be an exact absolute-minor sum, thereby breaking decimation, duality, and sign coherence. Around these sit the subtraction protocol — fitting independent bulk densities per side before extracting O(1) constants — and the positive doubled-Slater bridge, a sign-free Monte Carlo estimator of the exact absolute-minor difference used to rule out an order-unity selector artifact.","core_discovery":"The central claim is that, for the free-fermion BDI families studied, the O(1) boundary contribution to the α = 1/2 stabilizer-Rényi entropy, after independent removal of the bulk free-energy density on each side of a mass inversion, is a quantized response equal to R^SPT_{1/2} = |Δω| ln 2, with |Δω| the change in the winding number across the transition. Each Majorana channel that crosses the unit circle in the Laurent symbol leaves an imprint of ln 2 on the reduced boundary response, while modes common to both phases contribute nothing. The claim is established exactly for the shifted-decimated family z^r(h + z^d), where the whole finite-size crossover equals d times the transverse-field Ising crossover, and numerically for a primitive range-two deformation h − $z^{{−1}}$ + λ $z^{{−2}}$ in which the Pfaffian sign structure fails, bulk duality is broken, and local termination changes are applied; there the extrapolated endpoint lies within 7.4×10^−6 of ln 2, with a systematic envelope of 2×10^−5. The paper also argues that the response is not a property of the unprocessed total entropy: individual boundary constants are termination-dependent, and only the bulk-subtracted, same-termination, two-sided difference is universal.","pith_inferences":["If the quantization survives beyond free fermions, the α = 1/2 stabilizer-Rényi entropy could serve as a practical scalar diagnostic of fermionic SPT order accessible to Pauli-measurement protocols, since the entropy is assembled directly from Pauli expectation values.","The ln 2 value admits a boundary-state-counting reading — one transferred Majorana per end forms a nonlocal complex fermion with a two-dimensional occupation space — which the paper explicitly declines to identify with Affleck–Ludwig g-factors; connecting the two is a natural next step.","A testable extension is to map the full Rényi-index function R_α of the strict nested limit, which must interpolate between 2|Δω| ln 2 at α = 1/2 and 0 at α = ∞ and may develop structure at the boundary Rényi transition α = 4; the paper leaves this function open.","Because the generic anisotropic XY chain shows no resolved nonlinear boundary mismatch relative to TFI in the accessible fixed-x window, an analogous nested-limit analysis for the XY family — currently out of reach of exact enumeration — is a concrete place to test whether the quantization extends beyond the Pfaffian-solvable classes."],"forward_implications":["At Rényi index 1/2, the independently bulk-subtracted boundary response of the stabilizer-Rényi entropy across a BDI mass inversion equals |Δω| ln 2, so the change in the fermionic SPT index is carried by the O(1) boundary term rather than by the extensive volume law.","In the shifted-decimated family z^r(h + z^d), the entire finite-size crossover — not just its asymptotic value — is exactly d = |Δω| times the single-channel transverse-field Ising crossover, and is independent of the background winding r common to both phases.","Through the stabilizer–Shannon correspondence, the same quantization yields a computational-basis Shannon–Rényi boundary difference of 2|Δω| ln 2 in the half-filled SSH-type partner chain.","The exact min-entropy constraint M_∞ = L ln 2 forces the two-sided boundary difference to vanish at Rényi index infinity, so the response is genuinely index-dependent rather than a fixed constant.","The result is compatible with the reported vanishing topological-magic response in the Ising phases, because the present response probes the fermionic BDI boundary index carried by Jordan–Wigner Majoranas rather than nonlocally distributed bosonic SPT resources."],"supporting_citations":[{"why":"Defines the stabilizer Rényi entropy, the object whose boundary response is studied.","marker":"[20]"},{"why":"Establishes the stabilizer–Shannon correspondence and the absolute-minor/determinant-sum representation used at α = 1/2.","marker":"[46]"},{"why":"Supplies the minor-summation Pfaffian machinery for finite open chains that the paper extends to both physical ends away from criticality.","marker":"[60]"},{"why":"Identifies the open SSH chain as the number-conserving partner whose Shannon–Rényi boundary difference inherits the quantization.","marker":"[48]"},{"why":"Provides the boundary-CFT formulation of critical SREs whose open-boundary logarithm the exact family reproduces.","marker":"[44]"},{"why":"Supplies the open-boundary logarithmic coefficient (c/4 ln L) that the paper matches and extends to the massive side.","marker":"[45]"},{"why":"Provides the Pfaffian minor-summation identity underlying the single-Pfaffian reduction in Proposition S2.2.","marker":"[62]"},{"why":"Documents the vanishing topological-magic response in Ising phases, the contrasting construction the paper distinguishes from its own.","marker":"[53]"}],"fun_headline_variants":["Majorana channels imprint ln 2 on boundary stabilizer-Rényi entropy","Quantized boundary entropy: each changing Majorana channel adds ln 2","SPT boundary response: |Δω| ln 2 at α=1/2","Stabilizer-Rényi boundary term counts Majorana channels via |Δω| ln 2","Boundary entropy quantized: ln 2 per transferred Majorana channel"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that, in the non-solvable test family, the finitely many minors whose signs the Pfaffian formula gets wrong never accumulate enough weight to change the answer at large system sizes; the boundary-pattern analysis only bounds what it can detect, and the independent sign-free check has limited precision, so if hidden sign sectors contributed at order one, the fitted endpoint could move away from ln 2.","fun_headline_variants_meta":{"raw":{"variants":["Majorana channels imprint ln 2 on boundary stabilizer-Rényi entropy","Quantized boundary entropy: each changing Majorana channel adds ln 2","SPT boundary response: |Δω| ln 2 at α=1/2","Stabilizer-Rényi boundary term counts Majorana channels via |Δω| ln 2","Boundary entropy quantized: ln 2 per transferred Majorana channel"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000409,"raw_usage":{"total_tokens":2133,"prompt_tokens":969,"completion_tokens":1164,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":585,"completion_tokens_details":{"reasoning_tokens":1057}},"tokens_in":585,"tokens_out":1164,"duration_ms":5929,"temperature":1.0,"reasoning_tokens":1057,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:30:05.885457+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the exact absolute-minor difference for the primitive family at λ = 0.2 using the positive doubled-Slater bridge at larger sizes than the reported L = 120 and at smaller masses, and check whether Δabs_L − ΔPf_L grows with L beyond the ~3.6×10^−3 level already seen; a growth law in L, or a boundary refinement at window w ≥ 10 revealing negative sectors contributing at order 10^−10 or larger, would break the claimed quantization. Alternatively, build any finite-range BDI deformation with |Δω| = 1 that lies outside the coherent Pfaffian class and fit the same nine-mass global endpoint: a departure from ln 2 exceeding the stated systematic envelope would falsify the universality of the response.","supporting_citations":[{"cited_title":"Leone, S","cited_arxiv_id":null,"evidence_quote":"Defines the stabilizer Rényi entropy, the object whose boundary response is studied."},{"cited_title":"Hoshino, M","cited_arxiv_id":null,"evidence_quote":"Provides the boundary-CFT formulation of critical SREs whose open-boundary logarithm the exact family reproduces."},{"cited_title":"Quantized Stabilizer-Rényi Boundary Response across Fermionic SPT Transitions","cited_arxiv_id":null,"evidence_quote":"Provides the Pfaffian minor-summation identity underlying the single-Pfaffian reduction in Proposition S2.2."}],"review_version":1}