Multipartite entanglement in the random Ising chain

Quantifying entanglement of multiple subsystems is a challenging open problem in interacting quantum systems. Here, we focus on two subsystems of length l separated by a distance r = alpha l and quantify their entanglement negativity (epsilon) and mutual information (I) in critical random Ising chains. We find universal constant epsilon(alpha) and I(alpha) over any distances, using the asymptotically exact strong disorder renormalization group method. Our results are qualitatively different from both those in the clean Ising model and random spin chains of a singlet ground state, like the spin-1/2 random Heisenberg chain and the random XX chain. While for random singlet states I(alpha)/epsilon(alpha) = 2, in the random Ising chain this universal ratio is strongly alpha dependent. This deviation between systems contrasts with the behavior of the entanglement entropy of a single subsystem, for which the various random critical chains and clean models give the same qualitative behavior. The reason is that epsilon and I are sensitive to higher order correlations in the ground-state structure. Therefore, studying multipartite entanglement provides additional universal information in random quantum systems, beyond what we can learn from a single subsystem.