{"id":"b197a2a8-4f37-41c0-86a5-8545ccb201d9","arxiv_id":"2505.16776","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The NESS spectrum of boundary-driven integrable spin chains in the Zeno regime is expressed through coherent Bethe ansatz eigenstates with modified single-particle dispersions.","lead":"Researchers show that the steady state of a quantum spin chain with strong boundary dissipation can be described by the same quasiparticles as the closed system, but with a renormalized, dissipatively dressed dispersion relation. This maps a non-equilibrium problem onto a known Bethe ansatz structure, which may aid dissipative state preparation and the study of open integrable systems.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Additivity of dressed energies for M>1 in the chiral XXZ/XYZ cases is the load-bearing conjecture; numerical checks test final ratios, not the factorization itself.","rationale":"The reader's weakest assumption is exactly where the argument is least secure: the extension of the dressed dispersion from M=1 to M>1 in the non-U(1) cases. My reading confirms this. The main text's Section V is an honest and complete M=1 proof for the chiral XXZ case; the authors explicitly mark the multi-particle extension as a conjecture. The companion paper [6] supplies the missing proof for the U(1) case, so the gap is not uniform across models. The numerical checks in Appendix B (Tables 1-2) and Appendix C (Tables 3-4) are consistent with Eq. (36)/(43) at small N, and this is real evidence. However, those tables compare the final NESS population ratios with the product formula; they do not test the factorization property that distinguishes additivity from a coincidental fit. Because Eq. (20) is only numerically verified (Remark 2) and the derivation of (100) explicitly invokes BAE and non-orthogonality cancellations that become much harder for M>1, the risk is concrete rather than abstract. The right response is not rejection: the authors flag the conjecture, provide a proof in the companion for the U(1) case, and give reproducible numerical tables. CONDITIONAL is the correct verdict. My stress-test therefore leaves the reader's verdict unchanged.","tokens_in":21400,"tokens_out":11362,"duration_ms":84154,"concrete_test":"Derive the M=2 analogue of Eq. (100) in the chiral XXZ model by computing the two-kink form factors ⟨α|g_l|β⟩ with the coordinate Bethe ansatz of Ref. [17] and the non-orthogonal chiral basis (81), without assuming detailed balance. If the resulting population ratio ν_α/ν_β is exactly [∏_{j=1}^2 F(u_{j,α})]/[∏_{j=1}^2 F(u_{j,β})] with F(u)=|Δ+Δ^{-1}-2 cos p(u)|^2, then Eq. (36) for M=2 is proven and additivity is confirmed; if cross terms involving both rapidities survive the N-dependence cancellation that makes the M=1 derivation work, then the dressed energy is not additive and the central claim fails for M>1. Repeat with the N=6, M=3 data of Table 2 as a cross-check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (3)-(4) requires that for every multi-particle state |α⟩ in the chiral/XYZ sectors, log ν_α = Σ_j ε̃(u_{j,α}) with the same rapidities as H_D. This additivity is the load-bearing premise. For the U(1) sink-source case it is proven for all M in the companion paper [6], and for M=1 in the chiral XXZ case it is derived in Section V. But for M>1 in the chiral XXZ case the paper explicitly labels Eq. (36) a conjecture based on Appendix B, and for the XYZ model the same statement is Hypothesis 1 in Appendix C (Eq. A-10). The numerical tables in Appendices B and C verify the product formula at selected M=2/N=5 and M=3/N=6 points, but they do not isolate the additivity requirement: a generic symmetric function of the M Bethe rapidities can agree with the sum of single-particle terms on the tested states without factorizing. In addition, the Kolmogorov relation (20), which licenses the detailed-balance shortcut (21) used to obtain these ratios, is itself only verified numerically for these non-U(1) cases. Thus, if additivity fails at larger M or N (or if Kolmogorov fails outside the tested parameter window), the dressed-quasiparticle picture in the strongly anisotropic/non-diagonal models would not hold in the form claimed. The authors are transparent about the conjecture, which is why this is a condition for acceptance rather than a refutation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the nonequilibrium steady state (NESS) of boundary-driven integrable spin chains in the quantum Zeno limit, proposing that the NESS has the same eigenstates as the dissipation-projected Hamiltonian and that its eigenvalues are given by a Gibbs-like formula with a renormalized, 'dissipatively dressed' dispersion relation. The central claim, stated in Eqs. (3)-(4), is that for every Bethe state |α⟩ the NESS weight obeys log ν_α = Σ_j ε̃(u_{j,α}) with the same Bethe rapidities as the coherent model. Explicit dressed dispersions are derived for the XXX sink-source case (Eq. (26)), the diagonal XXZ case (Eq. (31)), the chiral XXZ case (Eq. (36)), and conjectured for the XYZ case (Eq. (43)). The derivation is rigorous for M=1 in the sink-source and chiral XXZ sectors and for all M in the U(1) sink-source case via the companion paper [6]; for M>1 in the chiral XXZ and XYZ cases the additivity of dressed energies is presented as a conjecture supported by numerical tables in Appendices B and C.","tokens_in":21693,"tokens_out":4717,"duration_ms":42508,"significance":"If the central claim holds, the paper establishes a striking and useful phenomenon: the quasiparticle content of an integrable system can survive strong boundary dissipation, with the only effect being a renormalization of the dispersion relation. The result would provide an analytic handle on the NESS spectrum of boundary-driven integrable chains beyond the U(1)-symmetric cases, and it has potential implications for dissipative state preparation and for understanding emergent integrability in open quantum circuits. The strengths of the manuscript are the exact algebraic Bethe ansatz derivations for the M=1 sectors, the explicit closed-form expressions for the dressed dispersions, and the transparent statement of which parts are conjectural. The numerical tables in the appendices provide supporting evidence, although, as discussed below, they do not yet fully test the specific additivity property that is load-bearing for the multi-particle claim.","major_comments":[{"comment":"The generalization of Eq. (36) to M>1 is explicitly labeled a conjecture, and the numerical evidence in Tables 1 and 2 does not actually display the predicted dressed-energy differences. To verify the central factorization log ν_α = Σ_j ε̃(p_j(α)), one must compare log(να/ν1) with Σ_j ε̃(p_j(α)) − Σ_j ε̃(p_j(1)); the tables list E and the NESS ratio but not this comparison. As written, the numerical check could be consistent with a more general symmetric function of the rapidities. Please add the predicted column, or better, test the additivity state-by-state by computing log ν_α − Σ_j ε̃(p_j(α)) for each α.","section":"Section V and Appendix B, Eqs. (36) and (A-1)"},{"comment":"For the XYZ model, the multi-particle formula is assumed as Hypothesis 1, and the a,b values in Table (A-12) are system-dependent parameters extracted from M=1 numerics. The later identification a=0, b=(1−τ)/2+iy_l is presented as a conclusion, but Tables 3 and 4 again verify final ratios rather than the additivity of dressed energies over rapidities. Since this is the only evidence for the XYZ dressed dispersion (43), the claim is not yet load-bearing; either provide a derivation or present a dedicated numerical test that separates the additivity hypothesis from the detailed-balance ratios.","section":"Appendix C, Hypothesis 1 and Eq. (A-10)"},{"comment":"The Kolmogorov relation (20) is 'observed numerically' and called 'postulated' in Remark 2, yet the detailed-balance relation (21) is used to obtain all NESS ratios, including the XYZ data in Appendix C. For the chiral XXZ case the relation is proved a posteriori from Eq. (100), but for the XYZ case no proof is given. If Eq. (20) fails for parameters outside the tested window, the ratios used in Appendix C would not be the NESS weights. Please either prove Eq. (20) for the XYZ case or provide a systematic scan over τ, η, x_l, y_l, N, and M that demonstrates its validity.","section":"Section II, Eqs. (20)-(21)"},{"comment":"The statement that Eq. (31) holds 'for all eigenstates' relies entirely on the companion paper [6], which is an arXiv preprint rather than a proof contained in this manuscript. Since the universal claim of the paper depends on this external result, please clarify whether [6] has been peer reviewed or provide the essential steps of the all-M proof here.","section":"Section VI and companion paper [6]"}],"minor_comments":[{"comment":"Reference [3] contains a typo: 'Rush. Math. Surveys' should be 'Russ. Math. Surveys'.","section":"References"},{"comment":"The phrase 'plain-wave like Bethe state' should read 'plane-wave like Bethe state'; the same typo appears in the main text as 'plane wave type'.","section":"Fig. 2 caption and Section IV"},{"comment":"The eigenvalues of the Kraus map are listed as λ_3 = λ_4 = √ϵ, but the eigenstates are indexed ψ_0,…,ψ_3; the notation should be λ_2 = λ_3 = √ϵ.","section":"Appendix A, Eq. (A-6)"},{"comment":"The reduction of the ratio w_{0α}/w_{α0} to |γ_1/γ_N|² skips an intermediate explanation of why the two terms in the denominator coincide; adding one sentence would improve readability.","section":"Section IV, Eqs. (48)-(49)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest and well structured, and the M=1 derivations are solid. The main concern is the gap between the broad claim of dissipative dressing and the proven cases: the multi-particle additivity in the chiral XXZ and XYZ cases is the load-bearing premise, and the current numerical evidence does not isolate it. The paper would be suitable for publication after the authors either strengthen the evidence with state-by-state additivity tests and a proof or systematic check of the Kolmogorov relation for the XYZ case, or clearly reclassify the non-U(1) multi-particle statements as conjectures in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plainly: the paper proves the dressed-quasiparticle picture for the one-particle sector of several integrable chains, and conjectures it for multi-particle sectors in the non-diagonal cases. The M=1 derivations are rigorous, the paper is honest about the conjectural parts, and the new formulas for chiral XXZ and XYZ are worth having.\n\nThe genuinely new content is the dressed dispersion for XXZ with non-diagonal boundary fields (Eq. 36) and for XYZ (Eq. 43), plus the careful derivation for M=1 in the chiral XXZ case. The sink-source XXX/XXZ results are already in the companion paper [6], but the present paper gives a more pedagogical algebraic-Bethe derivation. I appreciate that the authors clearly mark what is proved, what is conjectured, and where the Kolmogorov relation is verified numerically rather than derived.\n\nThe soft spot is exactly what the stress-test identifies. The central claim (3)-(4) requires additivity of the dressed energies over Bethe rapidities for every multi-particle state. That is proven for U(1) sink-source in [6] and for M=1 here, but for M>1 in the chiral XXZ and XYZ cases it is a conjecture. The numerical tables test the product formula at selected M=2/N=5 and M=3/N=6 points; they do not isolate the factorization step. A generic symmetric function could match those ratios without being a sum of single-particle terms. The Kolmogorov relation is also only numerically verified in these non-U(1) cases, though for chiral XXZ M=1 it follows from the derived rate ratio (Eq. 100). These are real limitations, but they are clearly stated, so the paper is conditional rather than wrong.\n\nOne smaller point: the abstract's 'substantially altering' overstates the effect for the chiral cases, where all Bethe roots are real and no localized state appears; the dramatic effect in Fig. 2 is specific to the sink-source setup.\n\nThis is a paper for researchers working on integrable open systems and NESS. It deserves a serious referee. The main request should be for stronger evidence or a proof for the M>1 additivity in the chiral/XYZ settings, or a sharper statement that the multi-particle dressed picture is currently a conjecture there.","headline":"Clear, honest paper proving the dressed-quasiparticle picture for one-particle sectors and conjecturing it for multi-particle sectors; deserves peer review with emphasis on the additivity conjecture.","tokens_in":22238,"tokens_out":2546,"would_cite":true,"duration_ms":19753,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82C23","82C10","81R12"],"pacs":["05.30.-d","75.10.Pq","03.65.Yz"],"model":"deepseek-v4-flash","headline":"This paper claims that in the quantum Zeno limit, the nonequilibrium steady state of a boundary-driven integrable spin chain is a Gibbs-like mixture over the same Bethe eigenstates as the coherent chain, with the quasiparticle dispersion…","keywords":["dissipative dressing","non-equilibrium steady state","integrable spin chains","quantum Zeno regime","Bethe ansatz","Lindblad master equation","boundary driving","quasiparticle dispersion"],"falsifier":"Compute the exact NESS populations $\\nu_\\alpha$ for $M=2$ in the chiral XXZ chain (or $M=2$ in the XYZ chain) by full diagonalization of the Zeno effective Lindblad dynamics for moderate $N$, and compare $-\\log\\nu_\\alpha$ with $\\tilde{\\epsilon}(u_1^{(\\alpha)})+\\tilde{\\epsilon}(u_2^{(\\alpha)})$ for every pair of Bethe roots; a single pair where $\\log(\\nu_\\beta/\\nu_\\alpha)$ differs from $\\tilde{\\epsilon}(u_1^{(\\alpha)})+\\tilde{\\epsilon}(u_2^{(\\alpha)})-\\tilde{\\epsilon}(u_1^{(\\beta)})-\\tilde{\\epsilon}(u_2^{(\\beta)})$ disproves additivity.","tokens_in":21137,"feed_emoji":"🧲","tokens_out":6969,"duration_ms":51178,"temperature":0.7,"pith_summary":"This paper argues that when an integrable spin chain is driven at its boundaries by strong dissipation, the resulting nonequilibrium steady state is not structureless: it is a Gibbs-like mixture over the same Bethe eigenstates as the coherent chain, with the quasiparticle dispersion $\\epsilon(u)$ replaced by a dissipatively dressed dispersion $\\tilde{\\epsilon}(u)$. The claim is that the steady-state \"energies\" are additive over Bethe rapidities exactly as in the coherent spectrum, so the quasiparticle content survives dissipation and only the dispersion relation is renormalized. Explicit analytic dressed dispersions are given for XXX and XXZ chains with diagonal boundary driving, and for XXZ and XYZ chains with non-diagonal chiral driving. If true, this gives an exact handle on far-from-equilibrium steady states of interacting many-body systems and links boundary-driven dissipative dynamics to Bethe ansatz integrability.","feed_headline":"Dissipation dresses quasiparticle energies in integrable spin chains","feed_subtitle":"Boundary-driven spin chains reach steady states that mirror equilibrium, only the quasiparticle energies are renormalized.","key_machinery":"The load-bearing construction is the dissipation-projected Hamiltonian $H_D = H_{\\rm bulk}+h_1+h_N$, whose boundary fields are fixed by the polarizations onto which the boundary spins are projected in the Zeno limit; its eigenstates form the common basis of the coherent and dissipative systems. The second ingredient is the classical Markov process for populations $\\nu_\\alpha$ with rates $w_{\\alpha\\beta}=|\\langle\\alpha|g_l|\\beta\\rangle|^2+|\\langle\\alpha|g_r|\\beta\\rangle|^2$, whose stationary solution is the NESS spectrum. The paper uses the Kolmogorov condition $w_{\\alpha\\beta}w_{\\beta\\gamma}w_{\\gamma\\alpha}=w_{\\alpha\\gamma}w_{\\gamma\\beta}w_{\\beta\\alpha}$, which implies detailed balance, so $\\nu_\\alpha/\\nu_\\beta$ reduces to a ratio of boundary correlation functions of Bethe eigenstates. Evaluating those ratios with the algebraic or chiral Bethe ansatz yields the dressed dispersion $\\tilde{\\epsilon}(u)$ in closed form.","core_discovery":"The central discovery is the one-to-one correspondence between the coherent spectrum and the NESS spectrum in the quantum Zeno regime: $\\rho_{\\rm NESS} = \\tilde{Z}^{-1}\\sum_\\alpha e^{-\\tilde{E}_\\alpha}|\\alpha\\rangle\\langle\\alpha|$, where $|\\alpha\\rangle$ are eigenstates of the dissipation-projected Hamiltonian $H_D$ and $\\tilde{E}_\\alpha = \\sum_j \\tilde{\\epsilon}(u_{j,\\alpha})$ uses the same Bethe rapidities as the coherent model. The paper derives the dressed dispersions from ratios of boundary correlation functions enforced by detailed balance: for the XXX and XXZ sink-source cases $\\tilde{\\epsilon}(u)=\\log|1-\\epsilon(u)|$ and $\\tilde{\\epsilon}(u)=\\log|1-\\Delta\\epsilon(u)|$, for the chiral XXZ case $\\tilde{\\epsilon}(u)=2\\log|\\cosh(u+i\\gamma/2)\\cosh(u-i\\gamma/2)/(\\sinh(u+i\\gamma/2)\\sinh(u-i\\gamma/2))|$, and for the XYZ case an elliptic expression. A further claim is that the dressing always adds a singularity to the dispersion, and in the boundary-localized sector this singularity produces exponentially large weights in the NESS.","pith_inferences":["If additivity holds beyond the proved cases, then the NESS is a genuine generalized Gibbs ensemble for the dissipation-projected Hamiltonian, and boundary dissipation could become a tool to prepare states with prescribed quasiparticle content by engineering the dispersion through boundary polarizations.","The Kolmogorov property of the rates may itself be a hidden integrable structure of the auxiliary Markov process; it would be worth checking whether it follows from the Yang-Baxter relation rather than being an independent assumption.","A quantum-circuit implementation with reset gates could test the predicted exponentially large weight of the localized Bethe state, which should appear as a sharp feature in the steady-state magnetization profile or structure factor.","The predicted level-order reversal in the chiral XXZ chain could be probed experimentally by quench spectroscopy on a small chain, since the NESS ordering should follow $\\tilde{\\epsilon}$ rather than $\\epsilon$."],"forward_implications":["In the U(1)-symmetric sink-source XXX/XXZ case, the dressed-dispersion formula is valid for all eigenstates, so the NESS spectrum is exactly computable for chains of any length in the Zeno limit.","In the XXX case, the Bethe root of the boundary-localized state sits exponentially close to the new singularity $u=3i/2$, giving that state a weight exponentially large in $N$ and making the NESS entropy subextensive.","In the chiral XXZ case, because all Bethe roots are real, no root approaches the imaginary-axis singularities, so the dressing reverses the ordering of level contributions without creating exponentially dominant states.","For the XYZ model, the conjectured elliptic dressed dispersion reduces to the chiral XXZ result in the appropriate trigonometric limit, so the dressing mechanism extends to the fully anisotropic chain.","In all cases the dressed dispersion has the form $\\tilde{\\epsilon}(u)\\sim\\log(\\epsilon(u)/f(u))$ with $f$ carrying an extra singularity, so the NESS spectrum is obtained from the coherent spectrum by a universal type of renormalization."],"supporting_citations":[{"why":"companion paper proving the dressed-dispersion result for the U(1)-symmetric sink-source case for all eigenstates","marker":"[6]"},{"why":"introduces the dissipation-projected Hamiltonian and the Zeno-limit condition underpinning the commuting NESS","marker":"[12]"},{"why":"provides the effective Lindblad master equation for internal spins and the rate-equation form used to find NESS populations","marker":"[13]"},{"why":"supplies the Sklyanin Bethe ansatz for the XXZ chain with boundary fields used in the diagonal-driving derivation","marker":"[14]"},{"why":"constructs the chiral invariant subspace for the XXZ model with non-diagonal boundary fields","marker":"[15]"},{"why":"gives the chiral coordinate Bethe ansatz and the Bethe equations for the kink eigenstates","marker":"[17]"},{"why":"constructs the chiral invariant subspace and boundary-field parametrization for the XYZ model","marker":"[19]"},{"why":"supplies the XYZ Bethe equations used for the elliptic dressed dispersion","marker":"[20]"}],"fun_headline_variants":["Dissipative dressing renormalizes quasiparticle energies in spin chains","How dissipation dresses quasiparticles in integrable spin chains","Boundary-driven spin chains: quasiparticles get dissipatively dressed","Dressed quasiparticles explain steady states of driven spin chains","Spin chains: dissipation adds a singularity to quasiparticle spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole picture rests on the claim that the steady-state \"energy\" of an $M$-quasiparticle state is the sum of $M$ single-particle dressed dispersions evaluated at the same Bethe rapidities; this additivity is proved for all $M$ only in the U(1) sink-source case, and for the chiral XXZ and XYZ cases it is a conjecture backed by numerics.","fun_headline_variants_meta":{"raw":{"variants":["Dissipative dressing renormalizes quasiparticle energies in spin chains","How dissipation dresses quasiparticles in integrable spin chains","Boundary-driven spin chains: quasiparticles get dissipatively dressed","Dressed quasiparticles explain steady states of driven spin chains","Spin chains: dissipation adds a singularity to quasiparticle spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000873,"raw_usage":{"total_tokens":3777,"prompt_tokens":944,"completion_tokens":2833,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":560,"completion_tokens_details":{"reasoning_tokens":2745}},"tokens_in":560,"tokens_out":2833,"duration_ms":21057,"temperature":1.0,"reasoning_tokens":2745,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:54:50.257278+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact NESS populations $\\nu_\\alpha$ for $M=2$ in the chiral XXZ chain (or $M=2$ in the XYZ chain) by full diagonalization of the Zeno effective Lindblad dynamics for moderate $N$, and compare $-\\log\\nu_\\alpha$ with $\\tilde{\\epsilon}(u_1^{(\\alpha)})+\\tilde{\\epsilon}(u_2^{(\\alpha)})$ for every pair of Bethe roots; a single pair where $\\log(\\nu_\\beta/\\nu_\\alpha)$ differs from $\\tilde{\\epsilon}(u_1^{(\\alpha)})+\\tilde{\\epsilon}(u_2^{(\\alpha)})-\\tilde{\\epsilon}(u_1^{(\\beta)})-\\tilde{\\epsilon}(u_2^{(\\beta)})$ disproves additivity.","supporting_citations":[{"cited_title":"Popkov, X","cited_arxiv_id":null,"evidence_quote":"companion paper proving the dressed-dispersion result for the U(1)-symmetric sink-source case for all eigenstates"},{"cited_title":"Zanardi and L","cited_arxiv_id":null,"evidence_quote":"introduces the dissipation-projected Hamiltonian and the Zeno-limit condition underpinning the commuting NESS"},{"cited_title":"Popkov, S","cited_arxiv_id":null,"evidence_quote":"provides the effective Lindblad master equation for internal spins and the rate-equation form used to find NESS populations"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the Sklyanin Bethe ansatz for the XXZ chain with boundary fields used in the diagonal-driving derivation"},{"cited_title":"Zhang, A","cited_arxiv_id":null,"evidence_quote":"constructs the chiral invariant subspace for the XXZ model with non-diagonal boundary fields"},{"cited_title":"Zhang, A","cited_arxiv_id":null,"evidence_quote":"gives the chiral coordinate Bethe ansatz and the Bethe equations for the kink eigenstates"},{"cited_title":"Yang and Y.-Z","cited_arxiv_id":null,"evidence_quote":"supplies the XYZ Bethe equations used for the elliptic dressed dispersion"}],"review_version":1}