We study the field dependence of the entanglement of formation in anisotropic S=1/2 antiferromagnetic chains displaying a T=0 field-driven quantum phase transition. The analysis is carried out via quantum Monte Carlo simulations. At zero temperature the entanglement estimators show abrupt changes at and around criticality, vanishing below the critical field, in correspondence with an exactly factorized state, and then immediately recovering a finite value upon passing through the quantum phase transition. At the quantum-critical point, a deep minimum in the pairwise-to-global entanglement ratio shows that multispin entanglement is strongly enhanced; moreover this signature represents a novel way of detecting the quantum phase transition of the system, relying entirely on entanglement estimators.

We study the pairwise entanglement close to separable ground states of a class of one dimensional quantum spin models. At T=0 we find that such ground states separate regions, in the space of the Hamiltonian parameters, which are characterized by qualitatively different types of entanglement, namely parallel and antiparallel entanglement; we further demonstrate that the range of the Concurrence diverges while approaching separable ground states, therefore evidencing that such states, with uncorrelated fluctuations, are reached by a long range reshuffling of the entanglement. We generalize our results to the analysis of quantum phase transitions occurring in bosonic and fermionic systems. Finally, the effects of finite temperature are considered: At T >0 we evidence the existence of a region where no pairwise entanglement survives, so that entanglement, if present, is genuinely multipartite.

Making use of exact results and quantum Monte Carlo data for the entanglement of formation, we show that the ground state of anisotropic two-dimensional S=1/2 antiferromagnets in a uniform field takes the classical-like form of a product state for a particular value and orientation of the field, at which the purely quantum correlations due to entanglement disappear. Analytical expressions for the energy and the form of such states are given, and a novel type of exactly solvable two-dimensional quantum models is therefore singled out. Moreover, we show that the field-induced quantum phase transition present in the models is unambiguously characterized by a cusp minimum in the pairwise-to-global entanglement ratio R, marking the quantum-critical enhancement of multipartite entanglement.