Pith. sign in

REVIEW 4 major objections 4 minor 23 references

Dissipatively dressed quasiparticles in boundary driven integrable spin chains

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read 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…

desk verdict 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. read the letter →

arxiv 2505.16776 v1 pith:YYDNYT6V submitted 2025-05-22 cond-mat.stat-mech math-phmath.MPquant-ph

Vladislav Popkov , Xin Zhang , Carlo Presilla , Tomaž Prosen This is my paper · ORCID
classification cond-mat.stat-mechmath-phmath.MPquant-ph MSC 82C2382C1081R12 PACS 05.30.-d75.10.Pq03.65.Yz
keywords dissipativedressingnon-equilibriumsteadystateintegrablespinchainsquantumZenoregimeBetheansatzLindbladmasterequationboundarydrivingquasiparticledispersion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

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.

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 (4)
  1. [Section V and Appendix B, Eqs. (36) and (A-1)] 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 α.
  2. [Appendix C, Hypothesis 1 and Eq. (A-10)] 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.
  3. [Section II, Eqs. (20)-(21)] 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.
  4. [Section VI and companion paper [6]] 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.
minor comments (4)
  1. [References] Reference [3] contains a typo: 'Rush. Math. Surveys' should be 'Russ. Math. Surveys'.
  2. [Fig. 2 caption and Section IV] 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'.
  3. [Appendix A, Eq. (A-6)] The eigenvalues of the Kraus map are listed as λ_3 = λ_4 = √ϵ, but the eigenstates are indexed ψ_0,…,ψ_3; the notation should be λ_2 = λ_3 = √ϵ.
  4. [Section IV, Eqs. (48)-(49)] 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.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular reduction: multi-particle additivity in chiral/XYZ cases is openly conjectural, and the companion-paper citation is prior work.

full rationale

The central derivation is not circular. For the U(1) XXX/XXZ sink-source case, the M=1 dressed dispersion is obtained from exact algebraic Bethe-ansatz matrix elements (Eqs. (48)-(67)), not by assuming the product form (3)-(4); the all-M proof is delegated to the companion paper [6], a prior separate derivation by the same authors, which is legitimate external support rather than a circular reduction. In Section V, Eq. (82) uses the detailed-balance shortcut (21), which initially rests on the numerically observed Kolmogorov relation (20); however, the independent derivation of the rate ratio (100) from the wavefunctions and BAE (87) allows Kolmogorov to be verified a posteriori, so the final ratio (101) and dressed dispersion (36) are not assumed. For M>1 in the chiral XXZ and XYZ cases, the additivity E_tilde_alpha = sum_j epsilon_tilde(u_j) is explicitly flagged as a conjecture (Sec. III after Eq. (36); Hypothesis 1, Appendix C) and supported by numerical tables that compare the conjectured ratio with 'other approaches'; this is an honest limitation, not a circular step. The XYZ parameter identification (a=0, b=(1-tau)/2+iy_l) is fitted on M=1 data and then extrapolated, which is a testable prediction, not a fitted quantity renamed as a result. Overall, no equation reduces by construction to its input; the nonzero score reflects only the reliance on companion work and on numerically verified (rather than fully proven) Kolmogorov relations.

Assumptions & free parameters 1 free parameters · 5 assumptions · 1 invented entities

The central derivation relies on standard Bethe ansatz machinery plus two ad hoc postulates: the Kolmogorov relation (20) and the additivity of dressed energies for M>1. The latter is only conjectured for the non-diagonal boundary field cases. No fitting parameters remain in the final analytic expressions once the model parameters are fixed; the XYZ conjecture involved fitting a and b during discovery, but the final formulas determine them from model inputs.

free parameters (1)
  • a and b in XYZ numerical hypothesis = a=0, b=(1-tau)/2 + i y_l (ultimately fixed by model parameters; initial fits e.g. b=0.5-0.125i in table A-12)
    In Appendix C, the parameters a and b in the conjectured ratio formula (A-11) were initially determined by fitting numerical NESS spectra, as shown in table (A-12). The final expression (A-13) sets a=0 and b=(1-tau)/2 + i y_l, which are model parameters, so no free parameter remains in the final result. This entry records the exploratory fitting step.
assumptions (5)
  • domain assumption Quantum Zeno expansion: the NESS commutes with the dissipation-projected Hamiltonian (Eq. 9) and the effective dynamics is Eq. (13).
    This is the standard Zeno-limit projection from [12,13], used in Section II to establish property (A). It assumes fast boundary relaxation with 1/Gamma as a small parameter.
  • ad hoc to paper Kolmogorov relation (Eq. 20) for the transition rates w_alpha beta of the auxiliary classical master equation.
    Equation (20) is postulated in Section II and Remark 2, with the statement that its validity can be checked or proved aposteriori in specific cases. It is load-bearing because it leads to detailed balance (21), which converts rate ratios into NESS weight ratios.
  • ad hoc to paper Additivity of dressed energies over Bethe rapidities for M>1, i.e. tilde E_alpha = sum_j tilde epsilon(u_j), in the non-diagonal boundary field cases.
    For the chiral XXZ case this is a conjecture based on numerics (Appendix B), and for XYZ it is Hypothesis 1 (Appendix C). For the U(1) sink-source case it is proved in the companion paper [6]. This is the central structural assumption needed to extend the one-particle result to all sectors.
  • standard math Bethe ansatz integrability of the boundary-field Hamiltonian (5) for arbitrary J_x, J_y, J_z and boundary fields h_1, h_N.
    Used throughout Section III to parametrize eigenstates by Bethe rapidities, citing Refs. [7,8] and Sklyanin [14].
  • domain assumption The chiral invariant subspace and the Bethe ansatz equations (39) for the XYZ model, as given in Refs. [19,20].
    Section III and Appendix C rely on the existence and structure of the chiral subspace G_M and the elliptic BAE (39) from the literature to parametrize the eigenstates.
invented entities (1)
  • Dissipatively dressed quasiparticle
    purpose: Conceptual bookkeeping device to express NESS weights as exp(-sum_j tilde epsilon(u_j)), i.e. the same quasiparticles as the coherent system with a renormalized dispersion.
    This is not a new physical particle but a mathematical construct. Its validity is supported by the derived one-particle results and numerical checks, but it has no directly measurable signature independent of the NESS spectrum itself.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Dissipatively dressed quasiparticles in boundary driven integrable spin chains." pith.science (2026). https://pith.science/paper/YYDNYT6V

@misc{pith2026250516776,
  author       = {Pith},
  title        = {Pith review of: Dissipatively dressed quasiparticles in boundary driven integrable spin chains},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YYDNYT6V}},
  note         = {Machine review of arXiv:2505.16776}
}
read the original abstract

The nonequilibrium steady state (NESS) of integrable spin chains experiencing strong boundary dissipation is accounted by introducing quasiparticles with a renormalized -- dissipatively dressed -- dispersion relation. This allows us to evaluate the spectrum of the NESS in terms of the Bethe ansatz equations for a related coherent system which has the same set of eigenstates, the so-called dissipation-projected Hamiltonian. We find explicit analytic expressions for the dressed energies of the XXX and XXZ models with effective, i.e., induced by the dissipation, diagonal boundary fields, which are U(1) invariant, as well as the XXZ and XYZ models with effective non-diagonal boundary fields. In all cases, the dissipative dressing generates an extra singularity in the dispersion relation, substantially altering the NESS spectrum with respect to the spectrum of the corresponding coherent model.

Figures

Figures reproduced from arXiv: 2505.16776 by the authors.

Figure 1
Figure 1. Schematic picture of the dissipative setup. The [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Explicitly, from (23) we find u1 − 3i/2 ≡ δ = 3i2 −2N−1 (1 + O(Nδ)). The corresponding dressed en￾ergy ˜ϵ(u1) is drastically renormalized by the singularity, acquiring a negative amplitude linearly growing with sys￾tem size, ˜ϵ(u1) ≈ −2(N +1) log 2+log 9, see (80). On the other hand, for plane wave type (real uj ) BAE solutions, the dressed and original energies are comparable. As a result, e −ϵ˜(u1) ≫ e −ϵ˜(uα) for… view at source ↗
Figure 3
Figure 3. Quasiparticle energies ϵ(uα) (blue joined points) and dissipatively dressed ones ˜ϵ(uα) (red joined points) in the XXZ model with chiral invariant subspace, in the block with one kink M = 1. Parameters: N = 18, ∆ = 0.7. There are no localized states: all the amplitudes are plain-wave like, data not shown. proof only for the U(1)-symmetric “sink and source” case (31), which is valid for all eigenstates, see the compa… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 19 canonical work pages

  1. [6]

    Popkov, X

    V. Popkov, X. Zhang, C. Presilla, and T. Prosen, arXiv:2408.09302

    work page arXiv

  2. [1]

    Bethe, Zeitschrift f¨ ur Physik71, 205 (1931)

    H. Bethe, Zeitschrift f¨ ur Physik71, 205 (1931)

    1931

  3. [2]

    R. J. Baxter,Exactly solved models in statistical mechan- ics(Courier Corporation, 2007)

    work page 2007

  4. [3]

    L. A. Takhtadzhan and L. D. Faddeev, Rush. Math. Sur- veys34, 11 (1979)

    work page 1979

  5. [4]

    E. K. Sklyanin, L. A. Takhtadzhyan, and L. D. Faddeev, Theor. Math. Phys.40, 688 (1980)

    work page 1980

  6. [5]

    E. K. Sklyanin, Progress of Theoretical Physics Supple- ment118, 35 (1995)

    work page 1995

  7. [7]

    J. Cao, S. Cui, W.-L. Yang, K. Shi, and Y. Wang, Nucl. Phys. B886, 185 (2014)

    work page 2014

  8. [8]

    Wang, W.-L

    Y. Wang, W.-L. Yang, J. Cao, and K. Shi,Off-Diagonal Bethe Ansatz for Exactly Solvable Models(Springer, 2016)

    work page 2016

Show all 23 references
  1. [9]

    Popkov, D

    V. Popkov, D. Karevski, and G. M. Sch¨ utz, Phys. Rev. E88, 062118 (2013)

    2013

  2. [10]

    G. T. Landi, E. Novais, M. J. de Oliveira, and D. Karevski, Phys. Rev. E90, 042142 (2014)

    2014

  3. [11]

    X. Mi, A. A. Michailidis, S. Shabani, K. C. Miao, Klimov, et al., Science383, 1332 (2024)

    2024

  4. [12]

    Zanardi and L

    P. Zanardi and L. Campos Venuti, Phys. Rev. Lett.113, 240406 (2014)

    2014

  5. [13]

    Popkov, S

    V. Popkov, S. Essink, C. Presilla, and G. Sch¨ utz, Phys. Rev. A98, 052110 (2018)

    2018

  6. [14]

    E. K. Sklyanin, J. Phys. A21, 2375 (1988)

    1988

  7. [15]

    Zhang, A

    X. Zhang, A. Kl¨ umper, and V. Popkov, Phys. Rev. B 103, 115435 (2021)

    2021

  8. [16]

    Popkov and M

    V. Popkov and M. Salerno, Europhysics Lett.140, 11004 (2022)

    2022

  9. [17]

    Zhang, A

    X. Zhang, A. Kl¨ umper, and V. Popkov, Phys. Rev. B 104, 195409 (2021)

    2021

  10. [18]

    E. T. Whittaker and G. N. Watson,A course of modern analysis(Cambridge University Press, 1950)

    1950

  11. [19]

    Zhang, A

    X. Zhang, A. Kl¨ umper, and V. Popkov, Phys. Rev. B 106, 075406 (2022)

    2022

  12. [20]

    Yang and Y.-Z

    W.-L. Yang and Y.-Z. Zhang, Nucl. Phys. B744, 312 (2006)

    2006

  13. [21]

    Popkov and T

    V. Popkov and T. Prosen, arXiv:2502.06731

    arXiv

  14. [22]

    Cao, H.-Q

    J. Cao, H.-Q. Lin, K.-J. Shi, and Y. Wang, Nucl. Phys. B663, 487 (2003)

    2003

  15. [23]

    0”( always in the state|↓⟩ ⟨↓|) + internal spins “1,2, . . . , N

    R. I. Nepomechie, J. Phys. A37, 433 (2003). 1 APPENDICES This Appendix part contains three sections. In Appendix A we sketch a derivation of the Lindblad Master equation starting from brickwall qubit circuit. Appendices B and C contain some details and Tables for XXZ and XYZ m...

    2003

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.