Metastability, chaos and spectrum tomography for Bose-Hubbard rings and chains
Pith reviewed 2026-05-15 00:42 UTC · model grok-4.3
The pith
The many-body spectrum of finite Bose-Hubbard rings and chains reflects underlying classical phase-space structures, explaining metastability through mixed regular-chaotic dynamics.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For finite-size Bose-Hubbard rings and chains, the many-body spectrum can be tomographically related to classical phase-space structures; this mapping accounts for metastability by showing how regular islands persist amid chaotic seas and how local Bogoliubov analysis captures the initial stability of condensates, while global ergodicity is limited by the mixed dynamics.
What carries the argument
Semiclassical tomographic perspective that maps the many-body spectrum onto classical phase-space structures to distinguish regular and chaotic regions.
If this is right
- Local stability of condensates follows directly from the Bogoliubov spectrum around fixed points in phase space.
- Global metastability persists because mixed regular-chaotic dynamics prevents full exploration of phase space on accessible timescales.
- Chaos is suppressed when the system is taken to the Gross-Pitaevskii limit, recovering integrable mean-field motion.
- Spectrum tomography provides a practical diagnostic for quantum ergodicity breaking in far-from-equilibrium Bose-Hubbard experiments.
Where Pith is reading between the lines
- The same tomographic relation may generalize to higher-dimensional lattices or to driven Bose-Hubbard systems, offering a route to predict lifetimes without full quantum simulation.
- Experiments could test the prediction by preparing condensates near regular islands and measuring the rate at which population leaks into chaotic regions.
- If the mapping holds, it suggests that controlled chaos engineering could be used to stabilize or destabilize specific many-body states in atomtronic devices.
Load-bearing premise
The semiclassical mapping between the finite many-body spectrum and classical phase-space structures remains accurate enough to explain metastability.
What would settle it
Numerical diagonalization of the Bose-Hubbard Hamiltonian for moderate particle numbers and lattice sizes that yields level statistics or eigenstate projections inconsistent with the expected regular-chaotic partitioning of classical phase space.
Figures
read the original abstract
We analyze the metastability of Bose-Hubbard condensates for finite-size one-dimensional ring lattices and open chains, using a semiclassical tomographic perspective that emphasizes the relation of the many-body spectrum to the underlying classical phase-space structures. This constitutes an arena for inspection of quantum ergodicity and localization, in far-from-equilibrium scenarios of experimental interest. Both local aspects (via Bogoliubov analysis) and global aspects (by inspecting the mixed regular-chaotic dynamics) are addressed. We also clarify how chaos is diminished in the limit of the Gross-Pitaevskii equation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes the metastability of Bose-Hubbard condensates for finite-size one-dimensional ring lattices and open chains. It employs a semiclassical tomographic perspective that relates the many-body spectrum to underlying classical phase-space structures. The work addresses quantum ergodicity and localization in far-from-equilibrium scenarios of experimental interest, covering local aspects via Bogoliubov analysis and global aspects through inspection of mixed regular-chaotic dynamics. It also clarifies the reduction of chaos in the Gross-Pitaevskii equation limit.
Significance. If the semiclassical tomographic mapping is accurate for these finite-size systems, the paper offers a coherent framework for connecting quantum spectra to classical phase-space features in Bose-Hubbard models. This contributes to studies of metastability, quantum chaos, and ergodicity in ultracold-atom platforms. The dual treatment of local (Bogoliubov) and global (mixed dynamics) aspects, together with the explicit discussion of the GPE limit, strengthens the manuscript's utility for both theory and experiment.
minor comments (3)
- The abstract states the central claims clearly but does not reference specific system sizes, interaction regimes, or quantitative measures (e.g., Lyapunov exponents or participation ratios) that appear in the full text; adding one or two such anchors would help readers gauge the scope immediately.
- Notation for the tomographic reconstruction (e.g., the precise definition of the phase-space projection operator) should be introduced with an equation number in the main text rather than only in the methods section, to improve readability for readers focused on the spectrum analysis.
- Figure captions for the mixed regular-chaotic spectra should explicitly state the classical energy or filling factor at which the Poincaré sections are computed, to allow direct comparison with the quantum level statistics shown in the same panels.
Simulated Author's Rebuttal
We thank the referee for the positive summary and significance assessment of our manuscript on metastability, chaos, and spectrum tomography in Bose-Hubbard rings and chains. We appreciate the recommendation for minor revision and will incorporate improvements to enhance clarity, particularly in the discussion of the GPE limit and experimental implications.
Circularity Check
No significant circularity detected
full rationale
The paper applies standard semiclassical tomography and Bogoliubov analysis to the Bose-Hubbard model on rings and chains. The central claim relates the many-body spectrum to classical phase-space structures via established techniques without reducing any prediction to a fitted parameter or self-referential definition. No load-bearing step collapses by construction to the inputs, and the derivation remains self-contained against external semiclassical benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Semiclassical tomographic perspective accurately relates the many-body spectrum to classical phase-space structures
Reference graph
Works this paper leans on
-
[1]
MODEL HAMIL TONIAN Consider a system withL s equally-spaced sitesx j, la- beled byj. The BHH is expressed in terms of Fock operatorsa j and associated site occupation operators nj =a † jaj. Given on-site potentialV j, inter-site hopping Kji, and on-site interactionU, it takes the form: H= LsX j=1 Vja† jaj + U 2 a† ja† jajaj − X ⟨ij⟩ Kji 2 a† jai.(1) We as...
-
[2]
SEMICLASSICAL LIMIT Formally, the BHH describes coupled non-linear oscil- lators, and therefore provides an ideal arena for the study of “quantum chaos” themes [26, 42]. In a semiclassical context, one defines action-angle coordinates via aj = √njeiφj ,(6) and the Hamiltonian expressed asH=H(φ,n), becomes similar (but not identical) to that of a supercond...
-
[3]
In the contin- uum limitv k = (1/M)kwith massM= (Kℓ 2)−1
GPE LIMIT The single-particle energies and the associated veloc- ities in a tight-binding chain that has lattice constantℓ areε k =−Kcos(ℓk) andv k =Kℓsin(ℓk). In the contin- uum limitv k = (1/M)kwith massM= (Kℓ 2)−1. The dimensionless semiclassical interaction parameter of the BHH, that controls the local non-linearity and the man- ifestation of chaos, i...
-
[4]
THE ENERGY LANDSCAPE The standard textbook procedure is to transform the BHH to the momentum basis using Eq.(3). Here, we treat the chain and the ring on equal footing and transform the finite-size BHH to the appropriate orbital basis. The Hamiltonian takes the form H= X k εkb† kbk + U 2L X C[k]b † k′′′′b† k′′′bk′′ bk′ (12) = X k εknk + U L N2 − 1 2 c U 2...
-
[5]
ZERO-ORDER BOGOLIUBOV ANAL YSIS Assuming condensation in orbital “o”, and using that the number of particles is a constant of motion, we can substitute in the first term of the Hamiltonian b† ko bko =N− P b† kbk. The zero-order Bogoliubov ap- proximation is obtained if we furthermore keep only the leading terms in the interaction that transferpairsfrom th...
-
[6]
THE DISPLACED CONDENSA TE There is a one-to-one correspondence between quan- tum and the classical Bogoliubov procedure, and in prac- tice, the latter is more transparent and more convenient. Using phase-space language, the first step is to identify the SP that supports the condensate. The DNLSE is de- rived from the BHH via canonical Hamilton’s equations...
-
[7]
EXACT BOGOLIUBOV ANAL YSIS This section is rather technical and can be skipped in the first reading. In the next section, we discuss the im- plied stability regimes for rings and chains, which provide a bridge to the tomographic inspection of the spectrum. For the purpose of stability analysis, one finds the Bo- goliubov frequenciesω q of the small oscill...
-
[8]
If all the Bogoliubov frequencies are real and positive, then the SP is a stable minimum (ES)
ST ABILITY REGIMES On the basis of the previous section, let us discuss the implications with regard to the stability of a given SP. If all the Bogoliubov frequencies are real and positive, then the SP is a stable minimum (ES). If some are nega- tive (but real), it indicates that the SP is a stable elliptic point that sits on a saddle in the energy landsc...
-
[9]
CLASSICAL TRAJECTORIES Having diagnosed whether an SP is stable or not, we now zoom out and want to know whether it is located within a stable island or embedded in a chaotic sea. For this purpose, we inspect long trajectories that are launched either in the vicinity of the SP or somewhere else, away from the SP, but at the same energy. Repre- sentative t...
-
[10]
SPECTRUM TOMOGRAPHY A tomographic image of the spectrum provides insight into the underlying phase-space structure that supports the quantum eigenstates. The following figures display quantum tomographic images of the spectrum (panels of the 3rd column), which are compared with classical tomography (panels of the 1st and 2nd columns). Fig.10 is for anL s=...
-
[11]
SPECTRAL FINGERPRINTS OF MET AST ABILITY We would like to inspect quantitatively whether or how the stability of the SP is reflected in the tomography of the quantum spectrum. For demonstration, we use the following procedure. Given a tomographic (⟨n o⟩, E) spectrum, the subscript “max” is used to indicate the many-body eigenstate|E ν⟩for which⟨n o⟩ν is m...
-
[12]
QUANTUM ERGODICITY A superficial indication of ergodicity is narrow dis- persion of the expectation value⟨n o⟩within an energy window in the tomographic image of the spectrum. A more refined procedure is to look on the full distribution. Fig.16 displays theP ν(no) of eigenstates that belong to a narrow window around the energyE SP of the inspected SP. The...
-
[13]
SUMMAR Y Systems that have two degrees of freedom, such as the Bose-Hubbard trimer, are in some sense special. We can say that they feature low-dimensional chaos: they have phase-space with distinct quasi-regular and chaotic re- gions. This is because the energy surface has 2d−1=3 di- mensions, while the KAM barriers haved=2 dimensions. Ford>2, the KAM ba...
-
[14]
Leggett,Bose-Einstein condensation in the alkali gases: Some fundamental concepts, Rev
A.J. Leggett,Bose-Einstein condensation in the alkali gases: Some fundamental concepts, Rev. Mod. Phys.73, 307 (2001)
work page 2001
-
[15]
A.J. Leggett,Superfluidity, Rev. Mod. Phys.71, S318 (1999)
work page 1999
-
[16]
A.J. Leggett,The Penrose–Onsager–Yang Approach to Superconductivity and Superfluidity, chapter in ”Dia- logues Between Physics and Mathematics”, Editors: Mo- Lin Ge, Yang-Hui He (Springer 2022)
work page 2022
- [17]
- [18]
- [19]
-
[20]
S. Levy, E. Lahoud, I. Shomroni, J. Steinhauer,The ac and dc josephson effects in a bose–einstein condensate, Nature (London)449, 579 (2007)
work page 2007
-
[21]
T. Opatrny, L. Richterek and M. Opatrny,Analogies of the classical Euler top with a rotor to spin squeezing and quantum phase transitions in a generalized Lipkin- Meshkov-Glick model, Scientific Reports8, 1984 (2018)
work page 1984
-
[22]
A. Trenkwalder, G. Spagnolli, G. Semeghini, S. Coop, M. Landini, P. Castilho, L. Pezz` e, G. Modugno, M. In- guscio, A. Smerzi, M. Fattori,Quantum phase transi- tions with parity-symmetry breaking and hysteresis, Na- ture Physics12, 826 (2016)
work page 2016
- [23]
-
[24]
S. Flach and V. Fleurov,Tunnelling in the nonintegrable trimer - a step towards quantum breathers, J. Phys.: Con- dens. Matter9, 7039 (1997) 22
work page 1997
-
[25]
K. Nemoto, C.A. Holmes, G.J. Milburn, and W.J. Munro,Quantum dynamics of three coupled atomic Bose- Einstein condensates, Phys. Rev. A63, 013604 (2000)
work page 2000
-
[26]
R. Franzosi, V. Penna,Chaotic behavior, collective modes, and self-trapping in the dynamics of three cou- pled Bose-Einstein condensates, Phys. Rev. E67, 046227 (2003)
work page 2003
-
[27]
M. Johansson,Hamiltonian Hopf bifurcations in the dis- crete nonlinear Schr¨ odinger trimer: oscillatory instabil- ities, quasi-periodic solutions and a new type of self- trapping transition, J. Phys. A: Math. Gen.37, 2201-2222 (2004)
work page 2004
- [28]
-
[29]
C. Lee, T. J. Alexander, and Y.S. Kivshar,Melting of Discrete Vortices via Quantum Fluctuations, Phys. Rev. Lett.97, 180408 (2006)
work page 2006
-
[30]
Kolovsky,Semiclassical Quantization of the Bogoli- ubov Spectrum, Phys
A.R. Kolovsky,Semiclassical Quantization of the Bogoli- ubov Spectrum, Phys. Rev. Lett.99, 020401 (2007)
work page 2007
-
[31]
P. Buonsante, V. Penna, and A. Vezzani,Quantum sig- natures of the self-trapping transition in attractive lattice bosons, Phys. Rev. A82, 043615 (2010)
work page 2010
-
[32]
T.F. Viscondi, K. Furuya,Dynamics of a Bose–Einstein condensate in a symmetric triple-well trap, J. Phys. A 44, 175301 (2011)
work page 2011
- [33]
-
[34]
L. Morales-Molina, S.A. Reyes, M. Orszag,Current and entanglement in a three-site Bose-Hubbard ring, Phys. Rev. A86, 033629 (2012)
work page 2012
-
[35]
A. Gallem´ ı, M. Guilleumas, J. Martorell, R. Mayol, A. Polls, B. Juli´ a-D´ ıaz,Fragmented condensation in Bose–Hubbard trimers with tunable tunnelling, New J. Phys.17, 073014 (2015)
work page 2015
- [36]
-
[37]
A. Brattley, H. Huang, K.K. Das,Quantum Scattering States in a Nonlinear Coherent Medium, Phys. Rev. A 108, 023314 (2023) [25–27]
work page 2023
- [38]
- [39]
- [40]
-
[41]
L. Amico et al.,Roadmap on atomtronics: State of the art and perspective, AVS Quantum Science3, 039201 (2021)
work page 2021
-
[42]
Fetter,Rotating trapped Bose-Einstein condensates, Rev
A.L. Fetter,Rotating trapped Bose-Einstein condensates, Rev. Mod. Phys.81, 647 (2009)
work page 2009
- [43]
-
[44]
T.P. Meyrath, F. Schreck, J.L. Hanssen, C.-S. Chuu, M.G. Raizen,Bose-Einstein condensate in a box, Phys. Rev. A71, 041604(R) (2005)
work page 2005
-
[45]
A.L. Gaunt, T.F. Schmidutz, I. Gotlibovych, R.P. Smith, Z. Hadzibabic,Bose-Einstein Condensation of Atoms in a Uniform PotentialPhys. Rev. Lett.110, 200406 (2013)
work page 2013
-
[46]
N. Navon, A.L. Gaunt, R.P. Smith, Z. Hadzibabic,Emer- gence of a turbulent cascade in a quantum gas, Nature 539, 72 (2016)
work page 2016
- [47]
-
[48]
N. Navon, R.P. Smith, Z. Hadzibabic,Quantum gases in optical boxes, Nature Physics17, 1334 (2021)
work page 2021
-
[49]
S.J. Garratt, C. Eigen, J. Zhang, P. Turz´ ak, R. Lopes, R.P. Smith, Z. Hadzibabic, N. Navon,From single parti- cle excitations to sound waves in a box-trapped atomic Bose-Einstein condensate, Phys. Rev. A99, 021601 (2019)
work page 2019
-
[50]
A.D. Garcia-Orozco1, L. Madeira, M.A. Moreno- Armijos1, A.R. Fritsch, P.E.S. Tavares, P.C.M. Castilho, A. Cidrim, G. Roati, V.S. Bagnato,Universal dynamics of a turbulent superfluid Bose gas, Phys. Rev. A106, 023314 (2022)
work page 2022
-
[51]
M.A. Moreno-Armijos, A.R. Fritsch, A.D. Garcia- Orozco, S. Sab, G. Telles, Y. Zhu, L. Madeira, S. Nazarenko, V.I. Yukalov, V.S. Bagnato,Observation of Relaxation Stages in a Nonequilibrium Closed Quantum System: Decaying Turbulence in a Trapped Superfluid, Phys. Rev. Lett.134, 023401 (2025)
work page 2025
-
[52]
Maximilian Schmidt, Sebastian Erne, Boris Nowak, Denes Sexty, Thomas Gasenzer,Non-thermal fixed points and solitons in a one-dimensional Bose gas, New Journal of Physics,14, 075005 (2012)
work page 2012
-
[53]
Y. Zhu, B. Semisalov, G. Krstulovic, and S. Nazarenko, Direct and Inverse Cascades in Turbulent Bose-Einstein Condensates, Phys. Rev. Lett.130, 133001 (2023)
work page 2023
-
[54]
Y. Zhu, B. Semisalov, G. Krstulovic, and S. Nazarenko, Self-similar evolution of wave turbulence in Gross- Pitaevskii system, Phys. Rev. E108, 064207 (2023)
work page 2023
-
[55]
A. R. Kolovsky,Bose-Hubbard hamiltonian: Quantum chaos approach, Int. J. Mod. Phys. B30, 1630009 (2016)
work page 2016
-
[56]
L. Pausch, E.G. Carnio, A. Rodriguez, A. Buchleitner, Chaos and Ergodicity across the Energy Spectrum of In- teracting Bosons, Phys. Rev. Lett.126, 150601 (2021)
work page 2021
- [57]
-
[58]
E.J. Mueller,Superfluidity and mean-field energy loops: Hysteretic behavior in Bose-Einstein condensates, Phys. Rev. A66, 063603 (2002)
work page 2002
- [59]
-
[60]
M. Machholm, C.J. Pethick, H. Smith,Band structure, elementary excitations, and stability of a Bose-Einstein condensate in a periodic potential, Phys. Rev. A67, 053613 (2003)
work page 2003
-
[61]
O. Fialko, M.-C. Delattre, J. Brand, A.R. Kolovsky,Nu- cleation in Finite Topological Systems During Continu- ous Metastable Quantum Phase Transitions, Phys. Rev. Lett.108, 250402 (2012)
work page 2012
-
[62]
S. Baharian, G. Baym,Bose-Einstein condensates in toroidal traps: Instabilities, swallow-tail loops, and self- trapping, Phys. Rev. A87, 013619 (2013)
work page 2013
-
[63]
A. J. Lichtenberg and M. A. Lieberman,in Regular and chaotic dynamics(Springer-Verlag, 1992)
work page 1992
-
[64]
B. V. Chirikov,A universal instability of many- 23 dimensional oscillator systems, Physics Reports52, 263 (1979)
work page 1979
-
[65]
B. V. Chirikov and V. V. Vecheslavov,Theory of Fast Arnold Diffusion in Many-Frequency Systems, Journal of Statistical Physics71(1–2), 243–258 (1993)
work page 1993
-
[66]
I. C. Percival,Regular and irregular spectra, Journal of Physics B: Atomic and Molecular Physics6, L229 (1973)
work page 1973
-
[67]
M. V. Berry,Regular and irregular semiclassical wave- functions, Journal of Physics A: Mathematical and Gen- eral10, 2083 (1977)
work page 2083
-
[68]
E. B. Stechel and E. J. Heller,Quantum ergodicity and spectral chaos, Annual Review of Physical Chemistry35, 563 (1984)
work page 1984
-
[69]
Srednicki,Chaos and quantum thermalization, Phys
M. Srednicki,Chaos and quantum thermalization, Phys. Rev. E50, 888 (1994)
work page 1994
-
[70]
V.V. Flambaum, F.M. Izrailev,Statistical theory of finite Fermi systems based on the structure of chaotic eigen- states, Phys. Rev. E56, 5144 (1997)
work page 1997
- [71]
- [72]
-
[73]
A. P. Luca D’Alessio, Yariv Kafri and M. Rigol,From quantum chaos and eigenstate thermalization to statisti- cal mechanics and thermodynamics, Advances in Physics 65, 239 (2016)
work page 2016
-
[74]
F. Borgonovi, F. Mattiotti, F.M. Izrailev,Temperature of a single chaotic eigenstate, Phys. Rev. E95, 042135 (2017)
work page 2017
-
[75]
J. M. Deutsch,Eigenstate thermalization hypothesis, Re- ports on Progress in Physics81, 082001 (2018)
work page 2018
- [76]
-
[77]
Amichay Vardi, Alba Ramos, Tsampikos Kottos, Nonconventional thermal states of interacting bosonic oligomers, Phys. Rev. Research 6, 043282 (2024)
work page 2024
- [78]
-
[79]
Steven Tomsovic, Denis Ullmo,Chaos-assisted tunneling, Phys. Rev. E50, 145 (1994)
work page 1994
-
[80]
Lars Hufnagel, Roland Ketzmerick, Marc-Felix Otto, and Holger Schanz,Eigenstates Ignoring Regular and Chaotic Phase-Space Structures, Phys. Rev. Lett.89, 154101 (2002)
work page 2002
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.