Wave function

In quantum physics, a wave function (or wavefunction) is a mathematical description of the quantum state of an isolated quantum system. The most common symbols for a wave function are the Greek letters ψ and Ψ (lower-case and capital psi, respectively). Wave functions are complex-valued. For example, a wave function might assign a complex number to each point in a region of space. The Born rule^[1][2][3] provides the means to turn these complex probability amplitudes into actual probabilities. In one common form, it says that the squared modulus of a wave function that depends upon position is the probability density of measuring a particle as being at a given place. The integral of a wavefunction's squared modulus over all the system's degrees of freedom must be equal to 1, a condition called normalization. Since the wave function is complex-valued, only its relative phase and relative magnitude can be measured; its value does not, in isolation, tell anything about the magnitudes or directions of measurable observables. One has to apply quantum operators, whose eigenvalues correspond to sets of possible results of measurements, to the wave function ψ and calculate the statistical distributions for measurable quantities.

Wave functions can be functions of variables other than position, such as momentum. The information represented by a wave function that is dependent upon position can be converted into a wave function dependent upon momentum and vice versa, by means of a Fourier transform. Some particles, like electrons and photons, have nonzero spin, and the wave function for such particles includes spin as an intrinsic, discrete degree of freedom; other discrete variables can also be included, such as isospin. When a system has internal degrees of freedom, the wave function at each point in the continuous degrees of freedom (e.g., a point in space) assigns a complex number for each possible value of the discrete degrees of freedom (e.g., z-component of spin). These values are often displayed in a column matrix (e.g., a 2 × 1 column vector for a non-relativistic electron with spin 1⁄2).

According to the superposition principle of quantum mechanics, wave functions can be added together and multiplied by complex numbers to form new wave functions and form a Hilbert space. The inner product between two wave functions is a measure of the overlap between the corresponding physical states and is used in the foundational probabilistic interpretation of quantum mechanics, the Born rule, relating transition probabilities to inner products. The Schrödinger equation determines how wave functions evolve over time, and a wave function behaves qualitatively like other waves, such as water waves or waves on a string, because the Schrödinger equation is mathematically a type of wave equation. This explains the name "wave function", and gives rise to wave–particle duality. However, the wave function in quantum mechanics describes a kind of physical phenomenon, as of 2023 still open to different interpretations, which fundamentally differs from that of classic mechanical waves.^{[4][5][6][7][8][9][10]}

Part of a series of articles about
Quantum mechanics
${\displaystyle i\hbar {\frac {d}{dt}}|\Psi \rangle ={\hat {H}}|\Psi \rangle }$ Schrödinger equation
Introduction Glossary History
Background Classical mechanics Old quantum theory Bra–ket notation Hamiltonian Interference
Fundamentals Complementarity Decoherence Entanglement Energy level Measurement Nonlocality Quantum number State Superposition Symmetry Tunnelling Uncertainty Wave function Collapse
Experiments Bell's inequality CHSH inequality Davisson–Germer Double-slit Elitzur–Vaidman Franck–Hertz Leggett inequality Leggett–Garg inequality Mach–Zehnder Popper Quantum eraser Delayed-choice Schrödinger's cat Stern–Gerlach Wheeler's delayed-choice
Formulations Overview Heisenberg Interaction Matrix Phase-space Schrödinger Sum-over-histories (path integral)
Equations Dirac Klein–Gordon Pauli Rydberg Schrödinger
Interpretations Bayesian Consistent histories Copenhagen de Broglie–Bohm Ensemble Hidden-variable Local Superdeterminism Many-worlds Objective-collapse Quantum logic Relational Transactional Von Neumann–Wigner
Advanced topics Relativistic quantum mechanics Quantum field theory Quantum information science Quantum computing Quantum chaos EPR paradox Density matrix Scattering theory Quantum statistical mechanics Quantum machine learning
Scientists Aharonov Bell Bethe Blackett Bloch Bohm Bohr Born Bose de Broglie Compton Dirac Davisson Debye Ehrenfest Einstein Everett Fock Fermi Feynman Glauber Gutzwiller Heisenberg Hilbert Jordan Kramers Lamb Landau Laue Moseley Millikan Onnes Pauli Planck Rabi Raman Rydberg Schrödinger Simmons Sommerfeld von Neumann Weyl Wien Wigner Zeeman Zeilinger
v t e

In 1900, Max Planck postulated the proportionality between the frequency ${\displaystyle f}$ of a photon and its energy ${\displaystyle E}$, ${\displaystyle E=hf}$,^[11][12] and in 1916 the corresponding relation between a photon's momentum ${\displaystyle p}$ and wavelength ${\displaystyle \lambda }$, ${\displaystyle \lambda ={\frac {h}{p}}}$,^[13] where ${\displaystyle h}$ is the Planck constant. In 1923, De Broglie was the first to suggest that the relation ${\displaystyle \lambda ={\frac {h}{p}}}$, now called the De Broglie relation, holds for massive particles, the chief clue being Lorentz invariance,^[14] and this can be viewed as the starting point for the modern development of quantum mechanics. The equations represent wave–particle duality for both massless and massive particles.

In the 1920s and 1930s, quantum mechanics was developed using calculus and linear algebra. Those who used the techniques of calculus included Louis de Broglie, Erwin Schrödinger, and others, developing "wave mechanics". Those who applied the methods of linear algebra included Werner Heisenberg, Max Born, and others, developing "matrix mechanics". Schrödinger subsequently showed that the two approaches were equivalent.^[15]

In 1926, Schrödinger published the famous wave equation now named after him, the Schrödinger equation. This equation was based on classical conservation of energy using quantum operators and the de Broglie relations and the solutions of the equation are the wave functions for the quantum system.^[16] However, no one was clear on how to interpret it.^[17] At first, Schrödinger and others thought that wave functions represent particles that are spread out with most of the particle being where the wave function is large.^[18] This was shown to be incompatible with the elastic scattering of a wave packet (representing a particle) off a target; it spreads out in all directions.^[1] While a scattered particle may scatter in any direction, it does not break up and take off in all directions. In 1926, Born provided the perspective of probability amplitude.^[1][2][19] This relates calculations of quantum mechanics directly to probabilistic experimental observations. It is accepted as part of the Copenhagen interpretation of quantum mechanics. There are many other interpretations of quantum mechanics. In 1927, Hartree and Fock made the first step in an attempt to solve the N-body wave function, and developed the self-consistency cycle: an iterative algorithm to approximate the solution. Now it is also known as the Hartree–Fock method.^[20] The Slater determinant and permanent (of a matrix) was part of the method, provided by John C. Slater.

Schrödinger did encounter an equation for the wave function that satisfied relativistic energy conservation before he published the non-relativistic one, but discarded it as it predicted negative probabilities and negative energies. In 1927, Klein, Gordon and Fock also found it, but incorporated the electromagnetic interaction and proved that it was Lorentz invariant. De Broglie also arrived at the same equation in 1928. This relativistic wave equation is now most commonly known as the Klein–Gordon equation.^[21]

In 1927, Pauli phenomenologically found a non-relativistic equation to describe spin-1/2 particles in electromagnetic fields, now called the Pauli equation.^[22] Pauli found the wave function was not described by a single complex function of space and time, but needed two complex numbers, which respectively correspond to the spin +1/2 and −1/2 states of the fermion. Soon after in 1928, Dirac found an equation from the first successful unification of special relativity and quantum mechanics applied to the electron, now called the Dirac equation. In this, the wave function is a spinor represented by four complex-valued components:^[20] two for the electron and two for the electron's antiparticle, the positron. In the non-relativistic limit, the Dirac wave function resembles the Pauli wave function for the electron. Later, other relativistic wave equations were found.

All these wave equations are of enduring importance. The Schrödinger equation and the Pauli equation are under many circumstances excellent approximations of the relativistic variants. They are considerably easier to solve in practical problems than the relativistic counterparts.

The Klein–Gordon equation and the Dirac equation, while being relativistic, do not represent full reconciliation of quantum mechanics and special relativity. The branch of quantum mechanics where these equations are studied the same way as the Schrödinger equation, often called relativistic quantum mechanics, while very successful, has its limitations (see e.g. Lamb shift) and conceptual problems (see e.g. Dirac sea).

Relativity makes it inevitable that the number of particles in a system is not constant. For full reconciliation, quantum field theory is needed.^[23] In this theory, the wave equations and the wave functions have their place, but in a somewhat different guise. The main objects of interest are not the wave functions, but rather operators, so called field operators (or just fields where "operator" is understood) on the Hilbert space of states (to be described next section). It turns out that the original relativistic wave equations and their solutions are still needed to build the Hilbert space. Moreover, the free fields operators, i.e. when interactions are assumed not to exist, turn out to (formally) satisfy the same equation as do the fields (wave functions) in many cases.

Thus the Klein–Gordon equation (spin 0) and the Dirac equation (spin 1⁄2) in this guise remain in the theory. Higher spin analogues include the Proca equation (spin 1), Rarita–Schwinger equation (spin 3⁄2), and, more generally, the Bargmann–Wigner equations. For massless free fields two examples are the free field Maxwell equation (spin 1) and the free field Einstein equation (spin 2) for the field operators.^[24] All of them are essentially a direct consequence of the requirement of Lorentz invariance. Their solutions must transform under Lorentz transformation in a prescribed way, i.e. under a particular representation of the Lorentz group and that together with few other reasonable demands, e.g. the cluster decomposition property,^[25] with implications for causality is enough to fix the equations.

This applies to free field equations; interactions are not included. If a Lagrangian density (including interactions) is available, then the Lagrangian formalism will yield an equation of motion at the classical level. This equation may be very complex and not amenable to solution. Any solution would refer to a fixed number of particles and would not account for the term "interaction" as referred to in these theories, which involves the creation and annihilation of particles and not external potentials as in ordinary "first quantized" quantum theory.

In string theory, the situation remains analogous. For instance, a wave function in momentum space has the role of Fourier expansion coefficient in a general state of a particle (string) with momentum that is not sharply defined.^[26]

This section does not cite any sources. Please help improve this section by adding citations to reliable sources. Unsourced material may be challenged and removed. (May 2021) (Learn how and when to remove this message)

Standing waves for a particle in a box, examples of stationary states.

Travelling waves of a free particle.

The real parts of position wave function Ψ(x) and momentum wave function Φ(p), and corresponding probability densities |Ψ(x)|² and |Φ(p)|², for one spin-0 particle in one x or p dimension. The colour opacity of the particles corresponds to the probability density (not the wave function) of finding the particle at position x or momentum p.

For now, consider the simple case of a non-relativistic single particle, without spin, in one spatial dimension. More general cases are discussed below.

According to the postulates of quantum mechanics, the state of a physical system, at fixed time ${\displaystyle t}$, is given by the wave function belonging to a separable complex Hilbert space.^[27][28] As such, the inner product of two wave functions Ψ₁ and Ψ₂ can be defined as the complex number (at time t)^{[nb 1]}

${\displaystyle (\Psi _{1},\Psi _{2})=\int _{-\infty }^{\infty }\,\Psi _{1}^{*}(x,t)\Psi _{2}(x,t)\,dx<\infty }$.

More details are given below. However, the inner product of a wave function Ψ with itself,

${\displaystyle (\Psi ,\Psi )=\|\Psi \|^{2}}$,

is always a positive real number. The number ‖Ψ‖ (not ‖Ψ‖²) is called the norm of the wave function Ψ. The separable Hilbert space being considered is infinite-dimensional,^{[nb 2]} which means there is no finite set of square integrable functions which can be added together in various combinations to create every possible square integrable function.

The state of such a particle is completely described by its wave function, $${\displaystyle \Psi (x,t)\,,}$$ where x is position and t is time. This is a complex-valued function of two real variables x and t.

For one spinless particle in one dimension, if the wave function is interpreted as a probability amplitude; the square modulus of the wave function, the positive real number $${\displaystyle \left|\Psi (x,t)\right|^{2}=\Psi ^{*}(x,t)\Psi (x,t)=\rho (x),}$$ is interpreted as the probability density for a measurement of the particle's position at a given time t. The asterisk indicates the complex conjugate. If the particle's position is measured, its location cannot be determined from the wave function, but is described by a probability distribution.

The probability that its position x will be in the interval a ≤ x ≤ b is the integral of the density over this interval: $${\displaystyle P_{a\leq x\leq b}(t)=\int _{a}^{b}\,|\Psi (x,t)|^{2}dx}$$ where t is the time at which the particle was measured. This leads to the normalization condition: $${\displaystyle \int _{-\infty }^{\infty }\,|\Psi (x,t)|^{2}dx=1\,,}$$ because if the particle is measured, there is 100% probability that it will be somewhere.

For a given system, the set of all possible normalizable wave functions (at any given time) forms an abstract mathematical vector space, meaning that it is possible to add together different wave functions, and multiply wave functions by complex numbers. Technically, wave functions form a ray in a projective Hilbert space rather than an ordinary vector space.

At a particular instant of time, all values of the wave function Ψ(x, t) are components of a vector. There are uncountably infinitely many of them and integration is used in place of summation. In Bra–ket notation, this vector is written $${\displaystyle |\Psi (t)\rangle =\int \Psi (x,t)|x\rangle dx}$$ and is referred to as a "quantum state vector", or simply "quantum state". There are several advantages to understanding wave functions as representing elements of an abstract vector space:

• All the powerful tools of linear algebra can be used to manipulate and understand wave functions. For example:
  • Linear algebra explains how a vector space can be given a basis, and then any vector in the vector space can be expressed in this basis. This explains the relationship between a wave function in position space and a wave function in momentum space and suggests that there are other possibilities too.
  • Bra–ket notation can be used to manipulate wave functions.
• The idea that quantum states are vectors in an abstract vector space is completely general in all aspects of quantum mechanics and quantum field theory, whereas the idea that quantum states are complex-valued "wave" functions of space is only true in certain situations.

The time parameter is often suppressed, and will be in the following. The x coordinate is a continuous index. The |x⟩ are called improper vectors which, unlike proper vectors that are normalizable to unity, can only be normalized to a Dirac delta function.^{[nb 3][nb 4][29]} $${\displaystyle \langle x'|x\rangle =\delta (x'-x)}$$ thus $${\displaystyle \langle x'|\Psi \rangle =\int \Psi (x)\langle x'|x\rangle dx=\Psi (x')}$$ and $${\displaystyle |\Psi \rangle =\int |x\rangle \langle x|\Psi \rangle dx=\left(\int |x\rangle \langle x|dx\right)|\Psi \rangle }$$ which illuminates the identity operator $${\displaystyle I=\int |x\rangle \langle x|dx\,.}$$which is analogous to completeness relation of orthonormal basis in N-dimensional Hilbert space.

Finding the identity operator in a basis allows the abstract state to be expressed explicitly in a basis, and more (the inner product between two state vectors, and other operators for observables, can be expressed in the basis).

The particle also has a wave function in momentum space: $${\displaystyle \Phi (p,t)}$$ where p is the momentum in one dimension, which can be any value from −∞ to +∞, and t is time.

Analogous to the position case, the inner product of two wave functions Φ₁(p, t) and Φ₂(p, t) can be defined as: $${\displaystyle (\Phi _{1},\Phi _{2})=\int _{-\infty }^{\infty }\,\Phi _{1}^{*}(p,t)\Phi _{2}(p,t)dp\,.}$$

One particular solution to the time-independent Schrödinger equation is $${\displaystyle \Psi _{p}(x)=e^{ipx/\hbar },}$$ a plane wave, which can be used in the description of a particle with momentum exactly p, since it is an eigenfunction of the momentum operator. These functions are not normalizable to unity (they are not square-integrable), so they are not really elements of physical Hilbert space. The set $${\displaystyle \{\Psi _{p}(x,t),-\infty \leq p\leq \infty \}}$$ forms what is called the momentum basis. This "basis" is not a basis in the usual mathematical sense. For one thing, since the functions are not normalizable, they are instead normalized to a delta function,^{[nb 4]} $${\displaystyle (\Psi _{p},\Psi _{p'})=\delta (p-p').}$$

For another thing, though they are linearly independent, there are too many of them (they form an uncountable set) for a basis for physical Hilbert space. They can still be used to express all functions in it using Fourier transforms as described next.

The x and p representations are $${\displaystyle {\begin{aligned}|\Psi \rangle =I|\Psi \rangle &=\int |x\rangle \langle x|\Psi \rangle dx=\int \Psi (x)|x\rangle dx,\\|\Psi \rangle =I|\Psi \rangle &=\int |p\rangle \langle p|\Psi \rangle dp=\int \Phi (p)|p\rangle dp.\end{aligned}}}$$

Now take the projection of the state Ψ onto eigenfunctions of momentum using the last expression in the two equations, $${\displaystyle \int \Psi (x)\langle p|x\rangle dx=\int \Phi (p')\langle p|p'\rangle dp'=\int \Phi (p')\delta (p-p')dp'=\Phi (p).}$$

Then utilizing the known expression for suitably normalized eigenstates of momentum in the position representation solutions of the free Schrödinger equation $${\displaystyle \langle x|p\rangle =p(x)={\frac {1}{\sqrt {2\pi \hbar }}}e^{{\frac {i}{\hbar }}px}\Rightarrow \langle p|x\rangle ={\frac {1}{\sqrt {2\pi \hbar }}}e^{-{\frac {i}{\hbar }}px},}$$ one obtains $${\displaystyle \Phi (p)={\frac {1}{\sqrt {2\pi \hbar }}}\int \Psi (x)e^{-{\frac {i}{\hbar }}px}dx\,.}$$

Likewise, using eigenfunctions of position, $${\displaystyle \Psi (x)={\frac {1}{\sqrt {2\pi \hbar }}}\int \Phi (p)e^{{\frac {i}{\hbar }}px}dp\,.}$$

The position-space and momentum-space wave functions are thus found to be Fourier transforms of each other.^[30] They are two representations of the same state; containing the same information, and either one is sufficient to calculate any property of the particle.

In practice, the position-space wave function is used much more often than the momentum-space wave function. The potential entering the relevant equation (Schrödinger, Dirac, etc.) determines in which basis the description is easiest. For the harmonic oscillator, x and p enter symmetrically, so there it does not matter which description one uses. The same equation (modulo constants) results. From this, with a little bit of afterthought, it follows that solutions to the wave equation of the harmonic oscillator are eigenfunctions of the Fourier transform in L².^{[nb 5]}

Following are the general forms of the wave function for systems in higher dimensions and more particles, as well as including other degrees of freedom than position coordinates or momentum components.

While Hilbert spaces originally refer to infinite dimensional complete inner product spaces they, by definition, include finite dimensional complete inner product spaces as well.^[31] In physics, they are often referred to as finite dimensional Hilbert spaces.^[32] For every finite dimensional Hilbert space there exist orthonormal basis kets that span the entire Hilbert space.

If the N-dimensional set ${\textstyle \{|\phi _{i}\rangle \}}$ is orthonormal, then the projection operator for the space spanned by these states is given by:

$${\displaystyle P=\sum _{i}|\phi _{i}\rangle \langle \phi _{i}|=I}$$where the projection is equivalent to identity operator since ${\textstyle \{|\phi _{i}\rangle \}}$ spans the entire Hilbert space, thus leaving any vector from Hilbert space unchanged. This is also known as completeness relation of finite dimensional Hilbert space.

The wavefunction is instead given by:

$${\displaystyle |\psi \rangle =I|\psi \rangle =\sum _{i}|\phi _{i}\rangle \langle \phi _{i}|\psi \rangle }$$where ${\textstyle \{\langle \phi _{i}|\psi \rangle \}}$, is a set of complex numbers which can be used to construct a wavefunction using the above formula.

If the set ${\textstyle \{|\phi _{i}\rangle \}}$ are eigenkets of a non-degenerate observable with eigenvalues ${\textstyle \lambda _{i}}$, by the postulates of quantum mechanics, the probability of measuring the observable to be ${\textstyle \lambda _{i}}$ is given according to Born rule as:^[33]

$${\displaystyle P_{\psi }(\lambda _{i})=|\langle \phi _{i}|\psi \rangle |^{2}}$$

For non-degenerate ${\textstyle \{|\phi _{i}\rangle \}}$ of some observable, if eigenvalues ${\textstyle \lambda }$ have subset of eigenvectors labelled as ${\textstyle \{|\lambda ^{(j)}\rangle \}}$, by the postulates of quantum mechanics, the probability of measuring the observable to be ${\textstyle \lambda }$ is given by:

$${\displaystyle P_{\psi }(\lambda )=\sum _{j}|\langle \lambda ^{(j)}|\psi \rangle |^{2}=|{\widehat {P}}_{\lambda }|\psi \rangle |^{2}}$$where ${\textstyle {\widehat {P}}_{\lambda }=\sum _{j}|\lambda ^{(j)}\rangle \langle \lambda ^{(j)}|}$ is a projection operator of states to subspace spanned by ${\textstyle \{|\lambda ^{(j)}\rangle \}}$. The equality follows due to orthogonal nature of ${\textstyle \{|\phi _{i}\rangle \}}$.

Hence, ${\textstyle \{\langle \phi _{i}|\psi \rangle \}}$ which specify state of the quantum mechanical system, have magnitudes whose square gives the probability of measuring the respective ${\textstyle |\phi _{i}\rangle }$ state.

While the relative phase has observable effects in experiments, the global phase of the system is experimentally indistinguishable. For example in a particle in superposition of two states, the global phase of the particle cannot be distinguished by finding expectation value of observable or probabilities of observing different states but relative phases can affect the expectation values of observables.

While the overall phase of the system is considered to be arbitrary, the relative phase for each state ${\textstyle |\phi _{i}\rangle }$ of a prepared state in superposition can be determined based on physical meaning of the prepared state and its symmetry. For example, the construction of spin states along x direction as a superposition of spin states along z direction, can done by applying appropriate rotation transformation on the spin along z states which provides appropriate phase of the states relative to each other.

An example of finite dimensional Hilbert space can be constructed using spin eigenkets of ${\textstyle s}$-spin particles which forms a ${\textstyle 2s+1}$ dimensional Hilbert space. However, the general wavefunction of a particle that fully describes its state, is always from an infinite dimensional Hilbert space since it involves a tensor product with Hilbert space relating to the position or momentum of the particle. Nonetheless, the techniques developed for finite dimensional Hilbert space are useful since they can either be treated independently or treated in consideration of linearity of tensor product.

Since the spin operator for a given ${\textstyle s}$-spin particles can be represented as a finite ${\textstyle (2s+1)^{2}}$ matrix which acts on ${\textstyle 2s+1}$ independent spin vector components, it is usually preferable to denote spin components using matrix/column/row notation as applicable.

For example, each |s_z⟩ is usually identified as a column vector:$${\displaystyle |s\rangle \leftrightarrow {\begin{bmatrix}1\\0\\\vdots \\0\\0\\\end{bmatrix}}\,,\quad |s-1\rangle \leftrightarrow {\begin{bmatrix}0\\1\\\vdots \\0\\0\\\end{bmatrix}}\,,\ldots \,,\quad |-(s-1)\rangle \leftrightarrow {\begin{bmatrix}0\\0\\\vdots \\1\\0\\\end{bmatrix}}\,,\quad |-s\rangle \leftrightarrow {\begin{bmatrix}0\\0\\\vdots \\0\\1\\\end{bmatrix}}}$$

but it is a common abuse of notation, because the kets |s_z⟩ are not synonymous or equal to the column vectors. Column vectors simply provide a convenient way to express the spin components.

Corresponding to the notation, the z-component spin operator can be written as:$${\displaystyle {\frac {1}{\hbar }}{\hat {S}}_{z}={\begin{bmatrix}s&0&\cdots &0&0\\0&s-1&\cdots &0&0\\\vdots &\vdots &\ddots &\vdots &\vdots \\0&0&\cdots &-(s-1)&0\\0&0&\cdots &0&-s\end{bmatrix}}}$$

since the eigenvectors of z-component spin operator are the above column vectors, with eigenvalues being the corresponding spin quantum numbers.

Corresponding to the notation, a vector from such a finite dimensional Hilbert space is hence represented as:

$${\displaystyle |\phi \rangle ={\begin{bmatrix}\langle s|\phi \rangle \\\langle s-1|\phi \rangle \\\vdots \\\langle -(s-1)|\phi \rangle \\\langle -s|\phi \rangle \\\end{bmatrix}}={\begin{bmatrix}\varepsilon _{s}\\\varepsilon _{s-1}\\\vdots \\\varepsilon _{-s+1}\\\varepsilon _{-s}\\\end{bmatrix}}}$$where ${\textstyle \{\varepsilon _{i}\}}$ are corresponding complex numbers.

In the following discussion involving spin, the complete wavefunction is considered as tensor product of spin states from finite dimensional Hilbert spaces and the wavefunction which was previously developed. The basis for this Hilbert space are hence considered: ${\displaystyle |\mathbf {r} ,s_{z}\rangle =|\mathbf {r} \rangle |s_{z}\rangle }$.

The position-space wave function of a single particle without spin in three spatial dimensions is similar to the case of one spatial dimension above: $${\displaystyle \Psi (\mathbf {r} ,t)}$$ where r is the position vector in three-dimensional space, and t is time. As always Ψ(r, t) is a complex-valued function of real variables. As a single vector in Dirac notation $${\displaystyle |\Psi (t)\rangle =\int d^{3}\!\mathbf {r} \,\Psi (\mathbf {r} ,t)\,|\mathbf {r} \rangle }$$

All the previous remarks on inner products, momentum space wave functions, Fourier transforms, and so on extend to higher dimensions.

For a particle with spin, ignoring the position degrees of freedom, the wave function is a function of spin only (time is a parameter); $${\displaystyle \xi (s_{z},t)}$$ where s_z is the spin projection quantum number along the z axis. (The z axis is an arbitrary choice; other axes can be used instead if the wave function is transformed appropriately, see below.) The s_z parameter, unlike r and t, is a discrete variable. For example, for a spin-1/2 particle, s_z can only be +1/2 or −1/2, and not any other value. (In general, for spin s, s_z can be s, s − 1, ..., −s + 1, −s). Inserting each quantum number gives a complex valued function of space and time, there are 2s + 1 of them. These can be arranged into a column vector

$${\displaystyle \xi ={\begin{bmatrix}\xi (s,t)\\\xi (s-1,t)\\\vdots \\\xi (-(s-1),t)\\\xi (-s,t)\\\end{bmatrix}}=\xi (s,t){\begin{bmatrix}1\\0\\\vdots \\0\\0\\\end{bmatrix}}+\xi (s-1,t){\begin{bmatrix}0\\1\\\vdots \\0\\0\\\end{bmatrix}}+\cdots +\xi (-(s-1),t){\begin{bmatrix}0\\0\\\vdots \\1\\0\\\end{bmatrix}}+\xi (-s,t){\begin{bmatrix}0\\0\\\vdots \\0\\1\\\end{bmatrix}}}$$

In bra–ket notation, these easily arrange into the components of a vector: $${\displaystyle |\xi (t)\rangle =\sum _{s_{z}=-s}^{s}\xi (s_{z},t)\,|s_{z}\rangle }$$

The entire vector ξ is a solution of the Schrödinger equation (with a suitable Hamiltonian), which unfolds to a coupled system of 2s + 1 ordinary differential equations with solutions ξ(s, t), ξ(s − 1, t), ..., ξ(−s, t). The term "spin function" instead of "wave function" is used by some authors. This contrasts the solutions to position space wave functions, the position coordinates being continuous degrees of freedom, because then the Schrödinger equation does take the form of a wave equation.

More generally, for a particle in 3d with any spin, the wave function can be written in "position–spin space" as: $${\displaystyle \Psi (\mathbf {r} ,s_{z},t)}$$ and these can also be arranged into a column vector $${\displaystyle \Psi (\mathbf {r} ,t)={\begin{bmatrix}\Psi (\mathbf {r} ,s,t)\\\Psi (\mathbf {r} ,s-1,t)\\\vdots \\\Psi (\mathbf {r} ,-(s-1),t)\\\Psi (\mathbf {r} ,-s,t)\\\end{bmatrix}}}$$ in which the spin dependence is placed in indexing the entries, and the wave function is a complex vector-valued function of space and time only.

All values of the wave function, not only for discrete but continuous variables also, collect into a single vector $${\displaystyle |\Psi (t)\rangle =\sum _{s_{z}}\int d^{3}\!\mathbf {r} \,\Psi (\mathbf {r} ,s_{z},t)\,|\mathbf {r} ,s_{z}\rangle }$$

For a single particle, the tensor product ⊗ of its position state vector |ψ⟩ and spin state vector |ξ⟩ gives the composite position-spin state vector $${\displaystyle |\psi (t)\rangle \!\otimes \!|\xi (t)\rangle =\sum _{s_{z}}\int d^{3}\!\mathbf {r} \,\psi (\mathbf {r} ,t)\,\xi (s_{z},t)\,|\mathbf {r} \rangle \!\otimes \!|s_{z}\rangle }$$ with the identifications $${\displaystyle |\Psi (t)\rangle =|\psi (t)\rangle \!\otimes \!|\xi (t)\rangle }$$ $${\displaystyle \Psi (\mathbf {r} ,s_{z},t)=\psi (\mathbf {r} ,t)\,\xi (s_{z},t)}$$ $${\displaystyle |\mathbf {r} ,s_{z}\rangle =|\mathbf {r} \rangle \!\otimes \!|s_{z}\rangle }$$

The tensor product factorization of energy eigenstates is always possible if the orbital and spin angular momenta of the particle are separable in the Hamiltonian operator underlying the system's dynamics (in other words, the Hamiltonian can be split into the sum of orbital and spin terms^[34]). The time dependence can be placed in either factor, and time evolution of each can be studied separately. Under such Hamiltonians, any tensor product state evolves into another tensor product state, which essentially means any unentangled state remains unentangled under time evolution. This is said to happen when there is no physical interaction between the states of the tensor products. In the case of non separable Hamiltonians, energy eigenstates are said to be some linear combination of such states, which need not be factorizable; examples include a particle in a magnetic field, and spin–orbit coupling.

The preceding discussion is not limited to spin as a discrete variable, the total angular momentum J may also be used.^[35] Other discrete degrees of freedom, like isospin, can expressed similarly to the case of spin above.

If there are many particles, in general there is only one wave function, not a separate wave function for each particle. The fact that one wave function describes many particles is what makes quantum entanglement and the EPR paradox possible. The position-space wave function for N particles is written:^[20] $${\displaystyle \Psi (\mathbf {r} _{1},\mathbf {r} _{2}\cdots \mathbf {r} _{N},t)}$$ where r_i is the position of the i-th particle in three-dimensional space, and t is time. Altogether, this is a complex-valued function of 3N + 1 real variables.

In quantum mechanics there is a fundamental distinction between identical particles and distinguishable particles. For example, any two electrons are identical and fundamentally indistinguishable from each other; the laws of physics make it impossible to "stamp an identification number" on a certain electron to keep track of it.^[30] This translates to a requirement on the wave function for a system of identical particles: $${\displaystyle \Psi \left(\ldots \mathbf {r} _{a},\ldots ,\mathbf {r} _{b},\ldots \right)=\pm \Psi \left(\ldots \mathbf {r} _{b},\ldots ,\mathbf {r} _{a},\ldots \right)}$$ where the + sign occurs if the particles are all bosons and − sign if they are all fermions. In other words, the wave function is either totally symmetric in the positions of bosons, or totally antisymmetric in the positions of fermions.^[36] The physical interchange of particles corresponds to mathematically switching arguments in the wave function. The antisymmetry feature of fermionic wave functions leads to the Pauli principle. Generally, bosonic and fermionic symmetry requirements are the manifestation of particle statistics and are present in other quantum state formalisms.

For N distinguishable particles (no two being identical, i.e. no two having the same set of quantum numbers), there is no requirement for the wave function to be either symmetric or antisymmetric.

For a collection of particles, some identical with coordinates r₁, r₂, ... and others distinguishable x₁, x₂, ... (not identical with each other, and not identical to the aforementioned identical particles), the wave function is symmetric or antisymmetric in the identical particle coordinates r_i only: $${\displaystyle \Psi \left(\ldots \mathbf {r} _{a},\ldots ,\mathbf {r} _{b},\ldots ,\mathbf {x} _{1},\mathbf {x} _{2},\ldots \right)=\pm \Psi \left(\ldots \mathbf {r} _{b},\ldots ,\mathbf {r} _{a},\ldots ,\mathbf {x} _{1},\mathbf {x} _{2},\ldots \right)}$$

Again, there is no symmetry requirement for the distinguishable particle coordinates x_i.

The wave function for N particles each with spin is the complex-valued function $${\displaystyle \Psi (\mathbf {r} _{1},\mathbf {r} _{2}\cdots \mathbf {r} _{N},s_{z\,1},s_{z\,2}\cdots s_{z\,N},t)}$$

Accumulating all these components into a single vector,

$${\displaystyle |\Psi \rangle =\overbrace {\sum _{s_{z\,1},\ldots ,s_{z\,N}}} ^{\text{discrete labels}}\overbrace {\int _{R_{N}}d^{3}\mathbf {r} _{N}\cdots \int _{R_{1}}d^{3}\mathbf {r} _{1}} ^{\text{continuous labels}}\;\underbrace {{\Psi }(\mathbf {r} _{1},\ldots ,\mathbf {r} _{N},s_{z\,1},\ldots ,s_{z\,N})} _{\begin{array}{c}{\text{wave function (component of }}\\{\text{ state vector along basis state)}}\end{array}}\;\underbrace {|\mathbf {r} _{1},\ldots ,\mathbf {r} _{N},s_{z\,1},\ldots ,s_{z\,N}\rangle } _{\text{basis state (basis ket)}}\,.}$$

For identical particles, symmetry requirements apply to both position and spin arguments of the wave function so it has the overall correct symmetry.

The formulae for the inner products are integrals over all coordinates or momenta and sums over all spin quantum numbers. For the general case of N particles with spin in 3-d, $${\displaystyle (\Psi _{1},\Psi _{2})=\sum _{s_{z\,N}}\cdots \sum _{s_{z\,2}}\sum _{s_{z\,1}}\int \limits _{\mathrm {all\,space} }d^{3}\mathbf {r} _{1}\int \limits _{\mathrm {all\,space} }d^{3}\mathbf {r} _{2}\cdots \int \limits _{\mathrm {all\,space} }d^{3}\mathbf {r} _{N}\Psi _{1}^{*}\left(\mathbf {r} _{1}\cdots \mathbf {r} _{N},s_{z\,1}\cdots s_{z\,N},t\right)\Psi _{2}\left(\mathbf {r} _{1}\cdots \mathbf {r} _{N},s_{z\,1}\cdots s_{z\,N},t\right)}$$ this is altogether N three-dimensional volume integrals and N sums over the spins. The differential volume elements d³r_i are also written "dV_i" or "dx_i dy_i dz_i".

The multidimensional Fourier transforms of the position or position–spin space wave functions yields momentum or momentum–spin space wave functions.

For the general case of N particles with spin in 3d, if Ψ is interpreted as a probability amplitude, the probability density is $${\displaystyle \rho \left(\mathbf {r} _{1}\cdots \mathbf {r} _{N},s_{z\,1}\cdots s_{z\,N},t\right)=\left|\Psi \left(\mathbf {r} _{1}\cdots \mathbf {r} _{N},s_{z\,1}\cdots s_{z\,N},t\right)\right|^{2}}$$

and the probability that particle 1 is in region R₁ with spin s_z1 = m₁ and particle 2 is in region R₂ with spin s_z2 = m₂ etc. at time t is the integral of the probability density over these regions and evaluated at these spin numbers:

${\displaystyle P_{\mathbf {r} _{1}\in R_{1},s_{z\,1}=m_{1},\ldots ,\mathbf {r} _{N}\in R_{N},s_{z\,N}=m_{N}}(t)=\int _{R_{1}}d^{3}\mathbf {r} _{1}\int _{R_{2}}d^{3}\mathbf {r} _{2}\cdots \int _{R_{N}}d^{3}\mathbf {r} _{N}\left|\Psi \left(\mathbf {r} _{1}\cdots \mathbf {r} _{N},m_{1}\cdots m_{N},t\right)\right|^{2}}$

In non-relativistic quantum mechanics, it can be shown using Schrodinger's time dependent wave equation that the equation:

$${\displaystyle {\frac {\partial \rho }{\partial t}}+\nabla \cdot \mathbf {J} =0}$$is satisfied, where ${\textstyle \rho (\mathbf {x} ,t)=|\psi (\mathbf {x} ,t)|^{2}}$ is the probability density and ${\textstyle \mathbf {J} (\mathbf {x} ,t)={\frac {\hbar }{2im}}(\psi ^{*}\nabla \psi -\psi \nabla \psi ^{*})={\frac {\hbar }{m}}{\text{Im}}(\psi ^{*}\nabla \psi )}$, is known as the probability flux in accordance with the continuity equation form of the above equation.

Using the following expression for wavefunction:$${\displaystyle \psi (\mathbf {x} ,t)={\sqrt {\rho (\mathbf {x} ,t)}}\exp {\frac {iS(\mathbf {x} ,t)}{\hbar }}}$$where ${\textstyle \rho (\mathbf {x} ,t)=|\psi (\mathbf {x} ,t)|^{2}}$ is the probability density and ${\textstyle S(\mathbf {x} ,t)}$ is the phase of the wavefunction, it can be shown that:

$${\displaystyle \mathbf {J} (\mathbf {x} ,t)={\frac {\rho \nabla S}{m}}}$$

Hence the spacial variation of phase characterizes the probability flux.

In classical analogy, for ${\textstyle \mathbf {J} =\rho \mathbf {v} }$, the quantity ${\textstyle {\frac {\nabla S}{m}}}$ is analogous with velocity. Note that this does not imply a literal interpretation of ${\textstyle {\frac {\nabla S}{m}}}$ as velocity since velocity and position cannot be simultaneously determined as per the uncertainty principle. Substituting the form of wavefunction in Schrodinger's time dependent wave equation, and taking the classical limit, ${\textstyle \hbar |\nabla ^{2}S|\ll |\nabla S|^{2}}$:

$${\displaystyle {\frac {1}{2m}}|\nabla S(\mathbf {x} ,t)|^{2}+V(\mathbf {x} )+{\frac {\partial S}{\partial t}}=0}$$

Which is analogous to Hamilton-Jacobi equation from classical mechanics. This interpretation fits with Hamilton–Jacobi theory, in which ${\textstyle \mathbf {P} _{\text{class.}}=\nabla S}$, where S is Hamilton's principal function.^[37]

For systems in time-independent potentials, the wave function can always be written as a function of the degrees of freedom multiplied by a time-dependent phase factor, the form of which is given by the Schrödinger equation. For N particles, considering their positions only and suppressing other degrees of freedom, $${\displaystyle \Psi (\mathbf {r} _{1},\mathbf {r} _{2},\ldots ,\mathbf {r} _{N},t)=e^{-iEt/\hbar }\,\psi (\mathbf {r} _{1},\mathbf {r} _{2},\ldots ,\mathbf {r} _{N})\,,}$$ where E is the energy eigenvalue of the system corresponding to the eigenstate Ψ. Wave functions of this form are called stationary states.

The time dependence of the quantum state and the operators can be placed according to unitary transformations on the operators and states. For any quantum state |Ψ⟩ and operator O, in the Schrödinger picture |Ψ(t)⟩ changes with time according to the Schrödinger equation while O is constant. In the Heisenberg picture it is the other way round, |Ψ⟩ is constant while O(t) evolves with time according to the Heisenberg equation of motion. The Dirac (or interaction) picture is intermediate, time dependence is places in both operators and states which evolve according to equations of motion. It is useful primarily in computing S-matrix elements.^[38]

The following are solutions to the Schrödinger equation for one non-relativistic spinless particle.

One of the most prominent features of wave mechanics is the possibility for a particle to reach a location with a prohibitive (in classical mechanics) force potential. A common model is the "potential barrier", the one-dimensional case has the potential $${\displaystyle V(x)={\begin{cases}V_{0}&|x|<a\\0&|x|\geq a\end{cases}}}$$ and the steady-state solutions to the wave equation have the form (for some constants k, κ) $${\displaystyle \Psi (x)={\begin{cases}A_{\mathrm {r} }e^{ikx}+A_{\mathrm {l} }e^{-ikx}&x<-a,\\B_{\mathrm {r} }e^{\kappa x}+B_{\mathrm {l} }e^{-\kappa x}&|x|\leq a,\\C_{\mathrm {r} }e^{ikx}+C_{\mathrm {l} }e^{-ikx}&x>a.\end{cases}}}$$

Note that these wave functions are not normalized; see scattering theory for discussion.

The standard interpretation of this is as a stream of particles being fired at the step from the left (the direction of negative x): setting A_r = 1 corresponds to firing particles singly; the terms containing A_r and C_r signify motion to the right, while A_l and C_l – to the left. Under this beam interpretation, put C_l = 0 since no particles are coming from the right. By applying the continuity of wave functions and their derivatives at the boundaries, it is hence possible to determine the constants above.

In a semiconductor crystallite whose radius is smaller than the size of its exciton Bohr radius, the excitons are squeezed, leading to quantum confinement. The energy levels can then be modeled using the particle in a box model in which the energy of different states is dependent on the length of the box.

The wave functions for the quantum harmonic oscillator can be expressed in terms of Hermite polynomials H_n, they are $${\displaystyle \Psi _{n}(x)={\sqrt {\frac {1}{2^{n}\,n!}}}\cdot \left({\frac {m\omega }{\pi \hbar }}\right)^{1/4}\cdot e^{-{\frac {m\omega x^{2}}{2\hbar }}}\cdot H_{n}{\left({\sqrt {\frac {m\omega }{\hbar }}}x\right)}}$$ where n = 0, 1, 2, ....

The wave functions of an electron in a Hydrogen atom are expressed in terms of spherical harmonics and generalized Laguerre polynomials (these are defined differently by different authors—see main article on them and the hydrogen atom).

It is convenient to use spherical coordinates, and the wave function can be separated into functions of each coordinate,^[39] $${\displaystyle \Psi _{n\ell m}(r,\theta ,\phi )=R(r)\,\,Y_{\ell }^{m}\!(\theta ,\phi )}$$ where R are radial functions and Y^m
_ℓ(θ, φ) are spherical harmonics of degree ℓ and order m. This is the only atom for which the Schrödinger equation has been solved exactly. Multi-electron atoms require approximative methods. The family of solutions is:^[40] $${\displaystyle \Psi _{n\ell m}(r,\theta ,\phi )={\sqrt {{\left({\frac {2}{na_{0}}}\right)}^{3}{\frac {(n-\ell -1)!}{2n[(n+\ell )!]}}}}e^{-r/na_{0}}\left({\frac {2r}{na_{0}}}\right)^{\ell }L_{n-\ell -1}^{2\ell +1}\left({\frac {2r}{na_{0}}}\right)\cdot Y_{\ell }^{m}(\theta ,\phi )}$$ where a₀ = 4πε₀ħ²/m_ee² is the Bohr radius, L^{2ℓ + 1}
_{n − ℓ − 1} are the generalized Laguerre polynomials of degree n − ℓ − 1, n = 1, 2, ... is the principal quantum number, ℓ = 0, 1, ..., n − 1 the azimuthal quantum number, m = −ℓ, −ℓ + 1, ..., ℓ − 1, ℓ the magnetic quantum number. Hydrogen-like atoms have very similar solutions.

This solution does not take into account the spin of the electron.

In the figure of the hydrogen orbitals, the 19 sub-images are images of wave functions in position space (their norm squared). The wave functions represent the abstract state characterized by the triple of quantum numbers (n, ℓ, m), in the lower right of each image. These are the principal quantum number, the orbital angular momentum quantum number, and the magnetic quantum number. Together with one spin-projection quantum number of the electron, this is a complete set of observables.

The figure can serve to illustrate some further properties of the function spaces of wave functions.

• In this case, the wave functions are square integrable. One can initially take the function space as the space of square integrable functions, usually denoted L².
• The displayed functions are solutions to the Schrödinger equation. Obviously, not every function in L² satisfies the Schrödinger equation for the hydrogen atom. The function space is thus a subspace of L².
• The displayed functions form part of a basis for the function space. To each triple (n, ℓ, m), there corresponds a basis wave function. If spin is taken into account, there are two basis functions for each triple. The function space thus has a countable basis.
• The basis functions are mutually orthonormal.

The concept of function spaces enters naturally in the discussion about wave functions. A function space is a set of functions, usually with some defining requirements on the functions (in the present case that they are square integrable), sometimes with an algebraic structure on the set (in the present case a vector space structure with an inner product), together with a topology on the set. The latter will sparsely be used here, it is only needed to obtain a precise definition of what it means for a subset of a function space to be closed. It will be concluded below that the function space of wave functions is a Hilbert space. This observation is the foundation of the predominant mathematical formulation of quantum mechanics.

A wave function is an element of a function space partly characterized by the following concrete and abstract descriptions.

• The Schrödinger equation is linear. This means that the solutions to it, wave functions, can be added and multiplied by scalars to form a new solution. The set of solutions to the Schrödinger equation is a vector space.
• The superposition principle of quantum mechanics. If Ψ and Φ are two states in the abstract space of states of a quantum mechanical system, and a and b are any two complex numbers, then aΨ + bΦ is a valid state as well. (Whether the null vector counts as a valid state ("no system present") is a matter of definition. The null vector does not at any rate describe the vacuum state in quantum field theory.) The set of allowable states is a vector space.

This similarity is of course not accidental. There are also a distinctions between the spaces to keep in mind.

Basic states are characterized by a set of quantum numbers. This is a set of eigenvalues of a maximal set of commuting observables. Physical observables are represented by linear operators, also called observables, on the vectors space. Maximality means that there can be added to the set no further algebraically independent observables that commute with the ones already present. A choice of such a set may be called a choice of representation.

• It is a postulate of quantum mechanics that a physically observable quantity of a system, such as position, momentum, or spin, is represented by a linear Hermitian operator on the state space. The possible outcomes of measurement of the quantity are the eigenvalues of the operator.^[18] At a deeper level, most observables, perhaps all, arise as generators of symmetries.^{[18][41][nb 6]}
• The physical interpretation is that such a set represents what can – in theory – simultaneously be measured with arbitrary precision. The Heisenberg uncertainty relation prohibits simultaneous exact measurements of two non-commuting observables.
• The set is non-unique. It may for a one-particle system, for example, be position and spin z-projection, (x, S_z), or it may be momentum and spin y-projection, (p, S_y). In this case, the operator corresponding to position (a multiplication operator in the position representation) and the operator corresponding to momentum (a differential operator in the position representation) do not commute.
• Once a representation is chosen, there is still arbitrariness. It remains to choose a coordinate system. This may, for example, correspond to a choice of x, y- and z-axis, or a choice of curvilinear coordinates as exemplified by the spherical coordinates used for the Hydrogen atomic wave functions. This final choice also fixes a basis in abstract Hilbert space. The basic states are labeled by the quantum numbers corresponding to the maximal set of commuting observables and an appropriate coordinate system.^{[nb 7]}

The abstract states are "abstract" only in that an arbitrary choice necessary for a particular explicit description of it is not given. This is the same as saying that no choice of maximal set of commuting observables has been given. This is analogous to a vector space without a specified basis. Wave functions corresponding to a state are accordingly not unique. This non-uniqueness reflects the non-uniqueness in the choice of a maximal set of commuting observables. For one spin particle in one dimension, to a particular state there corresponds two wave functions, Ψ(x, S_z) and Ψ(p, S_y), both describing the same state.

• For each choice of maximal commuting sets of observables for the abstract state space, there is a corresponding representation that is associated to a function space of wave functions.
• Between all these different function spaces and the abstract state space, there are one-to-one correspondences (here disregarding normalization and unobservable phase factors), the common denominator here being a particular abstract state. The relationship between the momentum and position space wave functions, for instance, describing the same state is the Fourier transform.

Each choice of representation should be thought of as specifying a unique function space in which wave functions corresponding to that choice of representation lives. This distinction is best kept, even if one could argue that two such function spaces are mathematically equal, e.g. being the set of square integrable functions. One can then think of the function spaces as two distinct copies of that set.

There is an additional algebraic structure on the vector spaces of wave functions and the abstract state space.

• Physically, different wave functions are interpreted to overlap to some degree. A system in a state Ψ that does not overlap with a state Φ cannot be found to be in the state Φ upon measurement. But if Φ₁, Φ₂, … overlap Ψ to some degree, there is a chance that measurement of a system described by Ψ will be found in states Φ₁, Φ₂, …. Also selection rules are observed apply. These are usually formulated in the preservation of some quantum numbers. This means that certain processes allowable from some perspectives (e.g. energy and momentum conservation) do not occur because the initial and final total wave functions do not overlap.
• Mathematically, it turns out that solutions to the Schrödinger equation for particular potentials are orthogonal in some manner, this is usually described by an integral $${\displaystyle \int \Psi _{m}^{*}\Psi _{n}w\,dV=\delta _{nm},}$$ where m, n are (sets of) indices (quantum numbers) labeling different solutions, the strictly positive function w is called a weight function, and δ_mn is the Kronecker delta. The integration is taken over all of the relevant space.

This motivates the introduction of an inner product on the vector space of abstract quantum states, compatible with the mathematical observations above when passing to a representation. It is denoted (Ψ, Φ), or in the Bra–ket notation ⟨Ψ|Φ⟩. It yields a complex number. With the inner product, the function space is an inner product space. The explicit appearance of the inner product (usually an integral or a sum of integrals) depends on the choice of representation, but the complex number (Ψ, Φ) does not. Much of the physical interpretation of quantum mechanics stems from the Born rule. It states that the probability p of finding upon measurement the state Φ given the system is in the state Ψ is $${\displaystyle p=|(\Phi ,\Psi )|^{2},}$$ where Φ and Ψ are assumed normalized. Consider a scattering experiment. In quantum field theory, if Φ_out describes a state in the "distant future" (an "out state") after interactions between scattering particles have ceased, and Ψ_in an "in state" in the "distant past", then the quantities (Φ_out, Ψ_in), with Φ_out and Ψ_in varying over a complete set of in states and out states respectively, is called the S-matrix or scattering matrix. Knowledge of it is, effectively, having solved the theory at hand, at least as far as predictions go. Measurable quantities such as decay rates and scattering cross sections are calculable from the S-matrix.^[42]

The above observations encapsulate the essence of the function spaces of which wave functions are elements. However, the description is not yet complete. There is a further technical requirement on the function space, that of completeness, that allows one to take limits of sequences in the function space, and be ensured that, if the limit exists, it is an element of the function space. A complete inner product space is called a Hilbert space. The property of completeness is crucial in advanced treatments and applications of quantum mechanics. For instance, the existence of projection operators or orthogonal projections relies on the completeness of the space.^[43] These projection operators, in turn, are essential for the statement and proof of many useful theorems, e.g. the spectral theorem. It is not very important in introductory quantum mechanics, and technical details and links may be found in footnotes like the one that follows.^{[nb 8]} The space L² is a Hilbert space, with inner product presented later. The function space of the example of the figure is a subspace of L². A subspace of a Hilbert space is a Hilbert space if it is closed.

In summary, the set of all possible normalizable wave functions for a system with a particular choice of basis, together with the null vector, constitute a Hilbert space.

Not all functions of interest are elements of some Hilbert space, say L². The most glaring example is the set of functions e^{2πip · x⁄h}. These are plane wave solutions of the Schrödinger equation for a free particle, but are not normalizable, hence not in L². But they are nonetheless fundamental for the description. One can, using them, express functions that are normalizable using wave packets. They are, in a sense, a basis (but not a Hilbert space basis, nor a Hamel basis) in which wave functions of interest can be expressed. There is also the artifact "normalization to a delta function" that is frequently employed for notational convenience, see further down. The delta functions themselves are not square integrable either.

The above description of the function space containing the wave functions is mostly mathematically motivated. The function spaces are, due to completeness, very large in a certain sense. Not all functions are realistic descriptions of any physical system. For instance, in the function space L² one can find the function that takes on the value 0 for all rational numbers and -i for the irrationals in the interval [0, 1]. This is square integrable,^{[nb 9]} but can hardly represent a physical state.

While the space of solutions as a whole is a Hilbert space there are many other Hilbert spaces that commonly occur as ingredients.

• Square integrable complex valued functions on the interval [0, 2π]. The set {e^int/2π, n ∈ Z} is a Hilbert space basis, i.e. a maximal orthonormal set.
• The Fourier transform takes functions in the above space to elements of l²(Z), the space of square summable functions Z → C. The latter space is a Hilbert space and the Fourier transform is an isomorphism of Hilbert spaces.^{[nb 10]} Its basis is {e_i, i ∈ Z} with e_i(j) = δ_ij, i, j ∈ Z.
• The most basic example of spanning polynomials is in the space of square integrable functions on the interval [–1, 1] for which the Legendre polynomials is a Hilbert space basis (complete orthonormal set).
• The square integrable functions on the unit sphere S² is a Hilbert space. The basis functions in this case are the spherical harmonics. The Legendre polynomials are ingredients in the spherical harmonics. Most problems with rotational symmetry will have "the same" (known) solution with respect to that symmetry, so the original problem is reduced to a problem of lower dimensionality.
• The associated Laguerre polynomials appear in the hydrogenic wave function problem after factoring out the spherical harmonics. These span the Hilbert space of square integrable functions on the semi-infinite interval [0, ∞).

More generally, one may consider a unified treatment of all second order polynomial solutions to the Sturm–Liouville equations in the setting of Hilbert space. These include the Legendre and Laguerre polynomials as well as Chebyshev polynomials, Jacobi polynomials and Hermite polynomials. All of these actually appear in physical problems, the latter ones in the harmonic oscillator, and what is otherwise a bewildering maze of properties of special functions becomes an organized body of facts. For this, see Byron & Fuller (1992, Chapter 5).

There occurs also finite-dimensional Hilbert spaces. The space Cⁿ is a Hilbert space of dimension n. The inner product is the standard inner product on these spaces. In it, the "spin part" of a single particle wave function resides.

• In the non-relativistic description of an electron one has n = 2 and the total wave function is a solution of the Pauli equation.
• In the corresponding relativistic treatment, n = 4 and the wave function solves the Dirac equation.

With more particles, the situations is more complicated. One has to employ tensor products and use representation theory of the symmetry groups involved (the rotation group and the Lorentz group respectively) to extract from the tensor product the spaces in which the (total) spin wave functions reside. (Further problems arise in the relativistic case unless the particles are free.^[44] See the Bethe–Salpeter equation.) Corresponding remarks apply to the concept of isospin, for which the symmetry group is SU(2). The models of the nuclear forces of the sixties (still useful today, see nuclear force) used the symmetry group SU(3). In this case, as well, the part of the wave functions corresponding to the inner symmetries reside in some Cⁿ or subspaces of tensor products of such spaces.

• In quantum field theory the underlying Hilbert space is Fock space. It is built from free single-particle states, i.e. wave functions when a representation is chosen, and can accommodate any finite, not necessarily constant in time, number of particles. The interesting (or rather the tractable) dynamics lies not in the wave functions but in the field operators that are operators acting on Fock space. Thus the Heisenberg picture is the most common choice (constant states, time varying operators).

Due to the infinite-dimensional nature of the system, the appropriate mathematical tools are objects of study in functional analysis.

Not all introductory textbooks take the long route and introduce the full Hilbert space machinery, but the focus is on the non-relativistic Schrödinger equation in position representation for certain standard potentials. The following constraints on the wave function are sometimes explicitly formulated for the calculations and physical interpretation to make sense:^[45][46]

• The wave function must be square integrable. This is motivated by the Copenhagen interpretation of the wave function as a probability amplitude.
• It must be everywhere continuous and everywhere continuously differentiable. This is motivated by the appearance of the Schrödinger equation for most physically reasonable potentials.

It is possible to relax these conditions somewhat for special purposes.^{[nb 11]} If these requirements are not met, it is not possible to interpret the wave function as a probability amplitude.^[47] Note that exceptions can arise to the continuity of derivatives rule at points of infinite discontinuity of potential field. For example, in particle in a box where the derivative of wavefunction can be discontinuous at the boundary of the box where the potential is known to have infinite discontinuity.

This does not alter the structure of the Hilbert space that these particular wave functions inhabit, but the subspace of the square-integrable functions L², which is a Hilbert space, satisfying the second requirement is not closed in L², hence not a Hilbert space in itself.^{[nb 12]} The functions that does not meet the requirements are still needed for both technical and practical reasons.^{[nb 13][nb 14]}

As has been demonstrated, the set of all possible wave functions in some representation for a system constitute an in general infinite-dimensional Hilbert space. Due to the multiple possible choices of representation basis, these Hilbert spaces are not unique. One therefore talks about an abstract Hilbert space, state space, where the choice of representation and basis is left undetermined. Specifically, each state is represented as an abstract vector in state space.^[48] A quantum state |Ψ⟩ in any representation is generally expressed as a vector $${\displaystyle |\Psi \rangle =\sum _{\boldsymbol {\alpha }}\int d^{m}\!{\boldsymbol {\omega }}\,\,\Psi _{t}({\boldsymbol {\alpha }},{\boldsymbol {\omega }})\,|{\boldsymbol {\alpha }},{\boldsymbol {\omega }}\rangle }$$ where

• |α, ω⟩ the basis vectors of the chosen representation
• d^mω = dω₁dω₂...dω_m a differential volume element in the continuous degrees of freedom
• ${\displaystyle {\boldsymbol {\Psi }}_{t}({\boldsymbol {\alpha }},{\boldsymbol {\omega }})}$ a component of the vector ${\displaystyle |\Psi \rangle }$, called the wave function of the system
• α = (α₁, α₂, ..., α_n) dimensionless discrete quantum numbers
• ω = (ω₁, ω₂, ..., ω_m) continuous variables (not necessarily dimensionless)

These quantum numbers index the components of the state vector. More, all α are in an n-dimensional set A = A₁ × A₂ × ... × A_n where each A_i is the set of allowed values for α_i; all ω are in an m-dimensional "volume" Ω ⊆ ℝ^m where Ω = Ω₁ × Ω₂ × ... × Ω_m and each Ω_i ⊆ R is the set of allowed values for ω_i, a subset of the real numbers R. For generality n and m are not necessarily equal.

Example:

The probability density of finding the system at time ${\displaystyle t}$ at state |α, ω⟩ is $${\displaystyle \rho _{\alpha ,\omega }(t)=|\Psi ({\boldsymbol {\alpha }},{\boldsymbol {\omega }},t)|^{2}}$$

The probability of finding system with α in some or all possible discrete-variable configurations, D ⊆ A, and ω in some or all possible continuous-variable configurations, C ⊆ Ω, is the sum and integral over the density,^{[nb 15]} $${\displaystyle P(t)=\sum _{{\boldsymbol {\alpha }}\in D}\int _{C}d^{m}\!{\boldsymbol {\omega }}\,\,\rho _{\alpha ,\omega }(t)}$$

Since the sum of all probabilities must be 1, the normalization condition $${\displaystyle 1=\sum _{{\boldsymbol {\alpha }}\in A}\int _{\Omega }d^{m}\!{\boldsymbol {\omega }}\,\,\rho _{\alpha ,\omega }(t)}$$ must hold at all times during the evolution of the system.

The normalization condition requires ρ d^mω to be dimensionless, by dimensional analysis Ψ must have the same units as (ω₁ω₂...ω_m)^−1/2.

Whether the wave function exists in reality, and what it represents, are major questions in the interpretation of quantum mechanics. Many famous physicists of a previous generation puzzled over this problem, such as Erwin Schrödinger, Albert Einstein and Niels Bohr. Some advocate formulations or variants of the Copenhagen interpretation (e.g. Bohr, Eugene Wigner and John von Neumann) while others, such as John Archibald Wheeler or Edwin Thompson Jaynes, take the more classical approach^[49] and regard the wave function as representing information in the mind of the observer, i.e. a measure of our knowledge of reality. Some, including Schrödinger, David Bohm and Hugh Everett III and others, argued that the wave function must have an objective, physical existence. Einstein thought that a complete description of physical reality should refer directly to physical space and time, as distinct from the wave function, which refers to an abstract mathematical space.^[50]

• Quantum Mechanics for Engineers
• Spin wave functions NYU
• Identical Particles Revisited, Michael Fowler
• The Nature of Many-Electron Wavefunctions
• Quantum Mechanics and Quantum Computation at BerkeleyX Archived 2013-05-13 at the Wayback Machine
• Einstein, The quantum theory of radiation