Dirac Spinors and Plane-Wave Solutions

Once the Dirac equation is in hand, the first thing any physicist does is solve it for a free particle — and the answer is where relativistic quantum mechanics stops being abstract and starts computing. A free electron of definite momentum is a plane wave, and because the Dirac equation is a four-component matrix equation, the plane wave carries a constant four-component passenger called a Dirac spinor. That spinor is not decoration: its structure is exactly what encodes spin-½, and its two energy branches are exactly what forces antimatter into the theory. The move that makes everything tractable is that a plane wave turns the differential Dirac operator into a plain number-valued matrix. Feed in ψ=u(p) e−ip⋅x\psi = u(p)\,e^{-ip\cdot x} and the derivatives just pull down factors of momentum, so "being a solution" collapses to a purely algebraic condition, (p̸−m) u=0(\not{p}-m)\,u=0. This is a matrix equation with a two-dimensional null space for each momentum — precisely the two spin states a spin-½ particle must have. The negative-energy branch, packaged into a second spinor v(p)v(p), gives two more solutions that become the positron. Master these spinors and you own the working language of every scattering calculation. The normalization uˉu=2m\bar{u}u=2m, the completeness relation ∑susuˉs=p̸+m\sum_s u^s\bar{u}^s=\not{p}+m, and the positive density ψ†ψ\psi^\dagger\psi are not formalities: they are the exact identities inserted when you square an amplitude and sum over spins in quantum electrodynamics. We build all four spinors explicitly, verify every identity as an honest 4×4 matrix statement, and watch spin and antimatter fall out of the algebra.

⚡ Key Takeaways
  • Free Dirac solutions factor into a phase times a constant spinor: ψ=us(p) e−ip⋅x\psi = u^s(p)\,e^{-ip\cdot x} (positive energy) and ψ=vs(p) e+ip⋅x\psi = v^s(p)\,e^{+ip\cdot x} (negative energy), with spin label s=1,2s=1,2.
  • The plane-wave ansatz turns the differential equation into the algebraic on-shell conditions (p̸−m)u=0(\not{p}-m)u=0 and (p̸+m)v=0(\not{p}+m)v=0, each with a two-dimensional solution space — the two spin states.
  • In the Dirac (standard) representation the positive-energy spinor is us(p)=(E+m ξs, σ⋅pE+m ξs) ⁣⊤u^s(p)=\bigl(\sqrt{E+m}\,\xi^s,\ \tfrac{\boldsymbol\sigma\cdot\mathbf p}{\sqrt{E+m}}\,\xi^s\bigr)^{\!\top}: a large upper block and a small lower block tied by σ⋅p/(E+m)\boldsymbol\sigma\cdot\mathbf p/(E+m).
  • Normalizations: the Lorentz scalar uˉrus=2m δrs\bar{u}^r u^s = 2m\,\delta^{rs} and the density ur†us=2E δrsu^{r\dagger}u^s = 2E\,\delta^{rs}; the completeness spin sums are ∑susuˉs=p̸+m\sum_s u^s\bar{u}^s=\not{p}+m and ∑svsvˉs=p̸−m\sum_s v^s\bar{v}^s=\not{p}-m.
  • The upper/lower split is the large/small (nonrelativistic) decomposition; the uu/vv split is the particle/antiparticle (energy-sign) decomposition — they are NOT the same thing.
  • The conserved Dirac density j0=ψ†ψ=∑a∣ψa∣2≥0j^0=\psi^\dagger\psi=\sum_a|\psi_a|^2\ge 0 is manifestly positive, curing the Klein–Gordon negative-density disease.
  • Key idea: a free Dirac plane wave is a phase e∓ip⋅xe^{\mp ip\cdot x} times a constant spinor obeying the algebraic condition (p̸∓m)w=0(\not{p}\mp m)w=0, whose two-dimensional null space per momentum encodes spin-½, with the two energy signs becoming electron and positron.
💡 The Core Insight

A free Dirac plane wave factorizes into a spacetime phase e∓ip⋅xe^{\mp ip\cdot x} times a constant four-spinor u(p)u(p) or v(p)v(p) obeying the purely algebraic on-shell condition (p̸∓m)w=0(\not{p}\mp m)w=0; its two-dimensional solution space per momentum encodes spin-½, and the positive- and negative-energy branches map to electrons and positrons.

The physical picture

Think of the Dirac equation as the "square root" of the Klein–Gordon equation. Taking a square root of a differential operator forces matrix-valued coefficients — the gamma matrices — and a matrix operator needs a multi-component object to act on. The four components of ψ\psi are the price of a first-order, Lorentz-covariant equation with a positive probability density. For a plane wave of definite four-momentum, all of spacetime dependence is squeezed into one universal phase e∓ip⋅xe^{\mp ip\cdot x}; the only thing left to determine is a single constant column of four numbers. That constant column lives or dies by a single matrix equation. The Dirac operator, acting on a plane wave, becomes the ordinary c-number matrix p̸∓m\not{p}\mp m, and "solving the Dirac equation" reduces to "find the null space of this 4×4 matrix." On the mass shell that null space is two-dimensional for each energy sign — exactly two spin states for a fermion, no more and no fewer. Spin is not put in by hand; it is the dimension count of a kernel.

Figure
Free Dirac plane wave: phase times spinor
ψ(x) = u(p) · e^(−ip·x) a constant 4-spinor × a universal spacetime phase phase e^(−ip·x): oscillation, all the x,t dependence large (∼ √(E+m) ξ) upper block small (∼ σ·p/(E+m)) lower block u(p): four numbers

A particle at rest makes the anatomy transparent. With p=0\mathbf p=0 the spinor splits cleanly: the positive-energy solution uu lives entirely in the upper two components (two spins), the negative-energy solution vv lives entirely in the lower two (two spins). Boosting to a moving frame mixes the blocks, so a moving electron grows nonzero lower components of order ∣p∣/(E+m)|\mathbf p|/(E+m). That ratio is the relativistic "smallness" parameter: tiny for slow electrons, order one when the kinetic energy rivals the rest mass. The negative-energy vv-solutions are not a pathology to be discarded — they are mandatory for a complete basis, just as they were in the Klein–Gordon theory of "Negative-Energy Solutions and the Interpretive Crisis." Reading the extra branch e+ip⋅xe^{+ip\cdot x} not as a sick electron but as an antiparticle is the conceptual hinge that turns a crisis into the prediction of the positron.

Figure
Rest frame vs boosted: where the components live
particle at rest (p = 0) u: large (ξ) u: large (ξ) small = 0 small = 0 boosted (p ≠ 0) large ∼ √(E+m) large ∼ √(E+m) small ∼ |p|/(E+m) small ∼ |p|/(E+m) the block split (large/small) is the NONRELATIVISTIC decomposition — NOT the particle/antiparticle (u vs v) split. A single electron is pure u.

The core equations

Everything begins with the Dirac equation, the first-order relativistic wave equation whose four-component solutions we are about to construct. We work throughout in natural units ℏ=c=1\hbar=c=1 and the mostly-minus metric gμν=diag(+,−,−,−)g_{\mu\nu}=\mathrm{diag}(+,-,-,-) in the Dirac (standard) representation of the gamma matrices, exactly the conventions pinned in "The Klein-Gordon Equation" and used by "Dirac Gamma Matrices and Clifford Algebra."

(iγμ∂μ−m) ψ(x)=0(1)(i\gamma^\mu\partial_\mu - m)\,\psi(x) = 0 \qquad (1)

where: γμ\gamma^\mu — the four 4×44\times4 Dirac matrices satisfying {γμ,γν}=2gμν1\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}\mathbb{1}; ψ\psi — the four-component Dirac spinor field; mm — the particle mass; ∂μ=∂/∂xμ\partial_\mu=\partial/\partial x^\mu. The equation is first order in all derivatives, which is exactly why its conserved density will turn out positive. Because the equation is linear with constant coefficients, a plane-wave ansatz separates the universal phase from a constant spinor. We keep the physical energy p0=Ep>0p^0=E_p>0 positive in both branches and put the energy-sign information into the choice of spinor and the sign of the exponent.

ψ(+)(x)=us(p) e−ip⋅x,ψ(−)(x)=vs(p) e+ip⋅x(2)\psi^{(+)}(x) = u^s(p)\,e^{-ip\cdot x}, \qquad \psi^{(-)}(x) = v^s(p)\,e^{+ip\cdot x} \qquad (2)

where: us(p)u^s(p) — the positive-energy momentum-space spinor with spin label s=1,2s=1,2; vs(p)v^s(p) — the negative-energy (antiparticle) spinor; p⋅x=p0x0−p⋅xp\cdot x = p^0x^0-\mathbf p\cdot\mathbf x with p0=Ep=p2+m2p^0=E_p=\sqrt{\mathbf p^2+m^2}. Writing the second branch as e+ip⋅xe^{+ip\cdot x} (not e−ip⋅xe^{-ip\cdot x} with p0<0p^0<0) keeps every energy positive and packages the negative-energy content into v(p)v(p) — the bookkeeping that later makes the positron manifest. Substituting the ansatz collapses the differential operator into a matrix. A derivative acting on e−ip⋅xe^{-ip\cdot x} pulls down −ipμ-ip_\mu, so i∂μ→pμi\partial_\mu\to p_\mu and the whole equation becomes purely algebraic.

(p̸−m) us(p)=0,(p̸+m) vs(p)=0,p̸≡γμpμ(3)(\not{p}-m)\,u^s(p)=0, \qquad (\not{p}+m)\,v^s(p)=0, \qquad \not{p}\equiv\gamma^\mu p_\mu \qquad (3)

where: p̸\not{p} — the Feynman slash γμpμ\gamma^\mu p_\mu; the sign flips for vv because e+ip⋅xe^{+ip\cdot x} gives i∂μ→−pμi\partial_\mu\to -p_\mu. Multiplying (p̸−m)u=0(\not p-m)u=0 by (p̸+m)(\not p+m) and using p̸ p̸=p2\not p\,\not p = p^2 forces p2=m2p^2=m^2: the mass shell falls out automatically, reproducing the relativistic energy–momentum relation with no extra input. Now the explicit solutions. In the Dirac representation the standard-basis positive-energy spinor has a large upper block and a small lower block linked by σ⋅p/(E+m)\boldsymbol\sigma\cdot\mathbf p/(E+m).

us(p)=(Ep+m  ξsσ⋅pEp+m ξs),vs(p)=(σ⋅pEp+m ηsEp+m  ηs)(4)u^s(p)=\begin{pmatrix}\sqrt{E_p+m}\,\,\xi^s\\[4pt]\dfrac{\boldsymbol{\sigma}\cdot\mathbf{p}}{\sqrt{E_p+m}}\,\xi^s\end{pmatrix},\qquad v^s(p)=\begin{pmatrix}\dfrac{\boldsymbol{\sigma}\cdot\mathbf{p}}{\sqrt{E_p+m}}\,\eta^s\\[4pt]\sqrt{E_p+m}\,\,\eta^s\end{pmatrix} \qquad (4)

where: ξs,ηs\xi^s,\eta^s — two orthonormal Pauli two-spinors, ξr†ξs=δrs\xi^{r\dagger}\xi^s=\delta^{rs}; σ\boldsymbol\sigma — the Pauli matrices. Note the role reversal: for the antiparticle spinor vv the large block is the lower one — the tell that vv realizes the negative-energy branch, with the particle content pushed downstairs. These spinors carry two distinct normalizations, and confusing them is a classic source of factor errors. The Lorentz scalar uˉu\bar u u is frame-independent; the density u†uu^\dagger u is the timelike current component and scales with energy.

uˉr(p)us(p)=2m δrs,vˉr(p)vs(p)=−2m δrs,uˉrvs=vˉrus=0(5)\bar{u}^r(p)u^s(p)=2m\,\delta^{rs},\qquad \bar{v}^r(p)v^s(p)=-2m\,\delta^{rs},\qquad \bar u^r v^s=\bar v^r u^s=0 \qquad (5)

where: uˉ=u†γ0\bar u = u^\dagger\gamma^0 is the Dirac adjoint; the scalar uˉu=2m\bar u u=2m is Lorentz invariant, while the sign flip vˉv=−2m\bar v v=-2m reflects that vv sits on the opposite energy branch. The particle and antiparticle spinors are orthogonal, so a general solution decomposes cleanly.

ur†(p)us(p)=2Ep δrs,vr†(p)vs(p)=2Ep δrs(6)u^{r\dagger}(p)u^s(p)=2E_p\,\delta^{rs},\qquad v^{r\dagger}(p)v^s(p)=2E_p\,\delta^{rs} \qquad (6)

where: u†u=uˉγ0uu^\dagger u=\bar u\gamma^0 u is the timelike component of the conserved current — the probability density — which is positive and grows as 2Ep2E_p, exactly the relativistic scaling that compensates the Lorentz contraction of a normalization volume. The single most useful identities in all of relativistic scattering are the spin sums, which say that summing the spinor dyadics over spin reconstructs a projector. These are the workhorses of the "Casimir trick."

∑s=1,2us(p) uˉs(p)=p̸+m,∑s=1,2vs(p) vˉs(p)=p̸−m(7)\sum_{s=1,2}u^s(p)\,\bar{u}^s(p)=\not{p}+m,\qquad \sum_{s=1,2}v^s(p)\,\bar{v}^s(p)=\not{p}-m \qquad (7)

where: p̸±m=2m Λ±\not p\pm m = 2m\,\Lambda_\pm are twice the positive/negative energy projectors Λ±=(±p̸+m)/2m\Lambda_\pm=(\pm\not p+m)/2m. When you square a QED amplitude and average over initial spins and sum over final ones, these collapse the spinor algebra into traces of gamma matrices — spinor bookkeeping becomes trace technology. Finally, the payoff that motivated the whole four-component apparatus: the Dirac density is manifestly non-negative, the exact defect the Klein–Gordon current suffered from.

jμ=ψˉγμψ,∂μjμ=0,j0=ψ†ψ=∑a=14∣ψa∣2≥0(8)j^\mu=\bar\psi\gamma^\mu\psi,\qquad \partial_\mu j^\mu=0,\qquad j^0=\psi^\dagger\psi=\sum_{a=1}^{4}|\psi_a|^2\ge 0 \qquad (8)

where: jμj^\mu — the conserved Dirac four-current; j0=ψ†ψj^0=\psi^\dagger\psi is a sum of squared moduli, so it can never be negative — the first-order structure of the Dirac equation buys a positive probability density where "Klein-Gordon Probability Current and Negative Densities" could not.

Derivation ABuilding the spinors at rest and boosting to reveal the small components

We construct the positive-energy spinor from scratch: solve the momentum-space equation (p̸−m)u=0(\not p-m)u=0 in the rest frame where it is trivial, then read off the general-momentum solution by keeping the full p̸\not p and solving block by block. The strategy is to split u=(ϕ,χ)⊤u=(\phi,\chi)^\top into upper (ϕ\phi) and lower (χ\chi) Pauli two-spinors and let the lower block of the matrix equation express χ\chi in terms of ϕ\phi; the watch-out is that the normalization must be fixed at the end so that u†u=2Eu^\dagger u=2E and uˉu=2m\bar u u=2m both come out right, and one must not confuse this large/small split with the particle/antiparticle split.

Step 1. In the rest frame pμ=(m,0)p^\mu=(m,\mathbf 0), so p̸=γ0m\not p=\gamma^0 m and the equation reads (γ0−1)m u=0(\gamma^0-1)m\,u=0. In the Dirac representation γ0=diag(1,−1)\gamma^0=\mathrm{diag}(\mathbb 1,-\mathbb 1), so this projects onto the +1+1 eigenspace: the upper two components. The rest-frame solution is u=2m (ξ,0)⊤u=\sqrt{2m}\,(\xi,0)^\top with ξ\xi any Pauli two-spinor. Why the upper block: γ0\gamma^0 has eigenvalue +1+1 on the top two slots and −1-1 on the bottom two, so demanding γ0u=u\gamma^0 u=u kills the lower block entirely at rest.

Step 2. For general momentum write p̸=γ0Ep−γ⋅p\not p = \gamma^0 E_p - \boldsymbol\gamma\cdot\mathbf p. With u=(ϕ,χ)⊤u=(\phi,\chi)^\top and γi=(0σi−σi0)\gamma^i=\begin{pmatrix}0&\sigma^i\\-\sigma^i&0\end{pmatrix}, the lower-block row of (p̸−m)u=0(\not p - m)u=0 becomes −σ⋅p ϕ−(Ep+m)χ=0-\boldsymbol\sigma\cdot\mathbf p\,\phi-(E_p+m)\chi=0.

Step 3. Solve for the small component: χ=σ⋅pEp+m ϕ\chi=\dfrac{\boldsymbol\sigma\cdot\mathbf p}{E_p+m}\,\phi. The lower components are locked to the upper ones; they vanish at rest and grow with momentum, with magnitude ratio ∣p∣/(Ep+m)|\mathbf p|/(E_p+m). Why "small": for ∣p∣≪m|\mathbf p|\ll m we have Ep+m≈2mE_p+m\approx 2m, so χ≈(σ⋅p/2m)ϕ\chi\approx (\boldsymbol\sigma\cdot\mathbf p/2m)\phi is suppressed by ∣p∣/2m≪1|\mathbf p|/2m\ll 1 — the nonrelativistic hierarchy.

Step 4. Fix the normalization. With ϕ=Ep+m ξ\phi=\sqrt{E_p+m}\,\xi one finds u†u=(Ep+m)[1+p2/(Ep+m)2]ξ†ξu^\dagger u=(E_p+m)\bigl[1+\mathbf p^2/(E_p+m)^2\bigr]\xi^\dagger\xi. Using Ep2−m2=p2E_p^2-m^2=\mathbf p^2 this simplifies to 2Ep2E_p, and the scalar to uˉu=(Ep+m)−p2/(Ep+m)=2m\bar u u=(E_p+m)-\mathbf p^2/(E_p+m)=2m, giving the boxed spinor of equation (4). Verified: at p→0\mathbf p\to 0 the spinor returns to 2m(ξ,0)⊤\sqrt{2m}(\xi,0)^\top; a symbolic check for pμ=(13,0,0,3), m=2p^\mu=(\sqrt{13},0,0,3),\,m=2 gives u†u=213=7.211u^\dagger u=2\sqrt{13}=7.211 and uˉu=4=2m\bar u u=4=2m exactly (Bjorken–Drell §3.2; Peskin–Schroeder §3.3).

Derivation BThe spin-sum completeness relation as a genuine matrix identity

We prove ∑susuˉs=p̸+m\sum_s u^s\bar u^s=\not p+m not by manipulating symbols but by building the explicit 4×4 dyadic and simplifying it with the Pauli-spinor completeness ∑sξsξs†=1\sum_s\xi^s\xi^{s\dagger}=\mathbb 1. The route is to form ∑susus†\sum_s u^s u^{s\dagger} block by block, collapse the spin sum in each block, and then right-multiply by γ0\gamma^0 to convert us†u^{s\dagger} into the Dirac adjoint uˉs\bar u^s; the payoff is the projector p̸+m\not p+m, and the trace of both sides gives an independent numerical cross-check.

Step 1. Form the outer-product sum. With us=(E+m ξs, σ⋅pE+mξs)⊤u^s=(\sqrt{E+m}\,\xi^s,\,\tfrac{\boldsymbol\sigma\cdot\mathbf p}{\sqrt{E+m}}\xi^s)^\top and ∑sξsξs†=1\sum_s\xi^s\xi^{s\dagger}=\mathbb 1, the upper-left block of ∑susus†\sum_s u^s u^{s\dagger} is (E+m)1(E+m)\mathbb 1, the lower-right is (σ⋅p)2E+m=p2E+m1=(E−m)1\tfrac{(\boldsymbol\sigma\cdot\mathbf p)^2}{E+m}=\tfrac{\mathbf p^2}{E+m}\mathbb 1=(E-m)\mathbb 1, and the off-diagonal blocks are σ⋅p\boldsymbol\sigma\cdot\mathbf p.

Step 2. Right-multiply by γ0=diag(1,−1)\gamma^0=\mathrm{diag}(\mathbb 1,-\mathbb 1) to build ∑susuˉs=(∑susus†)γ0\sum_s u^s\bar u^s=(\sum_s u^s u^{s\dagger})\gamma^0. This flips the sign of the two right-hand blocks, turning the lower-right into −(E−m)1=(m−E)1-(E-m)\mathbb 1=(m-E)\mathbb 1 and the lower-left into −σ⋅p-\boldsymbol\sigma\cdot\mathbf p.

Step 3. Assemble. The result is ((E+m)1−σ⋅pσ⋅p(m−E)1)\begin{pmatrix}(E+m)\mathbb 1 & -\boldsymbol\sigma\cdot\mathbf p\\ \boldsymbol\sigma\cdot\mathbf p & (m-E)\mathbb 1\end{pmatrix}, which is exactly p̸+m\not p+m in the Dirac representation, where p̸=γ0E−γ⋅p=(E 1−σ⋅pσ⋅p−E 1)\not p=\gamma^0 E-\boldsymbol\gamma\cdot\mathbf p=\begin{pmatrix}E\,\mathbb 1 & -\boldsymbol\sigma\cdot\mathbf p\\ \boldsymbol\sigma\cdot\mathbf p & -E\,\mathbb 1\end{pmatrix} and m=m 14m=m\,\mathbb 1_4. Why γ0\gamma^0 enters: the completeness relation is a statement about uuˉu\bar u, and uˉ=u†γ0\bar u=u^\dagger\gamma^0, so the γ0\gamma^0 is not optional — dropping it gives the (non-covariant) ∑uu†\sum u u^\dagger instead.

Step 4. Cross-check by trace: Tr(p̸+m)=Tr(p̸)+4m=0+4m=4m\mathrm{Tr}(\not p+m)=\mathrm{Tr}(\not p)+4m=0+4m=4m, and independently ∑suˉsus=∑s2m=4m\sum_s\bar u^s u^s=\sum_s 2m=4m. The two agree. Verified: for pμ=(13,0,0,3), m=2p^\mu=(\sqrt{13},0,0,3),\,m=2, direct 4×4 construction gives ∑susuˉs−(p̸+m)=0\sum_s u^s\bar u^s-(\not p+m)=0 to machine precision, and ∑svsvˉs=p̸−m\sum_s v^s\bar v^s=\not p-m likewise (Peskin–Schroeder Eq. 3.66–3.67).

Limiting cases

The spinor structure simplifies sharply in the natural kinematic regimes, and each limit isolates one physical idea.

CaseResultMeaning
Rest frame, p=0\mathbf p=0us=2m (ξs,0)⊤u^s=\sqrt{2m}\,(\xi^s,0)^\top, vs=2m (0,ηs)⊤v^s=\sqrt{2m}\,(0,\eta^s)^\topUpper and lower blocks decouple; ξs,ηs\xi^s,\eta^s are ordinary Pauli two-spinors. Particle sits "up," antiparticle "down."
Nonrelativistic, ∣p∣≪m|\mathbf p|\ll msmall block ≈σ⋅p2m \approx \tfrac{\boldsymbol\sigma\cdot\mathbf p}{2m}\,(large block)Lower components suppressed by ∣p∣/2m|\mathbf p|/2m; the large components obey the Pauli equation — nonrelativistic spin physics recovered.
Ultrarelativistic / massless, E≫mE\gg muˉu→0\bar u u\to 0; helicity =±=\pmchirality, PL,R=12(1∓γ5)P_{L,R}=\tfrac12(1\mp\gamma^5)The scalar normalization collapses; helicity becomes a good quantum number and left/right Weyl components decouple.
⚠️ The upper/lower blocks are NOT particle vs antiparticle

The two-component blocks of a Dirac spinor are the large and small components — a statement about the nonrelativistic limit, tied to the ratio ∣p∣/(E+m)|\mathbf p|/(E+m). They are not the electron and positron. A single electron of nonzero momentum has both blocks filled and is a pure uu-type solution. The particle/antiparticle distinction is the uu-versus-vv (energy-sign) split entirely. Conflating "small components" with "antiparticle content" is the single most common error in relativistic quantum mechanics, and it produces spurious factors and sign mistakes when you least expect them.

Worked examples

Example 1 — Construct an explicit spinor and check its density. Take a particle of mass m=2m=2 (natural units) moving along zz with pz=3p_z=3; build u1(p)u^1(p) and verify u†u=2Eu^\dagger u = 2E.

Step 1. On-shell energy: E=pz2+m2=9+4=13=3.6056E=\sqrt{p_z^2+m^2}=\sqrt{9+4}=\sqrt{13}=3.6056. The normalization factor is E+m=13+2=5.6056=2.3676\sqrt{E+m}=\sqrt{\sqrt{13}+2}=\sqrt{5.6056}=2.3676.

Step 2. Large (upper) block for spin-up ξ1=(1,0)⊤\xi^1=(1,0)^\top: E+m ξ1=(2.3676, 0)⊤\sqrt{E+m}\,\xi^1=(2.3676,\,0)^\top.

Step 3. Small (lower) block: σ⋅pE+mξ1\dfrac{\boldsymbol\sigma\cdot\mathbf p}{\sqrt{E+m}}\xi^1. Here σ⋅p=pzσz=3 σz\boldsymbol\sigma\cdot\mathbf p=p_z\sigma_z=3\,\sigma_z, and σzξ1=ξ1\sigma_z\xi^1=\xi^1, so the block is 32.3676(1,0)⊤=(1.2671, 0)⊤\dfrac{3}{2.3676}(1,0)^\top=(1.2671,\,0)^\top.

Step 4. Assemble: u1(p)=(2.3676, 0, 1.2671, 0)⊤u^1(p)=(2.3676,\,0,\,1.2671,\,0)^\top. Its density is u†u=2.36762+1.26712=5.6056+1.6055=7.2111u^\dagger u = 2.3676^2 + 1.2671^2 = 5.6056 + 1.6055 = 7.2111. The answer u1†u1=7.2111=213=2Epu^{1\dagger}u^1 = 7.2111 = 2\sqrt{13} = 2E_p, confirming the density normalization u†u=2Eu^\dagger u=2E exactly for a concrete spinor. The small-to-large ratio is 1.2671/2.3676=0.5352=∣p∣/(E+m)1.2671/2.3676 = 0.5352 = |\mathbf p|/(E+m), an order-one number — this particle is genuinely relativistic.

Example 2 — Verify the on-shell condition kills the spinor. For the same u1(p)u^1(p), check (p̸−m)u1=0(\not p - m)u^1 = 0.

Step 1. Build p̸=γ0E−γ3pz\not p = \gamma^0 E - \gamma^3 p_z. In block form p̸=(E 1−pzσzpzσz−E 1)\not p = \begin{pmatrix} E\,\mathbb 1 & -p_z\sigma_z\\ p_z\sigma_z & -E\,\mathbb 1\end{pmatrix} with E=13, pz=3E=\sqrt{13},\,p_z=3.

Step 2. Act on the large block. The upper two rows of (p̸−m)u1(\not p-m)u^1 give (E−m)⋅2.3676+(−pzσz)⋅(1.2671,0)⊤(E-m)\cdot 2.3676 + (-p_z\sigma_z)\cdot(1.2671,0)^\top. Numerically (E−m)=1.6056(E-m)=1.6056, so the first term is 1.6056×2.3676=3.80171.6056\times 2.3676=3.8017; the second is −3×1.2671=−3.8013-3\times 1.2671=-3.8013 (up to rounding). They cancel: 3.8017−3.8013≈03.8017-3.8013\approx 0.

Step 3. Act on the small block. The lower two rows give pzσz⋅(2.3676,0)⊤+(−E−m)⋅(1.2671,0)⊤=3×2.3676−5.6056×1.2671=7.1028−7.1028=0p_z\sigma_z\cdot(2.3676,0)^\top + (-E-m)\cdot(1.2671,0)^\top = 3\times 2.3676 - 5.6056\times 1.2671 = 7.1028 - 7.1028 = 0.

Step 4. Both blocks vanish, so (p̸−m)u1=0(\not p - m)u^1 = 0. The spinor sits exactly in the null space of p̸−m\not p - m — it is a genuine positive-energy solution. (The tiny residuals are decimal rounding; the exact symbolic check gives identically zero, verified in the accompanying script.)

The story behind the physics

The four-component spinor is the object Dirac was forced into when, in 1928, he demanded a wave equation first order in time whose square reproduced the relativistic energy relation. The two extra "negative-energy" solutions embarrassed him until, in 1930, he reinterpreted them as an infinite filled sea whose holes are antiparticles — a picture soon replaced by Feynman and Stückelberg's cleaner reading but vindicated in spirit when Anderson found the positron in a cloud chamber in 1932. The two-component Pauli spinors ξ\xi that sit inside the large block are the very objects Pauli had introduced in 1927 to describe nonrelativistic spin; Dirac's construction shows spin was relativity all along.

DS
dirac spinor
the four-component object of relativistic QM · 1928– · Cambridge · Zürich · Pasadena

The Dirac spinor is the four-component wavefunction on which the gamma matrices act. Its upper and lower blocks are not four particles but a spin-½ particle and its antiparticle partner; its structure encodes spin, the electron's magnetic moment, and — through the negative-energy branch reinterpreted via charge conjugation — the positron. It emerged from Paul Dirac's 1928 demand for a first-order relativistic equation, absorbed Pauli's 1927 two-component spin object into its large block, and was vindicated when Anderson's 1932 positron realized its negative-energy solutions.

⏳ History
1927
— Pauli introduces two-component spinors for nonrelativistic electron spin
1928
— Dirac writes the first-order relativistic equation, forcing four-component spinors
1930
— Dirac proposes the hole theory reinterpreting the negative-energy v-spinors
1932
— Anderson observes the positron in a cloud chamber, realizing the v-solutions
1948
— Feynman–Stückelberg reinterpretation recasts v-spinors as forward-in-time antiparticles

Where this lives today

Dirac spinors are the daily currency of high-energy physics. Every cross section computed at the LHC begins by writing external electrons, quarks, and their antiparticles as u(p)u(p) and v(p)v(p) spinors, then uses the spin sums ∑uuˉ=p̸+m\sum u\bar u=\not p+m to fold the spinor algebra into gamma-matrix traces that a computer algebra system evaluates automatically. The same spinors, quantized, become the mode functions of the electron field: the uu-spinors ride electron annihilation operators and the vv-spinors ride positron creation operators. In lattice QCD, in Standard Model event generators, and in searches for physics beyond it, the four-component plane-wave spinor built here is the irreducible starting point.

Going further

References. Griffiths, Introduction to Elementary Particles, 2nd ed., Ch. 7 · Sakurai & Napolitano, Modern Quantum Mechanics, Ch. 8 (Relativistic QM) · Shankar, Principles of Quantum Mechanics, 2nd ed., Ch. 20 · Peskin & Schroeder, An Introduction to Quantum Field Theory, §3.3 · Bjorken & Drell, Relativistic Quantum Mechanics, §3.2.

Provenance. All equations are standard relativistic QM. Independently verified in the accompanying script: (1) (p̸∓m)w=0(\not p\mp m)w=0 follows from the plane-wave ansatz and {γμ,γν}=2gμν\{\gamma^\mu,\gamma^\nu\}=2g^{\mu\nu}; (2) the Dirac-representation spinors give uˉu=2m\bar u u=2m, u†u=2Epu^\dagger u=2E_p for pμ=(13,0,0,3), m=2p^\mu=(\sqrt{13},0,0,3),\,m=2; (3) the spin sums ∑uuˉ=p̸+m\sum u\bar u=\not p+m, ∑vvˉ=p̸−m\sum v\bar v=\not p-m hold as exact 4×4 identities with Tr(p̸+m)=4m\mathrm{Tr}(\not p+m)=4m; (4) the projectors Λ±=(±p̸+m)/2m\Lambda_\pm=(\pm\not p+m)/2m are idempotent, complete, and orthogonal. Base repair: the source article stated the explicit spinors in the chiral (Weyl) basis while framing the derivation in the Dirac (standard) representation and quoting vˉv=−2m\bar v v=-2m — an internal basis inconsistency. This article pins everything to the Dirac representation of "The Klein-Gordon Equation," so every explicit component, block structure, and sign is self-consistent and machine-checked in one basis.

Builds on: The Dirac Equation · Dirac Gamma Matrices and Clifford Algebra · The Klein-Gordon Equation · Negative-Energy Solutions and the Interpretive Crisis. Leads to: Spin from the Dirac Equation · The Dirac Sea and Antiparticles · Charge Conjugation and the Positron · Bilinear Covariants and Dirac Adjoint · Nonrelativistic Reduction and the Pauli Equation.

Why does the plane-wave ansatz turn a differential equation into an algebraic one?

Because ∂μ\partial_\mu acting on e−ip⋅xe^{-ip\cdot x} just multiplies by −ipμ-ip_\mu. So iγμ∂μi\gamma^\mu\partial_\mu becomes γμpμ=p̸\gamma^\mu p_\mu=\not p, a constant matrix. The Dirac equation (iγμ∂μ−m)ψ=0(i\gamma^\mu\partial_\mu-m)\psi=0 then reduces to (p̸−m)u=0(\not p-m)u=0 — no derivatives left, just a 4×4 matrix acting on a constant column. This is the universal trick for constant-coefficient linear equations: plane waves diagonalize translations, converting calculus into linear algebra.

What fixes the on-shell condition p² = m²?

Multiply (p̸−m)u=0(\not p-m)u=0 on the left by (p̸+m)(\not p+m): you get (̸p2−m2)u=0(\not p^2-m^2)u=0. Since ̸p2=12{γμ,γν}pμpν=gμνpμpν=p2\not p^2=\tfrac12\{\gamma^\mu,\gamma^\nu\}p_\mu p_\nu=g^{\mu\nu}p_\mu p_\nu=p^2, this reads (p2−m2)u=0(p^2-m^2)u=0. For a nonzero spinor you need p2=m2p^2=m^2, i.e. E=±p2+m2E=\pm\sqrt{\mathbf p^2+m^2}. The mass shell is not an extra assumption; it is forced by the gamma algebra, and it is exactly why relativistic energy–momentum emerges automatically.

Why two spin states rather than one or four?

Once you are on shell, (p̸−m)(\not p-m) is a rank-2 matrix: its null space is two-dimensional. Equivalently, Λ+=(p̸+m)/2m\Lambda_+=(\not p+m)/2m is a projector with trace 2, so it projects onto a two-dimensional subspace. Those two dimensions are the two spin states. The other two dimensions of the four-component space are the null space of (p̸+m)(\not p+m) — the two vv-spinors. So 4=2 (spins)×2 (energy signs)4 = 2\,(\text{spins}) \times 2\,(\text{energy signs}).

Is the "small component" the same as the antiparticle part?

No — this is the central pitfall. The small (lower) component is present in every moving electron and is suppressed by ∣p∣/(E+m)|\mathbf p|/(E+m); it is a nonrelativistic-limit statement. The antiparticle is the entirely separate vv-spinor with the opposite energy sign. A pure electron is a uu-spinor with both large and small components; it contains zero antiparticle content. Watch the two decompositions carefully: upper/lower is large/small, while u/vu/v is particle/antiparticle.

Why does the density carry a factor 2E while the scalar carries 2m?

uˉu=u†γ0u\bar u u=u^\dagger\gamma^0 u is a Lorentz scalar, so it must be the same number in every frame — and at rest that number is 2m2m. The density u†uu^\dagger u is the time component of the four-vector current jμ=uˉγμuj^\mu=\bar u\gamma^\mu u, so it transforms and grows with energy as 2E2E. Physically, a relativistic normalization volume Lorentz-contracts by γ=E/m\gamma=E/m, so the density must scale up by E/mE/m to keep the total probability (density × volume) invariant. Hence 2m→2E2m\to 2E.

What happens to these spinors when m → 0?

The scalar normalization uˉu=2m→0\bar u u=2m\to 0, so it stops being useful; you switch to u†u=2Eu^\dagger u=2E. More importantly, chirality (eigenvalue of γ5\gamma^5) becomes conserved and coincides with helicity, so the left- and right-handed pieces PL,Ru=12(1∓γ5)uP_{L,R}u=\tfrac12(1\mp\gamma^5)u decouple and each satisfies a two-component Weyl equation. This is why massless (or ultrarelativistic) fermions are naturally described by helicity/chirality, and why the weak interaction — coupling only to left-handed fields — produces strongly polarized particles at high energy.

How do these plane waves become the quantized Dirac field?

You promote the expansion coefficients to operators. The free field is ψ(x)=∫d3p(2π)312Ep∑s(apsus(p)e−ip⋅x+bps†vs(p)e+ip⋅x)\psi(x)=\int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_p}}\sum_s(a^s_{\mathbf p}u^s(p)e^{-ip\cdot x}+b^{s\dagger}_{\mathbf p}v^s(p)e^{+ip\cdot x}). The uu-spinors multiply electron annihilation operators aa; the vv-spinors multiply positron creation operators b†b^\dagger. The plane-wave spinors built here are precisely the mode functions of second quantization, and the anticommutators of a,ba,b enforce Fermi statistics.

Before this

Where this leads