Skip to content

Momentum-Space Representation

The momentum-space wavefunction is the coordinate function of an abstract state in the generalized momentum basis. For a particle on a line,

ϕ(p)=⟨p∣ψ⟩.\phi(p)=\langle p|\psi\rangle.

The complex number ϕ(p)\phi(p) is the probability amplitude density associated with the momentum label pp. The function p↦ϕ(p)p\mapsto\phi(p) represents the same state vector as the position-space wavefunction ψ(x)=⟨x∣ψ⟩\psi(x)=\langle x|\psi\rangle; the two functions are related by a unitary Fourier transform.

The basic dictionary is:

  • position space makes XX multiplicative and PP differential;
  • momentum space makes PP multiplicative and XX differential;
  • free dynamics is diagonal in momentum;
  • spatially varying potentials generally couple different momenta;
  • probabilities and expectation values are representation independent when the state, operators, and measure are transformed consistently.

Momentum space is therefore both a conceptual second example of a continuous representation and a practical language adapted to translation symmetry, free motion, wave packets, and scattering.

Let PP denote the momentum operator. Its formal generalized eigenkets satisfy

P∣p⟩=p∣p⟩.P|p\rangle=p|p\rangle.

For momentum on the full line, the continuum normalization convention is

⟨p∣p′⟩=δ(p−p′)\langle p|p'\rangle = \delta(p-p')

and the formal completeness relation is

I=∫−∞∞dp ∣p⟩⟨p∣.I = \int_{-\infty}^{\infty} dp\,|p\rangle\langle p|.

Applying completeness to a normalizable state gives

∣ψ⟩=∫−∞∞dp ϕ(p)∣p⟩.|\psi\rangle = \int_{-\infty}^{\infty} dp\,\phi(p)|p\rangle.

These formulas are the continuous counterparts of

cn=⟨en∣ψ⟩,∣ψ⟩=∑ncn∣en⟩.c_n=\langle e_n|\psi\rangle, \qquad |\psi\rangle=\sum_n c_n|e_n\rangle.

An exact ∣p⟩|p\rangle is not normally a finite-norm vector in the physical Hilbert space. It is a generalized eigenvector used under pairings and integrals. The distributional interpretation belongs to Generalized Eigenvectors; the normalization comparison belongs to Normalization.

The projection

ϕ(p)=⟨p∣ψ⟩\phi(p)=\langle p|\psi\rangle

extracts the momentum component of ∣ψ⟩|\psi\rangle. Conversely, the complete function reconstructs the vector:

∣ψ⟩=∫dp ϕ(p)∣p⟩.|\psi\rangle = \int dp\,\phi(p)|p\rangle.

For any two normalizable states,

⟨χ∣ψ⟩=∫dp χ(p)∗ϕ(p).\langle\chi|\psi\rangle = \int dp\,\chi(p)^*\phi(p).

In particular,

∥ψ∥2=∫dp ∣ϕ(p)∣2.\|\psi\|^2 = \int dp\,|\phi(p)|^2.

Thus the momentum representation identifies the abstract one-particle Hilbert space with an appropriate L2L^2 space of momentum amplitudes. As in position space, functions that differ only on a set of momentum measure zero represent the same Hilbert-space vector.

The symbol ϕ\phi is conventional rather than compulsory. Many texts use ψ~(p)\widetilde\psi(p), ψ(p)\psi(p), or ψp(p)\psi_p(p). A trustworthy formula declares what its argument labels and which Fourier convention is being used.

The site fixes the one-dimensional overlap convention

⟨x∣p⟩=12πℏeipx/ℏ.\langle x|p\rangle = \frac{1}{\sqrt{2\pi\hbar}} e^{ipx/\hbar}.

Taking the adjoint gives

⟨p∣x⟩=12πℏe−ipx/ℏ.\langle p|x\rangle = \frac{1}{\sqrt{2\pi\hbar}} e^{-ipx/\hbar}.

These kernels are delta normalized. Formally,

⟨p∣p′⟩=∫dx ⟨p∣x⟩⟨x∣p′⟩=12πℏ∫dx ei(p′−p)x/ℏ=δ(p−p′).\begin{aligned} \langle p|p'\rangle &= \int dx\, \langle p|x\rangle \langle x|p'\rangle \\ &= \frac{1}{2\pi\hbar} \int dx\, e^{i(p'-p)x/\hbar} \\ &= \delta(p-p'). \end{aligned}

Similarly,

∫dp ⟨x∣p⟩⟨p∣x′⟩=δ(x−x′).\int dp\, \langle x|p\rangle \langle p|x'\rangle = \delta(x-x').

The normalization and phase of ⟨x∣p⟩\langle x|p\rangle determine the transform pair, the dimensions of ϕ(p)\phi(p), and the represented form of XX. Mixing an overlap kernel from one convention with a transform formula from another is a common source of missing signs and factors of 2π2\pi or ℏ\hbar.

Insert the position identity into ϕ(p)=⟨p∣ψ⟩\phi(p)=\langle p|\psi\rangle:

ϕ(p)=∫dx ⟨p∣x⟩⟨x∣ψ⟩=12πℏ∫−∞∞dx e−ipx/ℏψ(x).\begin{aligned} \phi(p) &= \int dx\, \langle p|x\rangle \langle x|\psi\rangle \\ &= \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} dx\,e^{-ipx/\hbar}\psi(x). \end{aligned}

The inverse relation follows by inserting the momentum identity:

ψ(x)=∫dp ⟨x∣p⟩⟨p∣ψ⟩=12πℏ∫−∞∞dp eipx/ℏϕ(p).\begin{aligned} \psi(x) &= \int dp\, \langle x|p\rangle \langle p|\psi\rangle \\ &= \frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty} dp\,e^{ipx/\hbar}\phi(p). \end{aligned}

The signs and symmetric normalization factors are fixed in Fourier Transform Conventions. The transform theorem, inversion conditions, and distributional extensions are canonical in Fourier Transform.

The two wavefunctions are not independent initial data. Either complete representative determines the other, up to the same almost-everywhere and global-phase equivalences that apply to the abstract state.

For suitable wavefunctions, direct substitution gives

∫dp χ(p)∗ϕ(p)=∫dx η(x)∗ψ(x),\int dp\,\chi(p)^*\phi(p) = \int dx\,\eta(x)^*\psi(x),

where η\eta and χ\chi represent one state in position and momentum space, and ψ\psi and ϕ\phi represent the other. Setting the two states equal gives

∫dp ∣ϕ(p)∣2=∫dx ∣ψ(x)∣2.\int dp\,|\phi(p)|^2 = \int dx\,|\psi(x)|^2.

This is the quantum-mechanical form of the Plancherel identity. It means that Fourier transformation is a unitary change of representation: it preserves inner products, norms, orthogonality, and therefore transition probabilities.

The theorem-level statement and proof are in Plancherel and Parseval Theorems. Here its physical role is the important point: normalization does not have to be reimposed after an exact representation change.

For a normalized pure state,

∫−∞∞dp ∣ϕ(p)∣2=1.\int_{-\infty}^{\infty} dp\,|\phi(p)|^2=1.

The momentum Born rule assigns

Pr⁡(P∈Δ)=∫Δdp ∣ϕ(p)∣2\Pr(P\in\Delta) = \int_\Delta dp\,|\phi(p)|^2

to a measurable momentum region Δ\Delta. Thus ∣ϕ(p)∣2|\phi(p)|^2 is a probability density with respect to dpdp, not a probability at one exact value.

In one dimension,

[∣ϕ(p)∣2]=[p]−1,[ϕ(p)]=[p]−1/2.[|\phi(p)|^2]=[p]^{-1}, \qquad [\phi(p)]=[p]^{-1/2}.

These dimensions ensure that ∣ϕ(p)∣2dp|\phi(p)|^2dp is dimensionless. They differ from the dimensions of the position wavefunction, which are L−1/2L^{-1/2} on the line. A unitary transform preserves norms, not the physical units of the coordinate functions.

The probability of obtaining one exact value is normally zero for an absolutely continuous distribution, even if ∣ϕ(p0)∣2|\phi(p_0)|^2 is nonzero. The full statement in terms of spectral projectors is in Born Rule for Continuous Spectra.

For the generalized momentum eigenket ∣p0⟩|p_0\rangle, the momentum representative is

⟨p∣p0⟩=δ(p−p0),\langle p|p_0\rangle = \delta(p-p_0),

while the position representative is the plane wave

⟨x∣p0⟩=12πℏeip0x/ℏ.\langle x|p_0\rangle = \frac{1}{\sqrt{2\pi\hbar}} e^{ip_0x/\hbar}.

Neither representative is an ordinary L2L^2 wavefunction. The delta distribution cannot be squared to form a normalizable density, and the plane wave has constant modulus on the full line.

A physical normalizable state is instead a wave packet:

ψ(x)=12πℏ∫dp eipx/ℏϕ(p),\psi(x) = \frac{1}{\sqrt{2\pi\hbar}} \int dp\, e^{ipx/\hbar}\phi(p),

with ϕ∈L2(R,dp)\phi\in L^2(\mathbb R,dp). A packet concentrated near p0p_0 has an approximately sharp momentum but remains a superposition of continuum eigenkets. Plane-wave normalization and box regularization are developed in Plane Waves and Delta Normalization.

Wave number kk and physical momentum pp are related by

p=ℏk.p=\hbar k.

Changing the spectral label changes the density and the normalized generalized basis. If

⟨k∣k′⟩=δ(k−k′),\langle k|k'\rangle=\delta(k-k'),

then in one dimension

∣k⟩=ℏ ∣p=ℏk⟩.|k\rangle = \sqrt{\hbar}\,|p=\hbar k\rangle.

For

ψ~(k)=⟨k∣ψ⟩,\widetilde\psi(k)=\langle k|\psi\rangle,

the amplitudes obey

ψ~(k)=ℏ ϕ(ℏk).\widetilde\psi(k) = \sqrt{\hbar}\,\phi(\hbar k).

Consequently,

∣ψ~(k)∣2dk=∣ϕ(p)∣2dp.|\widetilde\psi(k)|^2dk = |\phi(p)|^2dp.

The two functions are therefore not related by merely renaming the horizontal axis. Their Jacobian factor is required by probability conservation. In dd dimensions the corresponding factor is ℏd/2\hbar^{d/2}.

With the normalized kk basis,

⟨x∣k⟩=12πeikx,\langle x|k\rangle = \frac{1}{\sqrt{2\pi}}e^{ikx},

whereas the normalized pp basis carries (2πℏ)−1/2(2\pi\hbar)^{-1/2}. Declaring the spectral variable prevents dimensional mistakes.

For one spinless particle in three Cartesian dimensions,

ϕ(p)=⟨p∣ψ⟩\phi(\mathbf p) = \langle\mathbf p|\psi\rangle

and

⟨r∣p⟩=1(2πℏ)3/2eip⋅r/ℏ.\langle\mathbf r|\mathbf p\rangle = \frac{1}{(2\pi\hbar)^{3/2}} e^{i\mathbf p\cdot\mathbf r/\hbar}.

The transform pair is

ϕ(p)=1(2πℏ)3/2∫d3r e−ip⋅r/ℏψ(r)\phi(\mathbf p) = \frac{1}{(2\pi\hbar)^{3/2}} \int d^3r\, e^{-i\mathbf p\cdot\mathbf r/\hbar} \psi(\mathbf r)

and

ψ(r)=1(2πℏ)3/2∫d3p eip⋅r/ℏϕ(p).\psi(\mathbf r) = \frac{1}{(2\pi\hbar)^{3/2}} \int d^3p\, e^{i\mathbf p\cdot\mathbf r/\hbar} \phi(\mathbf p).

Normalization reads

∫d3p ∣ϕ(p)∣2=1,\int d^3p\,|\phi(\mathbf p)|^2=1,

so

[ϕ(p)]=[p]−3/2.[\phi(\mathbf p)]=[p]^{-3/2}.

For spherical momentum coordinates,

d3p=p2sin⁡θp dp dθp dφp.d^3p = p^2\sin\theta_p\, dp\,d\theta_p\,d\varphi_p.

As in position space, the probability density must be interpreted together with its measure.

If momentum does not form a complete label, additional indices remain. A spin-1/21/2 particle has components

ϕs(p)=⟨p,s∣ψ⟩,s∈{↑,↓}.\phi_s(\mathbf p) = \langle\mathbf p,s|\psi\rangle, \qquad s\in\{\uparrow,\downarrow\}.

When spin is not resolved, the momentum density is

ρ(p)=∑s∣ϕs(p)∣2.\rho(\mathbf p) = \sum_s|\phi_s(\mathbf p)|^2.

Normalization becomes

∑s∫d3p ∣ϕs(p)∣2=1.\sum_s \int d^3p\, |\phi_s(\mathbf p)|^2 =1.

For NN distinguishable particles, the momentum wavefunction is

Φ(p1,…,pN)=⟨p1,…,pN∣Ψ⟩.\Phi(\mathbf p_1,\ldots,\mathbf p_N) = \langle \mathbf p_1,\ldots,\mathbf p_N |\Psi\rangle.

Its domain is the 3N3N-dimensional momentum configuration space. Identical particles add exchange symmetry, and second quantization reorganizes the same information using modes and occupation numbers. Those many-body structures belong to Momentum-Space Representation in Many-Body Theory.

The delta normalization of ∣p⟩|p\rangle does not fix its pp-dependent phase. One may rephase the generalized basis as

∣p⟩′=eiα(p)∣p⟩.|p\rangle' = e^{i\alpha(p)}|p\rangle.

The same abstract state then has components

ϕ′(p)=e−iα(p)ϕ(p).\phi'(p) = e^{-i\alpha(p)}\phi(p).

The momentum density is unchanged, but overlap kernels and represented operators change. Under this rephasing, the standard position operator iℏ d/dpi\hbar\,d/dp becomes

X′=iℏddp−ℏα′(p).X' = i\hbar\frac{d}{dp} - \hbar\alpha'(p).

This is why the plane-wave overlap convention is part of the representation, not a disposable typographical choice. A basis rephasing is passive and leaves all predictions invariant when every represented object is transformed. By contrast, multiplying only the state amplitude by a pp-dependent phase is an active change of state.

In its own spectral representation, momentum acts by multiplication:

(Pϕ)(p)=pϕ(p).(P\phi)(p)=p\phi(p).

More generally, the spectral calculus gives

(f(P)ϕ)(p)=f(p)ϕ(p).(f(P)\phi)(p) = f(p)\phi(p).

Examples include

(P2ϕ)(p)=p2ϕ(p)(P^2\phi)(p)=p^2\phi(p)

and the nonrelativistic kinetic energy

(Tϕ)(p)=p22mϕ(p).(T\phi)(p) = \frac{p^2}{2m}\phi(p).

Diagonalization here does not mean that every momentum-space function is an eigenstate. It means the operator acts pointwise in the spectral coordinate. An exact eigenstate would be concentrated at one pp as a delta distribution.

Expectation values of suitable functions of momentum are ordinary weighted integrals:

⟨f(P)⟩ψ=∫dp ∣ϕ(p)∣2f(p).\langle f(P)\rangle_\psi = \int dp\, |\phi(p)|^2 f(p).

For an unbounded f(P)f(P), this expression also states a domain requirement: the integral must exist for the state under consideration.

Using the standard overlap phase,

⟨p∣X∣ψ⟩=iℏdϕdp.\langle p|X|\psi\rangle = i\hbar\frac{d\phi}{dp}.

Thus

(Xϕ)(p)=iℏdϕdp(X\phi)(p) = i\hbar\frac{d\phi}{dp}

on an appropriate domain. The canonical commutator is visible directly:

([X,P]ϕ)(p)=iℏddp(pϕ(p))−p iℏdϕdp=iℏϕ(p).\begin{aligned} ([X,P]\phi)(p) &= i\hbar\frac{d}{dp} \bigl(p\phi(p)\bigr) - p\,i\hbar\frac{d\phi}{dp} \\ &= i\hbar\phi(p). \end{aligned}

Therefore

[X,P]=iℏI[X,P]=i\hbar I

on vectors for which both compositions are defined.

The phrase “on an appropriate domain” is essential. A derivative formula does not by itself specify a self-adjoint operator. Regularity, endpoint behavior, and boundary conditions matter, especially when momentum space is compact or restricted. The general operator viewpoint is in Operator Representations.

An abstract operator AA has momentum-space kernel

A(p,p′)=⟨p∣A∣p′⟩.A(p,p') = \langle p|A|p'\rangle.

Its action is formally

(Aϕ)(p)=∫dp′ A(p,p′)ϕ(p′).(A\phi)(p) = \int dp'\, A(p,p')\phi(p').

Momentum itself has diagonal kernel

⟨p∣P∣p′⟩=p δ(p−p′),\langle p|P|p'\rangle = p\,\delta(p-p'),

and a function of momentum has

⟨p∣f(P)∣p′⟩=f(p)δ(p−p′).\langle p|f(P)|p'\rangle = f(p)\delta(p-p').

This continuum kernel notation is the integral analogue of matrix multiplication. Diagonal kernels act by multiplication; off-diagonal kernels mix amplitudes at different momentum labels.

Let V(X)V(X) act by multiplication with V(x)V(x) in position space, and define

V~(q)=12πℏ∫dx e−iqx/ℏV(x).\widetilde V(q) = \frac{1}{\sqrt{2\pi\hbar}} \int dx\,e^{-iqx/\hbar}V(x).

Then, whenever the expressions are ordinary functions or are interpreted distributionally,

(V(X)ϕ)(p)=12πℏ∫dp′ V~(p−p′)ϕ(p′).\bigl(V(X)\phi\bigr)(p) = \frac{1}{\sqrt{2\pi\hbar}} \int dp'\, \widetilde V(p-p')\phi(p').

Equivalently,

⟨p∣V(X)∣p′⟩=12πℏV~(p−p′).\langle p|V(X)|p'\rangle = \frac{1}{\sqrt{2\pi\hbar}} \widetilde V(p-p').

The argument p−p′p-p' is the momentum transferred by the potential. A constant potential is diagonal because its transform is proportional to a delta distribution. A rapidly varying potential contains broad momentum-transfer components and can couple widely separated momenta.

For real V(x)V(x),

V~(q)∗=V~(−q),\widetilde V(q)^* = \widetilde V(-q),

which is the kernel condition associated with Hermiticity. The full convolution calculation is developed in Momentum Representation as a Computational Tool.

For

H=P22m+V(X),H=\frac{P^2}{2m}+V(X),

the momentum-space Schrödinger equation is

iℏ∂ϕ(p,t)∂t=p22mϕ(p,t)+12πℏ∫dp′ V~(p−p′)ϕ(p′,t).\begin{aligned} i\hbar\frac{\partial\phi(p,t)}{\partial t} &= \frac{p^2}{2m}\phi(p,t) \\ &\quad+ \frac{1}{\sqrt{2\pi\hbar}} \int dp'\, \widetilde V(p-p')\phi(p',t). \end{aligned}

For a free particle, V=0V=0, so each momentum component evolves independently:

ϕ(p,t)=exp⁡(−ip2t2mℏ)ϕ(p,0).\phi(p,t) = \exp\left( -\frac{ip^2t}{2m\hbar} \right) \phi(p,0).

Consequently,

∣ϕ(p,t)∣2=∣ϕ(p,0)∣2.|\phi(p,t)|^2 = |\phi(p,0)|^2.

The free momentum distribution is time independent even though the position-space packet can translate and spread. The different momentum components acquire different phases, and their interference determines the later position profile. The canonical free-motion analysis belongs to Free Particle.

The unitary translation by aa is

T(a)=e−iaP/ℏ.T(a) = e^{-iaP/\hbar}.

Because PP is multiplicative in momentum space,

(T(a)ϕ)(p)=e−iap/ℏϕ(p).(T(a)\phi)(p) = e^{-iap/\hbar}\phi(p).

In position space the same active transformation is

(T(a)ψ)(x)=ψ(x−a).(T(a)\psi)(x) = \psi(x-a).

A spatial translation therefore leaves the momentum density unchanged:

∣e−iap/ℏϕ(p)∣2=∣ϕ(p)∣2.|e^{-iap/\hbar}\phi(p)|^2 = |\phi(p)|^2.

It changes relative phases across momentum components, which encode the packet’s position. Conversely, the momentum-displacement unitary

B(q)=eiqX/ℏB(q)=e^{iqX/\hbar}

acts as

(B(q)ϕ)(p)=ϕ(p−q).(B(q)\phi)(p)=\phi(p-q).

The generator and conservation-law interpretation is canonical in Translations and Momentum.

Translation Symmetry and Momentum Conservation

Section titled “Translation Symmetry and Momentum Conservation”

If the Hamiltonian is invariant under all spatial translations, then

[H,P]=0.[H,P]=0.

Momentum is conserved, and the Hamiltonian can be decomposed into momentum sectors. For one free particle this decomposition is completely diagonal:

H(p,p′)=p22mδ(p−p′).H(p,p') = \frac{p^2}{2m} \delta(p-p').

For systems with internal degrees of freedom, HH may be a finite matrix within each momentum sector. For interacting many-particle systems, translation invariance conserves total momentum even when individual particle momenta are exchanged.

On a lattice, continuous momentum is replaced by crystal momentum in the Brillouin zone, and momenta differing by a reciprocal-lattice vector are equivalent labels. Those are symmetry and many-body extensions, not changes to the basic representational principle.

Momentum need not have a continuous spectrum. On a periodic interval of length LL, normalized modes are

⟨x∣n⟩=1Leipnx/ℏ,pn=2πℏnL,\langle x|n\rangle = \frac{1}{\sqrt L} e^{ip_nx/\hbar}, \qquad p_n=\frac{2\pi\hbar n}{L},

with n∈Zn\in\mathbb Z. The momentum representation is then the discrete sequence

cn=⟨n∣ψ⟩c_n=\langle n|\psi\rangle

with

∑n∈Z∣cn∣2=1.\sum_{n\in\mathbb Z}|c_n|^2=1.

The spacing is

Δp=2πℏL.\Delta p=\frac{2\pi\hbar}{L}.

In a large periodic box, continuum and discrete amplitudes are related schematically by

cn≃Δp ϕ(pn),c_n \simeq \sqrt{\Delta p}\,\phi(p_n),

so that

∑n∣cn∣2⟶∫dp ∣ϕ(p)∣2.\sum_n|c_n|^2 \longrightarrow \int dp\,|\phi(p)|^2.

Boundary conditions determine which momentum operator and spectrum are available. One should not import the periodic-box spectrum into a hard-wall problem without checking the operator domain.

Consider the normalized position wavefunction

ψ(x)=1(2πσx2)1/4exp⁡[−(x−x0)24σx2+ip0xℏ].\psi(x) = \frac{1}{(2\pi\sigma_x^2)^{1/4}} \exp\left[ -\frac{(x-x_0)^2}{4\sigma_x^2} +\frac{ip_0x}{\hbar} \right].

Its momentum representative is

ϕ(p)=(2σx2πℏ2)1/4×exp⁡[−σx2(p−p0)2ℏ2−i(p−p0)x0ℏ].\begin{aligned} \phi(p) &= \left( \frac{2\sigma_x^2}{\pi\hbar^2} \right)^{1/4} \\ &\quad\times \exp\left[ -\frac{\sigma_x^2(p-p_0)^2}{\hbar^2} -\frac{i(p-p_0)x_0}{\hbar} \right]. \end{aligned}

The momentum density is centered at p0p_0:

∣ϕ(p)∣2=2σx2πℏ2exp⁡[−2σx2(p−p0)2ℏ2].|\phi(p)|^2 = \sqrt{\frac{2\sigma_x^2}{\pi\hbar^2}} \exp\left[ -\frac{2\sigma_x^2(p-p_0)^2}{\hbar^2} \right].

Its standard deviation is

Δp=ℏ2σx.\Delta p = \frac{\hbar}{2\sigma_x}.

The position shift x0x_0 appears only as a linear phase in momentum space, while the momentum shift p0p_0 appears as the center of the momentum density. This is the translation dictionary in a concrete normalizable state.

The product Δx Δp=ℏ/2\Delta x\,\Delta p=\hbar/2 is the Gaussian equality case of the uncertainty relation. Its canonical derivation is in Position–Momentum Uncertainty.

A general state is represented by a density operator ρ\rho. Its momentum kernel is

ρ(p,p′)=⟨p∣ρ∣p′⟩.\rho(p,p') = \langle p|\rho|p'\rangle.

The trace condition becomes

∫dp ρ(p,p)=1,\int dp\,\rho(p,p)=1,

and the momentum density is

w(p)=ρ(p,p).w(p)=\rho(p,p).

For a pure state,

ρψ(p,p′)=ϕ(p)ϕ(p′)∗.\rho_\psi(p,p') = \phi(p)\phi(p')^*.

For an ensemble,

ρ(p,p′)=∑jqjϕj(p)ϕj(p′)∗.\rho(p,p') = \sum_j q_j \phi_j(p)\phi_j(p')^*.

The diagonal determines momentum-measurement probabilities. Off-diagonal entries encode coherence between different momenta and matter for observables that are not diagonal in momentum, including position. A single momentum wavefunction does not represent a general mixed state.

Momentum representation is especially effective when:

  • the Hamiltonian is dominated by functions of momentum;
  • the system has translation symmetry;
  • free propagation is central;
  • a wave packet is prepared by its momentum distribution;
  • scattering is organized by incoming, outgoing, and transferred momentum;
  • convolution in position space becomes multiplication in momentum space;
  • derivatives in position space become algebraic powers of momentum;
  • total momentum sectors simplify a many-body problem.

Position representation is often more effective for localized potentials, spatial boundaries, and direct geometry. Neither representation is more physical in general. The useful choice is the one that makes the state, operator, symmetry, and boundary conditions easiest to express together.

A discrete Fourier transform does not directly produce the continuum function ϕ(p)\phi(p) without scale factors. A spatial grid with spacing Δx\Delta x and length L=NΔxL=N\Delta x determines a reciprocal grid with spacing approximately

Δp=2πℏL\Delta p = \frac{2\pi\hbar}{L}

and a finite Nyquist range set by Δx\Delta x. The exact ordering of positive and negative frequencies depends on the transform library.

Numerical checks should include:

  • discrete norm preservation with the chosen quadrature weights;
  • the factors of Δx\Delta x, Δp\Delta p, 2π2\pi, and ℏ\hbar;
  • consistent frequency ordering or an explicit shift operation;
  • enough position range to suppress wraparound;
  • enough spatial resolution to avoid aliasing high momenta;
  • boundary conditions compatible with the discrete transform.

The fast Fourier transform is an algorithm for a discrete transform. It does not remove the need to map discrete coefficients to dimensionful continuum amplitudes.

A correct position-to-momentum transformation preserves:

⟨χ∣ψ⟩,∥ψ∥,Pr⁡(physical event),⟨ψ∣A∣ψ⟩,spec⁡(A).\begin{gathered} \langle\chi|\psi\rangle, \qquad \|\psi\|, \qquad \Pr(\text{physical event}),\\ \langle\psi|A|\psi\rangle, \qquad \operatorname{spec}(A). \end{gathered}

Useful practical checks are:

  1. verify the declared Fourier convention;
  2. check that the exponent is dimensionless;
  3. check the physical units of the transformed amplitude;
  4. verify norm preservation;
  5. transform the operator as well as the state;
  6. distinguish an ordinary L2L^2 state from a delta-normalized eigenfunction;
  7. include Jacobians when relabelling pp by kk, energy, or angular variables;
  8. verify domain and boundary assumptions for derivative operators.
  • Treating ϕ(p)\phi(p) as a different physical state from ψ(x)\psi(x).
  • Calling ∣ϕ(p)∣2|\phi(p)|^2 a probability rather than a density with respect to dpdp.
  • Treating an exact momentum eigenket or delta distribution as a normalizable wavefunction.
  • Confusing momentum pp with wave number kk and omitting the Jacobian factor.
  • Dropping ℏ\hbar from eipx/ℏe^{ipx/\hbar} while still integrating over pp.
  • Combining forward and inverse transforms from different sign or normalization conventions.
  • Assuming that every momentum spectrum is continuous, regardless of boundary conditions.
  • Forgetting that a local potential becomes a momentum-transfer kernel rather than multiplication by V(p)V(p).
  • Using X=iℏ d/dpX=i\hbar\,d/dp without checking its domain or the phase convention of the momentum basis.
  • Interpreting a momentum-dependent phase as physically irrelevant merely because ∣ϕ(p)∣2|\phi(p)|^2 is unchanged.
  • Reading a raw FFT output as a continuum momentum wavefunction without grid weights and units.

This page owns the conceptual and operational momentum-space dictionary. It does not duplicate several neighboring developments:

  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958, Chapters II and III.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2021, Chapters 1 and 2.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994, Chapters 1, 4, and 5.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014, Chapters 2 and 3.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013, Chapters 3, 7, and 10.
  • G. B. Folland, Fourier Analysis and Its Applications, American Mathematical Society, 1992, Chapters 1 and 2.

Starting from ϕ(p)=⟨p∣ψ⟩\phi(p)=\langle p|\psi\rangle, derive the forward transform by inserting the position identity. Then derive the inverse transform by inserting the momentum identity into ψ(x)=⟨x∣ψ⟩\psi(x)=\langle x|\psi\rangle.

Solution

For the forward direction,

ϕ(p)=⟨p∣ψ⟩=∫dx ⟨p∣x⟩⟨x∣ψ⟩=12πℏ∫dx e−ipx/ℏψ(x).\begin{aligned} \phi(p) &= \langle p|\psi\rangle \\ &= \int dx\, \langle p|x\rangle \langle x|\psi\rangle \\ &= \frac{1}{\sqrt{2\pi\hbar}} \int dx\, e^{-ipx/\hbar}\psi(x). \end{aligned}

For the inverse direction,

ψ(x)=⟨x∣ψ⟩=∫dp ⟨x∣p⟩⟨p∣ψ⟩=12πℏ∫dp eipx/ℏϕ(p).\begin{aligned} \psi(x) &= \langle x|\psi\rangle \\ &= \int dp\, \langle x|p\rangle \langle p|\psi\rangle \\ &= \frac{1}{\sqrt{2\pi\hbar}} \int dp\, e^{ipx/\hbar}\phi(p). \end{aligned}

Both formulas use the same normalized overlap kernel.

Let p=ℏkp=\hbar k. If ϕ(p)\phi(p) is normalized with respect to dpdp, find the normalized amplitude ψ~(k)\widetilde\psi(k) with respect to dkdk. Verify the dimensions and probability equality.

Solution

Probability invariance requires

∣ψ~(k)∣2dk=∣ϕ(p)∣2dp.|\widetilde\psi(k)|^2dk = |\phi(p)|^2dp.

Since dp=ℏ dkdp=\hbar\,dk,

ψ~(k)=ℏ ϕ(ℏk)\widetilde\psi(k) = \sqrt\hbar\,\phi(\hbar k)

up to a convention-dependent phase. Therefore

∫dk ∣ψ~(k)∣2=∫dp ∣ϕ(p)∣2=1.\int dk\,|\widetilde\psi(k)|^2 = \int dp\,|\phi(p)|^2 =1.

In one dimension, [ϕ]=[p]−1/2[\phi]=[p]^{-1/2} and [ψ~]=[k]−1/2=L1/2[\widetilde\psi]=[k]^{-1/2}=L^{1/2}, as required by their respective measures.

Using

(Pϕ)(p)=pϕ(p),(Xϕ)(p)=iℏϕ′(p),(P\phi)(p)=p\phi(p), \qquad (X\phi)(p)=i\hbar\phi'(p),

verify [X,P]ϕ=iℏϕ[X,P]\phi=i\hbar\phi on a common invariant domain.

Solution

Compute

(XPϕ)(p)=iℏddp(pϕ(p))=iℏϕ(p)+iℏpϕ′(p),\begin{aligned} (XP\phi)(p) &= i\hbar\frac{d}{dp} \bigl(p\phi(p)\bigr) \\ &= i\hbar\phi(p) + i\hbar p\phi'(p), \end{aligned}

whereas

(PXϕ)(p)=iℏpϕ′(p).(PX\phi)(p) = i\hbar p\phi'(p).

Subtracting gives

([X,P]ϕ)(p)=iℏϕ(p).([X,P]\phi)(p) = i\hbar\phi(p).

The domain qualification ensures that both compositions and their difference are defined.

Let T(a)=e−iaP/ℏT(a)=e^{-iaP/\hbar}. Show that its momentum-space action is multiplication by a phase and recover (T(a)ψ)(x)=ψ(x−a)(T(a)\psi)(x)=\psi(x-a) by inverse Fourier transform.

Solution

Because PP acts by multiplication,

(T(a)ϕ)(p)=e−iap/ℏϕ(p).(T(a)\phi)(p) = e^{-iap/\hbar}\phi(p).

Transforming back,

(T(a)ψ)(x)=12πℏ∫dp eipx/ℏe−iap/ℏϕ(p)=12πℏ∫dp eip(x−a)/ℏϕ(p)=ψ(x−a).\begin{aligned} (T(a)\psi)(x) &= \frac{1}{\sqrt{2\pi\hbar}} \int dp\, e^{ipx/\hbar} e^{-iap/\hbar} \phi(p) \\ &= \frac{1}{\sqrt{2\pi\hbar}} \int dp\, e^{ip(x-a)/\hbar} \phi(p) \\ &= \psi(x-a). \end{aligned}

The phase changes the packet’s position while leaving its momentum density unchanged.

Solve the free momentum-space Schrödinger equation for arbitrary initial ϕ(p,0)\phi(p,0). Show that every moment of the momentum distribution that exists is time independent.

Solution

For each pp,

iℏ∂ϕ(p,t)∂t=p22mϕ(p,t),i\hbar\frac{\partial\phi(p,t)}{\partial t} = \frac{p^2}{2m}\phi(p,t),

so

ϕ(p,t)=e−ip2t/(2mℏ)ϕ(p,0).\phi(p,t) = e^{-ip^2t/(2m\hbar)} \phi(p,0).

The phase has unit modulus, hence

∣ϕ(p,t)∣2=∣ϕ(p,0)∣2.|\phi(p,t)|^2 = |\phi(p,0)|^2.

Therefore, whenever the integral exists,

⟨Pn⟩t=∫dp pn∣ϕ(p,t)∣2=⟨Pn⟩0.\langle P^n\rangle_t = \int dp\,p^n|\phi(p,t)|^2 = \langle P^n\rangle_0.

6. Momentum transfer from a cosine potential

Section titled “6. Momentum transfer from a cosine potential”

Let

V(x)=V0cos⁡(q0x/ℏ).V(x)=V_0\cos(q_0x/\hbar).

Show directly, without manipulating products of delta distributions, that the potential term couples ϕ(p)\phi(p) only to amplitudes displaced by ±q0\pm q_0.

Solution

Write

V(x)=V02(eiq0x/ℏ+e−iq0x/ℏ).V(x) = \frac{V_0}{2} \left( e^{iq_0x/\hbar} +e^{-iq_0x/\hbar} \right).

Multiplication by eiq0x/ℏe^{iq_0x/\hbar} in position space shifts the momentum amplitude by q0q_0, while multiplication by the opposite phase shifts it by −q0-q_0. Therefore

(V(X)ϕ)(p)=V02[ϕ(p−q0)+ϕ(p+q0)].\bigl(V(X)\phi\bigr)(p) = \frac{V_0}{2} \left[ \phi(p-q_0) + \phi(p+q_0) \right].

The periodic potential transfers momentum only in the two amounts ±q0\pm q_0.

Fourier transform

ψ(x)=1(2πσx2)1/4exp⁡[−(x−x0)24σx2+ip0xℏ]\psi(x) = \frac{1}{(2\pi\sigma_x^2)^{1/4}} \exp\left[ -\frac{(x-x_0)^2}{4\sigma_x^2} +\frac{ip_0x}{\hbar} \right]

and identify the momentum mean, variance, and phase associated with x0x_0.

Solution

Completing the square in the Gaussian integral gives

ϕ(p)=(2σx2πℏ2)1/4×exp⁡[−σx2(p−p0)2ℏ2−i(p−p0)x0ℏ].\begin{aligned} \phi(p) &= \left( \frac{2\sigma_x^2}{\pi\hbar^2} \right)^{1/4} \\ &\quad\times \exp\left[ -\frac{\sigma_x^2(p-p_0)^2}{\hbar^2} -\frac{i(p-p_0)x_0}{\hbar} \right]. \end{aligned}

Hence

∣ϕ(p)∣2=2σx2πℏ2exp⁡[−2σx2(p−p0)2ℏ2].|\phi(p)|^2 = \sqrt{\frac{2\sigma_x^2}{\pi\hbar^2}} \exp\left[ -\frac{2\sigma_x^2(p-p_0)^2}{\hbar^2} \right].

The mean is p0p_0 and the variance is

(Δp)2=ℏ24σx2.(\Delta p)^2 = \frac{\hbar^2}{4\sigma_x^2}.

The position center x0x_0 appears as the linear phase e−i(p−p0)x0/ℏe^{-i(p-p_0)x_0/\hbar} and does not alter the momentum density.

For a normalized ϕ(p)\phi(p), show that

ρψ(p,p′)=ϕ(p)ϕ(p′)∗\rho_\psi(p,p') = \phi(p)\phi(p')^*

has unit trace and is idempotent under kernel composition.

Solution

The trace is

Tr⁡ρψ=∫dp ρψ(p,p)=∫dp ∣ϕ(p)∣2=1.\operatorname{Tr}\rho_\psi = \int dp\,\rho_\psi(p,p) = \int dp\,|\phi(p)|^2 =1.

For the square,

(ρψ2)(p,p′)=∫dq ρψ(p,q)ρψ(q,p′)=ϕ(p)ϕ(p′)∗∫dq ∣ϕ(q)∣2=ρψ(p,p′).\begin{aligned} (\rho_\psi^2)(p,p') &= \int dq\, \rho_\psi(p,q) \rho_\psi(q,p') \\ &= \phi(p)\phi(p')^* \int dq\,|\phi(q)|^2 \\ &= \rho_\psi(p,p'). \end{aligned}

The kernel is therefore the momentum representation of the rank-one projector ∣ψ⟩⟨ψ∣|\psi\rangle\langle\psi|.