{"id":"07241152-dfa3-41b9-a74d-a30a417e106b","arxiv_id":"2505.24777","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A one-dimensional superconductor with two boundary magnetic impurities has a special 'supersymmetric' point where the nine degenerate low-energy boundary states form spl(2,1)⊗spl(2,1) representations and support zero energy modes.","lead":"Using an exact method called Bethe ansatz, the authors solve a one-dimensional superconducting wire whose two ends are coupled to magnetic impurities. They find a special coupling point where nine boundary states become degenerate and are connected by zero energy modes that realize a supersymmetry algebra.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on an unproven exhaustive classification of boundary Bethe strings; a second boundary string appears in App. B.4, so the nine-state count is not secured.","rationale":"The reader's weakest assumption—completeness of the Bethe ansatz and exhaustive string classification—is exactly the load-bearing point. The paper's Phase IV and Section V arguments are internally consistent within the assumed nine-dimensional subspace: the Clebsch-Gordan decomposition and the action of the Σ operators check out if and only if the nine listed states are the complete low-energy sector. The Appendix even supplies direct evidence that the string classification is fragile: a second short boundary string λ_bs'=±i(3/2−a_B) is introduced in App. B.4 without a general proof that no further short strings exist. No internal algebraic contradiction or sign error was found in the degeneracy calculation at a_A=a_B=1; the spin part of the bound-state energy E_spin,A=−Δ sin(a_A π) does vanish there, and the thermodynamic-limit caveat for the charging energy is stated honestly. Thus the appropriate verdict remains CONDITIONAL: the central claim is plausible and self-contained, but it is not secured until the completeness question is settled. A finite-N numerical solution of the Bethe equations, or a lattice DMRG/ED benchmark, would directly test the nine-state counting and would either remove or confirm the main uncertainty.","tokens_in":32231,"tokens_out":5899,"duration_ms":67984,"concrete_test":"Solve the full Bethe equations (19) numerically for finite N_e (e.g., N_e=20, 40, 80) at b>0 and a_A=a_B=1, allowing for arbitrary complex root configurations: scan for all solutions with distinct λ_α, discard λ=0, and classify by |Im λ| and energy from Eqs. (21) and (30). If any solution besides the all-real-root distribution, the single boundary strings λ=±i/2 at one boundary, and the two-boundary-string state lies below the gap Δ, the nine-state counting fails. A complementary DMRG or exact-diagonalization check on a lattice Gross-Neveu model with two boundary impurities would verify whether the low-energy sector at the supersymmetric point contains exactly nine states.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central result—nine degenerate boundary states at a_A=a_B=1 realizing [1/2,1/2]⊗[1/2,1/2]=[1,1]⊕[3/2,1/2] and the zero-energy operators (47)-(49)—depends on the assertion that the low-energy spectrum contains exactly the all-real-root state, the two states with one boundary string λ_bs=±i(a_m−1/2), and the state with two such strings. This classification is not proven. In App. B.4 the 'uniqueness' of λ_bs is stated 'by observation,' and immediately afterward a second boundary string λ_bs'=±i(3/2−a_B) is introduced to construct the S_z=0 unscreened states. The paper does not prove that these two exhaust the short-string sector, nor that additional complex root pairs in the bulk do not produce states within the gap. If even one extra low-energy state exists, the nine-dimensional Hilbert space and the spl(2,1)⊗spl(2,1) multiplet structure collapse, and the Σ operators are not zero modes of the full Hamiltonian. This is a standard but real gap in the coordinate Bethe ansatz for two boundaries: completeness is assumed, not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies a one-dimensional spin-singlet superconductor (Gross-Neveu model) coupled at both boundaries to spin-1/2 magnetic impurities via spin-exchange interactions. Using nested coordinate Bethe ansatz, the authors solve the model and obtain a boundary phase diagram with Kondo, Yu-Shiba-Rusinov (YSR), and unscreened phases at each edge. In the YSR-YSR phase the low-energy boundary Hilbert space contains nine states; at the point a_A=a_B=1, where the boundary string energy E_bound=-Delta sin(a pi) vanishes, these nine states become degenerate in the thermodynamic limit. The authors show that they realize the representation [1/2,1/2] tensor [1/2,1/2] = [1,1] direct sum [3/2,1/2] of spl(2,1) tensor spl(2,1), and construct operators Sigma^x_+, Sigma^x_-, and Sigma^z, Eqs. (47)-(49), that map between odd-parity states and are claimed to be exact zero-energy modes.","tokens_in":32476,"tokens_out":15175,"duration_ms":163550,"significance":"If the low-energy spectral classification is complete, this is a significant exact result: a two-boundary integrable model with dynamical boundary degrees of freedom exhibiting emergent boundary supersymmetry, an explicit zero-energy-mode construction, and a full phase diagram. The derivation of the boundary string energy, the integral-equation analysis, and the explicit Clebsch-Gordan decomposition are careful and internally consistent, and the paper is unusually self-contained. The central caveat is that the nine-state counting relies on an unproven completeness assumption for the boundary Bethe ansatz, so the result is conditional on closing that gap.","major_comments":[{"comment":"The nine-state counting at the YSR-YSR supersymmetric point, and hence the irreducible decomposition in Eq. (44) and the zero-energy operators in Eqs. (47)-(52), rests on the assumption that the low-energy spectrum contains exactly the states built from real Bethe roots plus the two boundary strings lambda_bs = +/-i(a_m-1/2) and lambda_bs' = +/-i(3/2-a_m). In Appendix B.2 the uniqueness of lambda_bs is asserted 'by observation', and in Appendix B.4 a second boundary string is introduced without a systematic derivation. No completeness proof for the two-boundary coordinate Bethe ansatz is provided, and it is not shown that other complex-root configurations (bulk strings, additional boundary strings) cannot produce states inside the superconducting gap. Since an extra low-energy state would change the dimension of the low-energy Hilbert space and destroy the claimed multiplet structure, this gap is load-bearing. I request a proof of completeness of the Bethe ansatz solution, or at least a systematic classification of all solutions to Eq. (19) in the thermodynamic limit, supplemented where possible by independent numerical evidence (for example DMRG) that the low-energy spectrum at a_A=a_B=1 is exactly nine-dimensional.","section":"Appendix B.4, Eq. (19)"},{"comment":"The paper calls the operators (47)-(52) 'exact zero energy modes', but their action is demonstrated only on the four odd-parity states, and the degeneracy of those states at the supersymmetric point is stated to hold only up to exponential accuracy e^{-L} in the system size. At finite L, the full Hamiltonian therefore does not strictly commute with these operators on the ground subspace, and no computation of [H,Sigma] is given. Please clarify in what sense the modes are 'exact' -- for example, by specifying that they are zero modes of the effective low-energy Hamiltonian in the thermodynamic limit -- or provide a proof that the commutator vanishes exactly in that limit.","section":"Secs. IV.B.1 and VI.A, Eqs. (47)-(52)"}],"minor_comments":[{"comment":"The definition of c_m in Eq. (A7) appears inverted relative to Eq. (13); it should read c_m = 2J_m/(1-3J_m^2/4), and the denominator '2J' should be '2J_m'.","section":"Eq. (A7)"},{"comment":"The caption contains the typo 'irreducbible'; it should be 'irreducible'.","section":"Fig. 2 caption"},{"comment":"The phase heading 'YSR-Kondo (YSR-K)' is inconsistent with the main-text notation 'Kondo-YSR (K-YSR)'; please unify the terminology.","section":"Appendix B.2 heading"},{"comment":"The phrase 'exponentially degenerated' should be 'exponentially degenerate'.","section":"Sec. IV.B.1"},{"comment":"Reference [8] duplicates reference [4], and the author list of [12] appears corrupted ('Glazman, and Glazman').","section":"References"},{"comment":"The operator B_A in Eq. (49) is not explicitly defined; please state that it is the baryon number operator at edge A.","section":"Eq. (49)"},{"comment":"The algebra name is typeset inconsistently ('spl(2,1)' in the text and 'SPL(2,1)' in Appendix C); please unify the notation.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the completeness of the boundary Bethe ansatz classification of low-energy states. If the authors can supply a completeness proof or independent numerical confirmation of the nine-dimensional low-energy subspace at the supersymmetric point, the paper would be a strong contribution. The second major comment is more about precision of the 'exact zero mode' claim than about the correctness of the underlying physics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. It solves the two-impurity Gross-Neveu superconductor with coordinate Bethe ansatz and finds a point where nine boundary states become degenerate and assemble into spl(2,1)⊗spl(2,1). That is new relative to the single-impurity literature, and the derivation of the boundary string energy E_bound=-Δ sin(aπ) is clean enough to be checked by hand. The supersymmetric point is not fitted; it follows from that formula. The Clebsch-Gordan decomposition is worked out explicitly. Genuine substance.\n\nThe soft spot is the one the stress-test note hits: the nine-state count is an assertion, not a theorem. The paper states by observation that the short boundary string λ_bs=±i(a_m-1/2) is unique, then introduces a second short string λ_bs'=±i(3/2-a_B) in App. B.4. No proof is given that these two exhaust the short-string sector, and no argument rules out complex root pairs in the bulk producing low-energy states. This is the standard completeness gap in coordinate BA, and it is real. If an extra low-energy state exists, the nine-dimensional Hilbert space and the supersymmetric degeneracy would not be complete. I don't think it is fatal—the same assumption underlies the accepted single-impurity solutions, and the phase picture is physically sensible—but it is load-bearing enough that a referee should push for a proof or at least a much more careful classification of boundary strings.\n\nMinor: the phrase \"exact zero energy modes\" overstates the finite-size situation. In a finite wire the odd-parity states are split by e^{-L} and the even-parity ones by 1/L; the paper acknowledges this but the abstract and Section VI say exact. In the thermodynamic limit it is fine.\n\nWho should read this: people who work on integrable boundary problems, YSR physics, or emergent supersymmetry in 1D. The paper deserves serious refereeing, with a referee who knows Bethe ansatz completeness issues. I'd accept it with the expectation that the string classification be tightened, perhaps in an appendix.","headline":"Read this one: a real two-impurity Bethe ansatz extension with a clean derivation of the supersymmetric point, but the nine-state counting rests on an unproven boundary-string classification that a referee should push to tighten.","tokens_in":33027,"tokens_out":2839,"would_cite":true,"duration_ms":32001,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B23","81R12","17B81"],"pacs":["74.20.-z","75.20.Hr"],"model":"deepseek-v4-flash","headline":"Coupling a one-dimensional spin-singlet superconductor to spin-1/2 impurities at both ends produces a nine-fold degenerate boundary spectrum at a special coupling point, organized by the supersymmetric algebra spl(2,1)⊗spl(2,1) and giving…","keywords":["Bethe ansatz","Gross-Neveu model","supersymmetry","spl(2,1) algebra","Yu-Shiba-Rusinov states","Kondo effect","zero-energy modes","one-dimensional superconductor"],"falsifier":"Solve the two-boundary Bethe equations (19) numerically in the regime 1/2 < a < 3/2 and enumerate all solutions with |Im λ| < 1; finding any short string other than ±i(a - 1/2) would add states beyond the nine and break the claimed degeneracy. Alternatively, diagonalize a lattice regularization of the model (for example with DMRG) at a_A = a_B = 1 and check whether exactly nine states become degenerate in the thermodynamic limit and whether the operators Σ^x_+, Σ^x_-, Σ^z of Eqs. (47)-(49) act as zero-energy modes with the stated action on those states.","tokens_in":31996,"feed_emoji":"⚛️","tokens_out":8506,"duration_ms":85098,"temperature":0.7,"pith_summary":"This paper claims that when a one-dimensional spin-singlet superconductor, described by the Gross-Neveu model, is coupled to spin-1/2 magnetic impurities at both ends via spin-exchange, the low-energy boundary degrees of freedom can become exactly supersymmetric. Solving the model with nested coordinate Bethe ansatz for arbitrary bulk and impurity couplings, the authors find three possible impurity phases: Kondo, Yu-Shiba-Rusinov (YSR), and unscreened. When both impurities are in the YSR phase, the boundary Hilbert space holds nine states; at the point where both RG-invariant parameters equal a_A=a_B=1, all nine become degenerate in the thermodynamic limit and organize into the representation [1/2,1/2]⊗[1/2,1/2] = [1,1]⊕[3/2,1/2] of the superalgebra spl(2,1)⊗spl(2,1). At that point the boundary bound-state energy vanishes and the system possesses exact zero-energy modes in the odd fermion parity sector, written explicitly in terms of the algebra's generators. If correct, this gives a concrete, exactly solvable example of supersymmetry emerging in a boundary Hilbert space, with zero-energy modes that survive in the thermodynamic limit and connect degenerate ground states.","feed_headline":"Exact zero modes emerge from boundary supersymmetry","feed_subtitle":"Nine boundary states of a 1D superconductor with two magnetic impurities become degenerate at one coupling point.","key_machinery":"The argument runs through the nested coordinate Bethe ansatz for the Gross-Neveu model with two boundary impurities. The single-particle momenta are quantized by boundary S-matrices satisfying reflection equations; the spin sector is governed by Bethe equations whose RG-invariant boundary parameters d_m (real or imaginary, d_m = i a_m) encode the impurity couplings. The low-energy states correspond to real Bethe roots plus a unique short boundary string λ_bs = ±i(a_m - 1/2) per edge, which represents a bound state of energy -Δ sin(a_m π). The screened and unscreened impurity configurations form the three-state representation [1/2,1/2] of the spl(2,1) superalgebra at each edge, and their tensor product splits as [1,1]⊕[3/2,1/2]; at a_A=a_B=1 the generators V±, W±, Q±, Q3, B of each edge combine into the exact zero-energy operators Σ^x_+, Σ^x_-, Σ^z.","core_discovery":"The central claim is that at the supersymmetric point a_A=a_B=1, the boundary bound state at each edge has spin energy E_spin,m = -Δ sin(a_m π), which vanishes exactly; consequently the nine low-energy states tabulated in Table I are degenerate in the thermodynamic limit. These nine states decompose as [1/2,1/2]⊗[1/2,1/2] = [1,1]⊕[3/2,1/2] under spl(2,1)⊗spl(2,1), and the operators Σ^x_+, Σ^x_-, and Σ^z in Eqs. (47)-(49) are exact zero-energy modes in the odd fermion parity sector, mapping the states of the [1,1] representation onto those of [3/2,1/2] and vice versa.","pith_inferences":["If the nine-fold degeneracy is robust to perturbations that preserve the spl(2,1) algebra, the boundary Hilbert space at the supersymmetric point could serve as a protected degenerate subspace, potentially useful for quantum information; the paper does not claim this.","The same algebra structure may appear in other exactly solvable models with dynamical boundaries, such as coupled Kondo impurities in gapped spin chains; checking for an analogous supersymmetric point in those models would test the generality of the mechanism.","A direct numerical check of Bethe-string completeness—enumerating all low-lying solutions of the two-boundary Bethe equations—would settle whether the supersymmetric degeneracy at a=1 is exact or an artifact of the assumed root configuration; the paper leaves this open.","The open question the authors raise, whether the nine degenerate ground states carry local fractional boundary quantum numbers, could be addressed by computing spin profiles and entanglement spectra in a lattice realization of the model."],"forward_implications":["At a_A=a_B=1, the ground-state manifold of the YSR-YSR phase is nine-fold degenerate, and the supersymmetry of the boundary degrees of freedom is restored in the thermodynamic limit.","The boundary bound-state energy E_spin = -Δ sin(a π) changes sign across a=1, so each impurity undergoes a first-order quantum phase transition between screened and unscreened ground states.","The three operators Σ^x_+, Σ^x_-, Σ^z are exact zero-energy modes in the odd fermion parity sector at the supersymmetric point.","Exact zero-energy modes require two boundaries: with a single impurity no exact ZEM exists, because removing the boundary bound state costs a charging energy of order 1/L."],"supporting_citations":[{"why":"Supplies the bulk Gross-Neveu Bethe ansatz solution used as the starting point for the two-boundary construction.","marker":"[38]"},{"why":"Provides the boundary algebraic Bethe ansatz transfer-matrix method used to diagonalize the two-boundary problem.","marker":"[42]"},{"why":"Gives the reflection equations that ensure the boundary S-matrices are consistent.","marker":"[43]"},{"why":"Supplies the off-diagonal Bethe ansatz framework and the exclusion of λ=0 roots used in the Bethe equations.","marker":"[44]"},{"why":"Establishes the single-boundary solution and the RG-invariant parameter d_m that the two-boundary phase diagram generalizes.","marker":"[41]"},{"why":"Identifies the YSR boundary string and its energy expression -Δ sin(aπ), which is central to the zero-energy condition at a=1.","marker":"[36]"},{"why":"Supplies the spl(2,1) representation theory and baryon-number notation used to label the boundary multiplets.","marker":"[46]"},{"why":"Supports the interpretation of boundary strings as poles of the dressed boundary S-matrix, the origin of the λ_bs solution.","marker":"[52]"}],"fun_headline_variants":["Boundary supersymmetry yields exact zero modes","Zero modes emerge from boundary SUSY in 1D superconductor","Exact zero modes at supersymmetric boundary coupling point","Degenerate nine states signal boundary supersymmetry","Boundary SUSY exact zero modes in one-dimensional superconductor"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the Bethe ansatz is complete: that only real Bethe roots plus the short boundary string λ_bs = ±i(a_m - 1/2) contribute to the low-energy spectrum, and no other string solutions or missing roots add states; the uniqueness of that boundary string is asserted by observation (Appendix B.4) rather than proven, and no completeness proof for the two-boundary Bethe ansatz is given.","fun_headline_variants_meta":{"raw":{"variants":["Boundary supersymmetry yields exact zero modes","Zero modes emerge from boundary SUSY in 1D superconductor","Exact zero modes at supersymmetric boundary coupling point","Degenerate nine states signal boundary supersymmetry","Boundary SUSY exact zero modes in one-dimensional superconductor"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000159,"raw_usage":{"total_tokens":1225,"prompt_tokens":936,"completion_tokens":289,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":208}},"tokens_in":552,"tokens_out":289,"duration_ms":3269,"temperature":1.0,"reasoning_tokens":208,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:13:54.177535+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the two-boundary Bethe equations (19) numerically in the regime 1/2 < a < 3/2 and enumerate all solutions with |Im λ| < 1; finding any short string other than ±i(a - 1/2) would add states beyond the nine and break the claimed degeneracy. Alternatively, diagonalize a lattice regularization of the model (for example with DMRG) at a_A = a_B = 1 and check whether exactly nine states become degenerate in the thermodynamic limit and whether the operators Σ^x_+, Σ^x_-, Σ^z of Eqs. (47)-(49) act as zero-energy modes with the stated action on those states.","supporting_citations":[{"cited_title":"We do naturally expect that the impurity is partially screened, but eventually be- comes unscreened deep in the unscreened phasea≫1","cited_arxiv_id":null,"evidence_quote":"Supplies the bulk Gross-Neveu Bethe ansatz solution used as the starting point for the two-boundary construction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the boundary algebraic Bethe ansatz transfer-matrix method used to diagonalize the two-boundary problem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the reflection equations that ensure the boundary S-matrices are consistent."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the YSR boundary string and its energy expression -Δ sin(aπ), which is central to the zero-energy condition at a=1."},{"cited_title":"Wang, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the spl(2,1) representation theory and baryon-number notation used to label the boundary multiplets."}],"review_version":1}