# Hilbert space

In mathematics, a Hilbert space (named for David Hilbert) generalizes the notion of Euclidean space. It extends the methods of linear algebra and calculus from the two-dimensional Euclidean plane and three-dimensional space to spaces with any finite or infinite number of dimensions. A Hilbert space is a vector space equipped with an inner product operation, which allows lengths and angles to be defined. Furthermore, Hilbert spaces are complete, which means that there are enough limits in the space to allow the techniques of calculus to be used. A Hilbert space is a special case of a Banach space.
Hilbert spaces were studied beginning in the first decade of the 20th century by David Hilbert, Erhard Schmidt, and Frigyes Riesz. They are indispensable tools in the theories of partial differential equations, quantum mechanics, Fourier analysis (which includes applications to signal processing and heat transfer), and ergodic theory (which forms the mathematical underpinning of thermodynamics). John von Neumann coined the term Hilbert space for the abstract concept that underlies many of these diverse applications. The success of Hilbert space methods ushered in a very fruitful era for functional analysis. Apart from the classical Euclidean vector spaces, examples of Hilbert spaces include spaces of square-integrable functions, spaces of sequences, Sobolev spaces consisting of generalized functions, and Hardy spaces of holomorphic functions.
Geometric intuition plays an important role in many aspects of Hilbert space theory. Exact analogs of the Pythagorean theorem and parallelogram law hold in a Hilbert space. At a deeper level, perpendicular projection onto a linear subspace plays a significant role in optimization problems and other aspects of the theory. An element of a Hilbert space can be uniquely specified by its coordinates with respect to an orthonormal basis, in analogy with Cartesian coordinates in classical geometry. When this basis is countably infinite, it allows identifying the Hilbert space with the space of the infinite sequences that are square-summable. The latter space is often in the older literature referred to as the Hilbert space.

## Definition and illustration
### Motivating example: Euclidean vector space

One of the most familiar examples of a Hilbert space is the Euclidean vector space consisting of three-dimensional vectors, denoted by R$^3$, and equipped with the dot product. The dot product takes two vectors x and y, and produces a real number x ⋅ y. If x and y are represented in Cartesian coordinates, then the dot product is defined by ${\begin{pmatrix}x_{1}\\x_{2}\\x_{3}\end{pmatrix}}\cdot {\begin{pmatrix}y_{1}\\y_{2}\\y_{3}\end{pmatrix}}=x_{1}y_{1}+x_{2}y_{2}+x_{3}y_{3}\,.$
The dot product satisfies the properties
* It is symmetric in x and y: x ⋅ y = y ⋅ x.
* It is linear in its first argument: (ax$_1$ + bx$_2$) ⋅ y = a(x$_1$ ⋅ y) + b(x$_2$ ⋅ y) for any scalars a, b, and vectors x$_1$, x$_2$, and y.
* It is positive definite: for all vectors x, x ⋅ x ≥ 0 , with equality if and only if x = 0.
An operation on pairs of vectors that, like the dot product, satisfies these three properties is known as a (real) inner product. A vector space equipped with such an inner product is known as a (real) inner product space. Every finite-dimensional inner product space is also a Hilbert space. The basic feature of the dot product that connects it with Euclidean geometry is that it is related to both the length (or norm) of a vector, denoted ‖x‖, and to the angle θ between two vectors x and y by means of the formula $\mathbf {x} \cdot \mathbf {y} =\left\|\mathbf {x} \right\|\left\|\mathbf {y} \right\|\,\cos \theta \,.$
Multivariable calculus in Euclidean space relies on the ability to compute limits, and to have useful criteria for concluding that limits exist. A mathematical series $\sum _{n=0}^{\infty }\mathbf {x} _{n}$ consisting of vectors in R$^3$ is absolutely convergent provided that the sum of the lengths converges as an ordinary series of real numbers: $\sum _{k=0}^{\infty }\|\mathbf {x} _{k}\|<\infty \,.$
Just as with a series of scalars, a series of vectors that converges absolutely also converges to some limit vector L in the Euclidean space, in the sense that ${\Biggl \|}\mathbf {L} -\sum _{k=0}^{N}\mathbf {x} _{k}{\Biggr \|}\to 0\quad {\text{as }}N\to \infty \,.$
This property expresses the completeness of Euclidean space: that a series that converges absolutely also converges in the ordinary sense.
Hilbert spaces are often taken over the complex numbers. The complex plane denoted by C is equipped with a notion of magnitude, the complex modulus |z|, which is defined as the square root of the product of z with its complex conjugate: $|z|^{2}=z{\overline {z}}\,.$
If z = x + iy is a decomposition of z into its real and imaginary parts, then the modulus is the usual Euclidean two-dimensional length: $|z|={\sqrt {x^{2}+y^{2}}}\,.$
The inner product of a pair of complex numbers z and w is the product of z with the complex conjugate of w: $\langle z,w\rangle =z{\overline {w}}\,.$
This is complex-valued. The real part of ⟨z, w⟩ gives the usual two-dimensional Euclidean dot product.
A second example is the space C$^2$ whose elements are pairs of complex numbers z = (z$_1$, z$_2$). Then an inner product of z with another such vector w = (w$_1$, w$_2$) is given by $\langle z,w\rangle =z_{1}{\overline {w_{1}}}+z_{2}{\overline {w_{2}}}\,.$
The real part of ⟨z, w⟩ is then the four-dimensional Euclidean dot product. This inner product is Hermitian symmetric, which means that the result of interchanging z and w is the complex conjugate: $\langle w,z\rangle ={\overline {\langle z,w\rangle }}\,.$

### Definition

A Hilbert space is a real or complex inner product space that is also a complete metric space with respect to the distance function induced by the inner product.
To say that a complex vector space H is a complex inner product space means that there is an inner product $\langle x,y\rangle$ associating a complex number to each pair of elements $x,y$ of H that satisfies the following properties:
* The inner product is conjugate symmetric; that is, the inner product of a pair of elements is equal to the complex conjugate of the inner product of the swapped elements: $\langle y,x\rangle ={\overline {\langle x,y\rangle }}\,.$ Importantly, this implies that $\langle x,x\rangle$ is a real number.
* The inner product is linear in its first argument. For all complex numbers $a$ and $b,$ $\langle ax_{1}+bx_{2},y\rangle =a\langle x_{1},y\rangle +b\langle x_{2},y\rangle \,.$
* The inner product of an element with itself is positive definite: ${\begin{alignedat}{4}\langle x,x\rangle >0&\quad {\text{ if }}x\neq 0,\\\langle x,x\rangle =0&\quad {\text{ if }}x=0\,.\end{alignedat}}$
It follows from properties 1 and 2 that a complex inner product is antilinear, also called conjugate linear, in its second argument, meaning that $\langle x,ay_{1}+by_{2}\rangle ={\bar {a}}\langle x,y_{1}\rangle +{\bar {b}}\langle x,y_{2}\rangle \,.$
A real inner product space is defined in the same way, except that H is a real vector space and the inner product takes real values. Such an inner product will be a bilinear map and $(H,H,\langle \cdot ,\cdot \rangle )$ will form a dual system.
The norm is the real-valued function $\|x\|={\sqrt {\langle x,x\rangle }}\,,$ and the distance $d$ between two points $x,y$ in H is defined in terms of the norm by $d(x,y)=\|x-y\|={\sqrt {\langle x-y,x-y\rangle }}\,.$
That this function is a distance function means firstly that it is symmetric in $x$ and $y,$ secondly that the distance between $x$ and itself is zero, and otherwise the distance between $x$ and $y$ must be positive, and lastly that the triangle inequality holds, meaning that the length of one leg of a triangle xyz cannot exceed the sum of the lengths of the other two legs: $d(x,z)\leq d(x,y)+d(y,z)\,.$
This last property is ultimately a consequence of the more fundamental Cauchy–Schwarz inequality, which asserts $\left|\langle x,y\rangle \right|\leq \|x\|\|y\|$ with equality if and only if $x$ and $y$ are linearly dependent.
With a distance function defined in this way, any inner product space is a metric space, and sometimes is known as a pre-Hilbert space. Any pre-Hilbert space that is additionally also a complete space is a Hilbert space.
The completeness of H is expressed using a form of the Cauchy criterion for sequences in H: a pre-Hilbert space H is complete if every Cauchy sequence converges with respect to this norm to an element in the space. Completeness can be characterized by the following equivalent condition: if a series of vectors $\sum _{k=0}^{\infty }u_{k}$ converges absolutely in the sense that $\sum _{k=0}^{\infty }\|u_{k}\|<\infty \,,$ then the series converges in H, in the sense that the partial sums converge to an element of H.
As a complete normed space, Hilbert spaces are by definition also Banach spaces. As such they are topological vector spaces, in which topological notions like the openness and closedness of subsets are well defined. Of special importance is the notion of a closed linear subspace of a Hilbert space that, with the inner product induced by restriction, is also complete (being a closed set in a complete metric space) and therefore a Hilbert space in its own right.

### Second example: sequence spaces

The sequence space l$^2$ consists of all infinite sequences z = (z$_1$, z$_2$, ...) of complex numbers such that the following series converges: $\sum _{n=1}^{\infty }|z_{n}|^{2}$
The inner product on l$^2$ is defined by: $\langle \mathbf {z} ,\mathbf {w} \rangle =\sum _{n=1}^{\infty }z_{n}{\overline {w_{n}}}\,,$
This second series converges as a consequence of the Cauchy–Schwarz inequality and the convergence of the previous series.
Completeness of the space holds provided that whenever a series of elements from l$^2$ converges absolutely (in norm), then it converges to an element of l$^2$. The proof is basic in mathematical analysis, and permits mathematical series of elements of the space to be manipulated with the same ease as series of complex numbers (or vectors in a finite-dimensional Euclidean space).

## History

Prior to the development of Hilbert spaces, other generalizations of Euclidean spaces were known to mathematicians and physicists. In particular, the idea of an abstract linear space (vector space) had gained some traction towards the end of the 19th century: this is a space whose elements can be added together and multiplied by scalars (such as real or complex numbers) without necessarily identifying these elements with "geometric" vectors, such as position and momentum vectors in physical systems. Other objects studied by mathematicians at the turn of the 20th century, in particular spaces of sequences (including series) and spaces of functions, can naturally be thought of as linear spaces. Functions, for instance, can be added together or multiplied by constant scalars, and these operations obey the algebraic laws satisfied by addition and scalar multiplication of spatial vectors.
In the first decade of the 20th century, parallel developments led to the introduction of Hilbert spaces. The first of these was the observation, which arose during David Hilbert and Erhard Schmidt's study of integral equations, that two square-integrable real-valued functions f and g on an interval [a, b] have an inner product
    $\langle f,g\rangle =\int _{a}^{b}f(x)g(x)\,\mathrm {d} x$
that has many of the familiar properties of the Euclidean dot product. In particular, the idea of an orthogonal family of functions has meaning. Schmidt exploited the similarity of this inner product with the usual dot product to prove an analog of the spectral decomposition for an operator of the form
    $f(x)\mapsto \int _{a}^{b}K(x,y)f(y)\,\mathrm {d} y$
where K is a continuous function symmetric in x and y. The resulting eigenfunction expansion expresses the function K as a series of the form
    $K(x,y)=\sum _{n}\lambda _{n}\varphi _{n}(x)\varphi _{n}(y)$
where the functions φ$_n$ are orthogonal in the sense that ⟨φ$_n$, φ$_m$⟩ = 0 for all n ≠ m. The individual terms in this series are sometimes referred to as elementary product solutions. However, there are eigenfunction expansions that fail to converge in a suitable sense to a square-integrable function: the missing ingredient, which ensures convergence, is completeness.
The second development was the Lebesgue integral, an alternative to the Riemann integral introduced by Henri Lebesgue in 1904. The Lebesgue integral made it possible to integrate a much broader class of functions. In 1907, Frigyes Riesz and Ernst Sigismund Fischer independently proved that the space L$^2$ of square Lebesgue-integrable functions is a complete metric space. As a consequence of the interplay between geometry and completeness, the 19th century results of Joseph Fourier, Friedrich Bessel and Marc-Antoine Parseval on trigonometric series easily carried over to these more general spaces, resulting in a geometrical and analytical apparatus now usually known as the Riesz–Fischer theorem.
Further basic results were proved in the early 20th century. For example, the Riesz representation theorem was independently established by Maurice Fréchet and Frigyes Riesz in 1907. John von Neumann coined the term abstract Hilbert space in his work on unbounded Hermitian operators. Although other mathematicians such as Hermann Weyl and Norbert Wiener had already studied particular Hilbert spaces in great detail, often from a physically motivated point of view, von Neumann gave the first complete and axiomatic treatment of them. Von Neumann later used them in his seminal work on the foundations of quantum mechanics, and in his continued work with Eugene Wigner. The name "Hilbert space" was soon adopted by others, for example by Hermann Weyl in his book on quantum mechanics and the theory of groups.
The significance of the concept of a Hilbert space was underlined with the realization that it offers one of the best mathematical formulations of quantum mechanics. In short, the states of a quantum mechanical system are vectors in a certain Hilbert space, the observables are hermitian operators on that space, the symmetries of the system are unitary operators, and measurements are orthogonal projections. The relation between quantum mechanical symmetries and unitary operators provided an impetus for the development of the unitary representation theory of groups, initiated in the 1928 work of Hermann Weyl. On the other hand, in the early 1930s it became clear that classical mechanics can be described in terms of Hilbert space (Koopman–von Neumann classical mechanics) and that certain properties of classical dynamical systems can be analyzed using Hilbert space techniques in the framework of ergodic theory.
The algebra of observables in quantum mechanics is naturally an algebra of operators defined on a Hilbert space, according to Werner Heisenberg's matrix mechanics formulation of quantum theory. Von Neumann began investigating operator algebras in the 1930s, as rings of operators on a Hilbert space. The kind of algebras studied by von Neumann and his contemporaries are now known as von Neumann algebras. In the 1940s, Israel Gelfand, Mark Naimark and Irving Segal gave a definition of a kind of operator algebras called C*-algebras that on the one hand made no reference to an underlying Hilbert space, and on the other extrapolated many of the useful features of the operator algebras that had previously been studied. The spectral theorem for self-adjoint operators in particular that underlies much of the existing Hilbert space theory was generalized to C*-algebras. These techniques are now basic in abstract harmonic analysis and representation theory.

## Examples
### Lebesgue spaces

Lebesgue spaces are function spaces associated to measure spaces (X, M, μ), where X is a set, M is a σ-algebra of subsets of X, and μ is a countably additive measure on M. Let L$^2$(X, μ) be the space of those complex-valued measurable functions on X for which the Lebesgue integral of the square of the absolute value of the function is finite, i.e., for a function f in L$^2$(X, μ), $\int _{X}|f|^{2}\mathrm {d} \mu <\infty \,,$ and where functions are identified if and only if they differ only on a set of measure zero.
The inner product of functions f and g in L$^2$(X, μ) is then defined as $\langle f,g\rangle =\int _{X}f(t){\overline {g(t)}}\,\mathrm {d} \mu (t)$ or $\langle f,g\rangle =\int _{X}{\overline {f(t)}}g(t)\,\mathrm {d} \mu (t)\,,$
where the second form (conjugation of the first element) is commonly found in the theoretical physics literature. For f and g in L$^2$, the integral exists because of the Cauchy–Schwarz inequality, and defines an inner product on the space. Equipped with this inner product, L$^2$ is in fact complete. The Lebesgue integral is essential to ensure completeness: on domains of real numbers, for instance, not enough functions are Riemann integrable.
The Lebesgue spaces appear in many natural settings. The spaces L$^2$(R) and L$^2$([0,1]) of square-integrable functions with respect to the Lebesgue measure on the real line and unit interval, respectively, are natural domains on which to define the Fourier transform and Fourier series. In other situations, the measure may be something other than the ordinary Lebesgue measure on the real line. For instance, if w is any positive measurable function, the space of all measurable functions f on the interval [0, 1] satisfying $\int _{0}^{1}{\bigl |}f(t){\bigr |}^{2}w(t)\,\mathrm {d} t<\infty$ is called the weighted L$^2$ space L$^2$ $_w$([0, 1]), and w is called the weight function. The inner product is defined by $\langle f,g\rangle =\int _{0}^{1}f(t){\overline {g(t)}}w(t)\,\mathrm {d} t\,.$
The weighted space L$^2$ $_w$([0, 1]) is identical with the Hilbert space L$^2$([0, 1], μ) where the measure μ of a Lebesgue-measurable set A is defined by $\mu (A)=\int _{A}w(t)\,\mathrm {d} t\,.$
Weighted L$^2$ spaces like this are frequently used to study orthogonal polynomials, because different families of orthogonal polynomials are orthogonal with respect to different weighting functions.

### Sobolev spaces

Sobolev spaces, denoted by H$^s$ or W$^{s,2}$, are Hilbert spaces. These are a special kind of function space in which differentiation may be performed, but that (unlike other Banach spaces such as the Hölder spaces) support the structure of an inner product. Because differentiation is permitted, Sobolev spaces are a convenient setting for the theory of partial differential equations. They also form the basis of the theory of direct methods in the calculus of variations.
For s a non-negative integer and Ω ⊂ R$^n$, the Sobolev space H$^s$(Ω) contains L$^2$ functions whose weak derivatives of order up to s are also L$^2$. The inner product in H$^s$(Ω) is $\langle f,g\rangle =\int _{\Omega }f(x){\bar {g}}(x)\,\mathrm {d} x+\int _{\Omega }Df(x)\cdot D{\bar {g}}(x)\,\mathrm {d} x+\cdots +\int _{\Omega }D^{s}f(x)\cdot D^{s}{\bar {g}}(x)\,\mathrm {d} x$ where the dot indicates the dot product in the Euclidean space of partial derivatives of each order. Sobolev spaces can also be defined when s is not an integer.
Sobolev spaces are also studied from the point of view of spectral theory, relying more specifically on the Hilbert space structure. If Ω is a suitable domain, then one can define the Sobolev space H$^s$(Ω) as the space of Bessel potentials; roughly, $H^{s}(\Omega )=\left\{(1-\Delta )^{-s/2}f\mathrel {\Big |} f\in L^{2}(\Omega )\right\}\,.$
Here Δ is the Laplacian and (1 − Δ)$^{−s/2}$ is understood in terms of the spectral mapping theorem. Apart from providing a workable definition of Sobolev spaces for non-integer s, this definition also has particularly desirable properties under the Fourier transform that make it ideal for the study of pseudodifferential operators. Using these methods on a compact Riemannian manifold, one can obtain for instance the Hodge decomposition, which is the basis of Hodge theory.

### Spaces of holomorphic functions
#### Hardy spaces

The Hardy spaces are function spaces, arising in complex analysis and harmonic analysis, whose elements are certain holomorphic functions in a complex domain. Let U denote the unit disc in the complex plane. Then the Hardy space H$^2$(U) is defined as the space of holomorphic functions f on U such that the means $M_{r}(f)={\frac {1}{2\pi }}\int _{0}^{2\pi }\left|f{\biglre^{i\theta }{\bigr )}\right|^{2}\,\mathrm {d} \theta$ remain bounded for r < 1. The norm on this Hardy space is defined by $\left\|f\right\|_{2}=\lim _{r\to 1}{\sqrt {M_{r}(f)}}\,.$
Hardy spaces in the disc are related to Fourier series. A function f is in H$^2$(U) if and only if $f(z)=\sum _{n=0}^{\infty }a_{n}z^{n}$ where $\sum _{n=0}^{\infty }|a_{n}|^{2}<\infty \,.$
Thus H$^2$(U) consists of those functions that are L$^2$ on the circle, and whose negative frequency Fourier coefficients vanish.

#### Bergman spaces

The Bergman spaces are another family of Hilbert spaces of holomorphic functions. Let D be a bounded open set in the complex plane (or a higher-dimensional complex space) and let L$^{2,h}$(D) be the space of holomorphic functions f in D that are also in L$^2$(D) in the sense that $\|f\|^{2}=\int _{D}|f(z)|^{2}\,\mathrm {d} \mu (z)<\infty \,,$ where the integral is taken with respect to the Lebesgue measure in D. Clearly L$^{2,h}$(D) is a subspace of L$^2$(D); in fact, it is a closed subspace, and so a Hilbert space in its own right. This is a consequence of the estimate, valid on compact subsets K of D, that $\sup _{z\in K}\left|f(z)\right|\leq C_{K}\left\|f\right\|_{2}\,,$ which in turn follows from Cauchy's integral formula. Thus convergence of a sequence of holomorphic functions in L$^2$(D) implies also compact convergence, and so the limit function is also holomorphic. Another consequence of this inequality is that the linear functional that evaluates a function f at a point of D is actually continuous on L$^{2,h}$(D). The Riesz representation theorem implies that the evaluation functional can be represented as an element of L$^{2,h}$(D). Thus, for every z ∈ D, there is a function η$_z$ ∈ L$^{2,h}$(D) such that $f(z)=\int _{D}f(\zeta ){\overline {\eta _{z}(\zeta )}}\,\mathrm {d} \mu (\zeta )$ for all f ∈ L$^{2,h}$(D). The integrand $K(\zeta ,z)={\overline {\eta _{z}(\zeta )}}$ is known as the Bergman kernel of D. This integral kernel satisfies a reproducing property $f(z)=\int _{D}f(\zeta )K(\zeta ,z)\,\mathrm {d} \mu (\zeta )\,.$
A Bergman space is an example of a reproducing kernel Hilbert space, which is a Hilbert space of functions along with a kernel K(ζ, z) that verifies a reproducing property analogous to this one. The Hardy space H$^2$(D) also admits a reproducing kernel, known as the Szegő kernel. Reproducing kernels are common in other areas of mathematics as well. For instance, in harmonic analysis the Poisson kernel is a reproducing kernel for the Hilbert space of square-integrable harmonic functions in the unit ball. That the latter is a Hilbert space at all is a consequence of the mean value theorem for harmonic functions.

## Applications

Many of the applications of Hilbert spaces exploit the fact that Hilbert spaces support generalizations of simple geometric concepts like projection and change of basis from their usual finite dimensional setting. In particular, the spectral theory of continuous self-adjoint linear operators on a Hilbert space generalizes the usual spectral decomposition of a matrix, and this often plays a major role in applications of the theory to other areas of mathematics and physics.

### Sturm–Liouville theory

In the theory of ordinary differential equations, spectral methods on a suitable Hilbert space are used to study the behavior of eigenvalues and eigenfunctions of differential equations. For example, the Sturm–Liouville problem arises in the study of the harmonics of waves in a violin string or a drum, and is a central problem in ordinary differential equations. The problem is a differential equation of the form $-{\frac {\mathrm {d} }{\mathrm {d} x}}\left[p(x){\frac {\mathrm {d} y}{\mathrm {d} x}}\right]+q(x)y=\lambda w(x)y$ for an unknown function y on an interval [a, b], satisfying general homogeneous Robin boundary conditions ${\begin{cases}\alpha y(a)+\alpha 'y'(a)&=0\\\beta y(b)+\beta 'y'(b)&=0\,.\end{cases}}$ The functions p, q, and w are given in advance, and the problem is to find the function y and constants λ for which the equation has a solution. The problem only has solutions for certain values of λ, called eigenvalues of the system, and this is a consequence of the spectral theorem for compact operators applied to the integral operator defined by the Green's function for the system. Furthermore, another consequence of this general result is that the eigenvalues λ of the system can be arranged in an increasing sequence tending to infinity.

### Partial differential equations

Hilbert spaces form a basic tool in the study of partial differential equations. For many classes of partial differential equations, such as linear elliptic equations, it is possible to consider a generalized solution (known as a weak solution) by enlarging the class of functions. Many weak formulations involve the class of Sobolev functions, which is a Hilbert space. A suitable weak formulation reduces to a geometrical problem, the analytic problem of finding a solution or, often what is more important, showing that a solution exists and is unique for given boundary data. For linear elliptic equations, one geometrical result that ensures unique solvability for a large class of problems is the Lax–Milgram theorem. This strategy forms the rudiment of the Galerkin method (a finite element method) for numerical solution of partial differential equations.
A typical example is the Poisson equation −Δu = g with Dirichlet boundary conditions in a bounded domain Ω in R$^2$. The weak formulation consists of finding a function u such that, for all continuously differentiable functions v in Ω vanishing on the boundary: $\int _{\Omega }\nabla u\cdot \nabla v=\int _{\Omega }gv\,.$
This can be recast in terms of the Hilbert space H$^1$ $_0$(Ω) consisting of functions u such that u, along with its weak partial derivatives, are square integrable on Ω, and vanish on the boundary. The question then reduces to finding u in this space such that for all v in this space $a(u,v)=b(v)$
where a is a continuous bilinear form, and b is a continuous linear functional, given respectively by $a(u,v)=\int _{\Omega }\nabla u\cdot \nabla v,\quad b(v)=\int _{\Omega }gv\,.$
Since the Poisson equation is elliptic, it follows from Poincaré's inequality that the bilinear form a is coercive. The Lax–Milgram theorem then ensures the existence and uniqueness of solutions of this equation.
Hilbert spaces allow for many elliptic partial differential equations to be formulated in a similar way, and the Lax–Milgram theorem is then a basic tool in their analysis. With suitable modifications, similar techniques can be applied to parabolic partial differential equations and certain hyperbolic partial differential equations.

### Ergodic theory

The field of ergodic theory is the study of the long-term behavior of chaotic dynamical systems. The protypical case of a field that ergodic theory applies to is thermodynamics, in which—though the microscopic state of a system is extremely complicated (it is impossible to understand the ensemble of individual collisions between particles of matter)—the average behavior over sufficiently long time intervals is tractable. The laws of thermodynamics are assertions about such average behavior. In particular, one formulation of the zeroth law of thermodynamics asserts that over sufficiently long timescales, the only functionally independent measurement that one can make of a thermodynamic system in equilibrium is its total energy, in the form of temperature.
An ergodic dynamical system is one for which, apart from the energy—measured by the Hamiltonian—there are no other functionally independent conserved quantities on the phase space. More explicitly, suppose that the energy E is fixed, and let Ω$_E$ be the subset of the phase space consisting of all states of energy E (an energy surface), and let T$_t$ denote the evolution operator on the phase space. The dynamical system is ergodic if every invariant measurable functions on Ω$_E$ is constant almost everywhere. An invariant function f is one for which $f(T_{t}w)=f(w)$ for all w on Ω$_E$ and all time t. Liouville's theorem implies that there exists a measure μ on the energy surface that is invariant under the time translation. As a result, time translation is a unitary transformation of the Hilbert space L$^2$(Ω$_E$, μ) consisting of square-integrable functions on the energy surface Ω$_E$ with respect to the inner product $\left\langle f,g\right\rangle _{L^{2}\left(\Omega _{E},\mu \right)}=\int _{E}f{\bar {g}}\,\mathrm {d} \mu \,.$
The von Neumann mean ergodic theorem states the following:
* If U$_t$ is a (strongly continuous) one-parameter semigroup of unitary operators on a Hilbert space H, and P is the orthogonal projection onto the space of common fixed points of U$_t$, {x ∈H | U$_t$x = x, ∀t > 0}, then $Px=\lim _{T\to \infty }{\frac {1}{T}}\int _{0}^{T}U_{t}x\,\mathrm {d} t\,.$
For an ergodic system, the fixed set of the time evolution consists only of the constant functions, so the ergodic theorem implies the following: for any function f ∈ L$^2$(Ω$_E$, μ), ${\underset {T\to \infty }{L^{2}-\lim }}{\frac {1}{T}}\int _{0}^{T}f(T_{t}w)\,\mathrm {d} t=\int _{\Omega _{E}}f(y)\,\mathrm {d} \mu (y)\,.$
That is, the long time average of an observable f is equal to its expectation value over an energy surface.

### Fourier analysis

One of the basic goals of Fourier analysis is to decompose a function into a (possibly infinite) linear combination of given basis functions: the associated Fourier series. The classical Fourier series associated to a function f defined on the interval [0, 1] is a series of the form $\sum _{n=-\infty }^{\infty }a_{n}e^{2\pi in\theta }$ where $a_{n}=\int _{0}^{1}f(\theta )\;\!e^{-2\pi in\theta }\,\mathrm {d} \theta \,.$
The example of adding up the first few terms in a Fourier series for a sawtooth function is shown in the figure. The basis functions are sine waves with wavelengths ⁠λ/n⁠ (for integer n) shorter than the wavelength λ of the sawtooth itself (except for n = 1, the fundamental wave).
A significant problem in classical Fourier series asks in what sense the Fourier series converges, if at all, to the function f. Hilbert space methods provide one possible answer to this question. The functions e$_n$(θ) = e$^{2πinθ}$ form an orthogonal basis of the Hilbert space L$^2$([0, 1]). Consequently, any square-integrable function can be expressed as a series $f(\theta )=\sum _{n}a_{n}e_{n}(\theta )\,,\quad a_{n}=\langle f,e_{n}\rangle$ and, moreover, this series converges in the Hilbert space sense (that is, in the L$^2$ mean).
The problem can also be studied from the abstract point of view: every Hilbert space has an orthonormal basis, and every element of the Hilbert space can be written in a unique way as a sum of multiples of these basis elements. The coefficients appearing on these basis elements are sometimes known abstractly as the Fourier coefficients of the element of the space. The abstraction is especially useful when it is more natural to use different basis functions for a space such as L$^2$([0, 1]). In many circumstances, it is desirable not to decompose a function into trigonometric functions, but rather into orthogonal polynomials or wavelets for instance, and in higher dimensions into spherical harmonics.
For instance, if e$_n$ are any orthonormal basis functions of L$^2$[0, 1], then a given function in L$^2$[0, 1] can be approximated as a finite linear combination $f(x)\approx f_{n}(x)=a_{1}e_{1}(x)+a_{2}e_{2}(x)+\cdots +a_{n}e_{n}(x)\,.$
The coefficients {a$_j$} are selected to make the magnitude of the difference ‖f − f$_n$‖$^2$ as small as possible. Geometrically, the best approximation is the orthogonal projection of f onto the subspace consisting of all linear combinations of the {e$_j$}, and can be calculated by $a_{j}=\int _{0}^{1}{\overline {e_{j}(x)}}f(x)\,\mathrm {d} x\,.$
That this formula minimizes the difference ‖f − f$_n$‖$^2$ is a consequence of Bessel's inequality and Parseval's formula.
In various applications to physical problems, a function can be decomposed into physically meaningful eigenfunctions of a differential operator (typically the Laplace operator): this forms the foundation for the spectral study of functions, in reference to the spectrum of the differential operator. A concrete physical application involves the problem of hearing the shape of a drum: given the fundamental modes of vibration that a drumhead is capable of producing, can one infer the shape of the drum itself? The mathematical formulation of this question involves the Dirichlet eigenvalues of the Laplace equation in the plane, that represent the fundamental modes of vibration in direct analogy with the integers that represent the fundamental modes of vibration of the violin string.
Spectral theory also underlies certain aspects of the Fourier transform of a function. Whereas Fourier analysis decomposes a function defined on a compact set into the discrete spectrum of the Laplacian (which corresponds to the vibrations of a violin string or drum), the Fourier transform of a function is the decomposition of a function defined on all of Euclidean space into its components in the continuous spectrum of the Laplacian. The Fourier transformation is also geometrical, in a sense made precise by the Plancherel theorem, that asserts that it is an isometry of one Hilbert space (the "time domain") with another (the "frequency domain"). This isometry property of the Fourier transformation is a recurring theme in abstract harmonic analysis (since it reflects the conservation of energy for the continuous Fourier Transform), as evidenced for instance by the Plancherel theorem for spherical functions occurring in noncommutative harmonic analysis.

### Quantum mechanics

In the mathematically rigorous formulation of quantum mechanics, developed by John von Neumann, the possible states (more precisely, the pure states) of a quantum mechanical system are represented by unit vectors (called state vectors) residing in a complex separable Hilbert space, known as the state space, well defined up to a complex number of norm 1 (the phase factor). In other words, the possible states are points in the projectivization of a Hilbert space, usually called the complex projective space. The exact nature of this Hilbert space is dependent on the system; for example, the position and momentum states for a single non-relativistic spin zero particle is the space of all square-integrable functions, while the states for the spin of a single proton are unit elements of the two-dimensional complex Hilbert space of spinors. Each observable is represented by a self-adjoint linear operator acting on the state space. Each eigenstate of an observable corresponds to an eigenvector of the operator, and the associated eigenvalue corresponds to the value of the observable in that eigenstate.
The inner product between two state vectors is a complex number known as a probability amplitude. During an ideal measurement of a quantum mechanical system, the probability that a system collapses from a given initial state to a particular eigenstate is given by the square of the absolute value of the probability amplitudes between the initial and final states. The possible results of a measurement are the eigenvalues of the operator—which explains the choice of self-adjoint operators, for all the eigenvalues must be real. The probability distribution of an observable in a given state can be found by computing the spectral decomposition of the corresponding operator.
For a general system, states are typically not pure, but instead are represented as statistical mixtures of pure states, or mixed states, given by density matrices: self-adjoint operators of trace one on a Hilbert space. Moreover, for general quantum mechanical systems, the effects of a single measurement can influence other parts of a system in a manner that is described instead by a positive operator valued measure. Thus the structure both of the states and observables in the general theory is considerably more complicated than the idealization for pure states.

### Probability theory

In probability theory, Hilbert spaces also have diverse applications. Here a fundamental Hilbert space is the space of random variables on a given probability space, having class $L^{2}$ (finite first and second moments). A common operation in statistics is that of centering a random variable by subtracting its expectation. Thus if $X$ is a random variable, then $X-E(X)$ is its centering. In the Hilbert space view, this is the orthogonal projection of $X$ onto the kernel of the expectation operator, which a continuous linear functional on the Hilbert space (in fact, the inner product with the constant random variable 1), and so this kernel is a closed subspace.
The conditional expectation has a natural interpretation in the Hilbert space. Suppose that a probability space $(\Omega ,P,{\mathcal {B}})$ is given, where ${\mathcal {B}}$ is a sigma algebra on the set $\Omega$, and $P$ is a probability measure on the measure space $(\Omega ,{\mathcal {B}})$. If ${\mathcal {F}}\leq {\mathcal {B}}$ is a sigma subalgebra of ${\mathcal {B}}$, then the conditional expectation $E[X|{\mathcal {F}}]$ is the orthogonal projection of $X$ onto the subspace of $L^{2}(\Omega ,P)$ consisting of the ${\mathcal {F}}$-measurable functions. If the random variable $X$ in $L^{2}(\Omega ,P)$ is independent of the sigma algebra ${\mathcal {F}}$ then conditional expectation $E(X|{\mathcal {F}})=E(X)$, i.e., its projection onto the ${\mathcal {F}}$-measurable functions is constant. Equivalently, the projection of its centering is zero.
In particular, if two random variables $X$ and $Y$ (in $L^{2}(\Omega ,P)$) are independent, then the centered random variables $X-E(X)$ and $Y-E(Y)$ are orthogonal. (This means that the two variables have zero covariance: they are uncorrelated.) In that case, the Pythagorean theorem in the kernel of the expectation operator implies that the variances of $X$ and $Y$ satisfy the identity: $\operatorname {Var} (X+Y)=\operatorname {Var} (X)+\operatorname {Var} (Y),$ sometimes called the Pythagorean theorem of statistics, and is of importance in linear regression. As Stapleton (1995) puts it, "the analysis of variance may be viewed as the decomposition of the squared length of a vector into the sum of the squared lengths of several vectors, using the Pythagorean Theorem."
The theory of martingales can be formulated in Hilbert spaces. A martingale in a Hilbert space is a sequence $x_{1},x_{2},\dots$ of elements of a Hilbert space such that, for each n, $x_{n}$ is the orthogonal projection of $x_{n+1}$ onto the linear hull of $x_{1},\dots ,x_{n}$. If the $x_{k}$ are random variables, this reproduces the usual definition of a (discrete) martingale: the expectation of $x_{n+1}$, conditioned on $x_{1},\dots ,x_{n}$, is equal to $x_{n}$.
Hilbert spaces are also used throughout the foundations of the Itô calculus. To any square-integrable martingale, it is possible to associate a Hilbert norm on the space of equivalence classes of progressively measurable processes with respect to the martingale (using the quadratic variation of the martingale as the measure). The Itô integral can be constructed by first defining it for simple processes, and then exploiting their density in the Hilbert space. A noteworthy result is then the Itô isometry, which attests that for any martingale M having quadratic variation measure $d\langle M\rangle _{t}$, and any progressively measurable process H: $E\left[\left(\int _{0}^{t}H_{s}dM_{s}\right)^{2}\right]=E\left[\int _{0}^{t}H_{s}^{2}d\langle M\rangle _{s}\right]$ whenever the expectation on the right-hand side is finite.
A deeper application of Hilbert spaces that is especially important in the theory of Gaussian processes is an attempt, due to Leonard Gross and others, to make sense of certain formal integrals over infinite dimensional spaces like the Feynman path integral from quantum field theory. The problem with integral like this is that there is no infinite dimensional Lebesgue measure. The notion of an abstract Wiener space allows one to construct a measure on a Banach space B that contains a Hilbert space H, called the Cameron–Martin space, as a dense subset, out of a finitely additive cylinder set measure on H. The resulting measure on B is countably additive and invariant under translation by elements of H, and this provides a mathematically rigorous way of thinking of the Wiener measure as a Gaussian measure on the Sobolev space $H^{1}([0,\infty ))$.

### Color perception

Any true physical color can be represented by a combination of pure spectral colors. As physical colors can be composed of any number of spectral colors, the space of physical colors may aptly be represented by a Hilbert space over spectral colors. Humans have three types of cone cells for color perception, so the perceivable colors can be represented by 3-dimensional Euclidean space. The many-to-one linear mapping from the Hilbert space of physical colors to the Euclidean space of human perceivable colors explains why many distinct physical colors may be perceived by humans to be identical (e.g., pure yellow light versus a mix of red and green light, see Metamerism).

## Properties
### Pythagorean identity

Two vectors u and v in a Hilbert space H are orthogonal when ⟨u, v⟩ = 0. The notation for this is u ⊥ v. More generally, when S is a subset in H, the notation u ⊥ S means that u is orthogonal to every element from S.
When u and v are orthogonal, one has $\|u+v\|^{2}=\langle u+v,u+v\rangle =\langle u,u\rangle +2\,\operatorname {Re} \langle u,v\rangle +\langle v,v\rangle =\|u\|^{2}+\|v\|^{2}\,.$
By induction on n, this is extended to any family u$_1$, ..., u$_n$ of n orthogonal vectors, $\left\|u_{1}+\cdots +u_{n}\right\|^{2}=\left\|u_{1}\right\|^{2}+\cdots +\left\|u_{n}\right\|^{2}.$
Whereas the Pythagorean identity as stated is valid in any inner product space, completeness is required for the extension of the Pythagorean identity to series. A series Σu$_k$ of orthogonal vectors converges in H if and only if the series of squares of norms converges, and ${\Biggl \|}\sum _{k=0}^{\infty }u_{k}{\Biggr \|}^{2}=\sum _{k=0}^{\infty }\left\|u_{k}\right\|^{2}\,.$ Furthermore, the sum of a series of orthogonal vectors is independent of the order in which it is taken.

### Parallelogram identity and polarization

By definition, every Hilbert space is also a Banach space. Furthermore, in every Hilbert space the following parallelogram identity holds: $\|u+v\|^{2}+\|u-v\|^{2}=2{\bigl\|u\|^{2}+\|v\|^{2}{\bigr )}\,.$
Conversely, every Banach space in which the parallelogram identity holds is a Hilbert space, and the inner product is uniquely determined by the norm by the polarization identity. For real Hilbert spaces, the polarization identity is $\langle u,v\rangle ={\tfrac {1}{4}}{\bigl\|u+v\|^{2}-\|u-v\|^{2}{\bigr )}\,.$
For complex Hilbert spaces, it is $\langle u,v\rangle ={\tfrac {1}{4}}{\bigl\|u+v\|^{2}-\|u-v\|^{2}+i\|u+iv\|^{2}-i\|u-iv\|^{2}{\bigr )}\,.$
The parallelogram law implies that any Hilbert space is a uniformly convex Banach space.

### Best approximation

This subsection employs the Hilbert projection theorem. If C is a non-empty closed convex subset of a Hilbert space H and x a point in H, there exists a unique point y ∈ C that minimizes the distance between x and points in C, $y\in C\,,\quad \|x-y\|=\operatorname {dist} (x,C)=\min {\bigl \\|x-z\|\mathrel {\big |} z\in C{\bigr \}}\,.$
This is equivalent to saying that there is a point with minimal norm in the translated convex set D = C − x. The proof consists in showing that every minimizing sequence (d$_n$) ⊂ D is Cauchy (using the parallelogram identity) hence converges (using completeness) to a point in D that has minimal norm. More generally, this holds in any uniformly convex Banach space.
When this result is applied to a closed subspace F of H, it can be shown that the point y ∈ F closest to x is characterized by $y\in F\,,\quad x-y\perp F\,.$
This point y is the orthogonal projection of x onto F, and the mapping P$_F$ : x → y is linear (see § Orthogonal complements and projections). This result is especially significant in applied mathematics, especially numerical analysis, where it forms the basis of least squares methods.
In particular, when F is not equal to H, one can find a nonzero vector v orthogonal to F (select x ∉ F and v = x − y). A very useful criterion is obtained by applying this observation to the closed subspace F generated by a subset S of H.
    A subset S of H spans a dense vector subspace if (and only if) the vector 0 is the sole vector v ∈ H orthogonal to S.

### Duality

The dual space H* is the space of all continuous linear functions from the space H into the base field. It carries a natural norm, defined by $\|\varphi \|=\sup _{\|x\|=1,x\in H}|\varphi (x)|\,.$ This norm satisfies the parallelogram law, and so the dual space is also an inner product space where this inner product can be defined in terms of this dual norm by using the polarization identity. The dual space is also complete so it is a Hilbert space in its own right. If e$_•$ = (e$_i$)$_{i∈I}$ is a complete orthonormal basis for H then the inner product on the dual space of any two $f,g\in H^{*}$ is $\langle f,g\rangle _{H^{*}}=\sum _{i\in I}f(e_{i}){\overline {g(e_{i})}}$ where all but countably many of the terms in this series are zero.
The Riesz representation theorem affords a convenient description of the dual space. To every element u of H, there is a unique element φ$_u$ of H*, defined by $\varphi _{u}(x)=\langle x,u\rangle$ where moreover, $\left\|\varphi _{u}\right\|=\left\|u\right\|.$
The Riesz representation theorem states that the map from H to H* defined by u ↦ φ$_u$ is surjective, which makes this map an isometric antilinear isomorphism. So to every element φ of the dual H* there exists one and only one u$_φ$ in H such that $\langle x,u_{\varphi }\rangle =\varphi (x)$ for all x ∈ H. The inner product on the dual space H* satisfies $\langle \varphi ,\psi \rangle =\langle u_{\psi },u_{\varphi }\rangle \,.$
The reversal of order on the right-hand side restores linearity in φ from the antilinearity of u$_φ$. In the real case, the antilinear isomorphism from H to its dual is actually an isomorphism, and so real Hilbert spaces are naturally isomorphic to their own duals.
The representing vector u$_φ$ is obtained in the following way. When φ ≠ 0, the kernel F = Ker(φ) is a closed vector subspace of H, not equal to H, hence there exists a nonzero vector v orthogonal to F. The vector u is a suitable scalar multiple λv of v. The requirement that φ(v) = ⟨v, u⟩ yields $u=\langle v,v\rangle ^{-1}\,{\overline {\varphi (v)}}\,v\,.$
This correspondence φ ↔ u is exploited by the bra–ket notation popular in physics. It is common in physics to assume that the inner product, denoted by ⟨x|y⟩, is linear on the right, $\langle x|y\rangle =\langle y,x\rangle \,.$ The result ⟨x|y⟩ can be seen as the action of the linear functional ⟨x| (the bra) on the vector |y⟩ (the ket).
The Riesz representation theorem relies fundamentally not just on the presence of an inner product, but also on the completeness of the space. In fact, the theorem implies that the topological dual of any inner product space can be identified with its completion. An immediate consequence of the Riesz representation theorem is also that a Hilbert space H is reflexive, meaning that the natural map from H into its double dual space is an isomorphism.

### Weakly convergent sequences

In a Hilbert space H, a sequence {x$_n$} is weakly convergent to a vector x ∈ H when $\lim _{n}\langle x_{n},v\rangle =\langle x,v\rangle$ for every v ∈ H.
For example, any orthonormal sequence {f$_n$} converges weakly to 0, as a consequence of Bessel's inequality. Every weakly convergent sequence {x$_n$} is bounded, by the uniform boundedness principle.
Conversely, every bounded sequence in a Hilbert space admits weakly convergent subsequences (Alaoglu's theorem). This fact may be used to prove minimization results for continuous convex functionals, in the same way that the Bolzano–Weierstrass theorem is used for continuous functions on R$^d$. Among several variants, one simple statement is as follows:
    If f : H → R is a convex continuous function such that f(x) tends to +∞ when ‖x‖ tends to ∞, then f admits a minimum at some point x$_0$ ∈ H.
This fact (and its various generalizations) are fundamental for direct methods in the calculus of variations. Minimization results for convex functionals are also a direct consequence of the slightly more abstract fact that closed bounded convex subsets in a Hilbert space H are weakly compact, since H is reflexive. The existence of weakly convergent subsequences is a special case of the Eberlein–Šmulian theorem.

### Banach space properties

Any general property of Banach spaces continues to hold for Hilbert spaces. The open mapping theorem states that a continuous surjective linear transformation from one Banach space to another is an open mapping meaning that it sends open sets to open sets. A corollary is the bounded inverse theorem, that a continuous and bijective linear function from one Banach space to another is an isomorphism (that is, a continuous linear map whose inverse is also continuous). This theorem is considerably simpler to prove in the case of Hilbert spaces than in general Banach spaces. The open mapping theorem is equivalent to the closed graph theorem, which asserts that a linear function from one Banach space to another is continuous if and only if its graph is a closed set. In the case of Hilbert spaces, this is basic in the study of unbounded operators (see Closed operator).
The (geometrical) Hahn–Banach theorem asserts that a closed convex set can be separated from any point outside it by means of a hyperplane of the Hilbert space. This is an immediate consequence of the best approximation property: if y is the element of a closed convex set F closest to x, then the separating hyperplane is the plane perpendicular to the segment xy passing through its midpoint.

## Operators on Hilbert spaces
### Bounded operators

The continuous linear operators A : H$_1$ → H$_2$ from a Hilbert space H$_1$ to a second Hilbert space H$_2$ are bounded in the sense that they map bounded sets to bounded sets. Conversely, if an operator is bounded, then it is continuous. The space of such bounded linear operators has a norm, the operator norm given by $\lVert A\rVert =\sup {\bigl \\|Ax\|\mathrel {\big |} \|x\|\leq 1{\bigr \}}\,.$
The sum and the composite of two bounded linear operators is again bounded and linear. For y in H$_2$, the map that sends x ∈ H$_1$ to ⟨Ax, y⟩ is linear and continuous, and according to the Riesz representation theorem can therefore be represented in the form $\left\langle x,A^{*}y\right\rangle =\langle Ax,y\rangle$ for some vector A*y in H$_1$. This defines another bounded linear operator A* : H$_2$ → H$_1$, the adjoint of A. The adjoint satisfies A** = A. When the Riesz representation theorem is used to identify each Hilbert space with its continuous dual space, the adjoint of A can be shown to be identical to the transpose $^t$A : H$_2$* → H$_1$* of A, which by definition sends $\psi \in H_{2}^{*}$ to the functional $\psi \circ A\in H_{1}^{*}.$
The set B(H) of all bounded linear operators on H (meaning operators H → H), together with the addition and composition operations, the norm and the adjoint operation, is a C*-algebra, which is a type of operator algebra.
An element A of B(H) is called 'self-adjoint' or 'Hermitian' if A* = A. If A is Hermitian and ⟨Ax, x⟩ ≥ 0 for every x, then A is called 'nonnegative', written A ≥ 0; if equality holds only when x = 0, then A is called 'positive'. The set of self adjoint operators admits a partial order, in which A ≥ B if A − B ≥ 0. If A has the form B*B for some B, then A is nonnegative; if B is invertible, then A is positive. A converse is also true in the sense that, for a non-negative operator A, there exists a unique non-negative square root B such that $A=B^{2}=B^{*}B\,.$
In a sense made precise by the spectral theorem, self-adjoint operators can usefully be thought of as operators that are "real". An element A of B(H) is called normal if A*A = AA*. Normal operators decompose into the sum of a self-adjoint operator and an imaginary multiple of a self adjoint operator $A={\frac {A+A^{*}}{2}}+i{\frac {A-A^{*}}{2i}}$ that commute with each other. Normal operators can also usefully be thought of in terms of their real and imaginary parts.
An element U of B(H) is called unitary if U is invertible and its inverse is given by U*. This can also be expressed by requiring that U be onto and ⟨Ux, Uy⟩ = ⟨x, y⟩ for all x, y ∈ H. The unitary operators form a group under composition, which is the isometry group of H.
An element of B(H) is compact if it sends bounded sets to relatively compact sets. Equivalently, a bounded operator T is compact if, for any bounded sequence {x$_k$}, the sequence {Tx$_k$} has a convergent subsequence. Many integral operators are compact, and in fact define a special class of operators known as Hilbert–Schmidt operators that are especially important in the study of integral equations. Fredholm operators differ from a compact operator by a multiple of the identity, and are equivalently characterized as operators with a finite dimensional kernel and cokernel. The index of a Fredholm operator T is defined by $\operatorname {index} T=\dim \ker T-\dim \operatorname {coker} T\,.$
The index is homotopy invariant, and plays a deep role in differential geometry via the Atiyah–Singer index theorem.

### Unbounded operators

Unbounded operators are also tractable in Hilbert spaces, and have important applications to quantum mechanics. An unbounded operator T on a Hilbert space H is defined as a linear operator whose domain D(T) is a linear subspace of H. Often the domain D(T) is a dense subspace of H, in which case T is known as a densely defined operator.
The adjoint of a densely defined unbounded operator is defined in essentially the same manner as for bounded operators. Self-adjoint unbounded operators play the role of the observables in the mathematical formulation of quantum mechanics. Examples of self-adjoint unbounded operators on the Hilbert space L$^2$(R) are:
* A suitable extension of the differential operator $(Af)(x)=-i{\frac {\mathrm {d} }{\mathrm {d} x}}f(x)\,,$ where i is the imaginary unit and f is a differentiable function of compact support.
* The multiplication-by-x operator: $(Bf)(x)=xf(x)\,.$
These correspond to the momentum and position observables, respectively. Neither A nor B is defined on all of H, since in the case of A the derivative need not exist, and in the case of B the product function need not be square integrable. In both cases, the set of possible arguments form dense subspaces of L$^2$(R).

## Constructions
### Direct sums

Two Hilbert spaces H$_1$ and H$_2$ can be combined into another Hilbert space, called the (orthogonal) direct sum, and denoted $H_{1}\oplus H_{2}\,,$
consisting of the set of all ordered pairs (x$_1$, x$_2$) where x$_i$ ∈ H$_i$, i = 1, 2, and inner product defined by ${\bigl \langle }(x_{1},x_{2}),(y_{1},y_{2}){\bigr \rangle }_{H_{1}\oplus H_{2}}=\left\langle x_{1},y_{1}\right\rangle _{H_{1}}+\left\langle x_{2},y_{2}\right\rangle _{H_{2}}\,.$
More generally, if H$_i$ is a family of Hilbert spaces indexed by i ∈ I, then the direct sum of the H$_i$, denoted $\bigoplus _{i\in I}H_{i}$ consists of the set of all indexed families $x=(x_{i}\in H_{i}\mid i\in I)\in \prod _{i\in I}H_{i}$ in the Cartesian product of the H$_i$ such that $\sum _{i\in I}\|x_{i}\|^{2}<\infty \,.$
The inner product is defined by $\langle x,y\rangle =\sum _{i\in I}\left\langle x_{i},y_{i}\right\rangle _{H_{i}}\,.$
Each of the H$_i$ is included as a closed subspace in the direct sum of all of the H$_i$. Moreover, the H$_i$ are pairwise orthogonal. Conversely, if there is a system of closed subspaces, V$_i$, i ∈ I, in a Hilbert space H, that are pairwise orthogonal and whose union is dense in H, then H is canonically isomorphic to the direct sum of V$_i$. In this case, H is called the internal direct sum of the V$_i$. A direct sum (internal or external) is also equipped with a family of orthogonal projections E$_i$ onto the ith direct summand H$_i$. These projections are bounded, self-adjoint, idempotent operators that satisfy the orthogonality condition $E_{i}E_{j}=0,\quad i\neq j\,.$
The spectral theorem for compact self-adjoint operators on a Hilbert space H states that H splits into an orthogonal direct sum of the eigenspaces of an operator, and also gives an explicit decomposition of the operator as a sum of projections onto the eigenspaces. The direct sum of Hilbert spaces also appears in quantum mechanics as the Fock space of a system containing a variable number of particles, where each Hilbert space in the direct sum corresponds to an additional degree of freedom for the quantum mechanical system. In representation theory, the Peter–Weyl theorem guarantees that any unitary representation of a compact group on a Hilbert space splits as the direct sum of finite-dimensional representations.

### Tensor products

If x$_1$, y$_1$ ∊ H$_1$ and x$_2$, y$_2$ ∊ H$_2$, then one defines an inner product on the (ordinary) tensor product as follows. On simple tensors, let $\langle x_{1}\otimes x_{2},\,y_{1}\otimes y_{2}\rangle =\langle x_{1},y_{1}\rangle \,\langle x_{2},y_{2}\rangle \,.$
This formula then extends by sesquilinearity to an inner product on H$_1$ ⊗ H$_2$. The Hilbertian tensor product of H$_1$ and H$_2$, sometimes denoted by H$_1$ ${\widehat {\otimes }}$ H$_2$, is the Hilbert space obtained by completing H$_1$ ⊗ H$_2$ for the metric associated to this inner product.
An example is provided by the Hilbert space L$^2$([0, 1]). The Hilbertian tensor product of two copies of L$^2$([0, 1]) is isometrically and linearly isomorphic to the space L$^2$([0, 1]$^2$) of square-integrable functions on the square [0, 1]$^2$. This isomorphism sends a simple tensor f$_1$ ⊗ f$_2$ to the function $(s,t)\mapsto f_{1}(s)\,f_{2}(t)$ on the square.
This example is typical in the following sense. Associated to every simple tensor product x$_1$ ⊗ x$_2$ is the rank one operator from H$^∗$ $_1$ to H$_2$ that maps a given x* ∈ H$^∗$ $_1$ as $x^{*}\mapsto x^{*}(x_{1})x_{2}\,.$
This mapping defined on simple tensors extends to a linear identification between H$_1$ ⊗ H$_2$ and the space of finite rank operators from H$^∗$ $_1$ to H$_2$. This extends to a linear isometry of the Hilbertian tensor product H$_1$ ${\widehat {\otimes }}$ H$_2$ with the Hilbert space HS(H$^∗$ $_1$, H$_2$) of Hilbert–Schmidt operators from H$^∗$ $_1$ to H$_2$.

## Orthonormal bases

The notion of an orthonormal basis from linear algebra generalizes over to the case of Hilbert spaces. In a Hilbert space H, an orthonormal basis is a family {e$_k$}$_{k∈B}$ of elements of H satisfying the conditions:
* Orthogonality: Every two different elements of B are orthogonal: ⟨e$_k$, e$_j$⟩ = 0 for all k, j ∈ B with k ≠ j.
* Normalization: Every element of the family has norm 1: ‖e$_k$‖ = 1 for all k ∈ B.
* Completeness: The linear span of the family e$_k$, k ∈ B, is dense in H.
A system of vectors satisfying the first two conditions basis is called an orthonormal system or an orthonormal set (or an orthonormal sequence if B is countable). Such a system is always linearly independent.
Despite the name, an orthonormal basis is not, in general, a basis in the sense of linear algebra (Hamel basis). More precisely, an orthonormal basis is a Hamel basis if and only if the Hilbert space is a finite-dimensional vector space.
Completeness of an orthonormal system of vectors of a Hilbert space can be equivalently restated as:
    for every v ∈ H, if ⟨v, e$_k$⟩ = 0 for all k ∈ B, then v = 0.
This is related to the fact that the only vector orthogonal to a dense linear subspace is the zero vector, for if S is any orthonormal set and v is orthogonal to S, then v is orthogonal to the closure of the linear span of S, which is the whole space.
Examples of orthonormal bases include:
* the set {(1, 0, 0), (0, 1, 0), (0, 0, 1)} forms an orthonormal basis of R$^3$ with the dot product;
* the sequence { f$_n$ | n ∈ Z} with f$_n$(x) = exp(2πinx) forms an orthonormal basis of the complex space L$^2$([0, 1]);
In the infinite-dimensional case, an orthonormal basis will not be a basis in the sense of linear algebra; to distinguish the two, the latter basis is also called a Hamel basis. That the span of the basis vectors is dense implies that every vector in the space can be written as the sum of an infinite series, and the orthogonality implies that this decomposition is unique.

### Sequence spaces

The space $\ell _{2}$ of square-summable sequences of complex numbers is the set of infinite sequences $(c_{1},c_{2},c_{3},\dots )$ of real or complex numbers such that $\left|c_{1}\right|^{2}+\left|c_{2}\right|^{2}+\left|c_{3}\right|^{2}+\cdots <\infty \,.$
This space has an orthonormal basis: ${\begin{aligned}e_{1}&=(1,0,0,\dots )\\e_{2}&=(0,1,0,\dots )\\&\ \ \vdots \end{aligned}}$
This space is the infinite-dimensional generalization of the $\ell _{2}^{n}$ space of finite-dimensional vectors. It is usually the first example used to show that in infinite-dimensional spaces, a set that is closed and bounded is not necessarily (sequentially) compact (as is the case in all finite dimensional spaces). Indeed, the set of orthonormal vectors above shows this: It is an infinite sequence of vectors in the unit ball (i.e., the ball of points with norm less than or equal one). This set is clearly bounded and closed; yet, no subsequence of these vectors converges to anything and consequently the unit ball in $\ell _{2}$ is not compact. Intuitively, this is because "there is always another coordinate direction" into which the next elements of the sequence can evade.
One can generalize the space $\ell _{2}$ in many ways. For example, if B is any set, then one can form a Hilbert space of sequences with index set B, defined by $\ell ^{2}(B)={\biggl \x:B\xrightarrow {x} \mathbb {C} \mathrel {\bigg |} \sum _{b\in B}\left|x(b)\right|^{2}<\infty {\biggr \}}\,.$
The summation over B is here defined by $\sum _{b\in B}\left|x(b)\right|^{2}=\sup \sum _{n=1}^{N}\left|x(b_{n})\right|^{2}$ the supremum being taken over all finite subsets of B. It follows that, for this sum to be finite, every element of l$^2$(B) has only countably many nonzero terms. This space becomes a Hilbert space with the inner product $\langle x,y\rangle =\sum _{b\in B}x(b){\overline {y(b)}}$
for all x, y ∈ l$^2$(B). Here the sum also has only countably many nonzero terms, and is unconditionally convergent by the Cauchy–Schwarz inequality.
An orthonormal basis of l$^2$(B) is indexed by the set B, given by
    $e_{b}(b')={\begin{cases}1&{\text{if }}b=b'\\0&{\text{otherwise.}}\end{cases}}$

### Bessel's inequality and Parseval's formula

Let f$_1$, ..., f$_n$ be a finite orthonormal system in H. For an arbitrary vector x ∈ H, let $y=\sum _{j=1}^{n}\langle x,f_{j}\rangle \,f_{j}\,.$
Then ⟨x, f$_k$⟩ = ⟨y, f$_k$⟩ for every k = 1, ..., n. It follows that x − y is orthogonal to each f$_k$, hence x − y is orthogonal to y. Using the Pythagorean identity twice, it follows that $\|x\|^{2}=\|x-y\|^{2}+\|y\|^{2}\geq \|y\|^{2}=\sum _{j=1}^{n}{\bigl |}\langle x,f_{j}\rangle {\bigr |}^{2}\,.$
Let {f$_i$}, i ∈ I, be an arbitrary orthonormal system in H. Applying the preceding inequality to every finite subset J of I gives Bessel's inequality: $\sum _{i\in I}{\bigl |}\langle x,f_{i}\rangle {\bigr |}^{2}\leq \|x\|^{2},\quad x\in H$ (according to the definition of the sum of an arbitrary family of non-negative real numbers).
Geometrically, Bessel's inequality implies that the orthogonal projection of x onto the linear subspace spanned by the f$_i$ has norm that does not exceed that of x. In two dimensions, this is the assertion that the length of the leg of a right triangle may not exceed the length of the hypotenuse.
Bessel's inequality is a stepping stone to the stronger result called Parseval's identity, which governs the case when Bessel's inequality is actually an equality. By definition, if {e$_k$}$_{k∈B}$ is an orthonormal basis of H, then every element x of H may be written as $x=\sum _{k\in B}\left\langle x,e_{k}\right\rangle \,e_{k}\,.$
Even if B is uncountable, Bessel's inequality guarantees that the expression is well-defined and consists only of countably many nonzero terms. This sum is called the Fourier expansion of x, and the individual coefficients ⟨x, e$_k$⟩ are the Fourier coefficients of x. Parseval's identity then asserts that $\|x\|^{2}=\sum _{k\in B}|\langle x,e_{k}\rangle |^{2}\,.$
Conversely, if {e$_k$} is an orthonormal set such that Parseval's identity holds for every x, then {e$_k$} is an orthonormal basis.

### Hilbert dimension

As a consequence of Zorn's lemma, every Hilbert space admits an orthonormal basis; furthermore, any two orthonormal bases of the same space have the same cardinality, called the Hilbert dimension of the space. For instance, since l$^2$(B) has an orthonormal basis indexed by B, its Hilbert dimension is the cardinality of B (which may be a finite integer, or a countable or uncountable cardinal number).
The Hilbert dimension is not greater than the Hamel dimension (the usual dimension of a vector space).
As a consequence of Parseval's identity, if {e$_k$}$_{k∈B}$ is an orthonormal basis of H, then the map Φ : H → l$^2$(B) defined by Φ(x) = ⟨x, e$_k$⟩$_{k∈B}$ is an isometric isomorphism of Hilbert spaces: it is a bijective linear mapping such that ${\bigl \langle }\Phi (x),\Phi (y){\bigr \rangle }_{l^{2}(B)}=\left\langle x,y\right\rangle _{H}$ for all x, y ∈ H. The cardinal number of B is the Hilbert dimension of H. Thus every Hilbert space is isometrically isomorphic to a sequence space l$^2$(B) for some set B.

### Separable spaces

By definition, a Hilbert space is separable provided it contains a dense countable subset. Along with Zorn's lemma, this means a Hilbert space is separable if and only if it admits a countable orthonormal basis. All infinite-dimensional separable Hilbert spaces are therefore isometrically isomorphic to the square-summable sequence space $\ell ^{2}.$
In the past, Hilbert spaces were often required to be separable as part of the definition.

#### In quantum field theory

Most spaces used in physics are separable, and since these are all isomorphic to each other, one often refers to any infinite-dimensional separable Hilbert space as "the Hilbert space" or just "Hilbert space". Even in quantum field theory, most of the Hilbert spaces are in fact separable, as stipulated by the Wightman axioms. However, it is sometimes argued that non-separable Hilbert spaces are also important in quantum field theory, roughly because the systems in the theory possess an infinite number of degrees of freedom and any infinite Hilbert tensor product (of spaces of dimension greater than one) is non-separable. For instance, a bosonic field can be naturally thought of as an element of a tensor product whose factors represent harmonic oscillators at each point of space. From this perspective, the natural state space of a boson might seem to be a non-separable space. However, it is only a small separable subspace of the full tensor product that can contain physically meaningful fields (on which the observables can be defined). Another non-separable Hilbert space models the state of an infinite collection of particles in an unbounded region of space. An orthonormal basis of the space is indexed by the density of the particles, a continuous parameter, and since the set of possible densities is uncountable, the basis is not countable.

## Orthogonal complements and projections

If S is a subset of a Hilbert space H, the set of vectors orthogonal to S is defined by $S^{\perp }=\left\{x\in H\mid \langle x,s\rangle =0\ {\text{ for all }}s\in S\right\}\,.$
The set S$^⊥$ is a closed subspace of H (can be proved easily using the linearity and continuity of the inner product) and so forms itself a Hilbert space. If V is a closed subspace of H, then V$^⊥$ is called the orthogonal complement of V. In fact, every x ∈ H can then be written uniquely as x = v + w, with v ∈ V and w ∈ V$^⊥$. Therefore, H is the internal Hilbert direct sum of V and V$^⊥$.
The linear operator P$_V$ : H → H that maps x to v is called the orthogonal projection onto V. There is a natural one-to-one correspondence between the set of all closed subspaces of H and the set of all bounded self-adjoint operators P such that P$^2$ = P. Specifically,
Theorem — The orthogonal projection P$_V$ is a self-adjoint linear operator on H of norm ≤ 1 with the property P$^2$ $_V$ = P$_V$. Moreover, any self-adjoint linear operator E such that E$^2$ = E is of the form P$_V$, where V is the range of E. For every x in H, P$_V$(x) is the unique element v of V that minimizes the distance ‖x − v‖.
This provides the geometrical interpretation of P$_V$(x): it is the best approximation to x by elements of V.
Projections P$_U$ and P$_V$ are called mutually orthogonal if P$_U$P$_V$ = 0. This is equivalent to U and V being orthogonal as subspaces of H. The sum of the two projections P$_U$ and P$_V$ is a projection only if U and V are orthogonal to each other, and in that case P$_U$ + P$_V$ = P$_{U+V}$. The composite P$_U$P$_V$ is generally not a projection; in fact, the composite is a projection if and only if the two projections commute, and in that case P$_U$P$_V$ = P$_{U∩V}$.
By restricting the codomain to the Hilbert space V, the orthogonal projection P$_V$ gives rise to a projection mapping π : H → V; it is the adjoint of the inclusion mapping $i:V\to H\,,$ meaning that $\left\langle ix,y\right\rangle _{H}=\left\langle x,\pi y\right\rangle _{V}$ for all x ∈ V and y ∈ H.
The operator norm of the orthogonal projection P$_V$ onto a nonzero closed subspace V is equal to 1: $\|P_{V}\|=\sup _{x\in H,x\neq 0}{\frac {\|P_{V}x\|}{\|x\|}}=1\,.$
Every closed subspace V of a Hilbert space is therefore the image of an operator P of norm one such that P$^2$ = P. The property of possessing appropriate projection operators characterizes Hilbert spaces:
* A Banach space of dimension higher than 2 is (isometrically) a Hilbert space if and only if, for every closed subspace V, there is an operator P$_V$ of norm one whose image is V such that P$^2$ $_V$ = P$_V$.
While this result characterizes the metric structure of a Hilbert space, the structure of a Hilbert space as a topological vector space can itself be characterized in terms of the presence of complementary subspaces:
* A Banach space X is topologically and linearly isomorphic to a Hilbert space if and only if, to every closed subspace V, there is a closed subspace W such that X is equal to the internal direct sum V ⊕ W.
The orthogonal complement satisfies some more elementary results. It is a monotone function in the sense that if U ⊂ V, then V$^⊥$ ⊆ U$^⊥$ with equality holding if and only if V is contained in the closure of U. This result is a special case of the Hahn–Banach theorem. The closure of a subspace can be completely characterized in terms of the orthogonal complement: if V is a subspace of H, then the closure of V is equal to V$^{⊥⊥}$. The orthogonal complement is thus a Galois connection on the partial order of subspaces of a Hilbert space. In general, the orthogonal complement of a sum of subspaces is the intersection of the orthogonal complements: ${\biggl\sum _{i}V_{i}{\biggr )}^{\perp }=\bigcap _{i}V_{i}^{\perp }\,.$
If the V$_i$ are in addition closed, then ${\overline {\sum _{i}V_{i}^{\perp }{\vphantom {\Big |}}}}={\biggl\bigcap _{i}V_{i}{\biggr )}^{\perp }\,.$

## Spectral theory

There is a well-developed spectral theory for self-adjoint operators in a Hilbert space, that is roughly analogous to the study of symmetric matrices over the reals or self-adjoint matrices over the complex numbers. In the same sense, one can obtain a "diagonalization" of a self-adjoint operator as a suitable sum (actually an integral) of orthogonal projection operators.
The spectrum of an operator T, denoted σ(T), is the set of complex numbers λ such that T − λ lacks a continuous inverse. If T is bounded, then the spectrum is always a compact set in the complex plane, and lies inside the disc |z| ≤ ‖T‖. If T is self-adjoint, then the spectrum is real. In fact, it is contained in the interval [m, M] where $m=\inf _{\|x\|=1}\langle Tx,x\rangle \,,\quad M=\sup _{\|x\|=1}\langle Tx,x\rangle \,.$
Moreover, m and M are both actually contained within the spectrum.
The eigenspaces of an operator T are given by $H_{\lambda }=\ker(T-\lambda )\,.$
Unlike with finite matrices, not every element of the spectrum of T must be an eigenvalue: the linear operator T − λ may only lack an inverse because it is not surjective. Elements of the spectrum of an operator in the general sense are known as spectral values. Since spectral values need not be eigenvalues, the spectral decomposition is often more subtle than in finite dimensions.
However, the spectral theorem of a self-adjoint operator T takes a particularly simple form if, in addition, T is assumed to be a compact operator. The spectral theorem for compact self-adjoint operators states:
* A compact self-adjoint operator T has only countably (or finitely) many spectral values. The spectrum of T has no limit point in the complex plane except possibly zero. The eigenspaces of T decompose H into an orthogonal direct sum: $H=\bigoplus _{\lambda \in \sigma (T)}H_{\lambda }\,.$ Moreover, if E$_λ$ denotes the orthogonal projection onto the eigenspace H$_λ$, then $T=\sum _{\lambda \in \sigma (T)}\lambda E_{\lambda }\,,$ where the sum converges with respect to the norm on B(H).
This theorem plays a fundamental role in the theory of integral equations, as many integral operators are compact, in particular those that arise from Hilbert–Schmidt operators.
The general spectral theorem for self-adjoint operators involves a kind of operator-valued Riemann–Stieltjes integral, rather than an infinite summation. The spectral family associated to T associates to each real number λ an operator E$_λ$, which is the projection onto the nullspace of the operator (T − λ)$^+$, where the positive part of a self-adjoint operator is defined by $A^{+}={\tfrac {1}{2}}{\Bigl{\sqrt {A^{2}}}+A{\Bigr )}\,.$
The operators E$_λ$ are monotone increasing relative to the partial order defined on self-adjoint operators; the eigenvalues correspond precisely to the jump discontinuities. One has the spectral theorem, which asserts $T=\int _{\mathbb {R} }\lambda \,\mathrm {d} E_{\lambda }\,.$
The integral is understood as a Riemann–Stieltjes integral, convergent with respect to the norm on B(H). In particular, one has the ordinary scalar-valued integral representation $\langle Tx,y\rangle =\int _{\mathbb {R} }\lambda \,\mathrm {d} \langle E_{\lambda }x,y\rangle \,.$
A somewhat similar spectral decomposition holds for normal operators, although because the spectrum may now contain non-real complex numbers, the operator-valued Stieltjes measure dE$_λ$ must instead be replaced by a resolution of the identity.
A major application of spectral methods is the spectral mapping theorem, which allows one to apply to a self-adjoint operator T any continuous complex function f defined on the spectrum of T by forming the integral $f(T)=\int _{\sigma (T)}f(\lambda )\,\mathrm {d} E_{\lambda }\,.$
The resulting continuous functional calculus has applications in particular to pseudodifferential operators.
The spectral theory of unbounded self-adjoint operators is only marginally more difficult than for bounded operators. The spectrum of an unbounded operator is defined in precisely the same way as for bounded operators: λ is a spectral value if the resolvent operator $R_{\lambda }=(T-\lambda )^{-1}$
fails to be a well-defined continuous operator. The self-adjointness of T still guarantees that the spectrum is real. Thus the essential idea of working with unbounded operators is to look instead at the resolvent R$_λ$ where λ is nonreal. This is a bounded normal operator, which admits a spectral representation that can then be transferred to a spectral representation of T itself. A similar strategy is used, for instance, to study the spectrum of the Laplace operator: rather than address the operator directly, one instead looks as an associated resolvent such as a Riesz potential or Bessel potential.
A precise version of the spectral theorem in this case is:
Theorem — Given a densely defined self-adjoint operator T on a Hilbert space H, there corresponds a unique resolution of the identity E on the Borel sets of R, such that $\langle Tx,y\rangle =\int _{\mathbb {R} }\lambda \,\mathrm {d} E_{x,y}(\lambda )$ for all x ∈ D(T) and y ∈ H. The spectral measure E is concentrated on the spectrum of T.
There is also a version of the spectral theorem that applies to unbounded normal operators.

## In popular culture

In Gravity's Rainbow (1973), a novel by Thomas Pynchon, one of the characters is called "Sammy Hilbert-Spaess", a pun on "Hilbert Space". The novel refers also to Gödel's incompleteness theorems.