Empirical Consequences of Emergent Mass

The Lagrangian that defines quantum chromodynamics (QCD), the strong interaction piece of the Standard Model, appears very simple. Nevertheless, it is responsible for an astonishing array of high-level phenomena with enormous apparent complexity, e.g., the existence, number and structure of atomic nuclei. The source of all these things can be traced to emergent mass, which might itself be QCD's self-stabilising mechanism. A background to this perspective is provided, presenting, inter alia, a discussion of the gluon mass and QCD's process-independent effective charge and highlighting an array of observable expressions of emergent mass, ranging from its manifestations in pion parton distributions to those in nucleon electromagnetic form factors.


Introduction
One might define emergent phenomena as those features of nature which do not readily admit an explanation solely in terms of known or conjectured mathematical rules. The concept is at least as old as Aristotle (384-322 BC), who argued that a compound item can have (emergent) properties in the whole which are not explicable merely through the independent actions of the item's constituent parts. His view is often represented by the statement "The whole is more than the sum of its parts." In this sense, emergence has its origins in the Greek "sunergos": "together" plus "working", which is the origin of our current concept of synergy, viz. things working together more effectively than could be anticipated from their independent actions in isolation. (Etymologically, the word "emergence" entered English in the mid 17th century, meaning "unforeseen occurrence" and derived from the medieval Latin "emergentia", itself from the Latin "emergere", meaning "bring to light" (Source: Oxford English Dictionary). Herein, "emergence" is seen as a larger notion than is typically expected from this literal connection.).
This perspective is typically contrasted with that described as reductionism; namely, the view that everything in nature can ultimately be viewed as no more complex in principle than, e.g., a (very good) watch, which is clearly a complex object; but, equally clearly, not more than the sum of its parts.
In developing such a contrast through debate, hard lines are sometimes drawn, with individuals deciding or being forced to choose one side or the other. This is taking the argument too far, however. During each epoch in history, there has always been a line dividing physics from metaphysics; but the location of that border is neither fixed nor impermeable. As time and humanity have progressed, more aspects of nature have seeped into the pool of physics. Notwithstanding such progress, it does not follow that fewer questions have been left beyond physics. Typically, as mathematics has succeeded in explaining more phenomena, new discoveries have been made, often emerging from attempts to test the These remarks implicitly highlight QCD, the quantum field theory formulated in four spacetime dimensions which defines what is arguably the SM's most important chapter. QCD is supposed to describe all of nuclear physics through the interactions between quarks (matter fields) that are mediated by gluons (gauge bosons). Yet, fifty years after the discovery of quarks [6][7][8], science is only just beginning to grasp how QCD moulds the basic bricks for nuclei: pions, neutrons, protons, etc., and it is far from understanding how QCD produces nuclei.
The natural scale for nuclear physics (strong interactions) is characterised by the proton mass: where m e is the electron mass, i.e., m p = 1.783 × 10 −27 kg. In the SM, m e is correctly attributed to the Higgs boson, but what is the cause of the prodigious enhancement required to produce m p ? Followed logically to its source, this question leads to an appreciation that our Universe's existence depends critically on, inter alia, the following empirical facts. (i) The proton mass is large, i.e., the mass-scale for QCD is very much larger than that of electromagnetism; (ii) the proton does not decay, despite being a compound state built from three valence-quarks; and (iii) the pion, which carries long-range interactions between nucleons (neutrons and protons), is abnormally light (not massless), having a lepton-like mass even though it is a strongly interacting object built from a valence-quark and valence antiquark. (The µ-lepton, discovered in 1936 [9], was initially mistaken for the pion. The pion was only found a decade later [10].) These qualities of Nature transport us to a reductionist definition of emergence. Namely, assuming it does describe strong interactions, then the one-line Lagrangian of QCD-a very simple low-level rule-must somehow produce high-level phenomena with enormous apparent complexity. At this point, it is worth studying the Lagrangian of chromodynamics, which appeared as the culmination of a distillation process applied to a large array of distinct ideas and discoveries [11,12]: Here, {q j } are the quark fields, with j their flavour label and m j their Higgs-generated current-quark masses, and {A a µ , a = 1, . . . , 8} are the gluon fields, with { 1 2 λ a } the generators of the SU(3) (chromo/colour) gauge-group in the fundamental representation. Comparing with quantum electrodynamics (QED), the solitary difference is the piece describing gluon self-interactions, marked as the underlined term in Equation (2b). Somehow, the origin, mass and extent of almost all visible matter in the Universe is attributable to L QCD -one line plus two definitions. That being true, then . . . QCD is quite possibly the most remarkable fundamental theory ever invented.
The only apparent energy scales in Equation (2) are the Higgs-generated current-quark masses, but focusing on the u (up) and d (down) quarks that define nucleons, this scale is more-than one-hundred-times smaller than m p [13]. No amount of "staring" at L QCD can reveal the source of that enormous amount of "missing mass"; yet, it must be there. (This is a stark contrast to QED wherein, e.g., the scale in the spectrum of the hydrogen atom is set by m e , which is a prominent feature of L QED that is generated by the Higgs boson.).
Models and effective field theories (EFTs) for nuclear physics typically assume existence of the m p ≈ 1 GeV mass-scale and build upon it. They also assume the reality of effectively pointlike nuclear constituents (nucleons) and force carriers (pions and, perhaps, other meson-like entities). Their task is not to elucidate the internal structure of such objects. Instead, they aim to develop systematically improvable techniques that can describe the number and properties of atomic nuclei. This may be seen as reductionism built on an emergent plateau. The basic reductionist question generated by this approach is "Can the plateau upon which the nuclear model/EFT paradigm is built be constructed from QCD?" If the answer is "yes", then all parameters used and fitted in such theories of nuclear structure will (some day) be confronted with ab initio predictions in a profound test of the SM.
Today there is a good case to be made for an affirmative answer to this question. Amongst the numerous supporting examples, one may list QCD-connected computations of the hadron spectrum. Depicted in Figure 1, one sees that two disparate approaches to solving QCD [14,15] produce a spectrum of ground-state hadrons in good agreement with experiment.  [14] and lattice (lQCD-circles, blue) [15] methods in QCD compared with experiment [13] (PDG-black bars, with decay-widths of unstable states shaded in grey). The continuum study did not include isospin symmetry breaking effects, which are evidently small, as highlighted by the empirically determined Σ-Λ mass difference (<7%). Importantly, however, neither calculation presented in Figure 1 was able to predict the size of the proton mass. Apart from the Higgs-generated current-quark masses, each has an undetermined mass-scale parameter, denoted hereafter as Λ QCD . Its value is chosen to fit one experimental number, e.g., the pion's leptonic decay constant, but after Λ QCD is fixed, all other results are predictions.

Whence Mass?
Herein, it is worth explaining the need for such a scale-setting procedure. Regarded classically, chromodynamics is a local non-Abelian gauge field theory, and defined in four spacetime dimensions, there is no mass-scale if Lagrangian masses for the fermions are omitted. (The absence of such masses defines the chiral limit.) Scale invariant theories do not support dynamics, only kinematics. Therefore, bound states are impossible; so, the Universe cannot exist.
Spontaneous symmetry breaking, as generated by the Higgs mechanism, does not resolve this issue. The masses of the neutron and proton, which lie at the heart of all visible matter, are two orders-of-magnitude greater than the Higgs-generated current-masses of the uand d-quarks, which are the defining constituents of protons and neutrons.
Consequently, questions like "How does a mass-scale appear?" and "Why does this scale have its observed value?" are consanguineous with the question "How did the Universe form?" Modern quantum field theories are built upon Poincaré invariance. In this connection, consider the energy-momentum tensor in classical chromodynamics, T µν . In field theory, conservation of energy and momentum is a consequence of spacetime translational invariance, one of the family of Poincaré transformations. Consequently, Consider now a global scale transformation in the Lagrangian of classical chromodynamics: where σ is spacetime-independent. The associated Noether current is viz. the dilation current. In the absence of fermion masses, the classical action is invariant under Equation (4), i.e., the theory is scale invariant; therefore, where the last equality follows from Equation (3). Evidently, the energy-momentum tensor must be traceless in a scale invariant theory. Classical chromodynamics is not a meaningful physical framework for many reasons; amongst them the fact, illustrated by Figure 1, that strong interactions are empirically known to be characterised by a large mass-scale, m p ≈ 1 GeV. In quantising the theory, a mass-scale is introduced by regularisation and renormalisation of ultraviolet divergences. This is "dimensional transmutation": all quantities, including the field operators themselves, become dependent on a mass-scale. This entails the violation of Equation (6), i.e., the appearance of the chiral-limit "trace anomaly" [16]: where β(α(ζ)) is QCD's β-function, α(ζ) is the associated running-coupling, and ζ is the renormalisation scale. Equation (7) indicates that a mass-scale linked to the resolving power of a given measurement is engendered via quantisation; to wit, the scale emerges as an indispensable part of the theory's quantum definition. This mass is exhibited in the gauge-boson vacuum polarization. In QED, the photon vacuum polarization does not possess an infrared mass-scale, and dimensional transmutation acts simply to generate the gentle evolution of the QED coupling, i.e., any dynamical breaking of QED's conformal features is small; therefore, the associated trace anomaly is normally negligible. In contrast, gauge sector dynamics drives a Schwinger mechanism in QCD [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34], so that the QCD trace anomaly expresses a mass-scale which is, empirically, very significant. This is discussed in Section 4.2.
There is another aspect of chromodynamics that should be highlighted, namely, the classical Lagrangian still defines a non-Abelian local gauge theory. Accordingly, the concept of local gauge invariance persists; but without a mass-scale, there is no notion of confinement. For instance, one can compose a colour-singlet combination of three quarks and colour rotations will preserve its colour neutrality; but the participating quarks need not be close together. In fact, it is meaningless to discuss proximity because, in a scale invariant theory, all lengths are equivalent. Accordingly, the question of "Whence mass?" is indistinguishable from "Whence a mass-scale?", is indistinguishable from "Whence a confinement scale?".
Evidently, one does not learn much from knowing that a trace anomaly exists. It only means there is a mass-scale. The central concern is whether the magnitude of that scale can be computed and understood.
The magnitude of the scale anomaly can definitely be measured, and this simply by considering the in-proton expectation value of the energy-momentum tensor (see, e.g., in [35]): where the equations-of-motion for a one-particle proton state were used to obtain the right-hand-side. Now, in the chiral limit p(P)|T µµ |p(P) = −P 2 = m 2 p = p(P)|Θ 0 |p(P) .
Therefore, there is a sound position from which one may conclude that gluons generate all the proton's mass: the measured value of the trace anomaly is large; and, logically, that feature owes to gluon self-interactions, the agent behind asymptotic freedom. This is a valid deduction because, ultimately, what else could account for a mass-scale in QCD? Gluon self-interactions are QCD's definitive feature, and it is these interactions that might enable a rigorous (non-perturbative) definition of the matrix element in Equation (9). Nevertheless, it is only reasonable to conclude this when the operator and wave function are computed at a resolving scale ζ m p , i.e., when one employs a parton-model basis [36].
There is another issue, too, which can be exposed by returning to Equation (8) and replacing the proton by the pion because the chiral-limit pion is a massless Nambu-Goldstone (NG) mode [37,38]. It is possible that Equation (10) means the scale anomaly is trivially zero in the pion, i.e., gluon self-interactions have no effect in the pion because each term required to express the operator vanishes separately. Yet, such a conclusion would sit uncomfortably with known QCD dynamics, which expresses both attraction and repulsion, but never passive inactivity. More likely, then, the final identity in Equation (10) results from cancellations between different terms in the complete operator matrix element. Naturally, such precise cancellation could not be accidental. It would require that some symmetry is broken in a very particular manner. (The mechanism is explained in Section 4.5.) Equations (9) and (10) present a quandary. They stress that any understanding of the proton's mass is incomplete unless it simultaneously explains Equation (10). Moreover, any discussion of confinement, fundamental to the proton's stability, is unreasonable before this conundrum is resolved. As will become clear, at least some of these features of Nature have a reductionist explanation grounded in the dynamics responsible for the emergence of m p as the natural mass-scale for nuclear physics; and one of the most important goals in modern science is to explain and elucidate the entire array of empirical consequences of this dynamics [39][40][41][42][43].

Confinement
It is a textbook result, with its origin in the Nobel Prize for the discovery of asymptotic freedom [1,44,45], that QCD is characterised by an interaction which becomes stronger as the participants try to separate. Remaining within those results that can be established using perturbation theory, one is led to contemplate some unusual possibilities: if the coupling strength rises rapidly with separation, then perhaps it is not bounded, and perhaps an infinite amount of energy is required to remove a gluon or quark from the interior of a hadron? Such thinking has produced The Confinement Hypothesis: Colour-charged particles cannot be isolated and therefore cannot be directly observed. They clump together in colour-neutral bound-states.
Confinement seems to be an empirical reality, but a credited mathematical proof is lacking. Partly as a result, the Clay Mathematics Institute proffered a "Millennium Problem" prize of $1 million for a proof that pure-glue QCD is mathematically well defined [46]. One necessary part of such a proof would be to establish whether the confinement postulate is correct in pure-gauge QCD i.e., the theory obtained from Equation (2) after omitting all terms containing quark fields.
There is a pitfall here, however: no reader of this material can be described within pure-glue QCD. Light quarks are essential to understanding all known visible matter. Thus, a proof of confinement in pure-glue QCD is chiefly irrelevant to our Universe. Life exists because nature has supplied two light quark flavours, combinations of which form the pion, and the pion is unnaturally light, thus very easily produced and capable of propagating over nuclear-size distances. Therefore, as noted previously by others [47,48], no explanation of SM confinement is empirically relevant unless it also describes the link between confinement and the emergence of mass, and so the existence and role of pions, i.e., pseudo-NG modes with m π m p . One piece of the Yang-Mills millennium problem [46] is to prove that a mass-gap, ∆ > 0, exists in pure-glue QCD. This conjecture is supported by some strong evidence, e.g., numerical studies of lattice-regularised QCD (lQCD) find ∆ 1.5 GeV [49]. This sharpens the conundrum described above: with ∆ 2 /m 2 π 100, can the mass-gap in pure Yang-Mills theory really play any role in understanding confinement when the emergence of mass, driven by kindred dynamics, ensures that our Universe supports almost-massless strongly-interacting excitations? Skirting the question, one can respond that any explanation of confinement must simultaneously describe its link to pion properties. From this position, pions are viewed as playing a critical role in any explanation of SM confinement, and a discussion that omits reference to pions is practically irrelevant.
These observations indicate that the potential between infinitely-heavy quarks computed in numerical simulations of quenched lQCD-the static potential [50], often associated with formation of an incredibly strong flux tube between the colour source and sink [51]-is detached from the question of confinement in our Universe. In fact, as light-particle annihilation and creation effects are essentially non-perturbative in QCD, it is impossible to calculate a quantum mechanical potential between two light quarks [52][53][54]. It follows that there is no discernible flux tube in a Universe with light quarks; so, the flux tube is not the correct paradigm for confinement.
As highlighted already, the emergence of mass is key here. It ensures the existence of pseudo-NG modes, and no flux tube linking a static colour source and sink can have an observable existence in the presence of these modes. To verify this statement, suppose that such a tube is stretched between a source and sink. The potential energy stored within the tube may only increase until it reaches the amount required to produce a particle-antiparticle pair of the system's pseudo-NG modes. Simulations of lQCD demonstrate [52,53] that the flux tube then disappears instantaneously, leaving two separated colour-singlet systems. The length-scale characterising this effect is r σ (1/3) fm; so, if any such string forms, it would disintegrate well within a hadron's interior.
It should be highlighted, however, that the appearance of such modifications when analysing some simplification of a given theory does not signify that the theory itself is truly confining: uncommon spectral properties can be introduced by approximations, yielding a truncated version of a theory which expresses confinement even though the complete theory does not, see, e.g., in [72,73]. Notwithstanding exceptions like these, in a veracious treatment of QCD the computed violation of reflection positivity by coloured functions does express confinement.

Dyson-Schwinger Equations
The appearance and size of the natural scale for nuclear physics (m p ≈ 1 GeV) and the confinement of gluons and quarks are emergent phenomena. They are not apparent in the QCD Lagrangian, yet they determine the character of QCD's spectrum, the structure of bound states and so forth. Given Equation (2), the natural question to ask is whether one can understand these features reductively, i.e., directly in terms of the degrees-of-freedom used to formulate QCD, or does the complexity of strong interaction phenomena make prediction and explanation impractical (impossible)? For instance, is it pointless to attempt the QCD-connected prediction of any nucleon structural property on a domain that is not yet empirically accessible?
If a reductive explanation is impossible, then science must rely on a tower of EFTs, each level developed for a different energy domain, in order, e.g., to express and understand the consequences of the emergence of mass and contingent effects, such as confinement, without identifying their source.
On the other hand, if a reductive approach is possible, then non-perturbative calculational methods must be developed to define and solve QCD. Prominent amongst such techniques today are (i) the numerical simulation of lQCD [74][75][76] and (ii) continuum Schwinger function methods (CSMs), viz.ã collection of models and schemes, each with varying degrees of connection to Equation (2). Currently, each of these two approaches has strengths and weaknesses, so the best way forward is to combine them to the fullest extent that is reasonably possible and exploit the synergies that emerge.
Among CSMs, the Dyson-Schwinger equations (DSEs) have proven useful [63,[77][78][79][80][81][82][83][84][85][86]. These quantum field theory generalisations of the Euler-Lagrange equations provide a continuum method for calculating Schwinger functions; namely, the same Euclidean-space Green functions that are computed using lQCD. Consequently, there are many opportunities for cross-fertilisation between DSE and lQCD studies, and this has been exploited to increasingly good effect during the past twenty years, especially at the level of propagators (2-point functions) and vertices (3-point functions) defined using QCD's elementary degrees-of-freedom.
The challenge to DSE analyses is found in the fact that the equation for any given n-point function is coupled with those for some higher-n-point functions, e.g., the gap equation for the quark 2-point function is coupled to those for the gluon 2-point function and the gluon-quark 3-point function. Therefore, truncations are necessary in order to define a tractable problem. Systematic, symmetry-preserving schemes have been developed [87][88][89][90][91][92][93], so that today DSE predictions can be distributed into three classes: (A) model-independent results in QCD; (B) illustrations of such results, using well-constrained model elements and possessing a recognisable connection to QCD; and (C) analyses that can reasonably be described as QCD-based, but whose elements have not been calculated using a truncation that maintains a systematically-improvable connection with QCD. Regarding the last two classes, comparisons between schemes and orders within schemes can be used to identify robust outcomes. Results can also be compared with lQCD predictions, when available, capitalising on the overlap domain of these two approaches, and, of course, predictions can be tested against experiment, which is, as always, the final arbiter in physics.

Gluon Mass
It is now possible and appropriate to return to the confinement hypothesis, introduced in the opening paragraph of Section 3 and expressed in the "Millennium Problem". Beginning with pioneering efforts roughly forty years ago [17,18], continuum and lattice studies of QCD's gauge sector have been growing in sophistication and reliability. The current state of understanding can be traced from an array of sources, see, e.g., in [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34]. Of specific interest is the property that the gluon propagator saturates at infrared momenta; to wit, where ∆(k 2 ) is the scalar function that characterises the dressed gluon propagator. This entails that the long-range propagation characteristics of gluons are markedly altered by their self-interactions. Significantly, one may associate a renormalisation-group-invariant (RGI) gluon mass-scale with this effect: m 0 ≈ m p /2, and summarise a large body of work by recording that gluons, although behaving as massless entities on the perturbative domain, actually possess a running mass, with m 0 characterising its value at infrared momenta. Asymptotic freedom guarantees that QCD's ultraviolet behaviour is tractable; but the emergence of a gluon mass reveals a new SM physics frontier because the existence of a running gluon mass, sizeable at infrared momenta, influences all analyses of the continuum bound-state problem. For instance, it could be a harbinger of gluon saturation [94,95].
Furthermore, m 0 > 0 entails that QCD dynamically generates its own infrared cut-off, so that gluons with wavelength λ σ := 1/m 0 ≈ 0.5 fm decouple from the strong interaction, hinting at a dynamical realisation of confinement [25,31,[56][57][58][59][60][61][62][63][64][65][66][67][68][69][70][71]. In this picture, once a gluon or quark is produced, it starts to propagate in spacetime, but following each "step" of average length σ an interaction occurs and the parton loses its identity, sharing it with others. Ultimately a parton cloud is produced, which fuses into colour-singlet final states. This physics is embodied in parton fragmentation functions (PFFs), which describe how QCD partons, (nearly) massless when produced in a high-energy event, transform into a cascade of massive hadrons. PFFs express the emergence of hadrons with mass from massless partons [96]. Such observations suggest that PFFs are the cleanest expression of dynamical confinement in QCD. This perspective can be explored at modern and anticipated facilities.

Effective Charge
Among the many other consequences of QCD's intricate non-perturbative gauge-sector dynamics is the generation of a process-independent (PI) running coupling,α(k 2 ), see, e.g., in [32][33][34]. Depicted as the solid (black) curve in Figure 2, this is a new type of effective charge. It is an analogue of the Gell-Mann-Low effective coupling in QED [97] because it is completely determined by the gauge-boson vacuum polarisation, when the problem is approached using the pinch technique [17,[98][99][100] and background field method [101]. The result in Figure 2 is a parameter-free Class-A DSE prediction, capitalising on analyses of QCD's gauge sector undertaken using continuum methods and informed by numerical simulations of lQCD.  [33]. (The grey band bordered by dashed curves indicates the uncertainty in the result arising from that in both continuum and lattice-QCD inputs and is detailed in [34].) For comparison, world data on the process-dependent charge, α g 1 , defined via the Bjorken sum rule, are also depicted. (The data sources are listed elsewhere [34]. For additional details, see e.g., in [102][103][104].) The k-axis scale is linear to the left of the vertical partition, and logarithmic otherwise. The vertical line, k = m G , marks the gauge sector screening mass, Equation (12).
The data in Figure 2 represent empirical information on α g 1 , a process-dependent effective-charge [105,106] determined from the Bjorken sum rule, a basic constraint on our knowledge of nucleon spin structure. Solid theoretical reasons underpin the almost precise agreement betweenα and α g 1 [32][33][34], so that the Bjorken sum serves as a window through which to gain empirical insight into QCD's effective charge. Figure 2 shows that QCD's unique effective coupling is everywhere finite, i.e., there is no Landau pole and the theory plausibly possesses an infrared-stable fixed point. Apparently, QCD is infrared-finite because a gluon mass-scale is dynamically generated. (A theory is said to possess a Landau pole at k 2 L if the effective charge diverges at that point. In QCD perturbation theory, such a pole exists at k 2 L = Λ 2 QCD . Were such a pole to persist in a complete treatment of QCD, it would signal an infrared failure of the theory. On the other hand, the absence of a Landau pole supports a view that QCD is alone amongst four-dimensional quantum field theories in being defined and internally consistent at all energy scales. This might have implications for attempts to develop an understanding of physics beyond the SM based upon non-Abelian gauge theories [71,[107][108][109][110][111][112][113]). In this case, the value of the PI charge at k 2 = Λ 2 QCD defines a screening mass [114,115]: m G ≈ 1.4 Λ QCD ≈ 0.71 GeV. As evident in Figure 2, m G marks a boundary: the running coupling alters character at k m G so that modes with k 2 m 2 G are screened from interactions and the theory enters a practically conformal domain. Evidently, the line k = m G draws a natural border between soft and hard physics; therefore, defines the "hadronic scale": This is the scale at which all the properties of a hadron are expressed by the dressed quasiparticles that form the DSE kernels and emerge as the self-consistent solutions of the associated equations.
As a unique PI effective charge,α appears in each of QCD's dynamical equations of motion, including the gauge-and matter-sector gap equations, setting the strength of all interactions. It therefore plays a critical role in settling the fate of chiral symmetry; to wit, the dynamical origin of light-quark masses in the SM even in the absence of a Higgs coupling.

Dynamical Chiral Symmetry Breaking
Just as a gluon mass-scale emerges dynamically in QCD, massless current-quarks become massive dressed-quarks through a phenomenon known as dynamical chiral symmetry breaking (DCSB) [116]. This effect is another of the critical emergent phenomena in QCD. It is expressed in hadron wave functions, not in vacuum condensates [81,[117][118][119][120], and modern theory indicates that more than 98% of the visible mass in the Universe can be attributed to DCSB. As classical massless-QCD is a scale-invariant theory (Section 2.2), this means that DCSB is fundamentally connected with the origin of mass from nothing.
DCSB is most readily apparent in the dressed-quark propagator which is obtained as the solution of a gap equation whose kernel is critically dependent uponα. M(p 2 ) in Equation (13) is the dressed-quark mass-function, depicted and explained in Figure 3.  (13): solid curves-gap equation solutions [121,122], "data"-numerical simulations of lQCD [123], available for current-quark masses m = 30, 70 MeV. QCD's current-quark evolves into a constituent-quark as its momentum becomes smaller. The constituent-quark mass arises from a cloud of low-momentum gluons attaching themselves to the current-quark. This is DCSB, the essentially non-perturbative effect that generates a quark mass from nothing; namely, it occurs even in the chiral limit. Notably, the size of M(0) is a measure of the magnitude of the QCD scale anomaly in n = 1-point Schwinger functions [124]. Moreover, experiments on Q 2 ∈ [0, 12] GeV 2 at the modern Thomas Jefferson National Accelerator Facility (JLab) will be sensitive to the momentum dependence of M(p) within a domain that is here indicated approximately by the shaded region.
One must insist that chiral symmetry breaking in the absence of a Higgs mechanism is dynamical. It is distinct from Higgs-induced spontaneous symmetry breaking because (a) nothing is added to QCD to catalyse this remarkable outcome and (b) no simple change of variables in the QCD action can reveal it. Instead, following quantisation of classical chromodynamics, with its massless gluons and quarks, a large mass-scale emerges in both the gauge-and matter-sectors.
DCSB is empirically revealed very clearly in properties of the pion, whose structure in QCD is described by a Bethe-Salpeter amplitude: where k is the relative momentum between the pion's valence-quark and -antiquark (defined such that the scalar functions in Equation (14) are even under k · P → −k · P) and P is their total momentum. Γ π (k; P) is directly related to an entity that, if a nonrelativistic limit were appropriate, would become the pion's Schrödinger wave function.
In chiral-limit QCD, if, and only if, chiral symmetry is dynamically broken, then [89,91,93] where f 0 π is the pion's leptonic decay constant evaluated in the chiral-limit. This identity is exceptional. It is independent of the renormalisation scheme and true in any covariant gauge, and it entails that the two-body problem is solved, nearly completely, once the solution to the one body problem is known. Furthermore, Equation (15) has many corollaries, e.g., it ensures that chiral-QCD generates a massless pion in the absence of a Higgs mechanism; predicts m 2 π ∝ m on m 0, where m is the current-quark mass; and entails that the chiral-limit leptonic decay constant vanishes for all excited-state pseudoscalar mesons with nonzero isospin [125,126]. It is also the keystone that supports the success of chiral effective field theories in nuclear physics.
Equation (15), which may be described as a quark-level Goldberger-Treiman relation, is the most basic statement in QCD of the Nambu-Goldstone theorem [37,38]. (Additional considerations applying to the η-η complex are described elsewhere [127,128].) The Nambu-Goldstone theorem is fundamentally an expression of equivalence between the one-body problem and the two-body problem in QCD's colour-singlet pseudoscalar channel.
It means that pion properties are practically a direct measure of the dressed-quark mass function rendered in Figure 3. Thus, the qualities of the nearly-massless pion are, enigmatically, the cleanest expression of the mechanism that is responsible for almost all visible mass in the Universe.
Reinstating the Higgs mechanism, so that the light-quarks possess their small current-masses, roughly commensurate with the electron mass, then DCSB is responsible for, inter alia, the physical size of the pion mass (m π ≈ 0.15 m p ); the large mass-splitting between the pion and its valence-quark spin-flip partner, the ρ-meson (m ρ > 5 m π ); and the natural scale of nuclear physics, m p ≈ 1 GeV. Curious things also happen to the kaon. Similar to a pion, except that a strange quark replaces one of the light quarks, the kaon comes to possess a mass m K ≈ 0.5 GeV. In this case, a competition is taking place between dynamical and Higgs-driven mass generation [114,115,[128][129][130][131].
Expanding upon these observations, it is worth highlighting that the physical size of m π is actually much bigger than that linked with the Higgs mechanism for light quarks. Empirically, the scale of the Higgs effect in the light-quark sector is ∼1 MeV [13]. As remarked above, m π = 0 without a Higgs mechanism, but the current-masses of the light quarks in the pion are the same as they are in nucleons. Therefore, the simple Higgs mechanism result is m π ≈ (m u + m d ), yielding a value which is only 5% of the physical mass. This physical mass emerges as the result of a hefty DCSB-induced enhancement factor, which multiplies the current-quark mass contribution [42]. However, the scale of DCSB is ∼ m p /3, i.e., the size of a typical u or d constituent-quark mass, and the special NG-character of the pion means that although it should have a mass ∼ (2/3)m p ≈ m ρ , most of that mass is cancelled by gluon binding effects owing to constraints imposed by DCSB [124]. The mechanism will now be explained.

Pion and the Trace Anomaly
At this point, it is possible to resolve the dichotomy expressed by Equations (9) and (10). These statements hold equally on a measurable neighbourhood of the chiral limit because each hadron's mass is a continuous function of the current-masses of its valence-quarks/antiquarks. So consider that for any meson, H 2 , constituted from a valence-quark with current-mass m q and a valence-antiquark with mass mq, where m Σ = (m q + mq), viz. the scalar form factor at zero momentum transfer measures the reaction of the meson's mass-squared to a variation in current-quark mass. It is merely a misleading convention to define the meson's σ-term as σ H 2 = s H 2 (0)/[2m H 2 ]. Notably, the pion (and any other NG mode) possesses the peculiar property that which is the statement that in the vicinity of the chiral limit, 100% of the pion mass-squared owes to the existence of the current-mass in L QCD . One should compare this result with that for the pion's valence-quark spin-flip partner, i.e., the ρ-meson [132]: indicating that just 6% of the ρ-meson's mass-squared is generated by the current-mass term in L QCD .
The remainder arises largely as a consequence of EHM [42], as suggested by the fact that m ρ ≈ 2M(0), where M(k 2 ) is the dressed-quark mass function in Figure 3. The key to understanding Equation (10) is Equation (15) and three associated Goldberger-Treiman-like relations, which are exact in chiral QCD. Utilising these identities when working with those DSEs needed to describe a pseudoscalar meson, an algebraic proof of the following statement can be constructed [89,91,93]: at every order in a symmetry-preserving analysis, the masses generated for the valence-quarks that constitute the system are exactly cancelled by the attraction produced by interactions between them. This cancellation ensures that the initial two-valence-parton system, which began massless, grows into a complex system, with a bound-state wave function tied to a pole in the scattering matrix at P 2 = 0. Stated simply, Equation (10) is obtained through cancellations between one-body dressing and two-body binding effects: with the sum being exactly zero if, and only if, chiral symmetry is broken dynamically in the Nambu pattern. (A full discussion is provided in Sec. 3.3 of Ref. [124].) Away from the chiral limit for NG modes and always in other channels, the cancellation is incomplete.
It is important to remark that such a transparent resolution of the conundrum expressed by Equations (9) and (10) is impossible if one insists on using a parton model basis, in which the trace anomaly operator is given simply by Equation (7). One must employ a modern quasiparticle formalism. Then, with dressed-particle propagators and bound-state wave functions, obtained at a low renormalisation scale, ζ m p , as solutions of coupled integral equations, each of which sums a countable infinity of diagrams, Equation (10) can be re-expressed: Equations (20) describe the transfigurement of the chiral-limit parton-basis expression of the trace anomaly's expectation value in the pion into a new form, written in terms of a non-perturbatively-dressed quasi-particle basis, with Q denoting dressed-quarks and G the dressed-gluon field strength tensor. Here, the first term expresses the one-body-dressing content of the trace anomaly. It is positive. Patently, a massless valence-quark (antiquark) gaining a large mass by way of interactions with its own gluon field is an effect of the trace-anomaly in what may be described as the one-quasiparticle subspace of a complete pion wave function. The second term expresses the two-particle-irreducible (2PI) interaction content of the forward scattering process represented by this trace anomaly matrix element. It is negative and acquires a scale because the gluon-and quark-propagators and the couplings in the 2PI processes have all gained a mass-scale.
The discussion of Equations (19) and (20) and their connection with Equation (10) bring a well-known adage to mind, quoted here from S. Weinberg in [133]: "You may use any degrees of freedom you like to describe a physical system, but if you use the wrong ones, you'll be sorry!"

Pion Wave Function
Empirical signals of the emergence of mass are ubiquitous, but their appearance typically changes from one system to another. This means that it is essential to study a wide variety of observables because each one is apt to expose different aspects of the underlying mechanisms. It is useful, therefore, to begin with the most obvious, whose importance derives from Equation (15), i.e., the pion's leading-twist two-particle parton distribution amplitude (PDA), ϕ π , the simplest component of its light-front wave function. Any framework that is capable of delivering a hadron's Poincaré-covariant bound-state amplitude can also be used to calculate its PDA, and in this case, ϕ π is given by a light-front projection of the pion's Bethe-Salpeter wave function [134] f π ϕ π (x; ζ) = tr CD Z 2 Λ dq δ(n · q + − x n · P) γ 5 γ · n S(q + )Γ(q; P)S(q − ) .
In Equation (21), the trace is over colour and spinor indices; Λ dq is a Poincaré-invariant regularisation of the four-dimensional integral, with Λ the ultraviolet regularization mass-scale; Z 2 (ζ, Λ) is the quark wavefunction renormalisation constant; n is a light-like four-vector, n 2 = 0; and P is the pion's four-momentum, P 2 = −m 2 π and n · P = −m π in the pion's rest frame. Two distinctively different truncations of QCD's DSEs [134] have been used to calculate the amplitude in Equation (21). Both agree: compared with the asymptotic profile, which is valid on Λ QCD /ζ 0, there is a marked broadening of ϕ π (x; ζ), which owes exclusively to DCSB, i.e., the emergence of mass as exhibited in Figure 3. This causal connection may be asserted because the PDA is calculated at a low renormalisation scale in the chiral limit, in which case the quark mass function owes entirely to DCSB via Equation (15). Moreover, the PDA's dilation is related to the rate at which a dressed-quark approaches the asymptotic bare-parton limit. The prediction determined using the most sophisticated kernel [114,115] is depicted in Figure 4: It can be verified empirically at JLab 12, e.g., in measurements of the pion's electromagnetic form factor (see Section 5.2 below).
It is worth remarking here that the authors of [134] chose to reconstruct the pion's DA from its Mellin moments using an order-α Gegenbauer expansion. This is useful as a first step because the procedure converges rapidly, so enables the qualitative feature of broadening driven by EHM in the SM to be exposed. Such broadening, however, need not and should not disturb the DA's endpoint  (22)-determined using the most sophisticated available symmetry preserving DSE kernels; dot-dashed green curve-original prediction from the work in [134]. These PDA results are consistent with contemporary lQCD results [135,136]. The dashed orange curve is ϕ asy π (x) = 6x(1 − x), the limiting form under QCD evolution [137][138][139]. The PDAs are dimensionless. behaviour, which QCD predicts to be linear in the neighbourhoods x 0, 1. Therefore, as a second step, the authors of [114,115] re-expressed the result from the work in [134] as the function in Equation (22). The first eleven moments agree at the level of 1.6 (1.4)%, i.e., well within any sensible estimate of uncertainty in the computation of high-order moments. Pointwise, the curves are practically indistinguishable, as highlighted in Figure 4.
It is important to emphasise that the computed PDAs in Figure 4 are concave functions. Such pointwise behaviour contrasts markedly with the "bimodal" or "double-humped" distributions that have been favoured in phenomenological applications by some authors [140]. It should be understood in this connection that a double-humped profile for the twist-two PDA places it in the class of distributions produced by a Bethe-Salpeter wave function which is suppressed at zero relative momentum, instead of maximal thereat. No ground-state solution of the pseudoscalar or vector meson Bethe-Salpeter equation exhibits corresponding behaviour [141][142][143][144][145]. Thus, a bimodal distribution cannot be an accurate pointwise rendering of the PDA for a ground-state meson. Notwithstanding that, if one employs such a distribution in a phenomenological application for which only its lowest Mellin moments are important, then some carefully-constrained bimodal distributions may supply useful approximations to the moments of a broad, concave PDA and thereby prove practically useful.
A question of more than thirty-years standing can be answered using Figure 4, namely, when does ϕ asy π (x) provide a good approximation to the pion PDA? Plainly, not at ζ H . The ERBL evolution equation [137][138][139] describes the ζ-evolution of ϕ π (x; ζ), and applied to ϕ π (x; ζ) in Figure 4, one finds [135,136] that ϕ asy π (x) is a poor approximation to the true result even at ζ = 10 GeV. Thus at energy scales accessible to experiment, the twist-two PDAs of ground-state hadrons are "fat and squat". Evidence supporting this picture had long been accumulating [146][147][148][149], and the dilation is confirmed by simulations of lQCD [135,136].

Pion Electromagnetic Form Factor
The cross sections for many hard exclusive hadronic reactions, i.e., processes involving a highly energetic probe that strikes a target and leaves it intact, can be written in terms of the PDAs of the hadrons involved. An example is the elastic electromagnetic form factor of the pion, for which the prediction can be stated thus [137][138][139]: where α(Q 2 ) is the running coupling, which is practically equivalent toα(Q 2 ) in Figure 2 on any domain within which perturbation theory is valid. The value of Q 0 is not predicted by perturbative QCD (pQCD), but it is computable in any non-perturbative framework that veraciously expresses EHM. It was anticipated that JLab would verify the fundamental prediction in Equation (23). Indeed, the first publication by the F π Collaboration [150] heralded the beginning of a new era in mapping the pion's internal structure. Ensuing measurements [151][152][153][154][155] confirmed the data's trend, leading to a common perception that, at Q 2 = 2.45 GeV 2 , one remains far from the resolution field wherein the pion acts like an elementary partonic quark-antiquark pair, i.e., far from validating Equation (23). This conclusion was based on the assumption that inserting ϕ asy π (x) into Equation (23) delivers a valid approximation at ζ 2 = Q 2 = 2.45 GeV 2 , so that The result in Equation (24) is a factor of 2.7 smaller than the empirical value quoted at Q 2 = 2.45 GeV 2 [154]: 0.41 +0.04 −0.03 , and a factor of three smaller than that computed at Q 2 = 4 GeV 2 in [156]. At the time, the authors of [156] supplied the lone prediction for the Q 2 -dependence of F π that was both applicable on the full spacelike region mapped reliably by experiment and confirmed thereby.
Here, the perception of a mismatch and a real discrepancy are not equivalent because, as elucidated above, Q 2 = 4 GeV 2 is not within the domain Λ 2 QCD /Q 2 0 upon which Equation (23) used with ϕ asy π (x) provides a valid approximation. This being so, and given the successful prediction in [156], it was natural to ask whether DSE methods could address the issue of the ultimate validity of Equation (23).
Developments within the past decade have made that possible and now a computation of the pion's electromagnetic form factor is available to arbitrarily large-Q 2 [131,157,158]. The result is the solid (black) curve in Figure 5. In addition, the study enables that result to be correlated with Equation (23) via the modern PDA computed in the same framework, which is the dashed (blue) curve in Figure 5. On the depicted domain, this leading-order, leading-twist QCD prediction, computed with a pion valence-quark PDA evaluated at a scale fitting the experiment, underestimates the full DSE result by a fairly uniform 15%. The mismatch is explained by a mix of higher-order, higher-twist corrections to Equation (23) in pQCD on the one hand and, on the other hand, shortcomings in the leading-order DSE truncation used in [131,157,158], which predicts the right power-law behaviour for the form factor but not precisely the correct anomalous dimension (logarithm exponent) in the strong-coupling calculation. It is now anticipated that the upgraded JLab facility will reveal a maximum in Q 2 F π (Q 2 ) at Q 2 ≈ 6 GeV 2 and an experiment at Q 2 = 9 GeV 2 will see an indication of parton model scaling and scaling violations for the first time in a hadron elastic form factor. While JLab's grip on these things might be tenuous, the reach of an electron ion collider would certainly enable validation of these predictions [39].
. Solid curve (black)-theoretical prediction [131,157,158]; dashed curve (blue)-pQCD prediction computed with the modern, dilated pion PDA described in Section 5.1; and dotted (red) curve-pQCD prediction computed with the asymptotic profile, ϕ asy π (x), which had previously been used to guide expectations for the asymptotic behaviour of Q 2 F π (Q 2 ). The filled-circles and -squares represent existing JLab data [154] and the filled diamonds and triangle, whose normalisations are arbitrary, indicate the projected Q 2 -reach and accuracy of forthcoming experiments [159,160].

Valence-Quark Distributions in the Pion
Given the pions' simple valence-quark content, i.e., one quark and one antiquark, a basic quantity in any discussion of their structure is the associated parton distribution function (PDF), q π (x; ζ). This density charts the probability that a valence q-quark in a pion carries a light-front fraction x of the system's total momentum [161], and one of the earliest predictions of pQCD is [162][163][164]: Moreover, the exponent evolves as ζ increases beyond ζ H , becoming 2 + γ, where γ 0 is an anomalous dimension that increases logarithmically with ζ [165][166][167][168]. (In the limit of exact G-parity symmetry, which is a good approximation in the SM, u π + (x) =d π + (x), etc. Therefore, it is only necessary to discuss one unique distribution.) As described in connection with Equations (19) and (20), for nature's pion there is near-complete cancellation between (a) strong-mass-generating dressing of the valence-quark and -antiquark and (b) binding attraction. Such effects distinguish this system from the more massive kaon, within which the cancellation is far less efficient because the s-quark current mass generated by the Higgs is so much larger than that of the uand d-quarks. Consequently, high-precision measurements of the valence-parton distributions in the pion and kaon are a high priority at existing and anticipated facilities [39,[169][170][171].
Such efforts are driven by ongoing progress in theory. Marking one significant class of advances, lQCD is beginning to yield results for the pointwise behaviour of the pion's valence-quark distribution [172][173][174][175][176], offering promise for information beyond the lowest few moments [177,178].
Extensions of the continuum analysis in [179] are also yielding new insights [180], leading to the first parameter-free predictions of the valence, glue and sea distributions within the pion [181,182], unifying them with, inter alia, electromagnetic pion elastic and transition form factors [128,131,157,158,183,184].
The analysis reveals that, like the pion's PDA in Figure 4, the valence-quark distribution function is hardened by DCSB, producing the following apportioning of light-front momentum at the scale ζ = ζ 2 = 2 GeV: x valence = 0.48(3) , x glue = 0.41(2) , x sea = 0.11 (2) . (26) A similar valence-quark momentum fraction was obtained in Ref. [185] by analysing data on π-nucleus Drell-Yan and leading neutron electroproduction [185]: 2x π q = 0.49(1) at ζ = ζ 2 . Even though this phenomenological PDF yields a compatible momentum fraction, its x-profile is different. In fact, the phenomenological PDF conflicts with the QCD constraint, Equation (25). Significantly, the analysis in [185] ignored threshold resummation effects, which are known to have a material impact at large x [186,187]. (Similar remarks apply to the analysis in [188].) Importantly, as illustrated in Figure 6, after evolution [165][166][167][168] to ζ = ζ 5 = 5.2 GeV, usingα(k 2 ) in Figure 2 to integrate the evolution equations, the continuum prediction for u π (x) from the works in [181,182] matches that obtained using lQCD [175]. Given that no parameters were varied in order to procure this or any other outcome in [181,182], a remarkable, modern confluence has been reached, which suggests that real strides are being made toward understanding pion structure and its relation to the emergence of mass. (Predictions for kaon PDFs are described in [114,115].) Figure 6. Pion valence-quark momentum distribution function, xu π (x; ζ 5 ): dot-dot-dashed (grey) curve within shaded band-lQCD result [175]; long-dashed (black) curve-early continuum analysis [179]; and solid (blue) curve within shaded band-modern, continuum calculation [181,182]. From the works in [181,182]: gluon momentum distribution in pion, xg π (x; ζ 5 )-dashed (green) curve within shaded band; and sea-quark momentum distribution, xS π (x; ζ 5 )-dot-dashed (red) curve within shaded band. (The shaded bands indicate the size of calculation-specific uncertainties, detailed in the source material [175,181,182].) Data (purple) from in [189], rescaled according to the analysis in [186].

Emergence of Diquark Correlations
The emergence of mass is also expressed in the properties of baryons, something which can readily be seen when the three valence-quark bound-state problem is studied with the same level of sophistication that is now typical for mesons. In this connection, DCSB, as displayed in the momentum-dependence of the dressed-quark mass- Figure 3, is as important to baryons as it is to mesons. In fact, one important consequence of DCSB is that any interaction able to create composite pseudo-NG modes from a light dressed-quark and -antiquark, and reproduce the observed value of their leptonic decay constants, will also generate tight colour-antitriplet correlations between any two of a nucleon's dressed-quarks. This conclusion is based upon evidence gathered in thirty years of studying two-and three-body bound-states in hadron physics [190].
Realistic diquark correlations are also soft and interacting. They all carry charge, scatter electrons, and have an electromagnetic size which is similar to that of the kindred mesonic system, see, e.g., in [194,204,205]: with r {uu} 1 + > r [ud] 0 + . As in the meson sector, these scales are set by that associated with DCSB. It should be emphasised that such fully dynamical diquark correlations are vastly different from the static, pointlike "diquarks" introduced originally [206] in an attempt to solve the so-called "missing resonance" problem [207], viz. the fact that quark models predict a far greater number of baryon states than were observed in the previous millennium [208]. Moreover, their existence enforces certain distinct interaction patterns for the singly-and doubly-represented valence-quarks within the proton, as exhibited elsewhere [81,84,[199][200][201][202][209][210][211][212][213].
The existence of these strong correlations between two dressed quarks is the path to converting the three valence-quark bound-state problem into the simpler Faddeev equation drawn in Figure 7, without loss of dynamical information [214]. The three-gluon vertex, a definitive feature of QCD's non-Abelian nature, is not an explicit part of the scattering kernel in this picture. Instead, one profits from the following features, phase-space factors materially enhance two-body interactions over n ≥ 3-body interactions and the primary three-body force, produced by a three-gluon vertex attaching once, and only once, to each valence quark, vanishes when projected into the colour-singlet channel; and subsequently capitalises on diquark dominance in the two-body subsystems. Then, while iterated, overlapping three-body terms might alter fine details of baryon structure, the primary effect of non-Abelian multi-gluon vertices is manifested in the formation of diquark correlations. (On the other hand, the three-gluon vertex appears to play a material and distinguishable role in the formation of hybrid hadrons [215] and glueballs [216]). Accordingly, the active kernel describes binding in the baryon through diquark breakup and reconstitution, mediated by exchange of a dressed-quark. The properties and interactions of such a baryon are chiefly determined by the quark+diquark structure evident in Figure 7. The spectrum of baryons produced by the Faddeev equation [197,198,211,217] is like that found in the three-constituent quark model and consistent with lQCD analyses [218]. Notably, modern data and recent analyses have reduced the number of missing resonances [219][220][221][222][223][224].

Proton Wave Function
After solving the Faddeev equation in Figure 7 for the proton's Faddeev amplitude, one can then compute the proton's dressed-quark leading-twist PDA [203]. The result for this proton "wave function" is depicted in Figure 8. Here, x 1 , x 2 , x 3 indicate the fraction of the bound-state's light-front momentum carried by the associated quark; naturally, x 1 + x 2 + x 3 = 1. The amplitudes are dimensionless; hence, the height is simply a real number. Table 1 lists the proton PDA's four lowest-order moments. They reveal important insights, e.g., when the proton is pictured as solely a quark+scalar-diquark correlation, x 2 u = x 3 d , because these two form the scalar quark+quark correlation, and the system is very asymmetric, with the PDA's peak being shifted strongly in favour of x 1 u > x 2 u . This outcome conflicts with lQCD results [225,226]. On the other hand, as explained above, realistic Faddeev equation calculations suggest that pseudovector diquark correlations are an essential piece of the proton's wave function. When these {uu} and {ud} diquarks are included, momentum is shared more evenly, shifting from the spectator u(x 1 ) quark into u(x 2 ), d(x 3 ). Adding these diquarks with the known weighting, the PDA's peak shifts back toward the centre, relocating to ([x]) = (0.55, 0.23, 0.22), and the computed values of the first moments line up with those computed using lQCD. This convergence delivers a more complete understanding of the lQCD simulations, which by that means are seen to confirm a picture of the proton as a bound-state with both tight scalar and pseudovector diquark correlations, wherein the scalar diquarks are responsible for ≈ 65% of the Faddeev amplitude's canonical normalisation.
Importantly, like ground-state S-wave mesons [134,136,[227][228][229][230], the leading-twist PDA of the ground-state nucleon is both broader than ϕ asy N ([x]) (defined by Table 1, Row 1) and decreases monotonically away from its maximum in all directions, i.e., the PDAs of these ground-state S-wave systems possess endpoint enhancements, but neither humps nor bumps. Models which produce such humped structures were previously considered reasonable [140]. However, it is now clear that pointwise behaviour of this type is in conflict with QCD. The models may nonetheless be viewed as possessing a qualitative truth, insofar as they represent a way by which endpoint enhancements can be expressed in hadron PDAs if one limits oneself to the polynomial basis characterising QCD in the neighbourhood Λ QCD /ζ 0.
As with the pion prediction described in Section 5.2, the veracity of this result for the proton PDA can be tested in future experiments. For instance, it can be used to provide the first realistic evaluation of the scale at which exclusive measurements involving proton targets may reasonably be compared with predictions based on pQCD hard scattering formulae. Analogous to the pion case, the value of this mass-scale is an empirical manifestation of the emergence of mass, here within the three-valence-quark proton bound-state. Table 1. Computed values of the first four moments of the proton and Roper-resonance PDAs. The error on f N , a dynamically-determined quantity which measures the proton's "wave function at the origin", reflects a nucleon scalar diquark content of 65 ± 5%, and values in rows marked with " ⊃ av" were obtained assuming the baryon is constituted solely from a scalar diquark. Including axial-vector diquark correlations, the mean absolute relative difference between continuum and lattice results improves by 36%. (All results listed at a renormalisation scale ζ = 2 GeV.)

Proton's First Radial Excitation
Almost immediately after discovery of the Roper resonance [231][232][233][234][235], the proton's first positive-parity excitation, questions were asked regarding the nature of like-parity excitations of ground-state positive-parity baryons. A lucid picture is now emerging following [84], (i) the acquisition and analysis of a large amount of very precise nucleon-resonance electroproduction data, involving single-and double-pion final states on a large domain of energy and momentum-transfer; (ii) development of an advanced dynamical reaction theory, able to simultaneously describe all partial waves extracted from available data; (iii) formulation and extensive application of the DSE approach to the continuum bound state problem in relativistic quantum field theory, which expresses diverse local and global impacts of EHM; and and (iv) the improvement of constituent quark models so that they can also qualitatively incorporate these features of strong QCD. In this picture, such states are primarily radial excitations of the associated ground-state baryon, comprised of a well-defined dressed-quark core supplemented by a meson cloud.
Using the Faddeev equation framework sketched above, properties of the Roper-resonance's dressed-quark core have been exposed [195,199]: it is found that the scalar functions in the Roper's Faddeev amplitude have a zero at quark+diquark relative momentum √ 2 ≈ 0.4 GeV≈ 1/[0.5 fm].
Working with this input, the authors of [203] delivered the associated leading-twist PDA, depicted in the rightmost panel of Figure 8 and whose first four moments are listed in Table 1. The prediction reveals some interesting features, e.g., the PDA of this excitation is not positive definite and there is a conspicuous locus of zeros in the lower-right section of the barycentric plot, both of which duplicate aspects of the wave function describing the first radial excitation of a quantum mechanical system that are also seen in the PDAs of meson radial excitations [236,237], and here the influence of pseudovector diquarks is contrary to that in the ground-state, viz. they shift momentum into u(x 1 ) from u(x 2 ), d(x 3 ). These observations highlight that the ground state is just one isolated member of a set of Hamiltonian eigenvectors with infinitely many elements: many Hamiltonians can possess practically equivalent ground states and yet lead to excited-state spectra that are vastly different. Moreover, masses alone, as infrared-dominated quantities, contain relatively little information. Different Hamiltonians may adequately reproduce known hadron spectra; but these same Hamiltonians can yield predictions that disagree markedly when used to compute structural properties. Such properties-like wave functions and the Q 2 -dependence of elastic and transition form factors-possess the greatest discriminating power. Therefore, study of the structure of nucleon resonances is a critical complement to that of ground-state nucleons and mesons because it is capable of revealing additional novel features of strong QCD. Modern theory must be deployed to compute observable properties of all these systems. Aspects of this effort are sketched elsewhere [41,43]. Paradigmatic shifts in our understanding of these things are being driven by new, precise form factor data. This is nowhere more evident than in analyses of experimental data acquired in the past vicennium, which have established a new ideal. Namely, despite its elementary valence-quark content, the nucleon's internal structure is very complex. For instance, there are measurable differences between the distributions of charge and magnetisation throughout the interiors [238] and between the way these qualities are carried by the different quark flavours [202,209,239]. The aim now is to explain the observations in terms of emergent features of the strong interaction.

Emergent Features of Nucleon Form Factors
In this connection, it is here worth highlighting that the behaviour of the proton's electric form factor on Q 2 4 GeV 2 is particularly sensitive to the running of the dressed-quark mass (displayed in Figure 3), especially the rate at which the dressed-quark mass runs between the infrared and ultraviolet domains [209,213,240].
The proton's momentum-space charge and magnetisation distributions are expressed in combinations of the two Poincaré-invariant elastic form factors that are needed to express the proton's electromagnetic current: where m N is the nucleon mass; Q = p − p, u(p) andū(p ) are, respectively, spinors describing the incident and scattered proton; and F 1,2 (Q 2 ) are the proton's Dirac and Pauli form factors. The combinations that feature in the electron-proton elastic scattering cross section are the charge and magnetisation distributions [241] Data available before the year 1999 led to a view that µ p G p M (Q 2 = 0) = 1; therefore, a conclusion that the distributions of charge and magnetisation within the proton are approximately identical [242]. At the time, the proton was viewed as a simple bound state, wherein quark+quark correlations and attendant orbital angular momentum played little role. Equation (31) is consistent with this picture.
The situation changed dramatically when the combination of high-energy, -current and -polarisation at JLab enabled a new type of experiment to be performed, viz. polarisation-transfer reactions [238], which are directly proportional to G E (Q 2 )/G M (Q 2 ) [243,244]. A series of these experiments [238,[245][246][247][248] has determined that where the constant is such that the ratio might become negative for Q 2 8 GeV 2 . This behaviour contrasts starkly with Equation (31); and since the proton's magnetic form factor is reliably known on a spacelike domain that extends to Q 2 ≈ 30 GeV 2 [249,250], the Q 2 -dependence of G p E /G p M exposes novel features of the proton's charge distribution, as expressed in G p E (Q 2 ). An understanding of the behaviour in Equation (32) is provided by the analyses in [209,213,240], and the answer lies largely with the proton's Pauli form factor. F p 2 is a gauge of the proton's magnetisation distribution. Ultimately, this magnetisation is borne by the dressed-quarks and influenced by the k 2 -dependence of (i) their propagators and (ii) the correlations amongst them. Both are expressed in the Faddeev wavefunction.
Thus, for the sake of argument, suppose that dressed-quarks are described by a momentumindependent mass-function, e.g., as obtained using a symmetry-preserving regularisation of a vector × vector contact interaction [251]. They then behave as Dirac particles with constant Dirac values for their magnetic moments. Consequently, the composite proton possesses a hard Pauli form factor and this produces a zero in G p E at Q 2 ≈ 4 GeV 2 [240,252]. Alternatively, suppose that the dressed-quarks have a momentum-dependent mass-function, like that depicted in Figure 3, which is large at small momenta but vanishes as their momentum increases. At infrared momenta the dressed-quarks will then act as constituent-like particles with a large magnetic moment, but as the probe momentum grows, their mass and magnetic moment will drop toward zero. (N.B. Massless fermions do not have a measurable magnetic moment [253], so that any significant magnetic moment for a constituent-like quark is an emergent feature.) Such dressed-quarks will produce a proton Pauli form factor that is large on Q 2 0 but drops rapidly on the domain of transition between the infrared and ultraviolet domains, to give a very small result at large-Q 2 . The proton's Dirac form factor is far less sensitive to spin-related effects; hence the interplay between the Dirac and Pauli form factors, expressed in Equation (30), entails that G E (Q 2 ) must have a zero at larger values of Q 2 when calculated with a running mass function than we computed with momentum-independent dressed-quark masses. The precise form of the Q 2 -dependence will depend on the evolving nature of the angular momentum correlations between the dressed-quarks [209,213].
The Class-C DSE analyses in [209,213] implement a dressed-quark mass function that is distinguished by a particular transition rate between the non-perturbative and perturbative domains. If that rate were increased, then the transformation to partonic quarks would become more rapid; hence the proton's Pauli form factor would drop even more quickly to zero. In this event, the quark angular momentum correlations, embodied in the diquark structure, remain but the separate dressed-quark magnetic moments diminish markedly. Thus a quicker transition will push the zero in µ p G p E /G p M to larger values of Q 2 . Moreover, there will be a rate of transformation beyond which the zero disappears completely: there is no zero in a theory whose mass function rapidly becomes partonic. (For instance, pQCD analyses cannot produce a zero.) These expectations are realised in explicit calculations, as illustrated in Figure 4 of Ref. [209].
It follows that the possible existence and location of a zero in the ratio of proton elastic form factors are a fairly direct measure of the character of EHM in the SM. Consequently, in pushing experimental measurements of this ratio, and thereby the proton's electric form factor, to larger momenta, i.e., in reliably determining the proton's charge distribution, there is an exceptional opportunity for a positive dialogue between experiment and theory. That feedback should assist greatly with contemporary efforts to reveal the character of the SM's strong interaction and its emergent phenomena.
Given such potential, the analysis in [209] was recently revisited [213], using improved algorithms to calculate the form factors on Q 2 10 m 2 p . The new study employed a statistical implementation of the Schlessinger point method (SPM) [31,210,[257][258][259][260] for the interpolation and extrapolation of smooth functions to deliver predictions for form factors on Q 2 > 9 m 2 p with a quantified uncertainty estimate. Results for the form factor ratios discussed herein are presented in Figure 9.
For the proton, the new analysis predicts This value is compatible with, although a little larger than, that obtained earlier [209]: Q 2 z p ≈ 9.5 GeV 2 ; and a more recent inference based on ρ-meson elastic form factors [261]: Q 2 z p ≈ 9.4(3) GeV 2 . Regarding the neutron, the authors of [209] predicted a peak in this ratio at Q 2 ≈ 6 GeV 2 , which is reproduced in [213]. Furthermore, it located a zero at Q 2 z n ≈ 12 GeV 2 . With the statistical SPM method for reaching to large-Q 2 : viz. at 1σ SPM confidence level, this ratio is likely to exhibit a zero, but it probably lies beyond the reach of 12 GeV beams at JLab. On the other hand, the prediction of a peak in R n EM = µ n G n E /G n M , which is a harbinger of the zero in this ratio, can be tested at the 12 GeV JLab.
The properties of the dressed-quark propagators and bound-state amplitudes which influence the appearance of a zero in R n EM are qualitatively the same as those described in connection with R p EM . However, because of the different electric charge weightings attached to the quark contributions in the neutron (1 valence u-quark and 2 valence d-quarks), the quantitative effect is opposite to that for the proton (2 valence u-quarks and 1 valence d-quark). Namely, when the transformation from dressed-quark to parton is accelerated, the zero occurs at smaller Q 2 . On the other hand, a model which generates a momentum-independent dressed-quark mass typically produces no zero in the neutron ratio [240].   [246]; black down-triangles [247]; and cyan diamonds [248]); right, compared with available lQCD results, drawn from Ref. [254]. Lower panels-Neutron. Left, comparison with data (blue circles [255] and green squares [256]); right, with available lQCD results, drawn from the work in [254]. In all panels, the 1σ band for the SPM approximants is shaded in light blue.
In order to understand the source of these features, consider that the strange quark contributes very little to nucleon electromagnetic form factors [262,263] must occur at a larger value of Q 2 than the zero in G p E , a prediction confirmed in Figure 9. , G E p Figure 10. Data: neutron electric form factor-blue circles [256], and proton electric form factor-green squares [264]. Associated curves depict least-squares fits obtained via a jackknife analysis of the depicted G n E , G p E data. Vertical black line within grey bands marks the boundary of the domain described by Equation (38).
A curious effect follows from the combination of faster-than-dipole decrease of the proton's electric form factor (and probable appearance of a zero) and steady increase of R n EM to a peak at Q 2 ≈ 6 GeV 2 . Namely, there must be a domain of larger Q 2 upon which the neutron's electric form factor is greater than that of the proton's; therefore, at some value of Q 2 , the electric form factor of the neutral compound fermion is actually larger than that of its positive counterpart. This was predicted in [209] and confirmed by the updated analyses in [213]. It is therefore worth testing with available data [256,264]. Figure 10 depicts G n E , G p E data along with a set of associated least-squares fits produced using a jackknife analysis. On the domain currently covered by data, G n E < G p E ; but extrapolation of the fits suggests that With 12 GeV operations under way at JLab, neutron electric form factor data to Q 2 = 10.2 GeV 2 will soon be available [265], enabling this prediction to be tested. It is worth reiterating that the existence and location of a zero in G p E -similarly for G n Eare sensitive to the rate at which the dressed-quark mass function transits from the strong to perturbative QCD domains: denote this rate by R sp M . Note, too, that in contrast to G p,u E , G p,d E evolves more slowly with changes in R sp M . This inertia derives from the d-quark being preferentially contained within a (soft) scalar diquark. Subject to these insights, consider Equation (37): with the location of a zero in G p,d E shifting slowly to larger values of Q 2 , but that in G p,u E moving rapidly (as noted above, the zero in G p,u E disappears if R sp M is sufficiently quick) one is subtracting from G p,d E (Q 2 ) a function whose domain of positive support is increasing in size. That operation typically shifts the zero in G n E to smaller values of Q 2 , potentially enabling a zero in G n E even when that in G p E has disappeared.

Epilogue
The Lagrangian that defines QCD, the strong interaction sector of the Standard Model (SM), appears very simple, yet it is responsible for a remarkable array of high-level phenomena with enormous apparent complexity. One can argue that the foundation for all such effects is laid by the emergence of hadronic mass (EHM), whose consequences are variously expressed in a wide range of empirical observables. This emergence itself is currently beyond a reductionist explanation, unless one is content to find that in the need for ultraviolet renormalisation of four-dimensional quantum field theory. Even supposing the latter, then the size of the associated mass-scale, Λ QCD , is not something which can be predicted from within the SM.
The simplest expression of EHM is to be found in the generation of a running mass for the gluon. Mechanically, this is driven by gluon self-interactions, and the scale of the effect is known once Λ QCD is fixed. It is interesting theoretically to explore the impact on observables induced by reducing/increasing Λ QCD . Perhaps there is a critical value (or set of values) at which the Universe would appear very different?
Gluon mass generation in the SM entails that quarks, which are massless in the absence of a Higgs mechanism, also acquire a running mass, whose scale at infrared momenta is roughly one-third of the proton mass. This effect is known as dynamical chiral symmetry breaking (DCSB). One of its corollaries is the emergence of pseudoscalar Nambu-Goldstone bosons-most recognisably, the pions; and their appearance and properties have a profound impact on the character of the physical Universe. Nuclear physics, e.g., the formation of the elements and the number of elements that may be formed, seems to be very sensitive to the pion mass; and the pion mass is a quantity which is uniquely sensitive to the interplay between emergent and explicit (Higgs) mass generating mechanisms.
A non-perturbative framework that can unify the emergence of gluon and quark masses, and express their impacts on solutions of the associated bound-state equations in quantum field theory and also in the matrix elements which describe complex observable processes, can provide reductive explanations for physical phenomena. No single such framework exists, but an amalgam of high-level approaches is bearing fruit.
Funding: This work was supported in part by Jiangsu Province Hundred Talents Plan for Professionals.
Acknowledgments: These notes are based on results obtained and insights developed through collaborations with many people, to all of whom I am greatly indebted.

Conflicts of Interest:
The author declares no conflicts of interest.

Abbreviations
The following abbreviations are used in this manuscript.