We study the zero-temperature phase diagram of the disordered spin-
1 Bose-Hubbard model in a two-dimensional square lattice. To this aim, we use a mean-field Gutzwiller ansatz and a probabilistic mean-field perturbation theory. The spin interaction induces two different regimes, corresponding to a ferromagnetic and antiferromagnetic order. In the ferromagnetic case, the introduction of disorder reproduces analogous features of the disordered scalar Bose-Hubbard model, consisting in the formation of a Bose glass phase between Mott insulator lobes. In the antiferromagnetic regime, the phase diagram differs more from the scalar case. Disorder in the chemical potential can lead to the disappearance of Mott insulator lobes with an odd-integer filling factor and, for sufficiently strong spin coupling, to Bose glass of singlets between even-filling Mott insulator lobes. Disorder in the spinor coupling parameter results in the appearance of a Bose glass phase only between the
n and the
n+1 lobes for
n odd. Disorder in the scalar Hubbard interaction inhibits Mott insulator regions for occupation larger than a critical value.
We study quantum phase transitions of Mott insulators between spin-glass and featureless paramagnetic phases. Specifically, we consider a triangular-lattice bilayer Heisenberg model with bond disorder. The clean system has two phases, a dimer quantum paramagnet and a non-collinear antiferromagnet, separated by a quantum critical point. Bond disorder destroys the antiferromagnetic phase via the interference of dipolar textures and results in spin-glass order, such that the system features a quantum phase transition between spin-glass and dimer phases. We study the vicinity of this transition using a variant of bond-operator theory and calculate thermodynamic observables as well as excitation spectra. Bond disorder leads to strong inhomogeneities near the transition which suppresses non-collinearities and leads to anomalously weak glassiness in the near-critical quantum spin glass. We characterize the low-energy excitations which have a strong tendency towards spatial localization, and we track the behavior of the amplitude (i.e. Higgs) mode across the glass phase whose spectral weight we find to be strongly suppressed due to interference effects.
We use quantum Monte Carlo simulations to study a disordered
S=1/2 Heisenberg quantum spin model with three different nearest-neighbor interactions, ${J}_{1}\ensuremath{\le}{J}_{2}\ensuremath{\le}{J}_{3}$, on the square lattice. We consider the regime in which
J1 represents weak bonds, and
J2 and
J3 correspond to two kinds of stronger bonds (dimers) which are randomly distributed on columns forming coupled two-leg ladders. When increasing the average intradimer coupling
(J2+J3)/2, the system undergoes a N\'eel to quantum glass transition of the ground state and later a second transition into a quantum paramagnet. The quantum glass phase is of the gapless Mott glass type (i.e., in boson language it is incompressible at temperature
T=0), and we find that the temperature dependence of the uniform magnetic susceptibility follows the stretched exponential form $\ensuremath{\chi}\ensuremath{\sim}exp(\ensuremath{-}b/{T}^{\ensuremath{\alpha}})$, with $0<\ensuremath{\alpha}<1$. At the N\'eel-glass transition, we observe the standard O(3) critical exponents, which implies that the Harris criterion for the relevance of the disorder is violated in this system.