I read on this Sobolev space with negative index question that the dirac delta distribution belongs to $H^{-1/2-\epsilon}(\mathbb{R})$. I was trying to figure out why that is the case but it wasn't clear to me. I know that it acts on test functions by $\langle \delta, f\rangle = f(0)$ Although I don't see why the test functions necessarily need to be in $H^{1/2}$.
Finally, I also read that $\delta^2$ is not well defined as a distribution. How would I see this? Is it simply because if we apply $\delta$ to a function once, we get a value in $\mathbb{R}$. But then why not treat $f(0)$ like a constant function that you can then apply $\delta$ to again?
I can see the issue with naively defining $\langle \delta^2,f\rangle = f(0)^2$, because here linearity breaks down. But is there a general proof for why there's no good way to define $\delta^2$?