We investigate the ground-state properties of ultracold atoms trapped in a two-leg ladder potential in the presence of an artificial magnetic field in a staggered configuration. We focus on the strongly interacting regime and use the Landau theory of phase transitions and a mean field Gutzwiller variational method to identify the stable superfluid phases and their boundaries with the Mott-insulator regime as a function of magnetic flux. In addition, we calculate the local and chiral currents of these superfluid phases, which show a staggered vortex-antivortex configuration. The analytical results are confirmed by numerical simulations using a cluster mean-field-theory 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 effects of artificial magnetic fields on the ground-state of the two-dimensional Bose-Hubbard model. Using an asymmetric Bose-Hubbard model, we demonstrate that anisotropic hopping amplitude localizes bosons and enlarges the insulating and supersolid regions. We show these increases in the presence of artificial gauge fields up to the symmetric point of the field. Moreover, our calculations exhibit real-space modulations of the superfluid and supersolid phases. The bosonic current exhibits vortices in these phases, whose configurations depend on the commensuration between the magnetic field and the lattice. Although the magnetic field explicitly breaks the translational symmetry of the square lattice, this symmetry is restored in the Mott insulator phase. The local densities shows a checkerboard pattern in the density-wave and supersolid phases, regardless of the magnetic field strength. We investigate thermal fluctuations and demonstrate the robustness of insulating and supersolid phases up to temperatures comparable to the interaction energy, which supports the feasibility of observing such phases in experiments.
We investigate the localization of a Bose-Einstein condensate trapped in a two-leg ladder in the presence of an artificial magnetic field via Bose-Hubbard model and discuss its abundant localized phenomena using variational approximation and direct numerical simulation. A two-leg bosonic ladder provides an ideal platform to study the complex dynamics of Bose-Einstein condensates in an optical lattice subject to an artificial magnetic field. As a main result, the system displays rich interesting localized states including self-trapping, soliton and breather. The coupling of the repulsive atomic interaction, artificial magnetic field, the rung-to-leg coupling ratio of the ladder, and the initial population difference of atoms between the two legs of the ladder results in the transitions between the localized states. The critical conditions for the transitions between the localized states are given analytically and are intuitively demonstrated in the phase diagram. Particularly, when the rung-to-leg coupling radio of the ladder satisfies a critical condition, we have the localized states with no local atomic currents between the two legs, otherwise, we observe the localized states with local atomic currents between the two legs. Finally, the results obtained via variational method are confirmed by the numerical simulations of the full discrete nonlinear Schrödinger equation describing the system.
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.
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.