We employ the stochastic series expansion quantum Monte Carlo method to investigate the
S=1/2 antiferromagnetic Heisenberg model on a bilayer square lattice with diluted interlayer couplings. Both regular and random dilution patterns are considered. In systems with regular dilution, tuning the interlayer interaction drives a quantum phase transition from a N\'eel-ordered phase to a quantum disordered phase, consistent with the
O(3) universality class. In contrast, random dilution gives rise to a two-step transition: from the N\'eel phase to an intermediate Mott glass (MG) phase, followed by a transition to the quantum disordered phase. Within the MG phase, the uniform magnetic susceptibility exhibits a stretched-exponential temperature dependence ${\ensuremath{\chi}}_{u}\ensuremath{\sim}exp(\ensuremath{-}b/{T}^{\ensuremath{\alpha}}), 0<\ensuremath{\alpha}<1$. At the N\'eel-to-glass transition, the quenched disorder modifies the critical exponents in a manner consistent with the Harris criterion. These findings provide new insights into disorder-driven quantum phase transitions and the emergence of glassy phases in diluted bilayer quantum magnets.
We present Monte Carlo simulations of a two-dimensional bilayer quantum Heisenberg antiferromagnet with random dimer dilution. In contrast with exotic scaling scenarios found in other random quantum systems, the quantum phase transition in this system is characterized by a finite-disorder fixed point with power-law scaling. After accounting for corrections to scaling, with a leading irrelevant exponent of omega approximately 0.48, we find universal critical exponents z=1.310(6) and nu=1.16(3). We discuss the consequences of these findings and suggest new experiments.