We establish the phase diagram of the strongly interacting Bose-Hubbard model defined on a two-leg ladder geometry in the presence of a homogeneous flux. Our work is motivated by a recent experiment [M. Atala et al., Nat. Phys. 10, 588 (2014)], which studied the same system, in the complementary regime of weak interactions. Based on extensive density matrix renormalization group simulations and a bosonization analysis, we fully explore the parameter space spanned by filling, interleg tunneling, and flux. As a main result, we demonstrate the existence of gapless and gapped Meissner and vortex phases, with the gapped states emerging in Mott-insulating regimes. We calculate experimentally accessible observables such as chiral currents and vortex patterns.
We perform a density-matrix renormalization-group study of strongly interacting bosons on a three-leg ladder in the presence of a homogeneous flux. Focusing on one-third filling, we explore the phase diagram in dependence of the magnetic flux and the inter-leg tunneling strength. We find several phases including a Meissner phase, vortex liquids, a vortex lattice, as well as a staggered-current (SC) phase. Moreover, there are regions where the chiral current reverses its direction, both in the Meissner and in the SC phase. While the reversal in the latter case can be ascribed to spontaneous breaking of translational invariance, in the first it stems from an effective flux increase in the rung direction. Interactions are a necessary ingredient to realize either type of chiral-current reversal.
The phase diagram of a half-filled hard-core boson two-leg ladder in a flux is investigated by means of numerical simulations based on the density matrix renormalization group (DMRG) algorithm and bosonization. We calculate experimentally accessible observables such as the momentum distribution, as well as rung current, density wave, and bond-order wave correlation functions, allowing us to identify the Mott Meissner and Mott vortex states. We follow the transition from commensurate Meissner to incommensurate vortex state at increasing interchain hopping until the critical value [Piraud et al. Phys. Rev. B 91, 140406 (2015)] above which the Meissner state is stable at any flux. For flux close to $\ensuremath{\pi}$, and below the critical hopping, we observe the formation of a second incommensuration in the Mott vortex state that could be detectable in current experiments.
We consider interacting bosonic particles on a two-leg triangular ladder in the presence of an artificial gauge field. We employ density matrix renormalization group numerical simulations and analytical bosonization calculations to study the rich phase diagram of this system. We show that the interplay between the frustration induced by the triangular lattice geometry and the interactions gives rise to multiple chiral quantum phases. Phase transitions between superfluid to Mott-insulating states occur, which can have Meissner or vortex character. Furthermore, a state that explicitly breaks the symmetry between the two legs of the ladder, the biased chiral superfluid, is found for values of the flux close to $\ensuremath{\pi}$. In the regime of hard-core bosons, we show that the extension of the bond order insulator beyond the case of the fully frustrated ladder exhibits Meissner-type chiral currents. We discuss the consequences of our findings for experiments in cold atomic systems.
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.
We study the quantum phases of bosons with repulsive contact interactions on a two-leg ladder in the presence of a uniform Abelian gauge field. The model realizes many interesting states, including Meissner phases, vortex fluids, vortex lattices, charge density waves, and the biased-ladder phase. Our work focuses on the subset of these states that breaks a discrete symmetry. We use density matrix renormalization group simulations to demonstrate the existence of three vortex-lattice states at different vortex densities and we characterize the phase transitions from these phases into neighboring states. Furthermore, we provide an intuitive explanation of the chiral-current reversal effect that is tied to some of these vortex lattices. We also study a charge-density-wave state that exists at 1/4 particle filling at large interaction strengths and flux values close to half a flux quantum. By changing the system parameters, this state can transition into a completely gapped vortex-lattice Mott-insulating state. We elucidate the stability of these phases against nearest-neighbor interactions on the rungs of the ladder relevant for experimental realizations with a synthetic lattice dimension. A charge-density-wave state at 1/3 particle filling can be stabilized for flux values close to half a flux quantum and for very strong on-site interactions in the presence of strong repulsion on the rungs. Finally, we analytically describe the emergence of these phases in the low-density regime, and, in particular, we obtain the boundaries of the biased-ladder phase, i.e., the phase that features a density imbalance between the legs. We make contact with recent quantum-gas experiments that realized related models and discuss signatures of these quantum states in experimentally accessible observables.
We introduce generic bosonic models exemplifying that chiral Meissner currents can persist in insulating phases of matter. We first consider interacting bosons on a two-leg ladder. The total density sector can be gapped in a bosonic Mott insulator at odd-integer filling, while the relative density sector remains superfluid due to interchain hopping. Coupling the relative density to gauge fields yields a pseudospin Meissner effect. We show that the same phase arises if the bosons are replaced by spinful fermions confined in Cooper pairs, and find a dual fermionic Mott insulator with spinon currents. We prove that, by tuning the mean density, the Mott insulator with Meissner currents turns into a low-dimensional bosonic $\ensuremath{\nu}=\frac{1}{2}$ Laughlin state for strong enough repulsive interactions across the ladder rungs. We finally discuss extensions to multileg ladders and bilayers in which spinon superfluids with Meissner currents become possible. We propose two experimental realizations, one with ultracold atoms in the setup of Atala et al. [Nat. Phys. 10, 588 (2014)] and another with Josephson junction arrays. We also address a Bose-Fermi mixture subject to a magnetic field in connection with the pseudogap phase of high-
Tc cuprates.