Pith. sign in

REVIEW 3 major objections 4 minor 37 references

Integrable 3-Site, Tilted, Extended Bose-Hubbard Model with Nearest-Neighbour Interactions

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

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

Pith's one-line read A three-site Bose-Hubbard model with nearest-neighbour interactions is shown to be integrable, even with a tilting potential.

desk verdict The paper's Bethe state ansatz is algebraically false, though the integrable model and the alternative polynomial method may still be salvageable. read the letter →

arxiv 2506.21789 v3 pith:VYLCTCSF submitted 2025-06-26 cond-mat.quant-gas nlin.SI

Jon Links This is my paper
classification cond-mat.quant-gasnlin.SI MSC 17B8081R12
keywords Bose-HubbardmodelquantumintegrabilityBetheansatzthree-sitedipolarbosonso(4)Liealgebratriplewellexactsolution
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

The paper constructs an integrable 3-site extended Bose–Hubbard model whose interactions are restricted to on-site and nearest-neighbour terms, and which remains integrable when a tilting potential is added. Prior integrable 3-site models required either periodic boundary conditions or couplings that connect sites 1 and 3 directly; here the model uses only nearest-neighbour tunneling while still admitting a Bethe ansatz solution. If correct, this is the first such integrable model with a tilting term that preserves integrability, making it a useful testbed for quantum dynamics in triple-well potentials with dipolar atoms.

What carries the argument

The central object is a realization of the o(4) Lie algebra in terms of the boson operators of three sites, with e1, h1, f1 forming one o(3) subalgebra and e2, h2, f2 a second. This algebra acts on the Fock space through a symmetry-adapted basis |N, m, n} where the Hamiltonian takes a simple quadratic form in h1 plus a linear coupling term. Bethe states are generated by applying products of (e1 − u_j) to lowest-weight vectors; the Bethe equations (3.6) arise from requiring the unwanted terms to cancel, and the ODE (B.8) with the distinct-roots argument provides the completeness proof.

What would settle it

Fix a small particle number (e.g., N = 2 or N = 4) and numerically diagonalize the Hamiltonian (2.1) under the constraints (2.2)–(2.4) for E = 0 and for E values approaching zero; compare with the energies predicted by (3.6)–(3.7). If at E = 0 the Bethe ansatz fails to reproduce the full spectrum, or if the enumeration of solutions to (3.6) becomes singular, the completeness claim fails at that point.

Watch

Extended reading notes

Core claim

Under the parameter constraints (2.2)–(2.4), the general tilted 3-site extended Bose–Hubbard Hamiltonian (2.1) can be expressed in terms of two commuting o(3) subalgebras of an o(4) Lie algebra, yielding H = U $N^{2}$ + (U/4) $h1^{2}$ − (mu/2) h1 + (E/$\sqrt$(2))(e1 + f1). This Hamiltonian conserves the Casimir invariants of o(4), giving it enough conserved quantities to be integrable. The paper then develops a Bethe ansatz: eigenstates take the product form (3.1) built from the o(3) raising operator e1, with rapidities satisfying the Bethe ansatz equations (3.6), and energies given by (3.7). A completeness proof in Appendix B shows that all eigenstates and all energy levels are obtained this way, with no spurious solutions, under the assumption that the tunneling strength E is nonzero.

Load-bearing premise

The completeness proof assumes the tunneling coupling E is nonzero and that the Bethe roots are distinct, because the recursion relations dividing by E break down at E = 0 where the spectrum becomes degenerate and the one-to-one correspondence is not established.

Editorial extensions

If this is right

  • The full spectrum and all eigenstates of the model are obtained analytically via Bethe roots, enabling exact studies of dynamics in a triple well with a tilt.
  • Because the model has no spurious Bethe solutions, every solution of (3.6) corresponds to a physical eigenstate, simplifying numerical and analytic work.
  • The formulation extends naturally to 4-site systems through o(3) ⊕ o(3) ⊕ o(3), as noted in the conclusion, which may be useful for atomtronic interferometry.
  • The integrability-preserving tilt offers a controlled way to tune between integrable and chaotic regimes in dipolar bosonic systems without breaking exact solvability.

Reading between the lines

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

  • The paper leaves the E = 0 limit unexamined; the completeness proof divides by E, and at E = 0 the Bethe equations become singular, so a separate treatment would be needed to claim complete integrability in that degenerate limit.
  • The existence of an integrable tilted model suggests that careful engineering of dipolar interaction strengths might realize integrable dynamics in existing triple-well cold-atom setups, since the required couplings (2.3) are in principle tunable by shaping the trapping ellipsoids.
  • The Bethe ansatz structure could allow closed-form expressions for observables like population imbalance or entanglement entropy in the resonant-tunnelling regime, which the paper identifies as a target for future work.
  • A direct numerical check for small particle numbers comparing the Bethe energies (3.7) with exact diagonalization for E close to zero would clarify the practical scope of the completeness claim.
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

3 major / 4 minor

Summary. The manuscript proposes an integrable three-site extended Bose-Hubbard model with nearest-neighbour interactions, obtained by imposing the parameter constraints (2.2)-(2.4). The Hamiltonian is expressed through two commuting su(2) algebras inside o(4), and conserved operators are identified. A Bethe ansatz is then formulated: eigenstates are claimed to have the product form (3.1), with roots satisfying the Bethe ansatz equations (3.6) and energies given by (3.7). Section 4 asserts that the solution is complete and free of spurious solutions, with Appendix B supplying an ODE-polynomial completeness proof.

Significance. If correct, the paper would provide a notable example of an integrable open-chain extended Bose-Hubbard model with strictly nearest-neighbour couplings and an integrability-preserving tilt. The algebraic embedding into o(4) and the explicit parameter constraints are attractive, and the construction is genuinely parameter-free rather than fitted. However, the central Bethe ansatz is incorrect: the operator identities on which it relies fail already in the first nontrivial sector, so the claimed exact eigenstates and the correspondence with the Bethe equations are not established. The o(4) construction may contain useful ideas, but the main result as stated is not valid.

major comments (3)
  1. [Section 3, Eqs. (3.2)-(3.7)] The claimed Bethe ansatz identities are false. Counterexample: take n=1, N=1, U=1, mu=0, and E/sqrt(2)=1. The roots u1=1, u2=-1 satisfy the Bethe equations (3.6): for k=1 the left side is 2*1^2/(1-(-1))=1 and the right side is 1+1-1=1, and similarly for k=2. Using the paper's su(2) actions e1|m>=(2-m)|m+1>, f1|m>=m|m-1>, h1|m>=2(m-1)|m>, the state (3.1) is |Psi>=2|2>-|0>. Direct computation gives e1|Psi>=-2|1>, whereas the right-hand side of (3.2) equals -4|1>, so (3.2) is false. Moreover, with H=1+(1/4)h1^2+e1+f1, one obtains H|Psi>=4|2>-2|0>+2|1>, which is not proportional to |Psi>, although formula (3.7) gives energy 2. Thus the pair (1,-1) is a spurious solution of (3.6), directly contradicting the assertion in Section 4 that spurious solutions cannot occur.
  2. [Appendix A] The derivation of (3.2)-(3.5) is invalid. In the x^m representation, the raising operator is e1=2nx-x^2 d/dx, not multiplication by x. The polynomial corresponding to the product state (3.1) is obtained by applying the differential operators (e1-u_j) to the constant function 1, which does not produce the polynomial product over j of (x-u_j). The manipulations in Appendix A treat e1 as multiplication by x and omit the derivative contributions; this is exactly why the n=1 counterexample in Major Comment 1 violates (3.2).
  3. [Appendix B] The completeness proof does not establish the claimed correspondence for the states (3.1). The polynomial Q(x) in (B.6) is defined through coefficients alpha_j/j!, whereas the product state (3.1) involves coefficients without the factorial weightings. Consequently, roots v_j of Q(x) satisfying (3.6) do not parametrize the state (3.1) in the manner claimed; the final one-to-one statement relates (B.6) to (3.6), not (3.1) to (3.6). The counterexample in Major Comment 1 shows that this distinction is not a minor technicality. In addition, the recursion (B.3)-(B.4) divides by the tunneling coupling E, leaving the E=0 case untreated, although that is a secondary issue relative to the mismatch with (3.1).
minor comments (4)
  1. [Notation] The symbol E is used both for the tunneling coupling in (2.2)-(2.5) and for the energy eigenvalue in (3.7); this notational clash should be fixed.
  2. [Section 2] The statement that the listed operators 'realise the o(4) Lie algebra' is under-specified: the cross-commutators such as [e1,e2], [h1,e2], and [h1,h2] are not listed. They should be stated explicitly to justify the o(3) direct-sum decomposition.
  3. [Section 4] The argument that spurious solutions cannot occur 'since boson creation operators do not admit a non-trivial kernel' is not applicable, because e1 is not a creation operator; the counterexample in Major Comment 1 is a spurious solution of the Bethe equations.
  4. [Appendix B] The treatment of the U=0 case as 'diagonalisable by a canonical transformation' is asserted without proof; if U=0 is meant to be included in the parameter domain, this claim should be substantiated.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the integrable Hamiltonian and Bethe ansatz are derived self-containedly from stated operator identities and parameter constraints.

full rationale

The paper does not fit any parameter and then rename a fit as a prediction. The coupling constraints (2.2)-(2.4) are explicit choices, and the Hamiltonian is then algebraically rewritten using the o(4) generators; the claim of integrability rests on exhibiting conserved Casimir operators rather than on importing a prior uniqueness theorem. The Bethe ansatz is derived from the action formulae (3.2)-(3.5), with the calculations placed in Appendix A, and the completeness argument in Appendix B is a self-contained polynomial/ODE construction that ends by showing the roots of Q(x) satisfy the same equations (3.6). Self-citations to [21,23] are citations to standard differential-operator Bethe ansatz techniques, not to the three-site result being claimed, and the external reference [30] is only invoked for established ODE analysis. There is no step in which a stated output is equivalent by definition to an input, no fitted input masquerading as a prediction, and no load-bearing self-citation chain. Even if the algebraic identities or completeness proof contained an error, that would be a correctness issue rather than circularity.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claim rests on standard su(2) representation theory, the Fock-space realization of bosons, and a hand-chosen parameter specialization. No parameters are fitted to data and no new entities are introduced; the only questionable input is the unproven distinctness of Bethe roots at zero tunneling, treated in the domain assumption.

assumptions (4)
  • standard math Bosonic canonical commutation relations and the Fock space representation define the Hilbert space and operators.
    Section 2, Eq. (2.1) and the basis |l,m,n>.
  • ad hoc to paper The parameter constraints E12=E23=E, U1=2U2=U3=U12=U23=2U, and mu3=-mu1=mu are imposed by hand.
    Section 2, Eqs. (2.2)-(2.4), chosen to make the Hamiltonian a function of o1(3) generators and N.
  • standard math The operators e1,h1,f1 form an su(2) algebra and provide a complete symmetry-adapted basis.
    Section 3, basis construction and counting argument.
  • domain assumption The Bethe roots are distinct for every eigenstate, which is required for the completeness proof.
    Appendix B, the distinctness argument; the proof divides by the tunneling coupling E and fails at E=0.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Integrable 3-Site, Tilted, Extended Bose-Hubbard Model with Nearest-Neighbour Interactions." pith.science (2026). https://pith.science/paper/VYLCTCSF

@misc{pith2026250621789,
  author       = {Pith},
  title        = {Pith review of: Integrable 3-Site, Tilted, Extended Bose-Hubbard Model with Nearest-Neighbour Interactions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VYLCTCSF}},
  note         = {Machine review of arXiv:2506.21789}
}
read the original abstract

Extended Bose-Hubbard models have been employed in the study of cold-atom systems with dipolar interactions. It is shown that, for a certain choice of the coupling parameters, there exists an integrable extended 3-site Bose-Hubbard model with nearest-neighbour interactions. A Bethe ansatz procedure is developed to obtain expressions for the energy spectrum and eigenstates.

Figures

Figures reproduced from arXiv: 2506.21789 by the authors.

Figure 1
Figure 1. Schematic representation for a system of dipolar bosons tunneling in a tilted 3-site potential. The narrow profile of the potential at site 2 leads to attractive on-site dipole interaction, such that U2 may be tuned to be a lower value than U1 and U3. where FN = span{|l, m, n⟩ | l +m+n = N }, and each non-zero element of FN is an eigenvector of N with eigenvalue N . The dimensions of the components in the decomposit… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

37 extracted references · 35 canonical work pages

  1. [1]

    Albertini G., Bethe-ansatz type equations for the Fateev–Zamolodchikov spin model,J. Phys. A25(1992), 1799–1813

    work page 1992

  2. [2]

    and Math

    Avdeev L.V., Vladimirov A.A., Exceptional solutions to the Bethe ansatz equations,Theoret. and Math. Phys.69(1986), 1071–1079

    work page 1986

  3. [3]

    Bennett L., Isaac P.S., Links J., NOON state measurement probabilities and outcome fidelities: a Bethe ansatz approach,J. Phys. A56(2023), 505202, 27 pages

    work page 2023

  4. [4]

    Chaotic and regular dynamics in the three-site Bose-Hubbard model

    Bychek A.A., Muraev P.S., Maksimov D.N., Kolovsky A.R., Bulgakov E.N., Chaotic and regular dynamics in the three-site Bose–Hubbard model,AIP Conf. Proc.2241(2020), 020007, 7 pages, arXiv:1910.12489

    work page Pith review arXiv 2020

  5. [5]

    Bethe ansatz solution of an integrable, non-Abelian anyon chain with D(D_3) symmetry

    Campbell C.W., Dancer K.A., Isaac P.S., Links J., Bethe ansatz solution of an integrable, non-abelian anyon chain withD(D 3) symmetry,Nuclear Phys. B836(2010), 171–185, arXiv:1003.3514

    work page Pith review arXiv 2010

  6. [6]

    Castro E.R., Ch´ avez-Carlos J., Roditi I., Santos L.F., Hirsch J.G., Quantum-classical correspondence of a system of interacting bosons in a triple-well potential,Quantum5(2021), 563–574, arXiv:2105.10515

    work page Pith review arXiv 2021

  7. [7]

    From integrability to chaos: the quantum-classical correspondence in a triple well bosonic model

    Castro E.R., Wittmann W. K., Ch´ avez-Carlos J., Roditi I., Foerster A., Hirsch J.G., Quantum-classical correspondence in a triple-well bosonic model: from integrability to chaos,Phys. Rev. A109(2024), 032225, 13 pages, arXiv:2311.13189

    work page Pith review arXiv 2024

  8. [8]

    Chanda T., Barbiero L., Lewenstein M., Mark M.J., Zakrzewski J., Recent progress on quantum simulations of non-standard Bose–Hubbard models,Rep. Progr. Phys.88(2025), 044501, 27 pages, arXiv:2405.07775

    work page Pith review arXiv 2025

Show all 37 references
  1. [9]

    Franzosi R., Penna V., Chaotic behavior, collective modes, and self-trapping in the dynamics of three coupled Bose–Einstein condensates,Phys. Rev. E67(2003), 046227, 37 pages, arXiv:cond-mat/0203509

    2003 arXiv

  2. [10]

    Giri P.R., Deguchi T., Singular eigenstates in the even (odd) length Heisenberg spin chain,J. Phys. A48 (2015), 175207, 26 pages, arXiv:1411.5839

    2015 arXiv

  3. [11]

    K., Ymai L.H., Links J., Foerster A., Protocol designs for NOON states,Commun

    Gr¨ un D.S., Wittmann W. K., Ymai L.H., Links J., Foerster A., Protocol designs for NOON states,Commun. Phys.5(2022), 36, 13 pages, arXiv:2102.02944

    2022 arXiv

  4. [12]

    K., Ymai L.H., Links J., Foerster A., Tonel A.P., Integrable atomtronic interfer- ometry,Phys

    Gr¨ un D.S., Wittmann W. K., Ymai L.H., Links J., Foerster A., Tonel A.P., Integrable atomtronic interfer- ometry,Phys. Rev. Lett.129(2022), 020401, 9 pages, arXiv:2004.11987

    2022 arXiv

  5. [13]

    Kirillov A.N., Sakamoto R., Singular solutions to the Bethe ansatz equations and rigged configurations, J. Phys. A47(2014), 205207, 20 pages, arXiv:1402.0651

    2014 arXiv

  6. [14]

    Koch T., Lahaye T., Metz J., Fr¨ ohlich B., Griesmaier A., Pfau T., Stabilization of a purely dipolar quantum gas against collapse,Nature Phys.4(2008), 218–222, arXiv:0710.3643

    2008 arXiv

  7. [15]

    Kollath C., Roux G., Biroli G., L¨ auchli A.M., Statistical properties of the spectrum of the extended Bose– Hubbard model,J. Stat. Mech. Theory Exp.2010(2010), P08011, 18 pages, arXiv:1004.2203

    2010 arXiv

  8. [16]

    Kolovsky A.R., Bose–Hubbard Hamiltonian: quantum chaos approach,Internat. J. Modern Phys. B30 (2016), 1630009, 21 pages, arXiv:1507.03413

    2016 arXiv

  9. [17]

    Lett.68(2004), 632–638, arXiv:cond-mat/0403213

    Kolovsky A.R., Buchleitner A., Quantum chaos in the Bose–Hubbard model,Europhys. Lett.68(2004), 632–638, arXiv:cond-mat/0403213

    2004 arXiv

  10. [18]

    Lahaye T., Menotti C., Santos L., Lewenstein M., Pfau T., The physics of dipolar bosonic quantum gases, Rep. Progr. Phys.72(2009), 126401, 71 pages, arXiv:0905.0386

    2009 arXiv

  11. [19]

    Lahaye T., Santos L., Pfau T., Mesoscopic ensembles of polar bosons in triple-well potentials,Phys. Rev. Lett.104(2010), 170404, 8 pages, arXiv:0911.5288

    2010 arXiv

  12. [20]

    Henri Poincar´ e7 (2006), 1591–1600, arXiv:cond-mat/0605486

    Links J., Foerster A., Tonel A.P., Santos G., The two-site Bose–Hubbard model,Ann. Henri Poincar´ e7 (2006), 1591–1600, arXiv:cond-mat/0605486

    2006 arXiv

  13. [21]

    Links J., Hibberd K.E., Bethe ansatz solutions of the Bose–Hubbard dimer,SIGMA2(2006), 095, 8 pages, arXiv:nlin.SI/0612063

    2006

  14. [22]

    Links J., Moghaddam A., Zhang Y.Z., BCS model with asymmetric pair scattering: a non-Hermitian, exactly solvable Hamiltonian exhibiting generalized exclusion statistics,J. Phys. A46(2013), 305205, 18 pages, arXiv:1304.5818

    2013 arXiv

  15. [23]

    Links J., Zhao S.Y., A Bethe ansatz study of the ground state energy for the repulsive Bose–Hubbard dimer, J. Stat. Mech. Theory Exp.2009(2009), P03013, 16 pages

    2009

  16. [24]

    Physics95(1975), 139–169

    Moshinsky M., Patera J., Sharp R.T., Winternitz P., Everything you always wanted to know about SU(3)⊃ O(3),Ann. Physics95(1975), 139–169. Integrable 3-Site, Extended Bose–Hubbard Model 11

    1975

  17. [25]

    Nakerst G., Haque M., Chaos in the three-site Bose–Hubbard model: classical versus quantum,Phys. Rev. E 107(2023), 024210, 14 pages, arXiv:2203.09953

    2023 arXiv

  18. [26]

    Nepomechie R.I., Wang C., Twisting singular solutions of Bethe’s equations,J. Phys. A47(2014), 505004, 9 pages, arXiv:1409.7382

    2014 arXiv

  19. [27]

    Oelkers N., Links J., Ground-state properties of the attractive one-dimensional Bose–Hubbard models,Phys. Rev. B75(2007), 115119, 17 pages, arXiv:cond-mat/0611510

    2007 arXiv

  20. [28]

    Pan F., Li A., Wu Y., Draayer J.P., An exact solution of the homogenous trimer Bose–Hubbard model, J. Stat. Mech. Theory Exp.2023(2023), 033101, 21 pages

    2023

  21. [29]

    Rovirola M., Briongos-Merino H., Juli´ a-D ´ ıaz B., Guilleumas M., Ultracold dipolar bosons trapped in atom- tronic circuits,Phys. Rev. A109(2024), 063331, 8 pages, arXiv:2403.11620

    2024 arXiv

  22. [30]

    Simon B., Sturm oscillation and comparison theorems, in Sturm-Liouville theory, Birkh¨ auser, Basel, 2005, 29–43, arXiv:math.SP/0311049

    2005

  23. [31]

    K., Foerster A., Links J., Entangled states of dipolar bosons generated in a triple-well potential,SciPost Phys.2(2020), 003, 20 pages, arXiv:1909.04815

    Tonel A.P., Ymai L.H., Wittmann W. K., Foerster A., Links J., Entangled states of dipolar bosons generated in a triple-well potential,SciPost Phys.2(2020), 003, 20 pages, arXiv:1909.04815

    2020 arXiv

  24. [32]

    K., Castro E.R., Foerster A., Santos L.F., Interacting bosons in a triple well: preface of many-body quantum chaos,Phys

    Wittmann W. K., Castro E.R., Foerster A., Santos L.F., Interacting bosons in a triple well: preface of many-body quantum chaos,Phys. Rev. E105(2022), 034204, 15 pages, arXiv:2111.13714

    2022 arXiv

  25. [33]

    K., Castro E.R., Roditi I., Foerster A., Hirsch J.G., Subtle nuances between quantum and classical regimes,Chaos35(2025), 043108, 9 pages, arXiv:2411.07373

    Wittmann W. K., Castro E.R., Roditi I., Foerster A., Hirsch J.G., Subtle nuances between quantum and classical regimes,Chaos35(2025), 043108, 9 pages, arXiv:2411.07373

    2025 arXiv

  26. [34]

    K., Ymai L.H., Barros B.H.C., Links J., Foerster A., Controlling entanglement in a triple-well system of dipolar atoms,Phys

    Wittmann W. K., Ymai L.H., Barros B.H.C., Links J., Foerster A., Controlling entanglement in a triple-well system of dipolar atoms,Phys. Rev. A108(2023), 033313, 11 pages, arXiv:2305.09754

    2023 arXiv

  27. [35]

    K., Ymai L.H., Foerster A., Tonel A.P., J

    Wittmann W. K., Ymai L.H., Foerster A., Tonel A.P., J. L., Control of tunneling in an atomtronic switching device,Commun. Phys.1(2018), 91, 10 pages, arXiv:1710.05831

    2018 arXiv

  28. [36]

    Ymai L.H., Tonel A.P., Foerster A., Links J., Quantum integrable multi-well tunneling models,J. Phys. A 50(2017), 264001, 13 pages, arXiv:1606.00816

    2017 arXiv

  29. [37]

    Zheng M., Qiao Y., Wang Y., Cao J., Chen S., Exact solution of the Bose–Hubbard model with unidirectional hopping,Phys. Rev. Lett.132(2024), 086502, 7 pages, arXiv:2305.00439

    2024 arXiv

Pith tools

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