Dirac function series expansion Dirac function is generalized function. However, is it possible to use taylor expansion on it? i.e., is it possible to write
$$\delta(x)=\sum_{n=0}^{\infty}A_{n}x^{n}$$?
I know that it can have Fourier transform and can be represented as Legendre polynomial, but I would like to know if it can written in terms of Taylor and corresponding $A_{n}$.
 A: Whether you use the limit approach suggested by @pointguardo or anything else, you need to provide the identity:
$$
\int_{-\infty}^\infty \delta(x) dx = 1.
$$
If $\delta$ admits a Taylor expansion, then
$$
\sum_{n = 0}^\infty A_n \int_{-\infty}^\infty x^n dx = 1.
$$
Because of the symmetry of $(-\infty, \infty)$ w.r.t. $x  = 0$, integrals with odd powers vanish and we have
$$
2\sum_{n = 0}^\infty A_{2n}  \int_0^\infty x^{2n} dx = 1.
$$
The integral is nothing else but the Mellin transform of $1$, which is equal to
$$
\int_0^\infty x^{2n} dx = \delta(2n).
$$
Thus, we arrive at
$$
2 \sum_{n = 0}^\infty A_{2n} \delta(2n) = 1.
$$
When $n \neq 0$, $\delta(2n) = 0$. Therefore, the only non-zero term is
$$
2 A_0 \delta(0) = 1, 
$$
which does not hold for any $A_0$.
Thus, the answe is: no.
A: The Taylor series for Dirac delta function:
$$
\delta(x) = \sum_{n=0}^{\infty} A_n x^n. 
$$
At any point $x \neq 0, ~ A_n \equiv 0$ for all $n \ge 0$. So, the series could be useful for the case when $x = 0$. Let's write down using straightforward definition of Taylor series.
$$
\delta(x) = \sum_{n=0}^{\infty} \frac{\delta^{(n)}(0)}{n!}x^n. 
$$
For Dirac function the following relation is true (try to prove it!): $x\delta^{(n)}(x) = -n \delta^{(n-1)}(x)$, so using this recurrent relation one finds:
$$
\delta(x) = \sum_{n=0}^{\infty} \frac{\delta^{(n)}(0)}{n!}x^n = \delta(0) + \sum_{n=1}^{\infty}  \underbrace{(-1)^n\left[\lim_{x\to 0} \frac{\delta(x)}{x^n}\right]}_{A_n} x^n 
$$
However, I don't see how this is going to help you, since you have $\delta(0)$ as a summand which is itself an infinity. 
Another approach that comes to mind (and potentially could be useful) is to try to use the limit definition of Dirac function, namely 
$$
\delta(x) = \lim_{a \to 0}\delta_a(x), 
$$
where $\delta_a(x) = \frac{1}{a\sqrt{\pi}}{e^{-x^2/a^2}}$. From here you can expand this smooth function in Taylor series and get the Taylor for $\delta(x)$ as a limit.
A: Probably it is not possible to find a Dirac delta function representing Taylor series in the real field however, as Jacques Hadamard once said, things change when you complexify i.e. transport to the complex domain the problem you are dealing with. Considering $\delta(x)$, $x\in\mathbb{R}$ as an analytic functional or as a residue current, it is possible to find its representative as a very simple Laurent series and an associated real Taylor series with complex coefficients, which overcomes the issues raised by the very basic observation by Asatur Khurshudyan.
Following Gel'fand and Shilov (1964) (ch. II, §1.5 pp. 161-162), calling $\mathscr{A}(\mathbb{R})$ is the spaces of all entire functions $\psi(s)$, $s=x+iy\in\mathbb{C}$ such that
$$
|s|^q|\psi(s)|\leq C_q e^{a|y|}\quad q\in\mathbb{N}
$$
we can define $\delta(x)$ by the following formula
$$
\psi(0)=\langle \delta,\psi\rangle=\frac{1}{2\pi}\left[\int\limits_{-\infty+ai}^{+\infty+ai}\frac{\psi(s)}{s}\mathrm{d}s-\int\limits_{-\infty-ai}^{+\infty-ai}\frac{\psi(s)}{s}\mathrm{d}s\right]=\frac{1}{2\pi}\int\limits_{|s|=1}\frac{\psi(s)}{s}\mathrm{d}s
\quad a>0\tag{1}\label{1}
$$
In this approach we can say that
$$
\delta(x)=\sum_{n=-\infty}^{+\infty}A_ns^n=s^{-1}=\frac{1}{x+iy}=\sum_{n=0}^{+\infty}\frac{(-1)^n}{(iy)^{n+1}}x^n\tag{2}\label{2}
$$
and already we see that the issue raised by Asatur Khurshudyan is avoided. Approaching the problem in $\mathscr{A}^\prime(\mathbb{R})$ could seem a limited way to deal with it: however, the same approach can be extended to  $\mathscr{D}^\prime(\mathbb{R})$ by understanding \eqref{1} in the sense of residue currents, as done by Berenstein et al. (1993) (ch. 1,§1, pp. 1-2). Following their approach, it is possible to define the following functional:
$$
I(\varphi, \epsilon)=\frac{1}{2\pi}\int\limits_{|s|=\epsilon}\frac{\varphi(s)}{s}\mathrm{d}s \qquad\forall\varphi\in\mathscr{D}(\mathbb{R})
$$
and by their proposition 1.1 (Berenstein et al. (1993), ch. 1, §1, p.2) we have that
$$
\varphi(0)=\langle \delta,\varphi\rangle=\lim_{\epsilon\to 0}I(\varphi, \epsilon)\triangleq\mathrm{res}_{x=0}\frac{\varphi(s)}{s}\qquad\forall\varphi\in\mathscr{D}(\mathbb{R})
$$
[1] Berenstein, C. A.; Gay, R.; Vidras, A.; Yger, A. (1993), Residue currents and Bezout identities, Progress in Mathematics, 114, Basel–Berlin–Boston: Birkhäuser Verlag, pp. xi+158, doi:10.1007/978-3-0348-8560-7, ISBN 3-7643-2945-9 / 0-8176-2945-9 / 978-3-0348-8560-7, MR 1249478, Zbl 0802.32001.
[2] Gel'fand, I. M.; Shilov, G. E. (1964) [1958], Generalized functions. Vol. I: Properties and operations, New York–London: Academic Press, pp. xviii+423, ISBN 978-0-12-279501-5, MR 0166596, Zbl 0115.33101. 
