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.
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.