Probability amplitude

Mathematical formulation

In a formal setup, the state of an isolated physical system in quantum mechanics is represented, at a fixed time ${\displaystyle t}$, by a state vector |Ψ⟩ belonging to a separable complex Hilbert space. Using bra–ket notation the relation between state vector and "position basis" ${\displaystyle \{|x⟩\}}$ of the Hilbert space can be written as

${\displaystyle \psi (x)=⟨x|\mathrm{\Psi }⟩}$.

Its relation with an observable can be elucidated by generalizing the quantum state ${\displaystyle \psi }$ to a measurable function and its domain of definition to a given σ-finite measure space ${\displaystyle (X,\mathcal{A},\mu )}$. This allows for a refinement of Lebesgue's decomposition theorem, decomposing μ into three mutually singular parts

${\displaystyle \mu ={\mu }_{\mathrm{a}\mathrm{c}}+{\mu }_{\mathrm{s}\mathrm{c}}+{\mu }_{\mathrm{p}\mathrm{p}}}$

where μ_ac is absolutely continuous with respect to the Lebesgue measure, μ_sc is singular with respect to the Lebesgue measure and atomless, and μ_pp is a pure point measure.

Continuous amplitudes

A usual presentation of the probability amplitude is that of a wave function ${\displaystyle \psi }$ belonging to the L² space of (equivalence classes of) square integrable functions, i.e., ${\displaystyle \psi }$ belongs to L²(X) if and only if

${\displaystyle \Vert \psi {\Vert }^2={\int }_X|\psi (x){|}^{2}dx<\mathrm{\infty }}$.

If the norm is equal to 1 and ${\displaystyle |\psi (x){|}^2\in {\mathbb{R}}_{\ge 0}}$ such that

${\displaystyle {\int }_X|\psi (x){|}^{2}dx\equiv {\int }_{X}d{\mu }_{ac}(x)=1}$,

then ${\displaystyle |\psi (x){|}^2}$ is the probability density function for a measurement of the particle's position at a given time, defined as the Radon–Nikodym derivative with respect to the Lebesgue measure (e.g. on the set R of all real numbers). As probability is a dimensionless quantity, |ψ(x)|² must have the inverse dimension of the variable of integration x. For example, the above amplitude has dimension [L^−1/2], where L represents length.

Whereas a Hilbert space is separable if and only if it admits a countable orthonormal basis, the range of a continuous random variable ${\displaystyle x}$ is an uncountable set (i.e. the probability that the system is "at position ${\displaystyle x}$" will always be zero). As such, eigenstates of an observable need not necessarily be measurable functions belonging to L²(X) (see normalization condition below). A typical example is the position operator ${\displaystyle \hat{\mathrm{x}}}$ defined as

${\displaystyle ⟨x|\hat{\mathrm{x}}|\mathrm{\Psi }⟩=\hat{\mathrm{x}}⟨x|\mathrm{\Psi }⟩=x_0\psi (x),x\in \mathbb{R},}$

whose eigenfunctions are Dirac delta functions

${\displaystyle \psi (x)=\delta (x-x_0)}$

which clearly do not belong to L²(X). By replacing the state space by a suitable rigged Hilbert space, however, the rigorous notion of eigenstates from spectral theorem as well as spectral decomposition is preserved.

Discrete amplitudes

Let ${\displaystyle {\mu }_{pp}}$ be atomic (i.e. the set ${\displaystyle A\subset X}$ in ${\displaystyle \mathcal{A}}$ is an atom); specifying the measure of any discrete variable x ∈ A equal to 1. The amplitudes are composed of state vector |Ψ⟩ indexed by A; its components are denoted by ψ(x) for uniformity with the previous case. If the ℓ²-norm of |Ψ⟩ is equal to 1, then |ψ(x)|² is a probability mass function.

A convenient configuration space X is such that each point x produces some unique value of the observable Q. For discrete X it means that all elements of the standard basis are eigenvectors of Q. Then ${\displaystyle \psi (x)}$ is the probability amplitude for the eigenstate |x⟩. If it corresponds to a non-degenerate eigenvalue of Q, then ${\displaystyle |\psi (x){|}^2}$ gives the probability of the corresponding value of Q for the initial state |Ψ⟩.

|ψ(x)| = 1 if and only if |x⟩ is the same quantum state as |Ψ⟩. ψ(x) = 0 if and only if |x⟩ and |Ψ⟩ are orthogonal. Otherwise the modulus of ψ(x) is between 0 and 1.

A discrete probability amplitude may be considered as a fundamental frequency in the probability frequency domain (spherical harmonics) for the purposes of simplifying M-theory transformation calculations. Discrete dynamical variables are used in such problems as a particle in an idealized reflective box and quantum harmonic oscillator.

Examples

An example of the discrete case is a quantum system that can be in two possible states, e.g. the polarization of a photon. When the polarization is measured, it could be the horizontal state ${\displaystyle |H⟩}$ or the vertical state ${\displaystyle |V⟩}$. Until its polarization is measured the photon can be in a superposition of both these states, so its state ${\displaystyle |\psi ⟩}$ could be written as

${\displaystyle |\psi ⟩=\alpha |H⟩+\beta |V⟩}$,

with ${\displaystyle \alpha }$ and ${\displaystyle \beta }$ the probability amplitudes for the states ${\displaystyle |H⟩}$ and ${\displaystyle |V⟩}$ respectively. When the photon's polarization is measured, the resulting state is either horizontal or vertical. But in a random experiment, the probability of being horizontally polarized is ${\displaystyle |\alpha {|}^2}$, and the probability of being vertically polarized is ${\displaystyle |\beta {|}^2}$.

Hence, a photon in a state $|\psi ⟩=\sqrt{\frac{1}{3}}|H⟩-i\sqrt{\frac{2}{3}}|V⟩$ would have a probability of $\frac{1}{3}$ to come out horizontally polarized, and a probability of $\frac{2}{3}$ to come out vertically polarized when an ensemble of measurements are made. The order of such results, is, however, completely random.

Another example is quantum spin. If a spin-measuring apparatus is pointing along the z-axis and is therefore able to measure the z-component of the spin (${\sigma }_z$), the following must be true for the measurement of spin "up" and "down":

${\displaystyle {\sigma }_z|u⟩=(+1)|u⟩}$

${\displaystyle {\sigma }_z|d⟩=(-1)|d⟩}$

If one assumes that system is prepared, so that +1 is registered in ${\sigma }_x$ and then the apparatus is rotated to measure ${\sigma }_z$, the following holds:

${\displaystyle \begin{array}{cc}⟨r|u⟩& =\left(\frac{1}{\sqrt{2}}⟨u|+\frac{1}{\sqrt{2}}⟨d|\right)\cdot |u⟩\\ & =\left(\frac{1}{\sqrt{2}}\left(\begin{array}{c}1\\ 0\end{array}\right)+\frac{1}{\sqrt{2}}\left(\begin{array}{c}0\\ 1\end{array}\right)\right)\cdot \left(\begin{array}{c}1\\ 0\end{array}\right)\\ & =\frac{1}{\sqrt{2}}\end{array}}$

The probability amplitude of measuring spin up is given by $⟨r|u⟩$, since the system had the initial state $|r⟩$. The probability of measuring $|u⟩$ is given by

${\displaystyle P(|u⟩)=⟨r|u⟩⟨u|r⟩={\left(\frac{1}{\sqrt{2}}\right)}^2=\frac{1}{2}}$

Which agrees with experiment.

Normalization

In the example above, the measurement must give either | H ⟩ or | V ⟩, so the total probability of measuring | H ⟩ or | V ⟩ must be 1. This leads to a constraint that α² + β² = 1; more generally the sum of the squared moduli of the probability amplitudes of all the possible states is equal to one. If to understand "all the possible states" as an orthonormal basis, that makes sense in the discrete case, then this condition is the same as the norm-1 condition explained above.

One can always divide any non-zero element of a Hilbert space by its norm and obtain a normalized state vector. Not every wave function belongs to the Hilbert space L²(X), though. Wave functions that fulfill this constraint are called normalizable.

The Schrödinger equation, describing states of quantum particles, has solutions that describe a system and determine precisely how the state changes with time. Suppose a wave function ψ(x, t) gives a description of the particle (position x at a given time t). A wave function is square integrable if

${\displaystyle \int |\psi (\mathbf{x},t){|}^2\mathrm{d}\mathbf{x}=a^2<\mathrm{\infty }.}$

After normalization the wave function still represents the same state and is therefore equal by definition to

${\displaystyle \psi (\mathbf{x},t):=\frac{\psi (\mathbf{x},t)}{a}.}$

Under the standard Copenhagen interpretation, the normalized wavefunction gives probability amplitudes for the position of the particle. Hence, ρ(x) = |ψ(x, t)|² is a probability density function and the probability that the particle is in the volume V at fixed time t is given by

${\displaystyle P_{\mathbf{x}\in V}(t)={\int }_V|\psi (\mathbf{x},t){|}^2\mathrm{d}\mathbf{x}={\int }_V\rho (\mathbf{x})\mathrm{d}\mathbf{x}.}$

The probability density function does not vary with time as the evolution of the wave function is dictated by the Schrödinger equation and is therefore entirely deterministic. This is key to understanding the importance of this interpretation: for a given particle constant mass, initial ψ(x, t₀) and potential, the Schrödinger equation fully determines subsequent wavefunctions. The above then gives probabilities of locations of the particle at all subsequent times.

In the context of the double-slit experiment

Probability amplitudes have special significance because they act in quantum mechanics as the equivalent of conventional probabilities, with many analogous laws, as described above. For example, in the classic double-slit experiment, electrons are fired randomly at two slits, and the probability distribution of detecting electrons at all parts on a large screen placed behind the slits, is questioned. An intuitive answer is that P(through either slit) = P(through first slit) + P(through second slit), where P(event) is the probability of that event. This is obvious if one assumes that an electron passes through either slit. When no measurement apparatus that determines through which slit the electrons travel is installed, the observed probability distribution on the screen reflects the interference pattern that is common with light waves. If one assumes the above law to be true, then this pattern cannot be explained. The particles cannot be said to go through either slit and the simple explanation does not work. The correct explanation is, however, by the association of probability amplitudes to each event. The complex amplitudes which represent the electron passing each slit (ψ_first and ψ_second) follow the law of precisely the form expected: ψ_total = ψ_first + ψ_second. This is the principle of quantum superposition. The probability, which is the modulus squared of the probability amplitude, then, follows the interference pattern under the requirement that amplitudes are complex: $${\displaystyle P={\left|{\psi }_{\text{first}}+{\psi }_{\text{second}}\right|}^2={\left|{\psi }_{\text{first}}\right|}^2+{\left|{\psi }_{\text{second}}\right|}^2+2\left|{\psi }_{\text{first}}\right|\left|{\psi }_{\text{second}}\right|\cos ({\phi }_1-{\phi }_2).}$$ Here, ${\displaystyle {\phi }_1}$ and ${\displaystyle {\phi }_2}$ are the arguments of ψ_first and ψ_second respectively. A purely real formulation has too few dimensions to describe the system's state when superposition is taken into account. That is, without the arguments of the amplitudes, we cannot describe the phase-dependent interference. The crucial term $2\left|{\psi }_{\text{first}}\right|\left|{\psi }_{\text{second}}\right|\cos ({\phi }_1-{\phi }_2)$ is called the "interference term", and this would be missing if we had added the probabilities.

However, one may choose to devise an experiment in which the experimenter observes which slit each electron goes through. Then, due to wavefunction collapse, the interference pattern is not observed on the screen.

One may go further in devising an experiment in which the experimenter gets rid of this "which-path information" by a "quantum eraser". Then, according to the Copenhagen interpretation, the case A applies again and the interference pattern is restored.

Conservation of probabilities and the continuity equation

Intuitively, since a normalised wave function stays normalised while evolving according to the wave equation, there will be a relationship between the change in the probability density of the particle's position and the change in the amplitude at these positions.

Define the probability current (or flux) j as

${\displaystyle \mathbf{j}=\frac{\hslash }{m}\frac{1}{2i}\left({\psi }^{\ast }\mathrm{\nabla }\psi -\psi \mathrm{\nabla }{\psi }^{\ast }\right)=\frac{\hslash }{m}\mathrm{Im}\left({\psi }^{\ast }\mathrm{\nabla }\psi \right),}$

measured in units of (probability)/(area × time).

Then the current satisfies the equation

${\displaystyle \mathrm{\nabla }\cdot \mathbf{j}+\frac{\mathrm{\partial }}{\mathrm{\partial }t}|\psi {|}^2=0.}$

The probability density is ${\displaystyle \rho =|\psi {|}^2}$, this equation is exactly the continuity equation, appearing in many situations in physics where we need to describe the local conservation of quantities. The best example is in classical electrodynamics, where j corresponds to current density corresponding to electric charge, and the density is the charge-density. The corresponding continuity equation describes the local conservation of charges.

Composite systems

For two quantum systems with spaces L²(X₁) and L²(X₂) and given states |Ψ₁⟩ and |Ψ₂⟩ respectively, their combined state |Ψ₁⟩ ⊗ |Ψ₂⟩ can be expressed as ψ₁(x₁) ψ₂(x₂) a function on X₁ × X₂, that gives the product of respective probability measures. In other words, amplitudes of a non-entangled composite state are products of original amplitudes, and respective observables on the systems 1 and 2 behave on these states as independent random variables. This strengthens the probabilistic interpretation explicated above.

Amplitudes in operators

The concept of amplitudes is also used in the context of scattering theory, notably in the form of S-matrices. Whereas moduli of vector components squared, for a given vector, give a fixed probability distribution, moduli of matrix elements squared are interpreted as transition probabilities just as in a random process. Like a finite-dimensional unit vector specifies a finite probability distribution, a finite-dimensional unitary matrix specifies transition probabilities between a finite number of states.

The "transitional" interpretation may be applied to L²s on non-discrete spaces as well.