# How to rigorously understand continuous bases?

In Quantum Mechanics it is quite common to see the idea of a continuous basis of a Hilbert space. In truth if $\mathcal{H}$ is the state space of a quantum system and if $X : U\subset \mathcal{H}\to \mathcal{H}$ is the position operator then it is quite usual to consider the set of eigenvectors $\{|x\rangle : x \in \mathbb{R}\}$ indexed by the eigenvalues $x\in \mathbb{R}$.

In that case, Quantum Mechanics books usually assume that one can write for any $|\psi\rangle \in \mathcal{H}$

$$|\psi\rangle = \int_{-\infty}^{\infty}|x\rangle\langle x|\psi\rangle dx.$$

And also assume we have the completeness relation

$$\int_{-\infty}^{\infty}|x\rangle \langle x| \ dx = I,$$

being $I$ the identity operator and considering $|x\rangle\langle x|$ to be the operator defined by

$$(|x\rangle \langle x|)|\psi\rangle = |x\rangle \langle x|\psi\rangle.$$

Although this works and gives the correct results, I can't understand what is going on here. One of the points is that the integral involved here is quite strange. We are integrating $f(x) = |x\rangle \langle x|\psi\rangle$ but $f : \mathbb{R}\to \mathcal{H}$ and I'm not sure how this integral is defined.

Also, I know that for bounded operators with discrete spectrum the set of eigenvectors does really form a basis in the usual sense, but here things are quite different.

How does one rigorously understand those basis? What is the correct definitions behind all of this?

• One can roughly think of a Hilbert space as the real line with a vector structure such that each point $x \in \mathbb{R}$ is a one dimensional vector space, and every $x$ is orthogonal to every $x' \neq x$. In this case, the completeness relation is just an infinite dimension version of $x x^t = I$ for an orthonormal basis $\{x\}$. Further, for the $f(x)$ bit: For every $x \in \mathbb{R}$, there exists a $|x \rangle \in \mathcal{H}$, giving you a bijection between $\mathbb{R}$ and $\mathcal{H}$. Then defining functions $\mathbb{R} \to \mathcal{H}$ shouldn't be a problem. Commented Oct 1, 2015 at 2:07
• For a very thorough account, see Mathematical surprises and Dirac's formalism in quantum mechanics. It's related to TrialAndError's answer. Commented Oct 1, 2015 at 10:50
• At the end of the day, I think this is all formal manipulation with spectral projections. Commented Oct 1, 2015 at 16:54

For example, consider a periodic problem on $[-\pi,\pi]$. The opertor $L=-\frac{d^{2}}{dx^{2}}$ has a two-dimensional eigenspace spanned by $\{ e^{inx},e^{-inx}\}$ for each eigenvalue $\lambda = 1^{2}, 2^{2}, 3^{2},\cdots$ and has a one-dimensional eigenspace spanned by $\{ 1 \}$ for eigenvalue $\lambda=0$. Ignoring the fact that the expansions cannot be written as a integral without some type of delta function, there is still the issue of dimension of the eigenspaces.
Dirac modeled his representation after the ODE representations for one-dimensional Sturm-Liouville operators $$Lf = -\frac{d}{dx}\left(p(x)\frac{df}{dx}\right)+qf.$$ For example, on a semi-infinite interval such as $[0,\infty)$ where $p$ is contnuously differentiable and non-vanishing, and where $q$ is continuous, but possibly unbounded as $x\rightarrow\infty$, you can impose an endpoint condition at $0$ of the form $Af(0)+Bf'(0)=0$ where $A$, $B$ are real and not both $0$, and the eigenfunctions $\varphi_{\lambda}(x)$ satisfying fixed endpoint condition $\varphi_{\lambda}(0)=-B$, $\varphi_{\lambda}'(0)=A$, and you'll get an expansion $$f(x) = \int_{-\infty}^{\infty}\left(\int_{0}^{\infty}f(y)\varphi_{\lambda}(y)dy\right)\varphi_{\lambda}(x)d\rho(\lambda).$$ Allowing for a general measure $\rho$ (or, equivalently, using a Riemann-Stieltjes integral with respect to a function $\rho$ of bounded variation on finite intervals), you do get what Dirac had in mind. This is what Dirac was working from. Further, you get Parseval's equality: $$\int_{0}^{\infty}|f(x)|^{2}dx = \int_{-\infty}^{\infty}\left|\int_{0}^{\infty}f(y)\varphi_{\lambda}(y)dy\right|^{2}d\rho(\lambda)$$ I think you can convince yourself that Dirac's notation is lifted from such a case, and that it can be interpreted as being correct in such a case. This formalism can be extended with a little care to all of $\mathbb{R}$ by patching representations on $(-\infty,0]$ with $[0,\infty)$. And there is no doubt why Dirac's notation gives you something correct for the cases where separation of variables applies and where the operators become classical differential operators.
For general selfadjoint operators, Dirac notation is not particularly good. A general selfadjoint operator $A$ can be represented as a spectral integral proposed by von Neumann in the form $$Af = \int_{-\infty}^{\infty}\lambda dE(\lambda)f.$$ Notice that von Neumann's notation (which is what I have given above) is related to Quantum as well, for which he was offering a rigorous alternative in the notation of Quantum. What is happening in this integral is that $E(\lambda)$ is an orthogonal projection for each $\lambda \in\mathbb{R}$, and the space onto which $E(\lambda)$ projects is increasing with $\lambda$. The underlying space is some complex Hilbert space $\mathcal{H}$, and $$E(\lambda_1)\mathcal{H} \subseteq E(\lambda_2)\mathcal{H},\;\;\;\lambda_1 \le \lambda_2.$$ The fact that $E(\lambda)$ is an orthogonal projection is equivalent to $$E(\lambda)^{2}=E(\lambda)^{\star} = E(\lambda).$$ Because $E$ is increasing, $E(\lambda_2)-E(\lambda_1)$ is an orthogonal projection, and $$A(E(\lambda+h)-E(\lambda))f \approx \lambda(E(\lambda+h)-E(\lambda))f$$ In infinitesimal notation: $$AdE(\lambda)f = \lambda dE(\lambda)f,$$ which is correct in an integral sense $$A\int_{-\infty}^{\infty}dE(\lambda)f=\int_{-\infty}^{\infty}\lambda dE(\lambda)f,\;\;\; f \in \mathcal{D}(A).$$ And you have rigorously correct Parseval identities: $$\left\|\int_{-\infty}^{\infty}dE(\lambda)f\right\|^{2} =\int_{-\infty}^{\infty}d\|E(\lambda)f\|^{2},\\ \|Af\|^{2} = \left\|\int_{-\infty}^{\infty}\lambda dE(\lambda)f\right\|^{2} = \int_{-\infty}^{\infty}\lambda^{2}d\|E(\lambda)f\|^{2}.$$ This version of the Spectral Theorem for selfadjoint linear operators holds very generally on any Hilbert space, and it dates back to the 1930's. I would encourage anyone serious about Quantum Mechanics to become familiar with it. Attractive as it may be, Dirac's formalism is not adequate for a foundational approach.
• @FrankScience : No matter what, you can't really have $$|\psi\rangle = \int_{-\infty}^{\infty}|x\rangle\langle x|\psi\rangle dx.$$ That's a one-dimensional representation without a discrete component and with uniform spectral density. You can overcome the uniform spectral density, and add in discrete sums, but you can't get a one-dimensional representation; you would need multiple such representations, and discrete components. With these mods for Pseudo Differential Operators, you might be able to make something work. Commented Nov 28, 2015 at 19:25