for illustrative purposes only; in quantum mechanics the wave function is generally
complex. It is vital to illustrate how the principle applies to relatively intelligible physical situations since it is indiscernible on the macroscopic scales that humans experience. Two alternative frameworks for quantum physics offer different explanations for the uncertainty principle. The
wave mechanics picture of the uncertainty principle is more visually intuitive, but the more abstract
matrix mechanics picture formulates it in a way that generalizes more easily. Mathematically, in wave mechanics, the uncertainty relation between position and momentum arises because the expressions of the wavefunction in the two corresponding
orthonormal bases in
Hilbert space are
Fourier transforms of one another (i.e., position and momentum are
conjugate variables). A nonzero function and its Fourier transform cannot both be sharply localized at the same time. A similar tradeoff between the variances of Fourier conjugates arises in all systems underlain by Fourier analysis, for example in sound waves: A pure tone is a
sharp spike at a single frequency, while its Fourier transform gives the shape of the sound wave in the time domain, which is a completely delocalized sine wave. In quantum mechanics, the two key points are that the position of the particle takes the form of a matter wave, and momentum is its Fourier conjugate, assured by the
de Broglie relation , where is the
wavenumber. In
matrix mechanics, the
mathematical formulation of quantum mechanics, any pair of non-
commuting self-adjoint operators representing
observables are subject to similar uncertainty limits. An eigenstate of an observable represents the state of the wavefunction for a certain measurement value (the eigenvalue). For example, if a measurement of an observable is performed, then the system is in a particular eigenstate of that observable. However, the particular eigenstate of the observable need not be an eigenstate of another observable : If so, then it does not have a unique associated measurement for it, as the system is not in an eigenstate of that observable.
Visualization The uncertainty principle can be visualized using the position- and momentum-space wavefunctions for one spinless
particle with mass in one dimension. The more localized the position-space wavefunction, the more likely the particle is to be found with the position coordinates in that region, and correspondingly the momentum-space wavefunction is less localized so the possible momentum components the particle could have are more widespread. Conversely, the more localized the momentum-space wavefunction, the more likely the particle is to be found with those values of momentum components in that region, and correspondingly the less localized the position-space wavefunction, so the position coordinates the particle could occupy are more widespread. These wavefunctions are
Fourier transforms of each other: mathematically, the uncertainty principle expresses the relationship between conjugate variables in the transform. of finding the particle with position
x or momentum component
p.
Top: If wavelength
λ is unknown, so are momentum
p, wave-vector
k and energy
E (de Broglie relations). As the particle is more localized in position space, Δ
x is smaller than for Δ
px.
Bottom: If
λ is known, so are
p,
k, and
E. As the particle is more localized in momentum space, Δ
p is smaller than for Δ
x.
Wave mechanics interpretation According to the
de Broglie hypothesis, every object in the universe is associated with a
wave. Thus every object, from an elementary particle to atoms, molecules and on up to planets and beyond are subject to the uncertainty principle. The time-independent wave function of a single-moded plane wave of wavenumber or momentum is \psi(x) \propto e^{ik_0 x} = e^{ip_0 x/\hbar} ~. The
Born rule states that this should be interpreted as a
probability density amplitude function in the sense that the probability of finding the particle between and is \operatorname P [a \leq X \leq b] = \int_a^b |\psi(x)|^2 \, \mathrm{d}x ~. In the case of the single-mode plane wave, |\psi(x)|^2 is 1 if X=x and 0 otherwise. In other words, the particle position is extremely uncertain in the sense that it could be essentially anywhere along the wave packet. On the other hand, consider a wave function that is a
sum of many waves, which we may write as \psi(x) \propto \sum_n A_n e^{i p_n x/\hbar}~, where represents the relative contribution of the mode to the overall total. The figures to the right show how with the addition of many plane waves, the wave packet can become more localized. We may take this a step further to the
continuum limit, where the wave function is an
integral over all possible modes \psi(x) = \frac{1}{\sqrt{2 \pi \hbar}} \int_{-\infty}^\infty \varphi(p) \cdot e^{i p x/\hbar} \, dp ~, with \varphi(p) representing the amplitude of these modes and is called the wave function in
momentum space. In mathematical terms, we say that \varphi(p) is the
Fourier transform of \psi(x) and that and are
conjugate variables. Adding together all of these plane waves comes at a cost, namely the momentum has become less precise, having become a mixture of waves of many different momenta. One way to quantify the precision of the position and momentum is the
standard deviation . Since |\psi(x)|^2 is a probability density function for position, we calculate its standard deviation. The precision of the position is improved, i.e. reduced , by using many plane waves, thereby weakening the precision of the momentum, i.e. increased . Another way of stating this is that and have an
inverse relationship or are at least bounded from below. This is the uncertainty principle, the exact limit of which is the Kennard bound.
Proof of the Kennard inequality using wave mechanics We are interested in the
variances of position and momentum, defined as \sigma_x^2 = \int_{-\infty}^\infty x^2 \cdot |\psi(x)|^2 \, dx - \left( \int_{-\infty}^\infty x \cdot |\psi(x)|^2 \, dx \right)^2 \sigma_p^2 = \int_{-\infty}^\infty p^2 \cdot |\varphi(p)|^2 \, dp - \left( \int_{-\infty}^\infty p \cdot |\varphi(p)|^2 \, dp \right)^2~.
Without loss of generality, we will assume that the
means vanish, which just amounts to a shift of the origin of our coordinates. (A more general proof that does not make this assumption is given below.) This gives us the simpler form \sigma_x^2 = \int_{-\infty}^\infty x^2 \cdot |\psi(x)|^2 \, dx \sigma_p^2 = \int_{-\infty}^\infty p^2 \cdot |\varphi(p)|^2 \, dp~. The function f(x) = x \cdot \psi(x) can be interpreted as a
vector in a
function space. We can define an
inner product for a pair of functions
u(
x) and
v(
x) in this vector space: \langle u \mid v \rangle = \int_{-\infty}^\infty u^*(x) \cdot v(x) \, dx, where the asterisk denotes the
complex conjugate. With this inner product defined, we note that the variance for position can be written as \sigma_x^2 = \int_{-\infty}^\infty |f(x)|^2 \, dx = \langle f \mid f \rangle ~. We can repeat this for momentum by interpreting the function \tilde{g}(p)=p \cdot \varphi(p) as a vector, but we can also take advantage of the fact that \psi(x) and \varphi(p) are Fourier transforms of each other. We evaluate the inverse Fourier transform through
integration by parts: \begin{align} g(x) &= \frac{1}{\sqrt{2 \pi \hbar}} \cdot \int_{-\infty}^\infty \tilde{g}(p) \cdot e^{ipx/\hbar} \, dp \\ &= \frac{1}{\sqrt{2 \pi \hbar}} \int_{-\infty}^\infty p \cdot \varphi(p) \cdot e^{ipx/\hbar} \, dp \\ &= \frac{1}{2 \pi \hbar} \int_{-\infty}^\infty \left[ p \cdot \int_{-\infty}^\infty \psi(\chi) e^{-ip\chi/\hbar} \, d\chi \right] \cdot e^{ipx/\hbar} \, dp \\ &= \frac{i}{2 \pi} \int_{-\infty}^\infty \left[ \cancel{ \left. \psi(\chi) e^{-ip\chi/\hbar} \right|_{-\infty}^\infty } - \int_{-\infty}^\infty \frac{d\psi(\chi)}{d\chi} e^{-ip\chi/\hbar} \, d\chi \right] \cdot e^{ipx/\hbar} \, dp \\ &= -i \int_{-\infty}^\infty \frac{d\psi(\chi)}{d\chi} \left[ \frac{1}{2 \pi}\int_{-\infty}^\infty \, e^{ip(x - \chi)/\hbar} \, dp \right]\, d\chi\\ &= -i \int_{-\infty}^\infty \frac{d\psi(\chi)}{d\chi} \left[ \delta\left(\frac{x - \chi }{\hbar}\right) \right]\, d\chi\\ &= -i \hbar \int_{-\infty}^\infty \frac{d\psi(\chi)}{d\chi} \left[ \delta\left(x - \chi \right) \right]\, d\chi\\ &= -i \hbar \frac{d\psi(x)}{dx} \\ &= \left( -i \hbar \frac{d}{dx} \right) \cdot \psi(x) , \end{align} where v=\frac{\hbar}{-ip}e^{-ip\chi/\hbar} in the integration by parts, the cancelled term vanishes because the wave function vanishes at both infinities and |e^{-ip\chi/\hbar}|=1, and then use the
Dirac delta function which is valid because \dfrac{d\psi(\chi)}{d\chi} does not depend on
p . The term -i \hbar \frac{d}{dx} is called the
momentum operator in position space. Applying
Plancherel's theorem, we see that the variance for momentum can be written as \sigma_p^2 = \int_{-\infty}^\infty |\tilde{g}(p)|^2 \, dp = \int_{-\infty}^\infty |g(x)|^2 \, dx = \langle g \mid g \rangle. The
Cauchy–Schwarz inequality asserts that \sigma_x^2 \sigma_p^2 = \langle f \mid f \rangle \cdot \langle g \mid g \rangle \ge |\langle f \mid g \rangle|^2 ~. The
modulus squared of any complex number
z can be expressed as |z|^{2} = \Big(\text{Re}(z)\Big)^{2}+\Big(\text{Im}(z)\Big)^{2} \geq \Big(\text{Im}(z)\Big)^{2} = \left(\frac{z-z^{\ast}}{2i}\right)^{2}. we let z=\langle f|g\rangle and z^{*}=\langle g\mid f\rangle and substitute these into the equation above to get |\langle f\mid g\rangle|^2 \geq \left(\frac{\langle f\mid g\rangle-\langle g \mid f \rangle}{2i}\right)^2 ~. All that remains is to evaluate these inner products. \begin{align} \langle f\mid g\rangle-\langle g\mid f\rangle &= \int_{-\infty}^\infty \psi^*(x) \, x \cdot \left(-i \hbar \frac{d}{dx}\right) \, \psi(x) \, dx - \int_{-\infty}^\infty \psi^*(x) \, \left(-i \hbar \frac{d}{dx}\right) \cdot x \, \psi(x) \, dx \\ &= i \hbar \cdot \int_{-\infty}^\infty \psi^*(x) \left[ \left(-x \cdot \frac{d\psi(x)}{dx}\right) + \frac{d(x \psi(x))}{dx} \right] \, dx \\ &= i \hbar \cdot \int_{-\infty}^\infty \psi^*(x) \left[ \left(-x \cdot \frac{d\psi(x)}{dx}\right) + \psi(x) + \left(x \cdot \frac{d\psi(x)}{dx}\right)\right] \, dx \\ &= i \hbar \cdot \int_{-\infty}^\infty \psi^*(x) \psi(x) \, dx \\ &= i \hbar \cdot \int_{-\infty}^\infty |\psi(x)|^2 \, dx \\ &= i \hbar \end{align} Plugging this into the above inequalities, we get \sigma_x^2 \sigma_p^2 \ge |\langle f \mid g \rangle|^2 \ge \left(\frac{\langle f\mid g\rangle-\langle g\mid f\rangle}{2i}\right)^2 = \left(\frac{i \hbar}{2 i}\right)^2 = \frac{\hbar^2}{4} and taking the square root \sigma_x \sigma_p \ge \frac{\hbar}{2}~. with equality if and only if
p and
x are linearly dependent. Note that the only
physics involved in this proof was that \psi(x) and \varphi(p) are wave functions for position and momentum, which are Fourier transforms of each other. A similar result would hold for
any pair of conjugate variables.
Matrix mechanics interpretation In matrix mechanics, observables such as position and momentum are represented by self-adjoint operators. When considering pairs of observables, an important quantity is the
commutator. For a pair of operators and \hat{B}, one defines their commutator as [\hat{A},\hat{B}]=\hat{A}\hat{B}-\hat{B}\hat{A}. In the case of position and momentum, the commutator is the
canonical commutation relation [\hat{x},\hat{p}]=i \hbar. The physical meaning of the non-commutativity can be understood by considering the effect of the commutator on position and momentum
eigenstates. Let |\psi\rangle be a right eigenstate of position with a constant eigenvalue . By definition, this means that \hat{x}|\psi\rangle = x_0 |\psi\rangle. Applying the commutator to |\psi\rangle yields [\hat{x},\hat{p}] | \psi \rangle = (\hat{x}\hat{p}-\hat{p}\hat{x}) | \psi \rangle = (\hat{x} - x_0 \hat{I}) \hat{p} \, | \psi \rangle = i \hbar | \psi \rangle, where is the
identity operator. Suppose, for the sake of
proof by contradiction, that |\psi\rangle is also a right eigenstate of momentum, with constant eigenvalue . If this were true, then one could write (\hat{x} - x_0 \hat{I}) \hat{p} \, | \psi \rangle = (\hat{x} - x_0 \hat{I}) p_0 \, | \psi \rangle = (x_0 \hat{I} - x_0 \hat{I}) p_0 \, | \psi \rangle=0. On the other hand, the above canonical commutation relation requires that [\hat{x},\hat{p}] | \psi \rangle=i \hbar | \psi \rangle \ne 0. This implies that no quantum state can simultaneously be both a position and a momentum eigenstate. When a state is measured, it is projected onto an eigenstate in the basis of the relevant observable. For example, if a particle's position is measured, then the state amounts to a position eigenstate. This means that the state is
not a momentum eigenstate, however, but rather it can be represented as a sum of multiple momentum basis eigenstates. In other words, the momentum must be less precise. This precision may be quantified by the standard deviations, \sigma_x=\sqrt{\langle \hat{x}^2 \rangle-\langle \hat{x}\rangle^2} \sigma_p=\sqrt{\langle \hat{p}^2 \rangle-\langle \hat{p}\rangle^2}. As in the wave mechanics interpretation above, one sees a tradeoff between the respective precisions of the two, quantified by the uncertainty principle.
Quantum harmonic oscillator stationary states Consider a one-dimensional
quantum harmonic oscillator. It is possible to express the position and momentum operators in terms of the
creation and annihilation operators: \hat x = \sqrt{\frac{\hbar}{2m\omega}}(a+a^\dagger) \hat p = i\sqrt{\frac{m \omega\hbar}{2}}(a^\dagger-a). Using the standard rules for creation and annihilation operators on the energy eigenstates, a^{\dagger}|n\rangle=\sqrt{n+1}|n+1\rangle a|n\rangle=\sqrt{n}|n-1\rangle, the variances may be computed directly, \sigma_x^2 = \frac{\hbar}{m\omega} \left( n+\frac{1}{2}\right) \sigma_p^2 = \hbar m\omega \left( n+\frac{1}{2}\right)\, . The product of these standard deviations is then \sigma_x \sigma_p = \hbar \left(n+\frac{1}{2}\right) \ge \frac{\hbar}{2}.~ In particular, the above Kennard bound is saturated for the
ground state , for which the probability density is just the
normal distribution.
Quantum harmonic oscillators with Gaussian initial condition In a quantum harmonic oscillator of characteristic angular frequency
ω, place a state that is offset from the bottom of the potential by some displacement
x0 as \psi(x)=\left(\frac{m \Omega}{\pi \hbar}\right)^{1/4} \exp{\left( -\frac{m \Omega (x-x_0)^2}{2\hbar}\right)}, where Ω describes the width of the initial state but need not be the same as
ω. Through integration over the
propagator, we can solve for the -dependent solution. After many cancelations, the probability densities reduce to |\Psi(x,t)|^2 \sim \mathcal{N}\left( x_0 \cos{(\omega t)} , \frac{\hbar}{2 m \Omega} \left( \cos^2(\omega t) + \frac{\Omega^2}{\omega^2} \sin^2{(\omega t)} \right)\right) |\Phi(p,t)|^2 \sim \mathcal{N}\left( -m x_0 \omega \sin(\omega t), \frac{\hbar m \Omega}{2} \left( \cos^2{(\omega t)} + \frac{\omega^2}{\Omega^2} \sin^2{(\omega t)} \right)\right), where we have used the notation \mathcal{N}(\mu, \sigma^2) to denote a normal distribution of mean
μ and variance
σ2. Copying the variances above and applying
trigonometric identities, we can write the product of the standard deviations as \begin{align} \sigma_x \sigma_p&=\frac{\hbar}{2}\sqrt{\left( \cos^2{(\omega t)} + \frac{\Omega^2}{\omega^2} \sin^2{(\omega t)} \right)\left( \cos^2{(\omega t)} + \frac{\omega^2}{\Omega^2} \sin^2{(\omega t)} \right)} \\ &= \frac{\hbar}{4}\sqrt{3+\frac{1}{2}\left(\frac{\Omega^2}{\omega^2}+\frac{\omega^2}{\Omega^2}\right)-\left(\frac{1}{2}\left(\frac{\Omega^2}{\omega^2}+\frac{\omega^2}{\Omega^2}\right)-1\right) \cos{(4 \omega t)}} \end{align} From the relations \frac{\Omega^2}{\omega^2}+\frac{\omega^2}{\Omega^2} \ge 2, \quad |\cos(4 \omega t)| \le 1, we can conclude the following (the right most equality holds only when ): \sigma_x \sigma_p \ge \frac{\hbar}{4}\sqrt{3+\frac{1}{2} \left(\frac{\Omega^2}{\omega^2}+\frac{\omega^2}{\Omega^2}\right)-\left(\frac{1}{2} \left(\frac{\Omega^2}{\omega^2}+\frac{\omega^2}{\Omega^2}\right)-1\right)} = \frac{\hbar}{2}.
Coherent states A coherent state is a right eigenstate of the
annihilation operator, \hat{a}|\alpha\rangle=\alpha|\alpha\rangle, which may be represented in terms of
Fock states as |\alpha\rangle =e^{-{|\alpha|^2\over2}} \sum_{n=0}^\infty {\alpha^n \over \sqrt{n!}}|n\rangle In the picture where the coherent state is a massive particle in a quantum harmonic oscillator, the position and momentum operators may be expressed in terms of the annihilation operators in the same formulas above and used to calculate the variances, \sigma_x^2 = \frac{\hbar}{2 m \omega}, \sigma_p^2 = \frac{\hbar m \omega}{2}. Therefore, every coherent state saturates the Kennard bound \sigma_x \sigma_p = \sqrt{\frac{\hbar}{2 m \omega}} \, \sqrt{\frac{\hbar m \omega}{2}} = \frac{\hbar}{2}. with position and momentum each contributing an amount \sqrt{\hbar/2} in a "balanced" way. Moreover, every
squeezed coherent state also saturates the Kennard bound although the individual contributions of position and momentum need not be balanced in general.
Particle in a box Consider a particle in a one-dimensional box of length L. The
eigenfunctions in position and momentum space are \psi_n(x,t) =\begin{cases} A \sin(k_n x)\mathrm{e}^{-\mathrm{i}\omega_n t}, & 0 and \varphi_n(p,t)=\sqrt{\frac{\pi L}{\hbar}}\,\,\frac{n\left(1-(-1)^ne^{-ikL} \right) e^{-i \omega_n t}}{\pi ^2 n^2-k^2 L^2}, where \omega_n=\frac{\pi^2 \hbar n^2}{8 L^2 m} and we have used the
de Broglie relation p=\hbar k. The variances of x and p can be calculated explicitly: \sigma_x^2=\frac{L^2}{12}\left(1-\frac{6}{n^2\pi^2}\right) \sigma_p^2=\left(\frac{\hbar n\pi}{L}\right)^2. The product of the standard deviations is therefore \sigma_x \sigma_p = \frac{\hbar}{2} \sqrt{\frac{n^2\pi^2}{3}-2}. For all n=1, \, 2, \, 3,\, \ldots, the quantity \sqrt{\frac{n^2\pi^2}{3}-2} is greater than 1, so the uncertainty principle is never violated. For numerical concreteness, the smallest value occurs when n = 1, in which case \sigma_x \sigma_p = \frac{\hbar}{2} \sqrt{\frac{\pi^2}{3}-2} \approx 0.568 \hbar > \frac{\hbar}{2}.
Constant momentum Assume a particle initially has a
momentum space wave function described by a normal distribution around some constant momentum
p0 according to \varphi(p) = \left(\frac{x_0}{\hbar \sqrt{\pi}} \right)^{1/2} \exp\left(\frac{-x_0^2 (p-p_0)^2}{2\hbar^2}\right), where we have introduced a reference scale x_0=\sqrt{\hbar/m\omega_0}, with \omega_0>0 describing the width of the distribution—cf.
nondimensionalization. If the state is allowed to evolve in free space, then the time-dependent momentum and position space wave functions are \Phi(p,t) = \left(\frac{x_0}{\hbar \sqrt{\pi}} \right)^{1/2} \exp\left(\frac{-x_0^2 (p-p_0)^2}{2\hbar^2}-\frac{ip^2 t}{2m\hbar}\right), \Psi(x,t) = \left(\frac{1}{x_0 \sqrt{\pi}} \right)^{1/2} \frac{e^{-x_0^2 p_0^2 /2\hbar^2}}{\sqrt{1+i\omega_0 t}} \, \exp\left(-\frac{(x-ix_0^2 p_0/\hbar)^2}{2x_0^2 (1+i\omega_0 t)}\right). Since \langle p(t) \rangle = p_0 and \sigma_p(t) = \hbar /(\sqrt{2}x_0), this can be interpreted as a particle moving along with constant momentum at arbitrarily high precision. On the other hand, the standard deviation of the position is \sigma_x = \frac{x_0}{\sqrt{2}} \sqrt{1+\omega_0^2 t^2} such that the uncertainty product can only increase with time as \sigma_x(t) \sigma_p(t) = \frac{\hbar}{2} \sqrt{1+\omega_0^2 t^2} ==Mathematical formalism==