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.

◆

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.

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 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 Research Career

Craig Darrian Roberts is one of the most influential theoretical physicists working in nonperturbative quantum chromodynamics. Over a career spanning nearly four decades, he has transformed the Dyson–Schwinger equations from a formal theoretical apparatus into the premier quantitative tool for understanding hadron structure, dynamical mass generation, and the emergence of visible matter in the Universe.

Adelaide and the early years (1985–1992)

Roberts received his PhD in theoretical physics from the Flinders University of South Australia on 10 May 1988. His doctoral work, in collaboration with Anthony G. Williams (later at the University of Adelaide), laid the groundwork for what would become a career-defining programme: applying the Dyson–Schwinger equations of QCD to nonperturbative phenomena including confinement and dynamical chiral symmetry breaking. A key early insight was that confinement in quantum field theory is expressed through the absence of real-axis mass poles in the quark propagator — replaced instead by complex-conjugate singularities, meaning quarks cannot exist as free particles.

Following a postdoctoral fellowship at the University of Melbourne (1987–1989), Roberts joined the Theory Group in the Physics Division at Argonne National Laboratory in 1989 as a Research Associate. His rise was rapid: Assistant Physicist (1991–1996), then Physicist (1996–2006), then Senior Physicist (2006–2019).

Building the DSE programme at Argonne (1992–2000)

In 1994, Roberts and Williams published “Dyson–Schwinger Equations and their Application to Hadronic Physics” in Progress in Particle and Nuclear Physics. This review — with over 1,280 citations — remains the most-cited article in hadro-particle theory by an Australian-born scientist. It systematically showed how the tower of coupled integral equations governing QCD’s Green functions could be truncated in a symmetry-preserving manner and solved to yield predictions for confinement and DCSB. It effectively established the DSE programme as a serious competitor to lattice QCD for nonperturbative calculations.

During this period Roberts began his extraordinarily productive collaboration with Peter Maris (then at Argonne, later Iowa State) and Peter C. Tandy (Kent State University). The trio produced landmark calculations of pseudoscalar meson properties. Their 1997 papers on pion mass and decay constant and on pi- and K-meson Bethe–Salpeter amplitudes demonstrated that DSE methods could simultaneously preserve the Goldstone theorem, reproduce the Gell-Mann–Oakes–Renner relation as a corollary of the axial-vector Ward–Takahashi identity, and yield quantitative predictions for the pion and kaon in excellent agreement with experiment. Each of these papers accumulated over 500 citations and became cornerstones of the field.

In 2000, Roberts and Sebastian M. Schmidt (Forschungszentrum Jülich) published another major review extending DSE applications to finite temperature and density, covering the QCD phase diagram, deconfinement transitions, and the properties of quark-gluon plasma. This too garnered over 640 citations.

Theory Group Leader

In 2001, Roberts was appointed Leader of the Theory Group at Argonne — the youngest person ever to hold that position. He served until 2017, the longest tenure of any holder. Under his leadership the group became synonymous with DSE-based hadron physics worldwide.

The Brodsky collaboration and in-hadron condensates (2008–2015)

The collaboration with Stanley J. Brodsky (SLAC) proved transformative. In a series of influential papers beginning around 2010, Roberts and Brodsky challenged the conventional picture that chiral symmetry breaking condensates fill all of spacetime. Their work demonstrated that the chiral condensate is properly understood as a property of hadrons themselves — an “in-hadron condensate” — rather than a space-filling order parameter. This was further elaborated in “Confinement Contains Condensates” (2012), which argued that if quark-hadron duality is real in QCD, then all condensates are wholly contained within hadrons.

The Brodsky–Roberts collaboration also produced the landmark paper “Imaging Dynamical Chiral Symmetry Breaking: Pion Wave Function on the Light Front” (2013, Physical Review Letters), with Lei Chang, Ian Cloët, Schmidt, and Tandy. This projected the pion’s Bethe–Salpeter wave function onto the light front and showed that DCSB produces a pion distribution amplitude that is a broad, concave function — significantly different from both the asymptotic form and the narrow, end-point-concentrated forms assumed in some perturbative QCD analyses.

Precision predictions and the baryon programme (2009–2018)

Working with Si-xue Qin and Jorge Segovia, Roberts developed a comprehensive baryon physics programme using a Poincaré-covariant quark-diquark Faddeev equation approach. The key insight: the same DCSB mechanism generating constituent quark masses also produces strong, non-pointlike diquark correlations within baryons. This framework yielded predictions for nucleon and Delta elastic form factors, nucleon-to-Delta and nucleon-to-Roper electromagnetic transition form factors, and the spectra of all flavour-SU(3) octet and decuplet baryons.

With Qin, Roberts developed the influential “Interaction model for the gap equation” (2011), providing a momentum-dependent kernel consistent with both DSE and lattice QCD results that became a standard tool in the field.

The collaboration with Volker D. Burkert (Jefferson Lab) on the Roper resonance was a landmark achievement. For fifty years after its discovery in 1963, the proton’s first radial excitation had defied explanation. Their 2019 Reviews of Modern Physics Colloquium synthesised a decade of Jefferson Lab CLAS data with DSE-based calculations to show that the Roper is fundamentally a dressed-quark core augmented by a meson cloud that reduces the core mass by approximately 20%.

The QCD effective charge (2012–present)

Working with Daniele Binosi (ECT*, Trento), Joannis Papavassiliou, and José Rodríguez-Quintero, Roberts established the concept of a “process-independent effective charge” for QCD — an analogue of the Gell-Mann–Low effective coupling in QED. Their work demonstrated that this effective charge saturates at infrared momenta, reflecting dynamical breaking of scale invariance and the emergence of a gluon mass-scale. Using lattice QCD configurations at physical pion mass, they obtained a parameter-free prediction achieving near-perfect agreement with the independently measured Bjorken sum rule effective charge.

The collaboration with Alexandre Deur (Jefferson Lab) and Brodsky produced two major reviews of the QCD running coupling in Progress in Particle and Nuclear Physics (2016 and 2023), connecting perturbative QCD at high energies with nonperturbative DSE predictions in the infrared.

Nanjing University and Emergent Hadron Mass (2019–present)

Roberts moved to Nanjing University in September 2019 as International Distinguished Professor and Head of the newly established Institute for Nonperturbative Physics (INP). Here he synthesised decades of DSE research into the unifying conceptual framework of Emergent Hadron Mass (EHM).

The Higgs boson generates less than 2% of visible mass; EHM produces the remaining 98%. Roberts identified three “pillars” supporting EHM: (1) a running gluon mass-scale, generated dynamically through the Schwinger mechanism; (2) a process-independent effective charge that saturates in the infrared; and (3) dressed-quark running masses that take constituent-like values at infrared momenta. The comprehensive review “Emergence of Hadron Mass and Structure” (2023, with Minghui Ding and Schmidt) is the definitive statement of this framework. EHM has become a central organising concept for the physics programme of the Electron–Ion Collider (EIC) at Brookhaven National Laboratory.

At Nanjing, Roberts continues to build a world-leading group — working with Ding, Zhu-Fang Cui, Khepani Raya, and others on pion and kaon parton distributions, generalised parton distributions, gravitational form factors, and the prediction of observables measurable at the EIC. He was instrumental in defining the science case for both the US EIC (contributing to the Yellow Report, his single most-cited paper with over 1,470 citations) and the proposed Chinese Electron–Ion Collider (EicC).

By the numbers

Publication record

    Over 280 peer-reviewed articles • 24,700+ citations • h-index 89 •
    2 Reviews of Modern Physics articles •
    6 Progress in Particle and Nuclear Physics articles (the first four ranking
    2nd, 9th, 69th, and 110th most-cited in that journal’s 43-year history) •
    13 Physical Review Letters articles •
    400+ presentations including 160+ invited conference talks and lectures at 26 international
    graduate schools.
  

Awards and recognition

2001

Fellow, American Physical Society

2003

Friedrich Wilhelm Bessel Research Award, Alexander von Humboldt Foundation

2009

Convocation Medal, Flinders University

2012

International Fellow, Helmholtz Association

2014

Distinguished Performance Award, University of Chicago / Argonne LLC Board of Governors

2015

International Distinguished Professor, Chinese Ministry of Education

2022

Envoy of People’s Friendship, Jiangsu Province

2023

Friendship Award, Jiangsu Province

2025

Doctoratum Honoris Causa, Universidad de Huelva (Spain)

2025

International Science and Technology Cooperation Award, Nanjing University

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 260 peer-reviewed publications. 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.

Papers

Selected Publications

62 papers across 10 themes

All of Roberts’ papers are freely available on the arXiv preprint server. The following is a curated selection of his most significant works, organised by theme. Citation counts are approximate and drawn from INSPIRE-HEP.

Early Foundational DSE Work (1990s)

Dyson–Schwinger Equations and their Application to Hadronic Physics

hep-ph/9403224 · 1994 · Prog. Part. Nucl. Phys. 33· ~1,287 citations

Established the DSE formalism as a quantitative framework for nonperturbative QCD, covering confinement, DCSB, and hadronic applications — the most-cited paper in hadro-particle physics by an Australian-born scientist.

On the Implications of Confinement

INSPIRE · 1992 · Int. J. Mod. Phys. A 7· ~195 citations

Demonstrated that confinement is expressed through the absence of real-axis mass poles in quark propagators, replaced by complex-conjugate singularities.

Goldstone Theorem and Diquark Confinement Beyond Rainbow–Ladder Approximation

nucl-th/9602012 · 1996 · Phys. Lett. B 380· ~438 citations

Showed that repulsive corrections beyond rainbow-ladder truncation eliminate diquark bound-state poles while preserving the Goldstone theorem for pseudoscalar mesons.

Pi- and K-meson Bethe–Salpeter Amplitudes

nucl-th/9708029 · 1997 · Phys. Rev. C 56· ~645 citations

Derived the exact relation between the pseudoscalar meson Bethe–Salpeter amplitude and the dressed-quark propagator, showing the Gell-Mann–Oakes–Renner relation is a corollary of the axial-vector Ward–Takahashi identity.

Pion Mass and Decay Constant

nucl-th/9707003 · 1997 · Phys. Rev. C 56· ~523 citations

First fully self-consistent Poincaré-covariant DSE calculation of the pion mass and leptonic decay constant, establishing quantitative benchmarks for the rainbow-ladder truncation.

Survey of Heavy Meson Observables

nucl-th/9812063 · 1998 · Int. J. Mod. Phys. E 12· ~194 citations

Extended DSE meson calculations to heavy-quark systems, computing masses, decay constants, and form factors across the full flavour range from pion to upsilon.

Electromagnetic Pion Form-factor and Neutral Pion Decay Width

hep-ph/9408233 · 1996 · Nucl. Phys. A 605· ~198 citations

Calculated the pion electromagnetic form factor and π⁰→γγ decay width using confining DSE quark propagators and dressed vertices.

Reviews and Overviews

DSEs: Density, Temperature and Continuum Strong QCD

nucl-th/0005064 · 2000 · Prog. Part. Nucl. Phys. 45· ~646 citations

Comprehensive review extending DSE applications to finite temperature and density, covering the QCD phase diagram, deconfinement, and quark-gluon plasma.

Dyson–Schwinger Equations: A Tool for Hadron Physics

nucl-th/0301049 · 2003 · Int. J. Mod. Phys. E 12· ~676 citations

Pedagogical review demonstrating that DSEs furnish a Poincaré-covariant framework with symmetry-preserving truncations enabling proof of exact results in hadron physics.

Hadron Properties and Dyson–Schwinger Equations

0712.0633 · 2008 · Prog. Part. Nucl. Phys. 61· ~205 citations

Overview of DSE phenomenology covering pseudoscalar and vector mesons, nucleon elastic and transition form factors, and the role of DCSB.

Distribution Functions of the Nucleon and Pion in the Valence Region

1002.4666 · 2010 · Rev. Mod. Phys. 82· ~217 citations

Reviews of Modern Physics article providing a combined experimental and theoretical perspective on nucleon and pion parton distribution functions at large Bjorken-x.

Collective Perspective on Advances in DSE QCD

1201.3366 · 2012 · Commun. Theor. Phys. 58· ~364 citations

Multi-author survey of contemporary DSE studies covering confinement, DCSB, the hadron spectrum, form factors, parton distributions, and heavy quarks.

Explanation and Prediction of Observables using Continuum Strong QCD

1310.2651 · 2014 · Prog. Part. Nucl. Phys. 77· ~317 citations

Major review highlighting progress in reliable computation and prediction of measurable hadron properties, introducing the paradigm of in-hadron condensates.

The Pion: An Enigma within the Standard Model

1602.04016 · 2016 · J. Phys. G 43· ~185 citations

Comprehensive review of the pion’s dual nature as both a Nambu–Goldstone boson and a quark-antiquark composite, surveying theory and experiment.

Pion Physics

Imaging Dynamical Chiral Symmetry Breaking: Pion Wave Function on the Light Front

1301.0324 · 2013 · Phys. Rev. Lett. 110· ~279 citations

Demonstrated that DCSB produces a pion light-front distribution amplitude that is a broad, concave function at the hadronic scale, significantly different from the asymptotic form.

Pion Electromagnetic Form Factor at Spacelike Momenta

1307.0026 · 2013 · Phys. Rev. C· ~215 citations

Unified the prediction of the pion electromagnetic form factor with a broad PDA, showing that DCSB hardening effects reconcile DSE predictions with Jefferson Lab data.

Neutral Pion and its Electromagnetic Transition Form Factor

1510.02799 · 2016 · Phys. Rev. D 93· ~130 citations

Computed the γ*γ→π⁰ transition form factor on the entire spacelike domain, unifying it with the PDA and elastic form factor predictions.

Symmetry, Symmetry Breaking, and Pion Parton Distributions

1905.05208 · 2020 · Chin. Phys. C 44· ~124 citations

Showed that DCSB causes hardening of the pion’s valence-quark distribution function, with the computed large-x behaviour matching the (1−x)² expectation.

Pion and K-meson Electromagnetic Form Factors

nucl-th/9905056 · 2000 · Phys. Rev. C 62

Definitive rainbow-ladder calculation of pion and kaon electromagnetic form factors in impulse approximation, achieving excellent agreement with experimental data.

Abelian Anomaly and Neutral Pion Production

1009.0067 · 2011 · Few-Body Syst.· ~166 citations

Explored the connection between the axial anomaly and neutral pion electroproduction, including the role of DCSB in the γ*γ*π⁰ vertex.

Dynamical Chiral Symmetry Breaking

Essence of the Vacuum Quark Condensate

1005.4610 · 2010 · Phys. Rev. C 82· ~154 citations

Showed that the chiral condensate is a property of hadrons (“in-hadron condensate”) rather than a space-filling order parameter, challenging the conventional vacuum picture.

Confinement Contains Condensates

1202.2376 · 2012 · J. Phys. G 39· ~136 citations

Argued that if quark-hadron duality is real in QCD, then all condensates are wholly contained within hadrons and do not fill spacetime.

Dynamical Chiral Symmetry Breaking and a Critical Mass

nucl-th/0605058 · 2007 · Phys. Rev. C 75

Demonstrated the existence of a critical current-quark mass above which DCSB solutions are destabilised, bounding the domain of validity of the chiral expansion.

Sketching the Bethe–Salpeter Kernel

0903.5461 · 2009 · Phys. Rev. Lett. 103· ~277 citations

Constructed the first Bethe–Salpeter kernel consistent with the gap equation kernel and capable of supporting lattice-QCD-verified DCSB.

Aspects and Consequences of a Dressed Quark–Gluon Vertex

nucl-th/0403012 · 2004 · Phys. Rev. C 70· ~207 citations

Investigated the impact of a dressed quark-gluon vertex on DCSB and meson properties, comparing DSE predictions with lattice data for the quark mass function.

Dynamical Chiral Symmetry Breaking and the Fermion–Gauge-Boson Vertex

1112.4847 · 2011 · Phys. Rev. C 84· ~103 citations

Showed that the dressed fermion–gauge-boson vertex must contain structures generated by DCSB that are essential for maintaining Ward–Takahashi identities.

Emergent Hadron Mass

Emergence of Hadron Mass and Structure

2211.07763 · 2023 · Particles 6· ~106 citations

Definitive review of the EHM framework, presenting the three pillars (running gluon mass, process-independent effective charge, running quark mass) and their measurable consequences.

Insights into Emergence of Mass from Pion and Kaon Structure

2102.01765 · 2021 · Prog. Part. Nucl. Phys. 120· ~202 citations

Reviewed the role of EHM in pion and kaon properties and laid out the experimental programme needed to validate EHM predictions.

Reflections upon the Emergence of Hadronic Mass

2006.08782 · 2020 · Eur. Phys. J. Spec. Top. 229

Discussed how the Higgs boson generates less than 2% of visible mass and how EHM, operating through DCSB and confinement, produces the remainder.

Empirical Consequences of Emergent Mass

2009.04011 · 2020 · Symmetry 12· ~93 citations

Articulated the experimental signatures of EHM measurable at existing and planned facilities. Winner of the Symmetry 2020 Best Paper Award.

On Mass and Matter

2101.08340 · 2021 · AAPPS Bull. 31

Concise overview of how the mass of visible matter emerges from QCD dynamics rather than the Higgs mechanism.

Perspective on the Origin of Hadron Masses

1606.03909 · 2017 · Few-Body Syst. 58· ~93 citations

Articulated the case that understanding hadron mass generation requires experiments at the EIC, JLab 12 GeV, and FAIR.

Hadron Structure: Perspective and Insights

2503.05984 · 2025

Recent review highlighting EHM’s role in shaping electromagnetic and gravitational form factors, nucleon resonances, and hadron parton distributions.

Insights into Meson and Baryon Structure using Continuum Schwinger Function Methods

2601.08046 · 2026

Latest unified explanation of pion and proton electromagnetic and gravitational form factors through the EHM framework.

Baryon Structure and Spectra

Colloquium: Roper Resonance — Toward a Solution to the Fifty Year Puzzle

1710.02549 · 2019 · Rev. Mod. Phys. 91· ~154 citations

Reviews of Modern Physics colloquium synthesising a decade of JLab data and DSE theory to explain the Roper resonance as a dressed-quark core augmented by a meson cloud.

Completing the Picture of the Roper Resonance

1504.04386 · 2015 · Phys. Rev. Lett. 115· ~142 citations

Showed that the Roper resonance is the proton’s first radial excitation with DCSB-generated diquark correlations, resolving a key feature of the fifty-year puzzle.

Toward Unifying the Description of Meson and Baryon Properties

0810.1222 · 2009 · Phys. Rev. C 79· ~159 citations

First unified DSE treatment of both mesons and baryons using identical interaction kernels, demonstrating consistency between the two sectors.

Nucleon and Δ Elastic and Transition Form Factors

1408.2919 · 2014 · Few-Body Syst. 55· ~132 citations

Comprehensive quark-diquark Faddeev equation calculations of nucleon and Delta elastic form factors and the N→Δ electromagnetic transition.

Spectrum of Hadrons with Strangeness

1204.2553 · 2013 · Phys. Rev. C 87· ~108 citations

Extended the quark-diquark Faddeev equation approach to strange baryons, computing masses of octet and decuplet states with strangeness.

Studies of Nucleon Resonance Structure in Exclusive Meson Electroproduction

1212.4891 · 2013 · Int. J. Mod. Phys. E 22· ~251 citations

White paper on the JLab/CLAS programme for nucleon resonance electroexcitation, connecting measurements to DSE predictions.

Spectrum of Fully-heavy Tetraquarks from a Diquark+Antidiquark Perspective

1911.00960 · 2020 · Eur. Phys. J. C 80· ~179 citations

Used a relativised diquark model to predict masses of fully-heavy tetraquark systems including cccc, bbbb, and mixed-flavour states.

Nucleon Form Factors

Survey of Nucleon Electromagnetic Form Factors

0812.0416 · 2009 · Few-Body Syst. 46· ~167 citations

Comprehensive survey including flavour separation of form factor contributions and prediction of a zero in the proton electric-to-magnetic form factor ratio.

Nucleon Electromagnetic Form Factors

nucl-th/0611050 · 2006 · Phys. Rev. C· ~241 citations

Calculated proton and neutron electromagnetic form factors with DSE-dressed quarks, predicting a zero in G_E^p/G_M^p at large Q².

Dressed-quark Anomalous Magnetic Moments

1009.3458 · 2013 · Phys. Rev. Lett. 111· ~206 citations

Showed that DCSB generates large dressed-quark anomalous magnetic moments that play a critical role in nucleon electromagnetic form factors.

Parton Distributions

Pion and Kaon Valence-quark Parton Distribution Functions

1102.2448 · 2011 · Phys. Rev. C 83· ~103 citations

First DSE calculation of pion and kaon valence-quark PDFs, with the ratio of kaon to pion u-quark distributions agreeing with Drell–Yan data.

Valence-quark Distribution Functions in the Pion

1602.01502 · 2016 · Phys. Rev. D 93· ~119 citations

Refined calculation of pion and kaon valence-quark distributions incorporating realistic momentum-dependent DSE kernels.

Kaon and Pion Parton Distributions

2006.14075 · 2020 · Eur. Phys. J. C 80· ~118 citations

First continuum QCD prediction of kaon parton distributions including all valence, glue, and sea contributions, revealing Higgs-boson modulation of EHM.

Revealing Pion and Kaon Structure via Generalised Parton Distributions

2109.11686 · 2022 · Chin. Phys. C 46· ~124 citations

First continuum QCD predictions for pion and kaon GPDs constrained by hadron-scale valence PDFs, providing windows onto EHM.

QCD Running Coupling

The QCD Running Coupling

1604.08082 · 2016 · Prog. Part. Nucl. Phys. 90

Major review of α_s across all momentum scales, connecting perturbative QCD at high energies with nonperturbative DSE and light-front holography predictions in the infrared.

QCD Running Couplings and Effective Charges

2303.00723 · 2023 · Prog. Part. Nucl. Phys. 134· ~90 citations

Updated comprehensive review covering process-independent effective charges, lattice QCD determinations, and connections to EHM.

Process-independent Strong Running Coupling

1612.04835 · 2017 · Phys. Rev. D 96· ~168 citations

Calculated a process-independent QCD effective charge analogous to the Gell-Mann–Low coupling, finding near-perfect agreement with the Bjorken sum rule effective charge.

Effective Charge from Lattice QCD

1912.08232 · 2020 · Chin. Phys. C 44· ~137 citations

Used lattice QCD configurations at physical pion mass to obtain a parameter-free prediction of QCD’s process-independent effective charge.

Confinement and the QCD Phase Diagram

Bridging a Gap between Continuum-QCD and Ab Initio Predictions of Hadron Observables

1412.4782 · 2015 · Phys. Lett. B 742· ~234 citations

Demonstrated that the renormalisation-group-invariant running interaction predicted by gauge-sector analyses coincides with that needed to describe ground-state hadron observables.

Phase Diagram and Critical Endpoint for Strongly-interacting Quarks

1011.2876 · 2011 · Phys. Rev. Lett. 106· ~230 citations

Located the critical endpoint in the QCD phase diagram using chiral susceptibility methods within the DSE framework.

Locating the Gribov Horizon

1706.04681 · 2018 · Phys. Rev. D 97· ~91 citations

Explored whether a tree-level gluon propagator is consistent with lattice QCD and located the boundary of the fundamental modular region of gauge-field configuration space.

Interaction Model for the Gap Equation

1108.0603 · 2011 · Phys. Rev. C 84· ~250 citations

Introduced a momentum-dependent interaction kernel for the DSE gap equation consistent with modern DSE and lattice QCD results, becoming a standard tool in the field.

Natural Constraints on the Gluon–Quark Vertex

1609.02568 · 2018 · Phys. Rev. C 97· ~118 citations

Identified constraints on the dressed gluon-quark vertex imposed by multiplicative renormalisability and BRST symmetry, necessary for reliable DSE truncations.

Analysis of a Quenched Lattice-QCD Dressed-quark Propagator

nucl-th/0304003 · 2003 · Phys. Rev. C 68· ~237 citations

Correlated quenched lattice QCD data on the dressed-quark propagator with the rainbow gap equation, identifying the need for vertex dressing.

Electron–Ion Collider Science Case

Science Requirements and Detector Concepts for the Electron–Ion Collider: EIC Yellow Report

2103.05419 · 2022 · Nucl. Phys. A 1026· ~1,474 citations

Comprehensive community white paper defining the physics case, detector requirements, and experimental programme for the US Electron–Ion Collider. Roberts’ single most-cited paper.

Electron–Ion Collider in China

2102.09222 · 2021 · Front. Phys. 16· ~554 citations

White paper for the proposed Chinese Electron–Ion Collider (EicC), covering science opportunities in hadron structure, EHM, and nuclear physics.

Pion and Kaon Structure at the Electron–Ion Collider

1907.08218 · 2019 · Eur. Phys. J. A 55· ~202 citations

Laid out the science case for measuring pion and kaon structure at the EIC, connecting the mass budget of pseudoscalar mesons to five key experimental measurements.

Revealing the Structure of Light Pseudoscalar Mesons at the Electron–Ion Collider

2102.11788 · 2021 · J. Phys. G 48· ~124 citations

Detailed programme for measuring pion and kaon electromagnetic form factors, PDFs, and GPDs at the EIC to expose EHM mechanisms.

Profiles and Databases

INSPIRE-HEP author page — complete publication list with citation metrics
Google Scholar profile
Institute for Nonperturbative Physics, Nanjing University
Personal website
Argonne profile (archived)