Skip to main content

Nuclear, particle & quantum-information sequence · Interactive lesson

Colour charge, SU(3)c and the QCD Lagrangian

Colour is the charge nobody has measured directly, and this is where you learn to trust it anyway: write the one line of Lagrangian that is the whole strong interaction, turn the Gell-Mann matrices into the three numbers every QCD amplitude reduces to, and watch three separate experiments count Nc = 3 without a colour meter.

UPV-18.07 · Lesson 7 of 9 in this unit. Position is not a completion score.
Starting model: Colour charge, SU(3)_c and the QCD LagrangianR = σ(e⁺e⁻→hadrons)/σ(e⁺e⁻→μ⁺μ⁻) against √s on a log axis, 2 to 50 GeV. Solid: N_c Σ Q_f² (1 + α_s/π) for N_c = 3 and α_s = 0.15, stepping up at the charm (3.7 GeV) and bottom (10.6 GeV) thresholds. Dashed: measured plateaus 2.2, 3.6 and 3.9. At √s = 30.0 GeV the model gives R = 3.84.solid: Nc Σ Qf² (1 + αₛ/π)dashed: measured R (PDG)R = 2R = 4R = 62charm 3.7bottom 10.650 GeV√s = 30.0 GeV: R = 3.84 → Nc = 3.05
Start with a picture. Change the model in the lab below.

Nuclear, particle & quantum-information sequence. These units continue into nuclear and particle physics and quantum information. Follow your university’s prerequisites, mathematical depth and assessment; this independent pathway is not a universal programme specification.

1

Build the model

Colour is a three-dimensional space, not three labels.

The key idea

Colour charge is the SU(3)c quantum number: each quark field is a three-component vector rotated by the 3×3 Gell-Mann generators, and the QCD Lagrangian is the renormalisable field theory that stays invariant under that rotation at every point of spacetime.

Example. The combination (λ⁴ + iλ⁵)/2 is the 3×3 matrix with a single 1 in row r, column b, so the gluon it multiplies inside Dμ turns a blue quark into a red one: that gluon carries colour r b̄, and there are six such colour-changing gluons plus two diagonal ones.

The complete physical model and assumptions

The strong interaction is one gauged SU(3) acting on a hidden index. Each quark flavour is a triplet q = (qᵣ, qg, qb) in the fundamental representation, and demanding that the Lagrangian survive an independent SU(3) rotation of that triplet at every point of spacetime forces the covariant derivative Dμ = ∂μ − i gₛ (λa/2) Gaμ, with eight gluon fields, one per Gell-Mann matrix. That is the entire model: L = Σf q̄f (iγμ Dμ − mf) qf − ¼ Gaμν Gaμν, a single coupling gₛ shared by all six flavours. Where it comes from is an antisymmetry problem — the Δ⁺⁺, three u quarks in a symmetric ground state, needs an extra index with at least three values to satisfy Pauli, and εᵢⱼₖ is the antisymmetric singlet that does it. What it costs is that colour is non-abelian. The eight charges Qa do not commute, so no quark has a definite value of all of them; gluons carry colour–anticolour and couple to one another through fabc; and every vertex is weighted by a colour factor — the Casimirs CF = 4/3 and CA = 3, or the index TF = ½ — rather than by a charge you can read off. Because only colour singlets propagate as far as a detector, Nc = 3 is never observed. It is inferred, from the R ratio, from the π⁰ → γγ rate and from the multiplicity of gluon jets: three measurements that agree on the same integer.

1Why three: the Δ⁺⁺ and the antisymmetric singlet

Colour was not observed; it was demanded. The Δ⁺⁺ is three u quarks with J = 3/2 in the ground state, so its spatial wavefunction is symmetric, its spin state ↑↑↑ is symmetric and its flavour content uuu is symmetric — a totally symmetric state of three identical fermions, which Pauli forbids. The repair is a hidden index i on each quark and a totally antisymmetric combination of it, εᵢⱼₖ qᵢ qⱼ qₖ/√6. That combination exists only if the index takes at least three values, and with exactly three it is unique. So quarks live in a three-dimensional internal space, and the group that rotates a complex triplet while preserving that antisymmetric singlet is SU(3): quarks in the fundamental 3, antiquarks in the conjugate 3̄, whose generators are −(λa)*/2. Everything else in the topic is what that one choice forces once the symmetry is gauged.

2Write the Lagrangian and count the gluons off the generators

Gauge the triplet and the covariant derivative is fixed: (Dμ)ᵢⱼ = δᵢⱼ ∂μ − i gₛ (λa/2)ᵢⱼ Gaμ, one gluon field per generator, so L = Σf q̄f (iγμDμ − mf) qf − ¼ Gaμν Gaμν with a single coupling for all six flavours. Why eight gluons and not nine? A gluon carries a colour and an anticolour, and 3 ⊗ 3̄ = 8 ⊕ 1. The six off-diagonal combinations r b̄, r ḡ, g r̄, g b̄, b r̄, b ḡ change a quark's colour; two diagonal ones, (r r̄ − g ḡ)/√2 and (r r̄ + g ḡ − 2b b̄)/√6, are λ³ and λ⁸. The ninth, (r r̄ + g ḡ + b b̄)/√3, is the identity — the extra generator of U(3). It commutes with the other eight, has no self-coupling, and would be a massless colour-singlet vector coupling to every quark alike: an unconfined, infinite-range strong force between nucleons, which nobody sees. So the group is SU(3), the gluons are eight, and the inputs are gₛ and six masses, plus a θ-term that the neutron's electric dipole moment bounds below 10⁻¹⁰.

3Trace identities turn colour into three numbers

Nobody sums eight 3×3 matrices inside a calculation; they use the Casimirs. Normalise tr(TaTb) = TF δab with TF = ½. Since Σₐ TaTa commutes with every generator, Schur's lemma makes it a multiple of the identity, CF 𝟙; trace both sides and (N² − 1) TF = N CF, so CF = (N² − 1)/2N = 8/6 = 4/3. Check it on a red quark directly: the four off-diagonal generators that touch row 1 — T¹, T², T⁴, T⁵ — each contribute |½|² = ¼, a total of 1; the diagonal ones add (½)² + (1/2√3)² = ¼ + 1/12 = 1/3; the sum is 4/3. In the adjoint the same construction reads facd fbcd = CA δab, and for a = b = 1 the six non-zero f¹cd (±1, ±½, ±½) square to 2(1 + ¼ + ¼) = 3, so CA = N = 3. These are the weights at every vertex once colour is summed: CF when a quark radiates a gluon, CA when a gluon radiates a gluon, TF when a gluon splits into q q̄.

4One gluon exchange: the channels that attract build singlets

One gluon exchanged between two colour charges gives a Coulomb-like potential V(r) = ⟨T₁⋅T₂⟩ αₛ ħc/r, and the colour factor comes from completing the square: (T₁ + T₂)² = T₁² + T₂² + 2T₁⋅T₂, so ⟨T₁⋅T₂⟩ = ½[C(R) − C(R₁) − C(R₂)] for a pair coupled to representation R. For q q̄ the singlet gives ½[0 − 8/3] = −4/3 and the octet ½[3 − 8/3] = +1/6; for q q the antisymmetric 3̄ gives ½[4/3 − 8/3] = −2/3 and the symmetric 6, with C = 10/3, gives +1/3. Put numbers on it: at r = 0.50 fm with αₛ = 0.30, αₛ ħc/r = 0.30 × 197.3/0.50 = 118 MeV, so the q q̄ singlet sits at −158 MeV, the octet at +20 MeV, the q q antitriplet at −79 MeV and the sextet at +39 MeV. The pattern is the point. The only attractive channels are the ones that build colour singlets: a q q̄ pair in a 1 is a meson, and two quarks in a 3̄ look like an antiquark to the third quark, which is how three quarks close into a singlet baryon. Octets and sextets repel and never bind — confinement's shape, before confinement is even invoked.

5Gluons carry colour, so fabc puts them in their own current

In QED the photon is neutral and the field strength is linear in A. In QCD the field strength Gaμν = ∂μGaν − ∂νGaμ + gₛ fabc Gbμ Gcν carries a term with no derivative at all, and −¼G² then holds vertices with three gluons (∝ gₛ fabc) and with four (∝ gₛ² fabe fcde). Physically, a gluon labelled r b̄ is itself a colour charge, so it radiates gluons exactly as a quark does — but with weight CA = 3 instead of CF = 4/3. That ratio, 9/4, is measurable: a gluon jet should end in more hadrons than a quark jet of the same energy, the multiplicity ratio tending to CA/CF = 2.25 at asymptotic energies and sitting near 1.5 at LEP's 91 GeV once finite-energy corrections are included. The same self-coupling is what makes the one-loop coefficient β₀ = (11/3)CA − (2/3)nf positive: the +11 exists only because fabc ≠ 0. An abelian theory of colour would have no gluon jets to compare, no asymptotic freedom, and no reason for quarks to be confined.

6Nc = 3 is inferred, three different ways

No detector registers colour, because only singlets reach it. What it registers are rates, and colour multiplies them. In e⁺e⁻ → hadrons the photon makes a q q̄ pair in any of Nc colours, so R = σ(hadrons)/σ(μ⁺μ⁻) = Nc Σ Qf²: with u, d, s that is Nc × 2/3, above the charm threshold Nc × 10/9, above bottom Nc × 11/9. For Nc = 3 the plateaus are 2, 3.33 and 3.67 before the (1 + αₛ/π) correction, and the data sit near 2.2, 3.6 and 3.9; for Nc = 1 they would be 0.67, 1.11 and 1.22. The π⁰ → γγ decay runs through a quark triangle whose amplitude is proportional to Nc(Qᵤ² − Qd²) = Nc/3, so the rate goes as (Nc/3)²: 7.8 eV predicted at Nc = 3 against 7.8 ± 0.1 eV measured, and a factor of nine too small at Nc = 1. And the τ decays to hadrons Nc times as often as to each lepton, giving Rτ = Nc (1 + αₛ/π + …) ≈ 3.6 against a measured 3.64. Three integers, one answer, none of them read off a colour meter.

2

Key equations

Read the units, assumptions and limits before substituting.

The QCD Lagrangian

L = Σf q̄fi (iγμ Dμ − mf)ᵢⱼ qfj − ¼ Gaμν Gaμν

Seven inputs — gₛ and six masses — and one line: everything from proton structure to jets is this Lagrangian plus gauge fixing.

i, j = r, g, b colour indices; f runs over u, d, s, c, b, t; a = 1…8; gₛ dimensionless, mf in GeV

Covariant derivative and generators

(Dμ)ᵢⱼ = δᵢⱼ ∂μ − i gₛ (λa/2)ᵢⱼ Gaμ, tr(λa λb) = 2δab

A 3×3 matrix acting on colour: the six off-diagonal λa swap a quark's colour when a gluon is emitted; λ³ and λ⁸ leave it alone and weigh it.

λa: the eight Gell-Mann matrices, 3×3, Hermitian, traceless; Ta = λa/2 obey [Ta, Tb] = i fabc Tc

Gluon field strength and self-coupling

Gaμν = ∂μ Gaν − ∂ν Gaμ + gₛ fabc Gbμ Gcν

Squared inside −¼G², the f-term gives three- and four-gluon vertices: gluons are coloured and couple directly to one another, which photons never do.

f¹²³ = 1; f¹⁴⁷ = f²⁴⁶ = f²⁵⁷ = f³⁴⁵ = ½; f¹⁵⁶ = f³⁶⁷ = −½; f⁴⁵⁸ = f⁶⁷⁸ = √3/2; all others zero up to antisymmetry

Casimirs from trace identities

tr(TaTb) = TF δabΣₐ TaTa = CF 𝟙facd fbcd = CA δab

Every colour-summed amplitude reduces to these: a quark radiates a gluon with weight CF, a gluon radiates a gluon with CA, and a gluon splits to q q̄ with TF.

SU(N): TF = ½, CF = (N² − 1)/2N, CA = N. For N = 3 that is ½, 4/3 and 3 — pure numbers, no units

One-gluon-exchange colour factor

V(r) = ⟨T₁⋅T₂⟩ αₛ ħc/r, ⟨T₁⋅T₂⟩ = ½[C(R) − C(R₁) − C(R₂)]

The attractive channels are exactly the singlet-building ones — why mesons are q q̄ singlets, and why two quarks in a 3̄ pull in a third.

q q̄: singlet −4/3, octet +1/6; q q: antitriplet −2/3, sextet +1/3; αₛ = gₛ²/4π, ħc = 197.3 MeV fm

Counting colours from measurement

R = Nc Σf Qf² (1 + αₛ/π + …)Γ(π⁰→γγ) = (Nc/3)² α² mπ³/(64π³ fπ²)

R steps 2 → 10/3 → 11/3 across the c and b thresholds only for Nc = 3, and the π⁰ rate would be nine times too small without it.

flavours with 2mf < √s only, Qf in units of e; α = 1/137.04, fπ = 92.1 MeV, mπ = 135.0 MeV
3

Interactive lab

Predict first. Explore the model. Compare the result.

Explore the model

3
1 to 5; steps of 1.
0.15
0 to 0.5; steps of 0.05.
30.0 GeV
2 to 50 GeV; steps of 0.5.

Try this. Leave αₛ at 0.15 and step Nc from 1 to 5: only Nc = 3 drops the top plateau onto the 3.9 line. Then slide √s down to 3 GeV and watch the implied Nc stay near 3 while that plateau wants αₛ nearer 0.3 — the coupling runs, and that is the next topic.

Live results

Σ Qf² BELOW √s
1.222
R PREDICTED
3.84
R MEASURED PLATEAU
3.9
Nc IMPLIED BY DATA
3.05

Results follow this lesson’s model and assumptions. Displayed values are rounded.

Interactive Physics V modelR = σ(e⁺e⁻→hadrons)/σ(e⁺e⁻→μ⁺μ⁻) against √s on a log axis, 2 to 50 GeV. Solid: N_c Σ Q_f² (1 + α_s/π) for N_c = 3 and α_s = 0.15, stepping up at the charm (3.7 GeV) and bottom (10.6 GeV) thresholds. Dashed: measured plateaus 2.2, 3.6 and 3.9. At √s = 30.0 GeV the model gives R = 3.84.solid: Nc Σ Qf² (1 + αₛ/π)dashed: measured R (PDG)R = 2R = 4R = 62charm 3.7bottom 10.650 GeV√s = 30.0 GeV: R = 3.84 → Nc = 3.05
R = σ(e⁺e⁻→hadrons)/σ(e⁺e⁻→μ⁺μ⁻) against √s on a log axis, 2 to 50 GeV. Solid: Nc Σ Qf² (1 + αₛ/π) for Nc = 3 and αₛ = 0.15, stepping up at the charm (3.7 GeV) and bottom (10.6 GeV) thresholds. Dashed: measured plateaus 2.2, 3.6 and 3.9. At √s = 30.0 GeV the model gives R = 3.84.

Predict → change one variable → check

Plan a change before running it. Every other input stays at the saved baseline.

Saved baseline. Colour count Nc: 3; Strong coupling αₛ: 0.15; Collision energy √s: 30.0 GeV

Baseline 3; choose 1 to 5, in steps of 1.
4. Compared with the saved baseline, this result will…

This experiment is kept only while this page is open. It does not change saved course progress.

4

Catch the common trap

The tempting picture is three copies of QED: red, green and blue as three conserved charges, each carried the way an electron carries −e, with gluons as neutral photons between them. It fails for one reason: the generators do not commute. Because [Ta, Tb] = i fabc Tc, the eight colour charges Qa cannot be diagonalised together — only T³ and T⁸ can — so no quark has a definite value of all eight, and 'red' is a basis vector in a three-dimensional space, not a measurable number. Gluon emission rotates that vector: absorbing an r b̄ gluon turns a blue quark red, which no photon does. And gluons sit in the adjoint, carrying colour–anticolour, so they radiate more strongly than quarks (CA = 3 against CF = 4/3) and couple to each other through fabc. Three commuting charges would give three neutral photons, no self-coupling, no asymptotic freedom and no confinement — not QCD.

5

Practice & worked examples

3 authored examples. Make an attempt, then reveal the reasoning one step at a time.

EasyExample 1

Using the Gell-Mann matrices with Ta = λa/2, compute Σₐ (TaTa)ᵣᵣ for a red quark (colour index 1) and confirm that it equals CF = (N² − 1)/2N for N = 3.

Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.

0 of 4 steps revealed.

MediumExample 2

Predict R = σ(e⁺e⁻ → hadrons)/σ(e⁺e⁻ → μ⁺μ⁻) at √s = 3.0 GeV, 7.0 GeV and 35 GeV for Nc = 3, including the leading QCD correction (1 + αₛ/π) with αₛ = 0.30, 0.20 and 0.14 respectively. Compare with the measured plateaus 2.2, 3.6 and 3.9, and state what Nc = 1 would have given.

Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.

0 of 5 steps revealed.

HardExample 3

The π⁰ → γγ rate through the quark triangle is Γ = (Nc/3)² α² mπ³/(64π³ fπ²), the Nc/3 arising as Nc(Qᵤ² − Qd²). Evaluate it for Nc = 3 with α = 1/137.04, mπ = 134.98 MeV and fπ = 92.1 MeV, convert to a lifetime using the branching fraction B(γγ) = 98.8%, compare with the measured τ = 8.43 × 10⁻¹⁷ s, and state the Nc = 1 result.

Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.

0 of 5 steps revealed.

6

Quick check

Select an answer and explain it before checking.

QCD has eight gluons rather than the nine colour–anticolour combinations a quark triplet can form. Which statement gives the actual reason?
7

Review & continue

Colour charge is the SU(3)c quantum number: each quark field is a three-component vector rotated by the 3×3 Gell-Mann generators, and the QCD Lagrangian is the renormalisable field theory that stays invariant under that rotation at every point of spacetime.

  • Explain the relationship without looking at the equation.
  • Check a sign, a unit and a model assumption in your working.
  • Return to the lab or worked example if a step is still unclear.

Extra challenging questions

Connect several ideas, explore supplied data, and compare your reasoning with a worked solution.