We introduce a generic bosonic model exemplifying that (spin) Meissner currents can persist in insulating phases of matter. We consider two species of interacting bosons on a lattice. Our model exhibits separation of charge (total density) and spin (relative density): the charge sector is gapped in a bosonic Mott insulator phase with total density one, while the spin sector remains superfluid due to interspecies conversion. Coupling the spin sector to the gauge fields yields a spin Meissner effect reflecting the long-range spin superfluid coherence. We investigate the resulting phase diagram and describe other possible spin phases of matter in the Mott regime possessing chiral currents as well as a spin-density wave phase. The model presented here is realizable in Josephson junction arrays and in cold atom experiments.
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 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 two-component bosons on the Harper-Hofstadter model with two legs. The synthetic magnetic fields for the two types of bosons point to either the same direction or opposite directions. The bosons have hardcore intraspecies interaction such that there can be no more than one boson of the same type on each lattice site. For certain filling factors in the absence of interspecies interaction, each component realizes a vortex Mott insulator with rung current or a Meissner superfluid without rung current. The system undergoes phase transitions to other phases as interspecies interaction is turned on, which are characterized numerically using the density matrix renormalization group method and supplemented with analytical studies when possible. The vortex Mott insulator transits to a gapped Meissner phase without rung current and the Meissner superfluid transits to a gapped vortex phase with rung current. In both cases, we observe gapped spin density wave states that break certain
Z2 symmetries.
Flux ladders constitute the minimal setup enabling a systematic understanding of the rich physics of interacting particles subjected simultaneously to strong magnetic fields and a lattice potential. In this paper, the ground-state phase diagram of a flux-ladder model is mapped out using extensive density-matrix renormalization-group simulations. The emphasis is put on parameters which can be experimentally realized exploiting the internal states of potassium atoms as a synthetic dimension. The focus is on accessible observables such as the chiral current and the leg-population imbalance. Considering a particle filling of one boson per rung, we report the existence of a Mott-insulating Meissner phase as well as biased-ladder phases on top of superfluids and Mott insulators. Furthermore, we demonstrate that quantum quenches from suitably chosen initial states can be used to probe the equilibrium properties in the transient dynamics. Concretely, we consider the instantaneous turning on of hopping matrix elements along the rungs or legs in the synthetic flux-ladder model, with different initial particle distributions. We show that clear signatures of the biased-ladder phase can be observed in the transient dynamics. Moreover, the behavior of the chiral current in the transient dynamics is discussed. The results presented in this paper provide guidelines for future implementations of flux ladders in experimental setups exploiting a synthetic dimension.