I know very well that Laplacian in bounded domain has a discrete spectrum. How about Laplacian in $\mathbb{R}^n$?(not in some fancy-shaped unbounded domain, but the whole domain)
Where can I find such results?
Moreover, is there a counterpart of Hilbert-Schmidt theorem for Laplacian in $\mathbb{R}^n$? Hilbert-Schmidt asserts there is a countable set of eigenfunctions $\phi_n$ so that $x=\sum \langle x,\phi_n\rangle \phi_n,\forall x\in H$.
Is there a similar theorem saying $x=\int_0^\infty \langle x,\phi_\lambda\rangle \phi_\lambda\,\mathrm{d}\lambda,\forall x\in H$ where $\phi_\lambda$ is the eigenfunction of Laplacian to spectral value $\lambda$?
Preliminaries: I will use the following sign convention for the Laplacian: $ \Delta u := - \sum_{j=1}^n \frac{\partial^2 u }{\partial x_j^2}.$ $\Delta$ is an unbounded operator on $L^2(\mathbb{R}^n)$. To define the domain of $\Delta$, recall that for $u \in C^\infty_0$, $$\mathcal{F}(\Delta u)(\xi) = 4\pi^2 |\xi|^2 \hat{u}(\xi).$$ (I use $\mathcal{F}(\phi)$ and $\hat{\phi}$ interchangeably to denote the Fourier transform of a function $\phi$.) Let's take the domain of $\Delta$ to be $$D(\Delta) := \{ u \in L^2 ~:~ 4\pi^2 |\xi|^2 \hat{u}(\xi) \in L^2 \}.$$ This makes $\Delta$ a closed unbounded operator.
Recall that the resolvent set of $\Delta$, denoted by $\rho(\Delta)$, is defined as the set of complex numbers $\lambda$ such that $\lambda I - \Delta$ is a bijection $D(\Delta) \to L^2$ (note: in this case, the closed graph theorem gives that the resolvent $(\lambda I - \Delta)^{-1}$ is necessarily bounded). The spectrum of $\Delta$, denoted by $\sigma(\Delta)$, is defined as the complement of $\rho(\Delta)$.
My goal is to sketch a proof of the following claim: \begin{equation} \sigma(\Delta) = [0, \infty) \end{equation}
First note that for $\lambda \in \mathbb{C} \backslash [0, \infty)$, $\lambda$ is in the resolvent set since we can easily invert $\lambda I - \Delta$ using the Fourier transform: $$(\lambda I - \Delta) u = f \iff (\lambda - 4\pi^2|\xi|^2)\hat{u} = \hat{f} \iff \hat{u} = (\lambda - 4\pi^2|\xi|^2)^{-1} \hat{f} .$$ Thus for $f \in L^2$, $(\lambda I - \Delta)^{-1} f = \mathcal{F}^{-1}\left((\lambda - 4\pi^2|\xi|^2)^{-1} \hat{f}\right)$. (Of course, this only works when $\lambda$ is not a nonnegative real number.) This proves that $\sigma(\Delta) \subset [0, \infty)$.
Now we'll show that $[0, \infty) \subset \sigma(\Delta)$. Let $\lambda \in [0, \infty)$. To prove that $\lambda \in \sigma(\Delta)$, it suffices to exhibit a sequence of functions $u_k$ in $D(\Delta)$ such that $$\frac{||u_k||_{L^2}}{||(\lambda I - \Delta)u_k||_{L^2}} \to \infty \text{ as } k \to \infty $$ as this shows that the resolvent $(\lambda I - \Delta)^{-1}$ cannot possibly be bounded. To that end, pick a point $x_0 \in \mathbb{R}^n$ such that $|x_0|^2 = \lambda$, and define u(x) by $$u(x) = e^{ix_0 \cdot x}.$$ (Note that $u$ is an eigenfunction for the Laplacian with eigenvalue $\lambda = |x_0|^2$, but $u$ is not in $L^2$.) Choose a sequence of cutoff functions $\phi_k \in C_0^\infty$ such that $0 \leq \phi_k \leq 1$, $\text{supp}(\phi_k) \subset B_{k+1}(0)$, and $\phi_k(x) \equiv 1$ for $x \in B_k(0)$. We can choose the $\phi_k$'s so that all their first two partial derivatives are uniformly bounded in $k$. Define the sequence $u_k \in C_0^\infty \subset L^2$ by $$u_k := \phi_k u.$$ Then $(\lambda I - \Delta)u_k$ is supported in the annulus $B_{k+1}(0) \backslash B_k(0)$, which has volume that is $O(k^{n-1})$ as $k \to \infty$. The set $\{ (\lambda I - \Delta)u_k \}$ is uniformly bounded in the sup norm, so there exists $C$ such that $$||(\lambda I - \Delta)u_k ||_{L^2}^2 \leq C k^{n-1}.$$ On the other hand, $|u_k (x)| \equiv 1$ for $x \in B_k(0)$, so we have the upper bound $$||u_k||_{L^2}^2 \geq \text{vol}(B_k(0)) = \omega_n k^n.$$ The last two inequalities together give the desired result.