Where Does Mass Come From? — The Physics of Craig D. Roberts

Chapter 1

The Mass Puzzle

Pick up a textbook on the Standard Model of particle physics and you will find a beautiful, precisely tested theory of quarks, leptons, and the forces that bind them. You will also find the Higgs field, celebrated since its experimental confirmation in 2012 as the origin of mass. And that story is true — but it is drastically incomplete.

The Higgs mechanism gives mass to the fundamental particles: the W and Z bosons that mediate the weak force, the electron, the muon, the quarks. For the electron, that is essentially the whole story — the electron's mass is its Higgs-generated mass. But for the proton, which accounts for most of the mass of every atom in your body, the situation is radically different.

A startling accounting

A proton is built from two up quarks and one down quark, bound together by gluons. The Higgs-generated ("current") masses of these quarks are tiny:

Quark masses vs. proton mass

Two up quarks and one down quark contribute roughly 9 MeV — less than 1% of the proton's mass. The gluons are strictly massless. So where does the other 99% come from?

~2%

~98% from QCD dynamics

The proton's mass budget. The Higgs mechanism (gold) accounts for barely 2%. The rest (blue) emerges from the dynamics of the strong force — Quantum Chromodynamics (QCD).

This is the mass puzzle. The mass of nearly everything visible in the universe — every star, planet, and person — is overwhelmingly not generated by the Higgs field. It is emergent: it arises from the collective behaviour of quarks and gluons interacting via QCD. Understanding how this happens is one of the central challenges of modern physics.

It is also the central achievement of Craig D. Roberts.

Craig D. Roberts in brief

International Distinguished Professor at Nanjing University and formerly Senior Physicist at Argonne National Laboratory. Over 260 peer-reviewed papers with more than 22,000 citations. The world's leading practitioner of nonperturbative continuum methods in QCD, and the architect of the Emergent Hadron Mass framework that explains how the strong force generates visible mass.

Chapter 2

QCD: The Strangely Beautiful Strong Force

Quantum Chromodynamics is the theory of the strong force — the force that holds quarks together inside protons and neutrons, and holds protons and neutrons together inside nuclei. It is, in a sense, the theory responsible for the existence of matter.

Like Quantum Electrodynamics (QED) before it, QCD is a gauge theory: its dynamics are dictated by a symmetry principle. In QED, the gauge symmetry is U(1) and the force carrier is the photon. In QCD, the gauge symmetry is $SU(3)_\text{colour}$ , and the force carriers are eight gluons.

The QCD Lagrangian (schematic)

The first term describes quarks (

) interacting with gluon fields through the covariant derivative $D_\mu$ . The second term describes the gluons themselves via their field-strength tensor $G^a_{\mu\nu}$ . The index

runs from 1 to 8 — one for each gluon. The mass

is the Higgs-generated current quark mass: small for up and down quarks, larger for strange, charm, bottom, and top.

Despite its compact form, QCD is vastly more complex than QED, for one critical reason: gluons carry colour charge. Photons are electrically neutral — they mediate the electromagnetic force but do not feel it themselves. Gluons, by contrast, are themselves colour-charged. They interact with each other. This self-interaction changes everything.

Asymptotic freedom and confinement

QCD exhibits two extraordinary properties that emerge from gluon self-interaction:

Asymptotic freedom (high energy)

At very high energies — equivalently, at very short distances — the QCD coupling becomes weak. Quarks and gluons behave almost like free particles. This is why perturbative calculations (Feynman diagrams, order by order in powers of the coupling) work brilliantly for high-energy collisions at the LHC. Discovered by Gross, Wilczek, and Politzer (2004 Nobel Prize).

Confinement (low energy)

At low energies — the scale of everyday matter — the coupling becomes strong. Quarks and gluons are never observed in isolation; they are permanently confined inside hadrons (protons, neutrons, pions, etc.). No one has ever detected a free quark. Understanding why remains one of the millennium prize problems in mathematics.

The running of the QCD coupling constant $\alpha_s(Q^2)$ — strong at low momentum transfer , weak at high — is the single most important feature of the theory:

Running coupling (one-loop, leading order)

Here

is the number of active quark flavours and $\Lambda_{\text{QCD}} \approx 250$ MeV is the scale at which the coupling becomes strong. Because gluons carry colour (contributing the "11" in $\beta_0$ , which dominates over the quark contribution " $\frac{2}{3}n_f$ "), the coupling decreases with energy — the opposite of QED. Below the scale $\Lambda_{\text{QCD}}$ , perturbation theory breaks down entirely.

Here lies the fundamental problem: the most interesting physics — the formation of protons, the generation of mass, the binding of nuclei — all happens in the strong-coupling regime where Feynman diagrams cannot be trusted. To understand how mass emerges, you need a different set of tools.

Deep dive: Why can't we just use lattice QCD?

Lattice QCD discretises spacetime onto a four-dimensional grid and uses enormous computers to evaluate the QCD path integral numerically. It is powerful and has produced landmark results — for example, calculating the proton mass to within 2% of experiment.

But lattice QCD has limitations. It works in Euclidean (imaginary-time) space, making it difficult to access real-time dynamics, scattering processes, and parton distribution functions directly. It is computationally expensive, limiting the range of observables that can be studied. And it produces numbers, not analytic insight — it can confirm that a proton has a certain mass, but it does not easily reveal the mechanism by which that mass is generated.

This is where continuum methods — and Craig Roberts' approach — become essential. They complement lattice QCD by providing analytic understanding of the underlying dynamics.

Chapter 3

Dyson–Schwinger Equations: The Equations of Motion of QCD

In the late 1940s, Freeman Dyson and Julian Schwinger independently derived a set of integral equations that are, in a precise sense, the equations of motion of any quantum field theory. Applied to QED, these Dyson–Schwinger equations (DSEs) reproduce every result you can get from Feynman diagrams — and more, because they are exact, not perturbative.

Applied to QCD, the DSEs become a system of coupled integral equations that relate every Green function (propagator, vertex, scattering amplitude) to every other. They are the native language of QCD. The challenge is that this system is infinite — you need to truncate it intelligently, preserving the symmetries that make QCD what it is.

This is Craig Roberts' core expertise. From the early 1990s onward, he has developed, refined, and applied symmetry-preserving truncation schemes for the DSEs that allow first-principles calculations of hadron properties in the nonperturbative domain.

The gap equation: where mass is born

The most important DSE is the quark gap equation — the equation of motion for the quark propagator. In QED, the electron propagator receives small corrections from virtual photons, slightly shifting the electron's mass. In QCD, the analogous corrections are so large that they transform the theory.

The quark gap equation

is the full (dressed) quark propagator — the object we want to solve for. On the right: the first term is the "bare" quark with its small Higgs mass $m_{\text{bm}}$ . The integral term represents the quark's interaction with the gluon field: $D_{\mu\nu}$ is the gluon propagator, $\Gamma^a_\nu$ is the quark-gluon vertex, and

is the coupling constant. This equation says: the quark is constantly emitting and reabsorbing gluons, and the cumulative effect of this dressing generates a large dynamical mass, even when the Higgs mass is negligibly small.

The dressed quark propagator can be decomposed as:

Dressed quark propagator

The mass function is M(p^2) = B(p^2)/A(p^2)

. At high momenta (short distances), M(p^2)

approaches the small current quark mass — the Higgs contribution. But at low momenta (large distances, hadronic scales), M(p^2)

grows to roughly 350–400 MeV — hundreds of times larger than the Higgs mass of the light quarks. This momentum-dependent, dynamically generated mass is the origin of most of the proton's mass.

M(p) — The quark mass function from DSE calculations. At low momenta (inside a hadron), the dynamical mass is roughly 350 MeV — more than 70 times the Higgs-generated mass. This is *emergent hadron mass* in action.

Deep dive: What does "symmetry-preserving truncation" mean?

The DSEs form an infinite tower of coupled equations. To solve them in practice, you must truncate — approximate the unknown higher-order Green functions in terms of the lower-order ones you are solving for.

A symmetry-preserving truncation is one that respects the Ward–Green–Takahashi identities of QCD. These identities are consequences of gauge symmetry and chiral symmetry, and they ensure relationships like the Goldberger–Treiman relation and the axial-vector Ward identity are satisfied. If your truncation violates these identities, your pions won't be Goldstone bosons and your results will be physically inconsistent.

Roberts and collaborators developed the "rainbow-ladder" truncation as the leading-order scheme, and then systematically improved it through the Qin–Chang interaction kernel and beyond-rainbow-ladder constructions — always preserving the essential symmetries. This is one of the most technically demanding aspects of the programme.

Chapter 4

Dynamical Chiral Symmetry Breaking: Mass from Nothing

If you could somehow switch off the Higgs field — set all current quark masses to zero — the QCD Lagrangian would possess an exact chiral symmetry. Left-handed and right-handed quarks would be completely independent, and no hadron made of light quarks could have mass.

But that is not what happens. Even with massless quarks, the gap equation still has a nontrivial solution: the dynamical mass function M(p^2) is still large at low momenta. The strong interaction, through the constant exchange of gluons, spontaneously generates mass where there was none. This is dynamical chiral symmetry breaking (DCSB).

DCSB in plain language

Imagine a perfectly symmetrical ball balanced on the peak of a Mexican hat. The system has rotational symmetry. But the ball must fall — and whichever direction it falls, the symmetry is broken. In QCD, the "ball" is the QCD vacuum, and the "falling" is the generation of a quark condensate $\langle\bar{q}q\rangle$ . The quarks inside hadrons acquire a large dynamical mass, roughly 350 MeV, even though their Lagrangian mass is essentially zero. This is a purely quantum-mechanical, nonperturbative phenomenon — it cannot be seen in any finite number of Feynman diagrams.

DCSB is the single most important emergent phenomenon in QCD. It is responsible for:

Consequence of DCSB	Significance
Generation of ~98% of visible mass	Explains why protons are heavy despite containing nearly massless quarks
Existence of pions as Goldstone bosons	Pions are anomalously light — the "signature" of DCSB
Quark condensate $\langle\bar{q}q\rangle$	An order parameter signalling the broken symmetry
Momentum-dependent quark mass	Quarks are "dressed" — their effective mass depends on the resolution scale
Connection between confinement and mass generation	These two features of QCD are deeply intertwined

Roberts was among the first to demonstrate rigorously, through DSE calculations, that DCSB occurs as a direct consequence of the strong gluon dressing of quarks, and to trace its quantitative consequences for the spectrum and structure of hadrons. His work showed that DCSB and confinement are not separate mysteries — they are intimately linked aspects of the same nonperturbative dynamics.

Deep dive: DCSB and the quark condensate

The quark condensate $\langle\bar{q}q\rangle$ was traditionally understood as a property of the QCD vacuum — a sea of virtual quark-antiquark pairs filling all of spacetime. This interpretation leads to an enormous vacuum energy density, contributing to the notorious "cosmological constant problem" (a mismatch of ~46 orders of magnitude between QCD predictions and the observed vacuum energy of the universe).

Roberts, together with Stanley Brodsky, proposed a radical reinterpretation: the condensate is not a property of empty space, but a property of hadrons. Quarks and gluons only fluctuate in and out of existence inside hadrons. The vacuum itself is much simpler — essentially trivial. This "in-hadron condensate" picture reduces the cosmological constant mismatch by a factor of $10^{45}$ , and is supported by light-front analyses of QCD.

Chapter 5

The Pion Paradox: The Lightest and the Deepest

The pion ( $\pi$ ) is the lightest hadron, with a mass of about 140 MeV — much lighter than the proton (938 MeV) despite being made of the same quarks. This is already surprising. But the real puzzle is deeper.

On one hand, the pion is a Nambu–Goldstone boson: the quantum of the field associated with DCSB, in the same way that phonons are the quanta of crystal lattice vibrations. Goldstone's theorem tells us that when a continuous symmetry is spontaneously broken, massless bosons appear. In the limit of zero current quark mass, the pion would be exactly massless. Its small physical mass (~140 MeV) is a direct measure of the small but nonzero Higgs masses of the up and down quarks.

On the other hand, the pion is a relativistic bound state of a quark and an antiquark, each with a dynamical mass of ~350 MeV. A naive expectation would be $m_\pi \approx 2 \times 350 = 700$ MeV. Yet the pion is five times lighter than that.

The pion mass relation (Gell-Mann–Oakes–Renner)

This exact relation links the pion mass ( $m_\pi$ ), the pion decay constant ( $f_\pi \approx 92$ MeV), the current quark mass ( m_q

), and the quark condensate ( $\langle\bar{q}q\rangle$ ). As $m_q \to 0$ , the pion mass vanishes: $m_\pi \to 0$ . The pion is simultaneously a tightly bound quark-antiquark system with strong internal dynamics and an almost-massless collective mode of the vacuum.

How can both things be true at once? Roberts' DSE framework resolves the paradox quantitatively. The Bethe–Salpeter equation — the relativistic bound-state equation derived from the DSEs — shows that the pion's wave function has a specific structure enforced by chiral symmetry. The enormous binding energy (two ~350 MeV quarks producing a ~140 MeV meson) is not mysterious: it is a precise, calculable consequence of the interplay between DCSB and the Goldstone nature of the pion.

In Roberts' famous phrase, the pion is "an enigma within the Standard Model" — the simplest hadron, yet the most profound, because its properties encode the mechanism of mass generation more directly than any other particle.

Deep dive: The Bethe–Salpeter equation

To describe bound states in quantum field theory, one solves the Bethe–Salpeter equation (BSE):

Bethe–Salpeter equation for a meson

$\Gamma(p;P)$ is the meson's Bethe–Salpeter amplitude (its quantum-field-theoretic wave function),

are dressed quark propagators (from the gap equation),

is the quark-antiquark scattering kernel, and

is the total meson momentum. The meson mass is the value of P^2

at which this equation has a solution. For the pion, the solution satisfies the axial-vector Ward–Takahashi identity — a non-trivial consistency check linking the gap equation and the BSE, and the mathematical expression of the pion's Goldstone nature.

Roberts and Peter Maris showed in their landmark 1997–98 papers that a single interaction, applied consistently to both the gap equation and the BSE, simultaneously predicts the correct pion mass, decay constant, and electromagnetic form factor. This was a pivotal demonstration that DSEs could be a precision tool for hadron physics.

Chapter 6

Emergent Hadron Mass: The Three Pillars

By the 2010s, Roberts' programme had matured into a coherent theoretical framework called Emergent Hadron Mass (EHM). This framework identifies three interconnected manifestations of mass generation in QCD — three "pillars" that can each be independently tested by experiment.

◆

Pillar 1: Running Gluon Mass

Gluons are massless in the QCD Lagrangian, but nonperturbative dynamics generate an effective gluon mass scale of ~500 MeV at low momenta — via the Schwinger mechanism. This tames the infrared behaviour of QCD and seeds mass generation.

◆

Pillar 2: Process-Independent Effective Charge

The QCD running coupling freezes to a finite value at low momenta — it does not diverge. This "infrared completion" of the coupling defines a process-independent effective charge that governs all nonperturbative QCD phenomena. Three independent approaches — Deur's Jefferson Lab measurements, Brodsky's light-front holography, and Roberts' DSE calculations — all converge on the same answer, without adjustable parameters.

◆

Pillar 3: Running Quark Mass

The quark mass function M(p^2)

— large at low momenta, falling to the current mass at high momenta — is the direct signature of DCSB. It is the "smoking gun" of emergent mass in the quark sector.

These three pillars are not independent: they are connected by the DSEs. The gluon mass scale (Pillar 1) feeds into the effective charge (Pillar 2), which drives the gap equation and generates the running quark mass (Pillar 3). Together, they explain why hadrons are heavy, why pions are light, and why quarks are confined.

The Schwinger mechanism for gluon mass

How can a massless particle — the gluon — acquire mass without breaking gauge symmetry? In 1962, Julian Schwinger showed that in certain gauge theories, the gauge boson can become massive through a purely dynamical mechanism, without a Higgs-like field. Roberts and collaborators (especially Binosi, Papavassiliou, and Rodríguez-Quintero) demonstrated that this is precisely what happens in QCD: a massless, colour-carrying gluon+gluon correlation emerges in the dressed three-gluon vertex, giving the gluon an effective mass of about 500 MeV.

Gluon propagator in the infrared

Instead of diverging as 1/q^2

like a massless particle's propagator, the gluon propagator saturates to a finite value at zero momentum. This is equivalent to the gluon having acquired a dynamical mass m_g

. Both DSE calculations and lattice QCD simulations agree on this result — a remarkable convergence of independent methods.

Illuminating "Terra Damnata": The Running Coupling Breakthrough

For half a century after the discovery of asymptotic freedom (2004 Nobel Prize), the infrared region of QCD — where the coupling $\alpha_s$ grows large and perturbation theory fails — was terra incognita, a domain so forbidding that physicists dubbed it "Terra Damnata" (cursed land). The standard short-distance calculation predicted that $\alpha_s$ would blow up to infinity at long range — a mathematical singularity called the Landau pole — but this was widely understood to be a failure of the method rather than a genuine physical prediction. The real question was: what does $\alpha_s$ actually do in this region?

A breakthrough emerged from the convergence of three entirely independent approaches:

1. Experiment (Alexandre Deur, Jefferson Lab): Using data from polarised deep-inelastic scattering at JLab's CEBAF facility, Deur extracted $\alpha_s$ via the Bjorken integral — a quantity that, by a stroke of luck, filters out the complex multi-quark effects that would contaminate most other measurements. His data showed $\alpha_s$ stops growing and freezes to a constant value in the infrared.

2. Light-Front Holography (Stanley Brodsky, SLAC/Stanford): Using a mathematical device known as the gauge/gravity duality, Brodsky and collaborators (de Téramond Peralta, Dosch) mapped the strong force from a five-dimensional gravitational theory to four-dimensional QCD. Their holographic prediction for $\alpha_s$ matched Deur's data with no adjustable parameters.

3. QCD Equations of Motion (Craig Roberts, DSEs): Roberts' lifelong programme of solving the Dyson–Schwinger equations provided the third independent route. Working with Binosi, Papavassiliou, Rodríguez-Quintero, Chang, and Mezrag, Roberts showed that the emergent gluon mass tames the infrared coupling. The DSE result was parameter-free and virtually indistinguishable from both Deur's measurements and Brodsky's holographic calculation.

The agreement among three completely different methods — experiment, holography, and continuum QCD — without any parameter tuning is one of the most striking convergences in modern particle physics. It demonstrates that there exists a unique, universal QCD effective charge that is well-defined at all momentum scales, from the perturbative ultraviolet down through Terra Damnata to zero momentum.

This discovery has profound implications. Because $\alpha_s$ is finite everywhere, QCD becomes the first full quantum field theory that predicts only finite quantities — unlike quantum electrodynamics, which develops its own Landau pole at ultra-high energies. Moreover, knowing the coupling at all scales opens the door to parameter-free predictions of hadron properties from first principles — predictions that underpin every other chapter of the EHM programme.

This convergence was described for a broad audience in S. J. Brodsky, A. Deur, and C. D. Roberts, "The Secret to the Strongest Force in the Universe," Scientific American, Vol. 330 No. 5 (May 2024), p. 32 (doi:10.1038/scientificamerican0524-32).

Deep dive: Higgs mass vs. emergent mass — a comparison

Feature	Higgs mechanism	Emergent Hadron Mass
Origin	Coupling to Higgs field	Nonperturbative QCD dynamics
Operates on	All fundamental particles	Quarks and gluons inside hadrons
Scale	~2% of proton mass	~98% of proton mass
Symmetry broken	Electroweak $SU(2) \times U(1)$	Chiral $SU(N_f)_L \times SU(N_f)_R$
Mechanism	Scalar field with nonzero VEV	Strong gluon self-interaction (Schwinger mechanism + DCSB)
Goldstone bosons	Eaten by W, Z	Pions (pseudo-Goldstone)
Experimental probe	LHC	JLab, EIC, AMBER, J-PARC

Chapter 7

Inside Hadrons: What Roberts' Framework Predicts

A theory is only as good as its predictions. Roberts' DSE/EHM programme has generated a vast number of predictions for measurable properties of hadrons — many of which have been confirmed by experiment, and many more that are now being tested.

Pion and kaon structure

The pion's internal structure has been a central focus. Roberts' group has predicted the pion's parton distribution function (PDF) — the probability of finding a quark carrying a fraction of the pion's momentum when probed at high energy. A key prediction is the large- behaviour:

Pion valence quark distribution at large x

This

power law is a direct consequence of DCSB. It differs from naive perturbative QCD predictions and from some model calculations. It has become a benchmark prediction that the Electron–Ion Collider (EIC), expected to begin operations in the 2030s, will test definitively.

The group has also predicted the pion's distribution amplitude — a related quantity describing how momentum is shared between quark and antiquark at a given resolution scale. Roberts' DSE calculations yield a distribution amplitude that is significantly broader than the asymptotic form, a prediction now confirmed by lattice QCD calculations.

Baryon structure: the quark-diquark picture

For baryons (three-quark systems like the proton and neutron), the full three-body Faddeev equation is computationally demanding. Roberts and collaborators developed a powerful simplification: the quark-diquark picture.

Inside a baryon, two of the three quarks form a tightly correlated pair — a diquark. This is not a fundamental particle but an emergent correlation, arising naturally from the same QCD dynamics that generate DCSB. The baryon is then described as a quark orbiting a diquark, bound by the exchange of dressed gluons.

The quark-diquark picture of the proton. Two quarks form a correlated diquark (green dashed ellipse), and the third quark (blue) interacts with it through gluon exchange (gold). This is not a model assumption — it emerges from solving the Faddeev equation.

Using this framework, Roberts' group has predicted electromagnetic form factors of the proton, neutron, and many other baryons; the spectrum of excited baryon states (including the enigmatic Roper resonance); nucleon-to-Delta transition form factors; and baryon parton distribution functions — all from the same underlying QCD interaction.

Form factors: mapping the shape of hadrons

An electromagnetic form factor describes how a hadron responds to being "photographed" by a virtual photon. It encodes the spatial distribution of charge and magnetism inside the hadron, as a function of the photon's resolution (momentum transfer Q^2 ).

Roberts' group has produced predictions for pion, kaon, nucleon, and Delta form factors that agree with Jefferson Lab data where available, and extend to the high- Q^2 domain that will be probed by the JLab 12 GeV upgrade. A key prediction: at sufficiently high Q^2 , the pion form factor $Q^2 F_\pi(Q^2)$ should approach a value set by the scale of DCSB, not by the perturbative asymptotic limit — a direct observable signal of emergent mass.

Chapter 8

Testing the Theory: From Jefferson Lab to the Electron–Ion Collider

A distinguishing feature of Roberts' programme is its tight coupling to experiment. Unlike some theoretical frameworks that make predictions only for asymptotic limits or unobservable quantities, the DSE/EHM approach produces predictions for real observables at real experimental facilities.

Jefferson Lab (JLab), Virginia, USA

The 12 GeV Continuous Electron Beam Accelerator Facility (CEBAF) at JLab is the world's premier facility for studying hadron structure. Its high-luminosity electron beam can probe the internal structure of protons, neutrons, and mesons with unprecedented precision. Many of Roberts' predictions — pion form factors, nucleon electromagnetic form factors, nucleon resonance transitions — are being directly tested by JLab experiments.

Electron–Ion Collider (EIC)

Under construction at Brookhaven National Laboratory, the EIC will collide polarised electrons with protons and ions, enabling the first detailed maps of the gluon content of hadrons, the three-dimensional structure of the proton, and the partonic structure of pions and kaons (accessed via the Sullivan process). Roberts has been a leading voice in defining the EIC science case, particularly for meson structure measurements that will directly test the EHM framework.

Five key measurements that will test EHM:

Measurement	What it tests	Facility
Pion elastic form factor at $Q^2 > 6$ GeV	Scale of DCSB in the pion	JLab 12 GeV
Pion valence quark PDF	Large- behaviour: prediction	EIC, AMBER
Kaon-to-pion form factor ratio	Higgs vs. emergent mass interplay	JLab, EIC
Pion and proton gravitational form factors	Mass and pressure distribution inside hadrons	JLab, EIC
Nucleon-to-Roper transition	Nature of baryon excited states and diquarks	JLab CLAS12

Chapter 9

A Career in Four Acts

Roberts' research trajectory mirrors the development of the field itself. His career can be understood as a sustained, systematic programme to answer a single question: How does QCD build matter?

1985–1995: Foundations

PhD at Flinders University (Adelaide, 1988) followed by a move to Argonne National Laboratory. Early work establishing the DSE framework as a viable nonperturbative tool for QCD. Key contributions: demonstrating DCSB in the gap equation, clarifying the relationship between confinement and the analytic structure of the quark propagator, and connecting DSE results to lattice QCD.

1996–2008: Hadron physics revolution

Collaboration with Peter Maris and Peter Tandy produces the foundational DSE calculations of pion properties (mass, decay constant, form factor). Development of the quark-diquark picture for baryons. Landmark review articles (2000, 2003) that established DSEs as the standard continuum approach to nonperturbative QCD. Joint work with Brodsky on in-hadron condensates and the cosmological constant problem.

2009–2018: Precision and prediction

Shift toward precision predictions for JLab experiments. The Qin–Chang interaction kernel (2011) provides a more sophisticated description of the quark-gluon interaction. Extensive calculations of nucleon form factors, transition form factors, baryon spectra, and parton distributions. Predictions for the pion's distribution amplitude and valence quark distribution become benchmark tests for lattice QCD and future experiments.

2019–present: Emergent Hadron Mass

Moves to Nanjing University as International Distinguished Professor and Head of the Institute for Nonperturbative Physics. The EHM framework crystallises: three pillars, each with clear experimental signatures. Active role in defining the science case for the EIC. Major reviews with Deur and Brodsky on the QCD running coupling, and with Ding and Schmidt on hadron mass and structure. In May 2024, co-authors (with Brodsky and Deur) a landmark Scientific American feature article explaining the convergence of three independent approaches to the QCD coupling — bringing the EHM story to the broadest possible audience. Over 112 papers now available as open-access preprints.

Chapter 10

The Bigger Picture: Why This Matters

The question "where does mass come from?" is not merely academic. It connects to some of the deepest open problems in physics.

The Standard Model is not complete

The Higgs mechanism explains how electroweak symmetry is broken and why the W, Z, and fundamental fermions have mass. But it says nothing about why the quark masses have the values they do — these are free parameters. And it accounts for only a sliver of the mass of composite matter. Understanding EHM completes the picture: the Standard Model generates mass through two mechanisms — the Higgs field for fundamental particles, and QCD dynamics for composite hadrons — and the interplay between them shapes the visible universe.

Nuclear physics from first principles

Ultimately, nuclear physics should be derivable from QCD. Roberts' framework provides the bridge: by computing hadron properties from QCD, and then using those properties as inputs to nuclear force models, one can work toward a first-principles understanding of nuclear structure, nuclear reactions, and the equation of state of dense matter in neutron stars.

Confinement remains unsolved

Despite decades of work, no one has rigorously proved that QCD confines quarks — it remains a Clay Millennium Prize problem. Roberts' approach provides deep physical insight: his calculations show that confinement is reflected in the analytic structure of the quark propagator (which has no real mass pole, meaning a free quark cannot propagate as a physical particle). The dressed quark is "confined" in the precise mathematical sense that it cannot appear in the asymptotic states of the theory.

A legacy of method and insight

Perhaps Roberts' most lasting contribution is methodological: he demonstrated that continuum quantum field theory, applied with care and symmetry-preserving rigour, can be a precision tool for the strong interaction. Before his programme, many physicists viewed the DSEs as unwieldy or unreliable. Today, the "continuum Schwinger function methods" approach is recognised alongside lattice QCD as one of the two pillars of nonperturbative QCD — and the two methods increasingly agree in their overlapping domain of applicability.

By the numbers

Over 369 peer-reviewed publications archived here, spanning 1993 to 2026. More than 22,000 citations. An h-index above 75. Collaborations spanning five continents. A generation of students and postdocs trained in nonperturbative QCD methods. And a theoretical framework — Emergent Hadron Mass — that is now shaping the scientific programme of the world's next great particle physics facility, the Electron–Ion Collider.

Chapter 11

369 Papers: The Full Scope

The archive assembled here contains 369 of Roberts' papers, spanning from 1993 to 2026. The sheer breadth of topics reveals something important: this is not a collection of isolated investigations but a single sustained programme, where each thread connects back to the central question of how QCD builds matter.

The corpus at a glance

Research theme	Papers	Years active
Baryon physics (spectrum, diquarks, form factors, Roper resonance)	~93	1996–2026
Pion and kaon properties (mass, form factors, PDFs, DAs, GPDs)	~72	1993–2026
Quark propagator, gap equation, and DCSB	~40	1993–2024
Heavy quark systems (charm, bottom, charmonia)	~15	1997–2026
Form factors (electromagnetic, gravitational, transition)	~11	1996–2025
Finite temperature and density / quark-gluon plasma	~9	1997–2016
Emergent Hadron Mass reviews	~7	2014–2024
Confinement and Schwinger functions	~6	1996–2024
Quark-gluon vertex and Ward–Takahashi identities	~5	1993–2016
Parton distributions (PDFs, GPDs, TMDs)	~4	2001–2024
QCD running coupling and effective charge	~3	2016–2023
Reviews, perspectives, and white papers	~9	2003–2026

What the newest papers reveal

The 2025–2026 papers show the programme at its most mature. Several recent highlights:

Gravitational form factors (2025)

For the first time, unified predictions for pion, kaon, and nucleon gravitational form factors from a single framework. A remarkable finding: the mass radius of each hadron is smaller than its charge radius, and the core pressures inside these particles exceed those at the centre of a neutron star. These quantities will be measurable at JLab and the EIC.

Orbital angular momentum in the pion and kaon (2025)

Even the pion — the simplest hadron — contains significant intrinsic orbital angular momentum. The pion's quark-antiquark pair is roughly a 50/50 mixture of zero and one units of orbital angular momentum. The kaon is 60/40, reflecting the heavier strange quark. This overturns the naive picture of the pion as a simple -wave bound state.

Fragmentation functions (2025)

Roberts' group showed that fragmentation functions — describing how a quark or gluon produced in a high-energy collision turns back into hadrons — can be derived from parton distribution functions using a crossing relation. This connects two apparently different experiments: deep inelastic scattering (probing hadron structure) and e^+e^- annihilation (producing hadrons from energy).

The Roper resonance — finally understood (2017–2025)

For fifty years, the Roper resonance N(1440) was a puzzle: the proton's first excited state appeared lighter than expected. Roberts' framework revealed it as the proton's first radial excitation — a dressed-quark core surrounded by a meson cloud that reduces its mass by about 20%. Multiple papers through 2025 now connect this picture to CLAS12 data, providing a quantitative resolution of one of hadron physics' oldest mysteries.

The most recent paper in the collection, from 2026, extends the continuum Schwinger methods to describe meson and baryon structure in a single unified treatment — a fitting capstone to three decades of work building a coherent picture of hadronic matter.

Chapter 12

A Physics Glossary for the Curious Reader

The language of quantum field theory can be forbidding. This glossary explains the key concepts in Roberts' work, grouped by theme rather than alphabetically, so each idea builds on the last.

The cast of characters

Quarks

Quarks are fundamental particles — they have no known internal structure. There are six types (called flavours): up (), down (), strange (), charm (), bottom (), and top (). Ordinary matter is built from up and down quarks. Each quark carries a fractional electric charge ( $+\frac{2}{3}$ for up, $-\frac{1}{3}$ for down) and a property called colour charge — the QCD analogue of electric charge. Colour comes in three types: red, green, and blue (these are just labels; they have nothing to do with visible light).

Gluons

Gluons are the force carriers of QCD — they mediate the strong force between quarks, just as photons mediate the electromagnetic force between charged particles. There are eight gluons, each carrying a combination of colour and anti-colour charge. Crucially, because gluons themselves carry colour, they interact with each other. This self-interaction is what makes QCD so much richer (and harder) than electromagnetism. Gluons are massless in the QCD Lagrangian, but acquire a dynamical mass of about 500 MeV through nonperturbative effects — one of the three pillars of Emergent Hadron Mass.

Hadrons, mesons, and baryons

Hadrons are composite particles made of quarks bound by gluons. They come in two main families. Mesons are made of one quark and one antiquark (e.g., the pion $\pi$ , the kaon ). Baryons are made of three quarks (e.g., the proton , the neutron , the Delta $\Delta$ ). All hadrons are colour-neutral: the colour charges of their constituents cancel out, which is why isolated colour is never observed — this is confinement.

The pion ( $\pi$ )

The pion deserves its own entry because it plays a unique role in QCD. It is the lightest hadron (~140 MeV) and the Nambu–Goldstone boson of dynamical chiral symmetry breaking (see below). In nuclear physics, pions are the primary carriers of the force between protons and neutrons. Roberts calls the pion "an enigma within the Standard Model" because its small mass conceals enormously complex internal dynamics: it is simultaneously an almost-massless collective mode and a tightly bound relativistic system whose constituent quarks each have a dynamical mass of ~350 MeV.

Diquarks

Diquarks are not fundamental particles but emergent correlations — two quarks inside a baryon that form a tightly bound pair. They arise naturally from the same QCD dynamics that cause DCSB. The baryon can then be understood as a quark orbiting a diquark, simplifying the three-body problem into a more tractable two-body one. Roberts' work showed that diquarks are essential for understanding baryon form factors, the baryon spectrum, and even exotic states like tetraquarks and pentaquarks. Importantly, diquarks are confined — they cannot exist as free particles, only inside hadrons.

Forces, symmetries, and how they break

Quantum Chromodynamics (QCD)

QCD is the quantum field theory of the strong force. "Chromo" refers to colour charge. It is part of the Standard Model and is described by the gauge symmetry group SU(3) . QCD is simple to write down (the Lagrangian fits on one line) but extraordinarily difficult to solve at low energies, where the coupling constant is large and perturbation theory fails. All of Roberts' work is essentially about finding ways to extract predictions from QCD in this strong-coupling regime.

Asymptotic freedom

At very high energies (or equivalently, very short distances), the strong force becomes weak. Quarks inside a proton that are probed with a very energetic photon behave almost as free particles. This property, called asymptotic freedom, was discovered by Gross, Wilczek, and Politzer (2004 Nobel Prize). It means perturbative methods (Feynman diagrams) work well at high energies — but fail completely at the energy scale of ordinary hadrons.

Confinement

No one has ever observed a free quark or gluon. They are permanently locked inside hadrons — this is confinement. Pull two quarks apart and the energy stored in the gluon field between them grows until it becomes energetically favourable to create a new quark-antiquark pair from the vacuum, producing new hadrons rather than free quarks. Proving confinement mathematically from QCD is one of the Clay Millennium Prize problems (with a $1 million reward). Roberts' work provides deep physical insight: the dressed quark propagator has no real mass pole, meaning a quark cannot propagate as a physical on-shell particle.

Chiral symmetry and its breaking

Chirality (from the Greek for "hand") refers to a particle's handedness. A massless spin- $\frac{1}{2}$ particle can be either left-handed or right-handed, and these two types are completely independent — this is chiral symmetry. If quarks were truly massless, QCD would possess exact chiral symmetry.

In reality, chiral symmetry is broken in two ways. First, the Higgs mechanism gives quarks small masses, breaking the symmetry explicitly (but only slightly, for up and down quarks). Second, and far more importantly, the strong interaction itself breaks chiral symmetry spontaneously — this is dynamical chiral symmetry breaking (DCSB). Even if Higgs masses were zero, the QCD vacuum would still generate a large effective quark mass (~350 MeV) through the constant exchange of gluons. DCSB is the mechanism responsible for ~98% of visible mass.

Nambu–Goldstone bosons

Goldstone's theorem states: whenever a continuous symmetry of a system is spontaneously broken, massless bosons appear — one for each broken symmetry direction. When chiral symmetry is spontaneously broken in QCD, the resulting (pseudo-)Goldstone bosons are the pions (and kaons, and the eta meson). They would be exactly massless if not for the small explicit breaking from Higgs-generated quark masses. The pion's small mass (~140 MeV) is thus a direct measure of how much explicit symmetry breaking the Higgs field contributes — it is a window into the interplay between the two mass-generation mechanisms.

The Schwinger mechanism

In 1962, Julian Schwinger showed that gauge bosons can acquire mass through purely dynamical effects, without a Higgs-like field, provided certain conditions on the theory's dynamics are met. In QCD, this mechanism operates on gluons: a massless gluon+gluon correlation emerges in the three-gluon vertex, effectively giving gluons a mass of ~500 MeV at low momenta. This "gluon mass scale" is the first pillar of Emergent Hadron Mass.

The theoretical tools

Dyson–Schwinger equations (DSEs)

The DSEs are the exact equations of motion of a quantum field theory, derived from the generating functional. For every propagator and vertex in the theory, there is a DSE that relates it to all the others, forming an infinite tower of coupled integral equations. In QCD, they play the role that Newton's F = ma plays in classical mechanics: they are the fundamental dynamical equations.

Because the system is infinite, it must be truncated — approximated in a controlled way. Roberts' core methodological contribution is the development of symmetry-preserving truncations that respect the Ward–Green–Takahashi identities of QCD, ensuring that the approximation does not destroy the physics (e.g., pions remain Goldstone bosons, gauge invariance is maintained).

The gap equation

The gap equation is the DSE for the quark propagator. It determines the mass function M(p^2) — how the effective mass of a quark depends on the momentum at which it is probed. Solving the gap equation is where DCSB happens mathematically: even with zero input mass, the equation has a solution where $M(0) \approx 350$ MeV. This "gap" between the trivial (massless) solution and the physical (massive) one is analogous to the energy gap in a superconductor — hence the name.

The Bethe–Salpeter equation (BSE)

The BSE is the relativistic equation for two-body bound states. In the DSE framework, it determines the properties of mesons: their masses, decay constants, and internal wave functions. The quark propagators that enter the BSE come from the gap equation, and the interaction kernel must be consistent with both — this is the essence of a symmetry-preserving approach. Solving the BSE for the pion, for example, simultaneously predicts its mass, its decay constant $f_\pi$ , and its Goldstone nature.

The Faddeev equation

The Faddeev equation extends the bound-state formalism to three-body systems — i.e., baryons. It is computationally more demanding than the BSE. In practice, Roberts' group often uses the quark-diquark simplification: two quarks form a diquark correlation, reducing the three-body problem to an effective two-body one (quark + diquark). This has been shown to reproduce results from the full three-body calculation to within a few percent for most observables.

Continuum Schwinger function methods (CSMs)

This is the modern umbrella term for the entire DSE-based approach, coined to emphasise that the methods work with the Schwinger functions (Green functions in Euclidean space) of QCD. The name also acknowledges the role of Julian Schwinger's foundational contributions to quantum field theory. When Roberts' recent papers refer to "CSMs," they mean the mature framework encompassing DSEs, BSEs, Faddeev equations, and symmetry-preserving truncations.

What experiments measure

Form factors

A form factor encodes the internal structure of a particle as "seen" by a particular probe. An electromagnetic form factor describes the distribution of electric charge and magnetism inside a hadron, measured by scattering electrons off it at various momentum transfers Q^2 . A gravitational form factor describes the distribution of mass, pressure, and angular momentum — probed indirectly through deeply virtual Compton scattering or meson production. Roberts' recent work shows that a hadron's mass radius (from gravitational form factors) is systematically smaller than its charge radius (from electromagnetic form factors) — a prediction now being tested at JLab.

Parton distribution functions (PDFs)

When a hadron is probed at high energy, its interior is revealed as a collection of partons — quarks, antiquarks, and gluons, each carrying a fraction of the hadron's total momentum. A PDF, q(x) , gives the probability of finding a parton of type carrying momentum fraction . Roberts' key prediction for the pion is that its valence quark PDF falls as (1-x)^2 at large — a direct signal of DCSB that differs from naive perturbative expectations and will be measured at the EIC.

Distribution amplitudes (DAs)

A distribution amplitude describes how momentum is shared between the constituents of a hadron at a given resolution scale. It is related to the hadron's light-front wave function. Roberts showed that the pion's DA is significantly broader than the "asymptotic" form predicted by perturbative QCD alone — it is dilated by DCSB. This was later confirmed by lattice QCD calculations.

Generalised parton distributions (GPDs) and TMDs

GPDs generalise both form factors and PDFs into a single framework, providing a three-dimensional picture of hadron structure in combined position and momentum space. Transverse momentum dependent distributions (TMDs) add information about parton motion perpendicular to the hadron's direction of travel. Both are among the primary observables planned for the Electron–Ion Collider. Roberts' group has provided predictions for pion and kaon GPDs and TMDs, including the recently calculated pion Boer–Mulders function (a TMD sensitive to the correlation between quark spin and orbital motion).

Fragmentation functions

When a quark or gluon is produced in a high-energy collision, it cannot exist freely — it must hadronise, forming new hadrons from the vacuum. A fragmentation function describes the probability that a parton of type produces a hadron of type carrying a fraction of the parton's momentum. Roberts' 2025 work demonstrated that fragmentation functions can be derived from parton distributions using a mathematical crossing relation — a striking unification of two apparently unrelated measurements.

Emergent Hadron Mass: the key concepts

Running masses and couplings

In quantum field theory, the "constants" of nature are not really constant — they run with the energy scale at which they are measured. The quark mass is large (~350 MeV) at low energies but shrinks to a few MeV at high energies. The strong coupling $\alpha_s$ is large at low energies but shrinks logarithmically at high energies. The gluon behaves as if it has a mass of ~500 MeV at low momentum but is effectively massless at high momentum. All three of these running quantities are pillars of Emergent Hadron Mass.

Effective charge

The process-independent effective charge is a particular definition of the QCD coupling that is well-defined at all momentum scales, including in the infrared (low-energy) regime — the "Terra Damnata" — where the conventional coupling diverges. Roberts, together with Brodsky and Deur, showed that this effective charge saturates to a finite value at zero momentum — it does not blow up. Three entirely independent approaches (Deur's Jefferson Lab measurements, Brodsky's light-front holography, and Roberts' DSE calculations) converge on the same result with no adjustable parameters. This infrared-finite coupling unifies perturbative and nonperturbative QCD, and enables parameter-free predictions of hadron properties.

The quark condensate and the vacuum

The quark condensate $\langle\bar{q}q\rangle$ is a measure of DCSB. Traditionally, it was interpreted as a property of empty space — a dense sea of virtual quark-antiquark pairs filling the vacuum. Roberts and Brodsky overturned this picture, arguing that condensates are properties of hadrons, not the vacuum. This "in-hadron condensate" reinterpretation has profound implications: it reduces the QCD contribution to the cosmological constant by a factor of $10^{45}$ , essentially resolving one of the largest numerical discrepancies in theoretical physics.

The Roper resonance

The Roper resonance N(1440) is the first excited state of the proton. For fifty years, it was a puzzle: constituent quark models predicted it should be heavier than the negative-parity N(1535) , but it is lighter. Roberts' framework explains this naturally: the Roper is the proton's first radial excitation (like the state of hydrogen compared to the ), with a dressed-quark core whose mass is reduced by ~20% through coupling to the pion cloud. This quantitative understanding was confirmed by Jefferson Lab CLAS data in 2017–2025.

Orbital angular momentum in hadrons

Angular momentum in quantum mechanics comes in two forms: spin (intrinsic rotation) and orbital (motion around a centre). Even the pion, long thought of as a simple -wave ( l=0 ) state, contains substantial orbital angular momentum. Roberts' 2025 calculations show the pion is roughly a 50/50 superposition of l=0 and l=1 components. This is a purely relativistic effect — it reflects the fact that quarks inside hadrons move at close to the speed of light, making nonrelativistic approximations qualitatively misleading.

Key Papers and Resources

All of Roberts' papers are freely available on the arXiv preprint server. Here are selected entry points, grouped by theme:

Reviews and overviews

nucl-th/0005064 — DSEs: Density, Temperature and Continuum Strong QCD (2000, ~1000 citations)
nucl-th/0301049 — DSEs: A Tool for Hadron Physics (2003)
1203.5341 — Strong QCD and DSEs (2012)
2102.01765 — Insights into Emergence of Mass from Pion and Kaon Structure (2021)
2303.00723 — QCD Running Couplings and Effective Charges (2023)
2211.07763 — Emergence of Hadron Mass and Structure (2023)

The pion

nucl-th/9707003 — Pion Mass and Decay Constant (1998)
1602.04016 — The Pion: An Enigma within the Standard Model (2016)
1301.0324 — Imaging DCSB: Pion Wave Function on the Light Front (2013, PRL)
2403.00629 — Pseudoscalar Mesons and Emergent Mass (2024)

Emergent Hadron Mass

2006.08782 — Reflections upon the Emergence of Hadronic Mass (2020)
2009.04011 — Empirical Consequences of Emergent Mass (2020)
2101.08340 — On Mass and Matter (2021)

Baryon physics

1902.00026 — Spectrum of Light- and Heavy-Baryons (2019)
1603.02722 — Baryons and the Borromeo (2016)
2008.07630 — Diquark Correlations in Hadron Physics (2020)

Popular science

Scientific American, May 2024 — "The Secret to the Strongest Force in the Universe" by Brodsky, Deur & Roberts (doi:10.1038/scientificamerican0524-32). A vivid account of the convergence of three independent approaches to the QCD running coupling, written by the three principal researchers.
2303.00723 — QCD Running Couplings and Effective Charges (2023, the technical companion)

Profiles and databases

INSPIRE-HEP author page
Institute for Nonperturbative Physics, Nanjing University
Argonne profile (archived)