everyone. I am wondering whether any one knows if the square of Dirac Delta function is defined somewhere. 

In the beginning, this question might look strange. But by restricting the space of the test functions, I think it is still possible. For example, in order to make sense of $\delta_0^2$, we can think that it is the limit of $\frac{e^{-x^2/t}}{2\pi t}$ as $t\rightarrow 0_+$. Now choose the test function $f(x)=x^2$. It is clear that
$$
\int_{-\infty}^{\infty} x^2 \frac{e^{-x^2/t}}{2\pi t} d x = \frac{1}{2\sqrt{\pi t}} \int_{-\infty}^{\infty} x^2 \frac{e^{-x^2/t}}{\sqrt{\pi t}} d x = \frac{1}{2\sqrt{\pi t}} \cdot \frac{t}{2} = \frac{\sqrt{t}}{4\sqrt{\pi}}\;.
$$
Then let $t$ tend to $0$, we get $<\delta_0^2,f>=0$ in this case. So we can restrict, for example, all test functions tend to 0 at the speed no less than $x^2$.

I don't want to invent the whole stuff if it already exists. Otherwise, I might take care of the every details. Thank you in advance for any hints.


EDIT:
Here are some references that I found to be useful.
1: Mikusiński, J. [On the square of the Dirac delta-distribution][1]. (Russian summary) 
Bull. Acad. Polon. Sci. Sér. Sci. Math. Astronom. Phys. 14 1966 511–513. 
44.40 (46.40) 

2:Ta Ngoc Tri, [The Colombeau theory of generalized functions][2] Master thesis, 2005

  [1]: http://www.ams.org/mathscinet/search/publdoc.html?arg3=&co4=AND&co5=AND&co6=AND&co7=AND&dr=all&pg4=AUCN&pg5=TI&pg6=PC&pg7=ALLF&pg8=ET&r=1&review_format=html&s4=&s5=on%2520the%2520square%2520of%2520the%2520dirac%2520delta&s6=&s7=&s8=All&vfpref=html&yearRangeFirst=&yearRangeSecond=&yrop=eq
  [2]: http://www.science.uva.nl/onderwijs/thesis/centraal/files/f59708808604356.pdf