pith. sign in

arxiv: 2602.22857 · v1 · submitted 2026-02-26 · ❄️ cond-mat.quant-gas

Macroscopic quantum self-trapping in bosonic Josephson junctions: an exact quantum treatment

Andrea Bardin , Anna Minguzzi , Luca Salasnich This is my paper

Pith reviewed 2026-05-15 19:24 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas
keywords Bose-Josephson junctionmacroscopic quantum self-trappingBose-Hubbard modelpopulation imbalanceexact quantum dynamicsfinite particle numbermean-field breakdown
0
0 comments X

The pith

Exact quantum dynamics causes macroscopic quantum self-trapping to break down after a finite time for any finite particle number.

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

The paper demonstrates that in a perfectly symmetric bosonic Josephson junction modeled by the two-mode Bose-Hubbard Hamiltonian, the exact quantum evolution of the population imbalance leads to breakdown of self-trapping after finite time, regardless of particle number or initial state. This occurs because the full many-body spectrum always allows the imbalance to return to zero eventually. A sympathetic reader cares because the result shows how mean-field nonlinear phenomena like persistent trapping arise only approximately for large particle numbers and are ultimately unstable in exact quantum mechanics. The authors prove the breakdown using Hamiltonian symmetries and identify eigenvalue-difference branching that creates two dynamical regimes, including a quasi-self-trapping window for large N.

Core claim

For any finite number of particles the exact quantum dynamics leads to the breakdown of macroscopic quantum self-trapping after a finite time, regardless of the initial state. Using the symmetries of the Bose-Hubbard Hamiltonian, we provide a mathematical demonstration of this result and analyze the spectral properties governing the dynamics. We identify a branching behavior in the eigenvalues differences and a nontrivial structure of the population-imbalance amplitudes. These features allow us to distinguish two clearly different dynamical regimes and to elucidate the mechanism leading to the emergence of a quasi-MQST regime for large particle numbers.

What carries the argument

Symmetries of the two-mode Bose-Hubbard Hamiltonian together with the branching structure of eigenvalue differences that controls the time evolution of population-imbalance amplitudes.

If this is right

  • Macroscopic quantum self-trapping is only a transient or approximate phenomenon for any finite particle number and must break down after sufficient time.
  • Two distinct dynamical regimes exist, separated by the branching of eigenvalue differences, with a quasi-self-trapping regime appearing only for large particle numbers.
  • The mean-field prediction of stable self-trapping is recovered only in the limit of infinite particle number.
  • The breakdown time is set by the smallest relevant eigenvalue spacing in the many-body spectrum.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • In cold-atom experiments the trapping lifetime should grow with particle number, offering a direct test of the predicted crossover to quasi-MQST.
  • The same spectral mechanism may limit self-trapping in other finite many-body systems that exhibit mean-field bistability or nonlinear oscillations.
  • Adding small asymmetries or higher-mode couplings could shorten or lengthen the breakdown time in a controllable way.

Load-bearing premise

The two-mode Bose-Hubbard Hamiltonian perfectly describes the symmetric junction without contributions from higher modes or imperfections.

What would settle it

An exact diagonalization or time-dependent simulation of the two-mode Bose-Hubbard model for finite N showing that the population imbalance remains trapped indefinitely without returning to zero.

Figures

Figures reproduced from arXiv: 2602.22857 by Andrea Bardin, Anna Minguzzi, Luca Salasnich.

Figure 1
Figure 1. Figure 1: FIG. 1. Eigenvalues difference [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Fully quantum (left panel) and mean-field (right [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
read the original abstract

We investigate the fully quantum evolution of the population imbalance in a perfectly symmetric Bose-Josephson junction modeled by a two-mode Bose-Hubbard Hamiltonian, focusing on the validity of macroscopic quantum self-trapping beyond the mean-field theory. We show that for any finite number of particles the exact quantum dynamics leads to the breakdown of macroscopic quantum self-trapping after a finite time, regardless of the initial state. Using the symmetries of the Bose-Hubbard Hamiltonian, we provide a mathematical demonstration of this result and analyze the spectral properties governing the dynamics. We identify a branching behavior in the eigenvalues differences and a nontrivial structure of the population-imbalance amplitudes. These features allow us to distinguish two clearly different dynamical regimes and to elucidate the mechanism leading to the emergence of a quasi-MQST regime for large particle numbers. These findings bridge the gap between mean-field predictions and exact quantum dynamics and provide insight into the emergence of classical nonlinear behavior from finite quantum many-body systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript investigates the exact quantum dynamics of population imbalance in a symmetric bosonic Josephson junction modeled by the two-mode Bose-Hubbard Hamiltonian. It claims that for any finite particle number N, macroscopic quantum self-trapping (MQST) breaks down after a finite time regardless of initial state. This is demonstrated mathematically via parity symmetry under well exchange, which implies that every energy eigenstate has vanishing expectation value of the imbalance operator Jz, forcing the time-dependent expectation value to have zero time average and thus cross zero. Spectral analysis of eigenvalue differences and population-imbalance amplitudes is used to identify two dynamical regimes and explain the emergence of a long-lived quasi-MQST regime at large N.

Significance. If the central symmetry argument holds, the result is significant because it supplies a parameter-free, exact demonstration that mean-field MQST cannot persist indefinitely in any finite-N quantum system. The identification of eigenvalue branching and nontrivial amplitude structure provides a concrete mechanism for the separation of timescales between short-time trapping and eventual breakdown, directly addressing how classical nonlinear behavior emerges from finite quantum many-body systems. The absence of fitted parameters or ad-hoc assumptions strengthens the claim.

minor comments (2)
  1. In the discussion of spectral properties (around the eigenvalue-difference branching), the notation for the two distinct regimes could be clarified by explicitly labeling the short-time and long-time scales in the relevant equations or figures.
  2. A brief remark on the validity range of the two-mode approximation (e.g., a sentence referencing the conditions under which higher modes remain negligible) would help readers assess the physical applicability without altering the mathematical result.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript and for recommending acceptance. The referee's summary correctly captures our central result on the breakdown of MQST due to parity symmetry and the identification of distinct dynamical regimes at large N.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper derives the breakdown of MQST for finite N directly from the parity symmetry of the two-mode Bose-Hubbard Hamiltonian: every eigenstate satisfies <Jz> = 0, so any time-dependent <Jz>(t) has strictly zero time average and must cross zero. This is a model-internal mathematical consequence with no fitted parameters, no self-referential definitions, and no load-bearing self-citations. The eigenvalue branching and amplitude analysis explain the quasi-MQST timescale at large N but are not required for the breakdown proof itself. The derivation is self-contained against the stated Hamiltonian.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard two-mode approximation for Josephson junctions and the use of Hamiltonian symmetries for exact solution; no free parameters or invented entities are introduced.

axioms (1)
  • domain assumption The system is described by the two-mode Bose-Hubbard Hamiltonian
    Standard model for bosonic Josephson junctions invoked throughout the abstract.

pith-pipeline@v0.9.0 · 5466 in / 1088 out tokens · 26211 ms · 2026-05-15T19:24:54.000042+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [1]

    ⟨n| ˆH|n∓1⟩=⟨n∓1| ˆH|n⟩, but also centrosymmetric since⟨n∓1| ˆH|n⟩=⟨−n±1| ˆH|−n⟩and⟨n| ˆH|n⟩= ⟨−n| ˆH|−n⟩

    The Hamiltonian ˆHis not only symmetric, i.e. ⟨n| ˆH|n∓1⟩=⟨n∓1| ˆH|n⟩, but also centrosymmetric since⟨n∓1| ˆH|n⟩=⟨−n±1| ˆH|−n⟩and⟨n| ˆH|n⟩= ⟨−n| ˆH|−n⟩. As shown in Refs. [21, 22], symmetric cen- trosymmetric (N+ 1)×(N+ 1) matrices with evenN have eigenvectors that can be ordered as alternating sym- metric and antisymmetric, beginning with the symmet- ric...

    work page 2022

  2. [2]

    B. D. Josephson, Phys. Lett.1, 251 (1962)

    work page 1962

  3. [3]

    E. L. Wolf, G. B. Arnold, M. A. Gurvitch, and J. F. Zasadzinski,Josephson Junctions: History, Devices, and Applications(Pan Stanford Publishing, Singapore, 2017)

    work page 2017

  4. [4]

    T. D. Ladd, F. Jelezko, R. Laflamme, Y. Nakamura, C. Monroe, and J. L. O’Brien, Nature464, 45 (2010)

    work page 2010

  5. [5]

    Buluta, S

    I. Buluta, S. Ashhab, and F. Nori, Rep. Prog. Phys.74, 104401 (2011)

    work page 2011

  6. [6]

    M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science269, 198 (1995)

    work page 1995

  7. [7]

    K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett.75, 3969 (1995)

    work page 1995

  8. [8]

    Smerzi, S

    A. Smerzi, S. Fantoni, S. Giovannazzi, and S.R. Shenoy, Phys. Rev. Lett.79, 4950 (1997)

    work page 1997

  9. [9]

    Albiez, R

    M. Albiez, R. Gati, J. F¨ olling, S. Hunsmann, M. Cris- tiani, and M. K. Oberthaler, Phys. Rev. Lett.95, 010402 (2005)

    work page 2005

  10. [10]

    Shin, G.-B

    Y. Shin, G.-B. Jo, M. Saba, T. A. Pasquini, W. Ketterle, and D. E. Pritchard, Phys. Rev. Lett.95, 170402 (2005)

    work page 2005

  11. [11]

    S. Levy, E. Lahoud, I. Shomroni, and J. Steinhauer., Na- ture449, 579 (2007)

    work page 2007

  12. [12]

    Raghavan, A

    S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, Phys. Rev. A59, 620 (1999)

    work page 1999

  13. [13]

    Ruostekoski and D

    J. Ruostekoski and D. F. Walls, Phys. Rev. A58, R50(R) (1998)

    work page 1998

  14. [14]

    Spagnolli, G

    G. Spagnolli, G. Semeghini, L. Masi, G. Ferioli, A. Trenkwalder, S. Coop, M. Landini, L. Pezz` e, G. Mod- ugno et al., Phys. Rev. Lett.118, 230403 (2017)

    work page 2017

  15. [15]

    Wimberger, G

    S. Wimberger, G. Manganelli, A. Brollo, and L. Salas- nich, Phys. Rev. A103, 023326 (2021)

    work page 2021

  16. [16]

    Bardin, F

    A. Bardin, F. Lorenzi, and L. Salasnich, New Journal of Physics26(1), 013021 (2024)

    work page 2024

  17. [17]

    G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, Phys. Rev. A55, 4318 (1997)

    work page 1997

  18. [18]

    Ferrini, A

    G. Ferrini, A. Minguzzi, and F. W. Hekking, Phys. Rev. A78, 023606 (2008)

    work page 2008

  19. [19]

    Links, A

    J. Links, A. Foerster, A. P. Tonel, and G. Santos, Ann. Henri Poincar´ e7, 1951 (2006)

    work page 1951

  20. [20]

    Susskind and J

    L. Susskind and J. Glogower, Physics Physique Fizika1, 49 (1964)

    work page 1964

  21. [21]

    See Supplemental Material for details, which include Ref. [21–25]

    work page

  22. [22]

    A. L. Andrew, Linear Algebra and Its Applications,7, 2 (1973)

    work page 1973

  23. [23]

    Cantoni, and P

    A. Cantoni, and P. Butler, Linear Algebra and its Appli- cations13, 3 (1976)

    work page 1976

  24. [24]

    W. W. Hager,Applied Numerical Linear Algebra (Prentice-Hall International Editions, 1988)

    work page 1988

  25. [25]

    F. T. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, Phys. Rev. A6, 2211 (1972)

    work page 1972

  26. [26]

    Zhang, D

    W-M. Zhang, D. H. Feng, and R. Gilmore, Rev. Mod. Phys.62, 867 (1990)

    work page 1990

  27. [27]

    B. N. Parlett,The symmetric eigenvalue problem(Soci- ety for Industrial and Applied Mathematics, 1997). 6 Supplemental Material EIGENV ALUES AND EIGENVECTORS OF A TRIDIAGONAL SYMMETRIC CENTROSYMMETRIC MA TRIX We are interested in solving the eigenvalue problem for the following tridiagonal symmetric centrosymmetric (N+ 1)×(N+ 1) matrix, withNeven ˆH=  ...

    work page 1997