of the lattice. Because one of the most interesting aspects of site removal is the possibility of enhancing magnetism[@Mendes17], we fix the hybridization at $V/t=1.2$ and the repulsive potential as $U_{f}/t=4$, so that we are in the spin-singlet phase of the undepleted PAM[@Vekic95; @Hu17].

![(Color online) Contour plot of the $f$-orbital local susceptibilities for $\beta t = 20$. The central black square marks the geometrical position of the site where the $f$-orbital was depleted.[]{data-label="Fig:local_chi_contour"}](local_Chi_beta20.pdf)

Figure \[Fig:local\_chi\] presents the behavior of the local susceptibility, Eq., of the $f$-orbitals as a function of the distance from the depletion site. At high temperature, $T/t=0.2$ ($\beta
t=5$, black diamonds), the magnetic response is large and positive, and almost homogeneous throughout the lattice. When the temperature is decreased, $T/t=0.10 - 0.05$, the local susceptibility increases on nearest neighbor (NN) sites. This is the opposite of what happens in the conventional PAM, where, in the singlet phase at $V/t=1.2$, the magnetic susceptibility goes to zero as a consequence of the spin gap in the ground state. Indeed, at $\beta t=20$ (red squares), the NN magnetic susceptibility is an order of magnitude larger than the undepleted PAM (dashed black line).

On the other hand, as $T$ is 