Entanglement microscopy reveals the true quantum correlations among the microscopic building blocks of many-body systems [Nat. Commun. 16, 96 (2025)]. Using this approach, we study the multipartite entanglement of the quantum Ising model in 1D, 2D, and 3D. We first obtain the full reduced density matrix (tomography) of subregions that have at most four sites via quantum Monte Carlo, exact diagonalization, and the exact solution in 1D. We then analyze both bipartite and genuine multipartite entanglement (GME) among the sites in the subregion. To do so, we use a variety of measures including the negativity, as well as a true measure of GME: the genuinely multipartite concurrence (or GME concurrence), and its computationally cheaper lower bound
I2. We provide a complete proof that
I2 bounds the GME concurrence, and show how the symmetries of the state simplify its evaluation. For adjacent sites, we find three- and four-spin GME present across large portions of the phase diagram, reaching a maximum near the quantum critical point. In 1D, we identify the singular scaling of the derivative
dI2/dh approaching the critical point. We observe a sharp decrease in GME with increasing dimensionality, consistent with the monogamous nature of entanglement. Furthermore, we find that GME vanishes for subregions consisting of nonadjacent sites in both 2D and 3D, offering a stark illustration of the short-ranged nature of entanglement in equilibrium quantum matter [arXiv:2402.06677]. Finally, we analyze the most collective form of entanglement by evaluating the GME concurrence among all spins in the lattice, which can be obtained from a simple observable: the single-site transverse magnetization. This global concurrence is larger in 1D compared to 2D and 3D, but it is relatively less robust against perturbations such as local measurements.
We investigate the emergence and properties of long-distance genuine multipartite entanglement, induced via three localized magnetic defects, in a one-dimensional transverse-field
XX spin-
1/2 chain. Using both analytical and numerical techniques, we determine the conditions for the existence of bound states localized at the defects. We find that the reduced density matrix (RDM) of the defects exhibits long-distance genuine multipartite entanglement (GME) across the whole range of the Hamiltonian parameter space, including regions where the two-qubit concurrence is zero. We quantify the entanglement by using numerical lower bounds for the GME concurrence, as well as by analytically deriving the GME concurrence in regions where the RDM is of rank 2. Our work provides insights into generating multipartite entanglement in many-body quantum systems via local control techniques.