Let $A$ be an unbounded self-adjoint operator on $L^2(\mathbb{R})$ and we are assuming the $L^2$ functions to be complex-valued.
We further assume (e.g. compactness of resolvent) that there exists an orthonormal eigenbasis $\{ \phi_n \}$ of $L^2(\mathbb{R})$ with real eigenvalues $\lambda_n$: \begin{equation} \int_{\mathbb{R}} \overline{\phi_n}(x) \phi_m(x) dx=\delta_{m,n} \text{ , } A\phi_n=\lambda_n \phi_n \end{equation}
Moreover, the completeness relation is found in physics liteature as follows: \begin{equation} \sum_{n}\phi_n(x) \overline{\phi_n}(y)=\delta(x-y) \end{equation} where $\delta(x-y)$ is the Dirac-delta distribution. My questions are as follows:
What is the most rigorous way I can interpret the completeness relation above? Perhaps in the sense of tempered distributions?
Since $A$ is self-adjoint, functional calculus allows us to define $e^{-A}$ as a bounded operator on $L^2(\mathbb{R})$. Then, is it possible to carry out a computation like the below? \begin{equation} e^{-{A_x}} \delta(x-y)=\text{some smooth function in }x,y = \sum_{n} e^{-{A_x}} \phi_n(x) \overline{\phi_n}(y)=\sum_n e^{-\lambda_n} \phi_n(x)\overline{\phi_n}(y) \end{equation}
where I used the notation $A_x$ to mean that the operator $A$ acts on the $x$ argument only.
- If the computation in item 2 is possible, in what sense should I understand it?
These seem quite subtle problems between math and physics.. Could anyone please clarify for me?