Equilibrium phases of dipolar lattice bosons in the presence of random diagonal disorder

Ultracold gases offer an unprecedented opportunity to engineer disorder and interactions in a controlled manner. In an effort to understand the interplay between disorder, dipolar interactions, and quantum degeneracy, we study two-dimensional hard-core dipolar lattice bosons in the presence of on-site bound disorder. Our results are based on large-scale path-integral quantum Monte Carlo simulations by the worm algorithm. We study the ground-state phase diagram at a fixed half-integer filling factor for which the clean system is either a superfluid at a lower dipolar interaction strength or a checkerboard solid at a larger dipolar interaction strength. We find that, even for weak dipolar interactions, superfluidity is destroyed in favor of a Bose glass at a relatively low disorder strength. Interestingly, in the presence of disorder, superfluidity persists for values of the dipolar interaction strength for which the clean system is a checkerboard solid. At a fixed disorder strength, as the dipolar interaction is increased, superfluidity is destroyed in favor of a Bose glass. As the interaction is further increased, the system eventually develops extended checkerboard patterns in the density distribution. Due to the presence of disorder, though, grain boundaries and defects, responsible for a finite residual compressibility, are present in the density distribution. Finally, we study the robustness of the superfluid phase against thermal fluctuations.