The effect of quenched disorder on the low-energy and low-temperature properties of various two- and three-dimensional Heisenberg models is studied by a numerical strong disorder renormalization-group method. For strong enough disorder we have identified two relevant fixed points, in which the gap exponent, $\ensuremath{\omega},$ describing the low-energy tail of the gap distribution $P(\ensuremath{\Delta})\ensuremath{\sim}{\ensuremath{\Delta}}^{\ensuremath{\omega}}$ is independent of disorder, the strength of couplings, and the value of the spin. The dynamical behavior of nonfrustrated random antiferromagnetic models is controlled by a singletlike fixed point, whereas for frustrated models the fixed point corresponds to a large spin formation and the gap exponent is given by $\ensuremath{\omega}\ensuremath{\approx}0.$ Another type of universality class is observed at quantum critical points and in dimerized phases but no infinite randomness behavior is found, in contrast to that of one-dimensional models.
We use a tensor network renormalization group method to study random
S=2 antiferromagnetic Heisenberg chains with alternating bond strength distributions. In the absence of randomness, an imposed dimerization with bond alternation induces two quantum critical points between the
S=2 Haldane phase, a partially dimerized phase, and a fully dimerized phase, depending on the strength of dimerization. These three phases, called ($\ensuremath{\sigma},4\ensuremath{-}\ensuremath{\sigma}$) = (2,2), (3,1), and (4,0) phases, are valence-bond solid (VBS) states characterized by $\ensuremath{\sigma}$ valence bonds (effective spin-1/2 singlet pairs) forming across even links and $4\ensuremath{-}\ensuremath{\sigma}$ valence bonds on odd links. Here we study the effects of bond randomness on the ground states of the dimerized spin chain, calculating disorder-averaged twist order parameters and spin correlations. We classify the types of random VBS phases depending on the strength of bond randomness
R and dimerization
D using the twist order parameter, which has a negative/positive sign for a VBS phase with odd/even $\ensuremath{\sigma}$. Our results demonstrate the existence of a multicritical point in the intermediate disorder regime with finite dimerization, where (2,2), (3,1), and (4,0) phases meet. This multicritical point is at the junction of three phase boundaries in the $R\text{\ensuremath{-}}D$ plane: the (2,2)-(3,1) and (3,1)-(4,0) boundaries that extend to zero randomness, and the (2,2)-(4,0) phase boundary that connects another multicritical point in the undimerized limit. The undimerized multicritical point separates a gapless Haldane phase and an infinite-randomness critical line with the diverging dynamic critical exponent in the large-
R limit at
D=0. Furthermore, we identify the (3,1)-(4,0) phase boundary as an infinite-randomness critical line even at small
R and find the signature of infinite randomness at the (2,2)-(3,1) phase boundary only in the vicinity of the multicritical point.
The effect of quenched disorder on the low-energy properties of various antiferromagnetic spin-ladder models is studied by a numerical strong disorder renormalization-group method and by density-matrix renormalization. For strong enough disorder the originally gapped phases with finite topological or dimer order become gapless. In these quantum Griffiths phases the scaling of the energy, as well as the singularities in the dynamical quantities, are characterized by a finite dynamical exponent, z, which varies with the strength of disorder. At the phase boundaries, separating topologically distinct Griffiths phases the singular behavior of the disordered ladders is generally controlled by an infinite randomness fixed point.
We consider S=1/2 antiferromagnetic Heisenberg chains with alternating bonds and quenched disorder, which represents a theoretical model of the compound CuCl_{2x}Br_{2(1-x)}(γ-{pic})_2. Using a numerical implementation of the strong disorder renormalization group method we study the low-energy properties of the system as a function of the concentration, x, and the type of correlations in the disorder. For perfect correlation of disorder the system is in the random dimer (Griffiths) phase having a concentration dependent dynamical exponent. For weak or vanishing disorder correlations the system is in the random singlet phase, in which the dynamical exponent is formally infinity. We discuss consequences of our results for the experimentally measured low-temperature susceptibility of CuCl_{2x}Br_{2(1-x)}(γ-{pic})_2.
Quantum phase transition of site-diluted and bond-diluted Heisenberg antiferromagnets on square lattices is studied. By using the novel continuous-time loop algorithm, we perform quantum Monte Carlo simulations on quite larger lattices at extremely lower temperatures than the previous numerical studies. It is found that the antiferromagnetic long-range order at T = 0 persists so long as a cluster of magnetic sites percolates, that is, the critical concentration is equal to the classical percolation threshold, in both of the site-diluted and bond-diluted cases. Furthermore, we find that some critical exponents, such as the magnetization exponent β, are non-classical and strongly depend on the spin size S. On the other hand, we show that the correlation-length exponent ν is universal and is equal to the classical value (ν=4/3).
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.