materials. In the heavy fermion CeCo(In$_{1-x}$Cd$_{x}$)$_{5}$, for example, an antiferromagnetic region appears around Cd impurities, with a size that can be tuned with pressure [@Urbano07; @Seo14]. In this particular situation, a more appropriate model than that of Eq. is one in which the moment on an impurity site has an altered hybridization $V$ to the conduction electrons[@Benali16]. However, magnetic domains around sites in which the moment is removed have also been studied[@martin97; @azuma97].

![(Color online) Contour plot of the conduction electron density of states, $N^{c}_{\mathbf{i}}(0)$, for $\beta t = 20$. []{data-label="Fig:local_DOS_Nc_contour"}](DOS_L160_Nc.pdf)

The preceding result is suggestive of the breaking of the local singlet state, an effect we will see in even more dramatic form when a collection of $f$ sites is removed. It is also worth examining the spectral properties of the system. We compute the local density of states (DOS) by analytic continuation of the imaginary-time dependent Green’s function, inverting the integral equation $$\begin{aligned}
G_{\mathbf{i}}(\mathbf{j}=0,\tau) = \int \mathrm{d}\omega\,
N_{\mathbf{i}}(\omega) \,
\frac{e^{-\omega \tau}}{e^{\beta \omega} + 1},
\label{eq:aw}\end{aligned}$$ where $\mathbf{i}$ denotes the site position, while $\mathbf{j}$ is the displacement between sites where the creation and annihilation operators of the Green’s function are applied. As discussed in Ref., for low temperatures the 