We investigate the ground-state properties of ultracold atoms with long-range interactions trapped in a two-leg ladder configuration in the presence of an artificial magnetic field. Using a Gross-Pitaevskii approach and a mean-field Gutzwiller variational method, we explore both the weakly interacting and strongly interacting regimes, respectively. We calculate the boundaries between the density-wave--supersolid and the Mott-insulator--superfluid phases as a function of magnetic flux and uncover regions of supersolidity. The mean-field results are confirmed by numerical simulations using a cluster mean-field approach.
We consider the Bose-Hubbard model on a two-leg ladder under an artificial magnetic field and investigate the superfluid--to--Mott insulator transition in this setting. Recently, this system has been experimentally realized [M. Atala et al., Nature Phys. 10, 588 (2014)], albeit in a parameter regime that is far from the Mott transition boundary. Depending on the strength of the magnetic field, the single-particle spectrum has either a single ground state or two degenerate ground states. The transition between these two phases is reflected in the many-particle properties. We first investigate these phases through the Bogoliubov approximation in the superfluid regime and calculate the transition boundary for weak interactions. For stronger interactions the system is expected to form a Mott insulator. We calculate the Mott transition boundary as a function of the magnetic field and interleg coupling with mean-field theory, strong-coupling expansion, and density matrix renormalization group (DMRG). Finally, using the DMRG, we investigate the particle-hole excitation gaps of this system at different filling factors and find peaks at simple fractions, indicating the possibility of correlated phases.
We investigate the ground-state phase diagram of repulsively interacting bosons on a two-leg ladder threaded by a uniform artificial magnetic flux, using the cluster Gutzwiller mean-field method. In the strong-rung-coupling regime, self-consistent calculations are performed on a
2×4 cluster. By analyzing the superfluid order parameter, leg-resolved currents, chiral current, the current ratio on adjacent legs, and the density imbalance between the two legs, we distinguish Mott-insulating from superfluid regimes and characterize the observed states as Meissner-like, vortex-like (superfluid or Mott insulating), or biased-ladder. In regions overlapping with previous DMRG studies, our results qualitatively agree with the established phase structure, demonstrating that the cluster Gutzwiller approach balances computational efficiency and physical accuracy. We then construct the first grand-canonical
t--
μ phase diagrams for this system, revealing how the magnetic flux modifies the shape, tilt, and extent of the Mott lobes. We further explore previously inaccessible regimes, including higher fillings
ρ≳1 and the intermediate interaction window
U/t∈[7.69,9.09]. Special attention is paid to
φ=π, where the effective triangular-ladder mapping becomes singular. Owing to the equivalence of
φ=π and
−π modulo
2π, a combined symmetry forbids net chiral currents, leading to a nonchiral Mott-insulating state, in contrast to the chiral-superfluid tendency expected away from
φ=π. Our results offer a computationally efficient route for mapping the global phase structure of bosonic flux ladders and provide guidance for future ultracold-atom experiments in artificial gauge fields.
In view of strongly interacting bosons in an optical lattice with a large gauge field, we study phase transitions in a two-dimensional Bose-Hubbard model with a staggered flux, on the basis of variational Monte Carlo calculations. Inthe trial states,besides typical onsite and intersite correlation factors, we introduce a configuration-dependent phase factor,which was recently found essential for treating current-carrying states. It is found that this phase factor is qualitativelyvitalfordescribing a Mott insulating (MI) state in the present model. Thereby, the Peierls phasesattached in relevant hopping processes are cancelled out. As a result, local currents completely suppressed in MI states, namely, a chiral Mott state does not appear for the square lattice, in contrast tothecorresponding two-leg ladder model. In addition, we discuss other features of the first-order superfluid-MI transition in this model.
Abstract We investigate the ground state properties of interacting bosons on a triangular ladder subject to an artificial magnetic field. By means of the mean-field approximation, rich phases are revealed in the limit of weak interaction and large bosons occupation per site, including the Meissner superfluid, vortex superfluid, biased chiral superfluid (C-SF), C-SF and insulating phases. The quantum phases can be precisely identified and the underlying physics of phase transitions can be clarified by utilizing order parameters based on the atomic density, the effective rung current, and the current asymmetry ratio. Meanwhile, the phase transition of the system can be effectively controlled by adjusting the rung-to-leg coupling ratio, the magnetic flux and the atomic interaction. Furthermore, using Bogoliubov theory, we have mapped the stability phase diagram across all parameters, and clarified the interdependence between superfluidity and the dispersion relation. The excitation spectrum, which effectively describes the system’s energetic stability, is governed by the rung-to-leg coupling ratio, the magnetic flux, and the atomic interaction. Numerical simulations validate the dynamical instability of the system. This work offers a theoretical foundation for identifying new quantum phases exhibited by ultracold atoms in optical lattices and new experimental approaches for achieving system superfluidity.
Motivated by experiments on Josephson junction arrays in a magnetic field and ultracold interacting atoms in an optical lattice in the presence of a ``synthetic'' orbital magnetic field, we study the ``fully frustrated'' Bose-Hubbard model and quantum
XY model with half a flux quantum per lattice plaquette. Using Monte Carlo simulations and the density matrix renormalization group method, we show that these kinetically frustrated boson models admit three phases at integer filling: a weakly interacting chiral superfluid phase with staggered loop currents which spontaneously break time-reversal symmetry, a conventional Mott insulator at strong coupling, and a remarkable ``chiral Mott insulator'' (CMI) with staggered loop currents sandwiched between them at intermediate correlation. We discuss how the CMI state may be viewed as an exciton condensate or a vortex supersolid, study a Jastrow variational wave function which captures its correlations, present results for the boson momentum distribution across the phase diagram, and consider various experimental implications of our phase diagram. Finally, we consider generalizations to a staggered flux Bose-Hubbard model and a two-dimensional (2D) version of the CMI in weakly coupled ladders.