Quantum Chromodynamics — The Violent Heart of Matter

Inside every proton, three quarks are bound together by a force so strong that separating them completely requires infinite energy.

This is not a metaphor, and not an approximation!

If you try to pull a quark out of a proton, the energy you put in creates new quark-antiquark pairs before the quark escapes.

You cannot free a quark. You can only create more bound quarks. Boom! :)

This is confinement, one of the most remarkable phenomena in physics, and it is the reason you(!) exist.

The protons and neutrons in your body are made of quarks. Yes! Those quarks are permanently imprisoned by the strong nuclear force, Quantum Chromodynamics (QCD).

QCD is the most mathematically sophisticated and the least experimentally accessible of the fundamental theories.

It is fully solved in the ultraviolet: at short distances, perturbation theory works beautifully.

It is unsolved in the infrared: at long distances, where confinement occurs, perturbation theory breaks down completely.

The theory is known. Its consequences, at the scales that matter most for nuclear physics, are still being computed.

This is the story of that theory.

I. The Quark Model and Color

Before QCD, there was the quark model.

In the early 1960s, the “particle zoo” had become a zoo in the pejorative sense: dozens of hadrons with no organizing principle.

Murray Gell-Mann and George Zweig (1964) independently proposed: hadrons are made of more fundamental constituents, quarks.

Three flavors: up (u), down (d), strange (s).

Mesons: $q\bar{q}$ pairs. Baryons: $qqq$ triplets.

The quark model organized the hadrons into elegant multiplets based on $SU(3)$ flavor symmetry. It predicted new particles, most famously the $\Omega^-$ baryon, which were subsequently found.

The problem: statistics.

The $\Delta^{++}$ baryon has:

Three up quarks
Spin 3/2 (all spins aligned)
Symmetric spatial wavefunction

The total wavefunction is symmetric under exchange of any two quarks.

But quarks are fermions: their wavefunction must be antisymmetric.

Contradiction.

The solution: color.

Oscar Greenberg (1964) proposed: quarks carry a new quantum number, color, with three values: red, green, blue.

The color wavefunction of baryons is antisymmetric:

\[|\Delta^{++}\rangle = \frac{1}{\sqrt{6}} \epsilon_{abc}|u^a u^b u^c\rangle\]

where $\epsilon_{abc}$ is the antisymmetric tensor.

The total wavefunction: (symmetric) × (symmetric) × (antisymmetric) = antisymmetric ✓

Color saved the quark model.

But color raised a question: if quarks have color, what is the force associated with color?

The answer is QCD.

II. The QCD Lagrangian

QCD is a non-Abelian gauge theory based on the gauge group $SU(3)_C$, the color gauge group.

The quark fields:

Quarks come in 6 flavors $f = u, d, s, c, b, t$ and 3 colors $a = 1, 2, 3$.

Each quark is a Dirac spinor: $\psi^a_f(x)$

transforming in the fundamental representation of $SU(3)_C$:

\[\psi^a \to U^a{}_b \psi^b, \quad U \in SU(3)\]

The gluon fields:

The gauge bosons of $SU(3)_C$ are gluons, 8 gluons corresponding to the 8 generators of $SU(3)$.

The gluon field: $A^A_\mu(x)$, $A = 1,...,8$

In matrix form: $A_\mu = A^A_\mu T^A$

where $T^A = \lambda^A/2$ are the generators ($\lambda^A$: Gell-Mann matrices).

The field strength tensor:

\[F^A_{\mu\nu} = \partial_\mu A^A_\nu - \partial_\nu A^A_\mu + g_s f^{ABC}A^B_\mu A^C_\nu\]

The crucial difference from QED: the $f^{ABC}A^B_\mu A^C_\nu$ term. Gluons carry color charge and therefore interact with each other.

This non-Abelian self-interaction is the source of QCD’s most remarkable properties.

The QCD Lagrangian:

\[\mathcal{L}_{QCD} = \sum_f \bar{\psi}^a_f (i\gamma^\mu D_\mu - m_f)\psi_{fa} - \frac{1}{4}F^A_{\mu\nu}F^{A\mu\nu}\]

where the covariant derivative:

\[D_\mu\psi^a = \partial_\mu\psi^a - ig_s A^A_\mu(T^A)^a{}_b\psi^b\]

This is one of the most compact and beautiful descriptions of a fundamental force ever written.

III. Asymptotic Freedom

In 1973, David Gross, Frank Wilczek, and David Politzer computed the one-loop beta function of QCD.

Their result changed physics.

The beta function:

\[\beta(g_s) = \mu\frac{dg_s}{d\mu} = -\frac{g_s^3}{16\pi^2} \left(\frac{11}{3}N_c - \frac{2}{3}N_f\right)\]

For $SU(3)$ with $N_c = 3$ colors and $N_f$ quark flavors:

\[\beta(g_s) = -\frac{g_s^3}{16\pi^2} \left(11 - \frac{2N_f}{3}\right)\]

For $N_f \leq 16$: $\beta < 0$.

With 6 quark flavors: $\beta(g_s) = -\frac{7g_s^3}{16\pi^2} < 0$

The running coupling:

\[\alpha_s(\mu) = \frac{g_s^2}{4\pi} = \frac{\alpha_s(\mu_0)} {1 + \frac{\alpha_s(\mu_0)}{2\pi} \left(\frac{11N_c - 2N_f}{3}\right) \ln(\mu/\mu_0)}\]

At high energy $\mu \to \infty$: $\alpha_s \to 0$

Asymptotic freedom:

The strong coupling becomes weak at high energies. Quarks inside a proton, probed at short distances, behave as nearly free particles.

This explains the Bjorken scaling observed in deep inelastic scattering: the quarks inside the proton appear as point-like free particles at high momentum transfer.

Gross, Wilczek, and Politzer received the Nobel Prize in 2004.

The two-loop running:

\[\alpha_s(\mu) = \frac{2\pi}{\beta_0\ln(\mu/\Lambda_{QCD})} \left[1 - \frac{\beta_1\ln\ln(\mu/\Lambda_{QCD})} {\beta_0^2\ln(\mu/\Lambda_{QCD})} + ...\right]\]

where: $\beta_0 = \frac{11N_c - 2N_f}{3}$ $\beta_1 = \frac{17N_c^2 - 5N_cN_f - 3C_F N_f}{3}$

The QCD scale: $\Lambda_{QCD} \approx 200-300 \text{ MeV}$

This is the fundamental scale of QCD, where the coupling becomes of order unity and perturbation theory breaks down.

At $\mu = M_Z = 91.2$ GeV: $\alpha_s(M_Z) \approx 0.118$

At $\mu = 1$ GeV: $\alpha_s \approx 0.5$

At $\mu \sim \Lambda_{QCD}$: $\alpha_s \sim 1$, and perturbation theory fails.

IV. Color Confinement

At low energies, at the scales of hadron physics, $\alpha_s \sim 1$.

Perturbation theory fails. The strongly coupled regime begins. And confinement emerges.

The linear potential:

Lattice QCD simulations show: the potential between a quark and antiquark at large separation $r$ grows linearly:

\[V(r) \xrightarrow{r\to\infty} \sigma r\]

where $\sigma \approx (420 \text{ MeV})^2 \approx 0.18 \text{ GeV}^2$ is the string tension.

This linear potential implies: separating the quarks requires energy proportional to the separation. The energy density is constant, like an elastic string.

To separate quarks to infinity: infinite energy required. Confinement.

The flux tube:

The color electric field between a quark and antiquark does not spread out as in QED.

It is collimated into a flux tube, a thin string of color field.

Diameter: $\sim 1$ fm $= 10^{-15}$ m String tension: $\sigma \approx 0.18$ GeV²

The flux tube is the physical realization of the linear potential.

Why does the flux tube form?

In QED: the photon is electrically neutral. Electric field lines spread out, the Coulomb potential.

In QCD: gluons carry color charge. They interact with each other. They attract each other.

The color field is squeezed into a tube by the mutual attraction of gluons.

This is the mechanism of confinement, qualitatively understood, not yet rigorously proven from first principles.

The Millennium Problem:

The Clay Mathematics Institute lists as one of its seven Millennium Problems:

Prove that Yang-Mills theory exists and has a mass gap.

The mass gap, the energy difference between the vacuum and the first excited state, is equivalent to confinement.

If QCD has a mass gap, isolated color charges are impossible.

This has not been proven mathematically.

A proof would be worth $1,000,000.

It would also be the most important result in mathematical physics since the 20th century.

V. The QCD Vacuum: Instantons and the $\theta$ Term

The QCD vacuum is extraordinarily rich.

Topological structure:

The gauge field configurations of QCD can be classified by a topological invariant, the winding number or topological charge:

\[Q = \frac{g_s^2}{32\pi^2} \int d^4x \, F^A_{\mu\nu}\tilde{F}^{A\mu\nu}\]

where $\tilde{F}^{A\mu\nu} = \frac{1}{2}\epsilon^{\mu\nu\rho\sigma}F^A_{\rho\sigma}$ is the dual field strength.

$Q \in \mathbb{Z}$, an integer.

The different topological sectors, labeled by $Q$, are separated by energy barriers.

Instantons:

An instanton is a classical solution to the Euclidean equations of motion with $Q = \pm 1$.

The BPST instanton (Belavin, Polyakov, Schwarz, Tyupkin, 1975):

\[A^a_\mu = \frac{2\eta^a_{\mu\nu}(x-x_0)_\nu} {(x-x_0)^2 + \rho^2}\]

where $\eta^a_{\mu\nu}$ is the ‘t Hooft symbol and $\rho$ is the instanton size.

The instanton represents a tunneling event, the gauge field tunneling between topologically inequivalent vacua.

The $\theta$ vacuum:

The true QCD vacuum is a superposition of all topological sectors:

\[|\theta\rangle = \sum_{n=-\infty}^{\infty} e^{in\theta}|n\rangle\]

This is the $\theta$ vacuum.

The $\theta$ parameter is a new free parameter of QCD, the coefficient of the topological term in the Lagrangian:

\[\mathcal{L}_\theta = \frac{\theta g_s^2}{32\pi^2} F^A_{\mu\nu}\tilde{F}^{A\mu\nu}\]

This term violates CP symmetry.

The strong CP problem:

The neutron electric dipole moment (nEDM) is proportional to $\theta$.

Experimental bound: $|d_n| < 3 \times 10^{-26} \text{ e·cm}$

This implies: $|\theta| < 10^{-10}$

Why is $\theta$ so small? Natural values are of order 1.

This is the strong CP problem, one of the most puzzling fine-tunings in particle physics.

The Peccei-Quinn solution:

Peccei and Quinn (1977) proposed a new $U(1)_{PQ}$ symmetry that dynamically relaxes $\theta \to 0$.

The pseudo-Goldstone boson of the spontaneously broken $U(1)_{PQ}$: the axion.

The axion is a dark matter candidate. Its mass: $m_a \sim \Lambda_{QCD}^2/f_{PQ}$

For $f_{PQ} \sim 10^{12}$ GeV: $m_a \sim 10^{-5}$ eV

Multiple experiments are searching for the axion: ADMX, HAYSTAC, ABRACADABRA.

No detection yet.

VI. Chiral Symmetry and Its Breaking

In the limit of massless quarks, QCD has a chiral symmetry:

\[SU(N_f)_L \times SU(N_f)_R\]

Left-handed and right-handed quarks transform independently.

For $N_f = 2$ (u, d quarks, approximately massless):

\[SU(2)_L \times SU(2)_R \cong SO(4)\]

Spontaneous chiral symmetry breaking:

The QCD vacuum breaks this symmetry to the vector subgroup:

\[SU(2)_L \times SU(2)_R \to SU(2)_V\]

The chiral condensate forms:

\[\langle\bar{q}q\rangle = \langle\bar{u}u + \bar{d}d\rangle \neq 0\]

Computed from lattice QCD: $\langle\bar{q}q\rangle \approx -(250 \text{ MeV})^3$

The Goldstone bosons:

3 generators broken → 3 Goldstone bosons: pions $\pi^+, \pi^-, \pi^0$

The pions are pseudo-Goldstone bosons (not exactly massless because quark masses are nonzero):

\[m_\pi^2 = -\frac{(m_u + m_d)\langle\bar{q}q\rangle}{f_\pi^2}\]

The Gell-Mann-Oakes-Renner relation.

This is one of the most important results in hadronic physics: it connects the pion mass to the quark masses and the chiral condensate.

Chiral Perturbation Theory (ChPT):

Below the chiral symmetry breaking scale ($\Lambda_\chi \sim 4\pi f_\pi \sim 1$ GeV), the effective theory is ChPT, an expansion in powers of $p/\Lambda_\chi$.

The degrees of freedom: pions and other pseudo-Goldstone bosons. The interactions: dictated by the chiral symmetry.

ChPT is the low-energy effective theory of QCD. It is perturbative in momenta, even though QCD itself is non-perturbative at these scales.

At leading order: $\mathcal{L}_{ChPT}^{(2)} = \frac{f_\pi^2}{4} \text{Tr}[\partial_\mu U^\dagger\partial^\mu U + m_\pi^2(U + U^\dagger)]$

where $U = e^{i\pi^a\tau^a/f_\pi} \in SU(2)$ is the pion field.

VII. The Proton Mass: Where Your Mass Comes From

You are made of protons and neutrons. Each proton has mass $m_p \approx 938$ MeV.

Where does this mass come from?

The quark masses: $m_u \approx 2.3$ MeV $m_d \approx 4.8$ MeV Sum: $\approx 7.1$ MeV

The quarks contribute $\sim 1\%$ of the proton mass.

The other $99\%$:

The trace anomaly:

The proton mass comes from the QCD trace anomaly, the non-conservation of the energy-momentum tensor due to quantum effects.

The QCD contribution to the proton mass:

\[M_p = \langle p|\hat{T}^\mu_\mu|p\rangle\]

The trace of the energy-momentum tensor:

\[T^\mu_\mu = \sum_f m_f\bar{\psi}_f\psi_f + \frac{\beta(g_s)}{2g_s}F^A_{\mu\nu}F^{A\mu\nu}\]

The second term, the gluon contribution, is proportional to the beta function.

This is the trace anomaly: the classical scale invariance of massless QCD is broken by quantum effects.

The symmetry breaking generates a scale: $\Lambda_{QCD} \approx 200$ MeV.

The proton mass: $M_p \sim \Lambda_{QCD}$

The proton mass is generated by the quantum breaking of scale invariance.

It is not the Higgs mechanism. It is not the quark masses.

It is the energy of the chromodynamic field, the gluons and their quantum fluctuations, confined within the proton by confinement.

$E = mc^2$ applied to the QCD vacuum:

\[M_p \approx \Lambda_{QCD}/c^2\]

Your mass is QCD. Your mass is the trace anomaly. Your mass is the energy of permanently confined color fields.

VIII. Lattice QCD: Computing the Uncomputable

Perturbation theory fails at low energies. How do we compute QCD predictions there?

Lattice QCD:

Discretize spacetime on a lattice: spacing $a$, volume $L^4$.

The path integral becomes a finite sum:

\[Z = \int \mathcal{D}A\mathcal{D}\psi\mathcal{D}\bar{\psi}\, e^{iS_{QCD}} \to \sum_{\text{configs}} e^{-S_E}\]

In Euclidean spacetime, the oscillating phase $e^{iS}$ becomes a damping factor $e^{-S_E}$ that can be sampled by Monte Carlo methods.

The fermion determinant:

\[\int \mathcal{D}\psi\mathcal{D}\bar{\psi}\, e^{-\bar{\psi}M\psi} = \det M\]

Computed exactly (expensive) or with pseudofermion approximations.

What lattice QCD can compute:

Hadron masses: computed to 1-2% precision
Quark masses: determined precisely
Decay constants: $f_\pi, f_K, f_D, f_B$
Form factors: electroweak matrix elements
The QCD phase diagram
The chiral condensate

The hadron spectrum:

The lattice QCD computation of the hadron spectrum (BMW collaboration, 2008):

Starting from the quark masses and $\alpha_s$, compute the masses of: proton, neutron, $\Delta$, $\Lambda$, $\Sigma$, $\Xi$, $\Omega$ baryons and several mesons.

Agreement with experiment at the 2-3% level, with the discrepancies attributable to finite lattice spacing and volume effects.

This was the first time the proton mass was computed from first principles.

The computational challenge:

Lattice QCD requires enormous computational resources.

The computational cost scales as: $\sim (L/a)^4 \times (m_\pi a)^{-3}$

For physical pion mass ($m_\pi = 140$ MeV) and lattice spacing $a \sim 0.05$ fm:

The ratio $1/(m_\pi a) \sim 28$.

Cost: $\sim 10^{10}$ floating point operations per configuration. Thousands of configurations needed.

World-leading lattice QCD calculations use the largest available supercomputers, millions of CPU-hours per calculation.

The calculations are getting better. Exascale computing is making previously impossible calculations possible.

IX. The QCD Phase Diagram

At extreme temperatures and densities, QCD undergoes phase transitions.

The QCD phase diagram:

Axes: temperature $T$ and baryon chemical potential $\mu_B$.

Hadronic phase (low $T$, low $\mu_B$): Quarks confined in hadrons. Chiral symmetry broken. Normal nuclear matter.

Quark-Gluon Plasma (QGP) (high $T$): Above $T_c \approx 155$ MeV: quarks and gluons deconfined. Chiral symmetry restored.

Color Superconductor (high $\mu_B$, low $T$): Dense quark matter with quark-quark Cooper pairing. Color symmetry spontaneously broken. Relevant to neutron star interiors.

The deconfinement transition:

Lattice QCD at zero chemical potential: the transition from hadronic to QGP phase is a crossover, not a sharp phase transition, at $T_c \approx 155$ MeV.

At finite $\mu_B$: the crossover may become a first-order phase transition. The critical endpoint, where the crossover meets the first-order line, is the subject of intense theoretical and experimental search.

The Quark-Gluon Plasma in experiment:

The QGP was created at:

RHIC (Relativistic Heavy Ion Collider, BNL)
LHC (ALICE experiment, CERN)

Heavy ion collisions create temperatures $\sim 300$ MeV $\sim 3.5 \times 10^{12}$ K, conditions similar to the universe $\sim 1$ microsecond after the Big Bang.

The QGP behaves as a perfect fluid, with the lowest viscosity-to-entropy ratio ever measured:

\[\eta/s \approx \frac{1}{4\pi} \frac{\hbar}{k_B}\]

This lower bound was derived from the AdS/CFT correspondence (Kovtun, Son, Starinets, 2005).

The QGP saturates it.

The most strongly coupled fluid in nature corresponds to the most strongly coupled gravitational system via holography.

X. OPE and the Structure of QCD Amplitudes

The Operator Product Expansion (OPE):

One of the most powerful tools in QCD: the OPE (Wilson, 1969).

When two operators approach each other:

\[\mathcal{O}_i(x)\mathcal{O}_j(0) = \sum_k C^k_{ij}(x)\mathcal{O}_k(0)\]

The coefficient functions $C^k_{ij}(x)$ are calculable in perturbation theory. The operators $\mathcal{O}_k$ encode the non-perturbative information.

Deep Inelastic Scattering (DIS):

In DIS, an electron scatters off a quark inside a proton via virtual photon exchange.

The cross section involves:

\[W^{\mu\nu}(q,p) = \frac{1}{2\pi} \int d^4x\, e^{iq\cdot x} \langle p|[J^\mu(x), J^\nu(0)]|p\rangle\]

The OPE factorizes this into:

Hard part: Wilson coefficients, calculable in perturbation theory
Soft part: parton distribution functions (PDFs), non-perturbative, extracted from data

DGLAP evolution:

The PDFs depend on the resolution scale $\mu$. Their evolution is governed by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations:

\[\mu^2\frac{\partial f_i(x,\mu^2)}{\partial\mu^2} = \sum_j \int_x^1 \frac{dz}{z} P_{ij}(z)\, f_j(x/z, \mu^2)\]

The splitting functions $P_{ij}(z)$ are computed in perturbation theory. Currently known to 4-loop (N³LO) accuracy.

The DGLAP equations describe how the quark and gluon content of the proton changes as you probe it at higher resolution.

At higher $\mu$: more gluons and sea quarks visible. At lower $\mu$: fewer, coarser structure.

XI. What QCD Does Not Explain

QCD is the correct theory of the strong force. It is also incomplete in ways that matter.

Confinement not proven:

Despite decades of work, confinement has not been derived analytically from the QCD Lagrangian.

Lattice QCD demonstrates it numerically. The linear potential is observed in simulation. But the mechanism, why the flux tube forms, why the gauge field topology produces confinement, is not analytically understood.

This is the Clay Millennium Problem.

The strong CP problem:

Why is $\theta < 10^{-10}$? The Peccei-Quinn solution predicts the axion. The axion has not been found.

The proton spin crisis:

The total spin of the proton is $\frac{1}{2}$. The quark spins contribute only $\sim 30\%$.

Where is the rest?

Gluon spin: $\sim 20-30\%$. Orbital angular momentum: $\sim 40-50\%$.

The proton spin puzzle, the decomposition of the proton spin among quarks, gluons, and orbital motion, is not fully resolved.

The decomposition itself is gauge-dependent: defining “quark spin” and “gluon spin” separately requires a gauge choice.

Hadron spectroscopy:

Lattice QCD can compute hadron masses. Many are well understood.

But: exotic hadrons, states beyond the simple $q\bar{q}$ and $qqq$ picture, are being discovered experimentally and are not fully understood theoretically.

Tetraquarks ($qq\bar{q}\bar{q}$): $X(3872)$, $Z_c(3900)$, $T_{cc}^+$

Pentaquarks ($qqqq\bar{q}$): $P_c^+(4380)$, $P_c^+(4450)$ (LHCb, 2015)

Hybrid mesons ($q\bar{q}g$): predicted, not definitively confirmed.

Glueballs ($gg$ or $ggg$): lightest predicted around $1.6-1.7$ GeV, candidates observed but not confirmed.

The hadron spectrum is richer and stranger than the simple quark model. QCD generates all of these states. Computing them from first principles remains a challenge.

XII. QCD as the Paradigm

QCD is not only the theory of the strong force.

It is the paradigm for what a fundamental theory can be.

Asymptotic freedom taught us: the same theory can be weakly coupled at short distances and strongly coupled at long distances. The effective description depends on scale. This insight transformed our understanding of quantum field theory.

Confinement showed something stranger: the physical spectrum of a theory can be completely different from the fields that appear in the Lagrangian. Quarks and gluons are in the Lagrangian. Protons and pions are in the spectrum. The relationship is non-perturbative and not yet fully understood.

The QCD phase diagram taught us: the same theory can have completely different phases, hadronic matter, quark-gluon plasma, color superconductor, depending on the conditions. Phase structure is rich even in theories we fully know.

Lattice QCD proved something methodological: non-perturbative quantum field theories can be computed numerically with controlled approximations. This is a methodology that extends beyond QCD to other strongly coupled theories.

The $\theta$ vacuum and instantons taught us: quantum field theories have topological structure, vacuum configurations that cannot be continuously deformed into each other, and this structure has physical consequences.

Everything we have learned from QCD, the running coupling, the phase structure, the topological vacuum, the effective field theory, the OPE, the lattice regularization, applies, in various forms, to every other quantum field theory.

QCD is not just one theory. It is a school.

Coda

At the center of every proton, in every atom of your body, three quarks are confined by gluons in a space smaller than $10^{-15}$ meters.

The gluons are massless. They travel at the speed of light. But they are permanently trapped: by their own mutual interactions, by the topology of the color field, by confinement.

The proton has existed for 13.8 billion years.

In all that time, no quark has ever escaped. Aha!!! :)

The force that holds them is, in the infrared, at the scales that matter, the least understood of all the fundamental forces.

We have the Lagrangian. We have asymptotic freedom. We have lattice QCD.

What we don’t have is an analytic proof that confinement follows from the Lagrangian. We don’t know why $\theta$ is small, we don’t fully understand the proton spin, and we don’t know the full phase diagram at high density.

The quarks in your protons have kept their secrets for 13.8 billion years.

We have been asking for fifty.

We are getting closer…

Essential reading:

Gross & Wilczek, “Ultraviolet Behavior of Non-Abelian Gauge Theories” (1973) Phys. Rev. Lett. 30, 1343

Politzer, “Reliable Perturbative Results for Strong Interactions?” (1973) Phys. Rev. Lett. 30, 1346 (the asymptotic freedom papers; Nobel Prize 2004)

Belavin, Polyakov, Schwarz & Tyupkin, “Pseudoparticle Solutions of the Yang-Mills Equations” (1975) Phys. Lett. B 59, 85 (instantons)

‘t Hooft, “Computation of the Quantum Effects Due to a Four-Dimensional Pseudoparticle” (1976) Phys. Rev. D 14, 3432

Peccei & Quinn, “CP Conservation in the Presence of Pseudoparticles” (1977) Phys. Rev. Lett. 38, 1440 (the axion solution to strong CP)

For reviews:

Shifman, “Quark-Hadron Duality” (2000) arXiv:hep-ph/0009131

Ioffe, “QCD at Low Energies” (2006) Prog. Part. Nucl. Phys. 56, 232

For textbooks:

Peskin & Schroeder, An Introduction to Quantum Field Theory (Chapters 16-18 for QCD)

Manohar & Wise, Heavy Quark Physics (2000) Cambridge University Press

Donoghue, Golowich & Holstein, Introduction to the Standard Model (1992)

For lattice QCD:

DeGrand & DeTar, Lattice Methods for Quantum Chromodynamics (2006) World Scientific

Have a question about this piece? Ask →