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.
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 train machine learning algorithms to infer the entanglement structure of disordered long-range interacting quantum spin chains by learning from the strong disorder renormalisation group (SDRG) method. The system consists of
S=1/2-quantum spins coupled by antiferromagnetic power-law interactions with decay exponent
α at random positions on a one-dimensional chain. Using SDRG as a physics-informed teacher, we compare a Random Forest classifier as a classical baseline with a graph neural network (GNN) that operates directly on the interaction graph and learns a bond-ranking rule mirroring the SDRG decimation policy. The GNN achieves a disorder-averaged pairing accuracy close to one and reproduces the entanglement entropy
S(ℓ) in excellent quantitative agreement with SDRG across all subsystem sizes and interaction exponents. RG flow heat maps confirm that the GNN learns the sequential decimation hierarchy rather than merely fitting final-state observables. Finite-temperature entanglement properties are incorporated via the SDRGX framework through a two-stage strategy, using the zero-temperature GNN to generate the RG flow and sampling thermal occupations from the canonical ensemble, yielding results in agreement with both numerical SDRGX and analytical predictions without retraining.
Using strong-disorder renormalization group, numerical exact diagonalization, and quantum Monte Carlo methods, we revisit the random antiferromagnetic
XXZ spin-
1∕2 chain focusing on the long-length and ground-state behavior of the average time-independent spin-spin correlation function $C(l)=\ensuremath{\upsilon}{l}^{\ensuremath{-}\ensuremath{\eta}}$. In addition to the well-known universal (disorder-independent) power-law exponent $\ensuremath{\eta}=2$, we find interesting universal features displayed by the prefactor $\ensuremath{\upsilon}={\ensuremath{\upsilon}}_{o}∕3$, if
l is odd, and $\ensuremath{\upsilon}={\ensuremath{\upsilon}}_{e}∕3$, otherwise. Although ${\ensuremath{\upsilon}}_{o}$ and ${\ensuremath{\upsilon}}_{e}$ are nonuniversal (disorder dependent) and distinct in magnitude, the combination ${\ensuremath{\upsilon}}_{o}+{\ensuremath{\upsilon}}_{e}=\ensuremath{-}1∕4$ is universal if
C is computed along the symmetric (longitudinal) axis. The origin of the nonuniversalities of the prefactors is discussed in the renormalization-group framework where a solvable toy model is considered. Moreover, we relate the average correlation function with the average entanglement entropy, whose amplitude has been recently shown to be universal. The nonuniversalities of the prefactors are shown to contribute only to surface terms of the entropy. Finally, we discuss the experimental relevance of our results by computing the structure factor whose scaling properties, interestingly, depend on the correlation prefactors.
Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\ensuremath{\alpha}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $\ensuremath{\alpha}>\frac{1}{2}$. A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\ensuremath{\alpha}$, corresponding to a critical behavior with an effective central charge
c=ln(2), independent of $\ensuremath{\alpha}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\ensuremath{\alpha}$. For a wide range of partition size
l, the EE fits a critical behavior with a central charge close to
c=1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\ensuremath{\alpha}=2$. Only for small $l\ensuremath{\ll}L$, in a range which increases with the number of spins
N, we find deviations which are rather consistent with the strong disorder fixed point central charge
c=ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. We also present results obtained with DMRG and find agreement with ED for sufficiently small $\ensuremath{\alpha}<2$, whereas for larger $\ensuremath{\alpha}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.
We study the scaling properties of entanglement entropy (EE) near the quantum critical points in interacting random antiferromagnetic (AFM) spin chains. Using density-matrix renormalization group, we compute the half-chain EE near the topological phase transition between the Haldane and random singlet phases in a disordered spin-1 chain. It is found to diverge logarithmically in system size with an effective central charge
ceff=1.17(4) at the quantum critical point (QCP). Moreover, a scaling analysis of EE yields the correlation length exponent $\ensuremath{\nu}=2.28(5)$. Our unbiased calculation establishes that the QCP is in the universality class of the infinite-randomness fixed point predicted by previous studies based on the strong disorder renormalization group technique. However, in the disordered spin-1/2 Majumdar-Ghosh chain, where a valence bond solid phase is unstable to disorder, the crossover length exponent obtained from a scaling analysis of EE disagrees with the expectation based on the Imry-Ma argument. We provide a possible explanation.