{"id":"a29da349-a8ed-4c3e-ad33-e05221823f5b","arxiv_id":"2508.19330","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A Bethe-ansatz analysis yields the impurity entropy across the four phases and predicts entropy overshoots above ln 2 when a midgap YSR bound state is thermally activated.","lead":"This paper derives the temperature-dependent free energy and entropy of a magnetic impurity sitting at the edge of a superconducting wire, an exactly solvable one-dimensional model. It predicts that at intermediate temperatures the impurity entropy climbs above the free-spin value ln 2 in the Yu-Shiba-Rusinov phases, a fingerprint of a midgap bound state that entropy-measurement experiments could detect.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tower-counting fractions in YSR-I/II are asserted without proof; if they are wrong, the partition-function decomposition and entropy overshoot are not established.","rationale":"The reader's verdict is CONDITIONAL for essentially this reason, and I agree with that assessment. I single out the tower-counting fractions as the most load-bearing concern because every quantitative result—including the entropy overshoot and the closed-form free-energy expressions—is a weighted sum over towers whose weights are the asserted fractions. The contour-shift issue flagged by the reader is a serious technical question, but the tower-counting gap is more fundamental: without the fractions, the Hilbert-space splitting itself is not established. The exact-diagonalization results in Sec. S.5 give qualitative support for the overshoot phenomenon in a related lattice model, but they do not verify the specific 3:1 or 4:3:1 weights. Therefore the central claim should remain conditional pending an explicit counting derivation or a numerical confirmation of the tower fractions.","tokens_in":27912,"tokens_out":15122,"duration_ms":162179,"concrete_test":"Independently count the number of Bethe-ansatz states in each tower for finite N using the root densities of Eqs. (S.15), (S.17), and (S.21): enumerate all hole/string configurations on |U⟩, |B⟩, and |Ũ⟩ for representative N (e.g., 3, 5, 7) at δ values in YSR-I (δ=0.75) and YSR-II (δ=1.3), classify each solution by the presence of the boundary roots λδ and λδ+i, and verify the N→∞ ratios 3:1 and 4:3:1. If the ratios differ, the tower decomposition and the derived free energy are incorrect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that the impurity partition function factorizes into towers with exactly known weights—depends on the asserted Hilbert-space counting fractions: T1:T2 = 3:1 in YSR-I and T1:T2:T3 = 4:3:1 in YSR-II (and analogous formulas for δ>3/2). These fractions are stated without derivation in the main text (after Fig. 2). The supplement (Sec. S.1C) merely says 'counting the states ... shows' that two towers do not exhaust the Hilbert space, then asserts the third tower; no counting calculation is shown for any phase. If the true fractions differ—e.g., if the boundary strings overcount or the base states |U⟩, |B⟩, |Ũ⟩ are not orthogonal—then Z=Σ_i e^{-βF_i} is not the impurity partition function, and the predicted entropy overshoot and saturation values (ln 2, ln 3/2, etc.) are not the model's thermodynamics. The ED in Sec. S.5 demonstrates non-monotonic entropy in a related lattice model, but it does not test the claimed tower weights or the closed-form free-energy expressions. This is a load-bearing, unproven step.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter studies the finite-temperature impurity thermodynamics of an integrable one-dimensional superconductor with a boundary magnetic impurity. Using the Bethe-ansatz solution of Refs. [4,5], the authors propose that the Hilbert space in each phase decomposes into distinct 'excitation towers' whose number changes across the four phases (Kondo, YSR I, YSR II, local moment). They derive TBA-based expressions for the impurity free energy and entropy in each phase, and report a monotonic Kondo entropy flow, non-monotonic YSR/local-moment entropies with an overshoot above ln 2, and closed-form saturation values. Numerical TBA solutions and exact diagonalization of a small related lattice model are presented in the supplementary material to support the qualitative entropy curves.","tokens_in":28324,"tokens_out":14653,"duration_ms":161362,"significance":"If the central derivation is correct, this is a substantial contribution to integrable impurity thermodynamics: it extends beyond conventional TBA to a gapped superconductor with boundary-bound modes, gives explicit tower-resolved free energies, and makes concrete, falsifiable predictions for impurity entropy overshoots and saturation values. The work is parameter-free, builds on previously solved Bethe-ansatz equations rather than fitting, and includes both numerical TBA and exact-diagonalization checks. These are real strengths. However, the manuscript currently contains a sign inconsistency in a central equation, an unproven tower-counting step on which the decomposition rests, and a contour-shift step in the supplementary that is not justified. These issues must be resolved before the central claims can be accepted.","major_comments":[{"comment":"The printed sign in Eq. (10) contradicts the text and the supplementary derivation. As written, F(T2)_imp = -T/4 ∫ Σ ln(1+η1)/cosh(π(λ+iυδ)); with η1(T→∞)→3 this gives S(T2)(∞)=+ln2, whereas the text and Fig. 1 require S(T2)(∞)=-ln2. The supplementary derivation in Eq. (S.69) has the opposite overall sign. Since the entropy overshoot decomposition and the values S(T2)(∞)=-ln2, S(T1)(∞)=ln(3/2) follow directly from these formulas, Eq. (10) must be corrected and the two equations reconciled.","section":"Main text, Eq. (10); Supplementary Eq. (S.69)"},{"comment":"The tower-counting fractions are asserted without proof. The main text states that in YSR-I the towers contain 3/4 and 1/4 of the states, and in YSR-II the fractions are 4/6, 1/4, and 1/12 (with analogous δ>3/2 formulas), but no counting calculation is shown. The supplement merely says that 'counting the states ... shows' that two towers do not exhaust the Hilbert space, and then asserts a third tower. The partition-function decomposition Z=Σ_i e^{-βF_Ti}, the saturation entropies, and the entropy overshoot all depend on these towers and their relative weights. The ED in Sec. S.5 tests a related lattice model and does not verify the tower weights or the closed-form free-energy expressions. A derivation or an explicit state-count in the supplement is needed.","section":"Main text after Fig. 2; Supplementary Sec. S.1C"},{"comment":"The analytic continuation/contour-shift step is not justified. Eq. (S.65) is stated to be valid only for |ζ|<1/2, but the passage from the first to the second expression in Eq. (S.69) replaces a shift δ-1 (which lies inside the strip for δ∈(1/2,1)) by δ (which lies outside the strip). The contour shift can pick up pole contributions from the 1/cosh kernel. The authors should either show these residues vanish, or identify them explicitly and show that they are already accounted for by the |Eδ| term in Eq. (10). As written, the derivation of the central YSR-I free energy is incomplete.","section":"Supplementary Eq. (S.65) and Eq. (S.69)"}],"minor_comments":[{"comment":"The abstract claims 'closed-form analytic expressions across the entire phase diagram.' At finite temperature, η_n(λ) are obtained by numerically solving the infinite TBA hierarchy (Supplementary Sec. S.3); only the asymptotics are literally closed form. Please qualify the wording, e.g., 'exact integral representations in terms of the TBA functions.'","section":"Abstract"},{"comment":"The phrase '∀d ∪ δ ∈ (0,1)' is imprecise; d is real or imaginary and δ is defined for d=iδ. Please restate the parameter ranges explicitly.","section":"Main text after Eq. (10)"},{"comment":"F(T2)_0 is introduced and then set to zero later; the reader should be told at the point of introduction that this is a constant ground-state energy and will be absorbed.","section":"Supplementary Sec. S.2B"},{"comment":"The numerical TBA solution is validated only by empirical convergence checks. The authors state no formal error bounds; this is acceptable for a numerical study, but it should be labeled as numerical evidence rather than part of the exact analytic derivation.","section":"Supplementary Sec. S.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the authors' own prior work (Refs. [4,5]) for the Bethe-ansatz solution and phase classification. This is legitimate and not circular, but for a Letter with this level of technical claim, the tower-counting proof and the contour-shift justification should be made self-contained in the supplement. The sign inconsistency in Eq. (10) is a concrete, fixable internal error that currently prevents a reader from reproducing the reported entropy curves."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a serious TBA paper with a genuinely new result, and it is also a paper whose central quantitative claim is not yet supported. The new content is the finite-temperature impurity free energy and entropy for the boundary mode of the 1D gapped superconductor, with the prediction of entropy overshoot above ln 2 in the YSR phases. That overshoot is novel and, as they note, measurable in quantum-dot entropy experiments. The derivation from the previously solved Bethe ansatz is elaborate and the qualitative ED on a lattice model in the supplement supports the non-monotonic entropy. So there is real substance.\n\nThe soft spots are concrete. First, Eq. (10) in the main text has the wrong sign for F^(T2)_imp. As printed it gives S^(T2)(∞)=+ln2, while the text states -ln2, and the supplement Eq. (S.69) uses the opposite sign. That's a mismatch between main text and supplement, and it propagates to the stated tower entropy values.\n\nSecond, and more serious: the tower-counting fractions—3/4 and 1/4 in YSR-I, 4/6, 1/4, 1/12 in YSR-II, and the analogous formulas for δ>3/2—are asserted without proof. The partition function decomposition Z=Σ e^{-βF_i} depends entirely on these fractions, and the entropy overshoot would change if they are wrong. The supplement says 'counting the states shows' but does not show the count anywhere. That is a load-bearing unproven step, not a stylistic detail.\n\nThird, the contour-shift step in the supplement moves arguments beyond the stated analyticity strip |ζ|<1/2, which can pick up pole contributions. This may be fixable, but it needs to be checked.\n\nFinally, the abstract's 'closed-form analytic expressions' overstates the case: the η_n(λ) are obtained numerically for each temperature; the closed forms are only in the T→0 and T→∞ limits.\n\nShould it be reviewed? Yes—this is exactly the kind of paper referees should see: a first exact finite-T thermodynamics for a boundary impurity in a gapped host, with a measurable prediction. But the referee should demand a proof of the counting and a corrected Eq. (10) before the overshoot prediction is accepted. As it stands, I would not build on the tower decomposition or cite the entropy overshoot quantitatively.","headline":"New TBA derivation yielding a measurable entropy overshoot, but the main text has a sign error and the tower-counting fractions are asserted, so the central quantitative claim isn't yet established.","tokens_in":28661,"tokens_out":5095,"would_cite":false,"duration_ms":45952,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B23","82D55"],"pacs":[],"model":"deepseek-v4-flash","headline":"A magnetic impurity at the edge of a one-dimensional superconductor can carry an entropy that overshoots the free-spin value ln 2 at intermediate temperatures.","keywords":["Yu–Shiba–Rusinov bound states","impurity entropy overshoot","thermodynamic Bethe ansatz","Hilbert space fragmentation into towers","Kondo impurity in a superconductor","boundary quantum phase transition","midgap bound states","exactly solvable impurity model"],"falsifier":"Compute the tower weights by exact diagonalization of the attractive Hubbard chain with an edge impurity used in the supplement, counting eigenstates in each boundary-root sector: if the numbers do not match the 3:1 ratio in YSR I and the 8:3:1 ratio in YSR II, the decomposition fails. Alternatively, measure impurity entropy of a quantum dot at the end of a superconducting nanowire via charge sensing: the central claim predicts a bump above ln 2 near T ≈ |Eδ| that should move and disappear as the gate-tuned phase boundary δ = 1 is crossed.","tokens_in":27846,"feed_emoji":"🧲","tokens_out":6996,"duration_ms":73692,"temperature":0.7,"pith_summary":"The paper studies a single spin-1/2 impurity at the open edge of a one-dimensional superconducting wire, a model that is exactly solvable by Bethe ansatz. It claims that the finite-temperature impurity entropy is organized by a splitting of the Hilbert space into distinct towers of excitations whose relative sizes are phase-dependent. In the Kondo phase the entropy flows monotonically from ln 2 at high temperature to 0 at low temperature, with the same critical exponents as the conventional Kondo model. In the Yu-Shiba-Rusinov (YSR) phases the entropy is non-monotonic and can overshoot ln 2, because thermal activation of a midgap bound state temporarily adds entropy; in the local-moment phase it stays near ln 2 with only small bumps that fade as δ grows. If correct, the paper provides closed-form analytic expressions for the impurity free energy and entropy across all four phases, including thermodynamic signatures that do not exist in a gapless host.","feed_headline":"Impurity entropy overshoots ln 2 in a superconducting wire","feed_subtitle":"Midgap Yu–Shiba–Rusinov states add an entropy bump near their energy, and the paper gives exact formulas for all phases.","key_machinery":"Excitation towers: distinguished families of Bethe-ansatz eigenstates built on a small set of base states—|K⟩, |U⟩, |B⟩, and |̃U⟩—that differ by purely imaginary boundary string solutions of the Bethe equations. The impurity partition function is the sum over towers of e^{−β F_imp^{(Ti)}}, and each tower free energy is a convolution of ln(1 + η_n) with 1/cosh kernels, where η_n solve the thermodynamic Bethe ansatz equations. The entropy decomposes into a weighted ensemble average over towers (the mixing term) plus an activation term peaked when the temperature matches the midgap bound-state energy Eδ.","core_discovery":"The central claim is that the spectrum of this boundary impurity problem organizes into excitation towers that change with the RG-invariant parameter δ, and that each tower carries its own free energy. In the Kondo regime all states form one tower built on the fully screened singlet |K⟩. For 1/2 < δ < 1 (YSR I) there are two towers built on the unscreened |U⟩ and screened |B⟩ states, with relative weights 3/4 and 1/4; for 1 < δ < 3/2 (YSR II) and δ > 3/2 (local moment) there are three towers, with weights (4/6, 1/4, 1/12) in YSR II and floor-function expressions in the local-moment regime. Summing tower Boltzmann weights gives the impurity free energy F_imp(T); differentiating gives the entr","pith_inferences":["I would expect the same entropy-overshoot mechanism to appear in other integrable boundary problems with boundary strings, such as an impurity at the edge of a gapped spin chain, where the floor-function tower weights could be checked by exact diagonalization of short chains.","A direct experimental test could use charge-sensing entropy measurements on a quantum dot attached to a superconductor-semiconductor nanowire: the bump should appear near T ≈ |Eδ| and move or disappear as the gate-tuned phase boundary δ = 1 is crossed.","The paper leaves open whether the tower structure survives away from integrability; if it does, the overshoot should be robust to weak non-integrable perturbations, which could be tested with tensor-network simulations of the lattice model.","The fixed-N variant's predicted negative entropy dips are a sharp falsifiable signature: a measurement protocol that fixes particle number should see a dip, while a grand-canonical setup should see only the overshoot."],"forward_implications":["The predicted entropy overshoot above ln 2 gives a measurable thermodynamic fingerprint of YSR midgap states, distinguishing them from the monotonic Kondo flow.","The Kondo phase retains the conventional critical exponents despite the gapped host, so the universal screening behavior is unchanged there.","The fixed-particle-number treatment in the supplement shifts tower energies by the hole mass m and can produce negative entropy dips near δ ≈ 0.5, making the grand-canonical and fixed-N predictions experimentally distinguishable.","The tower decomposition generalizes the TBA approach to any integrable impurity or defect that generates boundary bound states, not only YSR superconductors.","In the local-moment phase the impurity is only asymptotically decoupled: S_imp approaches ln 2 in both limits, but small intermediate-temperature bumps persist and vanish as δ → ∞."],"supporting_citations":[{"why":"Defines the integrable model of a Kondo impurity at the edge of a superconducting wire whose Bethe ansatz solution and phase structure this paper extends to thermodynamics.","marker":"[4]"},{"why":"Supplies the ground-state Bethe ansatz solution, the RG-invariant parameter d, and the classification into Kondo, YSR I, YSR II, and local-moment phases that the tower construction builds on.","marker":"[5]"},{"why":"Supplement containing the derivations of the tower-counting fractions, the free-energy integrals in each phase, the numerical TBA method, and the exact diagonalization checks on a lattice model.","marker":"[15]"},{"why":"Original Yang-Yang thermodynamic Bethe ansatz formalism underlying the TBA equations used in Eq. (5).","marker":"[16]"},{"why":"Standard reference for TBA equations and the string hypothesis that the tower decomposition assumes.","marker":"[17]"},{"why":"Provides the Bethe-ansatz/TBA solution of the conventional Kondo problem whose critical exponents and monotonic entropy flow are the comparison baseline.","marker":"[18]"},{"why":"Exact-results framework for magnetic alloys whose boundary free-energy approach the impurity contribution follows.","marker":"[19]"},{"why":"Establishes the numerical scheme for solving TBA equations with an explicit mass scale, used here to obtain η_n at finite temperature.","marker":"[20]"},{"why":"Defines the boundary entropy/ground-state degeneracy used to interpret S_imp at the RG fixed points.","marker":"[24]"},{"why":"Gives the Bethe ansatz solution of the attractive Hubbard chain, which is the lattice model used for the exact diagonalization support.","marker":"[30]"}],"fun_headline_variants":["Impurity entropy jumps past ln 2 in wire edge","Split Hilbert space reveals entropy overshoots","YSR towers add entropy bump beyond ln 2","Exact impurity entropy from tower splitting","Gapped host yields richer impurity phases"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The counting that each phase's Hilbert space is fully exhausted by the proposed towers (3/4 and 1/4 in YSR I; 4/6, 1/4, 1/12 in YSR II) and the analyticity of the contour shifts used to evaluate the tower free energies are the load-bearing premises; if either fails, the closed-form free energy and the entropy overshoot are not established.","fun_headline_variants_meta":{"raw":{"variants":["Impurity entropy jumps past ln 2 in wire edge","Split Hilbert space reveals entropy overshoots","YSR towers add entropy bump beyond ln 2","Exact impurity entropy from tower splitting","Gapped host yields richer impurity phases"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000195,"raw_usage":{"total_tokens":1273,"prompt_tokens":900,"completion_tokens":373,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":644,"completion_tokens_details":{"reasoning_tokens":304}},"tokens_in":644,"tokens_out":373,"duration_ms":5086,"temperature":1.0,"reasoning_tokens":304,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:48:54.657634+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the tower weights by exact diagonalization of the attractive Hubbard chain with an edge impurity used in the supplement, counting eigenstates in each boundary-root sector: if the numbers do not match the 3:1 ratio in YSR I and the 8:3:1 ratio in YSR II, the decomposition fails. Alternatively, measure impurity entropy of a quantum dot at the end of a superconducting nanowire via charge sensing: the central claim predicts a bump above ln 2 near T ≈ |Eδ| that should move and disappear as the gate-tuned phase boundary δ = 1 is crossed.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the integrable model of a Kondo impurity at the edge of a superconducting wire whose Bethe ansatz solution and phase structure this paper extends to thermodynamics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplement containing the derivations of the tower-counting fractions, the free-energy integrals in each phase, the numerical TBA method, and the exact diagonalization checks on a lattice model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original Yang-Yang thermodynamic Bethe ansatz formalism underlying the TBA equations used in Eq. (5)."},{"cited_title":"Takahashi et al., Thermodynamics of one-dimensional solvable models (Cambridge university press Cambridge, 1999)","cited_arxiv_id":null,"evidence_quote":"Standard reference for TBA equations and the string hypothesis that the tower decomposition assumes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the numerical scheme for solving TBA equations with an explicit mass scale, used here to obtain η_n at finite temperature."},{"cited_title":"ln 1 + η⌈2δ⌉(λ) cosh π λ + iυ 2 (2δ − ⌊2δ⌋) − ln 1 + η⌊2δ⌋(λ) cosh π λ + iυ 2 (⌈2δ⌉ −2δ) # , F (T2) imp = T 4 X υ=± Z dλ","cited_arxiv_id":null,"evidence_quote":"Gives the Bethe ansatz solution of the attractive Hubbard chain, which is the lattice model used for the exact diagonalization support."}],"review_version":1}