REVIEW 3 major objections 5 minor 41 references
Hierarchy of localized many-body bound states in an interacting open lattice
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper claims to identify the mechanism: boundary-localized many-body bound states arise from asymmetric string solutions of the Bethe-ansatz equations, and the same solutions predict a hierarchy of such states that is always present…
desk verdict Small-N boundary string solutions are new and numerically confirmed, but the all-N hierarchy rests on a singular L→∞ limit that is deferred to a missing supplement. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the asymmetric boundary string solution of the Bethe-ansatz equations under open boundary conditions: a set of rapidities that sits asymmetrically on the imaginary axis (for $N=3$, $i\frac{3}{2}\beta$, $i\frac{1}{2}\beta$, $-i\frac{1}{2}\beta$) rather than in complex-conjugate pairs, together with the $\gamma_j=e^{ik_j}$ values $\gamma_1=-U$, $\gamma_2=\infty$, $\gamma_3=0$ in the infinite-lattice limit. These solutions carry the argument because their wavefunctions factor into a hard nearest-neighbor constraint $\delta_{x_i,x_{i+1}-1}$ times exponential factors, so the imbalance of powers on the two sides of the delta immediately tells whether the state is localized and on which side. The recurrence $g(x)=-1/x-U$ (with companion $f(x)=-x-U$) generates the quasimomenta $\eta_j$ as particles are added one by one, and the energy formula $E_{lb,N,i}=-\sum_{l=1}^{i-1}(1/\eta_l+\eta_l)-\sum_{l=1}^{N-i-1}(1/\eta_l+\eta_l)$ closes the hierarchy.
What would settle it
Solve the finite-lattice Bethe-ansatz equations numerically for several values of $U<-2$ and $N=3,4,5$, track the roots as $L$ grows, and check whether they converge to $\gamma_2=\infty$, $\gamma_3=0$ while the corresponding eigenstates have size-independent inverse participation ratio and the predicted energy $1/U+2U$; any divergence of the norm or failure of size-independent localization would falsify the claim.
Extended reading notes
Core claim
The central discovery is that the open-boundary Bethe-ansatz equations of the spinless-fermion chain with nearest-neighbor interaction $U$ admit, in the $L\to\infty$ limit, solutions with total quasimomentum $e^{iK}=U$ (or $1/U$), giving the rapidity set $\gamma_1=-U$, $\gamma_2=\infty$, $\gamma_3=0$ for $N=3$. These 'boundary string solutions' have rapidities $\lambda_1=i\frac{3}{2}\beta$, $\lambda_2=i\frac{1}{2}\beta$, $\lambda_3=-i\frac{1}{2}\beta$, which are asymmetric about the real axis and have no counterpart under periodic boundary conditions. Their wavefunctions are $\delta_{x_1,x_2-1}e^{-ik_1x_3}$ and the mirror image $\delta_{x_2,x_3-1}e^{ik_1x_1}$: one pair of fermions is forced onto neighboring sites while the third decays exponentially away from the edge. For larger $N$, the quasimomenta obey the recurrence $\eta_{j+1}=-1/\eta_j-U$, the wavefunctions take the form of powers of these $\eta$'s multiplied by a single nearest-neighbor delta function, and the energy is a sum of terms $-(1/\eta_l+\eta_l)$. Counting the states gives $2\lfloor (N-1)/2\rfloor$ localized $N$-body bound states, and the recurrence proves they always exist for $U<-2$ and $N\ge 3$, with some parameter regimes placing their energies inside the continuum, i.e., bound states in the continuum.
Load-bearing premise
The derivation assumes that the infinite-lattice limit of the Bethe equations, in which two rapidities become zero and infinity, is a legitimate solution with a normalizable wavefunction; the paper defers the proof of this limit and of the boundary conditions to the Supplemental Material, and if that limit is invalid the hierarchy collapses.
Editorial extensions
If this is right
- For any attractive interaction stronger than $U=-2$, an open chain with $N\ge 3$ fermions always possesses at least one localized many-body bound state, regardless of the lattice length in the thermodynamic limit.
- The number of such states grows with particle number: there are $2\lfloor (N-1)/2\rfloor$ localized $N$-body bound states, and for $N\ge 5$ some of them are localized away from the edge rather than at it.
- In suitable parameter regions these localized states survive inside the continuous spectrum, so the system provides exact realizations of many-body bound states in the continuum.
- Because the fermion model maps by a Jordan-Wigner transformation to the XXZ spin chain, the same hierarchy of boundary-localized states appears in the spin chain.
- The recurrence $\eta_{j+1}=-1/\eta_j-U$ fixes the quasimomentum of every member of the hierarchy once the interaction strength is given, so the full family of localized states is determined by a single parameter.
Reading between the lines
- A natural extension, only implicit in the paper, is that similar asymmetric boundary string solutions may exist in other integrable open-boundary models, such as the Hubbard chain or hard-core Bose gases, and they could be searched for directly by looking for Bethe roots with complex total quasimomentum.
- The left/right independence noted for non-overlapping clusters suggests that for $N\ge 6$ the count of localized states may grow beyond $2\lfloor (N-1)/2\rfloor$ by combining a left-localized cluster with a right-localized cluster; this combinatorial structure is not enumerated in the paper.
- The exponentially decaying density tails and the pinned nearest-neighbor pair are concrete experimental signatures: quantum-gas microscopes that resolve individual sites could detect the fixed-distance pair and the exponential profile after a quench.
- The bound-states-in-continuum regime could be probed dynamically: because the localized states have no overlap with extended plane waves, preparing one should show persistent trapping while all other excitations radiate away.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies spinless fermions on an open one-dimensional lattice with nearest-neighbor interaction, equivalently an open XXZ chain, and claims to solve the Bethe-ansatz equations in the L→∞ limit. It identifies asymmetric string solutions, corresponding to boundary-localized many-body bound states, and proposes explicit wavefunctions and energies for N=3,4,5. It then states recurrence relations for general N, claims a hierarchy of 2×[(N−1)/2] localized N-body bound states, and concludes that localized bound states are always present for U<−2 and N≥3. Exact-diagonalization data are presented for N=3,4,5, including IPR scaling and density profiles.
Significance. If the existence proof were complete, the paper would be a valuable analytical mechanism for boundary-localized many-body bound states in an integrable open lattice, with parameter-free wavefunctions, energies, and a hierarchy, and it would explain earlier numerical observations. The strengths are the exact-diagonalization benchmarks for N=3,4,5, the absence of fitted parameters, the simple recurrence structure, and the IPR scaling evidence for localization. The central claim, however, currently rests on an unproved and singular L→∞ limit of the Bethe-ansatz equations, and the general-N hierarchy is supported mainly by pattern extrapolation; the paper therefore reads as a well-motivated and plausible conjecture rather than a fully established theorem.
major comments (3)
- [Model and solutions, Eqs. (3)-(8)] The existence of the asymmetric string solution is not established by the argument given. The text states that solving Eqs. (3) in the limit L→∞ gives e^{iK}=U and hence Eq. (5), but Eq. (3) is singular at γ=0 and γ=∞ because the factors 1+e^{i(kj+ki)}+Ue^{ik_i} on the right-hand side diverge, and e^{-i2kj(L+1)} has no well-defined limit for complex kj. No finite-L sequence of solutions of Eq. (3) approaching the singular configuration is exhibited. Since the wavefunctions (7)-(8) and all subsequent recurrence states are built on Eq. (5), the main existence claim is conditional unless a finite-L construction or an independent direct derivation of (7)-(8) as eigenstates of the open-chain Hamiltonian is supplied.
- [Model and solutions and finite-size checks, Eqs. (7)-(15)] The displayed delta-pinned wavefunctions are not shown to be exact eigenstates of the finite-L Hamiltonian (1). For u_{3,L,1}, the equation at the right boundary x3=L omits the hopping term to x3=L+1, so the boundary condition is satisfied only up to corrections of order |γ1|^{-L}; direct substitution also leaves the status of the other open boundary unclear. The paper should state precisely in which limit (semi-infinite chain, L→∞ with exponential convergence) these functions are eigenstates, and it should quantify the finite-size corrections that are visible, for example, in the IPR and energy comparisons in Figs. 1 and 2.
- [Hierarchy of localized many-body bound states, Eqs. (16)-(17) and final paragraph] The induction argument for the claim that localized bound states are always present for U<−2 and N≥3 assumes that the recurrence-generated wavefunctions are eigenstates, but this is not derived from Eq. (3) or Eq. (2) for general N. The displayed inequalities verify localization and normalizability conditions only; they do not prove that u_{N,L,i} and u_{N,R,i} solve the Schrödinger equation. The stated count 2×[(N−1)/2] also needs derivation and appears to exclude the composite left/right states mentioned for N≥6. Please provide a general proof, or clearly label the hierarchy as a conjecture supported by the N=3,4,5 numerics.
minor comments (5)
- [Throughout] There are several typographical errors: 'do not exit under PBC' should be 'do not exist under PBC'; reference [2] is mistranscribed; and reference [10] has 'specific hear' for 'specific heat'.
- [Eqs. (16)-(17)] The notation for the wavefunctions is hard to read because the subscripts and superscripts of η are not clearly distinguished; please typeset expressions such as η_2^{x_{i-2}} η_1^{x_{i-1}} and define the ranges of the indices explicitly.
- [General energy formula after Eq. (17)] The general energy formula appears to omit the interaction energy U associated with the pinned bond: for N=3,i=1 it gives U+1/U instead of the value 2U+1/U quoted in the text. Please correct the formula and check it against the N=4 and N=5 energies.
- [Counting statement] The bracket [(N−1)/2] should be replaced by an explicit floor function, and the sentence 'there are 2×[(N−1)/2] localized N-body bound states' should be reconciled with the additional left/right composite states mentioned for N≥6.
- [Supplemental Material] The derivation of the asymmetric string and of the general hierarchy is deferred to a Supplemental Material that is not available in the arXiv version; the main text should either include the key steps or make clear that the full derivation will be available with the published version.
Circularity Check
No significant circularity: the derivation is analytical from the Bethe-ansatz equations, and the hierarchy follows from recurrence relations plus an induction on localization inequalities, with no fitted parameters.
full rationale
The paper derives the boundary-localized bound states from the OBC Bethe-ansatz equations (3): the asymmetric string solution (5) with gamma1 = -U, gamma2 = infinity, gamma3 = 0 is presented as the L-infinity solution, the wavefunctions (7)-(8) and energy Elb,3,1 = 2U + 1/U follow from that solution, and the numerical exact-diagonalization results in Fig. 1 provide independent confirmation. For N >= 4 the solutions are obtained by solving the same BA equations (claimed for N = 4, 5) and organized by the recurrence eta_{j+1} = -1/eta_j - U (Fig. 3), with closed-form wavefunctions (16)-(17). The central claim, 'localized bound states are always present for U < -2 and N >= 3,' is supported by an induction that uses the recurrence relation and the localization/normalizability inequalities derived from the wavefunction form; it does not assume the target result as a premise. No parameter is fitted to data and then renamed a prediction. The two self-citations ([13] Hao-Zhang-Liang-Chen and [41] Liu-Chen) are contextual references (bound-state examples and BIC literature, respectively) and carry no load-bearing weight. No uniqueness theorem is imported from the authors' prior work. The main weakness - the treatment of gamma2 = infinity and gamma3 = 0 as legitimate limits of the finite-L BA equations, with the derivation deferred to the Supplemental Material - is a mathematical-rigor concern about whether the claimed solutions exist, not a circular reduction, since the paper never defines the boundary string in terms of the hierarchy it claims to predict.
Assumptions & free parameters
assumptions (4)
- domain assumption The Bethe-ansatz equations (3) correctly give the spectrum of the open-chain spinless fermion model (1).
- ad hoc to paper In the L→∞ limit, solutions of the BA equations with γ2=∞ and γ3=0 (and more generally, complex quasimomenta with unbounded imaginary parts) are legitimate and yield normalizable boundary-localized wavefunctions.
- domain assumption The recurrence maps g(x)=-1/x-U and f(x)=-x-U generate valid BA solutions for all particle numbers N.
- domain assumption Existence of localized solutions on left and right sides can be combined independently for N>6.
Cite this review
Pith. "Pith review of Hierarchy of localized many-body bound states in an interacting open lattice." pith.science (2026). https://pith.science/paper/V3KAWHKS
@misc{pith2026250523255,
author = {Pith},
title = {Pith review of: Hierarchy of localized many-body bound states in an interacting open lattice},
year = {2026},
howpublished = {\url{https://pith.science/paper/V3KAWHKS}},
note = {Machine review of arXiv:2505.23255}
}
read the original abstract
We unveil the mechanism for the formation of puzzled boundary-localized bound states in a spinless fermionic open lattice with nearest-neighbor interactions. By solving the Bethe-ansatz equation analytically, we uncover asymmetrical string solutions corresponding to the boundary-localized bound states, which emerge in systems with at least three particles. The localized bound states can become bound states in continuum in a suitable parameter region. When the number of particles increases to five or more, additional bound states away from the edge are also observed. Through rigorous analysis, we derive recurrence relations of the quasi-momentum of the localized states as a function of the number of particles, predicting the presence of hierarchy of localized many-body bound states in interacting open lattices.
Figures
Reference graph
Works this paper leans on
-
[1]
E. H. Lieb and F. Y. Wu, Absence of Mott transition in an exact solution of the short-range, one-band model in one dimension, Phys. Rev. Lett. 20, 1445 (1968)
work page 1968
-
[2]
H. A. Bethe, Entanglement in the Quantum Phase Tran- sition of the Half-Integer Spin One-Dimensional Heisen- berg Model, Z. Physik 71, 205 (1931)
work page 1931
-
[3]
L. A. Takhtadzhan and L. D. Faddeev, The Quantum method of the inverse problem and the Heisenberg XYZ model, Rush. Math. Surveys 34, 13 (1979)
work page 1979
-
[4]
Andrei, Diagonalization of the Kondo Hamiltonian, Phys
N. Andrei, Diagonalization of the Kondo Hamiltonian, Phys. Rev. Lett. 45, 379 (1980)
work page 1980
-
[5]
P. B. Wiegmann, Exact solution of s-d exchange model at T = 0, Sov.Phys. JETP Lett. 31, 392 (1980)
work page 1980
-
[6]
M. Takahashi, Thermodynamics of One-dimensional Solvable Models (Cambridge University Press, Cam- bridge, England, 1999)
work page 1999
-
[7]
Batchelor, and Chaohong Lee, Fermi gases in one dimension: From Bethe ansatz to ex- periments, Rev
Xi-Wen Guan, Murray T. Batchelor, and Chaohong Lee, Fermi gases in one dimension: From Bethe ansatz to ex- periments, Rev. Mod. Phys. 85, 1633 (2013)
work page 2013
-
[8]
Y.-P. Wang, W.-L. Yang, J. Cao, and K. Shi, Off- Diagonal Bethe Ansatz for Exactly Solvable Model (Springer, Berlin, Heidelberg, 2015)
work page 2015
Show all 41 references
-
[9]
Takahashi, M
M. Takahashi, M. Suzuki, One-dimensional anisotropic Heisenberg model at finite temperatures, Prog. Theor. Phys. 48, 2187 (1972)
1972
-
[10]
Takahashi, Low-temperature specific hear of spin 1/2 anisotropic Heisenberg ring, Prog
M. Takahashi, Low-temperature specific hear of spin 1/2 anisotropic Heisenberg ring, Prog. Theor. Phys. 50, 1519 (1973)
1973
-
[11]
E. H. Lieb and W. Liniger, Exact analysis of an inter- acting Bose gas. I. The general solution and the ground state, Phys. Rev. 130, 1605 (1963)
1963
-
[12]
J. B. McGuire, Study of exactly soluble one-dimensional N-body problems, J. Math. Phys. 5, 622 (1964)
1964
-
[13]
Y. Hao, Y. Zhang, J.-Q. Liang, and S. Chen, Ground- state properties of one-dimensional ultracold Bose gases in a hard-wall trap, Phys. Rev. A 73, 063617 (2006)
2006
-
[14]
R. Koch, A. Bastianello, and J.-S. Caux, Adiabatic for- mation of bound states in the one-dimensional Bose gas, Phys. Rev. B 103, 165121 (2021)
2021
-
[15]
W. Yang, J. Wu, S. Xu, Z. Wang, and C. Wu, One- dimensional quantum spin dynamics of Bethe string states, Phys. Rev. B 100, 184406 (2019)
2019
-
[16]
Z. Wang, J. Wu, W. Yang, A. K. Bera, D. Kamenskyi, A. T. M. N. Islam, S. Xu, J. M. Law, B. Lake, C. Wu, and A. Loidl, Experimental observation of Bethe strings, Nature 554, 219 (2018)
2018
-
[17]
A. K. Bera, J. Wu, W. Yang, R. Bewley, M. Boehm, J. Xu, M. Bartkowiak, O. Prokhnenko, B. Klemke, A. T. M. N. Islam, J. M. Law, Z. Wang, and B. Lake, Disper- sions of many-body Bethe strings, Nature Physics 16, 625 (2020)
2020
-
[18]
Horvath, A
M. Horvath, A. Bastianello, S. Dhar, R. Koch, Y. Guo, J.-S. Caux, M. Landini, H.-C. N¨ agerl, Observing Bethe strings in an attractive Bose gas far from equilibrium, arXiv:2505.10550
-
[19]
Kardar, G
M. Kardar, G. Parisi, and Y.-C. Zhang, Dynamic scaling of growing interfaces, Phys. Rev. Lett. 56, 889 (1986)
1986
-
[20]
Ljubotina, M
M. Ljubotina, M. Znidaric, and T. Prosen, Spin diffusion from an inhomogeneous quench in an integrable system, Nature Communications 8, 16117 (2017)
2017
-
[21]
Gopalakrishnan and R
S. Gopalakrishnan and R. Vasseur, Anomalous transport from hot quasiparticles in interacting spin chains, Re- ports on Progress in Physics 86, 036502 (2023)
2023
-
[22]
Ilievski, J
E. Ilievski, J. De Nardis, S. Gopalakrishnan, R. Vasseur, 6 and B. Ware, Superuniversality of superdiffusion, Phys. Rev. X 11, 031023 (2021)
2021
-
[23]
P. N. Jepsen, J. Amato-Grill, I. Dimitrova, W. W. Ho, E. Demler, and W. Ketterle, Spin transport in a tunable heisenberg model realized with ultracold atoms, Nature 588, 403 (2020)
2020
-
[24]
Scheie, N
A. Scheie, N. E. Sherman, M. Dupont, S. E. Nagler, M. B. Stone, G. E. Granroth, J. E. Moore, and D. A. Tennant, Detection of Kardar–Parisi–Zhang hydrodynamics in a quantum Heisenberg spin-1/2 chain, Nature Physics 17, 726 (2021)
2021
-
[25]
D. Wei, A. Rubio-Abadal, B. Ye, F. Machado, J. Kemp, K. Srakaew, S. Hollerith, J. Rui, S. Gopalakrishnan, N. Y. Yao, I. Bloch, and J. Zeiher, Quantum gas microscopy of Kardar-Parisi-Zhang superdiffusion, Science 376, 716 (2022)
2022
-
[26]
Rosenberg, T
E. Rosenberg, T. I. Andersen, R. Samajdar, et al., Dy- namics of magnetization at infinite temperature in a Heisenberg spin chain, Science 384, 48 (2024)
2024
-
[27]
R. A. Pinto, M. Haque, and S. Flach, Edge-localized states in quantum one-dimensional lattices, Phys. Rev. A 79, 052118 (2009)
2009
-
[28]
N. Sun, W. Zhang, H. Yuan and X. Zhang, Boundary- localized many-body bound states in the continuum, Communications Physics, 7, 299 (2024)
2024
-
[29]
C. N. Yang and C. P. Yang, One-Dimensional Chain of Anisotropic Spin-Spin Interactions. I. Proof of Bethe’s Hypothesis for Ground State in a Finite System, Phys. Rev. 150, 321 (1966)
1966
-
[30]
C. N. Yang and C. P. Yang, One-Dimensional Chain of Anisotropic Spin-Spin Interactions. II. Properties of the Ground-State Energy Per Lattice Site for an Infinite Sys- tem, Phys. Rev. 150, 327 (1966)
1966
-
[31]
C. N. Yang and C. P. Yang, One-Dimensional Chain of Anisotropic Spin-Spin Interactions. III. Applications, Phys. Rev. 151, 258 (1966)
1966
-
[32]
Takahashi, Half-filled Hubbard model at low temper- ature, J
M. Takahashi, Half-filled Hubbard model at low temper- ature, J. Phys. C 10, 1289 (1977)
1977
-
[33]
Jordan and E
P. Jordan and E. Wigner, ¨Uber das Paulische ¨Aquivalenzverbot, Z. Phys. 47, 631 (1928)
1928
-
[34]
(II) Details for the two- particle case
See Supplemental Material for details on: (I) The deriva- tion of Bethe-ansatz equations. (II) Details for the two- particle case. (III) Energies, string solutions and wave- functions of localized bound states for systems with N = 3 to 5
-
[35]
F. C. Alcaraz, M. N. Barber, M. T. Batchelor, R. J. Baxter and G. R. W. Quispel, Surface exponents of the quantum XXZ, Ashkin-Teller and Potts models, J. Phys. A: Math. Gen. 20 6397(1987)
1987
-
[36]
Skorik and H
S. Skorik and H. Saleur, Boundary bound states and boundary bootstrap in the sine-Gordon model with Dirichlet boundary conditions, J. Phys. A: Math. Gen. 28 6605 (1995)
1995
-
[37]
Huang, Y
B. Huang, Y. Ke, H. Zhong, Y. S. Kivshar, C. Lee, Interaction-induced multiparticle bound states in the continuum, Phys. Rev. Lett. 133, 140202 (2024)
2024
-
[38]
Sugimoto, Y
S. Sugimoto, Y. Ashida, and M. Ueda, Many-body bound states in the continuum, arXiv:2307.05456 (2023)
2023 arXiv
-
[39]
J. M. Zhang, D. Braak, and M. Kollar, Bound states in the continuum realized in the one-dimensional two- particle Hubbard model with an impurity, Phys. Rev. Lett. 109, 116405 (2012)
2012
-
[40]
J. M. Zhang, D. Braak, and M. Kollar, Bound states in the one-dimensional two-particle Hubbard model with an impurity, Phys. Rev. A 87, 023613 (2013)
2013
-
[41]
Liu and S
Y. Liu and S. Chen, Fate of Two-Particle Bound States in the Continuum in Non-Hermitian Systems, Phys. Rev. Lett. 133, 193001(2024)
2024
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.