Using quantum Monte Carlo simulations, we study effects of disorder on the
S=1/2 Heisenberg model with exchange constant
J on the square lattice supplemented by multispin interactions
Q. It was found recently [L. Liu et al., Phys. Rev. X 8, 041040 (2018)] that the ground state of this $J\ensuremath{-}Q$ model with random couplings undergoes a quantum phase transition from the N\'eel antiferromagnetic state into a randomness-induced spin-liquid-like state that is a close analog to the well known random-singlet (RS) state of the Heisenberg chain with random couplings. This 2D RS state arises from a spontaneously symmetry-broken fourfold degenerate columnar valence-bond solid that is broken up by the disorder into finite domains, with spinons localized at topological defects. The interacting spinons form a critical collective many-body state without magnetic long range order but with the mean spin-spin correlations decaying with distance
r as ${r}^{\ensuremath{-}2}$, as in the one-dimensional RS state. The dynamic exponent $z\ensuremath{\ge}2$, varying continuously with the model parameters. In this work, we further investigate the properties of the RS state in the $J\ensuremath{-}Q$ model with random
Q couplings. We study the temperature dependence of the specific heat and various susceptibilities for large enough systems to reach the thermodynamic limit. We also analyze the size dependence of the critical magnetic order parameter and its susceptibility in the ground state. For all these quantities, we find consistency with the conventional quantum-critical scaling laws when the condition implied by the ${r}^{\ensuremath{-}2}$ form of the spin correlations is imposed. In particular, all the different quantities can be explained by the same value of the dynamic exponent
z at fixed model parameters. We argue that the RS state identified in the $J\ensuremath{-}Q$ model corresponds to a generic renormalization group fixed point that can be reached in many quantum magnets with random couplings and that it has already been observed experimentally.
We introduce and implement a reformulation of the strong-disorder renormalization group method in real space, well suited to study bond-disordered antiferromagnetic power-law coupled quantum spin chains. We apply it to a chain of randomly placed spins coupled by power-law long-range interaction with power $\ensuremath{\alpha}$. First, keeping only interactions between adjacent spins, we derive the Master equation and confirm that it is solved by the infinite randomness fixed-point distribution. Then, we derive the master equation for power-law long-range interaction between all spins for any anisotropy $\ensuremath{\gamma}$ ranging from the XX limit, $\ensuremath{\gamma}=0$, to the isotropic Heisenberg limit, $\ensuremath{\gamma}=1,$ which corresponds to a tight-binding chain of disordered long-range interacting fermions with long-range hopping. We thereby show that the distribution function of couplings
J smaller than renormalization scale $\mathrm{\ensuremath{\Omega}}$ flows to the strong-disorder fixed-point distribution of finite width $2\ensuremath{\alpha}.$ We find only small corrections to that distribution, which depend on power exponent $\ensuremath{\alpha}$ and coupling anisotropy $\ensuremath{\gamma}.$ As a consequence, the low-temperature magnetic susceptibility diverges with an anomalous power law. The distribution of singlet lengths
l is found to decay as ${l}^{\ensuremath{-}2}$. The entanglement entropy of a subsystem of length
n increases in the ground state logarithmically for all $\ensuremath{\alpha}$ and $\ensuremath{\gamma}$. After a global quantum quench, the entanglement entropy increases with time logarithmically as $S(t)\ensuremath{\sim}ln(t)/(2\ensuremath{\alpha})$.
We study the infinite-temperature properties of an infinite sequence of random quantum spin chains using a real-space renormalization group approach, and demonstrate that they exhibit nonergodic behavior at strong disorder. The analysis is conveniently implemented in terms of SU(2)_{k} anyon chains that include the Ising and Potts chains as notable examples. Highly excited eigenstates of these systems exhibit properties usually associated with quantum critical ground states, leading us to dub them "quantum critical glasses." We argue that random-bond Heisenberg chains self-thermalize and that the excited-state entanglement crosses over from volume-law to logarithmic scaling at a length scale that diverges in the Heisenberg limit k→∞. The excited state fixed points are generically distinct from their ground state counterparts, and represent novel nonequilibrium critical phases of matter.
Low-dimensional spin systems, consisting of arrays of spins arranged in chains or ladders, have been investigated intensively in recent years using both exactly solvable theoretical models as well as a variety of experimental techniques. In case of random variations of the exchange couplings the renormalization group theory predicts the existence of a random-singlet (RS) state, corresponding to spins coupled at all possible distances and energy scales. The scarce availability of suitable random systems, however, has so far prevented the experimental identification of this peculiar magnetic ground state. In a combined effort using nuclear magnetic resonance (NMR), dc magnetometry and numerical simulations, we now find compelling evidence for the formation of a random-singlet state in this class of materials. Randomness seems to generate a broad distribution of local magnetic fields, in turn reflected in stretched exponential NMR relaxations. This distribution exhibits a progressive broadening with decreasing temperature, caused by the growing inequivalence of magnetic sites, as expected from the RS theory. Our findings suggest that NMR is the tool of choice for probing the low-energy physics also in other disordered magnets, where extended-scale excitations are dominant.
We study dynamical properties at finite temperature
(T) of Heisenberg spin chains with random antiferromagnetic exchange couplings, which realize the random singlet phase in the low-energy limit, using three complementary numerical methods: exact diagonalization, matrix-product-state algorithms, and stochastic analytic continuation of quantum Monte Carlo results in imaginary time. Specifically, we investigate the dynamic spin structure factor $S(q,\ensuremath{\omega})$ and its $\ensuremath{\omega}\ensuremath{\rightarrow}0$ limit, which are closely related to inelastic neutron scattering and nuclear magnetic resonance (NMR) experiments (through the spin-lattice relaxation rate
1/T1). Our study reveals a continuous narrow band of low-energy excitations in $S(q,\ensuremath{\omega})$, extending throughout the
q space, instead of being restricted to $q\ensuremath{\approx}0$ and $q\ensuremath{\approx}\ensuremath{\pi}$ as found in the uniform system. Close to $q=\ensuremath{\pi}$, the scaling properties of these excitations are well captured by the random-singlet theory, but disagreements also exist with some aspects of the predicted
q dependence further away from $q=\ensuremath{\pi}$. Furthermore we also find spin diffusion effects close to
q=0 that are not contained within the random-singlet theory but give non-negligible contributions to the mean
1/T1. To compare with NMR experiments, we consider the distribution of the local relaxation rates
1/T1. We show that the local
1/T1 values are broadly distributed, approximately according to a stretched exponential. The mean
1/T1 first decreases with
T, but below a crossover temperature it starts to increase and likely diverges in the limit of a small nuclear resonance frequency ${\ensuremath{\omega}}_{0}$. Although a similar divergent behavior has been predicted and experimentally observed for the static uniform susceptibility, this divergent behavior of the mean
1/T1 has never been experimentally observed. Indeed, we show that the divergence of the mean
1/T1 is due to rare events in the disordered chains and is concealed in experiments, where the typical
1/T1 value is accessed.
Using an adaptive strategy which enables the study of quenched disordered system via the density-matrix renormalization-group method, we compute the various ground-state spin-spin correlation measures of the spin-
21 antiferromagnetic Heisenberg chain with random coupling constants, namely, the mean values of the bulk and of the end-to-end correlations, the typical value of the bulk correlations, and the distribution of the bulk correlations. Our results are in agreement with the predictions of the strong-disorder renormalization-group method. We do not find any hint of logarithmic corrections either in the bulk average correlations, which were recently reported by Shu et al. [Phys. Rev. B 94, 174442 (2016)], or in the end-to-end average correlations. We report the existence of a logarithmic correction on the end-to-end correlations of the clean chain. Finally, we have determined that the distribution of the bulk correlations, when properly rescaled by an associated Lyapunov exponent, is a narrow and universal (disorder-independent) probability function.