pith. sign in

arxiv: 2604.11169 · v1 · submitted 2026-04-13 · ❄️ cond-mat.quant-gas · cond-mat.stat-mech· cond-mat.str-el· cond-mat.supr-con

Emergence of the unexpected charge-density-wave phase driven by artificial gauge field in three-leg Bose-Hubbard ladder

Takayuki Yokoyama , Yasuhiro Tada This is my paper

Pith reviewed 2026-05-10 15:49 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.stat-mechcond-mat.str-elcond-mat.supr-con
keywords Bose-Hubbard ladderartificial gauge fieldscharge density wavevortex phasesquantum phase transitionshard-core bosonsthree-leg ladder
0
0 comments X

The pith

Charge-density-wave phases emerge in three-leg bosonic ladders under artificial gauge flux even with only on-site interactions.

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

The paper studies hard-core bosons at half filling on a three-leg ladder subject to a uniform artificial gauge field. It finds that charge-density-wave order appears in flux regimes where vortex states with circulating currents are normally expected. This occurs despite the model containing only on-site repulsion, and the authors identify both connected and isolated CDW regions in the phase diagram. Increasing the gauge flux produces a reentrant sequence of transitions from CDW to vortex superfluid and back to CDW. The results are obtained by examining current patterns and correlation functions, revealing competition between density-wave order and vortex formation.

Core claim

In the three-leg Bose-Hubbard ladder at half filling under artificial gauge field, charge-density-wave phases emerge in the flux regime where vortex states are expected and realized nearby, with only on-site interactions present. Upon increasing the gauge flux, a reentrant sequence of quantum phase transitions occurs: CDW to vortex-superfluid to CDW. This reveals strong competition between vortex phases and density-wave order, with some behavior understandable from strong-coupling but an isolated CDW region not connected to such limits.

What carries the argument

Artificial gauge flux acting on the three-leg ladder geometry, distinguished from vortex states by analyzing current patterns and density-density correlation functions.

If this is right

  • Vortex phases are destabilized in favor of CDW order at specific flux values despite on-site interactions only.
  • The phase diagram contains multiple superfluid and insulating phases with reentrant behavior under increasing gauge flux.
  • Part of the CDW emergence follows from strong-coupling limits, while an isolated CDW region requires additional mechanisms tied to the ladder geometry.
  • Current patterns and correlation functions serve as reliable diagnostics to map the competition between density-wave order and vortex formation.

Where Pith is reading between the lines

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

  • The isolated CDW region may indicate a flux-induced stabilization mechanism that could appear in wider ladders or quasi-two-dimensional bosonic systems without requiring longer-range interactions.
  • Experimental tuning of artificial flux in ultracold-atom ladders could directly test whether the reentrant sequence persists at larger system sizes.
  • The results suggest that gauge fields can shift the balance toward density-wave order in bosonic systems more generally, potentially affecting proposals for fractionalized phases in related geometries.

Load-bearing premise

Numerical identification of phases via current patterns and correlation functions on finite ladders accurately distinguishes CDW from vortex states without significant finite-size or boundary effects, especially for the isolated CDW region.

What would settle it

Observation of long-range density correlations together with the predicted reentrant CDW-vortex-CD W sequence in an ultracold-atom experiment on a three-leg ladder with tunable artificial flux.

Figures

Figures reproduced from arXiv: 2604.11169 by Takayuki Yokoyama, Yasuhiro Tada.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic illustration of a three-leg ladder under an [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Schematic illustrations of current and density config [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Phase diagram of the hard-core bosonic three-leg [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Gauge flux dependence of the chiral current [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Superfluid correlation functions. [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Momentum-resolved rung-current structure factor [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Effective couplings obtained from perturbation the [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Flux dependence of the effective couplings the chiral [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Entanglement entrooy for representative parameters in different quantum phases. Panels (a)-(g) correspond to the [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
read the original abstract

We investigate hard-core bosons at half filling on a three-leg ladder under the uniform artificial gauge field. By analyzing current patterns and correlation functions, we uncover a rich quantum phase diagram containing multiple superfluid and insulating phases. In bosonic ladder systems, increasing the gauge flux typically destabilizes the Meissner phase and leads to vortex phases characterized by circulating currents. In the present system, however, we find that charge-density-wave (CDW) phases emerge precisely in such a flux regime despite the presence of only an on-site interaction, where vortex states are naturally expected and are indeed realized in nearby parameter regions. While part of this behavior can be qualitatively understood from a strong-coupling perspective, we also identify an isolated CDW region that cannot be connected to such limits. Furthermore, upon increasing the artificial gauge flux, we observe a reentrant sequence of quantum phase transitions, CDW $\to$ vortex-superfluid $\to$ CDW, revealing a strong competition between the vortex phase and the density-wave order.

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

2 major / 2 minor

Summary. The manuscript numerically investigates hard-core bosons at half filling on a three-leg Bose-Hubbard ladder under a uniform artificial gauge field. Using current patterns and correlation functions, it maps a rich phase diagram containing superfluid and insulating phases. The central finding is the emergence of charge-density-wave (CDW) phases in flux regimes where vortex states are expected, including an isolated CDW region not connected to strong-coupling limits, together with a reentrant sequence of transitions (CDW to vortex-superfluid to CDW) upon increasing the gauge flux.

Significance. If the results hold, the work is significant for demonstrating that artificial gauge fields can stabilize unexpected CDW order in a bosonic system with only on-site repulsion, contrary to the usual dominance of vortex phases in fluxed ladders. The reentrant behavior and isolated CDW lobe highlight a non-trivial competition between orders that may be realizable in ultracold-atom experiments with synthetic gauge fields. The numerical phase-diagram construction via DMRG-style methods and explicit analysis of currents and density correlations provides concrete, falsifiable predictions.

major comments (2)
  1. [§IV] §IV (Phase diagram and isolated CDW region): The central claim of an isolated CDW phase in an intermediate-flux window (disconnected from strong-coupling limits) is identified via density-density correlations and current patterns on finite open ladders. Open boundaries on three-leg systems generically produce staggered density modulations and localized currents that can be misidentified as bulk CDW order. No finite-size extrapolation of the CDW structure factor, correlation length, or order parameter with ladder length L, nor any comparison to periodic-boundary results, is presented to establish thermodynamic stability. This directly undermines the load-bearing distinction between the isolated CDW lobe and possible boundary artifacts.
  2. [§III] §III (Numerical method and phase identification): Phase boundaries and the reentrant CDW–vortex-SF–CDW sequence rest on finite-size data without reported convergence checks (bond dimension, truncation error, or L-scaling of the CDW and vortex order parameters). Because the isolated CDW region is the key unexpected feature, the absence of these controls leaves open the possibility that the reported reentrance is influenced by finite-size or boundary effects rather than intrinsic bulk physics.
minor comments (2)
  1. [Figures 2–5] Figure captions and legends should explicitly state the ladder lengths L and bond dimensions used for each panel to allow readers to assess finite-size effects.
  2. [§II] The definition of the CDW order parameter (e.g., the wave-vector and normalization of the density-density correlator) should be given in a dedicated equation rather than only in the text.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive criticism of our manuscript. The points raised regarding potential boundary artifacts and the need for explicit convergence checks are important, and we address them point by point below. We will incorporate additional finite-size analyses and numerical details in a revised version.

read point-by-point responses

  1. Referee: [§IV] §IV (Phase diagram and isolated CDW region): The central claim of an isolated CDW phase in an intermediate-flux window (disconnected from strong-coupling limits) is identified via density-density correlations and current patterns on finite open ladders. Open boundaries on three-leg systems generically produce staggered density modulations and localized currents that can be misidentified as bulk CDW order. No finite-size extrapolation of the CDW structure factor, correlation length, or order parameter with ladder length L, nor any comparison to periodic-boundary results, is presented to establish thermodynamic stability. This directly undermines the load-bearing distinction between the isolated CDW lobe and possible boundary artifacts.

    Authors: We acknowledge that open boundaries can induce local density modulations in ladder systems. Our phase identification, however, combines the density-density correlation functions (which exhibit long-range order away from the edges) with the spatial pattern of currents, which remain uniform and non-circulating in the reported CDW regions. To strengthen this, the revised manuscript will include finite-size scaling of the CDW structure factor S(π) for L = 20 to 60, showing that the peak height extrapolates to a finite value in the isolated lobe while vanishing in adjacent phases. We will also report the correlation length ξ extracted from exponential fits to the density correlations in the bulk, confirming ξ remains comparable to or larger than L/2 in the CDW phase. Periodic boundaries are technically challenging for arbitrary flux due to flux quantization on the three-leg geometry; we will add a brief discussion of this limitation and note consistency checks at commensurate fluxes where periodic data are feasible. revision: yes

  2. Referee: [§III] §III (Numerical method and phase identification): Phase boundaries and the reentrant CDW–vortex-SF–CDW sequence rest on finite-size data without reported convergence checks (bond dimension, truncation error, or L-scaling of the CDW and vortex order parameters). Because the isolated CDW region is the key unexpected feature, the absence of these controls leaves open the possibility that the reported reentrance is influenced by finite-size or boundary effects rather than intrinsic bulk physics.

    Authors: We agree that explicit convergence data are necessary. Although our DMRG calculations used bond dimensions up to χ = 400 with truncation errors kept below 10^{-8}, these details were not fully documented. In the revision we will add a dedicated subsection reporting (i) the dependence of the CDW order parameter and vortex current on χ, demonstrating convergence for χ ≥ 200, and (ii) L-scaling plots of both the central CDW amplitude (to suppress boundary effects) and the maximum vortex current for L up to 60. These plots will show that the reentrant CDW–vortex-SF–CDW sequence remains stable with increasing L, supporting that the transitions reflect bulk physics. revision: yes

Circularity Check

0 steps flagged

No circularity: numerical phase diagram from direct computation

full rationale

The paper performs a numerical study of the hard-core Bose-Hubbard model on a three-leg ladder with artificial gauge flux, identifying phases via computed current patterns and density-density correlations on finite systems. No analytical derivation chain exists that reduces predictions to inputs by construction, no parameters are fitted and then relabeled as predictions, and no load-bearing self-citations or ansatze are invoked to justify the central claims about CDW emergence. The isolated CDW region and reentrant transitions are reported as outcomes of the simulations themselves.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard Bose-Hubbard Hamiltonian for hard-core bosons plus uniform artificial gauge field on a three-leg ladder geometry; no free parameters, new axioms, or invented entities are mentioned in the abstract.

axioms (2)
  • domain assumption Hard-core boson constraint (infinite on-site repulsion) at half filling
    Stated in abstract as the model under study; required for the phase diagram.
  • domain assumption Uniform artificial gauge field threading flux through plaquettes
    Central control parameter whose effect on current patterns and correlations is analyzed.

pith-pipeline@v0.9.0 · 5495 in / 1425 out tokens · 38457 ms · 2026-05-10T15:49:58.870991+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    Gross and I

    C. Gross and I. Bloch, Quantum simulations with ultra- cold atoms in optical lattices, Science357, 995 (2017)

    work page 2017

  2. [2]

    Sch¨ afer, T

    F. Sch¨ afer, T. Fukuhara, S. Sugawa, Y. Takasu, and Y. Takahashi, Tools for quantum simulation with ultra- cold atoms in optical lattices, Nature Reviews Physics2, 411 (2020)

    work page 2020

  3. [3]

    Aidelsburger, M

    M. Aidelsburger, M. Atala, M. Lohse, J. T. Barreiro, B. Paredes, and I. Bloch, Realization of the hofstadter hamiltonian with ultracold atoms in optical lattices, Phys. Rev. Lett.111, 185301 (2013)

    work page 2013

  4. [4]

    Miyake, G

    H. Miyake, G. A. Siviloglou, C. J. Kennedy, W. C. Bur- ton, and W. Ketterle, Realizing the harper hamiltonian with laser-assisted tunneling in optical lattices, Phys. Rev. Lett.111, 185302 (2013)

    work page 2013

  5. [5]

    Giamarchi, Quantum Physics in One Dimension (Ox- ford University Press, Oxford, 2004)

    T. Giamarchi, Quantum Physics in One Dimension (Ox- ford University Press, Oxford, 2004)

    work page 2004

  6. [6]

    Dagotto and T

    E. Dagotto and T. M. Rice, Surprises on the way from one- to two-dimensional quantum magnets: The ladder materials, Science271, 618 (1996)

    work page 1996

  7. [7]

    Orignac and T

    E. Orignac and T. Giamarchi, Meissner effect in a bosonic ladder, Phys. Rev. B64, 144515 (2001)

    work page 2001

  8. [8]

    Piraud, F

    M. Piraud, F. Heidrich-Meisner, I. P. McCulloch, S. Greschner, T. Vekua, and U. Schollwock, Vortex and meissner phases of strongly interacting bosons on a two- leg ladder, Phys. Rev. B91, 140406 (2015)

    work page 2015

  9. [9]

    Di Dio, S

    M. Di Dio, S. De Palo, E. Orignac, R. Citro, and M.-L. Chiofalo, Persisting meissner state and incommensurate phases of hard-core boson ladders in a flux, Phys. Rev. 13 B92, 060506 (2015)

    work page 2015

  10. [10]

    Petrescu and K

    A. Petrescu and K. Le Hur, Bosonic mott insulator with meissner currents, Phys. Rev. Lett.111, 150601 (2013)

    work page 2013

  11. [11]

    Greschner, M

    S. Greschner, M. Piraud, F. Heidrich-Meisner, I. P. Mc- Culloch, U. Schollw¨ ock, and T. Vekua, Spontaneous in- crease of magnetic flux and chiral-current reversal in bosonic ladders: Swimming against the tide, Phys. Rev. Lett.115, 190402 (2015)

    work page 2015

  12. [12]

    Greschner, M

    S. Greschner, M. Piraud, F. Heidrich-Meisner, I. P. Mc- Culloch, U. Schollw¨ ock, and T. Vekua, Symmetry-broken states in a system of interacting bosons on a two-leg lad- der with a uniform abelian gauge field, Phys. Rev. A94, 063628 (2016)

    work page 2016

  13. [13]

    Sachdeva, M

    R. Sachdeva, M. Singh, and T. Busch, Extended bose- hubbard model for two-leg ladder systems in artificial magnetic fields, Phys. Rev. A95, 063601 (2017)

    work page 2017

  14. [14]

    S. S. Natu, Bosons with long-range interactions on two- leg ladders in artificial magnetic fields, Phys. Rev. A92, 053623 (2015)

    work page 2015

  15. [15]

    Tokuno and A

    A. Tokuno and A. Georges, Ground states of a bose–hubbard ladder in an artificial magnetic field: field- theoretical approach, New Journal of Physics16, 073005 (2014)

    work page 2014

  16. [16]

    Uchino, Analytical approach to a bosonic ladder sub- ject to a magnetic field, Phys

    S. Uchino, Analytical approach to a bosonic ladder sub- ject to a magnetic field, Phys. Rev. A93, 053629 (2016)

    work page 2016

  17. [17]

    Petrescu and K

    A. Petrescu and K. Le Hur, Chiral mott insulators, meiss- ner effect, and laughlin states in quantum ladders, Phys. Rev. B91, 054520 (2015)

    work page 2015

  18. [18]

    Orignac, R

    E. Orignac, R. Citro, M. Di Dio, and S. De Palo, Vortex lattice melting in a boson ladder in an artificial gauge field, Phys. Rev. B96, 014518 (2017)

    work page 2017

  19. [19]

    Calvanese Strinati, E

    M. Calvanese Strinati, E. Cornfeld, D. Rossini, S. Bar- barino, M. Dalmonte, R. Fazio, E. Sela, and L. Mazza, Laughlin-like states in bosonic and fermionic atomic syn- thetic ladders, Phys. Rev. X7, 021033 (2017)

    work page 2017

  20. [20]

    Calvanese Strinati, S

    M. Calvanese Strinati, S. Sahoo, K. Shtengel, and E. Sela, Pretopological fractional excitations in the two-leg flux ladder, Phys. Rev. B99, 245101 (2019)

    work page 2019

  21. [21]

    Halati and T

    C.-M. Halati and T. Giamarchi, Bose-hubbard triangu- lar ladder in an artificial gauge field, Phys. Rev. Res.5, 013126 (2023)

    work page 2023

  22. [22]

    Kolley, M

    F. Kolley, M. Piraud, I. P. McCulloch, U. Schollw¨ ock, and F. Heidrich-Meisner, Strongly interacting bosons on a three-leg ladder in the presence of a homogeneous flux, New Journal of Physics17, 092001 (2015)

    work page 2015

  23. [23]

    Barbiero, J

    L. Barbiero, J. Cabedo, M. Lewenstein, L. Tarruell, and A. Celi, Frustrated magnets without geometrical frustra- tion in bosonic flux ladders, Phys. Rev. Res.5, L042008 (2023)

    work page 2023

  24. [24]

    Zhang and S.-J

    D.-C. Zhang and S.-J. Yang, Three-leg bosonic triangular ladder in a staggered magnetic field, Physics Letters A 527, 130022 (2024)

    work page 2024

  25. [25]

    Baldelli, C

    N. Baldelli, C. R. Cabrera, S. Juli` a-Farr´ e, M. Aidels- burger, and L. Barbiero, Frustrated extended bose- hubbard model and deconfined quantum critical points with optical lattices at the antimagic wavelength, Phys. Rev. Lett.132, 153401 (2024)

    work page 2024

  26. [26]

    Dasgupta, M

    A. Dasgupta, M. Lkacki, H. Korbmacher, G. A. Dominguez-Castro, J. Zakrzewski, and L. Santos, Chi- ral phases and dynamics of dipoles in triangular optical ladders, Phys. Rev. A113, L031301 (2026)

    work page 2026

  27. [27]

    Atala, M

    M. Atala, M. Aidelsburger, M. Lohse, J. T. Barreiro, B. Paredes, and I. Bloch, Observation of chiral currents with ultracold atoms in bosonic ladders, Nature Physics 10, 588 (2014)

    work page 2014

  28. [28]

    B. K. Stuhl, H.-I. Lu, L. M. Aycock, D. Genkina, and I. B. Spielman, Visualizing edge states with an atomic bose gas in the quantum hall regime, Science349, 1514 (2015)

    work page 2015

  29. [29]

    Y. Li, H. Du, Y. Wang, J. Liang, L. Xiao, W. Yi, J. Ma, and S. Jia, Observation of frustrated chiral dynamics in an interacting triangular flux ladder, Nature Communi- cations14, 7560 (2023)

    work page 2023

  30. [30]

    Impertro, S

    A. Impertro, S. Huh, S. Karch, J. F. Wienand, I. Bloch, and M. Aidelsburger, Strongly interacting meissner phases in large bosonic flux ladders, Nature Physics21, 895 (2025)

    work page 2025

  31. [31]

    S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett.69, 2863 (1992)

    work page 1992

  32. [32]

    Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011), january 2011 Special Issue

    U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011), january 2011 Special Issue

    work page 2011

  33. [33]

    I. P. McCulloch, From density-matrix renormalization group to matrix product states, Journal of Statistical Me- chanics: Theory and Experiment2007, P10014 (2007)

    work page 2007

  34. [34]

    haldane gap

    M. Oshikawa, M. Yamanaka, and I. Affleck, Magneti- zation plateaus in spin chains: “haldane gap” for half- integer spins, Phys. Rev. Lett.78, 1984 (1997)

    work page 1984

  35. [35]

    Tada and M

    Y. Tada and M. Oshikawa, Projectification of point group symmetries with a background flux and the lieb-schultz- mattis theorem, Phys. Rev. B113, 094439 (2026)

    work page 2026

  36. [36]

    T. Sun, Y. Zheng, and S.-J. Yang, Chiral current reversal induced by a quadratic field in the three-leg magnetic lat- tice, Journal of Physics A: Mathematical and Theoretical 53, 455301 (2020)

    work page 2020

  37. [37]

    Fishman, S

    M. Fishman, S. R. White, and E. M. Stoudenmire, The ITensor Software Library for Tensor Network Calcula- tions, SciPost Phys. Codebases , 4 (2022)

    work page 2022

  38. [38]

    Calabrese and J

    P. Calabrese and J. Cardy, Entanglement entropy and quantum field theory, Journal of Statistical Mechanics: Theory and Experiment2004, P06002 (2004)

    work page 2004