# Momentum operator

In quantum mechanics, momentum (like all other physical variables) is defined as an operator, which "acts on" or pre-multiplies the wave function ψ(r, t) to extract the momentum eigenvalue from the wave function: the momentum vector a particle would have when measured in an experiment. The momentum operator is an example of a differential operator.

At the time quantum mechanics was developed in the 1920s, the momentum operator was found by many theoretical physicists, including Niels Bohr, Arnold Sommerfeld, Erwin Schrödinger, and Eugene Wigner.

## Origin from De Broglie plane waves

The momentum and energy operators can be constructed in the following way.[1]

### One dimension

Starting in one dimension, using the plane wave solution to Schrödinger's equation:

${\displaystyle \psi =e^{i(kx-\omega t)}}$

The first order partial derivative with respect to space is

${\displaystyle {\frac {\partial \psi }{\partial x}}=ike^{i(kx-\omega t)}=ik\psi }$

By expressing Template:Mvar from the De Broglie relation:

${\displaystyle p=\hbar k}$

the formula for the derivative of Template:Mvar becomes:

${\displaystyle {\frac {\partial \psi }{\partial x}}=i{\frac {p}{\hbar }}\psi }$

This suggests the operator equivalence:

${\displaystyle {\hat {p}}=-i\hbar {\frac {\partial }{\partial x}}}$

so the momentum value Template:Mvar is a scalar factor, the momentum of the particle and the value that is measured, is the eigenvalue of the operator.

Since the partial derivative is a linear operator, the momentum operator is also linear, and because any wavefunction can be expressed as a superposition of other states, when this momentum operator acts on the entire superimposed wave, it yields the momentum eigenvalues for each plane wave component, the momenta add to the total momentum of the superimposed wave.

### Three dimensions

The derivation in three dimensions is the same, except the gradient operator del is used instead of one partial derivative. In three dimensions, the plane wave solution to Schrödinger's equation is:

${\displaystyle \psi =e^{i({\mathbf {k}}\cdot {\mathbf {r}}-\omega t)}}$

{\displaystyle {\begin{aligned}\nabla \psi &={\mathbf {e}}_{x}{\frac {\partial \psi }{\partial x}}+{\mathbf {e}}_{y}{\frac {\partial \psi }{\partial y}}+{\mathbf {e}}_{z}{\frac {\partial \psi }{\partial z}}\\&=ik_{x}\psi {\mathbf {e}}_{x}+ik_{y}\psi {\mathbf {e}}_{y}+ik_{z}\psi {\mathbf {e}}_{z}\\&={\frac {i}{\hbar }}\left(p_{x}{\mathbf {e}}_{x}+p_{y}{\mathbf {e}}_{y}+p_{z}{\mathbf {e}}_{z}\right)\psi \\&={\frac {i}{\hbar }}{\mathbf {\hat {p}}}\psi \end{aligned}}}

where ex, ey and ez are the unit vectors for the three spatial dimensions, hence

${\displaystyle {\mathbf {\hat {p}}}=-i\hbar \nabla }$

This momentum operator is in position space because the partial derivatives were taken with respect to the spatial variables.

## Definition (position space)

For a single particle with no electric charge and no spin, the momentum operator can be written in the position basis as:[2]

${\displaystyle {\mathbf {\hat {p}}}=-i\hbar \nabla }$

where is the gradient operator, ħ is the reduced Planck constant, and i is the imaginary unit.

In one spatial dimension this becomes:

${\displaystyle {\hat {p}}={\hat {p}}_{x}=-i\hbar {\partial \over \partial x}.}$

This is a commonly encountered form of the momentum operator, though not the most general one. For a charged particle Template:Mvar in an electromagnetic field, described by the scalar potential Template:Mvar and vector potential A, the momentum operator must be replaced by:[3]

${\displaystyle {\mathbf {\hat {p}}}=-i\hbar \nabla -q{\mathbf {A}}}$

where the canonical momentum operator is the above momentum operator:

${\displaystyle {\mathbf {\hat {P}}}=-i\hbar \nabla }$

This is of course true for electrically neutral particles also, since the second term vanishes if q = 0 and the original operator appears.

## Properties

### Hermiticity

The momentum operator is always a Hermitian operator when it acts on physical (in particular, normalizable) quantum states.[4]

### Canonical commutation relation

One can easily show that by appropriately using the momentum basis and the position basis:

${\displaystyle \left[{\hat {x}},{\hat {p}}\right]={\hat {x}}{\hat {p}}-{\hat {p}}{\hat {x}}=i\hbar .}$

The Heisenberg uncertainty principle defines limits on how accurately the momentum and position of a single observable system can be known at once. In quantum mechanics, position and momentum are conjugate variables.

### Fourier transform

One can show that the Fourier transform of the momentum in quantum mechanics is the position operator. The Fourier transform turns the momentum-basis into the position-basis. The following discussion uses the bra–ket notation:

${\displaystyle \langle x|{\hat {p}}|\psi \rangle =-i\hbar {\frac {d}{dx}}\psi (x)}$

The same applies for the position operator in the momentum basis:

${\displaystyle \langle p|{\hat {x}}|\psi \rangle =i\hbar {\frac {d}{dp}}\psi (p)}$

and other useful relations:

${\displaystyle \langle p|{\hat {x}}|p'\rangle =i\hbar {\frac {d}{dp}}\delta (p-p')}$
${\displaystyle \langle x|{\hat {p}}|x'\rangle =-i\hbar {\frac {d}{dx}}\delta (x-x')}$

where Template:Mvar stands for Dirac's delta function.

## Derivation from infinitesimal translations

The translation operator is denoted T(ε), where Template:Mvar represents the length of the translation. It satisfies the following identity:

${\displaystyle T(\varepsilon )|\psi \rangle =\int dxT(\varepsilon )|x\rangle \langle x|\psi \rangle }$

that becomes

${\displaystyle \int dx|x+\varepsilon \rangle \langle x|\psi \rangle =\int dx|x\rangle \langle x-\varepsilon |\psi \rangle =\int dx|x\rangle \psi (x-\varepsilon )}$

Assuming the function Template:Mvar to be analytic (i.e. differentiable in some domain of the complex plane), one may expand in a Taylor series about Template:Mvar:

${\displaystyle \psi (x-\varepsilon )=\psi (x)-\varepsilon {\frac {d\psi }{dx}}}$

so for infinitesimal values of Template:Mvar:

${\displaystyle T(\varepsilon )=1-\varepsilon {d \over dx}=1-{i \over \hbar }\varepsilon \left(-i\hbar {d \over dx}\right)}$

As it is known from classical mechanics, the momentum is the generator of translation, so the relation between translation and momentum operators is:

${\displaystyle T(\varepsilon )=1-{i \over \hbar }\varepsilon {\hat {p}}}$

thus

${\displaystyle {\hat {p}}=-i\hbar {d \over dx}.}$

## 4-momentum operator

Inserting the 3d momentum operator above and the energy operator into the 4-momentum (as a 1-form with (+ − − −) metric signature):

${\displaystyle P_{\mu }=\left({\frac {E}{c}},-{\mathbf {p}}\right)}$

obtains the 4-momentum operator;

${\displaystyle {\hat {P}}_{\mu }=\left({\frac {1}{c}}{\hat {E}},-{\mathbf {\hat {p}}}\right)=i\hbar \left({\frac {1}{c}}{\frac {\partial }{\partial t}},\nabla \right)=i\hbar \partial _{\mu }}$

where μ is the 4-gradient, and the becomes + preceding the 3-momentum operator. This operator occurs in relativistic quantum field theory, such as the Dirac equation and other relativistic wave equations, since energy and momentum combine into the 4-momentum vector above, momentum and energy operators correspond to space and time derivatives, and they need to be first order partial derivatives for Lorentz covariance.

The Dirac operator and Dirac slash of the 4-momentum is given by contracting with the gamma matrices:

${\displaystyle \gamma ^{\mu }{\hat {P}}_{\mu }=i\hbar \gamma ^{\mu }\partial _{\mu }={\hat {P}}=i\hbar \partial }$

If the signature was (− + + +), the operator would be

${\displaystyle {\hat {P}}_{\mu }=\left(-{\frac {1}{c}}{\hat {E}},{\mathbf {\hat {p}}}\right)=-i\hbar \left({\frac {1}{c}}{\frac {\partial }{\partial t}},\nabla \right)=-i\hbar \partial _{\mu }}$