Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Interference with Non-Interacting Free Particles and a Special Type of Detector
Next Article in Special Issue
Masses of Compact (Neutron) Stars with Distinguished Cores
Previous Article in Journal
Constraints on Nuclear Symmetry Energy Parameters
Previous Article in Special Issue
Experimental Determination of the QCD Effective Charge αg1(Q)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Emergence of Hadron Mass and Structure

by
Minghui Ding
1,*,†,
Craig D. Roberts
2,3,*,† and
Sebastian M. Schmidt
1,4,*,†
1
Helmholtz-Zentrum Dresden-Rossendorf, Bautzner Landstraße 400, D-01328 Dresden, Germany
2
School of Physics, Nanjing University, Nanjing 210093, China
3
Institute for Nonperturbative Physics, Nanjing University, Nanjing 210093, China
4
III. Physikalisches Institut B, RWTH Aachen University, D-52074 Aachen, Germany
*
Authors to whom correspondence should be addressed.
All authors contributed equally to this work.
Particles 2023, 6(1), 57-120; https://doi.org/10.3390/particles6010004
Submission received: 21 November 2022 / Revised: 18 December 2022 / Accepted: 19 December 2022 / Published: 11 January 2023

Abstract

:
Visible matter is characterised by a single mass scale; namely, the proton mass. The proton’s existence and structure are supposed to be described by quantum chromodynamics (QCD); yet, absent Higgs boson couplings, chromodynamics is scale-invariant. Thus, if the Standard Model is truly a part of the theory of Nature, then the proton mass is an emergent feature of QCD; and emergent hadron mass (EHM) must provide the basic link between theory and observation. Nonperturbative tools are necessary if such connections are to be made; and in this context, we sketch recent progress in the application of continuum Schwinger function methods to an array of related problems in hadron and particle physics. Special emphasis is given to the three pillars of EHM—namely, the running gluon mass, process-independent effective charge, and running quark mass; their role in stabilising QCD; and their measurable expressions in a diverse array of observables.

Contents
 
1      Introduction2
 
2      Hadron Mass Budgets4
 
3      Gluons and the Emergence of Mass6
 
4      Process-Independent Effective Charge8
 
5      Confinement11
 
6      Spectroscopy16
 
7      Baryon Wave Functions22
 
8      Meson Form Factors26
 
9      Baryon Form Factors29
 
10      Transition Form Factors of Heavy+Light Mesons36
 
11      Distribution Functions42
 
12      Conclusions50
 
References51

1. Introduction

Our Universe exists; and even the small part that we occupy contains much which might be considered miraculous. Nevertheless, science typically assumes that the Universe’s evolution can be explained by some collection of equations—even a single equation, perhaps, which replaces distinct theories of many things with a single theory of everything. Choosing not to approach that frontier, then, within the current paradigm, the Standard Model of particle physics (SM) is given a central role; and it must account for a huge array of observable phenomena. Herein, we focus on one especially important aspect, viz. the fact that the mass of the vast bulk of visible material in the Universe is explained as soon as one understands why the proton is absolutely stable and how it comes to possess a mass m p 1 GeV. In elucidating this connection, we will argue that the theory of strong interactions may deliver far more than was originally asked of it.
We have evidently supposed that quantum gauge field theory is the correct paradigm for understanding Nature. In this connection, it is important to note that, in our tangible Universe, time and space give us four noncompact dimensions. Consider, therefore, that quantum gauge field theories in D 4 dimensions are characterised by an explicit, intrinsic mass scale: the basic couplings generated by minimal substitution are mass-dimensioned and set the scale for all calculated quantities. For D > 4 , such theories manifest uncontrollable ultraviolet divergences, making them of little physical use. In contrast, for D < 4 , they are super-convergent, but are afflicted with a hierarchy problem, viz. dynamical mass-generation effects are typically very small when compared with the theory’s explicit scale [1,2,3,4,5]. Hence, perhaps unsurprisingly, D = 4 is a critical point. Removing Higgs boson couplings, the classical gauge theory elements of the SM are scale-invariant. Taking the step to quantum theories, they are all (at least perturbatively) renormalisable; and that procedure introduces a mass scale. As we have noted, the scale for visible matter is m Nature m p 1 GeV. However, the size of this scale is not determined by the theory; so, whence does it come? Further, how much tolerance does Nature give us? Is the Universe habitable when m Nature ( 1 ± δ ) m Nature , with δ = 0.1 or 0.2 , etc.? It is comforting to imagine that our (ultimate?) theory of Nature will answer these questions, but the existence of such a theory is not certain.
Returning to concrete issues, strong interactions within the SM are described by quantum chromodynamics (QCD). Therefore, consider the classical Lagrangian density that serves as the starting point on the road to QCD:
L QCD = f = u , d , s , q ¯ f [ γ · + i g 1 2 λ a γ · A a + m f ] q f + 1 4 G μ ν a G μ ν a ,
G μ ν a = μ A ν a + ν A μ a g f a b c A μ b A ν c ,
where { q f | f = u , d , s , c , b , t } are fields associated with the six known flavours of quarks; { m f } are their current–masses, generated by the Higgs boson; { A μ a | a = 1 , , 8 } represent the gluon fields, whose matrix structure is encoded in { 1 2 λ a } , the generators of SU ( 3 ) in the fundamental representation; and g is the unique QCD coupling, using which one conventionally defines α = g 2 / [ 4 π ] . As remarked above, if one removes Higgs boson couplings into QCD, so that { m f 0 } in Equation (1), then the classical action associated with this Lagrangian is scale-invariant. A scale-invariant theory cannot produce compact bound states; indeed, scale-invariant theories do not support dynamics, only kinematics [6]. So, if Equation (1) is really capable of explaining, amongst other things, the proton’s mass, size, and stability, then remarkable features must emerge via the process of defining quantum chromodynamics.
This point is placed in stark relief when one appreciates that the gluon and quark fields used to express the one-line Lagrangian of QCD are not the degrees-of-freedom measured in detectors. This is an empirical manifestation of confinement. Amongst other things, a solution of QCD will reveal the meaning of confinement, predict the observable states, and explain how they are built from the Lagrangian’s gluon and quark partons. However, the search for a solution presumes that QCD is actually a theory. Effective theories are tools for use in obtaining a realistic description of phenomena perceived at a given scale. A true theory must be rigorously defined at all scales and unify phenomena perceived at vastly different energies. If QCD really is a well-defined quantum field theory, then it may serve as a paradigm for physics far beyond the SM.
Having raised this possibility, then it is appropriate to provide a working definition of “well-defined” in relation to quantum field theory. Aspects of the mathematical problem are discussed elsewhere [7,8]. Herein, we consider that a quantum (gauge) field theory is well-defined if its ultraviolet renormalisation can be accomplished with a finite number of renormalisation constants, { Z j | j = 1 , , N } , N 10 ,1 all of which can (a) be computed nonperturbatively and (b) remain bounded real numbers as any regularisation scale is removed. Further, that the renormalisation of ultraviolet divergences is sufficient to ensure that any/all infrared divergences are eliminated, i.e., the theory is infrared-complete.
Quantum electrodynamics (QED) is not well-defined owing to the existence of a Landau pole in the far ultraviolet (see, e.g., Reference [9] (Ch. 13) and References [10,11,12,13]). Furthermore, weak interactions are essentially perturbative because the inclusion of the Higgs scalar boson introduces an enormous infrared scale that suppresses all nonperturbative effects; moreover, the Higgs boson mass is quadratically divergent, making the theory non-renormalisable.
On the other hand, as we will explain herein, it is beginning to seem increasingly likely that QCD satisfies the tests listed above; hence, is the first well-defined quantum field theory that humanity has developed. QCD may thus stand alone as an internally consistent theory, so that after quantisation of Equation (1), with nothing further added, it is a genuinely predictive mathematical framework for the explanation of natural phenomena.
We have used a Euclidean metric and consistent Dirac matrices in writing Equation (1) because if there is any hope of arriving at a rigorous definition of QCD, then it is by formulating the theory in Euclidean space. There are many reasons for adopting this perspective. Amongst the most significant being the fact that a lattice regularisation of the theory is only possible in Euclidean space, where one can use the action associated with Equation (1) to define a probability measure [14] (Section 2.1). Notably, a choice must be made because any “Wick rotation” between Minkowski space and Euclidean space is a purely formal exercise, whose validity is only guaranteed for perturbative calculations [15,16]. If QCD really does (somehow) explain the emergence of hadron mass and structure, then nonlinear, nonperturbative dynamics must be crucial. Consequently, one cannot assume that any of the requirements necessary to mathematically justify a Wick rotation are satisfied when calculating and summing the necessarily infinite collection of processes associated with a given experimental observable.
One concrete example may serve to illustrate the point. Both continuum and lattice analyses of the gluon two-point Schwinger function (often called the Euclidean space gluon propagator) yield a result whose analytic properties are very different from those one would obtain in perturbation theory at any finite order [17]. As a consequence, the Minkowski space gluon gap equation that is obtained from the Euclidean form via the standard transcriptions used to implement the Wick rotation [16] (Section 2.3), whilst being similar in appearance, cannot possess the same solutions. Thus, to avoid confusion, one should begin with all such equations formulated in Euclidean space, where the solutions determined have a direct and unambiguous connection with results obtained using numerical simulations of the lattice-regularised theory. Anything else is an unnecessary and potentially misleading pretence. Furthermore, only those Schwinger functions corresponding to observable quantities need have a continuation to Minkowski space, and that can be accomplished following standard notions from constructive field theory ([15] (Sections 3 and 4), [16] (Section 2.3)).
We proceed then by supposing that QCD is defined by the Euclidean-space generating functional built using the Lagrangian density in Equation (1). Here, a new choice presents itself. One might attempt to solve the thus-quantised theory using a lattice regularisation [18,19]. Lattice-regularised QCD (lQCD) is a popular framework, which, owing to growth in computer power and algorithm improvements, is becoming more effective—see, e.g., Reference [20]. On the other hand, continuum Schwinger function methods (CSMs) are also available [14,16,21,22,23,24,25]. Much has been achieved using this approach, especially during the past decade [26,27,28,29,30,31] and particularly in connection with elucidating the origins and wide-ranging expressions of emergent hadron mass (EHM) [32,33,34,35,36,37]. It is upon those advances that we focus herein.

2. Hadron Mass Budgets

There is one generally recognised mass-generating mechanism in the SM; namely, that associated with Higgs boson couplings [38,39]. Insofar as QCD is concerned, there are six distinct such couplings, each of which generates the current–mass of a different quark flavour. Those current–quark masses exhibit a remarkable hierarchy of scales, ranging from an electron-like size for the u and d quarks up to a value five-orders-of-magnitude larger for the t quark ([40] (p. 32)). Faced with such discordance, we choose to begin our discussion of mass by considering the proton and its closest relatives, viz. the π - and ρ -mesons.
The proton is defined as the lightest state constituted from the valence quark combination u + u + d . π + is a pseudoscalar meson built from u + d ¯ valence quarks, and the ρ + is its kindred vector meson partner: in quark models, the π and ρ are identified as 1 S 0 and 3 S 1 states, respectively ([40] (Section 63)). Table 1 presents a breakdown of the masses of these states into three contributions: the simplest to count is that associated with the Higgs-generated current–masses of the valence quarks (HB); the least well understood is that part which has no connection with the Higgs boson (EHM); and the remainder is that arising from constructive interference between these two sources of mass (EHM+HB).
The information listed in Rows 1, 2, 5, and 6 of Table 1 is represented pictorially in Figure 1: plainly, there are significant differences between the upper and lower panels. Regarding the proton and ρ -meson, the HB-alone component of their masses is just 1% in each case. Notwithstanding that, their masses are large and remain so even in the absence of Higgs boson couplings into QCD, i.e., in the chiral limit. This overwhelmingly dominant component is a manifestation of EHM in the SM. It produces roughly 95% of the measured mass. Evidently, baryons and vector mesons are similar in these respects.
Conversely and yet still owing to EHM via its dynamical chiral symmetry breaking (DCSB) corollary, the pion is massless in the chiral limit—it is the SM’s Nambu–Goldstone (NG) mode [27,44,45,46,47,48,49,50,51,52]. Returning to the quark model picture, the only difference between ρ - and π -mesons is a spin-flip: in the ρ , the constituent quark spins are aligned, whereas they are antialigned in π . Yet, their mass budgets are fundamentally different: Figure 1B; cf. Figure 1C. An inability to explain this difference is a conspicuous failure of quark models: whilst it is easy to obtain a satisfactory mass for the ρ , a low-mass pion can only be obtained by fine-tuning the quark model’s potential. Nature, however, does not fine-tune the pion: in the absence of Higgs boson couplings, it is massless irrespective of the size of m ρ and, in fact, the mass of any other hadron.
The kaon mass budget is also drawn—see Figure 1D. In the chiral limit, then, like the π , the K-meson is an NG boson. However, with realistic values of Higgs boson couplings into QCD, the s quark current–mass is approximately 27-times the average of the u and d current–masses [40]: 2 m s 27 ( m u + m d ) . Consequently, the HB wedge in Figure 1D accounts for 20% of m K . The remaining 80% is generated by constructive EHM+HB interference. It follows that comparisons between π and K properties present good opportunities for studying Higgs boson modulation of EHM, because the HB mass fraction is four-times larger in kaons than in pions. Moreover, the array of images in Figure 1 highlights that additional, complementary information can be obtained from comparisons between baryons/vector mesons and the set of kindred pseudoscalar mesons. For instance, studies of spectra (Section 6), transitions between vector mesons and pseudoscalar mesons (Section 10), and comparative analyses of proton and pion parton distribution functions (DFs) (see Section 11). In all cases, predominantly EHM systems, on the one hand, are contrasted/overlapped with final states that possess varying degrees of EHM+HB interference.
These observations highlight that EHM—whatever it is—can be accessed via experiment. The task for theory is to identify and explain its source, then elucidate a broad range of observable consequences so that the origins and explanations can be validated.

3. Gluons and the Emergence of Mass

The requirement of gauge invariance ensures that the Higgs boson does not couple to gluons and precludes any other means of generating an explicit mass term for the gluon fields in Equation (1). Consequently, it is widely believed that gluons are massless; and this is recorded by the Particle Data Group (PDG) [40] (p. 25). (We stress that gluon partons are massless).
In QCD, this “gauge invariance” statement is properly translated into a property of the two-point gluon Schwinger function. Namely, using the class of covariant gauges as an illustrative tool, characterised by a gauge-fixing parameter ξ , the inverse of the gluon two-point function can be expressed in terms of a gluon vacuum polarisation (or self-energy):
D μ ν 1 ( k ) = δ μ ν k 2 k μ k ν ( 1 ξ ) + Π μ ν ( k ) = : 0 D μ ν 1 ( k ) + Π μ ν ( k ) + ξ k μ k ν ,
where k is the gluon momentum. (Regarding ξ , common choices in perturbation theory are ξ = 0 , 1 , viz. Landau and Feynman gauges, respectively.) Gauge invariance (BRST symmetry of the quantised theory [53] (Ch. II)) is expressed in the following Slavnov–Taylor identity [54,55]:
k μ Π μ ν ( k ) = 0 = Π μ ν ( k ) k ν .
This restrictive, yet generous, constraint states that interactions cannot affect the four-longitudinal component of the gluon two-point function, but leaves room for modifications of the propagation characteristics of the three four-transverse degrees-of-freedom.
Equation (3) means
Π μ ν ( k ) = [ δ μ ν k 2 k μ k ν ] Π ( k 2 ) = : T μ ν k k 2 Π ( k 2 ) ,
where Π ( k 2 ) is the dimensionless gluon self-energy; hence, the gauge invariance constraint entails
D μ ν ( k ) = T μ ν k 1 k 2 [ 1 + Π ( k 2 ) ] + ξ k μ k ν k 4 = : T μ ν k D ¯ ( k 2 ) + ξ k μ k ν k 4 .
This is the propagator of a massless vector boson, unless
Π ( k 2 ) = k 2 0 m J 2 k 2 ,
in the event of which both the dressed gluon acquires a mass and all symmetry constraints are preserved. That Equation (6) is possible in an interacting quantum gauge field theory was first shown in a study of two-dimensional QED [56,57], and the phenomenon is now known as the Schwinger mechanism of gauge boson mass generation. Three-dimensional QED supports a similar outcome [3,4,5,58,59], as does D = 3 QCD [60,61]; but, as already noted above, there is a difference between both these examples and QCD. Namely, whereas the Lagrangian couplings in D < 4 theories carry a mass dimension, which explicitly breaks scale invariance, this is not the case for D = 4 chromodynamics.
The existence of a Schwinger mechanism in QCD was first conjectured forty years ago [62]. The idea has subsequently been explored and refined [63,64,65,66,67], so that, today, a detailed picture is emerging, which unifies both the gauge and matter sectors [68]. The dynamical origin of the QCD Schwinger mechanism and its intimate connection with nonperturbative dynamics in the three-gluon vertex are elucidated elsewhere [36,37]. This is an area of continuing research, where synergies between continuum and lattice QCD are being exploited [69,70]. For our purposes, it is sufficient to know that Equation (6) is realised in QCD. Indeed, owing to their self-interactions, gluon partons transmogrify into gluon quasiparticles, whose propagation characteristics are determined by a momentum-dependent mass function. That mass function is power-law suppressed in the ultraviolet—hence, invisible in perturbation theory, yet large at infrared momenta, being characterised by a renormalisation-point-independent value [71]:
m 0 = 0.43 ( 1 ) GeV .
The renormalisation-group-invariant (RGI) gluon mass function is drawn in Figure 2.
Before closing this section, it is worth stressing the importance of Poincaré covariance in modern physics.2 If one chooses to formulate a problem in quantum field theory using a scheme that does not ensure Poincaré invariance of physical quantities, then artificial or “pseudodynamical” effects are typically encountered [77]. In connection with gauge theory Schwinger functions, Poincaré covariance very effectively limits the nature and number of independent amplitudes that are required for a complete representation. In contrast, analyses and quantisation procedures that violate Poincaré covariance engender a rapid proliferation in the number of such functions. For instance, the covariant-gauge gluon two-point function in Equation (5) is fully specified by one scalar function, whereas, in the class of axial gauges, two unconnected functions are required, and unphysical, kinematic singularities appear in the associated tensors [78,79]. This is why covariant gauges are normally employed for concrete calculations in both continuum and lattice-regularised QCD. In fact, Landau gauge, i.e., ξ = 0 in Equation (5), is often used because, amongst other things, it is a fixed point of the renormalisation group [53] (Ch. IV) and implemented readily in lQCD [80]. We typically refer to Landau gauge results herein. Naturally, gauge covariance of Schwinger functions ensures that expressions of EHM in physical observables are independent of the gauge used for their elucidation.
Equation (7) is the cleanest expression of EHM in Nature, being truly a manifestation of mass emerging from nothing: infinitely many massless gluon partons fuse together so that, for all intents and purposes, they behave as coherent quasiparticle fields with a long-wavelength mass, which is almost half that of the proton. The implications of this result are enormous and far-reaching, including, e.g., key steps toward the elimination of the problem of Gribov ambiguities [81], which were long thought to prevent a rigorous definition of QCD.

4. Process-Independent Effective Charge

In classical field theories, couplings and masses are constants. Typically, this is also true in quantum mechanics models of strong interaction phenomena. However, this is not the case in renormalisable quantum gauge field theories, as highlighted by the Gell-Mann–Low effective charge/running coupling in QED [82], which is a textbook case [9] (Ch. 13.1).
A highlight of Twentieth Century physics was the realisation that QCD, in particular, and non-Abelian gauge theories, in general, express asymptotic freedom ([83,84,85], [86] (Ch. 7.1)), i.e., the feature that the interaction between charge carriers in the theory becomes weaker as k 2 , the momentum-squared characterising the scattering process, becomes larger. Analysed perturbatively at one-loop order in the modified minimal subtraction renormalisation scheme, M S ¯ , the QCD running coupling takes the form
α M S ¯ ( k 2 ) = γ m π ln k 2 / Λ QCD 2 ,
where γ m = 12 / [ 33 2 n f ] , with n f the number of quark flavours whose mass does not exceed k 2 , and Λ QCD 0.2 GeV is the RGI mass parameter that sets the scale for perturbative analyses.
Asymptotic freedom comes with a “flip side”, which came to be known as infrared slavery ([87] (Section 3.1.2)). Namely, beginning with some k 2 Λ QCD 2 , then the interaction strength grows as k 2 is reduced, with the coupling diverging at k 2 = Λ QCD 2 . (This is the Landau pole.) Qualitatively, this statement is true at any finite order in perturbation theory: whilst the value of Λ QCD changes somewhat, the divergence remains. In concert with the area law demonstrated in Reference [18], which entails that the potential between any two infinitely massive colour sources grows linearly with their separation, many practitioners were persuaded that the complex dynamical phenomenon of confinement could simply be explained by an unbounded potential that grows with parton separation. As we shall see, that is not the case, but the notion is persistent.
Given the character of QCD’s perturbative running coupling, two big questions arise:
(a)
Does QCD possess a unique, nonperturbatively well-defined and calculable effective charge, viz. a veritable analogue of QED’s Gell-Mann–Low running coupling; and
(b)
Does Equation (8) express the large- k 2 behaviour of that charge?
If both questions can be answered in the affirmative, then great strides have been made toward verifying that QCD is truly a theory.
Following roughly forty years of two practically disjoint research efforts, one focused on QCD’s gauge sector [66,67,88] and another on its matter sector [21,23,24,25], a key step on the path to answering these questions was taken in Reference [68]. The two distinct efforts were designated therein as the top-down approach—ab initio computation of the interaction via direct analyses of gauge sector gap equations; and the bottom-up scheme—inferring the interaction by describing data within a well-defined truncation of those matter sector equations that are relevant to bound state properties. Reference [68] showed that the top-down and bottom-up approaches are unified when the RGI running interaction predicted by then-contemporary analyses of QCD’s gauge sector is used to explain ground state hadron observables using nonperturbatively improved truncations of the matter sector bound state equations. The first such truncation was introduced in Reference [89].
It was a short walk from this point to a realisation [90] that in QCD, by means of the pinch technique [65,91,92] and background field method [93], one can define and calculate a unique, process-independent (PI) and RGI analogue of the Gell-Mann–Low effective charge, now denoted α ^ ( k 2 ) . The analysis was refined in Reference [71], which combined modern results from continuum analyses of QCD’s gauge sector and lQCD configurations generated with three domain-wall fermions at the physical pion mass [73,74,75] to obtain a parameter-free prediction of α ^ ( k 2 ) . The resulting charge is drawn in Figure 3. It is reliably interpolated by writing
α ^ ( k 2 ) = γ m π ln K 2 ( k 2 ) / Λ QCD 2 , K 2 ( y = k 2 ) = a 0 2 + a 1 y + y 2 b 0 + y ,
with (in GeV 2 ): a 0 = 0.104 ( 1 ) , a 1 = 0.0975 , b 0 = 0.121 ( 1 ) . The curve was obtained using a momentum subtraction renormalisation scheme: Λ QCD = 0.52 GeV when n f = 4 .
Notably, α ^ ( k 2 ) is PI and RGI in any gauge; but, it is sufficient to know α ^ ( k 2 ) in Landau gauge, ξ = 0 in Equation (5), which is the choice both for easiest calculation and the result in Equation (9). This is because α ^ ( k 2 ) is form-invariant under gauge transformations, as may be shown using identities discussed elsewhere [122], and gauge covariance ensures that any such transformations can be absorbed into the Schwinger functions of the quasiparticles whose interactions are described by α ^ ( k 2 ) [123].
The following physical features of α ^ deserve to be highlighted because they expose a great deal about QCD.
  • Absence of a Landau pole. Whereas the perturbative running coupling, e.g., Equation (8), diverges at k 2 = Λ QCD 2 , revealing the Landau pole, the PI charge is a smooth function on k 2 0 : the Landau pole is eliminated owing to the appearance of a gluon mass scale, Equation (7).
    Implicit in the “screening function”, K ( k 2 ) , is a screening mass:
    ζ H = K ( k 2 = Λ QCD 2 ) 1.4 Λ QCD < m p ,
    at which point the perturbative coupling would diverge but the PI coupling passes through an inflection point on its way to saturation. On k 2 ζ H , the PI charge enters a new domain, upon which the running slows, practically ceasing on k 2 m 0 / 2 , so that QCD is once again effectively a conformal theory and the charge saturates to a constant infrared value α ^ ( k 2 = 0 ) = π × 0.97 ( 4 ) . This value is a prediction: within 3(4)%, the coupling saturates to a value of π at k 2 = 0 . It is not yet known whether this proximity to π has any deeper significance in Nature, but a potential explanation is provided in the next bullet.
    These features emphasise the role of EHM as expressed in Equation (7): the existence of m 0 m p / 2 guarantees that long-wavelength gluons are screened, so play no dynamical role. Consequently, ζ H marks the boundary between soft/nonperturbative and hard/perturbative physics. It is therefore a natural choice for the “hadron scale”, viz. the renormalisation scale at which valence quasiparticle degrees-of-freedom should be used to formulate and solve hadron bound state problems [71]. Implementing that notion, then those quasiparticles carry all hadron properties at ζ = ζ H . This approach is today being used to good effect in the prediction of hadron parton distribution functions (DFs)—see Section 11 and References [124,125,126,127,128,129,130,131,132,133,134,135,136].
  • Match with the Bjorken process-dependent charge. The theory of process-dependent (PD) charges was introduced in References [137,138]: “…to each physical quantity depending on a single scale variable is associated an effective charge, whose corresponding Stückelberg–Peterman–Gell-Mann–Low function is identified as the proper object on which perturbation theory applies.” PD charges have since been widely canvassed [94,139,140].
    One of the most fascinating things about the PI running coupling is highlighted by its comparison with the data in Figure 3, which express measurements of the PD effective charge, α g 1 ( k 2 ) , defined via the Bjorken sum rule [141,142]. The charge calculated in Reference [71] is an essentially PI charge. There are no parameters; and, prima facie, no reason to expect that it should match α g 1 ( k 2 ) . The almost precise agreement is a discovery, given more weight by new results on α g 1 ( k 2 ) [95], which now reach into the conformal window at infrared momenta.
    Mathematically, at least part of the explanation lies in the fact that the Bjorken sum rule is an isospin non-singlet relation, which eliminates many dynamical contributions that might distinguish between the two charges. It is known that the two charges are not identical; yet, equally, on any domain for which perturbation theory is valid, the charges are nevertheless very much alike:
    α g 1 ( k 2 ) α ^ ( k 2 ) = k 2 m 0 2 1 + 1 20 α MS ¯ ( k 2 ) ,
    where α MS ¯ is given in Equation (8). At the c quark current–mass, the ratio is 1.007 , i.e., indistinguishable from unity insofar as currently achievable precision is concerned. At the other extreme, in the far infrared, the Bjorken charge saturates to α g 1 ( k 2 = 0 ) = π ; hence,
    α g 1 ( k 2 ) α ^ ( k 2 ) = k 2 m 0 2 1.03 ( 4 ) .
    Evidently, the PD charge determined from the Bjorken sum rule is, for practical intents and purposes, indistinguishable from the PI charge generated by QCD’s gauge sector dynamics [71,90].
  • Infrared completion. Being process-independent, α ^ ( k 2 ) serves numerous purposes and unifies many observables. It is therefore a good candidate for that long-sought running coupling which describes QCD’s effective charge at all accessible momentum scales [139], from the deep infrared to the far ultraviolet, and at all scales in between, without any modification.
    Significantly, the properties of α ^ ( k 2 ) support the conclusion that QCD is actually a theory, viz. a well-defined D = 4 quantum gauge field theory. QCD therefore emerges as a viable tool for use in moving beyond the SM by giving substructure to particles that today seem elementary. A good example was suggested long ago; namely, perhaps all spin- J = 0 bosons may be [57] “…secondary dynamical manifestations of strongly coupled primary fermion fields and vector gauge fields …”. Adopting this position, the SM’s Higgs boson might also be composite, in which case, inter alia, the quadratic divergence of Higgs boson mass corrections would be eliminated.
Qualitatively equivalent remarks have been developed using light-front holographic models of QCD based on anti-de Sitter/conformal field theory (AdS/CFT) duality [143,144].
Returning to the two questions posed following Equation (8) in Items (a) and (b), it is now apparent that they are answerable in the affirmative: QCD does possess a unique, nonperturbatively well-defined and calculable effective charge whose large- k 2 behaviour connects smoothly with that in Equation (8). These facts provide strong support for the view that QCD is a well-defined 4D quantum gauge field theory.

5. Confinement

Confinement is much discussed, but little understood. In large part, both these things stem from the absence of a clear, agreed definition of confinement. With certainty, it is only known that nothing with quantum numbers matching those of the gluon or quark fields in Equation (1) has ever reached a detector.
An interpretation of confinement is included in the official description of the Yang–Mills Millennium Problem [145]. The simpler background statement is worth repeating:
Quantum Yang–Mills theory is now the foundation of most of elementary particle theory, and its predictions have been tested at many experimental laboratories, but its mathematical foundation is still unclear. The successful use of Yang–Mills theory to describe the strong interactions of elementary particles depends on a subtle quantum mechanical property called the ‘mass gap’: the quantum particles have positive masses, even though the classical waves travel at the speed of light. This property has been discovered by physicists from experiment and confirmed by computer simulations, but it still has not been understood from a theoretical point of view. Progress in establishing the existence of the Yang–Mills theory and a mass gap will require the introduction of fundamental new ideas both in physics and in mathematics.”
The formulation of this problem focuses entirely on quenched-QCD, i.e., QCD without quarks; so, its solution is not directly relevant to our Universe. Confinement in pure quantum SU ( 3 ) gauge theory and in QCD proper are probably very different because the pion exists and is unnaturally light on the hadron scale [27]. On the other hand, the remarks concerning the emergence of a “mass gap” relate directly to Figure 2 and Equation (7) herein. Whilst these properties of QCD may be considered proven by the canons of theoretical physics, such arguments do not meet the standards of mathematical physics and constructive field theory because they involve input from numerical analyses of QCD Schwinger functions. Hereafter, therefore, we will continue within the theoretical physics perspective.
As noted above, a mechanism for the total confinement of infinitely massive charge sources has been identified in the lattice-regularised treatment of quantum field theories using compact representations of Abelian or non-Abelian gauge fields [18], viz. the area law ≡ linear source–antisource potential. However, no treatment of the continuum meson bound state problem has yet been able to demonstrate how such an area law emerges as the masses of the meson’s valence degrees-of-freedom grow to infinity.
In the era of infrared slavery, it was widely assumed that some sort of nonperturbatively improved one-gluon exchange could simultaneously produce asymptotic freedom and a linearly rising potential between quarks; and many models were developed with just such features [146,147,148,149,150]. However, as highlighted by our discussion of QCD’s effective charge, ongoing developments in the study of mesons, using rigorous treatments of the Schwinger functions involved, do not support this picture of confinement via dressed one-gluon exchange. The path to an area law is far more complex.
One direction that deserves exploration is connected with the gluon “H-diagrams” drawn in Reference [151] (Figure 8) and reproduced in Figure 4A. Imagine a valence quark and antiquark scattering via such a process, as drawn in Figure 4B; then keep adding H-diagrams within H-diagrams, exploiting both gluon–quark and gluon self-couplings. Such H-diagram scattering processes produce an infrared divergence in the perturbative computation of a static quark potential [152], viz. a contribution that exhibits unbounded growth as the source–antisource separation increases. Nonperturbatively, that divergence is tamed because the effective charge saturates—Figure 3. On the other hand, there are infinitely many such contributions; and in the limit of static valence degrees-of-freedom, the entire unbounded sum of planar H-diagrams is contracted to a point connection of infinitely dense fisherman’s net/spider’s web diagrams on both the source and antisource. It is conceivable that the confluence of these effects could yield the long-sought area law via the Bethe–Salpeter equation [151,153].
Real-world QCD, however, is characterised by light degrees-of-freedom: u and d quarks with electron-size current–masses; s quarks with a mass roughly one order-of-magnitude larger, so still much less than m p . Pions and kaons are constituted from such valence degrees-of-freedom, and these mesons are light. In fact, the pion has a lepton-like mass [40]: m π m μ , where m μ is the mass of the μ lepton. Owing to the presence of such degrees-of-freedom, light particle annihilation and creation effects are essentially nonperturbative in QCD. Consequently, despite continuing dedicated efforts [154,155,156], it has thus far proven impossible to either define or calculate a static quantum mechanical potential between two light quarks.
This may be illustrated by apprehending that a potential which increases with separation can be described by a flux tube extending between the source and antisource. As the source–antisource separation increases, so does the potential energy stored in the flux tube. However, it can only increase until the stored energy matches that required to produce a particle+antiparticle pair of the theory’s lightest asymptotic states—in QCD, a π + π pair. Numerical simulations of lQCD reveal [157,158] that once the energy exceeds this critical value, the flux tube then dissolves along its entire length, leaving two isolated colour singlet systems. Given that m π = 0.14 GeV, then this disintegration must occur at source+antisource centre-of-mass separation r ( 1 / 3 ) fm [159], which is well within the interior of any hadron. This example assumes that the source and antisource are static. The situation is even more complex for real, dynamical quarks. Thus, at least in the u, d, s quark sector, confinement is manifested in features of Schwinger functions that are far more subtle than can be captured in typical potential models.
One non-static, i.e., dynamical, picture of confinement has emerged from studies of the analytic properties of the two-point Schwinger functions associated with the propagation of coloured gluon and quark quasiparticles—see, e.g., Figure 2. The development of this perspective may be traced back to its beginning almost forty years ago [15,160,161,162,163,164]. It has subsequently been carefully explored [17,22,33,51,81,163,165,166,167,168,169,170]; and in this connection, one may profitably observe that only Schwinger functions which satisfy the axiom of reflection positivity [7,171,172] can be connected with states that appear in the Hilbert space of observables.
The axioms referred to here are those first presented in References [171,172] and subsequently modified in [7], which identify the properties of Schwinger functions that are necessary and sufficient to ensure equivalence between the formulation of a given quantum field theory in Euclidean and Minkowski space. (A contemporary literature compilation is presented elsewhere [173].) In effect, this means that all and only those Schwinger functions which satisfy the five Osterwalder–Schrader axioms possess connections with elements in the Hilbert space of physical states. Regarding strong interactions, all physical states are colour singlets. Consequently, for QCD to be the theory of strong interactions, all its colour singlet Schwinger functions must satisfy the Osterwalder–Schrader axioms; equally, all its colour-nonsinglet functions must violate at least one.
Reflection positivity is a severe constraint. It requires that the Fourier transform of the momentum space Schwinger function, treated as a function of analytic, Poincaré-invariant arguments, is a positive-definite function. To illustrate, consider the gluon Schwinger function in Equation (5). A massless partonic gluon is described by D ¯ ( k 2 ) = 1 / k 2 ; and the 4 D Fourier transform of this function is
d 4 k ( 2 π ) 4 e i k · x 1 k 2 = 1 4 π 2 x 2 > 0 x 2 > 0 .
More generally, regarding two-point functions, viz. those connected with the propagation of elementary excitations in QCD, reflection positivity is satisfied if, and only if, the Schwinger function has a Källén–Lehmann representation. Returning to the gluon Schwinger function in Equation (5), this means one must be able to write
D ¯ ( k 2 ) = 0 d ζ ρ ( ς ) k 2 + ς 2 , ρ ( ς ) > 0 ς > 0 .
Plainly, ρ ( ς ) = δ ( ς ) yields D ¯ ( k 2 ) = 1 / k 2 , i.e., the two-point function for a bare gluon parton. Hence, according to Equations (13) and (14), absent dressing, the gluon parton could appear in the Hilbert space of physical states.
It is important to observe that any function which satisfies Equation (14) is positive-definite itself. Moreover, given Equation (14),
s g n [ d d k 2 ] n D ¯ ( k 2 ) = ( 1 ) n ;
consequently, inter alia, treated as a function of the analytic, Poincaré-invariant variable k 2 , no function with a Källén–Lehmann representation of the form written in Equation (14) can possess an inflection point. Conversely, any function that exhibits an inflection point or, more generally, has a second derivative that changes sign must violate the axiom of reflection positivity [22]; hence, the associated excitation cannot appear in the Hilbert space of observables.
Take another step, and consider the following configuration space Schwinger function ( τ = x 4 , = k 4 ):
Δ ( τ ) = d 3 x d 4 k ( 2 π ) 4 e i k · x D ¯ ( k 2 ) = 1 π 0 d cos ( τ ) D ¯ ( 2 ) .
Suppose that interactions generate a constant mass for the gluon parton, so that D ¯ ( k 2 ) = 1 / ( k 2 + μ 2 ) . Does that trigger confinement? The answer is “no” because this Schwinger function has a spectral representation with ρ ( ς ) = δ ( ς 2 μ 2 ) ; Equation (15) is satisfied; and so is positivity:
Δ ( τ ) = 1 2 μ e μ τ .
Suppose instead that interactions produce a momentum-dependent mass-squared function like that in Figure 2, which is 1 / k 2 suppressed in the ultraviolet:
m J 2 ( k 2 ) = μ 0 4 k 2 + μ 0 2 D ¯ ( k 2 ) = k 2 + μ 0 2 k 2 ( k 2 + μ 0 2 ) + μ 0 4 .
The mass function itself is a monotonically decreasing, concave-up function; yet, in this case, the Schwinger function has an inflection point at k 2 = 0.53 μ 0 2 . Hence, it does not have a Källén–Lehmann representation; so, the associated excitation cannot appear in the Hilbert space of observables. Furthermore, evaluation of the configuration space Schwinger function defined by Equation (16) yields [81]
Δ ( τ ) = 1 μ 0 e τ μ 0 3 2 cos μ 0 τ 2 = : Δ p ( τ ) cos μ 0 τ 2 ,
using which the curve in Figure 5 is drawn: plainly, the configuration space Schwinger function violates reflection positivity. (Notably, the algebraic calculation of Δ ( τ ) is often difficult and not always possible; so, uniform positivity of the second-derivative, Equation (15), is a much quicker means of testing for reflection positivity. Nevertheless, when it can be obtained, an explicit form of Δ ( τ ) does provide additional insights.)
It is interesting to extend Equation (18) using m J 2 ( k 2 ) = α μ 0 4 / [ k 2 + μ 0 2 ] , in which case
D ¯ ( k 2 = μ 0 2 y ) = 1 μ 0 2 1 + y y 2 + y + α = : 1 μ 0 2 d ( y ) .
This type of Schwinger function lies within the so-called “refined Gribov–Zwanziger” class [166]. For α > 1 2 , as in the example above, the function in Equation (20) exhibits an inflection point at some y > 0 ; and when α ( 1 4 , 1 2 ) , the inflection point is found at a location y ( 1 2 , 0 ) . With α ( 0 , 1 4 ) , on the other hand, the function in Equation (20) separates into a sum of two terms:
d ( y ) = n 1 y + p 1 n 2 y + p 2 ,
where p 1 , 2 R and s g n ( p 1 / p 2 ) = + 1 . In this case, there is no inflection point; nevertheless, the second derivative does change sign, switching ± as y passes through the pole locations. Consequently, any excitation whose propagation is described by a Schwinger function obtained with α > 0 in Equation (20) cannot appear in the Hilbert space of observable states.
Inserting Equation (20) into Equation (16), one finds, after some careful algebra [81] ( a = α 1 / 4 ):
Δ ( τ ) = 1 2 μ 0 a s φ / 2 e τ μ 0 a c φ / 2 × ( 1 + 1 a 2 ) s φ / 2 cos ( τ μ 0 a s φ / 2 ) ( 1 1 a 2 ) c φ / 2 sin ( τ μ 0 a s φ / 2 ) ,
where c φ / 2 = [ 1 2 + 1 4 a 2 ] 1 / 2 , s φ / 2 = [ 1 c φ / 2 2 ] 1 / 2 . (Equation (19) is a special case of this result.) Equation (22) reveals that the distance from τ = 0 of the first zero in the configuration space Schwinger function, τ z , increases with decreasing α , i.e., as the infrared value of the gluon mass-squared function is reduced. Thus, confinement would practically be lost if τ z were to become much greater than π / μ 0 . Considering realistic gluon two-point functions, however, one finds μ 0 ( 3 / 4 ) m p , α 1 [81]; so, τ z 1 fm, expressing a natural confinement length scale. (It is worth observing that if τ z were instead measured on the Å scale, then the notion of confinement would be lost because modern detectors are able to directly image targets of this size [174].)
This discussion is readily summarised. Owing to complex nonlinear dynamics in QCD, gluon and quark partons acquire momentum-dependent mass functions, as a consequence of which they emerge as quasiparticles whose propagation characteristics are described by two-point Schwinger functions that are incompatible with reflection positivity. Normally, the dynamical generation of running masses is alone sufficient to ensure this outcome. It follows that the dressed gluons and quarks cannot appear in the Hilbert space of physical states. In this sense, they are confined. The associated confinement length scale is τ z 1 fm. It is worth stressing that the use of such two-point functions in the calculation of colour singlet matrix elements ensures the absence of coloured particle+antiparticle production thresholds [175], thereby securing the empirical expression of real-QCD confinement.
Considering these quasiparticle Schwinger functions further, one may also define a parton persistence or fragmentation length, τ F , as the scale whereat the deviation of the Schwinger function from parton-like behaviour is 50%: Δ ( τ F ) / Δ p ( τ F ) = 0.5 . Referring to Equation (19), one reads τ F = ( 2 / 3 ) τ z . This result is also found using realistic gluon two-point functions [81]. (The value of 50% is merely a reasonable choice. At this level, 30% would also be acceptable, in which case τ F τ F = ( 1 / 2 ) τ z .)
A physical picture of dynamical confinement now becomes apparent [165]. Namely, once a gluon or quark parton is produced, it begins to propagate in spacetime. However, after traversing a spacetime distance characterised by τ F , interactions occur, causing the parton to lose its identity, sharing it with others. Ultimately, combining the effects on this parton with similar impacts on those produced along with it, a countable infinity of partons (a parton cloud) is produced, from which detectable colour singlet final states coalesce. This train of events is the physics expressed in parton fragmentation functions (FFs) [176]. Such distributions describe how the QCD partons in Equation (1), generated in a high-energy event and almost massless in perturbation theory, transform into a shower of massive hadrons, viz. they describe how hadrons with mass emerge from practically massless partons. It is natural, therefore, to view FFs as the cleanest expression of dynamical confinement in QCD. Furthermore, in the neighbourhood of their common boundary of support, DFs and FFs are related by crossing symmetry [177]: FFs are timelike analogues of DFs. Hence, an understanding of FFs and their deep connection with DFs can deliver fundamental insights into EHM. This picture of parton propagation, hadronisation, and confinement—of DFs and FFs—can be tested in experiments at modern and planned facilities [178,179,180,181,182,183,184,185,186]. A pressing demand on theory is delivery of predictions for FFs before such experiments are completed so as, for instance, to guide development of facilities and detectors. As yet, however, there are no realistic computations of FFs. In fact, even a formulation of this problem remains uncertain.
Before moving on, it is worth reiterating that confinement means different things to different people. Whilst some see confinement only in an area law for Wilson loops [18], our perspective stresses a dynamical picture, in which dynamically driven changes in the analytic structure of coloured Schwinger functions ensures the absence of colour-carrying objects from the Hilbert space of observable states. In time, perhaps, as strong QCD is better understood, it may be found that these two realisations are connected. The only certain thing is the necessity to keep an open mind on this subject.

6. Spectroscopy

Insofar as the spectrum of hadrons is concerned, results from nonrelativistic or somewhat relativised quark models [187,188,189] are still often cited as benchmarks. Indeed, a standard reference ([40] (Section 63)) includes the following assertions: “The spectrum of baryons and mesons exhibits a high degree of regularity. The organizational principle which best categorizes this regularity is encoded in the quark model. All descriptions of strongly interacting states use the language of the quark model.” This is despite the facts that neither the “quarks” nor the potentials in quark models have been shown to possess any mathematical link with Equation (1)—rigorous or otherwise; and, furthermore, the orbital angular momentum and spin used to label quark model states are not Poincaré-invariant (observable) quantum numbers.
In step with improvements in computer performance, lQCD is delivering interesting results for hadron spectra [190,191], amongst which one may highlight indications for the existence of hybrid and exotic hadrons [192,193,194,195]. Continuum studies in quantum field theory are lagging behind owing in part to impediments placed by the character of the Bethe–Salpeter equation; primarily the fact that it is impossible to write the complete Bethe–Salpeter kernel in a closed form.
A systematic approach to truncating the integral equations associated with bound state problems in QCD was introduced almost thirty years ago [196,197]. Amongst other things, the scheme highlighted the importance of preserving continuous and discrete symmetries when formulating bound state problems; enabled proof of Goldberger–Treiman identities and the Gell-Mann–Oakes–Renner relation in QCD [198,199]; and opened the door to symmetry-preserving, Poincaré-invariant predictions of hadron observables, including elastic and transition form factors and DFs [21,23,24,25,28,125,200,201,202,203,204,205,206,207]. Some of the more recent developments are sketched below.
An issue connected with the leading-order (rainbow ladder (RL)) term in the truncation scheme of References [196,197] is that it only serves well for those ground state hadrons which possess little rest-frame orbital angular momentum, L, between the dressed valence constituents [208,209,210,211,212,213,214,215,216,217,218]. This limitation can be traced to its inability to realistically express the impacts of EHM on hadron observables, a weakness that is not overcome at any finite order of elaboration [210]. Improved schemes, which express EHM in the kernels, have been identified [68,89,151,219,220,221]. They have shown promise in applications to ground state mesons constituted from u and d valence quarks and/or antiquarks. However, that is a small subset of the hadron spectrum; so, a recent extension to the spectrum and decay constants of u, d, and s meson ground and first-excited states is welcome [222].
Returning to quark models, it was long ago claimed [149] “…that all mesons—from the pion to the upsilon—can be described in a unified framework.” The context for this assertion was a model potential built using one-gluon-like exchange combined with an infrared slavery “confinement” term that increases linearly with colour source separation. The basic mass scales in such potential models are set by the constituent quark masses; and one might draw a qualitative link between those scales and the far-infrared values of the momentum-dependent dressed quark running masses ([35] (Figure 2.5)): M u , d ( 0 ) 0.41 GeV, M s ( 0 ) 0.53 GeV. Thereafter, mass splittings and level orderings are arranged by tuning details of the potential. Such a procedure can be quantitatively efficacious; however, it is qualitatively incorrect. This is readily seen by recalling the Gell-Mann–Oakes–Renner relation [46,198,223]: m π 2 m ^ , where m ^ is Nature’s explicit source of chiral symmetry breaking, generated by Higgs boson couplings to quarks in the SM. Such behaviour is impossible in a potential model [25,27], but natural in the CSM treatment of bound states—see, e.g., References [198,208], [21] (Figure 3.3), [224] (Figure 7), [89] (Figure 1A). Thus, whilst potential models might deliver a fit to hadron spectra, they do not provide an explanation.
That such challenges are surmounted when using CSMs to treat hadron bound state problems is further exemplified in Reference [222], which adapts a novel scheme for including EHM effects in the Bethe–Salpeter kernel [221] to simultaneously treat ground and first-excited states of u, d, and s quarks. As revealed in Figure 6, the empirical spectrum displays some curious features, e.g.: consistent with quark mass–scale counting, m ρ < m K * , but this ordering is reversed for the first excitations of these states; the first excited state of the π is lighter than that of the ρ , but the ordering is switched for K, K * ; and all axial vector mesons are nearly degenerate, with the larger mass of the s quarks appearing to have little or no impact. In delivering the first symmetry-preserving analysis of this collection of states to employ an EHM-improved kernel, Reference [222] supplies fresh insights into the dynamical foundations of the properties of lighter quark mesons.
In order to sketch that effort, we note that, using CSMs, the dressed propagator (two-point Schwinger function) for a quark with flavour g is obtained as the solution of the following gap equation:
S g 1 ( k ) = i γ · k A g ( k 2 ) + B g ( k 2 ) = [ i γ · k + M g ( k 2 ) ] / Z g ( k 2 ) ,
= Z 2 ( i γ · k + m g ) + Σ g ( k ) ,
Σ ( k ) = Z 1 d q Λ g 2 D μ ν ( k q ) λ a 2 γ μ S g ( q ) λ a 2 Γ ν g ( k , p ) ,
where M g ( k 2 ) is RGI;
D μ ν ( k ) = Δ ( k 2 ) 0 D μ ν ( k ) ;
Γ ν g is the quark-gluon vertex; d q Λ denotes a Poincaré-invariant regularisation of the four-dimensional integral, with Λ the regularisation mass scale; and Z 1 , 2 ( ζ 2 , Λ 2 ) are, respectively, the vertex and quark wave function renormalisation constants, with ζ the renormalisation point. (In such applications, when renormalisation is necessary, a mass-independent scheme is important, as discussed elsewhere [225].)
It was anticipated almost forty years ago [226,227], and confirmed more recently [76,228,229,230], that EHM engenders a large anomalous chromomagnetic moment (ACM) for the lighter quarks; and with development of the first EHM-improved kernels, it was shown that such an ACM has a big impact on the u, d meson spectrum [219,220,221].
The aim in Reference [222] was to extend References [219,220,221] and highlight additional impacts of an ACM on the spectrum of mesons constituted from u, d, and s quarks. An ACM emerges as a feature of the dressed gluon–quark vertex, a three-point Schwinger function. Its character and impacts can be exposed by writing ( l = k q )
Γ ν g ( q , k ) = γ ν + τ ν ( l ) , τ ν ( l ) = η κ ( l 2 ) σ l ν ,
σ l ν = σ ρ ν l ρ , κ ( l 2 ) = ( 1 / ω ) exp ( l 2 / ω 2 ) . Here, η > 0 is the strength of the ACM term, τ ν ( l ) ; and it is assumed that the vertex is flavour-independent, which is a sound approximation for the lighter quarks [224,231]. When considering Equation (25), one might remark that the complete gluon–quark vertex is far more complicated—potentially containing twelve distinct terms—and, in QCD, κ ( l 2 ) is power-law suppressed in the ultraviolet. Notwithstanding these things, illustrative purposes are well served by Equation (25).
ACM effects are most immediately felt by the dressed quark propagator. The presence of an ACM in the kernel of Equation (23) increases positive EHM-induced feedback on dynamical mass generation. Consequently, as shown elsewhere [76], realistic values of the dressed quark mass at infrared momenta are achieved using the PI effective charge in Figure 3. Such an outcome requires tuning when using the PI charge in a rainbow truncation of the gap equation; in fact, DCSB cannot be guaranteed in that case [68].
Following References [89,221,228] in continuing to emphasise clarity over numerical complexity, Reference [222] also simplified the kernel in Equation (23), writing
g 2 Δ ( k 2 ) = 4 π α ^ ( 0 ) D ( η ) ω 4 e k 2 / ω 2 ,
where ω = 0.8 GeV, a value matching that suggested by analyses of QCD’s gauge sector [68,71], and D ( η ) = D RL ( 1 + 0.27 η ) / ( 1 + 1.47 η ) , with ω D RL = ( 1.286 GeV ) 3 chosen to achieve m ρ = 0.77 GeV in RL truncation. The η -dependence of D ( η ) was fixed a posteriori by requiring that m ρ remain unchanged as η is increased. Since η > 0 adds EHM strength to the gap equation’s kernel, then D must become smaller as η grows in order to maintain a fixed value of m ρ . Following this procedure, m ρ becomes the benchmark against which all ACM-induced changes are measured.
It is worth noting that when one identifies ( g 2 / [ 4 π ] ) Δ ( k 2 = 0 ) = 1 / μ 0 2 , then μ 0 = 0.39 GeV in RL truncation and μ 0 = 0.57 GeV at η = 1.2 . Therefore, the interaction specified by Equation (26) is consistent with gluon mass generation, as it is described in Section 3. On the other hand, the large- k 2 behaviour of Equation (26) does not respect the renormalisation group flow of QCD. This would be an issue if one were using it, e.g., to calculate hadron form factors at large- Q 2 [200,201,203,204], where Q 2 is momentum transfer squared, or parton distribution functions and amplitudes near the endpoints of their support domains [124,125,126,127].3 However, it is far less important when calculating global, integrated properties, like hadron masses. In such applications, Equation (26) is satisfactory. Indeed, good results can even be obtained using a symmetry-preserving treatment of a momentum-independent interaction [232,233,234]. A key merit of Equation (26) lies in its elimination of the need for renormalisation, which simplifies analyses without materially affecting relevant results.
The scheme introduced in Reference [221] provides a direct route from any reasonable set of gap equation elements to closed-form Bethe–Salpeter kernels for meson bound state problems. Thus having specified physically reasonable gap equations via Equations (25) and (26), Reference [222] adapted that scheme to arrive at Bethe–Salpeter equations for each of the mesons identified in Figure 6, obtaining solutions in all cases. The image compares the experiment with the RL truncation results, also calculated in Reference [222], and predictions obtained using the EHM-improved kernel.4
It is first worth mentioning the RL truncation mass predictions in Figure 6. On the whole, the mean absolute relative difference, ard ¯ , between RL results and central experimental values is 13 ( 8 ) %. This is tolerable. However, there is substantial scatter and there are many qualitative discrepancies.
In contrast, compared with central experimental values, the EHM-improved masses in Figure 6 agree at the level of ard ¯ = 2.9 ( 2.7 ) %. This is a factor of 4.6 of improvement over the RL spectrum. Moreover, correcting RL truncation flaws and reproducing empirical results: m K > m π , m ρ > m π , m ρ m K * ; the mass splittings a 1 - ρ and b 1 - ρ match the empirical values because including the ACM in the kernel has markedly increased the masses of the a 1 and b 1 mesons, whilst m ρ was deliberately kept unchanged; m ϕ m ϕ agrees with experiment to within 2%; the K 1 + , K 1 + + -level order is correct; and quark+antiquark scalar mesons are heavy, providing room for the addition of strongly attractive resonant contributions to the bound state kernels [235,236].
Using the Bethe–Salpeter amplitudes obtained in solving for the meson spectrum, canonically normalised in the standard fashion ([237] (Section 3)), Reference [222] also delivered predictions for the entire array of associated leptonic decay constants, f H , including many that have not yet been measured. The predicted values are depicted in Figure 7, which also includes the few available empirical results. The ground state leptonic decay constants in Figure 7 were calculated directly on-shell, but extrapolation was necessary to obtain on-shell values for those of the excited states. For these observables, two extrapolation schemes were used and they yielded consistent results in all cases.
Given that Reference [222] used a simplified interaction, viz. Equation (26), then the Figure 7 comparison between predicted ground state decay constants and the few known empirical values is favourable, particularly because decay constants are sensitive to ultraviolet physics, which was omitted. There are also indications that the EHM-improved kernels deliver better agreement.
The decay constants of radially excited states are especially interesting. Quantum mechanics models of positronium-like systems produce a single zero in the radial wave function of n = 1 states. The decay constant of a first radial excitation is thus ( 1 / 8 ) -times that of the ground state. The predictions in Reference [222] are generally consistent with this pattern, except for J P = 0 mesons. In the pseudoscalar channel, as a corollary of EHM, QCD predicts f H n = 1 0 in the chiral limit [208,209,240,241]. The results in Reference [222] meet this requirement, whereas such outcomes cannot be achieved in quark models without tuning parameters. For this reason alone, the decay constant predictions of [222] warrant testing.
Notwithstanding the simplifications used in formulating the problem, Reference [222] delivered the first Poincaré-invariant analysis of the spectrum and decay constants of the u, d, and s meson ground and first-excited states. The results include predictions for masses of as-yet unseen mesons and many unmeasured decay constants. One may look forward to extensions of the approach to heavy+light mesons [239,242,243], hybrid/exotic mesons [213,244,245,246], and glueballs [247,248,249,250]. These directions are especially important owing to worldwide investments in studies of the former and searches for the latter [179,183,186,251,252,253].
Such progress with meson properties should not obscure the need to calculate the spectrum of baryons. Indeed, baryons are the most fundamental three-body systems in Nature; if we do not understand how QCD, a Poincaré-invariant quantum field theory, structures the spectrum of baryons, then we do not understand Nature. Within the context of the truncation scheme introduced in References [196,197], baryon masses and bound state amplitudes have been calculated using a Poincaré-covariant Faddeev equation that describes a six-point Schwinger function for three-quark three-quark scattering. The first solution of this problem for the nucleon (N) was presented in Reference [254]; continuing studies are reviewed elsewhere [28,31,181,255]; efforts are now under way to adapt the methods in References [221,222] to the formulation and solution of the baryon Faddeev equations.
Meanwhile, the quark + dynamical diquark approach to baryon properties, introduced in References [256,257,258,259], is also being pursued vigorously. This treatment begins with solutions of the equation illustrated in Figure 8. As sketched, e.g., in Reference [21] (Section 5.1), this is an approximation to the three-body Faddeev equation, whose kernel is constructed using dressed quark and nonpointlike diquark degrees-of-freedom. Binding energy is lodged within the diquark correlation and also produced by the exchange of a dressed quark, which, as drawn in Figure 8, emerges in the break-up of one diquark and propagates to be absorbed into formation of another. In the general case, five distinct diquark correlations are possible: isoscalar–scalar, ( I , J P = 0 , 0 + ) ; isovector–axial vector; isoscalar–pseudoscalar; isoscalar–vector; isovector–vector. Channel dynamics within a given baryon determines the relative strengths of these correlations therein.
Given the extensive coverage of the role of diquark correlations in hadron structure presented in Reference [255], herein, we will only draw some recent highlights from analyses of the baryon spectrum using the Faddeev equation in Figure 8, drawing largely from Reference [232]. That study was built upon a symmetry-preserving treatment of a vector×vector contact interaction (SCI), which was introduced a little over a decade ago [260] and has since been employed with success in numerous applications, some of which are reviewed in this volume [261]. Amongst the merits of the SCI are its algebraic simplicity; limited number of parameters; simultaneous applicability to many systems and processes; and potential for generating insights that connect and explain numerous phenomena.
Reference [232] used the SCI to calculate the ground state masses of J P = 0 ± , 1 ± ( f g ¯ ) mesons and J P = 1 / 2 ± , 3 / 2 ± ( f g h ) baryons, where f , g , h { u , d , s , c , b } . Using J P = 1 2 ± states as exemplars, Figure 9 highlights the level of quantitative accuracy.
Regarding the 33 mesons, then ard ¯ = 5 ( 6 ) % when comparing SCI predictions with empirical masses. In achieving this outcome, it was found that sound expressions of EHM were crucial. Turning to baryons, the SCI generated 88 distinct bound states, namely every possible three-quark 1 / 2 ± , 3 / 2 ± ground state. In this collection, 34 states have already been identified empirically, and lQCD results are available for another 30: for these 64 states, comparing SCI prediction with experiment, where available, or lQCD mass otherwise, ard ¯ = 1.4 ( 1.2 ) %. This level of agreement was only achieved through the implementation of EHM-induced effects associated with spin–orbit repulsion in 1 / 2 baryons. Notably, the same 88 ground states are also produced by a three-body Faddeev equation [218]: in comparison with those results, ard ¯ = 3.4 ( 3.0 ) %.
Overall, Reference [232] delivered SCI predictions for 164 distinct quantities, 114 of which have either been measured or calculated using lQCD: performing a comparison on this subset yields ard ¯ = 4.5 ( 7.1 ) %. Such quantitative success means that credibility should be given to the qualitative conclusions that follow from the SCI analysis. We list them here: (i) Nonpointlike, dynamical diquarks play a significant role in all baryons. Usually, the lightest allowed diquark is the most important part of a baryon’s Faddeev amplitude. This remains true, even if the lightest correlation is a (sometimes called bad) axial vector diquark, and also for baryons containing one or more heavy quarks. In the latter connection, this means one cannot safely assume that singly heavy baryons may realistically be described as two-body light diquark+heavy quark ( q q + Q ) bound states or that doubly heavy baryons ( q Q Q ) can be treated as two-body light quark+heavy diquark bound states, q + Q Q . Corresponding statements apply to the treatment of tetra- and penta-quark problems. (ii) Positive-parity diquarks dominate in positive-parity baryons. Axial vector diquarks are prominent in all states. (iii) Negative-parity diquarks play a minor role in positive-parity baryons. On the other hand, owing to EHM, they are significant and sometimes dominant in J = 1 / 2 baryons. (iv) Curiously, however, J = 3 / 2 baryons are built (almost) exclusively from J = 1 + diquark correlations. These conclusions are being checked using Faddeev equations with momentum-dependent exchange interactions; hence, a closer link to QCD. Where results are already available, the SCI conclusions have been confirmed [265,266,267,268].
Following more than fifty years of hadron spectroscopy based on quark models, we are beginning to see real progress with the use of bound state equations in quantum field theory. Poincaré-invariant, symmetry-preserving analyses that reveal the expressions of EHM in hadron masses and level orderings are becoming available. This increases the value of experimental hadron spectra measurements, making them a clearer window onto strong QCD.

7. Baryon Wave Functions

Concerning baryon structure, as noted when opening Section 6, quark models are still considered to provide a useful picture [187,188,189,269]. In such models, baryons built from combinations of u, d, and s valence quarks are grouped into multiplets of SU ( 6 ) O ( 3 ) . The multiplets are labelled by their flavour content—SU ( 3 ) , spin—SU ( 2 ) , and orbital angular momentum—O ( 3 ) . However, as has been emphasised, quark potential models do not have an explicit link with QCD, a Poincaré-invariant quantum gauge field theory.
For the lightest four ( I , J P = 1 2 , 1 2 ± ) baryons, with I denoting isospin, a comparison between quark model expectations and insights drawn from solutions of the Poincaré-covariant Faddeev equation is presented elsewhere [265]. Herein, we will illustrate the qualitative character of such comparisons by considering more recent studies of ( 3 2 , 3 2 ± ) , ( 1 2 , 3 2 ) baryons [267,268].
These systems were studied in Reference [217] using RL truncation and direct calculations of all primary Schwinger functions. With current algorithms, owing to singularities that enter the integration domains sampled by the Faddeev equations [199], this approach limits the ability to compute wave functions because the on-shell point for many systems is inaccessible. (Similar issues are encountered with meson structure studies—Section 6.) To circumvent this issue, References [267,268] employed the QCD-kindred framework introduced elsewhere [270], in which, instead of calculating all primary Schwinger functions, one uses physics-constrained algebraic representations of the Faddeev kernel elements. This weakens the connection with QCD, but that loss is well compensated because, with reliably informed choices for the representation functions, the expedient enables access to on-shell baryon wave functions. The QCD-kindred framework has widely been used with success—see, e.g., References [271,272,273,274,275,276].
Within the quark model framework and using standard spectroscopic notation, n 2 S + 1 L J , where n is the radial quantum number with “0” labelling the ground state, the lightest four ( 3 2 , 3 2 ± ) Δ baryons, constructed from isospin I = 3 2 combinations of three u and/or d quarks, are understood as follows [277]:
  • Δ ( 1232 ) 3 2 + 0 4 S 3 2 = S -wave ground state;
  • Δ ( 1600 ) 3 2 + 1 4 S 3 2 = S -wave radial excitation of Δ ( 1232 ) 3 2 + ;
  • Δ ( 1700 ) 3 2 0 2 P 3 2 = P -wave orbital angular momentum excitation of Δ ( 1232 ) 3 2 + ;
  • Δ ( 1940 ) 3 2 1 4 P 3 2 = P -wave excitation of Δ ( 1600 ) 3 2 + .
Analogously, the ( 1 2 , 3 2 ) states are interpreted thus [277]:
  • N ( 1520 ) 3 2 0 2 P 1 2 = P -wave ground state in this channel and an angular momentum coupling partner of N ( 1535 ) 1 2 ;
  • N ( 1700 ) 3 2 0 4 P 3 2 = P -wave angular momentum coupling partner of N ( 1520 ) 3 2 ;
  • N ( 1720 ) 3 2 + 0 2 D 3 2 = D -wave orbital angular momentum excitation of N ( 1520 ) 3 2 ;
  • N ( 1900 ) 3 2 + 0 4 D 3 2 = D -wave orbital angular momentum excitation of N ( 1700 ) 3 2 .
On the other hand, Poincaré-invariant quantum field theory does not readily admit such assignments. Instead, the states appear as poles in the six-point Schwinger functions associated with the given ( I , J P ) channels. Here, “ ( 1 3 ) ” and “ ( 2 4 ) ’ in each block above are related as parity partners. All differences between positive- and negative-parity states can be attributed to chiral symmetry breaking in quantum field theory. This is highlighted by the ρ - a 1 meson complex [219,220,221,222,278]. Regarding light quark hadrons, such symmetry breaking is almost exclusively dynamical [24,279,280,281,282,283,284]. As noted above, DCSB is a corollary of EHM [32,33,34,35,36,37]; hence, quite probably linked tightly with confinement—Section 5. These features imbue quantum field theory analyses of ( 3 2 , 3 2 ± ) , ( 1 2 , 3 2 ) baryons with particular interest; consequently, experiments that test predictions made for structural differences between parity partners are highly desirable.
Working with the Faddeev equation sketched in Figure 8, then, a priori, the ( 3 2 , 3 2 ± ) baryons are the simpler systems because they can only contain isovector–axial vector and isovector–vector diquarks, whereas the ( 1 2 , 3 2 ) systems may involve all five distinct types of diquarks: ( 0 , 0 + ) , ( 1 , 1 + ) , ( 0 , 0 ) , ( 1 , 1 ) , ( 0 , 1 ) . Nonetheless, the formulation of the bound state problems in both channels is practically identical, using the same dressed quark and diquark propagators; diquark correlation amplitudes; etc. This way, one guarantees a unified description of all states in the spectrum. The propagators are parametrised using entire functions [163,285,286]; hence, satisfy the confinement constraints described in Section 5. It is this feature that enables on-shell calculations for all baryons.
The calculated spectrum of states is displayed in Figure 10. As highlighted elsewhere [29,265,267,268,287,288], the kernel in Figure 8 does not include contributions that may be understood as meson–baryon final state interactions. These are the interactions that transform a bare baryon into the observed state, e.g., via dynamical coupled channels calculations [289,290,291,292,293]. The Faddeev amplitudes and masses calculated in References [265,267,268] should therefore be seen as describing the dressed quark core of the bound state, not the fully dressed, observable object [294,295,296]. That explains why the masses are uniformly too large. Evidently and importantly, in each sector, a single subtraction constant is sufficient to realign the masses and produce a good description of the spectrum.
Regarding ( 3 2 , 3 2 ± ) baryons, Reference [267] found that although these states may contain both ( 1 , 1 + ) and ( 1 , 1 ) quark+quark correlations, one can neglect the ( 1 , 1 ) diquarks because they have practically no impact on the masses or wave functions. After this simplification, the Poincaré-covariant wave functions involve eight independent matrix-valued terms, each multiplied by a scalar function of two variables: ( k 2 , k · Q ) , with k the quark+diquark relative momentum. Studying the properties of these functions, one may conclude that Δ ( 1600 ) 3 2 + is fairly interpreted as a radial excitation of the Δ ( 1232 ) 3 2 + , as suggested by the quark model. However, the wave functions of the Δ ( 1700 ) 3 2 and Δ ( 1940 ) 3 2 states are complicated and do not readily admit direct analogies with quark model pictures.
Projecting the Poincaré-covariant Faddeev wave functions of ( 3 2 , 3 2 ± ) baryons into their respective rest-frames, one arrives at a J = L + S separation which is comparable to that associated with quark models. (Here, L is quark–diquark orbital angular momentum.) Following this procedure, Reference [267] found that the angular momentum structure of all these states is much more complicated than is typically generated in quark models—see Figure 11.
Evidently, making a link to quark models, the Δ ( 1232 ) 3 2 + and Δ ( 1600 ) 3 2 + are characterised by a dominant S -wave component and the Δ ( 1700 ) 3 2 by a prominent P -wave. The Δ ( 1940 ) 3 2 , however, does not fit this picture: contrary to quark model expectations, indicated on page 79, this state is S -wave-dominated. Moreover, each state contains every admissible partial wave.
Combining all gathered information, Reference [267] furthermore concluded that the negative-parity Δ baryons are not merely orbital angular momentum excitations of positive parity ground states. In this observation, the results match those obtained earlier for ( 1 2 , 1 2 ± ) baryons [265].
Recalling now that the interpolating fields for positive- and negative-parity hadrons are related by chiral rotation of the quark spinors used in their construction, then the highlighted structural differences are largely generated by DCSB. Regarding the Δ ( 1940 ) 3 2 in particular, these novel structural predictions may be expected to encourage new experimental efforts aimed at extracting reliable information about this little-understood state from exclusive π + π p electroproduction data [297,298].
Observations upon similar features were made about ( 1 2 , 3 2 ± ) baryons in Reference [268]. To begin, despite the fact that such states may contain all possible diquark correlations, the analyses showed that a good approximation is obtained by keeping only ( 0 , 0 + ) , ( 1 , 1 + ) correlations. This runs counter to the nature of ( 1 2 , 1 2 ) systems, in which ( 0 , 0 ) , ( 0 , 1 ) diquarks are also important [265,299]. Projecting the Poincaré-covariant Faddeev wave functions into the baryon rest-frames and considering the baryon mass fraction contributed by each partial wave, this collection of states form a set related via orbital angular momentum excitation: as in quark models, the parity-negative states are primarily P -wave in nature, whereas the parity-positive states are D wave—see Reference [268] (Figure 4). However, looking with finer resolution, using charts of the canonical normalisation constant contributions from the various partial waves in the Poincaré-covariant wave functions, like Figure 11 herein, far greater L-complexity was observed than is usually found in quark models—see Reference [268] (Figure 7). Here, too, one may anticipate that these structural predictions can be tested using data from measurements of resonance electroexcitation at large momentum transfers. For the N ( 1520 ) 3 2 , data are already available [300,301,302,303,304,305], and calculations of the electroproduction form factors are underway. Large- Q 2 data on the other states is not available; so, the predictions in Reference [268] will also encourage new experimental efforts in this area.
Parity partner channels are identical when chiral symmetry is restored [14,30]. It is therefore interesting to note that the mass splitting between partner states does not exhibit a simple pattern, viz. empirically [40]:
states mass splitting / GeV N ( 1535 ) 1 2 N ( 940 ) 1 2 + 0.57 , N ( 1650 ) 1 2 N ( 1440 ) 1 2 + 0.29 , Δ ( 1700 ) 3 2 Δ ( 1232 ) 3 2 + 0.46 , Δ ( 1940 ) 3 2 Δ ( 1600 ) 3 2 + 0.44 , N ( 1720 ) 3 2 + N ( 1520 ) 3 2 0.17 , N ( 1900 ) 3 2 + N ( 1700 ) 3 2 0.22 .
This system dependence of the mass splitting is also linked to the deeper structural differences between these states that are expressed in their complex wave functions.
Using a familiar quantum mechanics framework, quark models produce baryon wave functions that have an appealing simplicity. However, far richer structures are found when quantum field theory is used to solve baryon bound state problems. The growing body of quantum field theory predictions can be tested, e.g., in modern and future large- Q 2 measurements of baryon elastic and transition form factors. In fact, the large- Q 2 character of such experiments is alone sufficient to demand the sort of Poincaré-invariant, symmetry-preserving treatment that only analyses in quantum field theory can deliver. One may, therefore, expect studies using the Faddeev equation approach outlined in this section to become steadily more widespread.

8. Meson Form Factors

The truncation scheme explained in References [196,197] has been used to calculate many meson elastic and transition form factors [21]; and modern algorithms have enabled predictions to be delivered on the entire domain of spacelike Q 2 [200,201,306,307,308,309], making it possible to draw connections with hard scattering formulae derived using QCD perturbation theory [310,311,312]. These new predictions, which unify the infrared and ultraviolet Q 2 domains, are providing the impetus for measurements at new-generation high-energy and high-luminosity facilities [35,180,182,184]. Recalling that QCD is not found in form factor scaling, but in scaling violations, then the goal of these new experiments is, of course, to discover the breakaway from scaling in a hard exclusive process and thus reveal the hand of QCD.
In connection with these new facilities, it has been argued that the interaction of a heavy vector meson, V = J / ψ , Υ , with a proton, p, may provide access to a QCD van der Waals interaction, produced by multiple gluon exchange [313,314] and/or the QCD trace anomaly [6,315,316,317]. The van der Waals interaction is of interest because it might relate to, amongst other things, the observation of hidden-charm pentaquark states [318]; and the trace anomaly is topical because of its connection with EHM.
Lacking vector meson beams, ongoing and anticipated experiments at electron (e) accelerators are based on an expectation that the desired V + p interactions can be accessed through the electromagnetic production of vector mesons from the proton, in reactions like e + p e + V + p [183,252,319]. This is because some practitioners imagine that single-pole vector meson dominance (VMD) [320,321,322] can reliably be employed to draw a clean link between e + p e + V + p and the desired V p V p cross-sections. In this picture (see Figure 12), the interaction is supposed to proceed via the following sequence of steps: (i) e e + γ ( * ) ( Q ) ; (ii) γ ( * ) ( Q ) V ; and (iii) V + p V + p . γ ( * ) ( Q ) is a virtual photon; and step (ii) expresses the VMD hypothesis. As commonly used, VMD assumes: (a) that a photon, which is generally spacelike, so that Q 2 > 0 , transforms into an on-shell vector meson, with timelike momentum Q 2 = m V 2 ; and (b) that the Q 2 > 0 strength and form of the transition in (ii) is the same as that measured in the real vector meson decay process, V γ * ( Q 2 = m V 2 ) e + + e . Property (b) means that γ γ V , the associated decay constant, is fixed at its meson on-shell value and acquires no momentum dependence:
γ γ V 2 = 4 π α em m V 2 f V 2 ,
where α em is the QED fine-structure constant, m V is the vector meson mass, and f V measures the strength of the meson’s Bethe–Salpeter wave function at the origin in configuration space [324] (Section IIB).
The fidelity of these VMD assumptions was recently subjected to scrutiny via analyses of the photon vacuum polarisation and photon–quark vertex [204]. Regarding the photon vacuum polarisation, it was shown that there is no vector meson contribution to this polarisation at the photoproduction point, Q 2 = 0 . Consequently, massless real photons cannot readily be linked with massive vector bosons, and the current field identity, Equation (28), typical of VMD implementations, should not be used literally because it entails violations of Ward–Green–Takahashi identities (breaking of symmetries) in QED.
Turning to the dressed photon–quark vertex, Γ ν γ ( k ; Q ) , this three-point Schwinger function exhibits a pole at the mass of any vector meson bound state. A physical property, it expresses the fact that the decay V e + e proceeds via a timelike virtual photon. The VMD hypothesis may thus be seen as a claim that Γ ν γ ( k ; Q ) Q 2 0 maintains a rigorous link, in both magnitude and momentum dependence, with the Bethe–Salpeter amplitude of an on-shell vector meson. However, as shown in Figure 13, that is not the case. The panels in Figure 13 depict the ratios:
R V ( k 2 ; Q 2 = 0 ) = 2 f V m V F 1 0 ( k 2 ; m V 2 ) G 1 0 ( k 2 ; Q 2 = 0 ) , R V B ( k 2 ; Q 2 = 0 ) = f V m V F 8 0 ( k 2 ; m V 2 ) G 3 0 ( k 2 ; Q 2 = 0 ) ,
where G 1 0 is the zeroth Chebyshev moment5 of the dominant amplitude in the photon–quark vertex, associated with the matrix structure γ · ϵ ( Q ) , where ϵ ( Q ) is the photon polarisation vector—in fact, using the vector Ward–Green–Takahashi identity and Equation (23), G 1 0 ( k 2 ; Q 2 = 0 ) = A g ( k 2 ) , where g is the flavour of the meson’s valence quark; F 1 0 is its analogue in the vector meson bound state amplitude; G 3 0 is the zeroth moment of that term in the photon–quark vertex, which is directly linked to the scalar piece of the dressed quark self-energy via the vector Ward–Green–Takahashi identity, i.e., G 3 0 ( k 2 ; Q 2 = 0 ) = 2 B g ( k 2 ) ; and F 8 0 is its analogue in the vector meson bound state amplitude. Were VMD to be a sound assumption, then all these curves would lie near the thin horizontal line drawn at unity in both panels of Figure 13. However, whilst one might discuss the case for lighter vector mesons, the VMD hypothesis is plainly false for heavy vector mesons: the momentum-dependence of the Q 2 = 0 photon–quark vertex is completely different from that of the vector meson Bethe–Salpeter amplitude.
One is consequently led to conclude that, insofar as heavy mesons are concerned, no extant attempt to link e + p e + V + p reactions with V + p V + p via VMD is reliable. A similar conclusion is drawn in Reference [323] using arguments within perturbative QCD. It is strengthened by Reference [325], which demonstrates that even if VMD were valid, then contributing relevant coupled channels processes would obscure connections between e + p e + V + p and V + p V + p reactions. These analyses make tenuous any interpretation of e + p e + V + p reactions as a path to hidden-charm pentaquark production or as a route to understanding the origin of the proton mass. In fact, the crudity of existing models prevents sound conclusions being drawn about the capacity of heavy meson production to reveal something about the SM [326].
On the other hand, numerous applications [28,31,32,33,327], including γ * γ π 0 , η , η , η c , and η b [201,308], show that a viable alternative to the VMD hypothesis exists in adapting the CSMs discussed herein to a direct analysis of processes like γ ( * ) + p V + p . Considering vector meson photoproduction/electroproduction from the proton, References [326,328,329] illustrate how one might proceed; and given developments in the past vicennium, improvements of such studies are now possible.
Following the first studies almost thirty years ago [330], the CSM treatment of meson elastic and transition form factors has now reached a mature level. Today, as sketched above, sound predictions with a traceable connection to QCD are being delivered. This enables new opportunities to be exploited, such as informative comparison with results from lQCD [200], weak transitions of heavy mesons (see References [205,206,207] and Section 10), calculations of form factors describing light-by-light scattering contributions to the muon anomalous magnetic moment [331], and a beginning to the analysis of meson gravitational form factors [332].

9. Baryon Form Factors

Advances in the calculation of meson form factors are complemented by progress with the parameter-free prediction of baryon elastic and transition form factors. The first direct RL truncation study of nucleon elastic form factors was presented in Reference [333]. It was extended to nucleon axial and pseudoscalar form factors in Reference [334], the elastic form factors of Δ and Ω baryons in Reference [335], and nucleon tensor charges in Reference [31].
Many other direct RL truncation analyses are in train. However, in all such studies, the challenge of quark propagator singularities moving into the complex integration domain must also be overcome [199]. Using existing algorithms, the singularities limit the Q 2 reach of baryon form factor calculations. Here, again, the QCD-kindred quark plus fully interacting, nonpointlike diquark picture of baryon structure—outlined in Section 7—can profitably be exploited. Some of the successes are summarised elsewhere [181,255]. Recently, the predictions made for γ + p Δ ( 1600 ) transition form factors [273] have been tested in analyses of π + π p electroproduction data collected at Jefferson Lab (JLab), with preliminary results supporting the quark+diquark picture [336].
Scanning the sources indicated above, it becomes apparent that, hitherto, nucleon properties have largely been probed in e + N scattering. Aspects of this field are reviewed elsewhere [337]. An entirely new window onto baryon structure is opened when one uses neutrino scattering. In fact, reliable predictions of nucleon and N Δ ( 1232 ) electroweak form factors are crucial for understanding new-generation long-baseline neutrino oscillation experiments [338,339,340,341,342,343,344]. Importantly, in this connection, recent developments within the framework of CSMs have enabled practitioners to deliver the first Poincaré-invariant parameter-free predictions for such form factors on a momentum transfer domain that extends to Q 2 = 10 GeV2 [274,275,276]. Extensions to even larger Q 2 are feasible. Where data are available, the predictions confirm the measurements. More significantly, the results are serving as motivation for new experiments at high-luminosity facilities.
The key step in the use of CSMs was the construction of a symmetry-preserving current that describes the coupling between axial vector and pseudoscalar probes and baryons whose structure is determined by the Faddeev equation in Figure 8. This current is illustrated in Figure 14 and explained by the legend in Table 2. The origins and characters of Diagrams 1–3 are obvious: the probe must interact with every “constituent” that carries a weak charge. Diagram 4 is a two-loop diagram, made necessary by the quark exchange nature of the kernel in Figure 8: the object exchanged in binding also carries a weak charge. Given the presence of Diagram 4, so-called seagull diagrams—Diagrams 5 and 6—are necessary to ensure symmetry preservation. The analogous contributions for baryon electromagnetic currents were derived in Reference [345], but it took more than twenty years before the seagull terms were derived for axial vector and pseudoscalar currents [275]. These contributions are both two-loop diagrams.
Axial vector interactions of a nucleon are described by two form factors—- G A ( Q 2 ) (axial) and G P ( Q 2 ) (induced pseudoscalar)—associated with the following matrix element:
J ^ 5 μ j ( K , Q ) : = N ( P f ) | A 5 μ j ( 0 ) | N ( P i )
= u ¯ ( P f ) τ j 2 γ 5 γ μ G A ( Q 2 ) + i Q μ 2 m N G P ( Q 2 ) u ( P i ) ,
where P i , f are, respectively, the initial and final nucleon momenta, with P i , f 2 = m N 2 , m N is the nucleon mass, and u ( P ) is the associated Euclidean spinor. (Associated conventions are specified elsewhere, e.g., Reference [271] (Appendix B).) The average momentum of the system is K = ( P i + P f ) / 2 , and Q = P f P i is the momentum transferred between the initial and final states. It is usual to consider the SU ( 2 ) F isospin limit m u = m d = : m q , with the flavour structure described using Pauli matrices { τ j | j = 1 , 2 , 3 } : τ 1 ± i 2 : = ( τ 1 ± i τ 2 ) / 2 correspond to the weak charged currents and τ 3 is the neutral current. Moreover, the isovector axial current operator is
A 5 μ j ( x ) = q ¯ ( x ) τ j 2 γ 5 γ μ q ( x ) , q = u d .
A third form factor is defined via the kindred pseudoscalar current, a matrix element of P 5 j ( x ) = q ¯ ( x ) τ j 2 γ 5 q ( x ) :
J ^ 5 j ( K , Q ) : = N ( P f ) | P 5 j ( 0 ) | N ( P i ) = u ¯ ( P f ) τ j 2 γ 5 G 5 ( Q 2 ) u ( P i ) .
Using the framework deployed to calculate the baryon wave functions discussed in Section 7, References [274,275] delivered the parameter-free prediction for G A ( Q 2 ) displayed in Figure 15: the lighter blue band expresses the impact of ± 5 % variations in the diquark masses: m 0 + = 0.80 ( 1 ± 0.05 ) ; m 1 + = 0.89 ( 1 ± 0.05 ) . The calculated values for the nucleon axial charge g A = G A ( Q 2 = 0 ) , associated axial charge radius, and axial mass are, respectively,
g A = 1.25 ( 3 ) , r A 2 1 / 2 m N = 3.25 ( 4 ) , m A / m N = 1.23 ( 3 ) .
For comparison, empirically, g A = 1.2756 ( 13 ) [40] (PDG) and [346] r A 2 1 / 2 m N = 3.23 ( 72 ) , m A / m N = 1.15 ( 08 ) . Evidently, the CSM predictions agree with extant data. Regarding the axial mass, we note that it is sometimes convenient, when comparing with other analyses, to use a dipole Ansatz as an approximation for the axial form factor:
G A ( Q 2 ) = g A 1 + Q 2 / m A 2 2 ;
therefore References [274,275] extracted m A using Equation (34) to interpolate the global Q 2 behaviour of G A on x [ 0 , 1.6 ] , in which case m A is not simply related to the axial radius. It is worth remarking that scalar and axial vector diquark mass variations interfere destructively, e.g., reducing m 0 + increases g A , whereas g A decreases with the same change in the axial vector mass.
Turning to the induced pseudoscalar form factor, muon capture experiments ( μ + p ν μ + n ) may be used to determine the induced pseudoscalar charge:
g p * = m μ 2 m N G P ( Q 2 = 0.88 m μ 2 ) .
The CSM prediction is g p * = 8.80 ( 23 ) . Compared with the recent MuCap Collaboration result, g p * = 8.06 ( 55 ) [348,349], it agrees within uncertainties, but is slightly larger. The CSM value is nicely aligned with the world average [350]: g p * = 8.79 ( 1.92 ) .
The pseudoscalar form factor, G 5 ( Q 2 ) , is of interest because, inter alia, it is used to define the pion–nucleon form factor G π N N ( Q 2 ) via
G 5 ( Q 2 ) = : m π 2 Q 2 + m π 2 f π m q G π N N ( Q 2 ) ,
where f π is the pion leptonic decay constant (appearing in Figure 7) and G π N N ( m π 2 ) = g π N N is the π N N coupling, a key input to nucleon+nucleon potentials. The CSM prediction is [274,275]: g π N N / m N = 14.02 ( 33 ) / GeV . This value overlaps with that inferred from pion–nucleon scattering [351] ( g π N N / m N = 13.97 ( 10 ) / GeV ) and compares favourably with a determination based on the Granada 2013 n p and p p scattering database [352] ( g π N N / m N = 14.11 ( 3 ) / GeV ) and a recent analysis of nucleon–nucleon scattering using effective field theory and related tools [353] ( g π N N / m N = 14.09 ( 4 ) / GeV ). All these results are compared in Figure 16, which also highlights their error-weighted average:
g π N N / m N = 14.10 ( 2 ) / GeV .
Continuing with an effort to inform nuclear physics potentials using hadron physics results, it is worth noting that on m π 2 < Q 2 < 2 m N 2 , a fair approximation to the CSM prediction for the pion–nucleon form factor is provided by ( x = Q 2 / m N 2 ) :
G π N N d ( x ) = 13.47 m N ( 1 + x / 0 . 845 2 ) 2 .
This corresponds to a π N N dipole mass Λ π N N = 0.845 m N = 0.79 GeV, viz. a soft form factor. (A commensurate value was obtained previously in a simpler quark+scalar-diquark model [354].) Being just ~ 20 % greater, the CSM prediction is qualitatively equivalent to the π N N dipole mass inferred in a dynamical coupled channels analysis of π N , γ N interactions [292]. Future such coupled channels studies may profit by implementing couplings and range parameters determined in CSM analyses.
An extension of this effort to the more general case of hyperon+nucleon potentials can be found in Reference [355]. Using an SCI, predictions were made therein for an array of meson+octet-baryon couplings. Comparing the results with extant phenomenological interactions [292,356,357], one finds a mean absolute relative deviation (P is the pseudoscalar meson absorbed in the baryon transition B B , and g P B B is the associated coupling):
δ r g : = { | g P B B SCI / g P B B phen . 1 | } = 0.18 ( 14 ) .
This result strengthens the case in favour of using CSM predictions for the couplings as new constraints in the development of baryon+baryon potentials.
A little algebra reveals that the proton axial charge can be interpreted as a measure of the valence quark contributions to the proton light-front helicity, e.g., Reference [358] (Equations (6)–(8)):
g A = 0 1 d x [ Δ u p ( x ; ζ H ) Δ d p ( x ; ζ H ) ] = : g A u g A d ,
where Δ q p ( x ; ζ H ) = q p ( x ; ζ H ) q p ( x ; ζ H ) is the light-front helicity DF for a quark q. Plainly, Δ q p is the difference between the light-front number density of quarks with helicity parallel and antiparallel to that of the proton. It is scale-dependent.
Equation (40) conveys additional significance to a flavour separation of the axial charge form factor:
G A ( Q 2 ) = G A u ( Q 2 ) G A d ( Q 2 ) .
A detailed analysis is presented in Reference [276] (Section 4), which reveals the following diagram contributions to the separate u, d axial form factors:
G A u = J q S J q A ̲ + J qq A A + 1 2 J qq { S A } + 2 J ex { S A } + 4 5 J ex A A ,
G A d = J q S 2 J q A ̲ + 1 2 J qq { S A } + J ex S S J ex { S A } + 1 5 J ex A A ,
where the nomenclature of Table 2 is used. Identified according to Equations (42), the calculated Q 2 = 0 contributions are listed in Table 3. It is worth stressing that Equations (42) express the fact that since a 0 + diquark cannot couple to an axial vector current, then Diagram 1 in Figure 14 only supplies a u quark contribution to the proton G A ( Q 2 ) , viz. J q S . It follows that in a scalar-diquark-only proton, a d quark contribution to G A ( Q 2 ) can only arise from Figure 14, Diagram 4, i.e., J ex S S ; and | J ex S S / J q S | 0.06 . Notably, many scalar-diquark-only models omit Diagram 4, in which case G A d ( Q 2 ) 0 . An extension of these observations to the complete array of octet baryons is described elsewhere [355] (Section IV).
Using the solution of the Faddeev equation (Figure 8), Reference [276] reports
g A u / g A = 0.76 ± 0.01 , g A d / g A = 0.24 ± 0.01 , g A d / g A u = 0.32 ± 0.02 .
In nonrelativistic quark models with uncorrelated wave functions, g A d / g A u = 1 / 4 . Hence, the relevant comparison reveals that the highly correlated wave function obtained by solving the Faddeev equation gives the valence d quark a markedly larger fraction of the proton’s light-front helicity than is found in simple quark models. Reviewing the discussion after Equation (42), it becomes apparent that this feature owes to the presence of axial vector diquarks in the proton: the current contribution arising from the { u u } correlation—underlined term in Equation (42b)—measuring the probe striking the valence d quark, is twice as strong as that from the { u d } correlation—underlined term in Equation (42a)—in which the probe strikes the valence u quark.
It is natural to enquire after the robustness of the results in Equation (43). Consider, therefore, that assuming SU ( 3 ) flavour symmetry in analyses of octet baryon axial charges, these charges are expressed in terms of two low-energy constants ([359] (Table 1)): D, F, with g A u = 2 F , g A d = F D . (This assumption is accurate to roughly 4%—see, e.g., Reference [355].) In this case, the values in Equation (43) correspond to
D = 0.78 ( 2 ) , F = 0.48 ( 1 ) , F / D = 0.61 ( 2 ) .
On the other hand, using available empirical information [40], one obtains D = 0.774 ( 26 ) , F = 0.503 ( 27 ) , and g A u / g A = 0.79 ( 4 ) , g A d / g A = 0.21 ( 3 ) , and g A d / g A u = 0.27 ( 4 ) , values that are consistent with the results in Equations (43) and (44).6 In addition, the SCI predicts [355]
D = 0.78 , F = 0.43 , F / D = 0.56 ,
and a covariant baryon chiral perturbation theory analysis yields D = 0.80 ( 1 ) , F = 0.47 ( 1 ) , F / D = 0.59 ( 1 ) [360].
Given the favourable realistic proton wave function comparisons presented above, the values in Equation (43) can be viewed as reliable. This is important because of the connection between flavour-separated axial charges and the so-called proton “spin crisis” [361,362]. At any given resolving scale, the singlet, triplet, and octet axial charges of the proton are, respectively:
a 0 = g A u + g A d + g A s , a 3 = g A = g A u g A d , a 8 = g A u + g A d 2 g A s .
If working at the hadron scale, ζ H , where dressed valence quasiparticles carry all proton properties [126,127,132,133,134,135], then g A s 0 a 0 = a 8 ; hence [276],
a 0 = 0.65 ( 2 ) .
In general, a 3 , 8 are conserved charges, viz. they are the same at all resolving scales, ζ . However, that is not true of the individual terms in their definitions: the separate valence quark charges g A u , g A d , g A s evolve with ζ [362]. Consequently, the value of a 0 , which is the fraction of the proton’s total J = 1 / 2 carried by valence degrees-of-freedom, changes with scale: the result in Equation (47) is the maximum value of a 0 , and the fraction falls slowly with increasing ζ .
Returning to expectations based on simple, nonrelativistic quark models, textbook-level algebra yields a 0 = 1 . Therefore, in such pictures, all the proton spin derives from that of the constituent quarks. On the other hand, the CSM analysis in Reference [276] predicts that proton dressed valence degrees-of-freedom carry only two-thirds of the spin. Since there are no other degrees-of-freedom at ζ H and the Poincaré-covariant proton wave function properly describes a J = 1 / 2 system, then the “missing” part of the total J must be associated with quark+diquark orbital angular momentum. Similar conclusions apply for all ground state octet baryons [355].
The study in Reference [276] delivered Poincaré-invariant parameter-free predictions for the proton axial form factor and its flavour separation out to Q 2 10 m N 2 . The axial form factor itself agrees with available data [363,364], which extends to Q 2 5 m N 2 —see Reference [276] (Figure 3). More importantly, the results will likely serve as motivation for new experiments aimed at exploring nucleon structure with axial vector probes instead of the photon, opening a new window onto hadron structure. It is worth highlighting here that a dipole fit to data is only a good approximation on the fitting domain. With increasing Q 2 , the dipole increasingly overestimates the actual result, being 56 ( 5 ) % too large at Q 2 = 10 m N 2 —see Reference [276] (Figure 4B). It therefore becomes an unsound tool for developing qualitative insights and quantitative cross-section estimates.
Furthermore, with flavour-separated form factors in hand on such a large- Q 2 -domain, Reference [276] was able to calculate and contrast the u and d quark contributions to the associated light-front transverse spatial density profiles:
ρ ^ A f ( | b ^ | ) = d 2 q ( 2 π ) 2 e i q · b ^ G A f ( x ) ,
with G A f ( x ) interpreted in a frame defined by Q 2 = m N 2 q 2 , m N q = ( q 1 , q 2 , 0 , 0 ) = ( Q 1 , Q 2 , 0 , 0 ) . These profiles are depicted in Figure 17. We note that | b ^ | and ρ ^ A f are dimensionless; so, the images drawn in Figure 17 can be mapped into physical units using:
ρ A f ( | b | = | b ^ | / m N ) = m N 2 ρ ^ A f ( | b ^ | ) ,
in which case | b ^ | = 1 corresponds to | b | 0.2 fm and ρ ^ 1 f = 0.1 ρ 1 f 2.3 / fm 2 .
The top row of Figure 17 provides two-dimensional renderings of the flavour-separated transverse density profiles calculated from a proton wave function which does not include axial vector diquarks, i.e., a scalar-diquark-only proton. In this case, the u quark profile is far more diffuse than that of the d quark, viz. its extent in the light-front transverse spatial plane is much greater. One may quantify this by reporting the associated radii: r A d = 0.24 fm, r A u = 0.48 fm, so the d / u ratio of radii is ≈ 0.5 .
The bottom row of Figure 17 was obtained using a realistic proton wave function, in which both scalar and axial vector diquarks are present with the strength determined by the Faddeev equation in Figure 8. In this realistic case, the d quark profile is not very different from that of the u quark: relative to the u quark profile, the intensity peak is only somewhat broader for the d quark; and comparing radii,
r A d = 0.43 fm , r A u = 0.49 fm , r A d / r A u 0.9 .
The CSM predictions for nucleon axial and pseudoscalar form factors discussed in this section complement those for the large- Q 2 behaviour of nucleon electromagnetic elastic and transition form factors reported, e.g., in References [272,365]. One may now anticipate that predictions for form factors characterising weak interaction induced N Δ ( 1232 ) and N N * ( 1535 ) transitions will soon become available. Each will shed new light on nucleon structure; and the former, calculated on a domain that stretches from low-to-large- Q 2 , will likely prove valuable in developing a more reliable understanding of neutrino scattering from nucleons and nuclei.

10. Transition Form Factors of Heavy+Light Mesons

Heavy mesons ( B q = c , s , u , D s , D) are special for many reasons; and their mass budgets and role in exposing constructive interference between Nature’s two known sources of mass are of particular interest herein. Consider, therefore, Figure 18 and contrast the images with those in Figure 1. Evidently, for heavy meson masses: (i) the HB component is largest in each case and its relative size grows as the current–masses of the valence constituents increase; (ii) all receive a significant EHM+HB interference component, but its relative strength diminishes with increasing current–masses; and (iii) for vector heavy mesons, but not pseudoscalar mesons, there is an EHM component, but its relative strength drops as the HB component grows.
Next consider semileptonic weak-interaction transitions between heavy and light mesons. Comparing Figure 1 and Figure 18, it is apparent that heavy-pseudoscalar to light-pseudoscalar transitions serve to probe the relative impacts of the strength of EHM+HB interference in the initial and final states, whereas heavy-pseudoscalar to light-vector transitions overlap systems in which HB mass is dominant with those whose mass is owed almost entirely to EHM. Both classes of transitions, therefore, present excellent opportunities for exposing the influence of Nature’s two known sources of mass on physical observables. These cases are of special interest, of course, because the transitions have long been used to place constraints on the values of the elements of the Cabibbo–Kobayashi–Maskawa (CKM) matrix, which parametrises quark flavour mixing in the SM.
A unifying analysis of both classes of transitions was recently completed using the SCI [233,366], yielding results that compare favourably with other reliable experimental or independent theory analyses. The SCI branching fraction predictions should therefore be a reasonable guide. This is important because predictions were made for several branching fraction ratios— R D ( s ) ( * ) , R J / ψ , R η c —whose values are direct tests of lepton universality in weak interactions: the SCI values confirm SM predictions, hence, speak for universality, as we discuss below. The analyses also used B ( s ) D ( s ) * transitions to predict the precursor functions which evolve into the universal Isgur–Wise function [367], obtaining results in agreement with empirical inferences—References [368] (Equations (177) and (181)) and [369] (Belle).
The SCI’s successes in these applications highlight the need for kindred studies using an interaction with a closer connection to QCD. The impediment has always been the large disparity in masses that typically exists between initial and final states. That mass imbalance requires, inter alia, that any approach to the problem be simultaneously able to deal with both chiral and heavy quark limits in quantum field theory. To date, compared with the SCI, no framework with a better link to QCD can directly surmount this difficulty. Nevertheless, following Reference [205], practicable and effective algorithms for continuum studies do now exist. They exploit the strengths of the statistical SPM as a tool for interpolating data (broadly defined) and therefrom deliver extrapolations with a rigorously defined and calculable uncertainty [238]. Namely, results are calculated on domains of current–quark mass for which transition form factors may straightforwardly be obtained. The SPM is then used to extrapolate those results and arrive at predictions for the physical processes of interest.
At present, CSM RL truncation predictions are available for the following transitions [206,207]: B c J / ψ , η c ; B ( s ) π ( K ) ; D s K ; D π , K ; and K π . The last process is something of a test for the approach because such K 3 transitions have long been of experimental and theoretical interest [40] (Section 62). The calculated branching fractions are gathered in Table 4, from which it will be seen that CSMs deliver sound results.
One of the key results in Reference [206] concerns a ratio of B c + J / ψ branching fractions measured for the first time by the LHCb Collaboration fairly recently [371]:
R J / ψ : = B B c + J / ψ τ ν B B c + J / ψ μ ν = 0.71 ± 0.17 ( stat ) ± 0.18 ( syst ) .
This value is plotted in Figure 19 and compared with the CSM prediction and other Standard Model calculations. Evidently, the LHCb measurement lies approximately 2 σ above the values predicted by reliable SM calculations. If future, precision experiments do not deliver a markedly lower central value, then one might begin to judge that lepton flavour universality is violated in B c J / ψ semileptonic decays. As yet, however, the experimental precision is insufficient to support such a claim. Furthermore, a compelling case would need to include information on B c η c semileptonic decays. The CSM prediction is R η c = 0.313 ( 22 ) —see Table 4; the SCI result is R η c = 0.25 [233]; and a mean value of 0.31 ( 4 ) is obtained from modern continuum analyses [372,373,374,375,376,377]. An experimental value is lacking.
The array of analyses in Reference [207] yields novel results in other areas. Of particular interest are the discussions of D 0 K transition form factors and the value of | V c s | . Two independent form factors characterise 0 + 0 + transitions, viz. vector and scalar, f + , 0 ( t ) , respectively, where t is the Mandelstam variable, whose value expresses the momentum transferred to the final state. The CSM predictions are plotted in Figure 20 and compared with available data [380,381,382,383]. The CSM result is largely consistent with this collection, although there may be a hint that it is too high at lower t values. Concerning branching fractions, form factor contributions from this domain are important. It is therefore notable that, within mutual uncertainties, the CSM value for f + D K ( 0 ) = 0.796 ( 9 ) agrees with the N f = 2 + 1 + 1 lQCD result in Reference [384]: 0.765 ( 31 ) .
With the form factors in Figure 20, Reference [207] obtained the D 0 K branching fractions listed in Table 4 when using the value of | V c s | listed in the caption: evidently, both the e + ν e and μ + ν μ fractions exceed their respective PDG values. On the other hand, the ratio agrees within 1.4 σ ; so, a common overall factor can remedy the mismatch. Adopting this perspective, then the value | V c s | = 0.937 ( 17 ) combined with the CSM form factors delivers branching fractions that match the PDG values, viz. 3.52 ( 18 ) % and 3.44 ( 18 ) %, respectively. Actually, referring to Reference [40] (Section 12.2.4), one sees that the inferred CSM value is both commensurate with and more precise than one of the two used to arrive at the PDG average listed in the caption of Table 4. With | V c s | = 0.937 ( 17 ) used instead to compute this average, one finds a slightly more precise central value that is 1 σ lower:
| V c s | = 0.974 ( 10 ) .
Predictions for semileptonic B ¯ 0 π + , B ¯ s 0 K + transition form factors and branching fractions were also delivered in Reference [207]. As emphasised above, such processes present challenges because π , K are Nambu–Goldstone bosons and there is a huge disparity between the masses of the initial and final states. Consequently, comparisons with data serve as a stringent test of the new CSM algorithms.
CSM predictions for the B ¯ 0 π + transition form factors are depicted in Figure 21A. Regarding f + B ¯ π , data have been collected by two collaborations [385,386,387,388]: within mutual uncertainties, the CSM predictions agree with these data. The data support a value
f + B ¯ π ( t = 0 ) = 0.27 ( 2 ) ,
which is consistent with the CSM prediction [207]: 0.29 ( 5 ) .
Using the form factors in Figure 21A, one obtains the B ¯ 0 π + branching fractions in Table 4. The PDG lists a result for the μ ν μ final state, which matches the CSM prediction within mutual uncertainties. Precise agreement is obtained using
| V u b | = 0.00374 ( 44 ) .
This value is commensurate with other analyses of B B ¯ 0 π + μ ν ¯ μ [40] (Section 76.3), thus increases tension with the higher value inferred from inclusive decays. No data are available on the τ ν τ final state; so, the τ : μ ratio is empirically unknown. Here, a N f = 2 + 1 -flavour lQCD study yields 0.69 ( 19 ) [390], which, within its uncertainty, matches the CSM result: 0.733 ( 2 ) .
Figure 21B displays CSM predictions for the B ¯ s K + form factors. Although the B s K transition was recently observed [391], with the measurement yielding the branching fraction:
B B s 0 K μ + ν μ = [ 0.106 ± 0 . 005 stat ± 0 . 008 syst ] × 10 3 ,
no form factor data are yet available. Comparisons are therefore made in Figure 21B with results obtained using N f = 2 + 1 -flavour lQCD [389,390]. Owing to difficulties encountered when using lattice methods to calculate form factors of heavy+light mesons, lQCD results are limited to a few points on the domain t 17 GeV2—see Figure 21B. Today, lattice analyses typically employ such results to construct a least-squares fit to the form factor points, using some practitioner-favoured functional form. That fit is then employed to define the form factor on the whole kinematically accessible domain: 0 t 25 GeV2 in this case. It is worth noting that, at this time, given the small number of points and their limited precision, the SPM cannot gainfully be used to develop function–form unbiased interpolations and extrapolations of the lQCD output.
The CSM form factors in Figure 21B yield the B ¯ s 0 K + branching fractions in Table 4. Figure 22A compares the μ ν ¯ μ value with the measurement in Equation (55) and also results obtained via various other means. Experiment and theory only agree because the theory uncertainty is large. The unweighted theory average is 0.141 ( 44 ) ‰, and the uncertainty-weighted mean is 0.139 ( 08 ) ‰. These values increase when Entries V–VI [390,392] are omitted: unweighted 0.159 ( 38 ) ‰ and uncertainty weighted 0.156 ( 10 ) ‰. The extrapolations employed in V–VI [390,392] lead to values of f + B ¯ s K ( 0 ) that are 50 % of those obtained in I–IV [389,393,394,395]: 0.148 ( 53 ) vs. 0.299 ( 86 ) . This can explain the difference in branching fractions: V–VI vs. I–IV in Figure 22A. Significantly, a different approach to fitting and extrapolating lQCD results, using the LHCb datum, Equation (55), as an additional constraint, produces [396]: f + B ¯ s K ( 0 ) = 0.211 ( 3 ) .
Reference [396] also infers f + B ¯ π ( 0 ) = 0.255 ( 5 ) , leading to f + B ¯ s K ( 0 ) / f + B ¯ π ( 0 ) < 1 . This outcome conflicts with the CSM prediction, which has f + B ¯ s K ( 0 ) / f + B ¯ π ( 0 ) > 1 at the 85% confidence level [207] (Equation (9)), and the results in a raft of other studies, e.g., References [393,394,397,398,399,400,401,402,403]. It is probable, therefore, that the value for f + B ¯ s K ( 0 ) in Reference [396] is too small. It is worth remarking that the SCI is unclear on the value of this ratio. It produces f + B ¯ s K ( 0 ) / f + B ¯ π ( 0 ) < 1 , but the individual t = 0 values are too large by a factor of two [233] (Table 3A). On the other hand, the t = 0 value of the kindred ratio of vector form factors in B ¯ s K * , B ¯ ρ transitions is greater than unity [366] (Table 1A). Notably, the ratio f + B ¯ s K ( 0 ) / f + B ¯ π ( 0 ) is a marker for SU ( 3 ) -flavour symmetry breaking and its modulation by EHM; so, it is worth reaching a sound conclusion on the value of the ratio. It is here relevant to observe that f K / f π = 1.2 > 1 .
Regarding the τ ν τ final state in B ¯ s K transitions, no empirical information is currently available, hence none on the | V u b | -independent ratio that would test lepton flavour universality. In Figure 22B, therefore, we compare the CSM prediction for this ratio, drawn from Table 4, with results obtained via other means. The unweighted average of theory results is 0.705 ( 87 ) , and the uncertainty weighted mean is 0.678 ( 03 ) . Omitting Entries V–VI [390,392], these values are: unweighted 0.653 ( 41 ) and uncertainty weighted 0.677 ( 03 ) . Within sound analyses, many uncertainties cancel in this ratio; so, the results should be more reliable than any calculation of either fraction alone. Nevertheless, the values are widely scattered, indicating that there is ample room for improving the precision of B ¯ s 0 K theory. Of course, measurements enabling extraction of B s 0 K form factors would be very useful, too, for refining both (a) comparisons with theory and between theory analyses and (b) making progress toward a more accurate value of | V u b | .
Such predictions for heavy-to-light meson electroweak transition form factors are a new branch of application for CSMs. They are far more sophisticated and robust than the Schwinger function parametrisation-based analyses in, e.g., References [324,403,404] and significantly improve upon earlier RL truncation studies of π 3 and K 3 transitions [405,406]. The keys to these advances are an improved understanding of RL truncation and the capacity to greatly expand its quark mass and mass splitting domains of applicability using the SPM. Fairly soon, one can expect these advances to be exploited in the study of kindred baryon transitions.

11. Distribution Functions

Hadron parton DFs are probability densities: each one describes the light-front fraction, x, of the hadron’s total momentum carried by a given parton species within the bound state [407]. They are a much-prized source of hadron structure information; and following the quark discovery experiments fifty years ago [408,409,410,411], measurements interpretable in terms of hadron DFs have been awarded a high priority. For much of this time, DFs were inferred from global fits to data, with the results viewed as benchmarks. Such fitting remains crucial, providing input for the conduct of a huge number of experiments worldwide; but the past decade has seen the dawn of a new theory era. Continuum and lattice studies of QCD are beginning to yield robust predictions for the pointwise behaviour of DFs; and these developments are exposing potential conflicts with the fitting results [35,129,132,133,134,135,136].
Notwithstanding the enormous expense of time and effort, much yet remains to be learnt before proton and pion structure may be judged as understood in terms of DFs. For instance and most simply, it is still unclear whether there are differences between the distributions of partons within the proton (Nature’s most fundamental bound state) and the pion (Nature’s most fundamental (near) Nambu–Goldstone boson). Plainly, if there are differences, then they must be explained. As we have stressed above, answering the question of similarity/difference between proton and pion DFs is particularly important today as science seeks to expose and explain EHM [32,33,34,35,36,37].
Regarding DFs measured in processes that do not resolve beam or target polarisation, practitioners experienced and involved with solving bound state problems in QCD have learnt that, at the hadron scale, ζ H < m p , valence quark DFs in the proton and pion behave as follows [132,133,134,412,413]:
d p ( x ; ζ H ) , u p ( x ; ζ H ) x 1 ( 1 x ) 3 , d ¯ π ( x ; ζ H ) , u π ( x ; ζ H ) x 1 ( 1 x ) 2 .
It subsequently follows from the DGLAP equations [177,414,415,416] that the large-x power on the related gluon DF is approximately one unit larger; and that for sea quark DFs is roughly two units larger. Moreover, as the resolving scale increases to ζ > ζ H , all these exponents grow logarithmically. However, fuelling controversy and leading some to question the veracity of QCD [132,133,417], these constraints are typically ignored in fits to the world’s data on deep inelastic scattering (DIS) and kindred processes [418,419,420,421,422]. Furthermore, largely because pion data are scarce [35] (Table 9.5), proton and pion data have never been considered simultaneously. Therefore, the unified body of results in Reference [135], which uses a single symmetry-preserving framework to predict the pointwise behaviour of all proton and pion DFs—valence, glue, and four-flavour-separated sea—is a significant advance.
In order to sketch this progress, it is necessary to recall that the modern approach to the CSM prediction of hadron DFs7 is based on a single proposition [131,132,133,134,135]:
P1
There is an effective charge, α 1 ( k 2 ) , which, when used to integrate the one-loop perturbative-QCD DGLAP equations, defines a DF evolution scheme that is all orders exact.
As noted in connection with Figure 3, charges of this sort are discussed in References [137,138,139]. They need not be process-independent (PI), hence not unique. Moreover, the results delivered are independent of the explicit form of α 1 ( k 2 ) . Notwithstanding these things, a suitable PI charge is available, viz. the coupling discussed in Section 4, which has proven efficacious. In being defined by an observable—in this instance, structure functions—each such α 1 ( k 2 ) is [94]: consistent with the renormalisation group and renormalisation scheme independent; everywhere analytic and finite; and, crucially, provides an infrared completion of any standard perturbative running coupling.
P1 was used in References [126,127,128,131,136] to deliver meson DFs with a flavour-symmetric sea. A generalisation, which expresses key quark current–mass effects in the evolution kernels, was introduced in Reference [135] and used for the proton and pion. It features a threshold function P q g ζ θ ( ζ δ f ) , which ensures that a given quark flavour only becomes active in DF evolution when the energy scale exceeds a value determined by the quark’s mass [35] (Figure 2.5): δ u , d 0 , δ s 0.1 GeV, δ c 0.9 GeV. The impact of this modification is readily anticipated. Supposing that all quark flavours are light, then each would be emitted with equal probability on ζ > ζ H ; so, evolution would produce a certain gluon momentum fraction in the hadron plus a sea quark fraction shared equally between all quark flavours. Considering mass differences between the quarks, with some flavours being heavier than the light quark threshold, then evolution on ζ > ζ H will generate a gluon momentum fraction that is practically unchanged from the all-light quark case and a sea quark fraction divided amongst the quarks in roughly inverse proportion to their mass.
It is worth reiterating here that ζ H is the scale at which the valence quasiparticle degrees-of-freedom carry all properties of a given hadron [124,125,126,127,128,129,130,131,132,133,134,135,136]. Moreover, the value of this scale is a prediction. Using the PI charge discussed in Section 4 to construct bound state kernels informed by References [68,212], then
ζ H = 0.331 ( 2 ) GeV .
The value in Equation (57) is the same for all hadrons.
Furthermore, combined with evolution according to P1, the character of ζ H ensures that all hadron DFs are intertwined at every scale ζ . Hence, this perspective suggests that it is incorrect to choose independent, uncorrelated functions to parametrise the DFs of different parton species when fitting data at any scale ζ > ζ H . If one nevertheless chooses to ignore the innate associations, then DFs with unphysical features may be obtained—see, e.g., Reference [132] [Figure 6].
CSM predictions for the ζ = ζ H proton and pion valence DFs are drawn in Figure 23A. The following points are significant.
(i) Each DF is consistent with the relevant large-x scaling law in Equation (56). Hence, from the outset, whilst the ζ = ζ H momentum sum rules for each hadron are necessarily saturated by valence degrees-of-freedom, viz.
x u p ζ H = 0.687 , x d p ζ H = 0.313 , x u π ζ H = 0.5 ,
The proton and pion valence DFs nevertheless have markedly different x-dependence. (Nature’s approximate G -parity symmetry [424] entails d ¯ π ( x ; ζ ) = u π ( x ; ζ ) .) (ii) Owing to DCSB [14,76,279,280,281,282,283], an important corollary of EHM, QCD dynamics simultaneously produce a dressed light quark mass function, M u , d ( k 2 ) , that is large at infrared momenta, M D : = M u , d ( k 2 0 ) 0.4 GeV, and an almost massless pion, m π 2 / M D 2 0.1 —see Reference [35] (Section 2). As a result, u π ( x ; ζ H ) is Nature’s most dilated hadron-scale valence DF. This is highlighted by Figure 23A and Refs. [126,127] and implicit in numerous other symmetry-preserving analyses, e.g., References [201,239,425,426].
Evolving the DFs in Figure 23A according to P1, one obtains the ζ = m J / ψ = : ζ 3 distributions in Figure 23B. Plainly, although the profiles change, the relative dilation of the DFs is preserved and is therefore a verifiable prediction of the EHM paradigm.
Given that Figure 23B depicts the first CSM predictions for proton valence quark DFs, one might question their reliability. That issue can partly be addressed through a comparison with lQCD results. The calculation of individual valence DFs using lQCD is problematic owing to difficulties in handling so-called disconnected contributions [427]. In the continuum limit, however, disconnected diagrams do not contribute to the isovector DF [ u p ( x ; ζ ) d p ( x ; ζ ) ] , so computations of this difference are available [428,429]. Both analyses use the quasidistribution approach [430], but the lattice algorithms and configurations are somewhat different. Their comparison with CSM predictions is depicted in Figure 24. The level of agreement is encouraging, especially because refinements of both continuum and lattice calculations may be anticipated. For instance, the CSM predictions were obtained using a simplified proton Faddeev amplitude, and the lattice studies must address issues with, inter alia, the pion masses used, lattice artefacts and systematic errors, and convergence of the boost expansion in the quasidistribution approach. (The last of these is a particular hindrance to lQCD extractions of DF endpoint behaviour [431].)
In typical evolution kernels, gluon splitting yields quark+antiquark pairs of all flavours with equal probability. However, it was long ago argued [432] that, because the proton contains two valence u quarks and one valence d quark, the Pauli exclusion principle should force gluon splitting to prefer d + d ¯ production over u + u ¯ . Consequently, when implementing evolution of proton singlet and glue DFs, Reference [135] followed Reference [134] and introduced a small Pauli blocking factor into the gluon splitting function. This correction preserves the baryon number, but shifts momentum into d + d ¯ from u + u ¯ , otherwise leaving the sum of sea quark momentum fractions unchanged. It vanishes with increasing ζ , in order to express the declining influence of valence quarks as the proton’s sea and glue content increases.
The resulting CSM predictions for ζ = ζ 3 proton and pion glue DFs are drawn in Figure 25A. The glue-in- π DF is directly related to the ζ = 2 GeV = : ζ 2 result discussed in Reference [129], which is drawn in Figure 25B: evidently, it agrees with a recent lQCD calculation of the glue-in- π DF [433]. Furthermore, reproducing the pattern seen with valence quark DFs in Figure 23 and Figure 25A reveals that the glue-in- π DF possesses significantly more support on the valence domain, x 0.1 , than the glue-in-p DF. Once again, this feature is a measurable expression of EHM.
The ζ = ζ 3 light quark sea DFs for the proton and pion are depicted in Figure 25C. The EHM-induced pattern is also apparent here, viz. the sea-in- π DF possesses greater support on x 0.1 than the kindred sea-in-p DFs. DFs of the heavier sea quarks are also generated via evolution, with the results drawn in Figure 25D. Interestingly, the ζ = ζ 3 s and c quark DFs are similar in size to those of the light quark sea DFs; and for these heavier quarks, as well, the pion DFs have significantly greater support on the valence domain, x 0.1 , than the related proton DFs.
An analysis of the endpoint exponents of all ζ = ζ 3 DFs is also contained in Reference [135] along with simple interpolations of each DF that can readily be used by any practitioner—Reference [135] (Table 1). It is worth reiterating the following remarks about the endpoint exponents.
(i) 
The power-laws express measurable effective exponents, obtained from separate linear fits to ln [ x p ( x ) ] on the domains 0 < x < 0.005 , 0.85 < x < 1 . (Here, p ( x ) denotes a generic DF.)
(ii) 
Within mutual uncertainties, proton and pion DFs have the same power-law behaviour on x 0 :
α valence p , π 0.22 , α glue p , π 1.6 , α sea p , π 1.5 .
(iii) 
On x 1 , the following relationships exist for and between pion and proton DF exponents:
β valence π 2.5 , β valence p β valence π + 1.6 ,
β glue p , π β valence p , π + 1.4 , β sea p , π β valence p , π + 2.4 .
(iv) 
Given (ii) and (iii), then the CSM predictions are consistent with the QCD expectations discussed in connection with Equation (56).
(v) 
Existing phenomenological fits to relevant scattering data typically arrive at DFs which are inconsistent with (ii) and (iii); hence, fail to meet many QCD-based expectations, e.g., References [419,421,422,434,435]. This point is also discussed elsewhere [132,133,436].
Owing to the Pauli blocking factor described above and as evident in Figure 25C, the DFs calculated in Reference [135] express an in-proton separation between d ¯ and u ¯ distributions. This entails a violation of the Gottfried sum rule [437,438], which has been seen in experiments [439,440,441,442,443]. Using the proton DFs in Figure 25C then, on the domain covered by the measurements in References [439,440], one obtains
0.004 0.8 d x [ d ¯ ( x ; ζ 3 ) u ¯ ( x ; ζ 3 ) ] = 0.116 ( 12 )
for the Gottfried sum rule discrepancy. This value matches that inferred from recent fits to a large sample of high-precision data ( ζ = 2 GeV) [419] (CT18), 0.110(80), and is far more precise.
The result in Equation (61) corresponds to a strength for the Pauli blocking term in the gluon splitting function that shifts just 25 % of the u quark sea momentum fraction into the d quark sea at ζ = ζ 2 . Changing the strength by ± 25 % leads to the uncertainty indicated in Equation (61). Data from the most recent experiment focused on the asymmetry of antimatter in the proton [443] (E906) are presented in Figure 26A. They may be compared with the CSM result obtained using the proton DFs in Figure 25C. Evidently, a modest Paul blocking effect in the gluon splitting function is sufficient to explain modern data.
Proton and, in principle, neutron structure functions— F 2 p , n —can be measured in the DIS of electrons from nucleons [408,409,410,411]. The ratio F 2 n ( x ) / F 2 p ( x ) is recognised as a sensitive measure of d p ( x ) / u p ( x ) on x 0.4 [407], and the latter ratio is important because it is a keen discriminator between pictures of proton structure [444,445,448]. The obstacle to an empirical result for F 2 n ( x ) / F 2 p ( x ) is the measurement of F 2 n : since isolated neutrons decay rather quickly, a suitable, effective “free neutron target” must be found. Following References [451,452], many experiments have used the deuteron. However, despite this being a weakly bound system, the representation dependence of proton–neutron interactions leads to large theory uncertainties in the extracted ratio on x 0.7 [453].
A more favourable approach is provided by DIS measurements on 3 H and 3 He. In this case, nuclear interaction effects cancel to a very large degree when extracting F 2 n ( x ) / F 2 p ( x ) from the 3 H: 3 He ratio of scattering rates [454,455]. Of course, 3 H is highly radioactive; so, careful planning and implementation are required to deliver a safe target. Recently, after years of development, all challenges were overcome, and such an experiment was completed [445]: the extracted data are drawn in Figure 26B. Importantly, within mutual uncertainties, the results from Reference [445] match those inferred from an analysis of nuclear DIS reactions, exploiting targets ranging from the deuteron to lead and accounting for the effects of short-range correlations in the nuclei [456]. This speaks in support of the reliability of the analyses in both cases.
As described in Section 7, the Faddeev equation in Figure 8 makes firm statements about proton structure. In particular, well-constrained studies predict that axial vector diquark correlations are responsible for approximately 40% of the proton’s charge—see Reference [267] (Figure 2)—e.g., this strength is confirmed in studies of nucleon axial vector and pseudoscalar currents—see Section 9. Consequently, this is the size of the axial vector diquark fraction in the nucleon Faddeev amplitudes used to calculate proton DFs in References [134,135]. Using the results therein, one may readily predict the neutron–proton structure function ratio:
F 2 n ( x ; ζ ) F 2 p ( x ; ζ ) = U ( x ; ζ ) + 4 D ( x ; ζ ) + Σ ( x ; ζ ) 4 U ( x ; ζ ) + D ( x ; ζ ) + Σ ( x ; ζ ) ,
where, in terms of quark and antiquark DFs, U ( x ; ζ ) = u p ( x ; ζ ) + u ¯ p ( x ; ζ ) , D ( x ; ζ ) = d p ( x ; ζ ) + d ¯ p ( x ; ζ ) , and Σ ( x ; ζ ) = s p ( x ; ζ ) + s ¯ p ( x ; ζ ) + c p ( x ; ζ ) + c ¯ p ( x ; ζ ) . Supposing that valence quarks dominate on x 1 , then the limiting cases d p ( x ) 0 and u p ( x ) 0 yield the Nachtmann bounds [457]:
1 / 4 F 2 n ( x ) / F 2 p ( x ) 4 on x 1 .
The ζ = 2.7 GeV = : ζ M CSM prediction for F 2 n ( x ) / F 2 p ( x ) is drawn in Figure 26B. Its comparison with modern data [445] (MARATHON) may be quantified by noting that the central curve yields χ 2 / degree-of-freedom = 1.3 . It is worth stressing that the x-dependence of the CSM prediction in Figure 26B was made without reference to any data. Consequently, the agreement with the results published in Reference [445] (MARATHON) is meaningful and should serve to allay any concerns that the associated data analysis omitted some systematic effect deriving from nuclear structure modelling.
Such heightened confidence in the MARATHON data adds impact to the model-independent SPM analysis of that data described in Reference [458]; so, it is worth recapitulating some of the material therein. The final results are highlighted by Figure 27, which compares the MARATHON-based SPM prediction for F 2 n / F 2 p x 1 with: the nuclear DIS value [456]; theory predictions [271,412,446,450]; and the phenomenological fit result in Reference [434]. The figure also marks the Nachtmann lower bound, Equation (63), which is saturated if valence d quarks play no significant role at x = 1 ; namely, when there are practically no valence d quarks in the proton: d p / u p x 1 = 0 . This outcome is characteristic of proton wave function models in which the valence d quark is (almost) always paired with one of the valence u quarks inside a scalar diquark [255,449,459]. Even allowing for the quark exchange dynamics in Figure 8, one still finds d p / u p x 1 0 if only scalar diquarks are retained [450].
The following observations serve as a summary of the analyses in Reference [458]:
Observation A
… Applied to MARATHON data, the SPM yields F 2 n / F 2 p x 1 = 0.437 ( 85 ) d p / u p x 1 = 0.227 ( 100 ) .8 The possibility d p / u p x 1 = 0 is thus excluded with a 98.7% level of confidence; hence, scalar-diquark-only models of proton structure are excluded with equal likelihood. On the other hand, with this same 98.7% level of confidence, the SPM analysis confirms the QCD parton model prediction [412,446]: d p ( x ) u p ( x ) on x 1 .
Observation B
… The value of F 2 n / F 2 p x 1 inferred from nuclear DIS [456] agrees with the SPM prediction; hence, they may be averaged to yield
F 2 n / F 2 p x 1 SPM & DIS A = 0.454 ± 0.047 .
This result is drawn on Row 3 in Figure 27. It corresponds to
lim x 1 d p ( x ) u p ( x ) = 0.230 ± 0.057
and entails that the probability that scalar-diquark-only models of proton structure might be consistent with available data is 1 / 141 , 000 . In fact, with a high level of confidence, one may discard any proton structure model that delivers a result for F 2 n / F 2 p x 1 that differs significantly from Equation (64). (As reviewed in Section 9, the ratio g A d / g A u places an even harder exclusion bound on scalar-diquark-only models.)
Observation C
… Within uncertainties, the result in Equation (64) agrees with both: (i) the value obtained by assuming an SU ( 4 ) -symmetric spin–flavour wave function for the proton and helicity conservation in high- Q 2 interactions [412,446]; and (ii) the prediction developed from proton Faddeev wave functions that contain both scalar and axial vector diquarks, with the axial vector contributing approximately 40% of the proton charge [271,450]. (Recall that this was the axial vector diquark fraction built into the analyses in References [134,135].)
Following common practice, Reference [135] (Table 2) lists low-order ζ = ζ 2 , ζ 3 Mellin moments of all proton and pion DFs. Given P1 and the character of the hadron scale, then comparable momentum fractions in the proton and pion are necessarily identical and, as anticipated above, the total sea quark momentum fraction is shared between the quark flavours in roughly inverse proportion to their infrared dressed mass, M f . Significantly, given that they are predictions, the calculated values of the proton DF moments in Reference [135] (Table 2) are in fair agreement with those produced by phenomenological fits—see e.g., Reference [419] (Table VI): referring to the CT18 column, the CSM results match at the level of 1.7 ( 1.5 ) σ . Furthermore, within mutual uncertainties, the pion valence quark DF moments agree with recent lQCD results [461,462].
It is worth emphasising that the quantitative similarities also extend to the c quark: Reference [135] predicts x c p ζ 2 = 1.32 ( 5 ) %, x c p ζ 3 = 1.82 ( 6 ) %, cf. 1.7 ( 4 ) , 2.5 ( 4 ) % in Reference [435] (Figure 60). Moreover, the CSM study predicts x c ζ = 1.5 GeV = 0.64 ( 3 ) % in both the pion and proton. Regarding the pion, nothing is known about this momentum fraction; and in the proton, phenomenological estimates are inconclusive, ranging from 0–2% [435] (Figure 59). Plainly, a significant c quark momentum fraction is obtained under P1 evolution without recourse to “intrinsic charm” [463]. This outcome, which is independent of the explicit form of α 1 ( k 2 ) , potentially challenges the findings in Reference [464]. Notwithstanding the size of these calculated fractions, we stress that, as apparent in Figure 25, S c π , p ( x ) have sea quark profiles.
Given these observations, one is led to re-evaluate what is meant by intrinsic charm in the proton or any other hadron. With ζ H being the scale at which valence quasiparticle degrees-of-freedom carry all measurable properties of a given hadron, then the Fock space components, which might be interpreted as intrinsic charm or intrinsic strangeness, etc., are sublimated into the nonperturbatively computed ζ = ζ H Schwinger functions that completely express the bound state’s structure. In this context, such Fock space components are interpreted as being members of the set of basis eigenvectors representing a free-field light-front Hamiltonian. One may exemplify this by noting that any true QCD solution for the dressed quark Schwinger function (propagator) must contain infinitely many and all possible contributing Fock space vectors. The putative “intrinsic” components within a bound state are then revealed by evolving the hadron-scale Schwinger functions to higher scales, whereat an interpretation of data in terms of a Fock space expansion is relevant and practicable. In one study or another, the actual expressions of the characters of intrinsic charm, strangeness, etc., will depend on the sophistication of the kernels used to calculate the hadron-scale Schwinger functions. Nevertheless, as highlighted by Reference [135], any reasonable kernels will predict that a measurable fraction of the proton’s light-front momentum is carried by the charm quark sea at all resolving scales for which data may be interpreted in terms of DFs.
With a symmetry-preserving framework in hand that has the demonstrated ability to provide simultaneous predictions for the entire array of proton and pion DFs, one is potentially in a position to bring a new order to the study of hadron structure functions. Stringent new tests of the approach, including P1—see page 98—and the general character of the hadron scale, Equation (57), will be found in, amongst other things, studies of helicity-dependent DFs. New insights into proton spin structure may then be forthcoming.

12. Conclusions

We sketched some recent advances in the use of continuum Schwinger function methods (CSMs) to link QCD and hadron observables. The connecting bridge from theory to observation is supported by the three pillars of emergent hadrons mass (EHM): (i) dynamical generation of a gluon mass scale, whose size is roughly one half the proton mass; (ii) existence of a unique, process-independent effective charge, α ^ ( k 2 ) , which runs to a finite value at infrared momenta, α ^ ( 0 ) / π 1 ; and (iii) emergence of a running quark mass in the chiral limit, whose infrared value matches that typically identified as the constituent quark mass. Our subsequent commentary stressed that the single phenomenon of EHM manifests itself differently in the diverse array of measurable quantities that define hadron physics. No single observable is alone sufficient to validate the EHM paradigm for understanding strong interactions within the Standard Model (SM). Instead, theory should identify a broad range of empirical consequences of EHM so that the order brought to a collective body of experimental results—existing and future—can be recognised as the signature of EHM.
Our developing understanding of EHM suggests that QCD is unique amongst known fundamental theories of natural phenomena. It might be the first well-defined four-dimensional quantum field theory ever contemplated. If so, then QCD could provide the archetype for theories that take physics beyond the SM.
Science has delivered theories of many things. The best of them remain a part of the grander theories developed in response to new observations. The basic question yet remains unanswered, viz.: Is there a theory of everything? Hadron physics and QCD might be pointing us toward an answer in exposing the special qualities of strong-coupling non-Abelian quantum gauge field theories.

Funding

Work supported by: National Natural Science Foundation of China (Grant No. 12135007); and Helmholtz-Zentrum Dresden Rossendorf under the High Potential Programme.

Data Availability Statement

Not applicable.

Acknowledgments

This contribution is based on results obtained and insights developed through collaborations with many people, to all of whom we are greatly indebted.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACManomalous chromomagnetic moment
AdS/CFT (duality)anti-de Sitter/conformal field theory (duality)
ard ¯ mean absolute relative difference
CKMCabibbo–Kobayashi–Maskawa (matrix)
CSMscontinuum Schwinger function methods
DCSBdynamical chiral symmetry breaking
DF(parton) distribution function
DISdeep inelastic scattering
DSEDyson–Schwinger equation
EHMemergent hadron mass
FF(parton) fragmentation function
JLabThomas Jefferson National Accelerator Facility
lQCDlattice-regularised quantum chromodynamics
NG (mode/boson)Nambu–Goldstone (mode/boson)
PD (charge)process-dependent (charge)
PDGParticle Data Group and associated publications
PI (charge)process-independent (charge)
QCDquantum chromodynamics
QEDquantum electrodynamics
RGIrenormalisation-group-invariant
RLrainbow ladder (truncation)
SCIsymmetry-preserving treatment of a vector×vector contact interaction
SMStandard Model of particle physics
SPMSchlessinger point method
VMDvector meson dominance
1
Here, the value “10” is arbitrary. More generally, the number should be small enough to ensure that predictive power is not lost through a need to fit too many renormalised observables to measured quantities.
2
When working with a Euclidean formulation, as we do, Poincaré covariance maps straightforwardly into Euclidean covariance, viz. valid Schwinger functions must transform covariantly under O ( 4 ) rotations and linear translations in R 4 . Owing to the simplicity of this connection, we avoid transliteration and speak of Poincaré covariance and invariance throughout.
3
This is well known and explains why the truncated interaction in Equation (26) was not used for any of the calculations described in Section 7, Section 8, Section 9, Section 10 and Section 11 below. All those studies are based on interactions that at least preserve QCD’s ultraviolet power-law behaviour, where more has not yet been achieved—Section 7 and Section 9—and also the one-loop logarithmic improvement, when the necessary algorithms are already available—Section 8, Section 10, and Section 11.
4
Dressed-quark propagators form an important part of the kernels of all bound state equations. As on-shell meson masses increase, poles in those propagators enter the complex plane integration domain sampled by the Bethe–Salpeter equation [199]. For such cases—here, meson excited states—a direct on-shell solution cannot be obtained using simple algorithms. Therefore, to obtain the masses of those mesons, Reference [222] used an extrapolation procedure based on Padé approximants. This is the origin of the uncertainty bar on the CSM predictions.
5
The Chebyshev (or hyperspherical) expansion of Poincaré-invariant functions of two scalar variables is discussed, e.g., in Reference [199] (IV.B).
6
If one eliminates axial vector diquarks from the proton wave function, then g A d / g A u = 0.054 ( 13 ) , a result disfavoured by experiments at the level of 5.1 σ , i.e., the probability that the scalar-diquark-only proton result could be consistent with data is 1 / 7 , 100 , 000 .
7
Contemporary continuum methods for obtaining light-front amplitudes and density distributions from Euclidean space Schwinger functions are detailed, e.g., in References [134,423]; [35] – Sections 3, 5; [125] – Section IV; [127] – Sections 2, 5.
8
Extrapolations based on [ 1 , 1 ] Padé fits to MARATHON data, obtained using a one-point jackknife procedure, yield F 2 n / F 2 p = 0.395 ( 3 ) on x 1 d p / u d = 0.169 ( 3 ) [182]. Another analysis [460], employing practitioner-chosen polynomials as the basis for extrapolation, obtains F 2 n / F 2 p = 0.37 ( 7 ) on x 1 d p / u d = 0.13 ( 8 ) . The latter is less precise, but both results are consistent with the function form unbiased SPM prediction.

References

  1. Appelquist, T.; Pisarski, R.D. High-Temperature Yang–Mills Theories and Three-Dimensional Quantum Chromodynamics. Phys. Rev. D 1981, 23, 2305. [Google Scholar] [CrossRef]
  2. Appelquist, T.W.; Bowick, M.J.; Karabali, D.; Wijewardhana, L.C.R. Spontaneous Chiral Symmetry Breaking in Three-Dimensional QED. Phys. Rev. D 1986, 33, 3704. [Google Scholar] [CrossRef] [PubMed]
  3. Bashir, A.; Raya, A.; Cloet, I.C.; Roberts, C.D. Regarding confinement and dynamical chiral symmetry breaking in QED3. Phys. Rev. C 2008, 78, 055201. [Google Scholar] [CrossRef] [Green Version]
  4. Bashir, A.; Raya, A.; Sánchez-Madrigal, S.; Roberts, C.D. Gauge invariance of a critical number of flavours in QED3. Few Body Syst. 2009, 46, 229–237. [Google Scholar] [CrossRef] [Green Version]
  5. Braun, J.; Gies, H.; Janssen, L.; Roscher, D. Phase structure of many-flavor QED3. Phys. Rev. D 2014, 90, 036002. [Google Scholar] [CrossRef] [Green Version]
  6. Roberts, C.D. Perspective on the origin of hadron masses. Few Body Syst. 2017, 58, 5. [Google Scholar] [CrossRef] [Green Version]
  7. Glimm, J.; Jaffee, A. Quantum Physics: A Functional Point of View; Springer: New York, NY, USA, 1981. [Google Scholar]
  8. Seiler, E. Gauge Theories as a Problem of Constructive Quantum Theory and Statistical Mechanics; Springer: New York, NY, USA, 1982. [Google Scholar]
  9. Itzykson, C.; Zuber, J.B. Quantum Field Theory; McGraw-Hill Inc.: New York, NY, USA, 1980. [Google Scholar]
  10. Rakow, P.E.L. Renormalization group flow in QED: An investigation of the Schwinger-Dyson equations. Nucl. Phys. B 1991, 356, 27–45. [Google Scholar] [CrossRef]
  11. Gockeler, M.; Horsley, R.; Linke, V.; Rakow, P.E.L.; Schierholz, G.; Stuben, H. Probing the continuum limit in noncompact QED: New results on large lattices. Nucl. Phys. B Proc. Suppl. 1995, 42, 660–662. [Google Scholar] [CrossRef] [Green Version]
  12. Reenders, M. On the nontriviality of Abelian gauged Nambu-Jona-Lasinio models in four dimensions. Phys. Rev. D 2000, 62, 025001. [Google Scholar] [CrossRef] [Green Version]
  13. Kızılersü, A.; Sizer, T.; Pennington, M.R.; Williams, A.G.; Williams, R. Dynamical mass generation in unquenched QED using the Dyson–Schwinger equations. Phys. Rev. D 2015, 91, 065015. [Google Scholar] [CrossRef]
  14. Roberts, C.D.; Schmidt, S.M. Dyson–Schwinger equations: Density, temperature and continuum strong QCD. Prog. Part. Nucl. Phys. 2000, 45, S1–S103. [Google Scholar] [CrossRef]
  15. Roberts, C.D.; Williams, A.G.; Krein, G. On the implications of confinement. Int. J. Mod. Phys. A 1992, 7, 5607–5624. [Google Scholar] [CrossRef]
  16. Roberts, C.D.; Williams, A.G. Dyson–Schwinger equations and their application to hadronic physics. Prog. Part. Nucl. Phys. 1994, 33, 477–575. [Google Scholar] [CrossRef] [Green Version]
  17. Binosi, D.; Tripolt, R.A. Spectral functions of confined particles. Phys. Lett. B 2020, 801, 135171. [Google Scholar] [CrossRef]
  18. Wilson, K.G. Confinement of quarks. Phys. Rev. D 1974, 10, 2445–2459. [Google Scholar] [CrossRef]
  19. Wilson, K.G. The origins of lattice gauge theory. Nucl. Phys. Proc. Suppl. 2005, 140, 3–19. [Google Scholar] [CrossRef] [Green Version]
  20. Cali, S.; Detmold, W.; Korcyl, G.; Korcyl, P.; Shanahan, P. Implementation of the conjugate gradient algorithm for heterogeneous systems. In Proceedings of the 38th International Symposium on Lattice Field Theory—PoS(LATTICE2021), Cambridge, MA, USA, 26–30 July 2021; Volume 396. [Google Scholar]
  21. Maris, P.; Roberts, C.D. Dyson–Schwinger equations: A tool for hadron physics. Int. J. Mod. Phys. E 2003, 12, 297–365. [Google Scholar] [CrossRef] [Green Version]
  22. Roberts, C.D. Hadron Properties and Dyson–Schwinger Equations. Prog. Part. Nucl. Phys. 2008, 61, 50–65. [Google Scholar] [CrossRef] [Green Version]
  23. Chang, L.; Roberts, C.D.; Tandy, P.C. Selected highlights from the study of mesons. Chin. J. Phys. 2011, 49, 955–1004. [Google Scholar]
  24. Bashir, A.; Chang, L.; Cloet, I.C.; El-Bennich, B.; Liu, Y.X.; Roberts, C.D.; Tandy, P.C. Collective perspective on advances in Dyson–Schwinger Equation QCD. Commun. Theor. Phys. 2012, 58, 79–134. [Google Scholar] [CrossRef]
  25. Roberts, C.D. Strong QCD and Dyson–Schwinger Equations. IRMA Lect. Math. Theor. Phys. 2015, 21, 356–458. [Google Scholar]
  26. Roberts, C.D. Three Lectures on Hadron Physics. J. Phys. Conf. Ser. 2016, 706, 022003. [Google Scholar] [CrossRef]
  27. Horn, T.; Roberts, C.D. The pion: An enigma within the Standard Model. J. Phys. G. 2016, 43, 073001. [Google Scholar] [CrossRef] [Green Version]
  28. Eichmann, G.; Sanchis-Alepuz, H.; Williams, R.; Alkofer, R.; Fischer, C.S. Baryons as relativistic three-quark bound states. Prog. Part. Nucl. Phys. 2016, 91, 1–100. [Google Scholar] [CrossRef] [Green Version]
  29. Burkert, V.D.; Roberts, C.D. Roper resonance: Toward a solution to the fifty-year puzzle. Rev. Mod. Phys. 2019, 91, 011003. [Google Scholar] [CrossRef] [Green Version]
  30. Fischer, C.S. QCD at finite temperature and chemical potential from Dyson–Schwinger equations. Prog. Part. Nucl. Phys. 2019, 105, 1–60. [Google Scholar] [CrossRef] [Green Version]
  31. Qin, S.X.; Roberts, C.D. Impressions of the Continuum Bound State Problem in QCD. Chin. Phys. Lett. 2020, 37, 121201. [Google Scholar] [CrossRef]
  32. Roberts, C.D.; Schmidt, S.M. Reflections upon the Emergence of Hadronic Mass. Eur. Phys. J. Spec. Top. 2020, 229, 3319–3340. [Google Scholar] [CrossRef]
  33. Roberts, C.D. Empirical Consequences of Emergent Mass. Symmetry 2020, 12, 1468. [Google Scholar] [CrossRef]
  34. Roberts, C.D. On Mass and Matter. AAPPS Bull. 2021, 31, 6. [Google Scholar] [CrossRef]
  35. Roberts, C.D.; Richards, D.G.; Horn, T.; Chang, L. Insights into the emergence of mass from studies of pion and kaon structure. Prog. Part. Nucl. Phys. 2021, 120, 103883. [Google Scholar] [CrossRef]
  36. Binosi, D. Emergent Hadron Mass in Strong Dynamics. Few Body Syst. 2022, 63, 42. [Google Scholar] [CrossRef]
  37. Papavassiliou, J. Emergence of mass in the gauge sector of QCD. Chin. Phys. C 2022, 46, 112001. [Google Scholar] [CrossRef]
  38. Englert, F. Nobel Lecture: The BEH mechanism and its scalar boson. Rev. Mod. Phys. 2014, 86, 843. [Google Scholar] [CrossRef] [Green Version]
  39. Higgs, P.W. Nobel Lecture: Evading the Goldstone theorem. Rev. Mod. Phys. 2014, 86, 851. [Google Scholar] [CrossRef] [Green Version]
  40. Workman, R.L.; Burkert, V.D.; Crede, V.; Klempt, E.; Thoma, U.; Tiator, L.; Agashe, K.; Aielli, G.; Allanach, B.C.; Amsler, C.; et al. Review of Particle Physics. Prog. Theor. Exp. Phys. 2022, 2022, 083C01. [Google Scholar]
  41. Flambaum, V.V.; Höll, A.; Jaikumar, P.; Roberts, C.D.; Wright, S.V. Sigma terms of light quark hadrons. Few Body Syst. 2006, 38, 31–51. [Google Scholar] [CrossRef] [Green Version]
  42. Ruiz de Elvira, J.; Hoferichter, M.; Kubis, B.; Meißner, U.G. Extracting the σ-term from low-energy pion–nucleon scattering. J. Phys. G 2018, 45, 024001. [Google Scholar] [CrossRef] [Green Version]
  43. Aoki, S.; Aoki, Y.; Bečirević, D.; Blum, T.; Colangelo, G.; Collins, S.; Della Morte, M.; Dimopoulos, P.; Dürr, S.; Fukaya, H.; et al. FLAG Review 2019. Eur. Phys. J. C 2020, 80, 113. [Google Scholar]
  44. Nambu, Y. Quasiparticles and Gauge Invariance in the Theory of Superconductivity. Phys. Rev. 1960, 117, 648–663. [Google Scholar] [CrossRef]
  45. Goldstone, J. Field Theories with Superconductor Solutions. Nuovo Cim. 1961, 19, 154–164. [Google Scholar] [CrossRef] [Green Version]
  46. Gell-Mann, M.; Oakes, R.J.; Renner, B. Behavior of current divergences under SU(3) x SU(3). Phys. Rev. 1968, 175, 2195–2199. [Google Scholar] [CrossRef] [Green Version]
  47. Casher, A.; Susskind, L. Chiral magnetism (or magnetohadrochironics). Phys. Rev. 1974, D9, 436–460. [Google Scholar] [CrossRef]
  48. Brodsky, S.J.; Shrock, R. Condensates in Quantum Chromodynamics and the Cosmological Constant. Proc. Natl. Acad. Sci. USA 2011, 108, 45–50. [Google Scholar] [CrossRef]
  49. Brodsky, S.J.; Roberts, C.D.; Shrock, R.; Tandy, P.C. New perspectives on the quark condensate. Phys. Rev. C 2010, 82, 022201. [Google Scholar] [CrossRef] [Green Version]
  50. Chang, L.; Roberts, C.D.; Tandy, P.C. Expanding the concept of in-hadron condensates. Phys. Rev. C 2012, 85, 012201. [Google Scholar] [CrossRef] [Green Version]
  51. Brodsky, S.J.; Roberts, C.D.; Shrock, R.; Tandy, P.C. Confinement contains condensates. Phys. Rev. C 2012, 85, 065202. [Google Scholar] [CrossRef] [Green Version]
  52. Cloet, I.C.; Roberts, C.D. Explanation and Prediction of Observables using Continuum Strong QCD. Prog. Part. Nucl. Phys. 2014, 77, 1–69. [Google Scholar] [CrossRef] [Green Version]
  53. Pascual, P.; Tarrach, R. QCD: Renormalization for the Practitioner, Lecture Notes in Physics; Springer: Berlin/Heidelberg, Germany, 1984; Volume 194. [Google Scholar]
  54. Taylor, J.C. Ward Identities and Charge Renormalization of the Yang–Mills Field. Nucl. Phys. B 1971, 33, 436–444. [Google Scholar] [CrossRef]
  55. Slavnov, A.A. Ward Identities in Gauge Theories. Theor. Math. Phys. 1972, 10, 99–107. [Google Scholar] [CrossRef]
  56. Schwinger, J.S. Gauge Invariance and Mass. Phys. Rev. 1962, 125, 397–398. [Google Scholar] [CrossRef]
  57. Schwinger, J.S. Gauge Invariance and Mass. 2. Phys. Rev. 1962, 128, 2425–2429. [Google Scholar] [CrossRef]
  58. Appelquist, T.; Nash, D.; Wijewardhana, L. Critical Behavior in (2+1)-Dimensional QED. Phys. Rev. Lett. 1988, 60, 2575. [Google Scholar] [CrossRef] [PubMed]
  59. Maris, P. The influence of the full vertex and vacuum polarization on the fermion propagator in QED3. Phys. Rev. 1996, D54, 4049–4058. [Google Scholar]
  60. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Nonperturbative gluon and ghost propagators for d=3 Yang–Mills. Phys. Rev. D 2010, 81, 125025. [Google Scholar] [CrossRef]
  61. Cornwall, J.M. Exploring dynamical gluon mass generation in three dimensions. Phys. Rev. D 2016, 93, 025021. [Google Scholar] [CrossRef] [Green Version]
  62. Cornwall, J.M. Dynamical Mass Generation in Continuum QCD. Phys. Rev. D 1982, 26, 1453. [Google Scholar] [CrossRef]
  63. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. Gluon and ghost propagators in the Landau gauge: Deriving lattice results from Schwinger-Dyson equations. Phys. Rev. D 2008, 78, 025010. [Google Scholar] [CrossRef] [Green Version]
  64. Boucaud, P.; Leroy, J.; Le Yaouanc, A.; Micheli, J.; Pene, O.; Rodríguez-Quintero, J. On the IR behaviour of the Landau-gauge ghost propagator. J. High Energy Phys. 2008, 0806, 099. [Google Scholar] [CrossRef] [Green Version]
  65. Binosi, D.; Papavassiliou, J. Pinch Technique: Theory and Applications. Phys. Rept. 2009, 479, 1–152. [Google Scholar] [CrossRef] [Green Version]
  66. Boucaud, P.; Leroy, J.P.; Le-Yaouanc, A.; Micheli, J.; Pene, O.; Rodríguez-Quintero, J. The Infrared Behaviour of the Pure Yang–Mills Green Functions. Few Body Syst. 2012, 53, 387–436. [Google Scholar] [CrossRef] [Green Version]
  67. Aguilar, A.C.; Binosi, D.; Papavassiliou, J. The Gluon Mass Generation Mechanism: A Concise Primer. Front. Phys. China 2016, 11, 111203. [Google Scholar] [CrossRef] [Green Version]
  68. Binosi, D.; Chang, L.; Papavassiliou, J.; Roberts, C.D. Bridging a gap between continuum-QCD and ab initio predictions of hadron observables. Phys. Lett. B 2015, 742, 183–188. [Google Scholar] [CrossRef] [Green Version]
  69. Aguilar, A.C.; Ferreira, M.N.; Papavassiliou, J. Exploring smoking-gun signals of the Schwinger mechanism in QCD. Phys. Rev. D 2022, 105, 014030. [Google Scholar] [CrossRef]
  70. Pinto-Gómez, F.; De Soto, F.; Ferreira, M.N.; Papavassiliou, J.; Rodríguez-Quintero, J. Lattice three-gluon vertex in extended kinematics: Planar degeneracy. arXiv 2022, arXiv:2208.01020. [Google Scholar]
  71. Cui, Z.F.; Zhang, J.L.; Binosi, D.; de Soto, F.; Mezrag, C.; Papavassiliou, J.; Roberts, C.D.; Rodríguez-Quintero, J.; Segovia, J.; Zafeiropoulos, S. Effective charge from lattice QCD. Chin. Phys. C 2020, 44, 083102. [Google Scholar] [CrossRef]
  72. Aguilar, A.C.; De Soto, F.; Ferreira, M.; Papavassiliou, J.; Rodríguez-Quintero, J.; Zafeiropoulos, S. Gluon propagator and three-gluon vertex with dynamical quarks. Eur. Phys. J. C 2020, 80, 154. [Google Scholar] [CrossRef] [Green Version]
  73. Blum, T.; Boyle, P.A.; Christ, N.H.; Frison, J.; Garron, N.; Hudspith, R.J.; Izubuchi, T.; Janowski, T.; Jung, C.; Jüttner, A.; et al. Domain wall QCD with physical quark masses. Phys. Rev. D 2016, 93, 074505. [Google Scholar] [CrossRef] [Green Version]
  74. Boyle, P.A.; Christ, N.H.; Garron, N.; Jung, C.; Jüttner, A.; Kelly, C.; Mawhinney, R.D.; McGlynn, G.; Murphy, D.J.; Ohta, S.; et al. Low energy constants of SU(2) partially quenched chiral perturbation theory from Nf=2+1 domain wall QCD. Phys. Rev. D 2016, 93, 054502. [Google Scholar] [CrossRef] [Green Version]
  75. Boyle, P.A.; Del Debbio, L.; Jüttner, A.; Khamseh, A.; Sanfilippo, F.; Tsang, J.T. The decay constants fD and fDs in the continuum limit of Nf = 2 + 1 domain wall lattice QCD. J. High Energy Phys. 2017, 12, 008. [Google Scholar] [CrossRef] [Green Version]
  76. Binosi, D.; Chang, L.; Papavassiliou, J.; Qin, S.X.; Roberts, C.D. Natural constraints on the gluon–quark vertex. Phys. Rev. D 2017, 95, 031501. [Google Scholar] [CrossRef] [Green Version]
  77. Brodsky, S.J.; Deur, A.; Roberts, C.D. Artificial dynamical effects in quantum field theory. Nat. Rev. Phys. 2022, 4, 489–495. [Google Scholar] [CrossRef]
  78. West, G.B. General Infrared and Ultraviolet Properties of the Gluon Propagator in Axial Gauge. Phys. Rev. D 1983, 27, 1878. [Google Scholar] [CrossRef] [Green Version]
  79. Brown, N.; Pennington, M.R. Studies of Confinement: How the Gluon Propagates. Phys. Rev. D 1989, 39, 2723. [Google Scholar] [CrossRef] [PubMed]
  80. Cucchieri, A.; Mendes, T.; Santos, E.M.S. Covariant gauge on the lattice: A New implementation. Phys. Rev. Lett. 2009, 103, 141602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Gao, F.; Qin, S.X.; Roberts, C.D.; Rodríguez-Quintero, J. Locating the Gribov horizon. Phys. Rev. D 2018, 97, 034010. [Google Scholar] [CrossRef]
  82. Gell-Mann, M.; Low, F.E. Quantum electrodynamics at small distances. Phys. Rev. 1954, 95, 1300–1312. [Google Scholar] [CrossRef] [Green Version]
  83. Politzer, H.D. The dilemma of attribution. Proc. Natl. Acad. Sci. USA 2005, 102, 7789–7793. [Google Scholar] [CrossRef] [Green Version]
  84. Wilczek, F. Asymptotic freedom: From paradox to paradigm. Proc. Natl. Acad. Sci. USA 2005, 102, 8403–8413. [Google Scholar] [CrossRef] [Green Version]
  85. Gross, D.J. The discovery of asymptotic freedom and the emergence of QCD. Proc. Natl. Acad. Sci. USA 2005, 102, 9099–9108. [Google Scholar] [CrossRef] [Green Version]
  86. Pickering, A. Constructing Quarks. A Sociological History of Particle Physics; University of Chicago Press: Chicago, IL, USA, 1999. [Google Scholar]
  87. Marciano, W.J.; Pagels, H. Quantum Chromodynamics: A Review. Phys. Rept. 1978, 36, 137. [Google Scholar] [CrossRef]
  88. Fischer, C.S. Infrared properties of QCD from Dyson–Schwinger equations. J. Phys. G 2006, 32, R253–R291. [Google Scholar] [CrossRef] [Green Version]
  89. Chang, L.; Roberts, C.D. Sketching the Bethe–Salpeter kernel. Phys. Rev. Lett. 2009, 103, 081601. [Google Scholar] [CrossRef] [Green Version]
  90. Binosi, D.; Mezrag, C.; Papavassiliou, J.; Roberts, C.D.; Rodríguez-Quintero, J. Process-independent strong running coupling. Phys. Rev. D 2017, 96, 054026. [Google Scholar] [CrossRef] [Green Version]
  91. Pilaftsis, A. Generalized pinch technique and the background field method in general gauges. Nucl. Phys. B 1997, 487, 467–491. [Google Scholar] [CrossRef] [Green Version]
  92. Cornwall, J.M.; Papavassiliou, J.; Binosi, D. The Pinch Technique and its Applications to Non-Abelian Gauge Theories; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  93. Abbott, L.F. Introduction to the Background Field Method. Acta Phys. Polon. B 1982, 13, 33. [Google Scholar]
  94. Deur, A.; Brodsky, S.J.; de Teramond, G.F. The QCD Running Coupling. Prog. Part. Nucl. Phys. 2016, 90, 1–74. [Google Scholar] [CrossRef]
  95. Deur, A.; Burkert, V.; Chen, J.P.; Korsch, W. Experimental determination of the QCD effective charge αg1(Q). Particles 2022, 5, 171. [Google Scholar] [CrossRef]
  96. Deur, A.; Burkert, V.; Chen, J.P.; Korsch, W. Experimental determination of the effective strong coupling constant. Phys. Lett. B 2007, 650, 244–248. [Google Scholar] [CrossRef] [Green Version]
  97. Deur, A.; Burkert, V.; Chen, J.P.; Korsch, W. Determination of the effective strong coupling constant αg1(s) from CLAS spin structure function data. Phys. Lett. B 2008, 665, 349–351. [Google Scholar] [CrossRef]
  98. Deur, A.; Prok, Y.; Burkert, V.; Crabb, D.; Girod, F.X.; Griffioen, K.A.; Guler, N.; Kuhn, S.E.; Kvaltine, N. High precision determination of the Q2 evolution of the Bjorken Sum. Phys. Rev. D 2014, 90, 012009. [Google Scholar] [CrossRef] [Green Version]
  99. The HERMES Collaboration; Ackerstaff, K. Measurement of the neutron spin structure function g 1 n with a polarized 3He internal target. Phys. Lett. B 1997, 404, 383–389. [Google Scholar] [CrossRef] [Green Version]
  100. Ackerstaff, K.; Airapetian, A.; Akopov, N.; Akushevich, I.; Amarian, M.; Aschenauer, E.; Avakian, H.; Avetissian, A.; Bains, B.; Barrow, S.; et al. Determination of the deep inelastic contribution to the generalized Gerasimov-Drell-Hearn integral for the proton and neutron. Phys. Lett. B 1998, 444, 531–538. [Google Scholar] [CrossRef] [Green Version]
  101. Airapetian, A.; Akopov, N.; Akushevich, I.; Amarian, M.; Aschenauer, E.; Avakian, H.; Avetissian, A.; Bains, B.; Baumgarten, C.; Beckmann, M.; et al. Measurement of the proton spin structure function g 1 p with a pure hydrogen target. Phys. Lett. B 1998, 442, 484–492. [Google Scholar] [CrossRef] [Green Version]
  102. Airapetian, A.; Akopov, N.; Akopov, Z.; Amarian, M.; Ammosov, V.V.; Aschenauer, E.C.; Avakian, H.; Avetissian, A.; Bailey, P.; Baturin, V.; et al. Evidence for quark hadron duality in the proton spin asymmetry A1. Phys. Rev. Lett. 2003, 90, 092002. [Google Scholar] [CrossRef] [Green Version]
  103. Airapetian, A.; Akopov, N.; Akopov, Z.; Andrus, A.; Aschenauer, E.C.; Augustyniak, W.; Tchuiko, B. Precise determination of the spin structure function g1 of the proton, deuteron and neutron. Phys. Rev. D 2007, 75, 012007. [Google Scholar] [CrossRef] [Green Version]
  104. NuTeV Collaboration. A Measurement of αs(Q2) from the Gross-Llewellyn Smith sum rule. Phys. Rev. Lett. 1998, 81, 3595–3598. [Google Scholar] [CrossRef]
  105. Alexakhin, V.Y.; Alexandrov, Y.; Alexeev, G.D.; Alexeev, M.; Amoroso, A.; Badełek, B.; Mikhailov, Y.V. The Deuteron Spin-dependent Structure Function g 1 d and its First Moment. Phys. Lett. B 2007, 647, 8–17. [Google Scholar] [CrossRef] [Green Version]
  106. Alekseev, M.G.; Alexakhin, V.Y.; Alexandrov, Y.; Alexeev, G.D.; Amoroso, A.; Austregesilo, A.; Padee, A. The Spin-dependent Structure Function of the Proton g 1 p and a Test of the Bjorken Sum Rule. Phys. Lett. B 2010, 690, 466–472. [Google Scholar] [CrossRef] [Green Version]
  107. COMPASS Collaboration; Bordalo, P.; Franco, C.; Nunes, A.S.; Quaresma, M.; Quintans, C.; Ramos, S.; Silva, L.; Stolarski, M. The spin structure function g 1 p of the proton and a test of the Bjorken sum rule. Phys. Lett. B 2016, 753, 18–28. [Google Scholar] [CrossRef]
  108. Anthony, P.L.; Arnold, R.G.; Band, H.R.; Borel, H.; Bosted, P.E.; Breton, V.; Cates, G.D.; Chupp, T.E.; Dietrich, F.S.; Dunne, J.; et al. Determination of the neutron spin structure function. Phys. Rev. Lett. 1993, 71, 959–962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. E143 Collaboration; Abe, K.; Akagi, T.; Anthony, P.L.; Antonov, R.; Arnold, R.G.; Averett, T. Precision measurement of the proton spin structure function g 1 p . Phys. Rev. Lett. 1995, 74, 346–350. [Google Scholar] [CrossRef] [Green Version]
  110. E143 Collaboration; Abe, K.; Akagi, T.; Anthony, P.L.; Antonov, R.; Arnold, R.G.; Averett, T. Precision measurement of the deuteron spin structure function g 1 d . Phys. Rev. Lett. 1995, 75, 25–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  111. Abe, K.; Akagi, T.; Anthony, P.L.; Antonov, R.; Arnold, R.G.; Averett, T.; Band, H.R.; Bauer, J.M.; Borel, H.; Bosted, P.E.; et al. Measurements of the proton and deuteron spin structure function g2 and asymmetry A2. Phys. Rev. Lett. 1996, 76, 587–591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Abe, K.; Akagi, T.; Anthony, P.; Antonov, R.; Arnold, R.; Averett, T.; Band, H.; Bauer, J.; Borel, H.; Bosted, P.; et al. Measurements of the Q2 dependence of the proton and deuteron spin structure functions g 1 p and g 1 d . Phys. Lett. B 1995, 364, 61–68. [Google Scholar] [CrossRef] [Green Version]
  113. Anthony, P.L.; Arnold, R.G.; Band, H.R.; Borel, H.; Bosted, P.E.; Breton, V.; Cates, G.D.; Chupp, T.E.; Dietrich, F.S.; Dunne, J.; et al. Deep inelastic scattering of polarized electrons by polarized 3He and the study of the neutron spin structure. Phys. Rev. D 1996, 54, 6620–6650. [Google Scholar] [CrossRef] [Green Version]
  114. Abe, K.; Akagi, T.; Anderson, B.D.; Anthony, P.L.; Arnold, R.G.; Averett, T.; Band, H.R.; Berisso, C.M.; Bogorad, P.; Borel, H.; et al. Precision determination of the neutron spin structure function g 1 n . Phys. Rev. Lett. 1997, 79, 26–30. [Google Scholar] [CrossRef]
  115. Abe, K.; Akagi, T.; Anderson, B.; Anthony, P.; Arnold, R.; Averett, T.; Band, H.; Berisso, C.; Bogorad, P.; Borel, H.; et al. Measurement of the neutron spin structure function g 2 n and asymmetry A 2 n . Phys. Lett. B 1997, 404, 377–382. [Google Scholar] [CrossRef] [Green Version]
  116. Abe, K.; Akagi, T.; Anderson, B.; Anthony, P.; Arnold, R.; Averett, T.; Band, H.; Berisso, C.; Bogorad, P.; Borel, H.; et al. Next-to-leading order QCD analysis of polarized deep inelastic scattering data. Phys. Lett. B 1997, 405, 180–190. [Google Scholar] [CrossRef] [Green Version]
  117. Abe, K.; Akagi, T.; Anthony, P.L.; Antonov, R.; Arnold, R.G.; Averett, T.; Band, H.R.; Bauer, J.M.; Borel, H.; Bosted, P.E.; et al. Measurements of the proton and deuteron spin structure functions g1 and g2. Phys. Rev. 1998, D58, 112003. [Google Scholar]
  118. E155 Collaboration; Anthony, P.L.; Arnold, R.G.; Averett, T.; Band, H.R.; Berisso, M.C.; Borelg, H.; Bosteda, P.E.; Bültmanns, S.L.; Buenerd, M.; et al. Measurement of the proton and deuteron spin structure functions g2 and asymmetry A2. Phys. Lett. 1999, B458, 529–535. [Google Scholar]
  119. Anthony, P.; Arnold, R.; Averett, T.; Band, H.; Berisso, M.; Borel, H.; Bosted, P.; Bültmann, S.; Buenerd, M.; Chupp, T.; et al. Measurement of the deuteron spin structure function g 1 d (x) for 1 (GeV/c)2 < Q2 < 40 (GeV/c)2. Phys. Lett. B 1999, 463, 339–345. [Google Scholar]
  120. Anthony, P.; Arnold, R.; Averett, T.; Band, H.; Berisso, M.; Borel, H.; Bosted, P.; Bültmann, S.; Buenerd, M.; Chupp, T.; et al. Measurements of the Q2 dependence of the proton and neutron spin structure functions g 1 p and g 1 n . Phys. Lett. B 2000, 493, 19–28. [Google Scholar] [CrossRef] [Green Version]
  121. Anthony, P.; Arnold, R.; Averett, T.; Band, H.; Benmouna, N.; Boeglin, W.; Borel, H.; Bosted, P.; Bültmann, S.; Court, G.; et al. Precision measurement of the proton and deuteron spin structure functions g2 and asymmetries A2. Phys. Lett. B 2003, 553, 18–24. [Google Scholar] [CrossRef]
  122. Binosi, D.; Quadri, A. Anti-BRST symmetry and background field method. Phys. Rev. D 2013, 88, 085036. [Google Scholar] [CrossRef] [Green Version]
  123. Aslam, M.J.; Bashir, A.; Gutierrez-Guerrero, L.X. Local Gauge Transformation for the Quark Propagator in an SU(N) Gauge Theory. Phys. Rev. D 2016, 93, 076001. [Google Scholar] [CrossRef] [Green Version]
  124. Ding, M.; Raya, K.; Binosi, D.; Chang, L.; Roberts, C.D.; Schmidt, S.M. Drawing insights from pion parton distributions. Chin. Phys. C (Lett.) 2020, 44, 031002. [Google Scholar] [CrossRef] [Green Version]
  125. Ding, M.; Raya, K.; Binosi, D.; Chang, L.; Roberts, C.D.; Schmidt, S.M. Symmetry, symmetry breaking, and pion parton distributions. Phys. Rev. D 2020, 101, 054014. [Google Scholar] [CrossRef] [Green Version]
  126. Cui, Z.F.; Ding, M.; Gao, F.; Raya, K.; Binosi, D.; Chang, L.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M. Higgs modulation of emergent mass as revealed in kaon and pion parton distributions. Eur. Phys. J. A (Lett.) 2021, 57, 5. [Google Scholar] [CrossRef]
  127. Cui, Z.F.; Ding, M.; Gao, F.; Raya, K.; Binosi, D.; Chang, L.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M. Kaon and pion parton distributions. Eur. Phys. J. C 2020, 80, 1064. [Google Scholar] [CrossRef]
  128. Han, C.; Xie, G.; Wang, R.; Chen, X. An Analysis of Parton Distribution Functions of the Pion and the Kaon with the Maximum Entropy Input. Eur. Phys. J. C 2021, 81, 302. [Google Scholar] [CrossRef]
  129. Chang, L.; Roberts, C.D. Regarding the distribution of glue in the pion. Chin. Phys. Lett. 2021, 38, 081101. [Google Scholar] [CrossRef]
  130. Xie, G.; Han, C.; Wang, R.; Chen, X. Tackling the kaon structure function at EicC. Chin. Phys. C 2022, 46, 064107. [Google Scholar] [CrossRef]
  131. Raya, K.; Cui, Z.F.; Chang, L.; Morgado, J.M.; Roberts, C.D.; Rodríguez-Quintero, J. Revealing pion and kaon structure via generalised parton distributions. Chin. Phys. C 2022, 46, 013105. [Google Scholar] [CrossRef]
  132. Cui, Z.F.; Ding, M.; Morgado, J.M.; Raya, K.; Binosi, D.; Chang, L.; Papavassiliou, J.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M. Concerning pion parton distributions. Eur. Phys. J. A 2022, 58, 10. [Google Scholar] [CrossRef]
  133. Cui, Z.F.; Ding, M.; Morgado, J.M.; Raya, K.; Binosi, D.; Chang, L.; De Soto, F.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M. Emergence of pion parton distributions. Phys. Rev. D 2022, 105, L091502. [Google Scholar] [CrossRef]
  134. Chang, L.; Gao, F.; Roberts, C.D. Parton distributions of light quarks and antiquarks in the proton. Phys. Lett. B 2022, 829, 137078. [Google Scholar] [CrossRef]
  135. Lu, Y.; Chang, L.; Raya, K.; Roberts, C.D.; Rodríguez-Quintero, J. Proton and pion distribution functions in counterpoint. Phys. Lett. B 2022, 830, 137130. [Google Scholar] [CrossRef]
  136. De Paula, W.; Ydrefors, E.; Nogueira Alvarenga, J.H.; Frederico, T.; Salmè, G. Parton distribution function in a pion with Minkowskian dynamics. Phys. Rev. D 2022, 105, L071505. [Google Scholar] [CrossRef]
  137. Grunberg, G. Renormalization Scheme Independent QCD and QED: The Method of Effective Charges. Phys. Rev. D 1984, 29, 2315. [Google Scholar] [CrossRef]
  138. Grunberg, G. On Some Ambiguities in the Method of Effective Charges. Phys. Rev. D 1989, 40, 680. [Google Scholar] [CrossRef]
  139. Dokshitzer, Y.L. Perturbative QCD theory (includes our knowledge of α(s)) - hep-ph/9812252. In Proceedings of the 29th International Conference High-Energy Physics (ICHEP’98), Vancouver, BC, Canada, 23–29 July 1998; Volume 1, pp. 305–324. [Google Scholar]
  140. Prosperi, G.M.; Raciti, M.; Simolo, C. On the running coupling constant in QCD. Prog. Part. Nucl. Phys. 2007, 58, 387–438. [Google Scholar] [CrossRef] [Green Version]
  141. Bjorken, J.D. Applications of the Chiral U(6) x (6) Algebra of Current Densities. Phys. Rev. 1966, 148, 1467–1478. [Google Scholar] [CrossRef]
  142. Bjorken, J.D. Inelastic Scattering of Polarized Leptons from Polarized Nucleons. Phys. Rev. D 1970, 1, 1376–1379. [Google Scholar] [CrossRef] [Green Version]
  143. Brodsky, S.J.; de Téramond, G.F. Light-front hadron dynamics and AdS/CFT correspondence. Phys. Lett. B 2004, 582, 211–221. [Google Scholar] [CrossRef] [Green Version]
  144. Brodsky, S.J.; de Teramond, G.F.; Deur, A. Nonperturbative QCD Coupling and its β-function from Light-Front Holography. Phys. Rev. D 2010, 81, 096010. [Google Scholar] [CrossRef] [Green Version]
  145. Carlson, J.; Jaffe, A.; Wiles, A. (Eds.) The Millennium Prize Problems; American Mathematical Society: Providence, RI, USA, 2006. [Google Scholar]
  146. Eichten, E.; Gottfried, K.; Kinoshita, T.; Lane, K.D.; Yan, T.M. Charmonium: The Model. Phys. Rev. D 1978, 17, 3090, Erratum: Phys. Rev. D 1980, 21, 313. [Google Scholar] [CrossRef]
  147. Richardson, J.L. The Heavy Quark Potential and the Upsilon, J/psi Systems. Phys. Lett. B 1979, 82, 272–274. [Google Scholar] [CrossRef]
  148. Buchmuller, W.; Tye, S.H.H. Quarkonia and Quantum Chromodynamics. Phys. Rev. D 1981, 24, 132. [Google Scholar] [CrossRef] [Green Version]
  149. Godfrey, S.; Isgur, N. Mesons in a Relativized Quark Model with Chromodynamics. Phys. Rev. D 1985, 32, 189–231. [Google Scholar] [CrossRef]
  150. Lucha, W.; Schöberl, F.F.; Gromes, D. Bound states of quarks. Phys. Rept. 1991, 200, 127–240. [Google Scholar] [CrossRef]
  151. Binosi, D.; Chang, L.; Qin, S.X.; Papavassiliou, J.; Roberts, C.D. Symmetry preserving truncations of the gap and Bethe–Salpeter equations. Phys. Rev. D 2016, 93, 096010. [Google Scholar] [CrossRef] [Green Version]
  152. Smirnov, A.V.; Smirnov, V.A.; Steinhauser, M. Three-loop static potential. Phys. Rev. Lett. 2010, 104, 112002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  153. Brodsky, S.J.; Deur, A.; de Téramond, G.F.; Dosch, H.G. Light-Front Holography and Superconformal Quantum Mechanics: A New Approach to Hadron Structure and Color Confinement. Int. J. Mod. Phys. Conf. Ser. 2015, 39, 1560081. [Google Scholar] [CrossRef] [Green Version]
  154. Brodsky, S.J.; de Teramond, G.F.; Dosch, H.G.; Erlich, J. Light-Front Holographic QCD and Emerging Confinement. Phys. Rept. 2015, 584, 1–105. [Google Scholar] [CrossRef] [Green Version]
  155. Reinhardt, H.; Burgio, G.; Campagnari, D.; Ebadati, E.; Heffner, J.; Quandt, M.; Vastag, P.; Vogt, H. Hamiltonian approach to QCD in Coulomb gauge - a survey of recent results. Adv. High Energy Phys. 2018, 2018, 2312498. [Google Scholar] [CrossRef]
  156. Hoyer, P. Journey to the Bound States; Springer Briefs in Physics; Springer: Berlin/Heidelberg, Germany, 2021. [Google Scholar]
  157. Bali, G.S.; Neff, H.; Duessel, T.; Lippert, T.; Schilling, K. Observation of string breaking in QCD. Phys. Rev. D 2005, 71, 114513. [Google Scholar] [CrossRef] [Green Version]
  158. Prkacin, Z.; Bali, G.S.; Dussel, T.; Lippert, T.; Neff, H.; Schilling, K. Anatomy of string breaking in QCD. High Energy Phys. Lattice 2006, LAT2005, 308. [Google Scholar]
  159. Chang, L.; Cloet, I.C.; El-Bennich, B.; Klähn, T.; Roberts, C.D. Exploring the light quark interaction. Chin. Phys. C 2009, 33, 1189–1196. [Google Scholar]
  160. Munczek, H.J.; Nemirovsky, A.M. The Ground State q q ¯ Mass Spectrum in QCD. Phys. Rev. D 1983, 28, 181–186. [Google Scholar] [CrossRef]
  161. Stingl, M. Propagation Properties and Condensate Formation of the Confined Yang–Mills Field. Phys. Rev. D 1986, 34, 3863–3881, [Erratum: Phys. Rev. D 36, 651 (1987)]. [Google Scholar] [CrossRef] [PubMed]
  162. Zwanziger, D. Local and Renormalizable Action From the Gribov Horizon. Nucl. Phys. B 1989, 323, 513–544. [Google Scholar] [CrossRef]
  163. Burden, C.J.; Roberts, C.D.; Williams, A.G. Singularity structure of a model quark propagator. Phys. Lett. B 1992, 285, 347–353. [Google Scholar] [CrossRef]
  164. Gribov, V.N. The theory of quark confinement. Eur. Phys. J. C 1999, 10, 91–105. [Google Scholar] [CrossRef]
  165. Stingl, M. A Systematic extended iterative solution for quantum chromodynamics. Z. Phys. A 1996, 353, 423–445. [Google Scholar] [CrossRef] [Green Version]
  166. Dudal, D.; Gracey, J.A.; Sorella, S.P.; Vandersickel, N.; Verschelde, H. A Refinement of the Gribov–Zwanziger approach in the Landau gauge: Infrared propagators in harmony with the lattice results. Phys. Rev. D 2008, 78, 065047. [Google Scholar] [CrossRef] [Green Version]
  167. Qin, S.X.; Rischke, D.H. Quark Spectral Function and Deconfinement at Nonzero Temperature. Phys. Rev. D 2013, 88, 056007. [Google Scholar] [CrossRef] [Green Version]
  168. Lucha, W.; Schöberl, F.F. Analytic Bethe–Salpeter Description of the Lightest Pseudoscalar Mesons. Phys. Rev. D 2016, 93, 056006. [Google Scholar] [CrossRef] [Green Version]
  169. Dudal, D.; Oliveira, O.; Roelfs, M.; Silva, P. Spectral representation of lattice gluon and ghost propagators at zero temperature. Nucl. Phys. B 2020, 952, 114912. [Google Scholar] [CrossRef]
  170. Fischer, C.S.; Huber, M.Q. Landau gauge Yang–Mills propagators in the complex momentum plane. Phys. Rev. D 2020, 102, 094005. [Google Scholar] [CrossRef]
  171. Osterwalder, K.; Schrader, R. Axioms for Euclidean Green’s Functions. Commun. Math. Phys. 1973, 31, 83–112. [Google Scholar] [CrossRef]
  172. Osterwalder, K.; Schrader, R. Axioms for Euclidean Green’s Functions. 2. Commun. Math. Phys. 1975, 42, 281. [Google Scholar] [CrossRef]
  173. Dedushenko, M. Snowmass White Paper: The Quest to Define QFT. arXiv 2022, arXiv:2203.08053. [Google Scholar]
  174. Stodolna, A.S.; Rouzée, A.; Lépine, F.; Cohen, S.; Robicheaux, F.; Gijsbertsen, A.; Jungmann, J.H.; Bordas, C.; Vrakking, M.J.J. Hydrogen Atoms under Magnification: Direct Observation of the Nodal Structure of Stark States. Phys. Rev. Lett. 2013, 110, 213001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  175. Bhagwat, M.; Pichowsky, M.; Tandy, P.C. Confinement phenomenology in the Bethe–Salpeter equation. Phys. Rev. D 2003, 67, 054019. [Google Scholar] [CrossRef] [Green Version]
  176. Field, R.D.; Feynman, R.P. A Parametrization of the Properties of Quark Jets. Nucl. Phys. B 1978, 136, 1–76. [Google Scholar] [CrossRef]
  177. Gribov, V.N.; Lipatov, L.N. Deep inelastic electron scattering in perturbation theory. Phys. Lett. B 1971, 37, 78–80. [Google Scholar] [CrossRef]
  178. Brodsky, S.J.; Deshpande, A.L.; Gao, H.; McKeown, R.D.; Meyer, C.A.; Meziani, Z.E.; Milner, R.G.; Qiu, J.; Richards, D.G.; Roberts, C.D. QCD and Hadron Physics. arXiv 2015, arXiv:1502.05728. [Google Scholar]
  179. Denisov, O. Letter of Intent (Draft 2.0): A New QCD facility at the M2 beam line of the CERN SPS. arXiv 2018, arXiv:1808.00848. [Google Scholar]
  180. Aguilar, A.C.; Ahmed, Z.; Aidala, C.; Ali, S.; Andrieux, V.; Arrington, J.; Bashir, A.; Berdnikov, V.; Binosi, D.; Chang, L.; et al. Pion and Kaon Structure at the Electron-Ion Collider. Eur. Phys. J. A 2019, 55, 190. [Google Scholar] [CrossRef]
  181. Brodsky, S.J.; Burkert, V.D.; Carman, D.S.; Chen, J.P.; Cui, Z.-F.; Döring, M.; Dosch, H.G.; Draayer, J.; Elouadrhiri, L.; Glazier, D.I.; et al. Strong QCD from Hadron Structure Experiments. Int. J. Mod. Phys. E 2020, 29, 2030006. [Google Scholar] [CrossRef]
  182. Chen, X.; Guo, F.K.; Roberts, C.D.; Wang, R. Selected Science Opportunities for the EicC. Few Body Syst. 2020, 61, 43. [Google Scholar] [CrossRef]
  183. Anderle, D.P.; Bertone, V.; Cao, X.; Chang, L.; Chang, N.; Chen, G.; Chen, X.; Chen, Z.; Cui, Z.; Dai, L.; et al. Electron-ion collider in China. Front. Phys. 2021, 16, 64701. [Google Scholar] [CrossRef]
  184. Arrington, J.; Gayoso, C.A.; Barry, P.C.; Berdnikov, V.; Binosi, D.; Chang, L.; Diefenthaler, M.; Ding, M.; Ent, R.; Frederico, T.; et al. Revealing the structure of light pseudoscalar mesons at the electron–ion collider. J. Phys. G 2021, 48, 075106. [Google Scholar] [CrossRef]
  185. Aoki, K.; Fujioka, H.; Gogami, T.; Hidaka, Y.; Hiyama, E.; Honda, R.; Hosaka, A.; Ichikawa, Y.; Ieiri, M.; Isaka, M.; et al. Extension of the J-PARC Hadron Experimental Facility: Third White Paper. arXiv 2021, arXiv:2110.04462. [Google Scholar]
  186. Quintans, C. The New AMBER Experiment at the CERN SPS. Few Body Syst. 2022, 63, 72. [Google Scholar] [CrossRef]
  187. Capstick, S.; Roberts, W. Quark models of baryon masses and decays. Prog. Part. Nucl. Phys. 2000, 45, S241–S331. [Google Scholar] [CrossRef] [Green Version]
  188. Giannini, M.M.; Santopinto, E. The hypercentral Constituent Quark Model and its application to baryon properties. Chin. J. Phys. 2015, 53, 020301. [Google Scholar]
  189. Plessas, W. The constituent-quark model — Nowadays. Int. J. Mod. Phys. A 2015, 30, 1530013. [Google Scholar] [CrossRef]
  190. Dudek, J.J.; Edwards, R.G.; Peardon, M.J.; Richards, D.G.; Thomas, C.E. Toward the excited meson spectrum of dynamical QCD. Phys. Rev. D 2010, 82, 034508. [Google Scholar] [CrossRef] [Green Version]
  191. Edwards, R.G.; Dudek, J.J.; Richards, D.G.; Wallace, S.J. Excited state baryon spectroscopy from lattice QCD. Phys. Rev. D 2011, 84, 074508. [Google Scholar] [CrossRef] [Green Version]
  192. Liu, L.; Moir, G.; Peardon, M.; Ryan, S.M.; Thomas, C.E.; Vilaseca, P.; Dudek, J.J.; Edwards, R.G.; Joo, B.; Richards, D.G. Excited and exotic charmonium spectroscopy from lattice QCD. J. High Energy Phys. 2012, 7, 126. [Google Scholar] [CrossRef] [Green Version]
  193. Dudek, J.J.; Edwards, R.G. Hybrid Baryons in QCD. Phys. Rev. D 2012, 85, 054016. [Google Scholar] [CrossRef] [Green Version]
  194. Ryan, S.M.; Wilson, D.J. Excited and exotic bottomonium spectroscopy from lattice QCD. J. High Energy Phys. 2021, 02, 214. [Google Scholar] [CrossRef]
  195. Woss, A.J.; Dudek, J.J.; Edwards, R.G.; Thomas, C.E.; Wilson, D.J. Decays of an exotic 1-+ hybrid meson resonance in QCD. Phys. Rev. D 2021, 103, 054502. [Google Scholar] [CrossRef]
  196. Munczek, H.J. Dynamical chiral symmetry breaking, Goldstone’s theorem and the consistency of the Schwinger-Dyson and Bethe–Salpeter Equations. Phys. Rev. D 1995, 52, 4736–4740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  197. Bender, A.; Roberts, C.D.; von Smekal, L. Goldstone Theorem and Diquark Confinement Beyond Rainbow- Ladder Approximation. Phys. Lett. B 1996, 380, 7–12. [Google Scholar] [CrossRef]
  198. Maris, P.; Roberts, C.D.; Tandy, P.C. Pion mass and decay constant. Phys. Lett. B 1998, 420, 267–273. [Google Scholar] [CrossRef] [Green Version]
  199. Maris, P.; Roberts, C.D. π and K meson Bethe–Salpeter amplitudes. Phys. Rev. C 1997, 56, 3369–3383. [Google Scholar] [CrossRef] [Green Version]
  200. Chen, M.; Ding, M.; Chang, L.; Roberts, C.D. Mass-dependence of pseudoscalar meson elastic form factors. Phys. Rev. D 2018, 98, 091505. [Google Scholar] [CrossRef] [Green Version]
  201. Ding, M.; Raya, K.; Bashir, A.; Binosi, D.; Chang, L.; Chen, M.; Roberts, C.D. γ*γη,η transition form factors. Phys. Rev. D 2019, 99, 014014. [Google Scholar] [CrossRef] [Green Version]
  202. Eichmann, G.; Fischer, C.S.; Williams, R. Kaon-box contribution to the anomalous magnetic moment of the muon. Phys. Rev. D 2020, 101, 054015. [Google Scholar] [CrossRef] [Green Version]
  203. Xu, Y.Z.; Binosi, D.; Cui, Z.F.; Li, B.L.; Roberts, C.D.; Xu, S.S.; Zong, H.S. Elastic electromagnetic form factors of vector mesons. Phys. Rev. D 2019, 100, 114038. [Google Scholar] [CrossRef] [Green Version]
  204. Xu, Y.Z.; Chen, S.; Yao, Z.Q.; Binosi, D.; Cui, Z.F.; Roberts, C.D. Vector-meson production and vector meson dominance. Eur. Phys. J. C 2021, 81, 895. [Google Scholar] [CrossRef]
  205. Yao, Z.Q.; Binosi, D.; Cui, Z.F.; Roberts, C.D.; Xu, S.S.; Zong, H.S. Semileptonic decays of D(s) mesons. Phys. Rev. D 2020, 102, 014007. [Google Scholar] [CrossRef]
  206. Yao, Z.Q.; Binosi, D.; Cui, Z.F.; Roberts, C.D. Semileptonic Bcηc,J/ψ transitions. Phys. Lett. B 2021, 818, 136344. [Google Scholar] [CrossRef]
  207. Yao, Z.Q.; Binosi, D.; Cui, Z.F.; Roberts, C.D. Semileptonic transitions: B(s)π(K); DsK; Dπ,K; and Kπ. Phys. Lett. B 2022, 824, 136793. [Google Scholar] [CrossRef]
  208. Höll, A.; Krassnigg, A.; Roberts, C.D. Pseudoscalar meson radial excitations. Phys. Rev. C 2004, 70, 042203. [Google Scholar] [CrossRef] [Green Version]
  209. Höll, A.; Krassnigg, A.; Maris, P.; Roberts, C.D.; Wright, S.V. Electromagnetic properties of ground and excited state pseudoscalar mesons. Phys. Rev. C 2005, 71, 065204. [Google Scholar] [CrossRef] [Green Version]
  210. Fischer, C.S.; Williams, R. Probing the gluon self-interaction in light mesons. Phys. Rev. Lett. 2009, 103, 122001. [Google Scholar] [CrossRef] [Green Version]
  211. Krassnigg, A. Survey of J=0,1 mesons in a Bethe–Salpeter approach. Phys. Rev. D 2009, 80, 114010. [Google Scholar] [CrossRef]
  212. Qin, S.X.; Chang, L.; Liu, Y.X.; Roberts, C.D.; Wilson, D.J. Interaction model for the gap equation. Phys. Rev. C 2011, 84, 042202. [Google Scholar] [CrossRef] [Green Version]
  213. Qin, S.X.; Chang, L.; Liu, Y.X.; Roberts, C.D.; Wilson, D.J. Investigation of rainbow-ladder truncation for excited and exotic mesons. Phys. Rev. C 2012, 85, 035202. [Google Scholar] [CrossRef] [Green Version]
  214. Blank, M.; Krassnigg, A. Bottomonium in a Bethe–Salpeter-equation study. Phys. Rev. D 2011, 84, 096014. [Google Scholar]
  215. Hilger, T.; Popovici, C.; Gomez-Rocha, M.; Krassnigg, A. Spectra of heavy quarkonia in a Bethe–Salpeter-equation approach. Phys. Rev. D 2015, 91, 034013. [Google Scholar] [CrossRef] [Green Version]
  216. Fischer, C.S.; Kubrak, S.; Williams, R. Spectra of heavy mesons in the Bethe–Salpeter approach. Eur. Phys. J. A 2015, 51, 10. [Google Scholar] [CrossRef] [Green Version]
  217. Eichmann, G.; Fischer, C.S.; Sanchis-Alepuz, H. Light baryons and their excitations. Phys. Rev. D 2016, 94, 094033. [Google Scholar] [CrossRef] [Green Version]
  218. Qin, S.X.; Roberts, C.D.; Schmidt, S.M. Spectrum of light- and heavy-baryons. Few Body Syst. 2019, 60, 26. [Google Scholar] [CrossRef]
  219. Chang, L.; Roberts, C.D. Tracing masses of ground state light quark mesons. Phys. Rev. C 2012, 85, 052201. [Google Scholar] [CrossRef] [Green Version]
  220. Williams, R.; Fischer, C.S.; Heupel, W. Light mesons in QCD and unquenching effects from the 3PI effective action. Phys. Rev. D 2016, 93, 034026. [Google Scholar] [CrossRef] [Green Version]
  221. Qin, S.X.; Roberts, C.D. Resolving the Bethe–Salpeter kernel. Chin. Phys. Lett. Express 2021, 38, 071201. [Google Scholar] [CrossRef]
  222. Xu, Z.N.; Yao, Z.Q.; Qin, S.X.; Cui, Z.F.; Roberts, C.D. Bethe–Salpeter kernel and properties of strange-quark mesons. arXiv 2022, arXiv:2208.13903. [Google Scholar]
  223. Qin, S.X.; Roberts, C.D.; Schmidt, S.M. Ward–Green–Takahashi identities and the axial-vector vertex. Phys. Lett. B 2014, 733, 202–208. [Google Scholar] [CrossRef] [Green Version]
  224. Bhagwat, M.S.; Höll, A.; Krassnigg, A.; Roberts, C.D.; Tandy, P.C. Aspects and consequences of a dressed quark-gluon vertex. Phys. Rev. C 2004, 70, 035205. [Google Scholar] [CrossRef] [Green Version]
  225. Chang, L.; Liu, Y.X.; Roberts, C.D.; Shi, Y.M.; Sun, W.M.; Zong, H.S. Chiral susceptibility and the scalar Ward identity. Phys. Rev. C 2009, 79, 035209. [Google Scholar] [CrossRef] [Green Version]
  226. Singh, J. Anomalous magnetic moment of light quarks and dynamical symmetry breaking. Phys. Rev. D 1985, 31, 1097–1108. [Google Scholar] [CrossRef]
  227. Bicudo, P.J.A.; Ribeiro, J.E.F.T.; Fernandes, R. The anomalous magnetic moment of quarks. Phys. Rev. C 1999, 59, 1107–1112. [Google Scholar] [CrossRef] [Green Version]
  228. Chang, L.; Liu, Y.X.; Roberts, C.D. Dressed-quark anomalous magnetic moments. Phys. Rev. Lett. 2011, 106, 072001. [Google Scholar] [CrossRef] [Green Version]
  229. Bashir, A.; Bermúdez, R.; Chang, L.; Roberts, C.D. Dynamical chiral symmetry breaking and the fermion–gauge-boson vertex. Phys. Rev. C 2012, 85, 045205. [Google Scholar] [CrossRef]
  230. Kızılersü, A.; Oliveira, O.; Silva, P.J.; Skullerud, J.I.; Sternbeck, A. Quark-gluon vertex from Nf=2 lattice QCD. Phys. Rev. D 2021, 103, 114515. [Google Scholar] [CrossRef]
  231. Williams, R. The quark-gluon vertex in Landau gauge bound state studies. Eur. Phys. J. A 2015, 51, 57. [Google Scholar] [CrossRef] [Green Version]
  232. Yin, P.L.; Cui, Z.F.; Roberts, C.D.; Segovia, J. Masses of positive- and negative-parity hadron ground states, including those with heavy quarks. Eur. Phys. J. C 2021, 81, 327. [Google Scholar] [CrossRef]
  233. Xu, Z.N.; Cui, Z.F.; Roberts, C.D.; Xu, C. Heavy + light pseudoscalar meson semileptonic transitions. Eur. Phys. J. C 2021, 81, 1105. [Google Scholar] [CrossRef]
  234. Gutiérrez-Guerrero, L.X.; Paredes-Torres, G.; Bashir, A. Mesons and baryons: Parity partners. Phys. Rev. D 2021, 104, 094013. [Google Scholar] [CrossRef]
  235. Höll, A.; Maris, P.; Roberts, C.D.; Wright, S.V. Schwinger functions and light quark bound states, and sigma terms. Nucl. Phys. Proc. Suppl. 2006, 161, 87–94. [Google Scholar] [CrossRef] [Green Version]
  236. Santowsky, N.; Eichmann, G.; Fischer, C.S.; Wallbott, P.C.; Williams, R. The σ-meson: Four-quark vs. two-quark components and decay width in a Bethe–Salpeter approach. Phys. Rev. D 2020, 102, 056014. [Google Scholar] [CrossRef]
  237. Nakanishi, N. A General survey of the theory of the Bethe–Salpeter equation. Prog. Theor. Phys. Suppl. 1969, 43, 1–81. [Google Scholar] [CrossRef] [Green Version]
  238. Cui, Z.F.; Binosi, D.; Roberts, C.D.; Schmidt, S.M. Hadron and light nucleus radii from electron scattering. Chin. Phys. C 2022, in press. [Google Scholar] [CrossRef]
  239. Binosi, D.; Chang, L.; Ding, M.; Gao, F.; Papavassiliou, J.; Roberts, C.D. Distribution Amplitudes of Heavy-Light Mesons. Phys. Lett. B 2019, 790, 257–262. [Google Scholar] [CrossRef]
  240. McNeile, C.; Michael, C. The decay constant of the first excited pion from lattice QCD. Phys. Lett. B 2006, 642, 244–247. [Google Scholar] [CrossRef] [Green Version]
  241. Ballon-Bayona, A.; Krein, G.; Miller, C. Decay constants of the pion and its excitations in holographic QCD. Phys. Rev. D 2015, 91, 065024. [Google Scholar] [CrossRef] [Green Version]
  242. Chen, M.; Chang, L. A pattern for the flavor dependent quark-antiquark interaction. Chin. Phys. C 2019, 43, 114103. [Google Scholar] [CrossRef] [Green Version]
  243. Qin, P.; Qin, S.X.; Liu, Y.X. Heavy-light mesons beyond the ladder approximation. Phys. Rev. D 2020, 101, 114014. [Google Scholar] [CrossRef]
  244. Burden, C.J.; Pichowsky, M.A. JPC exotic mesons from the Bethe–Salpeter equation. Few Body Syst. 2002, 32, 119–126. [Google Scholar] [CrossRef]
  245. Hilger, T.; Gomez-Rocha, M.; Krassnigg, A. Masses of JPC=1-+ exotic quarkonia in a Bethe–Salpeter-equation approach. Phys. Rev. D 2015, 91, 114004. [Google Scholar] [CrossRef] [Green Version]
  246. Xu, S.S.; Cui, Z.F.; Chang, L.; Papavassiliou, J.; Roberts, C.D.; Zong, H.S. New perspective on hybrid mesons. Eur. Phys. J. A (Lett.) 2019, 55, 113. [Google Scholar] [CrossRef] [Green Version]
  247. Meyers, J.; Swanson, E.S. Spin Zero Glueballs in the Bethe–Salpeter Formalism. Phys. Rev. D 2013, 87, 036009. [Google Scholar] [CrossRef] [Green Version]
  248. Souza, E.V.; Narciso Ferreira, M.; Aguilar, A.C.; Papavassiliou, J.; Roberts, C.D.; Xu, S.S. Pseudoscalar glueball mass: A window on three-gluon interactions. Eur. Phys. J. A (Lett.) 2020, 56, 25. [Google Scholar] [CrossRef] [Green Version]
  249. Kaptari, L.P.; Kämpfer, B. Mass spectrum of pseudo-scalar glueballs from a Bethe–Salpeter approach with the rainbow-ladder truncation. Few Body Syst. 2020, 61, 28. [Google Scholar] [CrossRef]
  250. Huber, M.Q.; Fischer, C.S.; Sanchis-Alepuz, H. Higher spin glueballs from functional methods. Eur. Phys. J. C 2021, 81, 1083. [Google Scholar] [CrossRef]
  251. Ablikim, M.; Achasov, M.N.; Adlarson, P.; Ahmed, S.; Albrecht, M.; Alekseev, M.; Amoroso, A.; An, F.F.; An, Q.; Bai, Y.; et al. Future Physics Programme of BESIII. Chin. Phys. C 2020, 44, 040001. [Google Scholar] [CrossRef]
  252. Khalek, R.A.; Accardi, A.; Adam, J.; Adamiak, D.; Akers, W.; Albaladejo, M.; Blin, A.H. Science Requirements and Detector Concepts for the Electron-Ion Collider: EIC Yellow Report. Nucl. Phys. A 2022, 1026, 122447. [Google Scholar] [CrossRef]
  253. Pauli, P. Accessing glue through photoproduction measurements at GlueX. Rev. Mex. Fis. Suppl. 2022, 3, 0308002. [Google Scholar] [CrossRef]
  254. Eichmann, G.; Alkofer, R.; Krassnigg, A.; Nicmorus, D. Nucleon mass from a covariant three-quark Faddeev equation. Phys. Rev. Lett. 2010, 104, 201601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  255. Barabanov, M.; Bedolla, M.; Brooks, W.; Cates, G.; Chen, C.; Chen, Y.; Cisbani, E.; Ding, M.; Eichmann, G.; Ent, R.; et al. Diquark Correlations in Hadron Physics: Origin, Impact and Evidence. Prog. Part. Nucl. Phys. 2021, 116, 103835. [Google Scholar] [CrossRef]
  256. Cahill, R.T.; Roberts, C.D.; Praschifka, J. Baryon structure and QCD. Austral. J. Phys. 1989, 42, 129–145. [Google Scholar] [CrossRef]
  257. Burden, C.J.; Cahill, R.T.; Praschifka, J. Baryon Structure and QCD: Nucleon Calculations. Austral. J. Phys. 1989, 42, 147–159. [Google Scholar] [CrossRef] [Green Version]
  258. Reinhardt, H. Hadronization of Quark Flavor Dynamics. Phys. Lett. B 1990, 244, 316–326. [Google Scholar] [CrossRef]
  259. Efimov, G.V.; Ivanov, M.A.; Lyubovitskij, V.E. Quark—Diquark approximation of the three quark structure of baryons in the quark confinement model. Z. Phys. C 1990, 47, 583–594. [Google Scholar] [CrossRef]
  260. Gutiérrez-Guerrero, L.X.; Bashir, A.; Cloet, I.C.; Roberts, C.D. Pion form factor from a contact interaction. Phys. Rev. C 2010, 81, 065202. [Google Scholar] [CrossRef] [Green Version]
  261. Gutiérrez-Guerrero, L.X.; Bashir, A. Contact Interaction Studies of Hadron Observables. Particles 2022, 5, 69–138. [Google Scholar]
  262. Brown, Z.S.; Detmold, W.; Meinel, S.; Orginos, K. Charmed bottom baryon spectroscopy from lattice QCD. Phys. Rev. D 2014, 90, 094507. [Google Scholar] [CrossRef] [Green Version]
  263. Mathur, N.; Padmanath, M.; Mondal, S. Precise predictions of charmed-bottom hadrons from lattice QCD. Phys. Rev. Lett. 2018, 121, 202002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  264. Bahtiyar, H.; Can, K.U.; Erkol, G.; Gubler, P.; Oka, M.; Takahashi, T.T. Charmed baryon spectrum from lattice QCD near the physical point. Phys. Rev. D 2020, 102, 054513. [Google Scholar] [CrossRef]
  265. Chen, C.; El-Bennich, B.; Roberts, C.D.; Schmidt, S.M.; Segovia, J.; Wan, S. Structure of the nucleon’s low-lying excitations. Phys. Rev. D 2018, 97, 034016. [Google Scholar] [CrossRef] [Green Version]
  266. Chen, C.; Krein, G.I.; Roberts, C.D.; Schmidt, S.M.; Segovia, J. Spectrum and structure of octet and decuplet baryons and their positive-parity excitations. Phys. Rev. D 2019, 100, 054009. [Google Scholar] [CrossRef] [Green Version]
  267. Liu, L.; Chen, C.; Lu, Y.; Roberts, C.D.; Segovia, J. Composition of low-lying J = 3 2 ± Δ-baryons. Phys. Rev. D 2022, 105, 114047. [Google Scholar] [CrossRef]
  268. Liu, L.; Chen, C.; Roberts, C.D. Wave functions of (I, JP) = ( 1 2 , 3 2 ) baryons. arXiv 2022, arXiv:2208.12353. [Google Scholar] [CrossRef]
  269. Crede, V.; Roberts, W. Progress towards understanding baryon resonances. Rept. Prog. Phys. 2013, 76, 076301. [Google Scholar] [CrossRef] [Green Version]
  270. Alkofer, R.; Höll, A.; Kloker, M.; Krassnigg, A.; Roberts, C.D. On nucleon electromagnetic form-factors. Few Body Syst. 2005, 37, 1–31. [Google Scholar] [CrossRef]
  271. Segovia, J.; Cloet, I.C.; Roberts, C.D.; Schmidt, S.M. Nucleon and Δ elastic and transition form factors. Few Body Syst. 2014, 55, 1185–1222. [Google Scholar] [CrossRef] [Green Version]
  272. Chen, C.; Lu, Y.; Binosi, D.; Roberts, C.D.; Rodríguez-Quintero, J.; Segovia, J. Nucleon-to-Roper electromagnetic transition form factors at large Q2. Phys. Rev. D 2019, 99, 034013. [Google Scholar] [CrossRef] [Green Version]
  273. Lu, Y.; Chen, C.; Cui, Z.F.; Roberts, C.D.; Schmidt, S.M.; Segovia, J.; Zong, H.S. Transition form factors: γ*+p→Δ(1232), Δ(1600). Phys. Rev. D 2019, 100, 034001. [Google Scholar] [CrossRef] [Green Version]
  274. Chen, C.; Fischer, C.S.; Roberts, C.D.; Segovia, J. Form Factors of the Nucleon Axial Current. Phys. Lett. B 2021, 815, 136150. [Google Scholar] [CrossRef]
  275. Chen, C.; Fischer, C.S.; Roberts, C.D.; Segovia, J. Nucleon axial-vector and pseudoscalar form factors and PCAC relations. Phys. Rev. D 2022, 105, 094022. [Google Scholar] [CrossRef]
  276. Chen, C.; Roberts, C.D. Nucleon axial form factor at large momentum transfers. Eur. Phys. J. A 2022, 58, 206. [Google Scholar] [CrossRef]
  277. Klempt, E.; Metsch, B.C. Multiplet classification of light quark baryons. Eur. Phys. J. A 2012, 48, 127. [Google Scholar] [CrossRef]
  278. Weinberg, S. Precise relations between the spectra of vector and axial vector mesons. Phys. Rev. Lett. 1967, 18, 507–509. [Google Scholar] [CrossRef] [Green Version]
  279. Lane, K.D. Asymptotic Freedom and Goldstone Realization of Chiral Symmetry. Phys. Rev. D 1974, 10, 2605. [Google Scholar] [CrossRef]
  280. Politzer, H.D. Effective Quark Masses in the Chiral Limit. Nucl. Phys. B 1976, 117, 397. [Google Scholar] [CrossRef]
  281. Pagels, H. Dynamical Chiral Symmetry Breaking in Quantum Chromodynamics. Phys. Rev. D 1979, 19, 3080. [Google Scholar] [CrossRef]
  282. Higashijima, K. Dynamical Chiral Symmetry Breaking. Phys. Rev. D 1984, 29, 1228. [Google Scholar] [CrossRef] [Green Version]
  283. Pagels, H.; Stokar, S. The Pion Decay Constant, Electromagnetic Form-Factor and Quark Electromagnetic Selfenergy in QCD. Phys. Rev. D 1979, 20, 2947. [Google Scholar] [CrossRef]
  284. Cahill, R.T.; Roberts, C.D. Soliton Bag Models of Hadrons from QCD. Phys. Rev. D 1985, 32, 2419. [Google Scholar] [CrossRef] [PubMed]
  285. Munczek, H.J. Infrared Effects on the Axial Gauge Quark Propagator. Phys. Lett. B 1986, 175, 215–218. [Google Scholar] [CrossRef]
  286. Efimov, G.V.; Ivanov, M.A. Confinement and quark structure of light hadrons. Int. J. Mod. Phys. A 1989, 4, 2031–2060. [Google Scholar] [CrossRef]
  287. Hecht, M.B.; Oettel, M.; Roberts, C.D.; Schmidt, S.M.; Tandy, P.C.; Thomas, A.W. Nucleon mass and pion loops. Phys. Rev. C 2002, 65, 055204. [Google Scholar] [CrossRef] [Green Version]
  288. Sanchis-Alepuz, H.; Fischer, C.S.; Kubrak, S. Pion cloud effects on baryon masses. Phys. Lett. B 2014, 733, 151–157. [Google Scholar] [CrossRef] [Green Version]
  289. Julia-Diaz, B.; Lee, T.S.H.; Matsuyama, A.; Sato, T. Dynamical coupled channel model of pi N scattering in the W≤2-GeV nucleon resonance region. Phys. Rev. C 2007, 76, 065201. [Google Scholar] [CrossRef] [Green Version]
  290. Suzuki, N.; Julia-Diaz, B.; Kamano, H.; Lee, T.S.H.; Matsuyama, A.; Sato, T. Disentangling the Dynamical Origin of P-11 Nucleon Resonances. Phys. Rev. Lett. 2010, 104, 042302. [Google Scholar] [CrossRef]
  291. Rönchen, D.; Döring, M.; Huang, F.; Haberzettl, H.; Haidenbauer, J.; Hanhart, C.; Krewald, S.; Meissner, U.G.; Nakayama, K. Coupled-channel dynamics in the reactions πNπN, ηN, KΛ, KΣ. Eur. Phys. J. A 2013, 49, 44. [Google Scholar] [CrossRef] [Green Version]
  292. Kamano, H.; Nakamura, S.X.; Lee, T.S.H.; Sato, T. Nucleon resonances within a dynamical coupled channels model of πN and γN reactions. Phys. Rev. C 2013, 88, 035209. [Google Scholar] [CrossRef] [Green Version]
  293. García-Tecocoatzi, H.; Bijker, R.; Ferretti, J.; Santopinto, E. Self-energies of octet and decuplet baryons due to the coupling to the baryon-meson continuum. Eur. Phys. J. A 2017, 53, 115. [Google Scholar] [CrossRef]
  294. Eichmann, G.; Alkofer, R.; Cloet, I.C.; Krassnigg, A.; Roberts, C.D. Perspective on rainbow-ladder truncation. Phys. Rev. C 2008, 77, 042202. [Google Scholar] [CrossRef] [Green Version]
  295. Eichmann, G.; Cloet, I.C.; Alkofer, R.; Krassnigg, A.; Roberts, C.D. Toward unifying the description of meson and baryon properties. Phys. Rev. C 2009, 79, 012202. [Google Scholar] [CrossRef]
  296. Roberts, H.L.L.; Chang, L.; Cloet, I.C.; Roberts, C.D. Masses of ground and excited state hadrons. Few Body Syst. 2011, 51, 1–25. [Google Scholar] [CrossRef]
  297. Isupov, E.L.; Burkert, V.D.; Carman, D.S.; Gothe, R.W.; Hicks, K.; Ishkhanov, B.S.; Mokeev, V.I.; Adhikari, K.P.; Adhikari, S.; Adikaram, D.; et al. Measurements of epeπ+πp Cross Sections with CLAS at 1.40 Gev < W < 2.0 GeV and 2.0 GeV2 < Q2 < 5.0 GeV2. Phys. Rev. C 2017, 96, 025209. [Google Scholar]
  298. Trivedi, A. Measurement of New Observables from the π+π p Electroproduction Off the Proton. Few Body Syst. 2019, 60, 5. [Google Scholar] [CrossRef]
  299. Raya, K.; Gutiérrez-Guerrero, L.X.; Bashir, A.; Chang, L.; Cui, Z.F.; Lu, Y.; Roberts, C.D.; Segovia, J. Dynamical diquarks in the γ(*)pN(1535) 1 2 transition. Eur. Phys. J. A 2021, 57, 266. [Google Scholar] [CrossRef]
  300. Aznauryan, I.G.; Burkert, V.D.; Biselli, A.S.; Egiyan, H.; Joo, K.; Kim, W.; Park, K.; Smith, L.C.; Ungaro, M.; Adhikari, K.P.; et al. Electroexcitation of nucleon resonances from CLAS data on single pion electroproduction. Phys. Rev. C 2009, 80, 055203. [Google Scholar] [CrossRef]
  301. Mokeev, V.I.; Burkert, V.D.; Carman, D.S.; Elouadrhiri, L.; Fedotov, G.V.; Golovatch, E.N.; Gothe, R.W.; Hicks, K.; Ishkhanov, B.S.; Isupov, E.L.; et al. New Results from the Studies of the N(1440)1/2+, N(1520)3/2, and Δ(1620)1/2 Resonances in Exclusive epepπ+π Electroproduction with the CLAS Detector. Phys. Rev. C 2016, 93, 025206. [Google Scholar] [CrossRef]
  302. Mokeev, V.I. Two Pion Photo- and Electroproduction with CLAS. EPJ Web Conf. 2020, 241, 03003. [Google Scholar] [CrossRef]
  303. Carman, D.; Joo, K.; Mokeev, V. Strong QCD Insights from Excited Nucleon Structure Studies with CLAS and CLAS12. Few Body Syst. 2020, 61, 29. [Google Scholar] [CrossRef]
  304. Mokeev, V.I.; Carman, D.S. New baryon states in exclusive meson photo-/electroproduction with CLAS. Rev. Mex. Fis. Suppl. 2022, 3, 0308024. [Google Scholar] [CrossRef]
  305. Mokeev, V.I.; Carman, D.S. Photo- and Electrocouplings of Nucleon Resonances. Few Body Syst. 2022, 63, 59. [Google Scholar] [CrossRef]
  306. Chang, L.; Cloet, I.C.; Roberts, C.D.; Schmidt, S.M.; Tandy, P.C. Pion electromagnetic form factor at spacelike momenta. Phys. Rev. Lett. 2013, 111, 141802. [Google Scholar] [CrossRef] [Green Version]
  307. Raya, K.; Chang, L.; Bashir, A.; Cobos-Martinez, J.J.; Gutiérrez-Guerrero, L.X.; Roberts, C.D.; Tandy, P.C. Structure of the neutral pion and its electromagnetic transition form factor. Phys. Rev. D 2016, 93, 074017. [Google Scholar] [CrossRef] [Green Version]
  308. Raya, K.; Ding, M.; Bashir, A.; Chang, L.; Roberts, C.D. Partonic structure of neutral pseudoscalars via two photon transition form factors. Phys. Rev. D 2017, 95, 074014. [Google Scholar] [CrossRef] [Green Version]
  309. Gao, F.; Chang, L.; Liu, Y.X.; Roberts, C.D.; Tandy, P.C. Exposing strangeness: Projections for kaon electromagnetic form factors. Phys. Rev. D 2017, 96, 034024. [Google Scholar] [CrossRef] [Green Version]
  310. Lepage, G.P.; Brodsky, S.J. Exclusive Processes in Quantum Chromodynamics: Evolution Equations for Hadronic Wave Functions and the Form-Factors of Mesons. Phys. Lett. B 1979, 87, 359–365. [Google Scholar] [CrossRef]
  311. Efremov, A.V.; Radyushkin, A.V. Factorization and Asymptotical Behavior of Pion Form- Factor in QCD. Phys. Lett. B 1980, 94, 245–250. [Google Scholar] [CrossRef]
  312. Lepage, G.P.; Brodsky, S.J. Exclusive Processes in Perturbative Quantum Chromodynamics. Phys. Rev. D 1980, 22, 2157–2198. [Google Scholar] [CrossRef] [Green Version]
  313. Brodsky, S.J.; Schmidt, I.A.; de Teramond, G.F. Nuclear-bound quarkonium. Phys. Rev. Lett. 1990, 64, 1011. [Google Scholar] [CrossRef] [Green Version]
  314. Tarrús Castellà, J.; Krein, G. Effective field theory for the nucleon-quarkonium interaction. Phys. Rev. D 2018, 98, 014029. [Google Scholar] [CrossRef] [Green Version]
  315. Luke, M.E.; Manohar, A.V.; Savage, M.J. A QCD Calculation of the interaction of quarkonium with nuclei. Phys. Lett. B 1992, 288, 355–359. [Google Scholar] [CrossRef] [Green Version]
  316. Kharzeev, D. Quarkonium interactions in QCD. Proc. Int. Sch. Phys. Fermi 1996, 130, 105–131. [Google Scholar]
  317. Krein, G.; Peixoto, T.C. Femtoscopy of the Origin of the Nucleon Mass. Few Body Syst. 2020, 61, 49. [Google Scholar] [CrossRef]
  318. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Observation of J/ψp Resonances Consistent with Pentaquark States in Λ b 0 J/ψKp Decays. Phys. Rev. Lett. 2015, 115, 072001. [Google Scholar] [CrossRef] [Green Version]
  319. Ali, A.; Amaryan, M.; Anassontzis, E.G.; Austregesilo, A.; Baalouch, M.; Barbosa, F.; Barlow, J.; Barnes, A.; Barriga, E.; Beattie, T.D.; et al. First Measurement of Near-Threshold J/ψ Exclusive Photoproduction off the Proton. Phys. Rev. Lett. 2019, 123, 072001. [Google Scholar] [CrossRef] [Green Version]
  320. Sakurai, J.J. Theory of strong interactions. Annals Phys. 1960, 11, 1–48. [Google Scholar] [CrossRef]
  321. Sakurai, J.J. Currents and Mesons; University of Chicago Press: Chicago, IL, USA, 1969. [Google Scholar]
  322. Fraas, H.; Schildknecht, D. Vector-meson dominance and diffractive electroproduction of vector mesons. Nucl. Phys. B 1969, 14, 543–565. [Google Scholar] [CrossRef]
  323. Sun, P.; Tong, X.B.; Yuan, F. Near threshold heavy quarkonium photoproduction at large momentum transfer. Phys. Rev. D 2022, 105, 054032. [Google Scholar] [CrossRef]
  324. Ivanov, M.A.; Kalinovsky, Y.L.; Roberts, C.D. Survey of heavy-meson observables. Phys. Rev. D 1999, 60, 034018. [Google Scholar] [CrossRef] [Green Version]
  325. Du, M.L.; Baru, V.; Guo, F.K.; Hanhart, C.; Meißner, U.G.; Nefediev, A.; Strakovsky, I. Deciphering the mechanism of near-threshold J/ψ photoproduction. Eur. Phys. J. C 2020, 80, 1053. [Google Scholar] [CrossRef]
  326. Lee, T.S.H.; Sakinah, S.; Oh, Y. Models of J/Ψ photo-production reactions on the nucleon. arXiv 2022, arXiv:2210.02154. [Google Scholar]
  327. Eichmann, G.; Fischer, C.S.; Heupel, W.; Santowsky, N.; Wallbott, P.C. Four-quark states from functional methods. Few Body Syst. 2020, 61, 38. [Google Scholar] [CrossRef]
  328. Pichowsky, M.A.; Lee, T.S.H. Exclusive diffractive processes and the quark substructure of mesons. Phys. Rev. D 1997, 56, 1644–1662. [Google Scholar] [CrossRef] [Green Version]
  329. Shi, C.; Xie, Y.P.; Li, M.; Chen, X.; Zong, H.S. Light front wave functions and diffractive electroproduction of vector mesons. Phys. Rev. D 2021, 104, L091902. [Google Scholar] [CrossRef]
  330. Roberts, C.D. Electromagnetic pion form-factor and neutral pion decay width. Nucl. Phys. A 1996, 605, 475–495. [Google Scholar] [CrossRef] [Green Version]
  331. Miramontes, A.; Bashir, A.; Raya, K.; Roig, P. Pion and Kaon box contribution to a μ HLbL . Phys. Rev. D 2022, 105, 074013. [Google Scholar] [CrossRef]
  332. Rodríguez-Quintero, J. Gravitational form factors of Nambu–Goldstone bosons. In Proceedings of the ECT* Workshop: Revealing Emergent Mass through Studies of Hadron Spectra and Structure, Online, 12–16 September 2022. [Google Scholar]
  333. Eichmann, G. Nucleon electromagnetic form factors from the covariant Faddeev equation. Phys. Rev. D 2011, 84, 014014. [Google Scholar] [CrossRef] [Green Version]
  334. Eichmann, G.; Fischer, C.S. Nucleon axial and pseudoscalar form factors from the covariant Faddeev equation. Eur. Phys. J. A 2012, 48, 9. [Google Scholar] [CrossRef] [Green Version]
  335. Sanchis-Alepuz, H.; Williams, R.; Alkofer, R. Delta and Omega electromagnetic form factors in a three-body covariant Bethe–Salpeter approach. Phys. Rev. D 2013, 87, 095015. [Google Scholar] [CrossRef]
  336. Mokeev, V.I. Insight into EHM from results on electroexcitation of Δ(1600)3/2+ resonance. In Proceedings of the Workshop on Receiving the Emergence of Hadron Mass through AMBER @ CERN-VII, Geneva, Switzerland, 10–13 May 2022. [Google Scholar]
  337. Perdrisat, C.F.; Punjabi, V.; Vanderhaeghen, M. Nucleon electromagnetic form factors. Prog. Part. Nucl. Phys. 2007, 59, 694–764. [Google Scholar] [CrossRef] [Green Version]
  338. Mosel, U. Neutrino Interactions with Nucleons and Nuclei: Importance for Long-Baseline Experiments. Ann. Rev. Nucl. Part. Sci. 2016, 66, 171–195. [Google Scholar] [CrossRef] [Green Version]
  339. Alvarez-Ruso, L.; Athar, M.S.; Barbaro, M.B.; Cherdack, D.; Christy, M.; Coloma, P.; Donnelly, T.; Dytman, S.; de Gouvêa, A.; Hill, R.; et al. NuSTEC White Paper: Status and challenges of neutrino–nucleus scattering. Prog. Part. Nucl. Phys. 2018, 100, 1–68. [Google Scholar] [CrossRef] [Green Version]
  340. Hill, R.J.; Kammel, P.; Marciano, W.J.; Sirlin, A. Nucleon Axial Radius and Muonic Hydrogen — A New Analysis and Review. Rept. Prog. Phys. 2018, 81, 096301. [Google Scholar] [CrossRef] [Green Version]
  341. Gysbers, P.; Hagen, G.; Holt, J.D.; Jansen, G.R.; Morris, T.D.; Navrátil, P.; Papenbrock, T.; Quaglioni, S.; Schwenk, A.; Stroberg, S.R.; et al. Discrepancy between experimental and theoretical β-decay rates resolved from first principles. Nat. Phys. 2019, 15, 428–431. [Google Scholar] [CrossRef] [Green Version]
  342. Lovato, A.; Carlson, J.; Gandolfi, S.; Rocco, N.; Schiavilla, R. Ab initio study of (ν,) and ( ν ¯ ,+) inclusive scattering in 12C: Confronting the MiniBooNE and T2K CCQE data. Phys. Rev. X 2020, 10, 031068. [Google Scholar]
  343. King, G.B.; Andreoli, L.; Pastore, S.; Piarulli, M.; Schiavilla, R.; Wiringa, R.B.; Carlson, J.; Gandolfi, S. Chiral Effective Field Theory Calculations of Weak Transitions in Light Nuclei. Phys. Rev. C 2020, 102, 025501. [Google Scholar] [CrossRef]
  344. Simons, D.; Steinberg, N.; Lovato, A.; Meurice, Y.; Rocco, N.; Wagman, M. Form factor and model dependence in neutrino-nucleus cross section predictions. arXiv 2022, arXiv:2210.02455. [Google Scholar]
  345. Oettel, M.; Pichowsky, M.; von Smekal, L. Current conservation in the covariant quark–diquark model of the nucleon. Eur. Phys. J. A 2000, 8, 251–281. [Google Scholar] [CrossRef] [Green Version]
  346. Meyer, A.S.; Betancourt, M.; Gran, R.; Hill, R.J. Deuterium target data for precision neutrino-nucleus cross sections. Phys. Rev. D 2016, 93, 113015. [Google Scholar] [CrossRef] [Green Version]
  347. Jang, Y.C.; Gupta, R.; Yoon, B.; Bhattacharya, T. Axial Vector Form Factors from Lattice QCD that Satisfy the PCAC Relation. Phys. Rev. Lett. 2020, 124, 072002. [Google Scholar] [CrossRef] [Green Version]
  348. Andreev, V.A.; Banks, T.I.; Carey, R.M.; Case, T.A.; Clayton, S.M.; Crowe, K.M.; Deutsch, J.; Egger, J.; Freedman, S.J.; Ganzha, V.A.; et al. Measurement of Muon Capture on the Proton to 1% Precision and Determination of the Pseudoscalar Coupling gP. Phys. Rev. Lett. 2013, 110, 012504. [Google Scholar] [CrossRef]
  349. Andreev, V.A.; Banks, T.I.; Carey, R.M.; Case, T.A.; Clayton, S.M.; Crowe, K.M.; Deutsch, J.; Egger, J.; Freedman, S.J.; Ganzha, V.A.; et al. Measurement of the Formation Rate of Muonic Hydrogen Molecules. Phys. Rev. C 2015, 91, 055502. [Google Scholar] [CrossRef] [Green Version]
  350. Bernard, V.; Elouadrhiri, L.; Meissner, U.G. Axial structure of the nucleon: Topical Review. J. Phys. G 2002, 28, R1–R35. [Google Scholar] [CrossRef] [Green Version]
  351. Baru, V.; Hanhart, C.; Hoferichter, M.; Kubis, B.; Nogga, A.; Phillips, D.R. Precision calculation of threshold πd scattering, πN scattering lengths, and the GMO sum rule. Nucl. Phys. A 2011, 872, 69–116. [Google Scholar] [CrossRef] [Green Version]
  352. Navarro Pérez, R.; Amaro, J.E.; Ruiz Arriola, E. Precise Determination of Charge Dependent Pion-Nucleon-Nucleon Coupling Constants. Phys. Rev. C 2017, 95, 064001. [Google Scholar] [CrossRef] [Green Version]
  353. Reinert, P.; Krebs, H.; Epelbaum, E. Precision determination of pion–nucleon coupling constants using effective field theory. Phys. Rev. Lett. 2021, 126, 092501. [Google Scholar] [CrossRef]
  354. Bloch, J.C.R.; Roberts, C.D.; Schmidt, S.M. Selected nucleon form-factors and a composite scalar diquark. Phys. Rev. C 2000, 61, 065207. [Google Scholar] [CrossRef] [Green Version]
  355. Cheng, P.; Serna, F.E.; Yao, Z.Q.; Chen, C.; Cui, Z.F.; Roberts, C.D. Contact interaction analysis of octet baryon axial-vector and pseudoscalar form factors. Phys. Rev. D 2022, 106, 054031. [Google Scholar] [CrossRef]
  356. Haidenbauer, J.; Meissner, U.G. The Julich hyperon-nucleon model revisited. Phys. Rev. C 2005, 72, 044005. [Google Scholar] [CrossRef] [Green Version]
  357. Rijken, T.A.; Nagels, M.M.; Yamamoto, Y. Baryon-baryon interactions: Nijmegen extended-soft-core models. Prog. Theor. Phys. Suppl. 2010, 185, 14–71. [Google Scholar] [CrossRef] [Green Version]
  358. Chang, L.; Roberts, C.D.; Schmidt, S.M. Dressed-quarks and the nucleon’s axial charge. Phys. Rev. C 2013, 87, 015203. [Google Scholar] [CrossRef]
  359. Cabibbo, N.; Swallow, E.C.; Winston, R. Semileptonic hyperon decays. Ann. Rev. Nucl. Part. Sci. 2003, 53, 39–75. [Google Scholar] [CrossRef] [Green Version]
  360. Ledwig, T.; Martin Camalich, J.; Geng, L.S.; Vicente Vacas, M.J. Octet-baryon axial-vector charges and SU(3)-breaking effects in the semileptonic hyperon decays. Phys. Rev. D 2014, 90, 054502. [Google Scholar] [CrossRef] [Green Version]
  361. Aidala, C.A.; Bass, S.D.; Hasch, D.; Mallot, G.K. The Spin Structure of the Nucleon. Rev. Mod. Phys. 2013, 85, 655–691. [Google Scholar] [CrossRef] [Green Version]
  362. Deur, A.; Brodsky, S.J.; De Téramond, G.F. The Spin Structure of the Nucleon. Rept. Prog. Phys. 2019, 82, 655. [Google Scholar] [CrossRef] [Green Version]
  363. Del Guerra, A.; Giazotto, A.; Giorgi, M.A.; Stefanini, A.; Botterill, D.R.; Montgomery, H.E.; Norton, P.R.; Matone, G. Threshold π+ electroproduction at high momentum transfer: A determination of the nucleon axial vector form-factor. Nucl. Phys. B 1976, 107, 65–81. [Google Scholar] [CrossRef]
  364. Park, K.; Gothe, R.W.; Adhikari, K.P.; Adikaram, D.; Anghinolfi, M.; Baghdasaryan, H.; Ball, J.; Battaglieri, M.; Batourine, V.; Bedlinskiy, I.; et al. Measurement of the generalized form factors near threshold via γ*p+ at high Q2. Phys. Rev. C 2012, 85, 035208. [Google Scholar] [CrossRef] [Green Version]
  365. Cui, Z.F.; Chen, C.; Binosi, D.; de Soto, F.; Roberts, C.D.; Rodríguez-Quintero, J.; Schmidt, S.M.; Segovia, J. Nucleon elastic form factors at accessible large spacelike momenta. Phys. Rev. D 2020, 102, 014043. [Google Scholar] [CrossRef]
  366. Xing, H.Y.; Xu, Z.N.; Cui, Z.F.; Roberts, C.D.; Xu, C. Heavy + heavy and heavy + light pseudoscalar to vector semileptonic transitions. Eur. Phys. J. C 2022, 82, 889. [Google Scholar] [CrossRef]
  367. Isgur, N.; Wise, M.B. Weak Transition Form Factors between Heavy Mesons. Phys. Lett. B 1990, 237, 527–530. [Google Scholar] [CrossRef]
  368. Amhis, Y.; Banerjee, S.; Ben-Haim, E.; Bernlochner, F.U.; Bona, M.; Bozek, A.; Bozzi, C.; Brodzicka, J.; Chrzaszcz, M.; Dingfelder, J.; et al. Averages of b-hadron, c-hadron, and τ-lepton properties as of 2018. Eur. Phys. J. C 2021, 81, 226. [Google Scholar]
  369. Glattauer, R. Measurement of the decay BDℓν in fully reconstructed events and determination of the Cabibbo–Kobayashi–Maskawa matrix element |Vcb|. Phys. Rev. D 2016, 93, 032006. [Google Scholar] [CrossRef]
  370. Ablikim, M.; Achasov, M.N.; Ahmed, S.; Albrecht, M.; Alekseev, M.; Amoroso, A.; An, F.F.; An, Q.; Bai, Y.; Bakina, O.; et al. First Measurement of the Form Factors in D s + K0e+νe and D s + K*0e+νe Decays. Phys. Rev. Lett. 2019, 122, 061801. [Google Scholar] [CrossRef] [Green Version]
  371. LHCb Collaboration; Alfonso Albero, A.; Badalov, A.; Camboni, A.; Coquereau, S.; Fernandez, G.; Garrido Beltrán, L. Measurement of the ratio of branching fractions B( B c + J/ψτ+ντ)/B( B c + J/ψμ+νμ). Phys. Rev. Lett. 2018, 120, 121801. [Google Scholar]
  372. Tran, C.T.; Ivanov, M.A.; Körner, J.G.; Santorelli, P. Implications of new physics in the decays Bc→(J/ψ,ηc)τν. Phys. Rev. D 2018, 97, 054014. [Google Scholar] [CrossRef] [Green Version]
  373. Issadykov, A.; Ivanov, M.A. The decays BcJ/ψ+ ¯ ν and BcJ/ψ+π(K) in covariant confined quark model. Phys. Lett. B 2018, 783, 178–182. [Google Scholar] [CrossRef]
  374. Wang, W.; Zhu, R. Model independent investigation of the RJ/ψ,ηc and ratios of decay widths of semileptonic Bc decays into a P-wave charmonium. Int. J. Mod. Phys. A 2019, 34, 1950195. [Google Scholar] [CrossRef]
  375. Leljak, D.; Melic, B.; Patra, M. On lepton flavour universality in semileptonic Bcηc,J/ψ decays. J. High Energy Phys. 2019, 5, 094. [Google Scholar] [CrossRef] [Green Version]
  376. Hu, X.Q.; Jin, S.P.; Xiao, Z.J. Semileptonic decays Bc→(ηc,J/ψ)l ν ¯ l in the “PQCD + Lattice” approach. Chin. Phys. C 2020, 44, 023104. [Google Scholar] [CrossRef]
  377. Zhou, T.; Wang, T.; Jiang, Y.; Tan, X.Z.; Li, G.; Wang, G.L. Relativistic calculations of R(D(*)), R( D s (*) ), R(ηc) and R(J/ψ). Int. J. Mod. Phys. A 2020, 35, 2050076. [Google Scholar] [CrossRef]
  378. Harrison, J.; Davies, C.T.H.; Lytle, A. R(J/ψ) and Bc-→J/ψℓ ν ¯ Lepton Flavor Universality Violating Observables from Lattice QCD. Phys. Rev. Lett. 2020, 125, 222003. [Google Scholar] [CrossRef] [PubMed]
  379. Harrison, J.; Davies, C.T.H.; Lytle, A. BcJ/ψ form factors for the full q2 range from lattice QCD. Phys. Rev. D 2020, 102, 094518. [Google Scholar] [CrossRef]
  380. Widhalm, L.; Abe, K.; Abe, K.; Adachi, I.; Aihara, H.; Arinstein, K.; Asano, Y.; Aushev, T.; Bakich, A.M.; Balagura, V.; et al. Measurement of D0πlν(Klν) Form Factors and Absolute Branching Fractions. Phys. Rev. Lett. 2006, 97, 061804. [Google Scholar] [CrossRef] [PubMed]
  381. Besson, D.; Pedlar, T.K.; Xavier, J.; Cronin-Hennessy, D.; Gao, K.Y.; Hietala, J.; Kubota, Y.; Klein, T.; Poling, R.; Scott, A.W.; et al. Improved measurements of D meson semileptonic decays to π and K mesons. Phys. Rev. D 2009, 80, 032005. [Google Scholar] [CrossRef] [Green Version]
  382. Ablikim, M.; Achasov, M.N.; Ahmed, S.; Albrecht, M.; Alekseev, M.; Amoroso, A.; An, F.F.; An, Q.; Bai, J.Z.; Bai, Y.; et al. Study of Dynamics of D0Ke+νe and D0πe+νe Decays. Phys. Rev. D 2015, 92, 072012. [Google Scholar] [CrossRef] [Green Version]
  383. Ablikim, M.; Achasov, M.N.; Ahmed, S.; Ai, X.C.; Albayrak, O.; Albrecht, M.; Ambrose, D.J.; Amoroso, A.; An, F.F.; An, Q.; et al. Analysis of D+ K ¯ 0e+νe and D+π0e+νe semileptonic decays. Phys. Rev. D 2017, 96, 012002. [Google Scholar] [CrossRef] [Green Version]
  384. Lubicz, V.; Riggio, L.; Salerno, G.; Simula, S.; Tarantino, C. Scalar and vector form factors of Dπ(K)ℓν decays with Nf=2+1+1 twisted fermions. Phys. Rev. D 2017, 96, 054514, Erratum: Phys. Rev. D 2019, 99, 099902. [Google Scholar] [CrossRef]
  385. Del Amo Sanchez, P.; Lees, J.P.; Poireau, V.; Prencipe, E.; Tisserand, V.; Garra Tico, J.; Grauges, E.; Martinelli, M.; Palano, A.; Pappagallo, M.; et al. Study of Bπℓν and Bρℓν Decays and Determination of |Vub|. Phys. Rev. D 2011, 83, 032007. [Google Scholar] [CrossRef] [Green Version]
  386. Ha, H.; Won, E.; Adachi, I.; Aihara, H.; Aziz, T.; Bakich, A.M.; Balagura, V.; Barberio, E.; Bay, A.; Belous, K.; et al. Measurement of the decay B0π+ν and determination of |Vub|. Phys. Rev. D 2011, 83, 071101. [Google Scholar] [CrossRef] [Green Version]
  387. Lees, J.P.; Poireau, V.; Tisserand, V.; Garra Tico, J.; Grauges, E.; Palano, A.; Eigen, G.; Stugu, B.; Brown, D.N.; Kerth, L.T.; et al. Branching fraction and form-factor shape measurements of exclusive charmless semileptonic B decays, and determination of |Vub|. Phys. Rev. D 2012, 86, 092004. [Google Scholar] [CrossRef] [Green Version]
  388. Sibidanov, A.; Varvell, K.E.; Adachi, I.; Aihara, H.; Asner, D.M.; Aulchenko, V.; Aushev, T.; Bakich, A.M.; Bala, A.; Bozek, A.; et al. Study of Exclusive BXuℓν Decays and Extraction of ‖Vub‖ using Full Reconstruction Tagging at the Belle Experiment. Phys. Rev. D 2013, 88, 032005. [Google Scholar] [CrossRef] [Green Version]
  389. Bouchard, C.M.; Lepage, G.P.; Monahan, C.; Na, H.; Shigemitsu, J. BsKℓν form factors from lattice QCD. Phys. Rev. D 2014, 90, 054506. [Google Scholar] [CrossRef] [Green Version]
  390. Flynn, J.M.; Izubuchi, T.; Kawanai, T.; Lehner, C.; Soni, A.; Van de Water, R.S.; Witzel, O. Bπℓν and BsKℓν form factors and |Vub| from 2+1-flavor lattice QCD with domain wall light quarks and relativistic heavy quarks. Phys. Rev. D 2015, 91, 074510. [Google Scholar] [CrossRef]
  391. Aaij, R.; Abellán Beteta, C.; Ackernley, T.; Adeva, B.; Adinolfi, M.; Afsharnia, H.; Aidala, C.A.; Aiola, S.; Ajaltouni, Z.; Akar, S.; et al. First observation of the decay B s 0 Kμ+νμ and Measurement of |Vub|/|Vcb|. Phys. Rev. Lett. 2021, 126, 081804. [Google Scholar] [CrossRef]
  392. Bazavov, A.; Bernard, C.; DeTar, C.; Du, D.; El-Khadra, A.X.; Freeland, E.D.; Gámiz, E.; Gelzer, Z.; Gottlieb, S.; Heller, U.M.; et al. BsKℓν decay from lattice QCD. Phys. Rev. D 2019, 100, 034501. [Google Scholar] [CrossRef] [Green Version]
  393. Wu, Y.L.; Zhong, M.; Zuo, Y.B. B(s),D(s)π,K,η,ρ,K*,ω,ϕ Transition Form Factors and Decay Rates with Extraction of the CKM parameters |Vub|, |Vcs|, |Vcd|. Int. J. Mod. Phys. A 2006, 21, 6125–6172. [Google Scholar] [CrossRef] [Green Version]
  394. Faustov, R.N.; Galkin, V.O. Charmless weak Bs decays in the relativistic quark model. Phys. Rev. D 2013, 87, 094028. [Google Scholar] [CrossRef] [Green Version]
  395. Xiao, Z.J.; Fan, Y.Y.; Wang, W.F.; Cheng, S. The semileptonic decays of B/Bs meson in the perturbative QCD approach: A short review. Chin. Sci. Bull. 2014, 59, 3787–3800. [Google Scholar] [CrossRef] [Green Version]
  396. Gonzàlez-Solís, S.; Masjuan, P.; Rojas, C. Padé approximants to Bπℓν and BsKℓν and determination of |Vub|. arXiv 2021, arXiv:2110.06153. [Google Scholar] [CrossRef]
  397. Melikhov, D. Dispersion approach to quark binding effects in weak decays of heavy mesons. Eur. Phys. J. Direct C 2002, 4, 1–154. [Google Scholar] [CrossRef] [Green Version]
  398. Faessler, A.; Gutsche, T.; Ivanov, M.A.; Korner, J.; Lyubovitskij, V.E. The Exclusive rare decays B→ K(K*) ¯ and Bc→ D(D*) ¯ in a relativistic quark model. Eur. Phys. J. Direct C 2002, 4, 18. [Google Scholar]
  399. Ebert, D.; Faustov, R.N.; Galkin, V.O. Weak decays of the Bc meson to Bs and B mesons in the relativistic quark model. Eur. Phys. J. C 2003, 32, 29–43. [Google Scholar] [CrossRef] [Green Version]
  400. Ball, P.; Zwicky, R. New results on Bπ,K,η decay form factors from light-cone sum rules. Phys. Rev. D 2005, 71, 014015. [Google Scholar] [CrossRef] [Green Version]
  401. Khodjamirian, A.; Mannel, T.; Offen, N. Form-factors from light-cone sum rules with B-meson distribution amplitudes. Phys. Rev. D 2007, 75, 054013. [Google Scholar] [CrossRef]
  402. Lü, C.D.; Wang, W.; Wei, Z.T. Heavy-to-light form factors on the light cone. Phys. Rev. D 2007, 76, 014013. [Google Scholar] [CrossRef] [Green Version]
  403. Ivanov, M.A.; Körner, J.G.; Kovalenko, S.G.; Roberts, C.D. B- to light-meson transition form factors. Phys. Rev. D 2007, 76, 034018. [Google Scholar] [CrossRef] [Green Version]
  404. El-Bennich, B.; Ivanov, M.A.; Roberts, C.D. Strong D* -> D+pi and B* -> B+pi couplings. Phys. Rev. C 2011, 83, 025205. [Google Scholar] [CrossRef] [Green Version]
  405. Ji, C.R.; Maris, P. Kℓ3 transition form-factors. Phys. Rev. 2001, D64, 014032. [Google Scholar]
  406. Chen, C.; Chang, L.; Roberts, C.D.; Wan, S.L.; Schmidt, S.M.; Wilson, D.J. Features and flaws of a contact interaction treatment of the kaon. Phys. Rev. C 2013, 87, 045207. [Google Scholar] [CrossRef] [Green Version]
  407. Holt, R.J.; Roberts, C.D. Distribution Functions of the Nucleon and Pion in the Valence Region. Rev. Mod. Phys. 2010, 82, 2991–3044. [Google Scholar] [CrossRef] [Green Version]
  408. Taylor, R.E. Deep inelastic scattering: The Early years. Rev. Mod. Phys. 1991, 63, 573–595. [Google Scholar] [CrossRef] [Green Version]
  409. Kendall, H.W. Deep inelastic scattering: Experiments on the proton and the observation of scaling. Rev. Mod. Phys. 1991, 63, 597–614. [Google Scholar] [CrossRef]
  410. Friedman, J.I. Deep inelastic scattering: Comparisons with the quark model. Rev. Mod. Phys. 1991, 63, 615–629. [Google Scholar] [CrossRef] [Green Version]
  411. Friedman, J.I.; Kendall, H.W.; Taylor, R.E. Deep inelastic scattering: Acknowledgements. Rev. Mod. Phys. 1991, 63, 629. [Google Scholar] [CrossRef]
  412. Brodsky, S.J.; Burkardt, M.; Schmidt, I. Perturbative QCD constraints on the shape of polarized quark and gluon distributions. Nucl. Phys. B 1995, 441, 197–214. [Google Scholar] [CrossRef]
  413. Yuan, F. Generalized parton distributions at x→1. Phys. Rev. D 2004, 69, 051501. [Google Scholar] [CrossRef] [Green Version]
  414. Dokshitzer, Y.L. Calculation of the Structure Functions for Deep Inelastic Scattering and e+ e- Annihilation by Perturbation Theory in Quantum Chromodynamics. Sov. Phys. JETP 1977, 46, 641–653. (In Russian) [Google Scholar]
  415. Lipatov, L.N. The parton model and perturbation theory. Sov. J. Nucl. Phys. 1975, 20, 94–102. [Google Scholar]
  416. Altarelli, G.; Parisi, G. Asymptotic Freedom in Parton Language. Nucl. Phys. B 1977, 126, 298–318. [Google Scholar] [CrossRef]
  417. Aicher, M.; Schäfer, A.; Vogelsang, W. Soft-Gluon Resummation and the Valence Parton Distribution Function of the Pion. Phys. Rev. Lett. 2010, 105, 252003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  418. Ball, R.D.; Nocera, E.R.; Rojo, J. The asymptotic behaviour of parton distributions at small and large x. Eur. Phys. J. C 2016, 76, 383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  419. Hou, T.-J.; Gao, J.; Hobbs, T.J.; Xie, K.; Dulat, S.; Guzzi, M.; Huston, J.; Nadolsky, P.; Pumplin, J.; Schmidt, C.; et al. New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC. Phys. Rev. D 2021, 103, 014013. [Google Scholar] [CrossRef]
  420. Bailey, S.; Cridge, T.; Harland-Lang, L.A.; Martin, A.D.; Thorne, R.S. Parton distributions from LHC, HERA, Tevatron and fixed target data: MSHT20 PDFs. Eur. Phys. J. C 2021, 81, 341. [Google Scholar] [CrossRef]
  421. Novikov, I.; Abdolmaleki, H.; Britzger, D.; Cooper-Sarkar, A.; Giuli, F.; Glazov, A.; Kusina, A.; Luszczak, A.; Olness, F.; Starovoitov, P.; et al. Parton Distribution Functions of the Charged Pion Within The xFitter Framework. Phys. Rev. D 2020, 102, 014040. [Google Scholar] [CrossRef]
  422. Barry, P.C.; Ji, C.R.; Sato, N.; Melnitchouk, W. Global QCD Analysis of Pion Parton Distributions with Threshold Resummation. Phys. Rev. Lett. 2021, 127, 232001. [Google Scholar] [CrossRef]
  423. Chang, L.; Cloet, I.C.; Cobos-Martinez, J.J.; Roberts, C.D.; Schmidt, S.M.; Tandy, P.C. Imaging dynamical chiral symmetry breaking: Pion wave function on the light front. Phys. Rev. Lett. 2013, 110, 132001. [Google Scholar] [CrossRef] [Green Version]
  424. Lee, T.D.; Yang, C.N. Charge Conjugation, a New Quantum Number G, and Selection Rules Concerning a Nucleon Anti-nucleon System. Nuovo Cim. 1956, 10, 749–753. [Google Scholar] [CrossRef]
  425. Gao, F.; Chang, L.; Liu, Y.X.; Roberts, C.D.; Schmidt, S.M. Parton distribution amplitudes of light vector mesons. Phys. Rev. D 2014, 90, 014011. [Google Scholar] [CrossRef] [Green Version]
  426. Lu, Y.; Binosi, D.; Ding, M.; Roberts, C.D.; Xing, H.Y.; Xu, C. Distribution amplitudes of light diquarks. Eur. Phys. J A (Lett) 2021, 57, 115. [Google Scholar] [CrossRef]
  427. Alexandrou, C.; Constantinou, M.; Drach, V.; Hatziyiannakou, K.; Jansen, K.; Kallidonis, C.; Koutsou, G.; Leontiou, T.; Vaquero, A. Nucleon Structure using lattice QCD. Nuovo Cim. C 2013, 036, 111–120. [Google Scholar]
  428. Lin, H.W.; Chen, J.W.; Zhang, R. Lattice Nucleon Isovector Unpolarized Parton Distribution in the Physical-Continuum Limit. arXiv 2020, arXiv:2011.14971. [Google Scholar]
  429. Alexandrou, C.; Constantinou, M.; Hadjiyiannakou, K.; Jansen, K.; Manigrasso, F. Flavor decomposition of the nucleon unpolarized, helicity, and transversity parton distribution functions from lattice QCD simulations. Phys. Rev. D 2021, 104, 054503. [Google Scholar] [CrossRef]
  430. Ji, X.; Liu, Y.S.; Liu, Y.; Zhang, J.H.; Zhao, Y. Large-momentum effective theory. Rev. Mod. Phys. 2021, 93, 035005. [Google Scholar] [CrossRef]
  431. Xu, S.S.; Chang, L.; Roberts, C.D.; Zong, H.S. Pion and kaon valence quark parton quasidistributions. Phys. Rev. D 2018, 97, 094014. [Google Scholar] [CrossRef] [Green Version]
  432. Field, R.D.; Feynman, R.P. Quark Elastic Scattering as a Source of High Transverse Momentum Mesons. Phys. Rev. D 1977, 15, 2590–2616. [Google Scholar] [CrossRef] [Green Version]
  433. Fan, Z.; Lin, H.W. Gluon parton distribution of the pion from lattice QCD. Phys. Lett. B 2021, 823, 136778. [Google Scholar] [CrossRef]
  434. Accardi, A.; Brady, L.T.; Melnitchouk, W.; Owens, J.F.; Sato, N. Constraints on large-x parton distributions from new weak boson production and deep-inelastic scattering data. Phys. Rev. D 2016, 93, 114017. [Google Scholar] [CrossRef] [Green Version]
  435. Ball, R.D.; Bertone, V.; Carrazza, S.; Del Debbio, L.; Forte, S.; Groth-Merrild, P.; Guffanti, A.; Hartland, N.P.; Kassabov, Z.; Latorre, J.I.; et al. Parton distributions from high-precision collider data. Eur. Phys. J. C 2017, 77, 663. [Google Scholar] [CrossRef] [PubMed]
  436. Courtoy, A.; Nadolsky, P.M. Testing momentum dependence of the nonperturbative hadron structure in a global QCD analysis. Phys. Rev. D 2021, 103, 054029. [Google Scholar] [CrossRef]
  437. Gottfried, K. Sum rule for high-energy electron - proton scattering. Phys. Rev. Lett. 1967, 18, 1174. [Google Scholar] [CrossRef]
  438. Brock, R.; Collins, J.C.; Huston, J.; Kuhlmann, S.; Mishra, S.; Morfin, J.G.; Olness, F.I.; Owens, J.F.; Pumplin, J.; Qiu, J.-W.; et al. Handbook of perturbative QCD: Version 1.0. Rev. Mod. Phys. 1995, 67, 157–248. [Google Scholar]
  439. Amaudruz, P.; Arneodo, M.; Arvidson, A.; Badelek, B.; Baum, G.; Beaufays, J.; Bird, I.G.; Botje, M.; Broggini, C.; Bruckner, W.; et al. The Gottfried sum from the ratio F 2 n / F 2 p . Phys. Rev. Lett. 1991, 66, 2712–2715. [Google Scholar] [CrossRef] [Green Version]
  440. Arneodo, M.; Arvidson, A.; Badelek, B.; Ballintijn@f, M.; Baum@f, G.; Beaufays, J.; Bird, I.G.; Björkholm, P.; Botje, M.; Broggini, C.; et al. A Reevaluation of the Gottfried sum. Phys. Rev. D 1994, 50, R1–R3. [Google Scholar] [CrossRef] [Green Version]
  441. Baldit, A.; Barrière, C.; Castor, J.; Chambon, T.; Devaux, A.; Espagnon, B.; Fargeix, J.; Force, P.; Landaud, G.; Saturnini, P.; et al. Study of the isospin symmetry breaking in the light quark sea of the nucleon from the Drell-Yan process. Phys. Lett. B 1994, 332, 244–250. [Google Scholar] [CrossRef] [Green Version]
  442. Towell, R.S.; McGaughey, P.L.; Awes, T.C.; Beddo, M.E.; Brooks, M.L.; Geesaman, D.F.; Kaufman, S.B.; Makins, N.; Mueller, B.A.; Reimer, P.E.; et al. Improved measurement of the d ¯ / u ¯ asymmetry in the nucleon sea. Phys. Rev. D 2001, 64, 052002. [Google Scholar] [CrossRef] [Green Version]
  443. Dove, J.; Kerns, B.; McClellan, R.E.; Miyasaka, S.; Morton, D.H.; Nagai, K.; Prasad, S.; Sanftl, F.; Scott, M.B.C.; Tadepalli, A.S.; et al. The asymmetry of antimatter in the proton. Nature 2021, 590, 561–565. [Google Scholar] [CrossRef]
  444. Tkachenko, S.; Baillie, N.; Kuhn, S.E.; Zhang, J.; Arrington, J.; Bosted, P.; Bültmann, S.; Christy, M.E.; Dutta, D.; Ent, R.; et al. Measurement of the structure function of the nearly free neutron using spectator tagging in inelastic 2H(e, e’p)X scattering with CLAS. Phys. Rev. C 2014, 89, 045206, Addendum: Phys. Rev. C 2014, 90, 059901. [Google Scholar] [CrossRef] [Green Version]
  445. Abrams, D.; Albataineh, H.; Aljawrneh, B.S.; Alsalmi, S.; Aniol, K.; Armstrong, W.; Arrington, J.; Atac, H.; Averett, T.; Ayerbe Gayoso, C.; et al. Measurement of the Nucleon F 2 n / F 2 p Structure Function Ratio by the Jefferson Lab MARATHON Tritium/Helium-3 Deep Inelastic Scattering Experiment. Phys. Rev. Lett. 2022, 128, 132003. [Google Scholar] [CrossRef]
  446. Farrar, G.R.; Jackson, D.R. Pion and Nucleon Structure Functions Near x=1. Phys. Rev. Lett. 1975, 35, 1416. [Google Scholar] [CrossRef] [Green Version]
  447. Brodsky, S.J.; Lepage, G.P. Perturbative Quantum Chromodynamics. Prog. Math. Phys. 1979, 4, 255–422. [Google Scholar]
  448. Roberts, C.D.; Holt, R.J.; Schmidt, S.M. Nucleon spin structure at very high x. Phys. Lett. B 2013, 727, 249–254. [Google Scholar] [CrossRef]
  449. Close, F.E.; Thomas, A.W. The Spin and Flavor Dependence of Parton Distribution Functions. Phys. Lett. B 1988, 212, 227. [Google Scholar] [CrossRef]
  450. Xu, S.S.; Chen, C.; Cloet, I.C.; Roberts, C.D.; Segovia, J.; Zong, H.S. Contact-interaction Faddeev equation and, inter alia, proton tensor charges. Phys. Rev. D 2015, 92, 114034. [Google Scholar] [CrossRef] [Green Version]
  451. Bodek, A.; Breidenbach, M.; Dubin, D.L.; Elias, J.E.; Friedman, J.I.; Kendall, H.W.; Poucher, J.S.; Riordan, E.M.; Sogard, M.R.; Coward, D.H. Comparisons of Deep Inelastic ep and en Cross-Sections. Phys. Rev. Lett. 1973, 30, 1087. [Google Scholar] [CrossRef] [Green Version]
  452. Poucher, J.S.; Breidenbach, M.; Ditzler, R.; Friedman, J.I.; Kendall, H.W.; Bloom, E.D.; Taylor, R.E. High-Energy Single-Arm Inelastic ep and ed Scattering at 6-Degrees and 10-Degrees. Phys. Rev. Lett. 1974, 32, 118. [Google Scholar] [CrossRef]
  453. Whitlow, L.W.; Riordan, E.M.; Dasu, S.; Rock, S.; Bodek, A. Precise measurements of the proton and deuteron structure functions from a global analysis of the SLAC deep inelastic electron scattering cross-sections. Phys. Lett. B 1992, 282, 475–482. [Google Scholar] [CrossRef]
  454. Afnan, I.R.; Bissey, F.R.P.; Gomez, J.; Katramatou, A.T.; Melnitchouk, W.; Petratos, G.G.; Thomas, A.W. Neutron structure function and A = 3 mirror nuclei. Phys. Lett. B 2000, 493, 36–42. [Google Scholar] [CrossRef] [Green Version]
  455. Pace, E.; Salme, G.; Scopetta, S.; Kievsky, A. Neutron structure function F 2 n (x) from deep inelastic electron scattering off few nucleon systems. Phys. Rev. C 2001, 64, 055203. [Google Scholar] [CrossRef] [Green Version]
  456. Segarra, E.; Schmidt, A.; Kutz, T.; Higinbotham, D.; Piasetzky, E.; Strikman, M.; Weinstein, L.; Hen, O. Neutron Valence Structure from Nuclear Deep Inelastic Scattering. Phys. Rev. Lett. 2020, 124, 092002. [Google Scholar] [CrossRef] [Green Version]
  457. Nachtmann, O. Positivity constraints for anomalous dimensions. Nucl. Phys. B 1973, 63, 237–247. [Google Scholar] [CrossRef]
  458. Cui, Z.F.; Gao, F.; Binosi, D.; Chang, L.; Roberts, C.D.; Schmidt, S.M. Valence quark ratio in the proton. Chin. Phys. Lett. Express 2022, 39, 041401. [Google Scholar] [CrossRef]
  459. Anselmino, M.; Predazzi, E.; Ekelin, S.; Fredriksson, S.; Lichtenberg, D.B. Diquarks. Rev. Mod. Phys. 1993, 65, 1199–1234. [Google Scholar] [CrossRef]
  460. Pace, E.; Rinaldi, M.; Salmè, G.; Scopetta, S. The European Muon Collaboration effect in Light-Front Hamiltonian Dynamics. arXiv 2022, arXiv:2206.05485. [Google Scholar]
  461. Joó, B.; Karpie, J.; Orginos, K.; Radyushkin, A.V.; Richards, D.G.; Sufian, R.S.; Zafeiropoulos, S. Pion valence structure from Ioffe-time parton pseudodistribution functions. Phys. Rev. D 2019, 100, 114512. [Google Scholar] [CrossRef] [Green Version]
  462. Alexandrou, C.; Bacchio, S.; Cloet, I.; Constantinou, M.; Hadjiyiannakou, K.; Koutsou, G.; Lauer, C. Pion and kaon 〈x3〉 from lattice QCD and PDF reconstruction from Mellin moments. Phys. Rev. D 2021, 104, 054504. [Google Scholar] [CrossRef]
  463. Brodsky, S.J.; Hoyer, P.; Peterson, C.; Sakai, N. The Intrinsic Charm of the Proton. Phys. Lett. B 1980, 93, 451–455. [Google Scholar] [CrossRef]
  464. Ball, R.D.; Candido, A.; Cruz-Martinez, J.; Forte, S.; Giani, T.; Hekhorn, F.; Kudashkin, K.; Magni, G.; Rojo, J. Evidence for intrinsic charm quarks in the proton. Nature 2022, 608, 483–487. [Google Scholar]
Figure 1. Mass budgets: (A) proton; (B) ρ -meson; (C) pion; (D) kaon. Each is drawn using a Poincaré-invariant decomposition and the numerical values listed in Table 1 (separation at ζ = 2 G e V , calculated using information from References [40,41,42,43]).
Figure 1. Mass budgets: (A) proton; (B) ρ -meson; (C) pion; (D) kaon. Each is drawn using a Poincaré-invariant decomposition and the numerical values listed in Table 1 (separation at ζ = 2 G e V , calculated using information from References [40,41,42,43]).
Particles 06 00004 g001
Figure 2. Renormalisation-group-invariant dressed gluon mass function (solid blue curve) calculated, following the method in Reference [72], from a gluon two-point function obtained using the lQCD configurations in References [73,74,75]. The mass-squared curve is plotted, normalised by its k = 0 value, and compared with the kindred chiral-limit dressed quark mass function drawn from Reference [76] (dotted-dashed green curve). It is this pair of curves that is 1 / k 2 -suppressed in the ultraviolet, each with additional logarithmic corrections.
Figure 2. Renormalisation-group-invariant dressed gluon mass function (solid blue curve) calculated, following the method in Reference [72], from a gluon two-point function obtained using the lQCD configurations in References [73,74,75]. The mass-squared curve is plotted, normalised by its k = 0 value, and compared with the kindred chiral-limit dressed quark mass function drawn from Reference [76] (dotted-dashed green curve). It is this pair of curves that is 1 / k 2 -suppressed in the ultraviolet, each with additional logarithmic corrections.
Particles 06 00004 g002
Figure 3. Process-independent effective charge, α ^ ( k ) / π , obtained by combining modern results from continuum and lattice analyses of QCD’s gauge sector [71]. Existing data on the process-dependent charge α g 1 [94,95], defined via the Bjorken sum rule, is shown for comparison—see References [95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121]. (Image courtesy of D. Binosi).
Figure 3. Process-independent effective charge, α ^ ( k ) / π , obtained by combining modern results from continuum and lattice analyses of QCD’s gauge sector [71]. Existing data on the process-dependent charge α g 1 [94,95], defined via the Bjorken sum rule, is shown for comparison—see References [95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121]. (Image courtesy of D. Binosi).
Particles 06 00004 g003
Figure 4. Panel (A): primitive gluon H-diagram. Panel (B): one-H-diagram contribution to quark+antiquark scattering. Legend: Dressed gluon two-point function—spring with open circle insertion; dressed quark two-point function—straight line with open circle insertion; blue-filled circle at 3-gluon junction—dressed 3-gluon vertex; red-filled circle at gluon–quark junction—dressed gluon–quark vertex. Repeated insertion of H-diagrams within H-diagrams, exploiting both gluon–quark and gluon self-couplings, leads to a plaquette-like area-filling structure, reminiscent of the planar summation of elementary squares in Reference [18].
Figure 4. Panel (A): primitive gluon H-diagram. Panel (B): one-H-diagram contribution to quark+antiquark scattering. Legend: Dressed gluon two-point function—spring with open circle insertion; dressed quark two-point function—straight line with open circle insertion; blue-filled circle at 3-gluon junction—dressed 3-gluon vertex; red-filled circle at gluon–quark junction—dressed gluon–quark vertex. Repeated insertion of H-diagrams within H-diagrams, exploiting both gluon–quark and gluon self-couplings, leads to a plaquette-like area-filling structure, reminiscent of the planar summation of elementary squares in Reference [18].
Particles 06 00004 g004
Figure 5. Δ ( τ ) μ 0 exp ( μ 0 τ 3 2 ) computed from the Schwinger function in Equation (18) using Equation (16).
Figure 5. Δ ( τ ) μ 0 exp ( μ 0 τ 3 2 ) computed from the Schwinger function in Equation (18) using Equation (16).
Particles 06 00004 g005
Figure 6. Empirical spectrum [40] (PDG summary tables): blue circles (bars)—u, d states; green diamonds (bar)—systems with s and/or s ¯ quarks. Little is known about K ( 1460 ) , which is therefore drawn as an open red diamond. Gold six-pointed stars—spectrum of low-lying u, d, and s mesons predicted by the EHM-improved Bethe–Salpeter kernel developed in Reference [222]; black five-pointed stars—same spectrum computed using RL truncation.
Figure 6. Empirical spectrum [40] (PDG summary tables): blue circles (bars)—u, d states; green diamonds (bar)—systems with s and/or s ¯ quarks. Little is known about K ( 1460 ) , which is therefore drawn as an open red diamond. Gold six-pointed stars—spectrum of low-lying u, d, and s mesons predicted by the EHM-improved Bethe–Salpeter kernel developed in Reference [222]; black five-pointed stars—same spectrum computed using RL truncation.
Particles 06 00004 g006
Figure 7. Leptonic decay constants for all states whose masses are reported in Figure 6: ground states, n = 0 ; lowest lying radial excitations, n = 1 . For the excited states, two extrapolation results are presented for each state, viz. one obtained with Padé approximants and the other employing the Schlessinger point method (SPM) [206,207,238,239]—the distinct approaches yield consistent outcomes. Results inferred from data are also plotted, where available [40] (PDG).
Figure 7. Leptonic decay constants for all states whose masses are reported in Figure 6: ground states, n = 0 ; lowest lying radial excitations, n = 1 . For the excited states, two extrapolation results are presented for each state, viz. one obtained with Padé approximants and the other employing the Schlessinger point method (SPM) [206,207,238,239]—the distinct approaches yield consistent outcomes. Results inferred from data are also plotted, where available [40] (PDG).
Particles 06 00004 g007
Figure 8. Figurative representation of the integral equation satisfied by the Poincaré-covariant matrix-valued function Ψ , viz. the Faddeev amplitude for a baryon with total momentum Q = q + d = k q + k d built from three valence quarks, two of which are always participants in a nonpointlike, interacting diquark correlation. Ψ describes the sharing of relative momentum between the dressed quarks and diquarks. Shaded rectangle—Faddeev kernel. Legend: single line—dressed quark propagator, S ( q ) ; Γ J P ( k ; K ) —diquark correlation amplitude; double line—diquark propagator, Δ J P ( K ) .
Figure 8. Figurative representation of the integral equation satisfied by the Poincaré-covariant matrix-valued function Ψ , viz. the Faddeev amplitude for a baryon with total momentum Q = q + d = k q + k d built from three valence quarks, two of which are always participants in a nonpointlike, interacting diquark correlation. Ψ describes the sharing of relative momentum between the dressed quarks and diquarks. Shaded rectangle—Faddeev kernel. Legend: single line—dressed quark propagator, S ( q ) ; Γ J P ( k ; K ) —diquark correlation amplitude; double line—diquark propagator, Δ J P ( K ) .
Particles 06 00004 g008
Figure 9. Upper panel (A): SCI-calculated masses of ground state flavour-SU ( 5 ) J P = 1 / 2 + baryons [232] compared with either experiment (first 15) [40] or lQCD (last 9) [262,263]. Lower panel (B): SCI masses of ground state flavour-SU ( 5 ) J P = 1 / 2 baryons in [232] compared with experiment [40] (green bars), lQCD [264] (gold triangles), or three-body Faddeev equation results [218] (orange asterisks). Analogous plots for J P = 3 / 2 ± baryons are presented elsewhere ([232] (Figures 4B and 5B)).
Figure 9. Upper panel (A): SCI-calculated masses of ground state flavour-SU ( 5 ) J P = 1 / 2 + baryons [232] compared with either experiment (first 15) [40] or lQCD (last 9) [262,263]. Lower panel (B): SCI masses of ground state flavour-SU ( 5 ) J P = 1 / 2 baryons in [232] compared with experiment [40] (green bars), lQCD [264] (gold triangles), or three-body Faddeev equation results [218] (orange asterisks). Analogous plots for J P = 3 / 2 ± baryons are presented elsewhere ([232] (Figures 4B and 5B)).
Particles 06 00004 g009
Figure 10. Real part of empirical pole position for each designated baryon [40] (red circle) compared with: calculated masses in Reference [268] (gold asterisks) after subtracting δ MB N 3 / 2 = 0.13 GeV from each; Reference [265] (teal diamonds) after subtracting δ MB N 1 / 2 = 0.30 GeV; and Reference [267] (green five-pointed stars) after subtracting δ MB Δ 3 / 2 = 0.17 GeV. All Faddeev equation predictions are drawn with an uncertainty that reflects a ± 5 % change in diquark masses.
Figure 10. Real part of empirical pole position for each designated baryon [40] (red circle) compared with: calculated masses in Reference [268] (gold asterisks) after subtracting δ MB N 3 / 2 = 0.13 GeV from each; Reference [265] (teal diamonds) after subtracting δ MB N 1 / 2 = 0.30 GeV; and Reference [267] (green five-pointed stars) after subtracting δ MB Δ 3 / 2 = 0.17 GeV. All Faddeev equation predictions are drawn with an uncertainty that reflects a ± 5 % change in diquark masses.
Particles 06 00004 g010
Figure 11. Rest-frame quark+axial vector–diquark orbital angular momentum content of ( 3 2 , 3 2 ± ) baryons considered in Reference [267], as measured by the contribution from the various components to the associated canonical normalisation constant: A Δ ( 1232 ) 3 2 + ; B Δ ( 1600 ) 3 2 + ; C Δ ( 1700 ) 3 2 ; D Δ ( 1940 ) 3 2 . The overall positive normalisations receive both positive (above plane) and negative (below plane) contributions. The central image provides the legend for the interpretation of the other panels, identifying interference between the various identified orbital angular momentum basis components in the baryon rest-frame.
Figure 11. Rest-frame quark+axial vector–diquark orbital angular momentum content of ( 3 2 , 3 2 ± ) baryons considered in Reference [267], as measured by the contribution from the various components to the associated canonical normalisation constant: A Δ ( 1232 ) 3 2 + ; B Δ ( 1600 ) 3 2 + ; C Δ ( 1700 ) 3 2 ; D Δ ( 1940 ) 3 2 . The overall positive normalisations receive both positive (above plane) and negative (below plane) contributions. The central image provides the legend for the interpretation of the other panels, identifying interference between the various identified orbital angular momentum basis components in the baryon rest-frame.
Particles 06 00004 g011
Figure 12. Electroproduction of a vector meson from the proton: e + p e + V + p , often interpreted as providing access to V + p V + p using a vector meson dominance model. The γ ( * ) ( Q ) V VMD transition is indicated by the crossed circle. It is usually assumed to occur with a momentum-independent strength γ γ V [320,321,322]. However, regarding heavy mesons, at least, the VMD hypothesis is unsound, as discussed below and in, e.g., References [204,323].
Figure 12. Electroproduction of a vector meson from the proton: e + p e + V + p , often interpreted as providing access to V + p V + p using a vector meson dominance model. The γ ( * ) ( Q ) V VMD transition is indicated by the crossed circle. It is usually assumed to occur with a momentum-independent strength γ γ V [320,321,322]. However, regarding heavy mesons, at least, the VMD hypothesis is unsound, as discussed below and in, e.g., References [204,323].
Particles 06 00004 g012
Figure 13. Left panel (A): First ratio in Equation (29) computed using matched solutions of the gap and Bethe–Salpeter equations for V = ρ , ϕ , J / ψ , Υ . Right panel (B): Second ratio in Equation (29), computed similarly. In cases where the VMD hypothesis were sound, all these curves would lie close to the thin horizontal line drawn at unity.
Figure 13. Left panel (A): First ratio in Equation (29) computed using matched solutions of the gap and Bethe–Salpeter equations for V = ρ , ϕ , J / ψ , Υ . Right panel (B): Second ratio in Equation (29), computed similarly. In cases where the VMD hypothesis were sound, all these curves would lie close to the thin horizontal line drawn at unity.
Particles 06 00004 g013
Figure 14. Axial/pseudoscalar current that ensures a symmetry-preserving interaction with an on-shell baryon described by a Faddeev amplitude obtained from the equation depicted in Figure 8: single line, dressed quark propagator; undulating line, axial/pseudoscalar current; Γ , diquark correlation amplitude; double line, diquark propagator; χ , seagull terms. A detailed legend is provided in Table 2.
Figure 14. Axial/pseudoscalar current that ensures a symmetry-preserving interaction with an on-shell baryon described by a Faddeev amplitude obtained from the equation depicted in Figure 8: single line, dressed quark propagator; undulating line, axial/pseudoscalar current; Γ , diquark correlation amplitude; double line, diquark propagator; χ , seagull terms. A detailed legend is provided in Table 2.
Particles 06 00004 g014
Figure 15. Predicted result for G A ( x ) in Reference [275] (Chen et al.) (blue curve within lighter blue uncertainty band), compared with lQCD results from Reference [347] (Jang et al.) (green) diamonds. With respect to the CSM central results, this comparison may be quantified by reporting the mean- χ 2 value, which is 0.27.
Figure 15. Predicted result for G A ( x ) in Reference [275] (Chen et al.) (blue curve within lighter blue uncertainty band), compared with lQCD results from Reference [347] (Jang et al.) (green) diamonds. With respect to the CSM central results, this comparison may be quantified by reporting the mean- χ 2 value, which is 0.27.
Particles 06 00004 g015
Figure 16. Comparison of the CSM prediction for g π N N / m N (blue asterisk) [275] with values extracted from pion–nucleon scattering [351] (Baru) (black star), the Granada 2013 n p and p p scattering database [352] (Navarro) (purple circle), and nucleon–nucleon scattering [353] (Reinert) (red triangle). The vertical grey band marks the estimated uncertainty in the CSM prediction. The error-weighted average of the depicted results, Equation (37), is drawn as the gold line within the like-coloured band.
Figure 16. Comparison of the CSM prediction for g π N N / m N (blue asterisk) [275] with values extracted from pion–nucleon scattering [351] (Baru) (black star), the Granada 2013 n p and p p scattering database [352] (Navarro) (purple circle), and nucleon–nucleon scattering [353] (Reinert) (red triangle). The vertical grey band marks the estimated uncertainty in the CSM prediction. The error-weighted average of the depicted results, Equation (37), is drawn as the gold line within the like-coloured band.
Particles 06 00004 g016
Figure 17. Transverse density profiles, Equation (48), calculated in Reference [276] from flavour-separated proton axial form factors. Top row: scalar-diquark-only proton. Panel A: two-dimensional plot of ρ ^ A d ( | b ^ | ) / g A d ; Panel B: similar plot of ρ ^ A u ( | b ^ | ) / g A u . Removing the 1 / g A f normalisation, the b = 0 values are ρ ^ A d ( 0 ) = 0.009 , ρ A u ( 0 ) = 0.097 . Bottom row: realistic proton Faddeev amplitude, including axial vector diquarks as predicted by the Faddeev equation in Figure 8. Panel C: ρ ^ A d ( | b ^ | ) / g A d ; Panel D: ρ ^ A u ( | b ^ | ) / g A u . Removing the 1 / g A f normalisation, the b = 0 values of these profiles are ρ ^ A d ( 0 ) = 0.038 , ρ ^ A u ( 0 ) = 0.12 . N.B. d 2 b ^ ρ ^ A f ( | b ^ | ) / g A f = 1 , f = u , d .
Figure 17. Transverse density profiles, Equation (48), calculated in Reference [276] from flavour-separated proton axial form factors. Top row: scalar-diquark-only proton. Panel A: two-dimensional plot of ρ ^ A d ( | b ^ | ) / g A d ; Panel B: similar plot of ρ ^ A u ( | b ^ | ) / g A u . Removing the 1 / g A f normalisation, the b = 0 values are ρ ^ A d ( 0 ) = 0.009 , ρ A u ( 0 ) = 0.097 . Bottom row: realistic proton Faddeev amplitude, including axial vector diquarks as predicted by the Faddeev equation in Figure 8. Panel C: ρ ^ A d ( | b ^ | ) / g A d ; Panel D: ρ ^ A u ( | b ^ | ) / g A u . Removing the 1 / g A f normalisation, the b = 0 values of these profiles are ρ ^ A d ( 0 ) = 0.038 , ρ ^ A u ( 0 ) = 0.12 . N.B. d 2 b ^ ρ ^ A f ( | b ^ | ) / g A f = 1 , f = u , d .
Particles 06 00004 g017
Figure 18. Mass budgets: (A): D * meson; (B): B * -meson; (C): D meson; (D): B meson. Each is drawn using a Poincaré-invariant decomposition and the numerical values listed in Table 1. (Separation at ζ = 2 G e V , calculated using information from References [35,40]).
Figure 18. Mass budgets: (A): D * meson; (B): B * -meson; (C): D meson; (D): B meson. Each is drawn using a Poincaré-invariant decomposition and the numerical values listed in Table 1. (Separation at ζ = 2 G e V , calculated using information from References [35,40]).
Particles 06 00004 g018
Figure 19. Ratio R J / ψ in Equation (51)—red circle, empirical result [371] (LHCb); blue asterisk—CSM prediction [206]; green star—SCI prediction [366]; grey circle—lQCD result [378,379]; gold band—unweighted mean of central values from several calculations [372,373,374,375,376,377].
Figure 19. Ratio R J / ψ in Equation (51)—red circle, empirical result [371] (LHCb); blue asterisk—CSM prediction [206]; green star—SCI prediction [366]; grey circle—lQCD result [378,379]; gold band—unweighted mean of central values from several calculations [372,373,374,375,376,377].
Particles 06 00004 g019
Figure 20. D K transition form factors. f + —solid blue curve; f 0 —long-dashed purple curve; like-coloured shaded bands—associated SPM uncertainty in each case. Data: cyan squares [380]; green diamonds [381]; dark-blue up-triangles [382]; indigo down-triangles [383].
Figure 20. D K transition form factors. f + —solid blue curve; f 0 —long-dashed purple curve; like-coloured shaded bands—associated SPM uncertainty in each case. Data: cyan squares [380]; green diamonds [381]; dark-blue up-triangles [382]; indigo down-triangles [383].
Particles 06 00004 g020
Figure 21. CSM predictions for B ( s ) semileptonic transition form factors: f + —solid blue curve; f 0 —long-dashed purple curve; SPM uncertainty in each—like-coloured shaded bands. Left panel (A): B ¯ 0 π + . Data (green stars): reconstructed from the average [368] (Table 81) of data reported in References [385,386,387,388]. Right panel (B): B ¯ s 0 K + . lQCD results: f + —indigo open up-triangles [389] and green open boxes [390]; f 0 —brown open down-triangles [389] and red open circles [390].
Figure 21. CSM predictions for B ( s ) semileptonic transition form factors: f + —solid blue curve; f 0 —long-dashed purple curve; SPM uncertainty in each—like-coloured shaded bands. Left panel (A): B ¯ 0 π + . Data (green stars): reconstructed from the average [368] (Table 81) of data reported in References [385,386,387,388]. Right panel (B): B ¯ s 0 K + . lQCD results: f + —indigo open up-triangles [389] and green open boxes [390]; f 0 —brown open down-triangles [389] and red open circles [390].
Particles 06 00004 g021
Figure 22. Left panel (A): Branching fraction B B ¯ s 0 K + μ ν ¯ μ computed in Reference [207], “Y”, compared with the value in Equation (55), “E”, viz. a measurement of B B s 0 K μ + ν μ [391], and some results obtained using other approaches: continuum I–III [393,394,395]; lattice IV–VI [389,390,392]. Right panel (B): Branching fraction ratio B B ¯ s 0 K + τ ν ¯ τ / B B ¯ s 0 K + μ ν ¯ μ computed in Reference [207] compared with some results obtained using other approaches. The legend matches Panel A, except “P(E)” is the result from [396], i.e., an estimate constrained by the datum in Reference [391]. Both panels: Grey line: unweighted mean of theory results. Pink dotted-dashed line: unweighted mean of theory results, omitting V–VI. Green dashed line: uncertainty-weighted average of theory results. Like-coloured bands mark associated uncertainties in each case.
Figure 22. Left panel (A): Branching fraction B B ¯ s 0 K + μ ν ¯ μ computed in Reference [207], “Y”, compared with the value in Equation (55), “E”, viz. a measurement of B B s 0 K μ + ν μ [391], and some results obtained using other approaches: continuum I–III [393,394,395]; lattice IV–VI [389,390,392]. Right panel (B): Branching fraction ratio B B ¯ s 0 K + τ ν ¯ τ / B B ¯ s 0 K + μ ν ¯ μ computed in Reference [207] compared with some results obtained using other approaches. The legend matches Panel A, except “P(E)” is the result from [396], i.e., an estimate constrained by the datum in Reference [391]. Both panels: Grey line: unweighted mean of theory results. Pink dotted-dashed line: unweighted mean of theory results, omitting V–VI. Green dashed line: uncertainty-weighted average of theory results. Like-coloured bands mark associated uncertainties in each case.
Particles 06 00004 g022
Figure 23. Left panel (A); Proton and pion hadron scale valence parton DFs: x u p ( x ; ζ H ) —solid red curve; x d p ( x ; ζ H ) —dotted-dashed blue curve; x u π ( x ; ζ H ) —dashed green curve. Right panel (B): Valence DFs in Panel A evolved to ζ 3 = m J / ψ = 3.097 GeV according to P1. The band surrounding each CSM curve expresses the response to a ± 5 % variation in ζ H .
Figure 23. Left panel (A); Proton and pion hadron scale valence parton DFs: x u p ( x ; ζ H ) —solid red curve; x d p ( x ; ζ H ) —dotted-dashed blue curve; x u π ( x ; ζ H ) —dashed green curve. Right panel (B): Valence DFs in Panel A evolved to ζ 3 = m J / ψ = 3.097 GeV according to P1. The band surrounding each CSM curve expresses the response to a ± 5 % variation in ζ H .
Particles 06 00004 g023
Figure 24. Isovector distribution [ u p ( x ; ζ ) d p ( x ; ζ ) ] . CSM prediction—solid purple curve, ζ = ζ 3 ; lQCD result from Reference [428]—dashed grey curve, ζ = ζ 3 ; lQCD result from Reference [429]—dotted-dashed brown curve, ζ = ζ 2 . Like-coloured band bracketing each curve indicates associated uncertainty.
Figure 24. Isovector distribution [ u p ( x ; ζ ) d p ( x ; ζ ) ] . CSM prediction—solid purple curve, ζ = ζ 3 ; lQCD result from Reference [428]—dashed grey curve, ζ = ζ 3 ; lQCD result from Reference [429]—dotted-dashed brown curve, ζ = ζ 2 . Like-coloured band bracketing each curve indicates associated uncertainty.
Particles 06 00004 g024
Figure 25. Left upper panel (A): Glue DFs— x 2 g , in the proton (solid purple curve) and pion (dashed green curve) at ζ = ζ 3 . Right upper panel (B): Comparison between continuum [127] (Cui 2020) and lattice [433] (Fan 2021) results for the glue-in-pion DF at ζ = 2 GeV. Left lower panel (C): Proton and pion light quark sea DFs: x 2 S u p ( x ; ζ 3 ) —solid red curve; x 2 S d p ( x ; ζ 3 ) —dashed blue curve; x 2 S u π ( x ; ζ 3 ) —dotted-dashed green curve. Right lower panel (D): Proton and pion c- and s quark sea DFs: x 2 S s p ( x ; ζ 3 ) —solid red curve; x 2 S s π ( x ; ζ 3 ) —dashed green curve; x 2 S c p ( x ; ζ 3 ) —dotted-dashed blue curve; x 2 S c π ( x ; ζ 3 ) —long-dashed orange curve. The band surrounding each CSM curve expresses the response to a ± 5 % variation in ζ H . The uncertainty in the lQCD result is similarly indicated in Panel B.
Figure 25. Left upper panel (A): Glue DFs— x 2 g , in the proton (solid purple curve) and pion (dashed green curve) at ζ = ζ 3 . Right upper panel (B): Comparison between continuum [127] (Cui 2020) and lattice [433] (Fan 2021) results for the glue-in-pion DF at ζ = 2 GeV. Left lower panel (C): Proton and pion light quark sea DFs: x 2 S u p ( x ; ζ 3 ) —solid red curve; x 2 S d p ( x ; ζ 3 ) —dashed blue curve; x 2 S u π ( x ; ζ 3 ) —dotted-dashed green curve. Right lower panel (D): Proton and pion c- and s quark sea DFs: x 2 S s p ( x ; ζ 3 ) —solid red curve; x 2 S s π ( x ; ζ 3 ) —dashed green curve; x 2 S c p ( x ; ζ 3 ) —dotted-dashed blue curve; x 2 S c π ( x ; ζ 3 ) —long-dashed orange curve. The band surrounding each CSM curve expresses the response to a ± 5 % variation in ζ H . The uncertainty in the lQCD result is similarly indicated in Panel B.
Particles 06 00004 g025
Figure 26. Left panel (A): Ratio of light antiquark DFs. Data: Reference [443] (E906). CSM result (solid purple curve) obtained from valence quark DFs in Figure 23 after evolution to ζ 2 = ζ SQ 2 = 30 GeV 2 [135]. The shaded band expresses the impact of a ± 25 % variation in the strength of Pauli blocking. Right panel (B): Neutron-to-proton structure function ratio. Data: [444] (BoNuS)—open grey squares; [445] (MARATHON)—gold asterisks. Contemporary CSM results (solid purple curve) obtained from valence quark DFs in Figure 23 after evolution to ζ = ζ M = 2.7 GeV [134,135]. Other predictions: green star—helicity conservation in the QCD parton model [412,446,447]; red diamond—large-x estimate based on Faddeev equation solutions [448]; retaining only scalar diquarks in the proton wave function, which produces a large-x value for this ratio that lies in the neighbourhood of the filled circle [449,450]. The band surrounding the CSM curve expresses the response to a ± 5 % variation in the size of axial vector diquark contributions to the proton charge. It is only noticeable on the valence quark domain.
Figure 26. Left panel (A): Ratio of light antiquark DFs. Data: Reference [443] (E906). CSM result (solid purple curve) obtained from valence quark DFs in Figure 23 after evolution to ζ 2 = ζ SQ 2 = 30 GeV 2 [135]. The shaded band expresses the impact of a ± 25 % variation in the strength of Pauli blocking. Right panel (B): Neutron-to-proton structure function ratio. Data: [444] (BoNuS)—open grey squares; [445] (MARATHON)—gold asterisks. Contemporary CSM results (solid purple curve) obtained from valence quark DFs in Figure 23 after evolution to ζ = ζ M = 2.7 GeV [134,135]. Other predictions: green star—helicity conservation in the QCD parton model [412,446,447]; red diamond—large-x estimate based on Faddeev equation solutions [448]; retaining only scalar diquarks in the proton wave function, which produces a large-x value for this ratio that lies in the neighbourhood of the filled circle [449,450]. The band surrounding the CSM curve expresses the response to a ± 5 % variation in the size of axial vector diquark contributions to the proton charge. It is only noticeable on the valence quark domain.
Particles 06 00004 g026
Figure 27. lim x 1 F 2 n ( x ) / F 2 p ( x ) . SPM prediction derived from Reference [445] (MARATHON) compared with results inferred from: nuclear DIS [456]; large-x estimate from CSM studies [271,450]; quark counting (helicity conservation) [446]; and a phenomenological fit (CJ15) [434]. The red vertical line indicates the Nachtmann lower limit, Equation (63), which is saturated if valence d quarks play no material role at x = 1 ; Row 3 is the average in Equation (64).
Figure 27. lim x 1 F 2 n ( x ) / F 2 p ( x ) . SPM prediction derived from Reference [445] (MARATHON) compared with results inferred from: nuclear DIS [456]; large-x estimate from CSM studies [271,450]; quark counting (helicity conservation) [446]; and a phenomenological fit (CJ15) [434]. The red vertical line indicates the Nachtmann lower limit, Equation (63), which is saturated if valence d quarks play no material role at x = 1 ; Row 3 is the average in Equation (64).
Particles 06 00004 g027
Table 1. Mass budgets of a collection of hadrons, with each panel ordered according to the contribution from Higgs boson couplings into QCD (HB) and including the component that is entirely unrelated to the Higgs (EHM) and that arising from constructive interference between these two mass sources (EHM+HB) (separation at ζ = 2 G e V , produced using information from References [35,40,41,42,43]).
Table 1. Mass budgets of a collection of hadrons, with each panel ordered according to the contribution from Higgs boson couplings into QCD (HB) and including the component that is entirely unrelated to the Higgs (EHM) and that arising from constructive interference between these two mass sources (EHM+HB) (separation at ζ = 2 G e V , produced using information from References [35,40,41,42,43]).
Mass Fraction (%)
Hadron (Mass/GeV)HBEHM+HBEHM
p ( 0.938 ) 1 6 93
ρ ( 0.775 ) 1 2 97
D * ( 2.010 ) 63 30 7
B * ( 5.325 ) 78 21 1
π ( 0.140 ) 5 95 0
K ( 0.494 ) 20 80 0
D ( 1.870 ) 68 32 0
B ( 5.279 ) 79 21 0
Table 2. Enumeration of terms in the baryon axial vector/pseudoscalar current, drawn in Figure 14.
Table 2. Enumeration of terms in the baryon axial vector/pseudoscalar current, drawn in Figure 14.
  • Diagram 1, two terms: J q S —probe strikes dressed quark with scalar diquark spectator; J q A —probe strikes dressed quark with axial vector diquark spectator.
  • Diagram 2: J qq A A —probe strikes axial vector diquark with dressed quark spectator.
  • Diagram 3: J qq { S A } —probe mediates transition between scalar and axial vector diquarks, with dressed quark spectator.
  • Diagram 4, three terms: J ex S S —probe strikes dressed quark “in-flight” between one scalar diquark correlation and another; J ex { S A } —dressed quark “in-flight” between a scalar diquark correlation and an axial vector correlation; J ex A A —“in-flight” between one axial vector correlation and another.
  • Diagrams 5 and 6, seagull diagrams describing the probe coupling into the diquark correlation amplitudes: J sg . There is one contribution from each diagram to match every term in Diagram (4).
Table 3. Flavour and diagram—Figure 14—separation of the proton axial charge: g A u = G A u ( 0 ) , g A d = G A d ( 0 ) ; g A u g A d = 1.25 ( 3 ) . The listed uncertainties express the effect of ± 5 % variations in the diquark masses, e.g., 0 . 88 6 0.88 0.06 .
Table 3. Flavour and diagram—Figure 14—separation of the proton axial charge: g A u = G A u ( 0 ) , g A d = G A d ( 0 ) ; g A u g A d = 1.25 ( 3 ) . The listed uncertainties express the effect of ± 5 % variations in the diquark masses, e.g., 0 . 88 6 0.88 0.06 .
J q S J q A J qq AA J qq { SA } J ex SS J ex { SA } J ex AA
g A u 0 . 88 6 0 . 08 0 ± 0 . 03 0 ± 0 . 08 0 0≈0 0 . 03 1 ±
g A d 0 0 . 16 0 ± 0 0 . 08 0 0 . 05 1 ± ≈0 0 . 01 0 ±
Table 4. CSM predictions for pseudoscalar meson semileptonic branching fractions [206,207]—each such fraction is to be multiplied by 10 3 . The column labelled “ratio” is the ratio of the preceding two entries in the row, so no factor of 10 3 is applied in this column. (A 1 σ SPM uncertainty is listed for the CSM predictions.) Reference [40] (PDG) lists the following values for the CKM matrix elements: | V u s | = 0.2245 ( 8 ) , | V c d | = 0.221 ( 4 ) , | V c s | = 0.987 ( 11 ) | V u b | = 0.00382 ( 24 ) ; and the following lifetimes (in seconds): τ K + = 1.2379 ( 21 ) × 10 8 , τ D 0 = 4.10 × 10 13 , τ D s ± = 5.04 × 10 13 , τ B ¯ 0 = 1.519 × 10 12 , τ B ¯ s 0 = 1.515 × 10 12 , τ B c ± = 0.51 × 10 12 .
Table 4. CSM predictions for pseudoscalar meson semileptonic branching fractions [206,207]—each such fraction is to be multiplied by 10 3 . The column labelled “ratio” is the ratio of the preceding two entries in the row, so no factor of 10 3 is applied in this column. (A 1 σ SPM uncertainty is listed for the CSM predictions.) Reference [40] (PDG) lists the following values for the CKM matrix elements: | V u s | = 0.2245 ( 8 ) , | V c d | = 0.221 ( 4 ) , | V c s | = 0.987 ( 11 ) | V u b | = 0.00382 ( 24 ) ; and the following lifetimes (in seconds): τ K + = 1.2379 ( 21 ) × 10 8 , τ D 0 = 4.10 × 10 13 , τ D s ± = 5.04 × 10 13 , τ B ¯ 0 = 1.519 × 10 12 , τ B ¯ s 0 = 1.515 × 10 12 , τ B c ± = 0.51 × 10 12 .
B I F ( ν ) References [206,207]PDG [40] or Other, If Indicated
e + ν e μ + ν μ ratio e + ν e μ + ν μ ratio
K + π 0 50.0 ( 9 ) 00 33.0 ( 6 ) 00 0.665 ( 70 ) 50.7 ( 6 ) 00 33.5 ( 3 ) 22 0.661 ( 07 )
D 0 π 1 2.70 ( 12 ) 2 2.66 ( 12 ) 0.987 ( 02 ) 0 2.91 ( 4 ) 1 0 2.67 ( 12 ) 0.918 ( 40 )
D s + K 0 1 2.73 ( 12 ) 2 2.68 ( 12 ) 0.982 ( 01 ) 00 3.25 ( 36 ) [370]
D 0 K 39.0 ( 1.7 ) 38.1 ( 1.7 ) 0.977 ( 01 ) 0 35.41 ( 34 ) 1 34.1 ( 4 ) 22 0.963 ( 10 )
μ ν ¯ μ τ ν ¯ τ ratio μ ν ¯ μ τ ν ¯ τ ratio
B ¯ 0 π + 11 0.162 ( 44 ) 11 0.120 ( 35 ) 0.733 ( 02 ) 01 0.150 ( 06 )
B ¯ s 0 K + 11 0.186 ( 53 ) 11 0.125 ( 37 ) 0.667 ( 09 )
B c η c 00 8.10 ( 45 ) 0 00 2.54 ( 10 ) 0 0.31 ( 2 ) 00
B c J / ψ 0 17.2 ( 1.9 ) 0 0 4.17 ( 66 ) 0.24 ( 5 ) 00
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ding, M.; Roberts, C.D.; Schmidt, S.M. Emergence of Hadron Mass and Structure. Particles 2023, 6, 57-120. https://doi.org/10.3390/particles6010004

AMA Style

Ding M, Roberts CD, Schmidt SM. Emergence of Hadron Mass and Structure. Particles. 2023; 6(1):57-120. https://doi.org/10.3390/particles6010004

Chicago/Turabian Style

Ding, Minghui, Craig D. Roberts, and Sebastian M. Schmidt. 2023. "Emergence of Hadron Mass and Structure" Particles 6, no. 1: 57-120. https://doi.org/10.3390/particles6010004

Article Metrics

Back to TopTop