Skip to content

A.2

Fourier analysis reference

Transform conventions (including the ħ version for momentum space), a table of standard transform pairs, Parseval and convolution theorems, the δ function and Gaussian integrals — the complete bridge between the position and momentum representations.

Recommended first

  • Calculus (definite integrals and integration by parts)

After this section you should be able to

  • Distinguish the mathematical and quantum-mechanical (ħ) Fourier conventions and convert between them
  • Look up and use the standard transform pairs (Gaussian, square pulse, plane wave, δ function)
  • Use Parseval's theorem to explain why "normalised in position ⇔ normalised in momentum"
  • Evaluate Gaussian integrals with linear and quadratic terms (completing the square)

The physical content of the Fourier transform in one sentence: decompose a function into a superposition of waves of different frequencies. In quantum mechanics, “frequency” is momentum (p=kp=\hbar k), so the Fourier transform is precisely the change-of-basis formula between the position and momentum representations.

Where the main text needs this: 2.10 The free particle uses it to build wave packets; 3.3 Hilbert space interprets it as an expansion in a continuous basis; 7.6 Fermi’s golden rule and the Chapter 11 path integral lean heavily on Gaussian integrals; appendix C.2 The split-operator method uses the FFT to hop back and forth between the two representations.

Transform conventions

The same “Fourier transform” differs between books by a constant or a sign in the exponent, so let us nail down this book’s convention first.

Mathematical convention (symmetric version, kk is the wavenumber):

f~(k)=12πf(x)eikxdx,f(x)=12πf~(k)eikxdk(A.2.1)\tilde f(k)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} f(x)\,\ee^{-\ii kx}\,\dd x, \qquad f(x)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} \tilde f(k)\,\ee^{\ii kx}\,\dd k\tag{A.2.1}

Quantum-mechanical convention (p=kp=\hbar k, momentum wavefunction written ϕ(p)\phi(p)):

ϕ(p)=12πψ(x)eipx/dx,ψ(x)=12πϕ(p)eipx/dp(A.2.2)\phi(p)=\frac{1}{\sqrt{2\pi\hbar}}\int_{-\infty}^{\infty}\psi(x)\,\ee^{-\ii px/\hbar}\,\dd x, \qquad \psi(x)=\frac{1}{\sqrt{2\pi\hbar}}\int_{-\infty}^{\infty}\phi(p)\,\ee^{\ii px/\hbar}\,\dd p\tag{A.2.2}

The prefactor 1/2π1/\sqrt{2\pi\hbar} is there so that both sides are normalised simultaneously (see Parseval below). In Dirac notation this is nothing but inserting a complete set of momentum eigenstates: ϕ(p)=p|ψ\phi(p)=\braket{p}{\psi}, with x|p=eipx//2π\braket{x}{p}=\ee^{\ii px/\hbar}/\sqrt{2\pi\hbar}.

The δ function

The Dirac δ function is not an ordinary function but a “value picker”: whenever it appears inside an integral, it pulls out the value of the integrand at the position of the spike:

f(x)δ(xx0)dx=f(x0)(A.2.3)\int_{-\infty}^{\infty} f(x)\,\delta(x-x_0)\,\dd x=f(x_0)\tag{A.2.3}

Picture it as a limiting spike whose width goes to zero while its area stays 1. Its most important representation is a superposition of plane waves:

δ(x)=12πeikxdk(A.2.4)\delta(x)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\ee^{\ii kx}\,\dd k\tag{A.2.4}

— plane waves of every wavenumber added with equal weight, cancelling everywhere and piling up in phase only at x=0x=0. This formula is the origin of the continuous-basis orthonormality relations x|x=δ(xx)\braket{x}{x'}=\delta(x-x') and p|p=δ(pp)\braket{p}{p'}=\delta(p-p').

Frequently used properties:

PropertyFormula
Scalingδ(ax)=δ(x)/a\delta(ax)=\delta(x)/\lvert a\rvert
Composite argumentδ(g(x))=iδ(xxi)/g(xi)\delta(g(x))=\sum_i \delta(x-x_i)/\lvert g'(x_i)\rvert, where the xix_i are the simple roots of gg
Derivativef(x)δ(xx0)dx=f(x0)\int f(x)\,\delta'(x-x_0)\,\dd x=-f'(x_0)
Step functiondΘ/dx=δ(x)\dd\Theta/\dd x=\delta(x)

Table of standard transform pairs

Listed in the mathematical (symmetric) convention. To switch to the ħ convention, set kp/k\to p/\hbar and multiply the result by 1/1/\sqrt{\hbar}.

f(x)f(x)f~(k)\tilde f(k)Remarks
ex2/2σ2\ee^{-x^2/2\sigma^2}σeσ2k2/2\sigma\,\ee^{-\sigma^2k^2/2}Gaussian stays Gaussian; widths are reciprocal, σxσk=1\sigma_x\sigma_k=1
eik0x\ee^{ik_0x}2πδ(kk0)\sqrt{2\pi}\,\delta(k-k_0)Plane wave = a single momentum
δ(xx0)\delta(x-x_0)eikx0/2π\ee^{-\ii kx_0}/\sqrt{2\pi}A position eigenstate contains every momentum
Square pulse (equal to 1 for x<a\lvert x\rvert<a)2/πsinkak\sqrt{2/\pi}\,\dfrac{\sin ka}{k}sinc-shaped; the mathematical prototype of single-slit diffraction
eax\ee^{-a\lvert x\rvert}2/πaa2+k2\sqrt{2/\pi}\,\dfrac{a}{a^2+k^2}Lorentzian; exponential decay ⇔ power-law tails
f(xx0)f(x-x_0)eikx0f~(k)\ee^{-\ii kx_0}\tilde f(k)Translation ⇔ attaching a phase
f(x)f'(x)ikf~(k)\ii k\tilde f(k)Differentiation ⇔ multiplying by ik\ii k; the algorithmic basis of C.2

The first row deserves a second look: the narrower a Gaussian wave packet (small σ\sigma), the wider its momentum distribution — this is the uncertainty principle showing up at the Fourier level, and the Gaussian saturates the bound exactly, σxσp=/2\sigma_x\sigma_p=\hbar/2.

Parseval’s theorem and the convolution theorem

Parseval’s (Plancherel’s) theorem: the “total amount” is unchanged by the transform,

f(x)2dx=f~(k)2dk(A.2.5)\int_{-\infty}^{\infty}\lvert f(x)\rvert^2\,\dd x=\int_{-\infty}^{\infty}\lvert\tilde f(k)\rvert^2\,\dd k\tag{A.2.5}

Physical translation: if ψ(x)\psi(x) is normalised, then ϕ(p)\phi(p) is automatically normalised too — total probability is 1 in the position representation and still 1 after switching to momentum. That is why ϕ(p)2\lvert\phi(p)\rvert^2 can be directly interpreted as the probability density for momentum.

Convolution theorem: the convolution of two functions, (fg)(x)=f(y)g(xy)dy(f*g)(x)=\int f(y)\,g(x-y)\,\dd y (shift gg everywhere and stack the copies weighted by ff), becomes a product after transforming:

fg~(k)=2πf~(k)g~(k)(A.2.6)\widetilde{f*g}(k)=\sqrt{2\pi}\,\tilde f(k)\,\tilde g(k)\tag{A.2.6}

The use: convolutions are hard, products are easy, so “transform — multiply — transform back” is the standard play. Numerically this is exactly how FFT convolution works, and it is what the C.2 split-operator method does at every step.

Gaussian integrals

The most frequently occurring family of integrals in quantum mechanics. The basic formula (Rea>0\operatorname{Re}a>0):

eax2dx=πa(A.2.7)\int_{-\infty}^{\infty}\ee^{-ax^2}\,\dd x=\sqrt{\frac{\pi}{a}}\tag{A.2.7}

Section 100 of 106 · use to turn the page