Dirac adjoint & bilinear covariants

Write down the Dirac equation and you immediately own a four-component object ψ\psi — but four components is not the same as a four-vector, and that is the trap the whole chapter turns on. Spinors rotate into themselves under a representation S(Λ)S(\Lambda) of the Lorentz group that is not unitary: boosts are generated by anti-Hermitian matrices, so S†S≠1S^\dagger S \ne 1. The consequence is blunt — the naïve probability density ψ†ψ\psi^\dagger\psi, which served Schrödinger theory faithfully, is no longer a Lorentz scalar. It is the time slot of something larger, and it changes when you run past the particle. The fix is one factor of γ0\gamma^0. Define the Dirac adjoint ψˉ≡ψ†γ0\bar\psi \equiv \psi^\dagger\gamma^0, and the metric-aware γ0\gamma^0 exactly undoes the non-unitarity of SS: where ψ→Sψ\psi \to S\psi, the adjoint runs the other way, ψˉ→ψˉ S−1\bar\psi \to \bar\psi\,S^{-1}. Sandwich any fixed matrix Γ\Gamma between ψˉ\bar\psi and ψ\psi and the two SS's meet in the middle, leaving a clean tensor law. These sandwiches are the bilinear covariants — exactly sixteen of them, a scalar, a pseudoscalar, a four-vector, an axial vector and an antisymmetric tensor. They are the alphabet from which every Lorentz-invariant fermion coupling is spelled: the mass term mψˉψm\bar\psi\psi, the electromagnetic current ψˉγμψ\bar\psi\gamma^\mu\psi, the weak V ⁣− ⁣AV\!-\!A current ψˉγμ(1−γ5)ψ\bar\psi\gamma^\mu(1-\gamma^5)\psi, the Pauli moment ψˉσμνψ Fμν\bar\psi\sigma^{\mu\nu}\psi\,F_{\mu\nu}. Get the adjoint right and relativistic covariance takes care of itself.

Key takeaways
  • The Dirac adjoint is ψˉ≡ψ†γ0\bar\psi \equiv \psi^\dagger\gamma^0; the extra γ0\gamma^0 is what makes spinor products transform as Lorentz tensors.
  • The spinor representation is non-unitary (S†S≠1S^\dagger S\ne1 for boosts), so ψ†ψ\psi^\dagger\psi is not invariant — it is the time component of the four-vector current ψˉγμψ\bar\psi\gamma^\mu\psi.
  • The master identity γ0S†γ0=S−1\gamma^0 S^\dagger\gamma^0 = S^{-1} gives ψˉ→ψˉ S−1\bar\psi\to\bar\psi\,S^{-1}, so ψˉψ\bar\psi\psi is a scalar and ψˉγμψ\bar\psi\gamma^\mu\psi a four-vector.
  • There are exactly 1616 bilinear covariants: 11 scalar + 11 pseudoscalar + 44 vector + 44 axial + 66 tensor, a complete basis of 4×44\times4 Dirac matrices.
  • The conserved current jμ=ψˉγμψj^\mu=\bar\psi\gamma^\mu\psi satisfies ∂μjμ=0\partial_\mu j^\mu=0; its five siblings build every fermion interaction in the Standard Model.
  • Key idea: spinors do not transform like vectors, so the naïve product ψ†ψ\psi^\dagger\psi has no clean Lorentz behaviour. Inserting the metric-aware factor γ0\gamma^0 — using ψˉ=ψ†γ0\bar\psi=\psi^\dagger\gamma^0 — exactly compensates for the spinor transformation, turning ψˉ Γ ψ\bar\psi\,\Gamma\,\psi into objects with definite tensor character.
The core insight

Spinors do not transform like vectors, so naïve ψ†ψ\psi^\dagger\psi-type products have no clean Lorentz behaviour. Inserting the metric-aware factor γ0\gamma^0 — using ψˉ=ψ†γ0\bar\psi=\psi^\dagger\gamma^0 instead of ψ†\psi^\dagger — exactly compensates for the spinor transformation, turning every ψˉ Γ ψ\bar\psi\,\Gamma\,\psi into an object with a definite Lorentz-tensor transformation law.

Where this lives today

The bilinear covariants are the literal grammar of every calculation in particle physics. In any Standard Model Lagrangian each fermion interaction is one of the five families: the QED vertex is the vector current ψˉγμψ\bar\psi\gamma^\mu\psi coupled to the photon, the mass and Higgs–Yukawa terms are the scalar ψˉψ\bar\psi\psi, and the charged-current weak interaction is the parity-violating V ⁣− ⁣AV\!-\!A form ψˉγμ(1−γ5)ψ\bar\psi\gamma^\mu(1-\gamma^5)\psi coupled to the WW boson — tested to sub-percent precision at the LHC and in beta-decay experiments. The tensor ψˉσμνψ\bar\psi\sigma^{\mu\nu}\psi has a vivid modern life: as the Pauli term ψˉσμνψ Fμν\bar\psi\sigma^{\mu\nu}\psi\,F_{\mu\nu} it encodes the anomalous magnetic moment, which the Fermilab Muon g ⁣− ⁣2g\!-\!2 experiment (final result 2023) measures for the muon to about 0.20.2 ppm, where a tension with theory hunts for new physics. Pseudoscalar bilinears ψˉγ5ψ\bar\psi\gamma^5\psi govern pion couplings and drive electric-dipole-moment searches — a nonzero neutron or electron EDM would signal CP-violating scalar–pseudoscalar interference. Even lattice QCD extracts hadron structure by evaluating these five currents between quark fields. The adjoint Dirac introduced to make a probability come out positive is the alphabet in which the whole matter sector is written.

The physical picture

The whole subject sits inside one thought experiment: take an electron at rest, then run past it. The density ψ†ψ\psi^\dagger\psi is a charge density, and density is length-contracted — sprint past the particle and you sweep up more per unit volume, so ψ†ψ\psi^\dagger\psi grows. That growth is the signature that ψ†ψ\psi^\dagger\psi was never a scalar: it is the time component of the current jμ=ψˉγμψj^\mu=\bar\psi\gamma^\mu\psi, and like any four-vector its pieces reshuffle under a boost while one combination — its Minkowski length — stays nailed down. That fixed length is what the scalar ψˉψ\bar\psi\psi measures. The first figure follows a real electron through a boost.

Figure
ψ̄ψ is the same number in every frame while ψ†ψ = 2E is the time component of the current four-vector jᵘ = ū γᵘ u = 2pᵘ; boosting an electron from rest E = m = 0.511 MeV to E = 2m = 1.022 MeV grows u†u from 1.022 to 2.044 MeV but leaves ū u = 2m = 1.022 MeV fixed
ψ̄ψ is invariant; ψ†ψ is the time-slot of the current four-vector jᵘ = 2pᵘ rest frame p = 0, E = m = 0.511 MeV boosted E = 2m = 1.022 MeV space time j⁰ = u†u = 2m = 1.022 length = ū u = 2m = 1.022 space j⁰ = u†u = 2E = 2.044 j³ = 2p = 1.771 invariant length = ū u = 2m = 1.022

Why sixteen and not more? Because ψˉ Γ ψ\bar\psi\,\Gamma\,\psi is built from one 4×44\times4 matrix Γ\Gamma, and the space of 4×44\times4 matrices is sixteen-dimensional. A chosen basis of those matrices — the identity, γ5\gamma^5, the four γμ\gamma^\mu, the four γμγ5\gamma^\mu\gamma^5, the six σμν\sigma^{\mu\nu} — sorts the sixteen sandwiches into five Lorentz families that never mix under proper transformations. The second figure maps that alphabet, worth memorising because every fermion coupling you will write is one row of it.

Figure
The complete set of 16 bilinear covariants ψ̄Γψ splits into five Lorentz types — scalar 1, pseudoscalar 1, vector 4, axial-vector 4, antisymmetric tensor 6 — which sum to 16 and exhaust the space of 4×4 Dirac matrices; the last column shows the sign each picks up under a parity flip
The 16 bilinear covariants: 1 + 1 + 4 + 4 + 6 = 16 type Γ components parity scalar ψ̄ψ 1 + pseudoscalar ψ̄γ⁵ψ 1 − vector ψ̄γᵘψ 4 V⁰ + , Vⁱ − axial vector ψ̄γᵘγ⁵ψ 4 A⁰ − , Aⁱ + tensor ψ̄σᵘᵛψ 6 antisymmetric total 16 σᵘᵛ = (i/2)[γᵘ , γᵛ] , γ⁵ = i γ⁰γ¹γ²γ³ dim of 4×4 complex matrices = 16 — the covariants are a complete basis

The core equations

Everything begins with a definition. The Dirac adjoint is the row spinor you get by taking the Hermitian conjugate and reinserting one factor of γ0\gamma^0; that factor is the entire trick, because γ0\gamma^0 is the object that knows the Minkowski metric.

ψˉ≡ψ†γ0(1)\bar\psi \equiv \psi^\dagger\gamma^0 \qquad (1)

where: ψ\psi — the four-component Dirac spinor (a column); ψ†\psi^\dagger — its Hermitian conjugate (a row of complex conjugates); γ0\gamma^0 — the timelike Dirac matrix, Hermitian with (γ0)2=1(\gamma^0)^2 = 1; ψˉ\bar\psi — the Dirac adjoint, a row spinor. Read it as "conjugate-transpose, then reflect through the time axis." The reason γ0\gamma^0 is the right insertion is a two-line matrix identity and how it propagates to the finite transformation S(Λ)S(\Lambda). The hermiticity conditions γ0†=γ0\gamma^{0\dagger}=\gamma^0, γk†=−γk\gamma^{k\dagger}=-\gamma^k package into one covariant statement that forces the adjoint to transform inversely to ψ\psi.

(γμ)†=γ0γμγ0,γ0S†γ0=S−1    ⟹    ψˉ→ψˉ S−1(2)(\gamma^\mu)^\dagger = \gamma^0\gamma^\mu\gamma^0, \qquad \gamma^0 S^\dagger\gamma^0 = S^{-1} \;\;\Longrightarrow\;\; \bar\psi \to \bar\psi\,S^{-1} \qquad (2)

where: S=S(Λ)S=S(\Lambda) — the 4×44\times4 spinor representation of a Lorentz transformation, satisfying S−1γμS=ΛμνγνS^{-1}\gamma^\mu S = \Lambda^\mu{}_\nu\gamma^\nu, with S†S≠1S^\dagger S\ne1 for boosts. The middle identity is the workhorse — γ0\gamma^0 turns a dagger into an inverse: since ψ→Sψ\psi\to S\psi, ψˉ=ψ†γ0→ψ†S†γ0=ψ†γ0S−1=ψˉ S−1\bar\psi=\psi^\dagger\gamma^0\to\psi^\dagger S^\dagger\gamma^0=\psi^\dagger\gamma^0 S^{-1}=\bar\psi\,S^{-1}. With the adjoint transforming as ψˉ S−1\bar\psi\,S^{-1}, any fixed matrix between ψˉ\bar\psi and ψ\psi yields a definite tensor. A complete basis of 4×44\times4 matrices sorts the sixteen sandwiches into five families.

ψˉψ,ψˉγ5ψ,ψˉγμψ,ψˉγμγ5ψ,ψˉσμνψ(3)\bar\psi\psi,\quad \bar\psi\gamma^5\psi,\quad \bar\psi\gamma^\mu\psi,\quad \bar\psi\gamma^\mu\gamma^5\psi,\quad \bar\psi\sigma^{\mu\nu}\psi \qquad (3)

where: ψˉψ\bar\psi\psi — scalar (11 component); ψˉγ5ψ\bar\psi\gamma^5\psi — pseudoscalar (11), with γ5≡iγ0γ1γ2γ3\gamma^5\equiv i\gamma^0\gamma^1\gamma^2\gamma^3, (γ5)2=1(\gamma^5)^2=1, {γ5,γμ}=0\{\gamma^5,\gamma^\mu\}=0; ψˉγμψ\bar\psi\gamma^\mu\psi — vector (44); ψˉγμγ5ψ\bar\psi\gamma^\mu\gamma^5\psi — axial vector (44); ψˉσμνψ\bar\psi\sigma^{\mu\nu}\psi — antisymmetric tensor (66), with σμν≡i2[γμ,γν]\sigma^{\mu\nu}\equiv\tfrac{i}{2}[\gamma^\mu,\gamma^\nu]. The component counts 1+1+4+4+6=161+1+4+4+6=16 exhaust the sixteen dimensions of the Dirac-matrix space. The most important covariant is the vector one: the conserved probability-and-charge current, whose conservation is now manifestly Lorentz-invariant. The scalar covariant supplies the mass term of the Dirac Lagrangian.

jμ=ψˉγμψ,∂μjμ=0,L=ψˉ (iγμ∂μ−m) ψ(4)j^\mu=\bar\psi\gamma^\mu\psi,\qquad \partial_\mu j^\mu=0,\qquad \mathcal{L}=\bar\psi\,(i\gamma^\mu\partial_\mu-m)\,\psi \qquad (4)

where: jμj^\mu — the conserved four-current, with density j0=ψ†ψ≥0j^0=\psi^\dagger\psi\ge0 and flux j=ψ†αψ\mathbf{j}=\psi^\dagger\boldsymbol\alpha\psi; ∂μjμ=0\partial_\mu j^\mu=0 — the continuity equation, a scalar and hence frame-independent; L\mathcal{L} — the Dirac Lagrangian, whose mass term −mψˉψ-m\bar\psi\psi is a scalar. Coupling jμj^\mu to the photon, −e ψˉγμψ Aμ-e\,\bar\psi\gamma^\mu\psi\,A_\mu, is the QED interaction.

Derivation AWhy the adjoint carries a factor of gamma-zero

The goal is to show, with no hand-waving, that ψ†ψ\psi^\dagger\psi fails to be a Lorentz scalar while ψˉψ=ψ†γ0ψ\bar\psi\psi=\psi^\dagger\gamma^0\psi succeeds, and that the difference is forced by one algebraic fact — the spinor representation is non-unitary. We start from the spinor transformation law, establish the master identity γ0S†γ0=S−1\gamma^0 S^\dagger\gamma^0=S^{-1} from the hermiticity of the gammas, and feed it through the two simplest sandwiches to read off a scalar and a four-vector.

Step 1. Under xμ→Λμνxνx^\mu\to\Lambda^\mu{}_\nu x^\nu the spinor transforms as ψ→S(Λ)ψ\psi\to S(\Lambda)\psi, with SS fixed (up to sign) by demanding the Dirac equation keep its form, which requires S−1γμS=ΛμνγνS^{-1}\gamma^\mu S=\Lambda^\mu{}_\nu\gamma^\nu. For a rotation the generator is Hermitian and SS unitary, but for a boost it is anti-Hermitian, so S†=SS^\dagger=S and S†S=S2≠1S^\dagger S=S^2\ne1. Why: this non-unitarity is the whole reason ψ†ψ\psi^\dagger\psi misbehaves — a unitary SS would leave ψ†ψ\psi^\dagger\psi invariant, but boosts are not unitary in this finite-dimensional representation.

Step 2. Establish the master identity. The conditions γ0†=γ0\gamma^{0\dagger}=\gamma^0, γk†=−γk\gamma^{k\dagger}=-\gamma^k combine into (γμ)†=γ0γμγ0(\gamma^\mu)^\dagger=\gamma^0\gamma^\mu\gamma^0. Since SS is built by exponentiating σμν\sigma^{\mu\nu} and (σμν)†=γ0σμνγ0(\sigma^{\mu\nu})^\dagger=\gamma^0\sigma^{\mu\nu}\gamma^0 follows from the same relation, exponentiation gives γ0S†γ0=S−1\gamma^0 S^\dagger\gamma^0=S^{-1}, i.e. S†γ0=γ0S−1S^\dagger\gamma^0=\gamma^0 S^{-1}. Why: this is the exact statement that sandwiching with γ0\gamma^0 converts a Hermitian conjugate into an inverse — precisely what cancels the offending S†S^\dagger.

Step 3. Transform the adjoint. Using Step 2, ψˉ=ψ†γ0→ψ†S†γ0=ψ†γ0S−1=ψˉ S−1\bar\psi=\psi^\dagger\gamma^0\to\psi^\dagger S^\dagger\gamma^0=\psi^\dagger\gamma^0 S^{-1}=\bar\psi\,S^{-1}, whereas the bare conjugate ψ†→ψ†S†\psi^\dagger\to\psi^\dagger S^\dagger keeps the uncancelled non-unitary factor. Why: the adjoint runs inversely to ψ\psi, so the two transformations annihilate each other — the single fact that makes every bilinear covariant work.

Step 4. Read off two covariants. The scalar sandwich gives ψˉψ→ψˉ S−1S ψ=ψˉψ\bar\psi\psi\to\bar\psi\,S^{-1}S\,\psi=\bar\psi\psi — invariant. The vector sandwich gives ψˉγμψ→ψˉ S−1γμS ψ=Λμν ψˉγνψ\bar\psi\gamma^\mu\psi\to\bar\psi\,S^{-1}\gamma^\mu S\,\psi=\Lambda^\mu{}_\nu\,\bar\psi\gamma^\nu\psi — a genuine four-vector, using S−1γμS=ΛμνγνS^{-1}\gamma^\mu S=\Lambda^\mu{}_\nu\gamma^\nu from Step 1. Meanwhile ψ†ψ→ψ†S†Sψ≠ψ†ψ\psi^\dagger\psi\to\psi^\dagger S^\dagger S\psi\ne\psi^\dagger\psi. Why: the same S−1γμSS^{-1}\gamma^\mu S that makes the Dirac equation covariant promotes the current to a four-vector — covariance of the equation and tensor character of the bilinears are one fact seen twice.

Go deeper · Derivation B · 1 minThe conserved probability current
Derivation BThe conserved probability current

The aim is to prove that jμ=ψˉγμψj^\mu=\bar\psi\gamma^\mu\psi is conserved, ∂μjμ=0\partial_\mu j^\mu=0, working directly from the Dirac equation and its adjoint, and to see that conservation now holds as a Lorentz-scalar statement rather than the frame-dependent one of Schrödinger theory. We derive the equation obeyed by ψˉ\bar\psi, differentiate the current, and watch the mass terms cancel — which is what lets j0=ψ†ψj^0=\psi^\dagger\psi serve as a positive-definite density with a consistent flux in every frame.

Step 1. Start from (iγμ∂μ−m)ψ=0(i\gamma^\mu\partial_\mu-m)\psi=0 and take the Hermitian conjugate: −i ∂μψ†(γμ)†−mψ†=0-i\,\partial_\mu\psi^\dagger(\gamma^\mu)^\dagger-m\psi^\dagger=0. Why: to build ∂μ(ψˉγμψ)\partial_\mu(\bar\psi\gamma^\mu\psi) we need an equation of motion for ψˉ\bar\psi, obtained by conjugating the equation for ψ\psi.

Step 2. Right-multiply by γ0\gamma^0 and use (γμ)†γ0=γ0γμ(\gamma^\mu)^\dagger\gamma^0=\gamma^0\gamma^\mu (a rearrangement of the master identity), turning ψ†(γμ)†γ0\psi^\dagger(\gamma^\mu)^\dagger\gamma^0 into ψˉγμ\bar\psi\gamma^\mu. This gives the adjoint Dirac equation i ∂μψˉ γμ+mψˉ=0i\,\partial_\mu\bar\psi\,\gamma^\mu+m\bar\psi=0, i.e. ∂μψˉ γμ=imψˉ\partial_\mu\bar\psi\,\gamma^\mu=im\bar\psi. Why: the factor of γ0\gamma^0 converts the awkward (γμ)†(\gamma^\mu)^\dagger into a clean γμ\gamma^\mu acting from the right — again the adjoint is engineered to make the algebra close.

Step 3. Differentiate the current and substitute both equations of motion: ∂μjμ=(∂μψˉ)γμψ+ψˉγμ(∂μψ)=(imψˉ)ψ+ψˉ(−imψ)=0\partial_\mu j^\mu=(\partial_\mu\bar\psi)\gamma^\mu\psi+\bar\psi\gamma^\mu(\partial_\mu\psi)=(im\bar\psi)\psi+\bar\psi(-im\psi)=0. Why: the two mass terms are equal and opposite by construction, so the current is conserved for any solution — and because jμj^\mu is a four-vector, ∂μjμ\partial_\mu j^\mu is a scalar, so conservation in one frame holds in all frames.

Go deeper · Derivation C · 1 minSixteen matrices, sixteen covariants
Derivation CSixteen matrices, sixteen covariants

The goal here is to explain the number sixteen — why there are exactly that many independent bilinear covariants, why they organise into the five families scalar, pseudoscalar, vector, axial and tensor, and how a parity flip is the quickest way to tell the "true" tensors from the "pseudo" ones. We count the dimension of the matrix space, exhibit a basis, use trace-orthogonality to prove independence, and act with parity to split each pair — the backbone guaranteeing the alphabet is complete and non-redundant.

Step 1. A bilinear ψˉ Γ ψ\bar\psi\,\Gamma\,\psi is linear in a single 4×44\times4 matrix Γ\Gamma, and that complex vector space has dimension 1616, so there can be at most sixteen independent covariants — a basis of the matrix space yields a basis of covariants. Why: the counting is pure linear algebra; the physics only enters in choosing a basis whose members transform simply.

Step 2. Take the sixteen matrices { I, γ5, γμ, γμγ5, σμν }\{\,I,\ \gamma^5,\ \gamma^\mu,\ \gamma^\mu\gamma^5,\ \sigma^{\mu\nu}\,\} (μ<ν\mu<\nu in the tensor), counts 1+1+4+4+6=161+1+4+4+6=16. They are linearly independent because trace-orthogonal: Tr(ΓAΓB)=4 δAB\mathrm{Tr}(\Gamma^A\Gamma_B)=4\,\delta^A{}_B, so no nontrivial combination vanishes. Why: trace-orthogonality is the fermionic analogue of an orthonormal basis — any 4×44\times4 matrix expands as M=14∑ATr(ΓAM) ΓAM=\tfrac14\sum_A \mathrm{Tr}(\Gamma^A M)\,\Gamma_A, confirming completeness.

Step 3. Assign Lorentz character via S−1ΓSS^{-1}\Gamma S: the identity gives a scalar, γμ\gamma^\mu a vector, σμν\sigma^{\mu\nu} an antisymmetric rank-2 tensor; the two "pseudo" families differ from their partners by one γ5\gamma^5. Why: γ5\gamma^5 commutes with SS (a scalar under proper transformations), so multiplying by it does not change boost/rotation behaviour — the distinction only appears under parity.

Step 4. Act with parity P:ψ(t,x)→γ0ψ(t,−x)P:\psi(t,\mathbf{x})\to\gamma^0\psi(t,-\mathbf{x}). Then ψˉψ→+ψˉψ\bar\psi\psi\to+\bar\psi\psi (scalar) but ψˉγ5ψ→−ψˉγ5ψ\bar\psi\gamma^5\psi\to-\bar\psi\gamma^5\psi (pseudoscalar), since γ0γ5γ0=−γ5\gamma^0\gamma^5\gamma^0=-\gamma^5. The vector current has j0→+j0, ji→−jij^0\to+j^0,\ j^i\to-j^i (polar) while the axial current flips oppositely, A0→−A0, Ai→+AiA^0\to-A^0,\ A^i\to+A^i. Why: parity separates the five families into their observable roles — scalar mass term, pseudoscalar pion coupling, polar electromagnetic current, and the axial current whose interference with the vector current is the fingerprint of weak parity violation.

Limiting cases

CaseResultMeaning
Rest frame (p=0\mathbf p=0)ψ†ψ=ψˉψ=2m\psi^\dagger\psi=\bar\psi\psi=2mdensity and scalar coincide; the current has no spatial part
General boostψ†ψ=2E\psi^\dagger\psi=2E, ψˉψ=2m\bar\psi\psi=2mdensity grows as EE; the scalar is pinned at 2m2m — the Minkowski length of jμ=2pμj^\mu=2p^\mu
Massless / ultra-relativisticψˉψ, ψˉγ5ψ→0\bar\psi\psi,\ \bar\psi\gamma^5\psi\to0scalar and pseudoscalar vanish; only VV, AA, TT survive and chirality becomes exact
Non-relativistic limitjμ→(ρ, ρv)j^\mu\to(\rho,\ \rho\mathbf v)upper two components dominate; the current reduces to the Schrödinger probability current
Parity flip PPS +S\,+, P −P\,-, VV polar, AA axialthe five families sort by their sign under x→−x\mathbf x\to-\mathbf x
Charge conjugation CCψˉγμψ→−ψˉγμψ\bar\psi\gamma^\mu\psi\to-\bar\psi\gamma^\mu\psithe vector current flips sign — particle and antiparticle carry opposite charge

Worked examples

1 · The two invariants of a plane-wave electron (fully numerical)

Take an electron of mass m=0.511 MeVm=0.511\ \mathrm{MeV} prepared in a momentum eigenstate with E=2mE=2m. Using the Dirac-representation spinor us(p)=E+m (ξs, σ⋅pE+mξs) ⁣⊤u^s(p)=\sqrt{E+m}\,\big(\xi^s,\ \tfrac{\boldsymbol\sigma\cdot\mathbf p}{E+m}\xi^s\big)^{\!\top} with ξ†ξ=1\xi^\dagger\xi=1 and the normalisation u†u=2Eu^\dagger u=2E, evaluate u†uu^\dagger u, uˉu\bar u u, and the current uˉγμu\bar u\gamma^\mu u. Confirm which one is frame-independent.

Step 1. Kinematics. With E=2mE=2m, ∣p∣=E2−m2=3 m=1.732×0.511=0.885 MeV|\mathbf p|=\sqrt{E^2-m^2}=\sqrt{3}\,m=1.732\times0.511=\mathbf{0.885\ MeV}. Point p\mathbf p along zz, so (σ⋅p)2=∣p∣2=E2−m2(\boldsymbol\sigma\cdot\mathbf p)^2=|\mathbf p|^2=E^2-m^2.

Step 2. The density u†uu^\dagger u. Squaring, u†u=(E+m)[1+E2−m2(E+m)2]=(E+m)+(E−m)=2Eu^\dagger u=(E+m)\big[1+\tfrac{E^2-m^2}{(E+m)^2}\big]=(E+m)+(E-m)=2E. Numerically u†u=2E=2.044 MeVu^\dagger u=2E=\mathbf{2.044\ MeV}. This is j0j^0 — a density, depending on the frame through EE.

Step 3. The scalar uˉu\bar u u. Insert γ0=diag(I,−I)\gamma^0=\mathrm{diag}(I,-I), flipping the sign of the lower components: uˉu=u†γ0u=(E+m)[1−E2−m2(E+m)2]=(E+m)−(E−m)=2m\bar u u=u^\dagger\gamma^0 u=(E+m)\big[1-\tfrac{E^2-m^2}{(E+m)^2}\big]=(E+m)-(E-m)=2m. Numerically uˉu=2m=1.022 MeV\bar u u=2m=\mathbf{1.022\ MeV} — and crucially it carries no EE, so it is the same in every frame. The scalar uˉu=2m\bar u u=2m is frame-independent; the density u†u=2Eu^\dagger u=2E is not.

Step 4. The current uˉγμu\bar u\gamma^\mu u. A short calculation (the Gordon identity for equal spins) gives uˉγμu=2pμ\bar u\gamma^\mu u=2p^\mu. The time slot is uˉγ0u=u†u=2E=2p0\bar u\gamma^0 u=u^\dagger u=2E=2p^0; the spatial slot is uˉγ3u=u†α3u=2∣p∣=2×0.885=1.771 MeV=2p3\bar u\gamma^3 u=u^\dagger\alpha^3 u=2|\mathbf p|=2\times0.885=\mathbf{1.771\ MeV}=2p^3. So jμ=(2.044, 0, 0, 1.771) MeVj^\mu=(2.044,\,0,\,0,\,1.771)\ \mathrm{MeV}, with Minkowski length jμjμ=2.0442−1.7712=1.041=1.020≈2m\sqrt{j^\mu j_\mu}=\sqrt{2.044^2-1.771^2}=\sqrt{1.041}=1.020\approx2m. The vector's components shuffle with the frame, but its invariant length equals 2m2m — the very number the scalar uˉu\bar u u reports.

2 · Scalar versus pseudoscalar under a mirror (sign-tracking)

Show that ψˉψ\bar\psi\psi and ψˉγ5ψ\bar\psi\gamma^5\psi behave oppositely under a parity flip, and explain why this forces the pion and a hypothetical scalar meson to couple to nucleons through different bilinears.

Step 1. Parity acts as P:ψ(t,x)→γ0ψ(t,−x)P:\psi(t,\mathbf x)\to\gamma^0\psi(t,-\mathbf x), so the adjoint becomes ψˉ ′=(γ0ψ)†γ0=ψ†γ0†γ0=ψˉ γ0\bar\psi\,'=(\gamma^0\psi)^\dagger\gamma^0=\psi^\dagger\gamma^{0\dagger}\gamma^0=\bar\psi\,\gamma^0 (at −x-\mathbf x). Keep the γ0\gamma^0's explicit through the sandwiches.

Step 2. The scalar. ψˉψ→ψˉγ0⋅γ0ψ=ψˉ(γ0)2ψ=ψˉψ\bar\psi\psi\to\bar\psi\gamma^0\cdot\gamma^0\psi=\bar\psi(\gamma^0)^2\psi=\bar\psi\psi (at −x-\mathbf x). Even under parity — a true scalar.

Step 3. The pseudoscalar. ψˉγ5ψ→ψˉγ0γ5γ0ψ\bar\psi\gamma^5\psi\to\bar\psi\gamma^0\gamma^5\gamma^0\psi. Because {γ5,γ0}=0\{\gamma^5,\gamma^0\}=0, we have γ0γ5γ0=−γ5(γ0)2=−γ5\gamma^0\gamma^5\gamma^0=-\gamma^5(\gamma^0)^2=-\gamma^5, so ψˉγ5ψ→−ψˉγ5ψ\bar\psi\gamma^5\psi\to-\bar\psi\gamma^5\psi. Odd under parity — a pseudoscalar.

Step 4. Physics. A pion has parity JP=0−J^P=0^-, so a parity-conserving pion–nucleon coupling must pair the pion with a pseudoscalar bilinear, g ψˉγ5ψ ϕπg\,\bar\psi\gamma^5\psi\,\phi_\pi; a genuine scalar meson (0+0^+) couples through g ψˉψ ϕg\,\bar\psi\psi\,\phi. The one factor of γ5\gamma^5 is the difference between a 0+0^+ and a 0−0^- coupling — the bilinears carry the quantum numbers. The same VV-versus-AA bookkeeping applied to the currents makes the weak ψˉγμ(1−γ5)ψ\bar\psi\gamma^\mu(1-\gamma^5)\psi violate parity maximally.

Physics of the portrait

The adjoint spinor and its bilinears are so woven into field theory that it is easy to forget they were one person's inventions, arrived at in a rush between 1928 and the early 1930s while the meaning of a relativistic wave equation was still molten.

Origins · 1 minPaul Dirac and how this was discovered
PD
Paul Dirac
1902–1984 · relativistic quantum mechanics and the equation of the electron

Paul Adrien Maurice Dirac, trained as an electrical engineer in Bristol before turning to physics, wrote down in 1928 the relativistic wave equation that bears his name — and was forced to invent nearly everything in this article at once. Demanding a first-order equation consistent with special relativity produced four-component wavefunctions and the four anticommuting matrices γμ\gamma^\mu; insisting the theory carry a positive probability density forced the adjoint ψˉ=ψ†γ0\bar\psi=\psi^\dagger\gamma^0 and the current ψˉγμψ\bar\psi\gamma^\mu\psi that reduces to Schrödinger's in the slow limit. Dirac's temperament — spare, guided by an almost aesthetic conviction that fundamental equations must be beautiful — is stamped on the covariant form (iγμ∂μ−m)ψ=0(i\gamma^\mu\partial_\mu-m)\psi=0, which hides a mountain of structure behind one line. The "extra" negative-energy solutions troubled him until, in 1930–31, he reinterpreted them as antiparticles and predicted the positron, found by Anderson in 1932; he shared the 1933 Nobel Prize with Schrödinger. The classification of the sixteen bilinear covariants grew directly out of his formalism, and when parity violation surfaced in the 1950s it was the vector–axial distinction among his bilinears that let Feynman, Gell-Mann, Sudarshan and Marshak write down the weak interaction. Every fermion line in every Feynman diagram is, at bottom, a statement about Dirac adjoints.

History
1928
— Dirac writes the relativistic electron equation; the four-current ψ̄γᵘψ appears in the same paper
1929
— the adjoint spinor ψ̄ = ψ†γ⁰ and the covariant bilinears are formalized by Dirac and Darwin
1937
— Fierz derives the identities relating the five families of bilinear covariants
1957
— Wu's cobalt-60 experiment reveals parity violation, physically splitting the vector from the axial current
1958
— Feynman–Gell-Mann and Sudarshan–Marshak fix the weak current as the V minus A form ψ̄γᵘ(1−γ⁵)ψ

Going further

Texts. Peskin & Schroeder, An Introduction to Quantum Field Theory, §3.2–3.4 (spinor bilinears, discrete symmetries) · Griffiths, Introduction to Elementary Particles, 2nd ed., §7.3 (bilinear covariants) · Halzen & Martin, Quarks and Leptons, §5.4–5.6 (the sixteen Γ\Gamma-matrices and their transformations) · Thomson, Modern Particle Physics, §4.5–4.7 · Schwartz, Quantum Field Theory and the Standard Model, Ch. 10–11 (spinors and Lorentz representations); Bjorken & Drell, Relativistic Quantum Mechanics, §3–4.

Provenance. The covariant Dirac equation and γ0≡β\gamma^0\equiv\beta, γk≡βαk\gamma^k\equiv\beta\alpha^k are recalled from aq-pp-dirac-equation and aq-pp-gamma-matrices-clifford-algebra; the density–current pair ρ=ψ†ψ\rho=\psi^\dagger\psi, j=ψ†αψ\mathbf j=\psi^\dagger\boldsymbol\alpha\psi from aq-pp-dirac-probability-current; the index handling of Λμν\Lambda^\mu{}_\nu from aq-pp-four-vector-notation. A verifier should confirm: (1) γ0S†γ0=S−1\gamma^0 S^\dagger\gamma^0=S^{-1} from (γμ)†=γ0γμγ0(\gamma^\mu)^\dagger=\gamma^0\gamma^\mu\gamma^0; (2) ∂μ(ψˉγμψ)=0\partial_\mu(\bar\psi\gamma^\mu\psi)=0; (3) 1+1+4+4+6=161+1+4+4+6=16 and Tr(ΓAΓB)=4δAB\mathrm{Tr}(\Gamma^A\Gamma_B)=4\delta^A{}_B; (4) the numerics u†u=2E=2.044u^\dagger u=2E=2.044, uˉu=2m=1.022\bar u u=2m=1.022, uˉγμu=2pμ\bar u\gamma^\mu u=2p^\mu for E=2mE=2m, m=0.511 MeVm=0.511\ \mathrm{MeV}.

Prerequisites. The Dirac Equation · Gamma Matrices and the Clifford Algebra · The Dirac Probability Current · Four-Vector Notation. Leads to. The Covariant Dirac Equation · Feynman Rules for QED · Trace Techniques · Chirality and γ5\gamma^5.

Ask the tutor: Ask me to boost a spinor and watch ψ†ψ\psi^\dagger\psi change while ψˉψ\bar\psi\psi stays fixed, or to build any Standard-Model coupling from the five bilinear covariants.

Why not just use ψ† instead of ψ̄ = ψ†γ⁰?

Because ψ†\psi^\dagger transforms with the wrong factor. Under a boost ψ→Sψ\psi\to S\psi with SS non-unitary, so ψ†→ψ†S†\psi^\dagger\to\psi^\dagger S^\dagger and the product ψ†S†Sψ\psi^\dagger S^\dagger S\psi carries an uncancelled S†S≠1S^\dagger S\ne1. Inserting γ0\gamma^0 exploits S†γ0=γ0S−1S^\dagger\gamma^0=\gamma^0 S^{-1}, so the adjoint runs as ψˉ→ψˉ S−1\bar\psi\to\bar\psi\,S^{-1} and the S−1S^{-1} cancels the SS from ψ\psi. In short, ψ†ψ\psi^\dagger\psi is the time component of a four-vector (it grows under a boost), while ψˉψ\bar\psi\psi is a true scalar. The γ0\gamma^0 is not decoration; it is the metric doing its job in spinor space.

Are the gamma matrices a four-vector because they carry an index μ?

No — this is the classic trap. The four γμ\gamma^\mu are fixed, frame-independent constant matrices, the same numbers in every inertial frame. The relation S−1γμS=ΛμνγνS^{-1}\gamma^\mu S=\Lambda^\mu{}_\nu\gamma^\nu looks like a transformation law but states the opposite: it is a consistency condition among constant matrices. The genuine four-vector is the number-valued bilinear jμ=ψˉγμψj^\mu=\bar\psi\gamma^\mu\psi; its components mix under Λ\Lambda because ψ\psi and ψˉ\bar\psi transform, not because the gammas do. Think of γμ\gamma^\mu as fixed "machinery" and ψˉ(⋯ )ψ\bar\psi(\cdots)\psi as the tensor it manufactures.

Where does the number sixteen come from?

From pure linear algebra. A bilinear ψˉ Γ ψ\bar\psi\,\Gamma\,\psi is linear in one 4×44\times4 matrix Γ\Gamma, and the space of 4×44\times4 complex matrices has exactly 1616 dimensions. Any basis of that space gives a basis of bilinears; the useful basis {I,γ5,γμ,γμγ5,σμν}\{I,\gamma^5,\gamma^\mu,\gamma^\mu\gamma^5,\sigma^{\mu\nu}\} has 1+1+4+4+6=161+1+4+4+6=16 members that sort into five clean Lorentz families. Independence follows from trace-orthogonality, Tr(ΓAΓB)=4δAB\mathrm{Tr}(\Gamma^A\Gamma_B)=4\delta^A{}_B. So "sixteen" is not deep physics — it is the dimension of the matrix algebra — but the sorting into scalar, pseudoscalar, vector, axial and tensor is.

What is the difference between a scalar and a pseudoscalar, physically?

Both are single numbers invariant under rotations and boosts, so under proper Lorentz transformations they are indistinguishable. The difference appears only under a mirror reflection (parity): a scalar ψˉψ\bar\psi\psi is unchanged, a pseudoscalar ψˉγ5ψ\bar\psi\gamma^5\psi flips sign. Physically this labels the intrinsic parity of the particle a field couples to: a 0+0^+ scalar meson couples through ψˉψ\bar\psi\psi, a 0−0^- pion through ψˉγ5ψ\bar\psi\gamma^5\psi. In a parity-conserving theory you must match the parity of the coupling to that of the particle — mixing them would let a mirror experiment tell left from right for free.

What does the axial vector current have to do with the weak force?

Everything. The charged weak current is the combination ψˉγμ(1−γ5)ψ\bar\psi\gamma^\mu(1-\gamma^5)\psi — a vector minus an axial vector, the famous V ⁣− ⁣AV\!-\!A structure. Because the vector part flips one way under parity and the axial the opposite way, their sum has no definite parity, so the weak interaction violates parity — maximally, since the two enter with equal weight. This is why only left-handed fermions (and right-handed antifermions) feel the charged weak force. Once Wu's 1957 experiment showed parity is violated, the question became which mix of Dirac's bilinears nature uses; the answer, V ⁣− ⁣AV\!-\!A, was pinned down within a year.

Is ψ̄ψ the same thing as the probability density?

No, and conflating them is a common slip. The probability (and charge) density is j0=ψ†ψj^0=\psi^\dagger\psi, the time component of jμ=ψˉγμψj^\mu=\bar\psi\gamma^\mu\psi; it is positive-definite and frame-dependent, growing under a boost like any density. The scalar ψˉψ=ψ†γ0ψ\bar\psi\psi=\psi^\dagger\gamma^0\psi is a different object — Lorentz-invariant, equal to 2m2m for a plane wave in any frame, and not sign-definite (negative for antiparticle solutions). They coincide numerically only in the rest frame. In the Lagrangian ψ†ψ\psi^\dagger\psi never appears alone precisely because it is not invariant, whereas ψˉψ\bar\psi\psi is the mass term.

How does the tensor bilinear σᵘᵛ connect to spin and magnetism?

The tensor covariant ψˉσμνψ\bar\psi\sigma^{\mu\nu}\psi, with σμν=i2[γμ,γν]\sigma^{\mu\nu}=\tfrac{i}{2}[\gamma^\mu,\gamma^\nu], is the relativistic spin operator in disguise: its spatial components reduce to the Pauli spin matrices non-relativistically. Coupled to the field strength it forms the Pauli moment term ψˉσμνψ Fμν\bar\psi\sigma^{\mu\nu}\psi\,F_{\mu\nu}, whose coefficient is the magnetic moment. The Dirac equation alone predicts g=2g=2; radiative corrections add a small σμν\sigma^{\mu\nu} coupling, the anomalous moment a=(g−2)/2a=(g-2)/2, which the Fermilab Muon g ⁣− ⁣2g\!-\!2 experiment measures to sub-ppm precision. The sixth family of bilinears is where spin, magnetism, and precision QED all meet.

Before this

Where this leads