The local response to a uniform field around vacancies in the two-dimensional (2D) spin-1/2 Heisenberg antiferromagnet is determined by numerical quantum Monte Carlo simulations as a function of temperature. It is possible to separate the Knight shifts into uniform and staggered contributions on the lattice which are analyzed and understood in detail. The contributions show interesting long and short range behavior that may be of relevance in NMR and susceptibility measurements. For more than one impurity remarkable non-linear enhancement and cancellation effects take place. We predict that the Curie impurity susceptibility will be observable for a random impurity concentration even in the thermodynamic limit.PACS numbers: 75.10. Jm, 74.25.Nf, 75.20.Hr, 75.40.Mg The deliberate use of substitutional impurities in strongly correlated electron systems has become a valuable tool for controlled studies of the underlying correlations [1,2,3,4,5]. Doping of non-magnetic Znions in antiferromagnetic CuO 2 planes is known to lead to a surprisingly large reduction of T c [1] and induces local staggered magnetic moments around the impurity sites [2]. It has now become quite standard to analyze the local magnetic moments around static magnetic and non-magnetic impurities in an antiferromagnetic background for a deeper understanding of the correlated states [3,4,5].Theoretically the Knight shifts around impurities and vacancies have been studied for many low-dimensional antiferromagnets [6,7,8,9,10,11,12], which typically show a strong enhancement of the antiferromagnetic order. The detailed behavior of the staggered magnetic moments can even give rather exotic results such as an increasing amplitude as a function of distance [7] and show renormalization effects [11,12]. The magnetization pattern can often be interpreted as the interference of incoming and scattered quasiparticle excitations [11] or in terms of a pruned valence bond basis [8].The sum of all Knight shifts adds up to the total susceptibility χ 1 . The resulting impurity susceptibility χ imp = χ 1 − χ 0 has been studied in detail in more recent theoretical studies [13,14,15,16,17,18,19,20]. In systems with long range magnetic order the leading temperature dependence of the impurity susceptibility from one single vacancy is given by a classical Curie spin which has been confirmed by numerical simulations for the 2D Heisenberg antiferromagnet [14,15], where a subleading logarithmic term has also been established [13,14,15,16,17] Here χ 0 is the total susceptibility without any impurities, ρ s is the spin stiffness, and the limit of large correlation length ξ(T ) and system size ξ(T ) > L → ∞ is assumed.We now want to examine how this impurity susceptibility is distributed on the lattice by considering the linear local response to a small uniform magnetic field (Knight shift) at each siteWe find that long and short range patterns of the Knight shifts contribute to the impurity susceptibility differently. We are able to study the interference of several impuritie...