thumb|upright=1.25|[[Quantum harmonic oscillator|Quantum harmonic oscillators for a single spinless particle. The oscillations have no trajectory, but are instead represented each as waves; the vertical axis shows the real part (blue) and imaginary part (red) of the wave function. Panels A–D show four different standing-wave solutions of the Schrödinger equation. Panels E–F show two different wave functions that are solutions of the Schrödinger equation but not standing waves.]]

thumb|The wave function of an initially very localized free particle

In quantum mechanics, 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).

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 of 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, whether the wave function in quantum mechanics describes a kind of physical phenomenon is still open to different interpretations, fundamentally differentiating it from classic mechanical waves.

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 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 column vector for a non-relativistic electron with spin ).

Historical background

In 1900, Max Planck postulated the proportionality between the frequency <math>f</math> of a photon and its energy

and in 1916 the corresponding relation between a photon's momentum <math>p</math> and wavelength

where <math>h</math> is the Planck constant. In 1923, De Broglie was the first to suggest that the relation now called the De Broglie relation, holds for massive particles, the chief clue being Lorentz invariance, 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.

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. However, no one was clear on how to interpret it. 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. 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.

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. 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. 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.

In 1927, Pauli phenomenologically found a non-relativistic equation to describe spin-1/2 particles in electromagnetic fields, now called the Pauli equation. 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: 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.

Wave functions and wave equations in modern theories

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.

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 ) and the Dirac equation (spin ) in this guise remain in the theory. Higher spin analogues include the Proca equation (spin ), Rarita–Schwinger equation (spin ), and, more generally, the Bargmann–Wigner equations. For massless free fields two examples are the free field Maxwell equation (spin ) and the free field Einstein equation (spin ) for the field operators.

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,

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.

Relativistic particles: classical vs quantum

Source:

Classical particles, like bullets, rockets and planets are spatially localized. Classical waves like sound and ocean waves expand over large spatial regions. Quantum wave-particle duality is a unique characteristic of quantum particles. For instance, spatially localized atoms are stable due to the existence of waves associate to the electrons in the atom. There is no consensus about the nature of the waves associated to quantum particles; however, there is consensus on how to calculate the mathematical expressions corresponding to these waves.

:<math>( \Psi_1 , \Psi_2 ) = \int_{-\infty}^\infty \, \Psi_1^*(x, t)\Psi_2(x, t)\,dx < \infty</math>.

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

:<math>(\Psi,\Psi) = \|\Psi\|^2</math>,

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, 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.

Position-space wave functions

The state of such a particle is completely described by its wave function, <math display="block">\Psi(x,t)\,,</math> where is position and is time. This is a complex-valued function of two real variables and .

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

<math display="block"> \left|\Psi(x, t)\right|^2 = \Psi^*(x, t)\Psi(x, t) = \rho(x), </math>

is interpreted as the probability density for a measurement of the particle's position at a given time . 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.

Normalization condition

The probability that its position will be in the interval is the integral of the density over this interval:

<math display="block">P_{a\le x\le b} (t) = \int_a^b \,|\Psi(x,t)|^2 dx </math>

where is the time at which the particle was measured. This leads to the normalization condition:

<math display="block">\int_{-\infty}^\infty \, |\Psi(x,t)|^2dx = 1\,,</math>

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.

Quantum states as vectors

At a particular instant of time, all values of the wave function 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

<math display="block">|\Psi(t)\rangle = \int\Psi(x,t) |x\rangle dx </math>

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 coordinate is a continuous index. The are called improper vectors which, unlike proper vectors that are normalizable to unity, can only be normalized to a Dirac delta function.

<math display="block">(\Psi_{p},\Psi_{p'}) = \delta(p - p').</math>

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.

Relations between position and momentum representations

The and representations are

<math display="block">\begin{align}

|\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{align}</math>

Now take the projection of the state onto eigenfunctions of momentum using the last expression in the two equations,

<math display="block">\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).</math>

Then utilizing the known expression for suitably normalized eigenstates of momentum in the position representation solutions of the free Schrödinger equation

<math display="block">\langle x | p \rangle = p(x) = \frac{1}{\sqrt{2\pi\hbare^{\frac{i}{\hbar}px} \Rightarrow \langle p | x \rangle = \frac{1}{\sqrt{2\pi\hbare^{-\frac{i}{\hbar}px},</math>

one obtains

<math display="block">\Phi(p) = \frac{1}{\sqrt{2\pi\hbar\int \Psi(x)e^{-\frac{i}{\hbar}px}dx\,.</math>

Likewise, using eigenfunctions of position,

<math display="block">\Psi(x) = \frac{1}{\sqrt{2\pi\hbar\int \Phi(p)e^{\frac{i}{\hbar}px}dp\,.</math>

The position-space and momentum-space wave functions are thus found to be Fourier transforms of each other. 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, and 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 .

Definitions (other cases)

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.

Finite dimensional Hilbert space

While Hilbert spaces originally refer to infinite dimensional complete inner product spaces they, by definition, include finite dimensional complete inner product spaces as well.

In physics, they are often referred to as finite dimensional Hilbert spaces. For every finite dimensional Hilbert space there exist orthonormal basis kets that span the entire Hilbert space.

If the -dimensional set <math display="inline">\{ |\phi_i\rangle \}</math> is orthonormal, then the projection operator for the space spanned by these states is given by:

<math display="block">P = \sum_i |\phi_i\rangle\langle \phi_i | = I </math>where the projection is equivalent to identity operator since <math display="inline">\{ |\phi_i\rangle \}</math> 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:

<math display="block">|\psi\rangle = I|\psi\rangle = \sum_i |\phi_i\rangle\langle \phi_i |\psi\rangle </math>where <math display="inline">\{ \langle \phi_i |\psi\rangle \} </math>, is a set of complex numbers which can be used to construct a wavefunction using the above formula.

Probability interpretation of inner product

If the set <math display="inline">\{ |\phi_i\rangle \}</math> are eigenkets of a non-degenerate observable with eigenvalues <math display="inline">\lambda_i</math>, by the postulates of quantum mechanics, the probability of measuring the observable to be <math display="inline">\lambda_i</math> is given according to Born rule as:

<math display="block">P_\psi(\lambda_i) = |\langle \phi_i|\psi \rangle|^2 </math>

For non-degenerate <math display="inline">\{ |\phi_i\rangle \}</math> of some observable, if eigenvalues <math display="inline">\lambda</math> have subset of eigenvectors labelled as <math display="inline">\{ |\lambda^{(j)}\rangle \}</math>, by the postulates of quantum mechanics, the probability of measuring the observable to be <math display="inline">\lambda</math> is given by:

<math display="block">P_\psi(\lambda) =\sum_j |\langle \lambda^{(j)}|\psi \rangle|^2 = |\widehat P_\lambda |\psi \rangle |^2 </math>where <math display="inline">\widehat P_\lambda =\sum_j|\lambda^{(j)}\rangle\langle\lambda^{(j)}| </math> is a projection operator of states to subspace spanned by <math display="inline">\{ |\lambda^{(j)}\rangle \}</math>. The equality follows due to orthogonal nature of <math display="inline">\{ |\phi_i\rangle \}</math>.

Hence, <math display="inline">\{ \langle \phi_i |\psi\rangle \} </math> which specify state of the quantum mechanical system, have magnitudes whose square gives the probability of measuring the respective <math display="inline">|\phi_i\rangle </math> state.

Physical significance of relative phase

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 <math display="inline">|\phi_i\rangle </math> 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.

Application to include spin

An example of finite dimensional Hilbert space can be constructed using spin eigenkets of <math display="inline">s</math>-spin particles which forms a <math display="inline">2s+1</math> 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 <math display="inline">s</math>-spin particles can be represented as a finite <math display="inline">(2s+1)^2 </math> matrix which acts on <math display="inline">2s+1</math> independent spin vector components, it is usually preferable to denote spin components using matrix/column/row notation as applicable.

For example, each is usually identified as a column vector:<math display="block">|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}</math>

but it is a common abuse of notation, because the kets 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:<math display="block">\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} </math>

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:

<math display="block">|\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} </math>where <math display="inline"> \{ \varepsilon_i \} </math> 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: <math> |\mathbf{r}, s_z\rangle = |\mathbf{r}\rangle |s_z\rangle </math>.

One-particle states in 3d position space

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: <math display="block">\Psi(\mathbf{r},t)</math> where is the position vector in three-dimensional space, and is time. As always is a complex-valued function of real variables. As a single vector in Dirac notation

<math display="block">|\Psi(t)\rangle = \int d^3\! \mathbf{r}\, \Psi(\mathbf{r},t) \,|\mathbf{r}\rangle </math>

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);

<math display="block">\xi(s_z,t)</math>

where is the spin projection quantum number along the axis. (The axis is an arbitrary choice; other axes can be used instead if the wave function is transformed appropriately, see below.) The parameter, unlike and , is a discrete variable. For example, for a spin-1/2 particle, can only be or , and not any other value. (In general, for spin , can be ). Inserting each quantum number gives a complex valued function of space and time, there are of them. These can be arranged into a column vector

<math display="block">\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} </math>

In bra–ket notation, these easily arrange into the components of a vector:

<math display="block">|\xi (t)\rangle = \sum_{s_z=-s}^s \xi(s_z,t) \,| s_z \rangle </math>

The entire vector is a solution of the Schrödinger equation (with a suitable Hamiltonian), which unfolds to a coupled system of ordinary differential equations with solutions . 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:

<math display="block">\Psi(\mathbf{r},s_z,t)</math>

and these can also be arranged into a column vector

<math display="block">\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}</math>

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

<math display="block">|\Psi(t)\rangle =

\sum_{s_z}\int d^3\!\mathbf{r} \,\Psi(\mathbf{r},s_z,t)\, |\mathbf{r}, s_z\rangle </math>

For a single particle, the tensor product of its position state vector and spin state vector gives the composite position-spin state vector

<math display="block">|\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 </math>

with the identifications

<math display="block">|\Psi (t)\rangle = |\psi(t)\rangle

\!\otimes\!

|\xi(t)\rangle </math>

<math display="block">\Psi(\mathbf{r},s_z,t) = \psi(\mathbf{r},t)\,\xi(s_z,t) </math>

<math display="block">|\mathbf{r},s_z \rangle= |\mathbf{r}\rangle \!\otimes\! |s_z\rangle </math>

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). 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. Other discrete degrees of freedom, like isospin, can expressed similarly to the case of spin above.

Many-particle states in 3d position space

upright=1.4|thumb|Traveling waves of two free particles, with two of three dimensions suppressed. Top is position-space wave function, bottom is momentum-space wave function, with corresponding probability densities.

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 particles is written:

<math display="block">\Psi(\mathbf{r}_1,\mathbf{r}_2 \cdots \mathbf{r}_N,t)</math>

where is the position of the -th particle in three-dimensional space, and is time. Altogether, this is a complex-valued function of 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. This translates to a requirement on the wave function for a system of identical particles:

<math display="block">\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 )</math>

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. 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 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 and others distinguishable (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 only:

<math display="block">\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 )</math>

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

The wave function for N particles each with spin is the complex-valued function

<math display="block">\Psi(\mathbf{r}_1, \mathbf{r}_2 \cdots \mathbf{r}_N, s_{z\,1}, s_{z\,2} \cdots s_{z\,N}, t)</math>

Accumulating all these components into a single vector,

<math display="block">| \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 \; \underbrace{| \mathbf{r}_1, \ldots, \mathbf{r}_N , s_{z\,1} , \ldots , s_{z\,N} \rangle }_{\text{basis state (basis ket)\,.</math>

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 particles with spin in 3-d,

<math display="block"> ( \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 ) </math>

this is altogether three-dimensional volume integrals and sums over the spins. The differential volume elements are also written "" or "".

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

Probability interpretation

For the general case of particles with spin in 3d, if is interpreted as a probability amplitude, the probability density is

<math display="block">\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</math>

and the probability that particle 1 is in region with spin and particle 2 is in region with spin etc. at time is the integral of the probability density over these regions and evaluated at these spin numbers:

:<math>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</math>

Physical significance of phase

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

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

Using the following expression for wavefunction:<math display="block">\psi(\mathbf x,t)= \sqrt{\rho(\mathbf x,t)}\exp{\frac{iS(\mathbf x,t )}{\hbar </math>where <math display="inline">\rho(\mathbf x,t) = | \psi(\mathbf x,t)|^2 </math> is the probability density and <math display="inline">S(\mathbf x,t) </math> is the phase of the wavefunction, it can be shown that:

<math display="block">\mathbf J(\mathbf x,t) = \frac{\rho \nabla S}{m} </math>

Hence the spacial variation of phase characterizes the probability flux.

In classical analogy, for <math display="inline">\mathbf J = \rho \mathbf v </math>, the quantity <math display="inline"> \frac{\nabla S}{m} </math> is analogous with velocity. Note that this does not imply a literal interpretation of <math display="inline"> \frac{\nabla S}{m} </math> 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, <math display="inline"> \hbar |\nabla^2 S| \ll |\nabla S|^2 </math>:

<math display="block">\frac{1}{2m} |\nabla S(\mathbf x, t)|^2 + V(\mathbf x) + \frac{\partial S}{\partial t} = 0 </math>

Which is analogous to Hamilton-Jacobi equation from classical mechanics. This interpretation fits with Hamilton–Jacobi theory, in which <math display="inline"> \mathbf{P}_{\text{class. = \nabla S </math>, where ' is Hamilton's principal function.

Time dependence

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 particles, considering their positions only and suppressing other degrees of freedom,

<math display="block">\Psi(\mathbf{r}_1,\mathbf{r}_2,\ldots,\mathbf{r}_N,t) = e^{-i Et/\hbar} \,\psi(\mathbf{r}_1,\mathbf{r}_2,\ldots,\mathbf{r}_N)\,,</math>

where 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 , in the Schrödinger picture changes with time according to the Schrödinger equation while is constant. In the Heisenberg picture it is the other way round, is constant while 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.

Non-relativistic examples

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

Finite potential barrier

thumb|Scattering at a finite potential barrier of height . The amplitudes and direction of left and right moving waves are indicated. In red, those waves used for the derivation of the reflection and transmission amplitude. for this illustration.

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

<math display="block">V(x)=\begin{cases}V_0 & |x|<a \\ 0 & | x | \geq a\end{cases}</math>

and the steady-state solutions to the wave equation have the form (for some constants )

<math display="block">\Psi (x) = \begin{cases}

A_{\mathrm{re^{ikx}+A_{\mathrm{le^{-ikx} & x<-a, \\

B_{\mathrm{re^{\kappa x}+B_{\mathrm{le^{-\kappa x} & |x|\le a, \\

C_{\mathrm{re^{ikx}+C_{\mathrm{le^{-ikx} & x>a.

\end{cases}</math>

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 ): setting corresponds to firing particles singly; the terms containing and signify motion to the right, while and – to the left. Under this beam interpretation, put 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.

thumb|upright=1.3|3D confined electron wave functions in a quantum dot. Here, rectangular and triangular-shaped quantum dots are shown. Energy states in rectangular dots are more s-type and p-type. However, in a triangular dot the wave functions are mixed due to confinement symmetry. (Click for animation)|link=File:QuantumDot_wf.gif

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.

Quantum harmonic oscillator

The wave functions for the quantum harmonic oscillator can be expressed in terms of Hermite polynomials , they are

<math display="block"> \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)} </math>

where .

thumb|upright=1.5|The electron probability density for the first few [[hydrogen atom electron orbitals shown as cross-sections. These orbitals form an orthonormal basis for the wave function of the electron. Different orbitals are depicted with different scale.]]

Hydrogen atom

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,

<math display="block"> \Psi_{n\ell m}(r,\theta,\phi) = R(r)\,\,Y_\ell^m\!(\theta, \phi)</math>

where are radial functions and are spherical harmonics of degree and order . 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:

<math display="block"> \Psi_{n\ell m}(r,\theta,\phi) = \sqrt \int d^m\!\boldsymbol{\omega}\,\,

\Psi_t(\boldsymbol\alpha,\boldsymbol\omega)\,

|\boldsymbol\alpha,\boldsymbol\omega\rangle</math>

where

  • the basis vectors of the chosen representation
  • a differential volume element in the continuous degrees of freedom
  • <math>\boldsymbol{\Psi}_t(\boldsymbol\alpha, \boldsymbol\omega)</math> a component of the vector <math>|\Psi\rangle</math>, called the wave function of the system
  • dimensionless discrete quantum numbers
  • continuous variables (not necessarily dimensionless)

These quantum numbers index the components of the state vector. More, all are in an -dimensional set where each is the set of allowed values for ; all are in an -dimensional "volume" where and each is the set of allowed values for , a subset of the real numbers . For generality and are not necessarily equal.

Example:

The probability density of finding the system at time <math>t</math> at state is

<math display="block">\rho_{\alpha, \omega} (t)= |\Psi(\boldsymbol{\alpha},\boldsymbol{\omega},t)|^2</math>

The probability of finding system with in some or all possible discrete-variable configurations, , and in some or all possible continuous-variable configurations, , is the sum and integral over the density,

<math display="block">P(t)=\sum_{\boldsymbol{\alpha}\in D}\int_C d^m\!\boldsymbol{\omega}\,\,\rho_{\alpha, \omega}(t)</math>

Since the sum of all probabilities must be 1, the normalization condition

<math display="block">1=\sum_{\boldsymbol{\alpha}\in A}\int_{\Omega}d^m\!\boldsymbol{\omega}\,\,\rho_{\alpha, \omega} (t)</math>

must hold at all times during the evolution of the system.

The normalization condition requires to be dimensionless, by dimensional analysis must have the same units as .

Ontology

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 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.

See also

  • Boson
  • De Broglie–Bohm theory
  • Double-slit experiment
  • Faraday wave
  • Fermion
  • Phase-space formulation
  • Schrödinger equation
  • Wave function collapse
  • Wave packet

Notes

Remarks

Citations

References

  • Online copy (French) Online copy (English)
  • Online copy

Further reading

  • 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
  • Einstein, The quantum theory of radiation