Pseudoscalar Mesons and Emergent Mass

Despite its role in the continuing evolution of the Universe


Introduction
When studying the many facets of the Universe, one must be aware of the length/energy scales involved.For instance, the laws that seem to govern everyday phenomena might not be applicable at astronomical or subatomic scales.This apparent shift in perspective is probably only a reflection of our current limitations in comprehending the Universe.Despite the challenges, we have been able to separate the mass-energy budget of the universe into three sources: dark energy, dark matter and visible matter [1].While dark energy and dark matter contribute the largest portion (71% and 24%, respectively) [2], their effects on our daily life and accurate description are practically unimportant.Visible matter, on the other hand, constitutes only 5% of the whole; yet, it is the source of almost everything that is tangible.In terms of fundamental particles and their interactions, this small percentage is described by the so-called Standard Model (SM) of particle physics [3][4][5][6][7].The SM is remarkably successful in explaining observations; so much so that there are few empirical indications that it needs improvement -see, e.g., Refs.[8,9].Thus, before attempting to go beyond, it is worth resolving outstanding issues within the SM.
One such challenge is to explain the origin of visible mass.In this connection, the Higgs boson typically comes to mind [10,11].However, although it is vital to the character of our Universe, the Higgs only contributes about 1-2% of the visible mass.The remaining ≳ 98% owes to strong interactions in the SM, i.e., quantum chromodynamics (QCD) [12,13], through a dynamical source that is now often dubbed emergent hadron mass (EHM).These same strong interactions are also responsible for the formation of protons and neutrons (nucleons) and, from them, the atomic nuclei that constitute almost the entirety of visible material.
QCD is the Poincaré-invariant quantum non-Abelian gauge field theory that describes the color-charge interactions between gluon and quark partons.However, these partons cannot be studied in isolation.Insofar as QCD is concerned, only compact (fm-size) color-neutral objects can be detected.One says that the the gluons and quarks are confined, but the meaning of this statement is vigorously debated [14,Sec. 5]: it can be argued that confinement and EHM are two sides of the same coin.Our perspective stresses a dynamical picture of confinement, wherein nonperturbative interactions drive changes in the analytic structure of colored Schwinger functions that ensure the absence of color-carrying objects from the Hilbert space of observable states [14,Sec. 5].Such changes are evident in both continuum and lattice analysessee, e.g., Refs.[15][16][17][18].
Being the most abundant and stable hadron in the Universe, the proton has played a major role in the scientific endeavor.It debuted a little over a century ago, in Rutherford's experiments that exposed the substructure of the atom [19][20][21][22]; and later in the 1950s, when it was revealed that the proton also has a finite size [23].Subsequently, deep inelastic scattering measurements on the proton provided empirical evidence for the existence of quarks [24].
In addressing EHM and confinement, the proton has also been prominent.For instance, it is argued that the proton mass, m p ≈ 1 GeV, about 2000-times that of the electron, is a natural mass scale for visible matter [25,26].To this value, by itself, the Higgs boson (HB) only contributes ≈ 1%, a fact highlighted in Fig. 1.The largest Fig. 1: Mass budgets of the proton, ρ meson and different pseudoscalars.Three types of contributions play a role: EHM, HB and the EHM + HB interference.Clearly, proton and ρ are dominated by EHM (above 90%), heavy quarkonia by HB (85−90%), and finally, pions and kaons by the interplay between EHM and HB mass generation (95% and 80%, respectively).These budgets are drawn at a resolving scale of ζ = 2 GeV, but the images are little changed if one chooses instead to use renormalizationpoint-independent quantities.fraction, some 94%, is a definitive expression of EHM; and the remaining ∼ 5% results from constructive EHM + HB interference [27].
In contrast to protons, π-and K-mesons, pions and kaons, are light, despite also being hadron bound states, and somewhat more compact [29].At some level of approximation, they carry the strong force between nucleons [41]; and the existence of our known Universe requires that they be light compared to m p .At first glance, pions and kaons are the simplest QCD bound states; in fact, they are often drawn as two-body systems, viz. a quark and antiquark, somehow held together.However, this picture is simplistic [42,43].
In reality, pions and kaons are Nature's most fundamental Nambu-Goldstone (NG) bosons, which emerge as a consequence of dynamical chiral symmetry breaking (DCSB) in the SM [44,45].DCSB is a corollary of EHM.In the absence of HB couplings into QCD, i.e., in the chiral limit, π-and K-mesons would be massless and structurally indistinguishable.
Restoring HB couplings, pions and kaons acquire their unusually low masses: m π ≈ 0.15 m p and m K ≈ 0.53 m p ; and exhibit structural dissimilarities.As displayed in Fig. 1, alone, HB couplings into QCD generate only a small fraction of the π and K masses.In these cases, however, the bulk of the bound-state masses are generated by constructive EHM + HB interference.This contrasts starkly with the ρ-meson: constituted from the same valence degrees-of-freedom as the pion, m ρ ≈ 0.82 m p ≈ 5.6 m π and the associated mass budget is qualitatively indistinguishable from that of the proton -Fig. 1.
It should now be clear that in attempting to discover the origin of visible mass, a focus solely on proton structure is inadequate.It is imperative to broaden the goal and provide a simultaneous and unifying understanding of Nature's most fundamental Nambu-Goldstone bosons, viz.pions and kaons.Furthermore, every pseudoscalar meson would be a Nambu-Goldstone boson if it were not for HB couplings into QCD.This means that exploring and explaining the properties of the entire collection of pseudoscalar mesons offers a unique opportunity for elucidating the observable expressions and impacts of Nature's two known mass generating mechanisms and the interference between them -compare, e.g., Figs. 1 and 18 in Ref. [14].No wonder, then, that with high-energy, high-luminosity facilities becoming a reality, much experimental and theoretical attention is shifting toward the study of pseudoscalar mesons [30][31][32][33][34][35][36][37][38].
Hereafter, we sketch possibilities offered by studies of pseudoscalar mesons to provide insights into the emergence of hadron mass and structure.The vehicle for this discussion is provided by continuum Schwinger function methods (CSMs), an approach that has delivered significant progress, especially in the past decade [14,26,43,[46][47][48][49][50][51][52][53][54].The manuscript is organized as follows: Sec. 2 discusses emergent phenomena in QCD and their study using CSMs; pseudoscalar meson distribution amplitudes and functions, and their connection with EHM are addressed in Sec.3; a complementary perspective, obtained via electromagnetic and gravitational form factors, is provided in Sec.4; Sec. 5 follows, using light-front wave functions (LFWFs) and generalized parton distributions (GPDs) to address some questions concerning the three-dimensional structure of pseudoscalar mesons; and finally, a perspective is drawn in Sec. 6.

Emergent phenomena in QCD
Recall the QCD Lagrangian density: here q f denote the f -flavored quark fields, with current quark masses m f ; A a µ are the gluon fields and λ a the generators of SU(3) in the fundamental representation; and g is the unique QCD coupling [54].Plainly, L QCD looks very similar to the Lagrangian density of quantum electrodynamics (QED), except for the underlined piece.This special term gives QCD its non-Abelian character and underlies asymptotic freedom [55][56][57].It also plays a key role in all non-perturbative facets of QCD, including confinement and EHM.
Decades of analyses have crystallized in a now widely accepted picture of dynamical gluon mass generation in QCD, which owes to a Schwinger mechanism [52][53][54].Amongst other important implications, this translates into the physical image of a gluon developing a running mass, m g (k), which is large in the infrared and decreases monotonically toward its perturbative (massless) limit as k 2 increases [58].Likewise, in the matter sector, the nonperturbative phenomenon of DCSB generates a dressed lightquark mass that is equally large in the infrared, whether there is a Higgs mechanism or not [59,60].Thus, as illustrated in Fig. 2, dynamical mass generation is manifest in both the gauge and matter sectors.In Fig. 2, the left panel highlights similarities and differences between m g (k) and the chiral limit (m f ≡ 0) quark mass function; notably, both saturate in the infrared at a value ∼ 400 MeV.The right panel displays the influence of DCSB on the different quark flavors: although its effects are visually more conspicuous in the light sector, the magnitude of the vertical displacement above the current-mass at infrared momenta is roughly the same in each case.
The generation of mass in the gauge and matter sectors is inextricably linked with the behavior of the QCD running coupling [54] and recent advances in experiment and theory are leading to a solid understanding of this fundamental quantity.Indeed, exploiting continuum advances in analyses of QCD's gauge sector, a unique QCD analogue of the Gell-Mann-Low charge [61], used widely in QED, is now available [62].Moreover, using the best existing continuum and lattice results for low-order gauge-sector Schwinger functions, Ref. [58] delivered a parameter-free prediction for this process-independent (PI) charge, α(k 2 ) -see Fig. 3.Even at far-infrared momenta, the result possesses an uncertainty that is just 4%: α(k 2 = 0) = 0.97(4)π.Such accuracy is quite remarkable in an ab initio calculation of an essentially nonperturbative quantity in QCD.
The characteristics and virtues of α(k 2 ) are canvassed elsewhere -see, e.g., Refs.[14,54], and they include: the absence of a Landau pole; a match with the processdependent charge built using the Bjorken sum rule [54,63,64]; and an unambiguous infrared extension (completion) of the perturbative running coupling.Fig. 3: Process-independent charge [58] compared with empirical extractions [54,63] of the process-dependent charge defined via the Bjorken sum rule [64].
Each of the unique features described above contributes to the formation of observable color-neutral bound states and determination of their properties.In this connection, as already noted, charts of pseudoscalar meson properties provide clear windows onto such emergent phenomena in QCD.We now proceed to illustrate this fact using CSMs, implemented via QCD's Dyson-Schwinger equations (DSEs) [14,26,43,[46][47][48][49][50][51][52][53][54].The DSEs may be described as QCD's quantum equations of motion and have been instrumental in exposing emergent phenomena in the strong interaction.
Regarding mesons, all structural information is encoded in the Poincaré-covariant Bethe-Salpeter wave function (BSWF): χ P (k; P ) = S(k + P/2)Γ P (k; P )S(k − P/2) . (4) Here, S(p) refers to the dressed quark propagator and Γ P to the meson Bethe-Salpeter amplitude (flavor indices omitted); P is the total momentum of the meson P, and k is the relative momentum between the valence-quarks.(In a Poincaré-invariant treatment, the choice of k is practitioner dependent and no observable can depend on its definition.) The quark propagator can be expressed as follows: where Z(p 2 ) is the quark wave function renormalization function and M (p 2 ) is the quark (running) mass function.(Herein, we largely omit a discussion of renormalization.It is nevertheless worth noting that all renormalization point dependence of the quark propagator is contained in that of Z(p 2 ).M (p 2 ) is renormalization-point independent.)The mass functions associated with a physical range of quark current-masses are drawn in Fig. 2.
The structure of the Bethe-Salpeter amplitude (BSA) depends on the quantum numbers of the meson.For a pseudoscalar meson: Γ P (q; P ) = γ 5 [iE P (q; P ) + γ • P F P (q; P ) + γ • qG P (q; P ) + q µ σ µν P ν H P (q; P )] .(6) The BSWF has an analogous decomposition, in terms of functions we will write as E χ , F χ , G χ , H χ .In the meson rest-frame, G χ , H χ correspond to L = 1 components.At pion and kaon masses, one measure of their strength indicates that L = 1 components provide ≈ 20% of the canonical normalization [81].Under Poinaré transformations, there is mixing between the terms in the wave function analogue of Eq. ( 6).Consequently, pseudoscalar mesons contain nonzero quark orbital angular momentum in any reference frame.
The 1-body kernel, K (1) , is connected to higher-order Schwinger functions (such as the quark-gluon vertex) ad infinitum, as well as with the two-body kernel, K (2) .Thus, these DSEs form an infinite system of coupled integral equations [82].Consequently, a sound treatment of QCD bound states demands a systematic and symmetry-preserving truncation scheme [83][84][85].
Any such truncation will ensure compliance with the Goldstone theorem [86,87], which is most fundamentally expressed as an equivalence between the one-body quark propagator problem and the two-body meson problem in the pseudoscalar channel [44,45]: where f 0 P is the pseudoscalar meson decay constant and the super/subscript "0" indicates that the quantities are obtained in the chiral limit.In addition to exhibiting the equivalence between the one-and two-body problems, Eq. ( 9) states that Nambu-Goldstone bosons exist if, and only if, DCSB is realized.The leading-order symmetry-preserving truncation, dubbed rainbow-ladder (RL), is sufficient to guarantee the relevant symmetry principles [83,84], and thus describe the properties of ground-state pseudoscalar mesons.Ways to systematically improve upon this leadingorder kernel are known; so, there are increasingly more sophisticated truncationssee, e.g., [88][89][90][91][92][93][94].
3 Distribution amplitudes and functions

Parton distribution amplitudes
Parton distribution amplitudes (DAs) characterize the probability that a nominated parton within a hadron carries a light-front fraction x of that hadron's total momentum, P .Each is a particular one-dimensional projection of the hadron's light-front wave function (LFWF), which itself is the closest analogue in a quantum field theory to the wave function familiar from quantum mechanics [118].Such DAs play a crucial role in the description of hard exclusive processes [119][120][121].
For a pseudoscalar meson P = q h, where q, h specify valence-quark flavor, a nonperturbative extension of the leading-twist DA, φ q P (x; ζ), may be obtained by projecting the meson's Poincaré covariant BSWF onto the light-front [99]: where with n a light-like four vector satisfying n 2 = 0 and n • P = −m P in the meson's rest frame; tr CD indicates a trace taken over color and spinor indices; and we have explicitly indicated the regularization scale, Λ, and ζ, the point whereat the quark propagators and meson Bethe-Salpeter amplitude are renormalized.As written in Eq. ( 10), φ u P is dimensionless and normalized to unity, i.e., 1 0 dxφ q P (x; ζ) = 1.Notably, in a framework that preserves the multiplicative renormalizability of QCD, both f P , φ q P (x; ζ) are independent of the renormalization point.This means that ζ need not be identified with the ERBL evolution scale [119][120][121].That scale, ζ ev , is a somewhat amorphous measure of the degree of collinearity of the partons in the relevant Fock components of the LFWF, i.e., an upper bound on the k 2 ⊥ values of partons contributing to the subject DA.
The kinematics in Eq. ( 10) are set such that the expression refers to the q quark inside the meson P (q-in-P).The corresponding h DA is readily obtained as This guarantees momentum conservation and entails φ h P (x; ζ) = φ q P (x; ζ) in the isospin symmetric limit, viz.m q = m h .We assume the limit m u = m d hereafter.
The Mellin moments of the DA are defined as usual: and, following from Eq. ( 10), these can be computed using The asymptotic behavior of the leading-twist DA has been known for more than forty years [119][120][121]: Over the years, however, it has become clear that this form is unrealistic, in the case of π and K, on any domain accessible to terrestrial experiments.It is thus imperative to compute the DA at experimentally accessible scales.
Such efforts, using continuum and lattice methods, are canvassed elsewhere [43].They have concluded that, at all scales accessible to existing and foreseeable experiments, φ π,K (x; ζ) is a broad concave function, viz.strongly dilated and flattened in comparison with φ asy (x).This is illustrated in the left-panel of Fig. 4, which depicts CSM predictions for the π, K DAs [66] and contrasts them with φ asy (x).
The dilation of φ u π,K (x; ζ) is a manifestation of EHM; and the DA of the pion -Nature's lightest hadron -is the most dilated of all.Owing to flavor symmetry breaking, expressed in the QCD Lagrangian by a large difference between the s and u quark current-masses, the maximum of φ u K (x; ζ) is shifted slightly to x ≈ 0.4.This 20% relocation is a statement that EHM dominates, but HB modulations are beginning to be felt, as may be seen by considering that Turning attention to heavy q q pseudoscalar mesons, the pattern is reversed.Owing to the large HB-induced quark current masses, the DAs of such systems are compressed/contracted [122], something seen in the right-panel of Fig. 4.These DAs become increasingly narrow, with greater x = 1/2 peak height, as the current masses become larger.Indeed, in the limit of infinitely heavy quarks: Considering, just for illustration, a pure ss =: η s pseudoscalar meson, with m ηs ≈ 0.69 GeV, one finds φ s ηs (x, ζ 2 ) ≈ φ asy (x).Thus, ss systems define a boundary, whereat EHM and HB mass generating effects are of roughly equal importance.Numerical results from the simulations of lattice-regularized QCD (lQCD) confirm this CSM prediction [123].Moreover, CSM studies of the η-η ′ complex [105] show that the lightand s-quark component DAs of the η ′ -meson, especially, match this expectation.
The analysis of heavy-light mesons provides additional information about the interplay between strong and weak mass generation [101].The DAs of D and B mesons, compared with the π, K cases, are shown in Fig. 5. Clearly, as the size of flavor symmetry breaking increases and, thus, the HB impact on the heavier quark -see Ref. [14,Fig. 18], the DA distortion becomes more pronounced.At fixed light-quark mass, then with increasing heavy-quark mass, the location of the DA peak moves toward a minimum value x min ∈ (0, 0.5).That x min ̸ = 0 is a special feature of heavy + light systems [124].One may quantify this peak relocation by computing values of 1 0 dx xφ(x; ζ 2 ), which is a DA-weighted momentum fraction: π, 0.5; K, 0.48; D, 0.32; B, 0.19; and 0.12(1) in the case of one infinitely heavy quark.
For ground-state pseudoscalar mesons with masses m P ≤ m asy ≈ m ηs , the following parametrization is efficacious: Refs.[66,101].
whereas this one works well on m P ≥ m asy : In both cases, n P ensures unit normalization and ρ P , γ P are interpolation parameters.
These simple forms enable one to express the endpoint behavior predicted by QCD; dilation/compression of the DAs, via ρ P ; and skewing in flavor asymmetric systems, via γ P .Notably, both forms reproduce φ asy (x) in appropriate circumstances, viz.ρ P → ∞, γ P → 0, and that is why they have a common boundary of applicability.At the other extreme, i.e., ρ P → 0, γ P → 0, one finds φ light P (x; ζ) → φ SCI (x) = 1 and φ heavy P (x; ζ) → φ ∞ (x).The former is the broadest possible distribution and corresponds to that produced using a symmetry-preserving treatment of a vector ⊗ vector contact interaction [125], whereas the latter is the narrowest.Finally, the combination of numerous analyses, Refs.[66,101,105,106] enables us to determine the interpolation parameters shown in the Table 1.Associated low-order moments are also listed.

Parton distribution functions
Complementing DAs, parton distribution functions (DFs) play a key role in the description of hard inclusive processes.A given DF, p P (x; ζ), is a number density, so that p P (x; ζ)dx is the number of p partons carrying a light-front fraction between x and x + dx of the total momentum of hadron h at a resolving scale ζ [126].
The DFs of P are accessible via the associated forward Compton amplitude, γP → γP [126].Detailed considerations of that amplitude have led to the following expression Table 1: Interpolation parameters for the pseudoscalar meson DAs, to be used in Eqs. ( 16), (17), as appropriate.Entries above the horizontal line correspond to mesons whose masses are less than m asy .The complement lies below this line.In the η − η ′ case, the superscript refers to the light l = u/d and strange s components of its wave function [105].for the q-in-P valence-quark DF [102][103][104]: where The result from Eq. ( 18) does not depend on the scale at which the quark propagator and meson Bethe-Salpeter amplitude are renormalized [66]; so, what is the meaning of ζ H ?
In developing the answer to this question, it is important to note that the following identities are readily verified: The statements in Eq. (19a) express baryon number conservation.They must be valid, independent of the value and meaning of ζ H . On the other hand, those in Eq. (19b) mean that valence-quark degrees-of-freedom carry all the hadron's light-front momentum at the scale ζ H .This is a principal reason behind the identification of ζ H as the hadron scale [103].We will see that Λ QCD ≲ ζ H < m p .A hadron's LFWF provides a bridge between its valence-quark DAs and DFs.This connection is best introduced by illustration; so, suppose one has a two-body system  25) (dot-dashed green curve).Dashed cyan curve: scale-free function in Eq. ( 22).described by the following LFWF: where M is a mass whose size is assumed to be set by EHM and n ψ is the normalisation constant.For δ = 0, this LFWF exhibits the large-k 2 ⊥ scaling behaviour of a leadingtwist two-body wave function in QCD [121,Eq. (2.15)].(True QCD wave functions also include ln k 2 ⊥ scaling violations.This is mimicked by δ ≳ 0.) The valence-quark DF is obtained as This is a scale-free function, i.e., it is not explicitly dependent on the mass-scale in the LFWF.However, the presence of that scale and, hence, EHM is manifested in the denominator structure 1 + x(1 − x), which, as evident Fig. 6, introduces a dilation with-respect-to the numerator function alone: It is worth noting that one may recover this EHM-insensitive function by changing 20) and then taking ϵ → 0.
The meson's DA is obtained via where extraction of the leptonic decay constant, f P , means 1 0 dx φ q P (x, ζ H ) = 1.With δ = 0, this integral is ln-divergent and that explains the renormalization constant in Eq. (10).For the purpose of our illustration, we define the result by expressing the integral as a series in δ on δ ≃ 0, discarding the 1/δ piece that characterizes the ln-divergence, then taking the limit δ → 0. This procedure yields Here, the only explicit dependence on M is contained in the decay constant; so, the DA is seemingly independent of this mass-scale.However, as with the DF, EHM is expressed in the 1 + x(1 − x) term.(To make these things readily apparent, we have replaced special-functions at fixed arguments by their numerical values.)Once again, akin to the DF, one may recover the asymptotic DA in Eq. ( 14) by making the replacement in Eq. ( 20) and then taking ϵ → 0. Consider now the following rescaled DA: Figure 6 compares [ φq P (x, ζ H )] 2 with the DF derived from the same LFWF.The two functions are practically indistinguishable: mathematically, the L 1 -difference between these curves is just 2.7%.Following the same procedure with the asymptotic DA in Eq. ( 14) yields exactly the scale-free DF in Eq. ( 22).This highlights a simple fact.Namely, M 2 [1+x(1−x)] → M 2 in Eq. ( 20) produces a factorized LFWF: ; and whenever this is a good approximation for quantities obtained by integration -it need not be pointwise precise, then Such factorized representations are known [69] to be a good approximation for groundstate mesons in which either the valence quarks are mass degenerate or EHM leads to significant suppression of HB-induced flavor symmetry violation, e.g., kaons.Its value in treating other mesons, such as excited and heavy + light states, is currently being explored.

Hadron scale distributions
The reliability of Eq. ( 26) for pions and kaons has been exploited to good effect [65,66,72,73,78,127]. The first step is to associate the DA obtained using Eq. ( 10) with [left] π, K distributions at ζ H , in constrast with the scale-free profile q sf [66].[right] Analogous comparison for u π u p , d p [73].The π, K uncertainty bands stem from that associated with the corresponding DAs, whereas those of the proton express a 10% variation in the values of the low-order moments ⟨x − x q 0 ⟩ q P and ⟨(x − x q 0 ) 2 ⟩ q P , where x u,d 0 ≈ 0.35, 0.32 corresponds to the maximum of the xu p and xd p distributions, respectively.In each panel, the vertical lines indicate the maximum of a given u-quark DF.
the hadron scale, ζ H .This places it on the same level as the DF calculated using Eq. ( 18).(Recall that both expressions produce distributions that are independent of the propagator and Bethe-Salpeter amplitude renormalization scale.) Following this procedure and working with the DAs drawn in Fig. 4, one obtains the π and K DFs drawn in Fig. 7 -left panel.They are noticeably dilated in comparison with q sf (x).The DFs produce the following low-order Mellin moments: N.B.At ζ H , (u K +s K )/2 is a symmetric distribution, so its moments can be compared directly with those of the pion DF: evidently, HB modulation does not significantly affect the dilation of the kaon DF.Using the scale-free DF, these low-order moments are, respectively, 1 2 , 2 7 = 0.286.It is worth stressing that, despite the EHM-induced dilation, each DF is compatible with QCD constraints [126,128,129].
An interpolating DF parametrization that simultaneously expresses the EHMinduced dilation, soft endpoint behavior, and skewing, when present, is provided by the following function: The DFs in Fig. 7 -left panel are reproduced with For some purposes, such as the calculation of kaon fragmentation functions [78], one may use an alternative, practically equivalent form: ρK = 0.062, γK = 13.83.Calculation of proton DFs is described in Ref. [72,73,80].Such analyses yield the following light-front momentum fractions: Contrary to the pion, and despite also being composed of light quarks, the proton momentum distributions associated with the different flavors are not the same, even accounting for the 2 × u : 1 × d ratio.This is a manifestation of SU (4) spin-flavor symmetry breaking in the proton wave function, which may be attributed to the emergence of strong, nonpointlike diquark correlations [130].Furthermore, as highlighted by Fig. 7 -right panel, the valence parton DFs in the pion and proton have markedly different profiles.This owes partly to the difference in the number of valence degreesof-freedom; but it is also an expression of EHM, with the pion valence-quark DF being the most dilated amongst all hadrons.For instance, the η ′ meson DF profile is much less dilated.

Evolved distributions
In order for data to be connected with DFs, the experiments should involve energy and/or momentum transfers (far) in excess of m p , i.e., be conducted on a kinematic domain for which QCD factorization is valid.Consequently, before comparisons with such data can be made, the hadron scale DFs must be evolved to the energy scale appropriate to a given experiment.That can be accomplished using the all-orders (AO) evolution scheme explained elsewhere [76], which is a particular realization of DGLAP evolution [131][132][133][134] that has proved efficacious in numerous applicationssee, e.g., Refs.[66-69, 72-74, 77-80, 135].
Regarding DF Mellin moments, the AO scheme provides closed algebraic relations between them.For instance, the moments of any given valence-parton DF are related as follows: where γ n qq are the appropriate 1-loop anomalous dimensions [131][132][133][134].This identity states that all Mellin moments of the DF at any ζ > ζ H are completely determined by the valence-quark momentum fraction at this scale, so long as all moments are known at the hadron scale.Recall that, by definition, ζ H is the scale at which valence degrees-of-freedom carry all the hadron's properties.This entails that glue and sea DFs are identically zero at the hadron scale: It has thus far been found that the same value of ζ H serves well for all hadrons.Its actual value is immaterial and need not be specified.Notwithstanding that, practical analyses of lQCD results indicate that ζ H ≈ 0.35 GeV [79]; and using the PI charge described above, one predicts [66, Eq. ( 15)]: It is worth stressing that, under evolution, glue and sea DFs are nonzero ∀ζ > ζ H ; moreover, even on ζ ≃ m p , a significant fraction of a given hadron's light-front momentum is lodged with glue and sea.This entails that, even without recourse to "intrinsic charm" [136], roughly 1.5% of the hadron's momentum is lodged with the c-quark sea at ζ ≃ 1.5 m p .These features are illustrated for the pion in Fig. 8.The left panel displays the valence, glue, and sea DFs calculated in Ref. [66].Notably, the parameter-free CSM prediction for the glue DF agrees well with a recent lQCD computation [127,137].At this scale, ζ = ζ 2 = 2 GeV, referred to the light-front [66]: valence degrees-of-freedom carry 48(4)% of the pion's momentum; glue, 41(2)%, and four-flavor sea, 11(2)%.(In Ref. [66], quark current-mass effects were not included in the evolution equations.More recent analyses have introduced mass thresholds [73,76].) Figure 8 -right depicts the CSM predictions at ζ 5 = 5.2 GeV, i.e., the scale of Ref. [139, E615].There is excellent agreement with the analysis of that data described [left] Ratio p K /p π for different parton species: u valence-quark (solid blue), gluon (dashed green), and sea (dot-dashed red).Data [141] and lQCD (grey band) [142] correspond to the ratio valence ratio u K /u π .[right] CSM predictions for sK (x), u K (x) and u π (x).The error bands associated with the kaon correspond to a ±10% variation of ζ H . in Ref. [140].At this scale, referred to the light-front [66]: valence degrees-of-freedom carry 41(4)% of the pion's momentum; glue, 45(2)%, and four-flavor sea, 14(2)%.The lQCD computation from Ref. [138] is also displayed in Fig. 8 -right.That study used a novel techique for extracting DF pointwise behavior from a Euclidean lattice.Further discussion of these and related points can be found in Refs.[70,71,79].
Regarding the kaon, empirical information is scarce.Only eight points are available and those relate solely to the valence-quark ratio u K (x; ζ 5 )/u π (x; ζ 5 ) [141]. Figure 9 left displays the CSM prediction for this ratio compared with experimental data and a lQCD result [142]: plainly, all results are compatible.Actually, compared individually, the CSM and lQCD results for u K (x; ζ 5 ), u K (x; ζ 5 ) are quite different.Evidently, therefore, the ratio is a forgiving measure and data on the individual DFs would provide a far keener tool for discriminating between pictures of kaon structure.This panel also displays CSM predictions for analogous glue and sea ratios.The right panel of Fig. 9 contrasts the valence-quark DFs within the kaon and pion.In this calculation [66], recognizing that gluon splitting must produce less heavy s + s pairs than light u + ū pairs and heavy s quarks should produce less gluons via bremsstrahlung, quark current-mass threshold factors were included in the evolution kernel.As a consequence, when compared with mass-independent evolution results, low-order Mellin moments of the s-in-K valence DF are increased by 4.8(8)% and the glue moments are commensurately smaller.
Before closing this section, we list CSM predictions for low-order Mellin moments of valence parton DFs in the pion and kaon: These predictions may be viewed as benchmarks for phenomenology.Existing approaches to fitting relevant data typically see such phenomenology place too much momentum in the sea with the cost paid by the valence fraction [70].
Finally, comparisons between pion and proton DFs are drawn in Ref. [73].Notably, in all cases, QCD-connected CSM predictions for evolved DFs comply with SM constraints on the large-x behavior.

Electromagnetic and gravitational form factors 4.1 Hard-scattering formulae and scaling violations
Electromagnetic form factors (EFFs) also provide opportunities for the examination of diverse aspects of hadron internal structure.Obvious examples are the charge and magnetization distributions, but it goes much further than that.For instance, such EFFs of pseudoscalar mesons present an ideal platform for testing fundamental QCD predictions, since rigorous connections have been drawn between them and the DAs discussed above [119][120][121].
Focusing first on pseudoscalar meson elastic and transition electromagnetic form factors (EFFs and TFFs), QCD predicts the following behavior.
Transition: γ * (Q)γ → P 0 , considering any q q component of P 0 , where: f q P is the q q-component contribution to the pseudovector projection of the meson's wave function onto the origin in configuration space, i.e., a leptonic decay constant; and e q is the electric charge of quark q.The complete transition form factor is obtained as a sum over the various q q subcomponent contributions: where ψ q P is a flavour weighting factor originating in the meson's wave function.It is made plain by Eq. ( 36) that QCD is not seen in EFF Q 2 -scaling, but in the violations of scaling that reveal the character of the running coupling and evolution of the DA.Regarding TFFs, scaling violations are also evident in the leading order result, Eq. ( 38), through the evolution of the DA.Importantly, the absolute magnitude of either the EFF or TFF on the ultraviolet domain is set by the leptonic decay constant of the meson involved.This quantity is an order parameter for DCSB; hence, a measure of EHM.In addition, when considering the neutral pion TFF in the neighbourhood of the chiral limit, one has [143][144][145][146] Thus EHM sets the infrared scale as well; and deviations from this result for other (heavier) mesons are a measure of EHM + HB interference.

The longstanding question is:
At what value of Q 2 ≫ m 2 p do Eqs.( 36), (38) begin to serve as good approximations, i.e., how hard is hard for exclusive processes?Forty years of experiment and theory have shown that if φ asy (x) is used in these equations, then that domain is beyond the reach of terrestrial experiments.This opens up the possibility that existing and foreseeable EFF and TFF measurements might be interpreted as placing constraints on the 1/x moment of meson DAs.Such potential is the subject of much debate and analysis -see, e.g., Refs.[34,36,[108][109][110]147].The η and η ′ TFFs are of additional interest because they can expose observable consequences of the non-Abelian anomaly and topological effects within hadrons [105,148,149].
with a similar expression for F h P (Q 2 ), where Q is the incoming photon momentum, o,i = −m 2 P , m P is the meson mass.The calculation also requires quark propagators, S f , f = q, h, which, consistent with Eq. (41b), should be obtained from the rainbow-truncation gap equation; the meson Bethe-Salpeter amplitude, Γ P , computed in rainbow-ladder truncation; and consistent unamputated dressed-quark-photon vertices, χ f µ .Parameter-free CSM predictions for the π + and K + EFFs [107,114], calculated from Eq. ( 41), are displayed in Fig. 10.Plainly, where precise data are available [152,153], they agree with the CSM predictions.Importantly, unlike many other approaches, CSMs deliver QCD-connected predictions on the entire domain of spacelike momentum transfer.Timelike momenta are also accessible -see, e.g., Ref. [158], but this requires careful treatment of resonance contributions to the photon-quark vertex.
Considering the low-Q 2 domain, the CSM results deliver predictions for charged pion and kaon radii via the usual definition: Using the curves in Fig. 10, one obtains r π = 0.64(2) fm, r K ≈ 0.9 r π .These values are consistent with modern determinations [159].
Given the range accessible to CSM analyses, it is possible to develop an answer to the question "how hard is hard?"Focusing first on the pion, consider Fig. 10 left.Precise higher-Q 2 data are available [152,153]: in comparison with the CSM prediction, one finds χ 2 /datum = 1.0.Extending beyond the range of extant data, the CSM prediction follows a monopole function, whose scale is fixed by the pion charge radius, until Q 2 ≈ 6 GeV 2 .Thereafter, the CSM prediction breaks away from the simple scaling result, trending below with a separation that grows as Q 2 increases.This onset of scaling violation is the signal for QCD in hard exclusive scattering; and based on the anticipated precision of forthcoming data, experiments that probe above Q 2 ≈ 9 GeV 2 will be sensitive to this signal [34,36,43].
Figure 10 -left also includes the result obtained from Eq. ( 36) by using a dilated pion DA of the type displayed in Fig. 4 (dashed purple curve) and that produced by φ asy (x) in Eq. ( 14) (dotted black curve).Evidently, the broadened DA provides semiquantitative agreement with the CSM prediction.Indeed, the quantitative difference between these curves may be explained by a combination of higher-order and -twist corrections to Eq. ( 36) in perturbative QCD on the one hand, and shortcomings in the rainbow-ladder truncation, which predicts the correct power-law behaviour for the form factor but not precisely the right anomalous dimension in the strong coupling calculation, on the other hand.Empirical support for the broadened pion DA is also found in analyses of pion + proton Drell-Yan data [77].
The CSM prediction for the charged kaon EFF is drawn in Fig. 10 -right.At this time, on the entire domain, precise data are lacking [29,159], but that is expected to change in the foreseeable future [34,36,43].Nevertheless, the figure reveals that scaling violations should be visible in the charged kaon EFF on Q 2 ≳ 6 GeV 2 .
transition form factors [105,110].Hence, the following comparisons have weight: sources Refs.[160][161][162] Refs.[160,161,163] Refs.[ Plainly, the BaBar Collaboration data [162] are not compatible with the CSM prediction, whereas the Belle Collaboration data [163] match well.Going further and focusing on data at Q 2 > 10 GeV 2 , one finds χ 2 /datum = 4.14 [162,BaBar] and χ 2 /datum = 0.64 [163,Belle].These comparisons suggest to us that Eq. ( 38) is confirmed by the bulk of existing data and, hence, such data support a picture of the pion DA as a broad, concave function at experimentally accessible probe momenta.This perspective may be tested by new data [30].The γ * (Q)γ → {η c , η b } TFFs are drawn in Fig. 11 -right and compared with their π 0 analogue, all normalized to unity at Q 2 = 0.The parameter-free CSM prediction agrees well with the Ref. [171,BaBar] data.There is no data on γ * (Q)γ → η b , but a result is available from an analysis made using a non-relativistic QCD effective field theory [170].This result and the CSM prediction are in excellent agreement.
The impact of HB couplings into QCD on these TFFs is readily apparent in Fig. 11 -right.As the HB generated quark current-mass increases, the associated TFF falls more slowly.This effect is manifest in the TFF radii: r π 0 ≈ 0.66(2) fm, r ηc /r π 0 ≈ 0.25, r η b /r π 0 ≈ 0.06.Consequently, with increasing current-mass, the boundary of the domain upon which the hard-QCD limit, Eq. ( 38), delivers a good approximation is pushed ever deeper into the spacelike region.
Notably, the successes of CSMs in describing all neutral pseudoscalar meson TFFs [105,[109][110][111], have made the approach a credible contributor in matters related to the anomalous magnetic moment of the muon [8].

Gravitational form factors
The interaction of a pseudoscalar meson with a J = 2 probe is characterized by the current: which corresponds to the in-meson expectation value of the QCD energy-momentum tensor [172].Here θ P 2,1 are the meson gravitational form factors (GFFs) associated with the mass and pressure distributions.Symmetries impose the following relations: The first identity is a statement of mass normalization; the second is connected with a soft-pion theorem [172,173], a corollary of EHM; and the third is a direct consequence of energy-momentum conservation, i.e., The GFFs defined in Eq. ( 46) are the result of adding up the individual contributions of each type of parton.Working at ζ H , however, only fully dressed valence quarks play a role and F P = F q P + F h P (F P = θ P 2,1 , cP ).Consequently, the GFFs maybe be extracted from an analogue of Eq. ( 41), obtained by replacing the photon-quark vertex by the "graviton"-quark vertex [114].with F π : θ 2 shows the least rapid decay with Q 2 , F π falls faster, but θ 1 decays most rapidly.Naturally, therefore, the corresponding radii are ordered as follows: in other words, the mass distribution is more compact than the charge distribution (r π θ2 /r π F ≈ 0.74), which is, in turn, tighter than the pressure distribution (r π F /r π θ1 = 0.79).These patterns are also consistent with available lQCD results [174] -see Fig. 12 -right, and have been confirmed empirically [75,175].
It is worth noting here that by exploiting general physical constraints on the form of the pion DF [70,71,79], the following bounds can be established: The results above are consistent with these limits.Kaon FFs exhibit similar profiles [114].Naturally, being heavier than the pion, the kaon is more compact when judged by the measures discussed herein.Quantitatively, averaging all relevant radii, one finds r K /r π = 0.85 (6).
Gravitational form factors are also accessible via generalized parton distributions (GPDs) [117,176,177], which themselves will be discussed below.However, following that route, the pressure distribution, θ 1 , suffers from the so-called D-term ambiguity [178].This is avoided when one calculates the form factors directly from the analogue of Eq. ( 41) and this makes θ 1 especially interesting.The CSM analysis in Ref. [114] predicts θ π 1 (0) = 0.97 , θ K 1 (0) = 0.77 (10) , both in agreement with estimates made using chiral effective field theory [172].The deviations from Eq. ( 47) are a measure of HB modulation of EHM in the systems considered: the magnitudes match expectations.
Importantly, too, pressure and shear profiles -p P (r) and s P (r), respectively, can be obtained from θ 1 as follows [172]: where 2E(∆) = 4m 2 P + ∆ 2 and j 2 (z) is a spherical Bessel function.The pressure is shown in Fig. 13 -left: it is large and positive at small separations; but at some critical value, r c , it changes sign.This point may be interpreted as marking the beginning of the domain on which confinement forces become dominant.The locations are (in fm): r π c = 0.39(1), r K c = 0.26(1), r Ku c = 0.30(1), r Ks c = 0.25 (1).Notably, this qualitative change occurs when the deformation forces, shown in Fig. 13 -right, are maximal.Finally, it is worth highlighting that the meson core pressures are commensurate with those in neutron stars [69,179].
5 Toward a 3-dimensional picture

Light-front wave functions
Light-front wave functions provide probability amplitudes for different parton configurations within a hadron.The following light-front projection of a pseudoscalar meson BSWF connects with a leading-twist LFWF: Associations with the discussion in Sec. 3 are readily made.
Owing to the open dependence on k 2 ⊥ , explicit calculation of ψ u P (x, k 2 ⊥ ; ζ H ) from Eq. ( 52) is more complicated than that of the DA via Eq.(10).The challenges can be overcome by developing perturbation theory integral representations (PTIRs) of the numerical solutions for χ P (k − P ; P P ) -see, e.g., Ref. [115].
Another, simpler approach, which is nevertheless insightful and realistic, can be found in Refs.[69,180,181].In character, it exploits PTIRs, but the representations are less complex.For instance, one may write where M q,h are constituent-quark mass-scales, associated with the infrared size of the dressed-quark mass function -see Fig.  spectral density; and n P ensures canonical normalization.Introducing two Feynman parameters, (α, v), the BSWF can be re-expressed as with Inserting these expressions into Eq.( 52), then after a series of algebraic manipulations, the Mellin moments can be expressed as: with σ ⊥ = k 2 ⊥ + Ω 2 P .Capitalizing now on uniqueness properties of Mellin moments, one may immediately conclude that Plainly, the spectral function, ρ P (w), must play a key role in determining the LFWF's profile and, hence, the related DAs and DFs.Alternately, as highlighted by the discussion in Sec. 3, if one has reliable results for a given meson's DA and/or DF, then those results can be used to determine the spectral density.This is the procedure followed in Refs.[69,181] and the resulting LFWFs are displayed in Fig. 14.
Looking deeper into Eqs.( 54)-( 58), some interesting and useful features are revealed.Working in the chiral limit (M q = M h , m 2 P = 0) and setting Λ P = M q , then one finds Ω 2 P = M 2 q .In this case, direct evaluation of all integrals in Eq. ( 56) is possible and, using Eqs.(21a), ( 23), (25), one arrives at: viz. the x-k ⊥ dependences factorize and Eq. ( 26) is recovered.Focusing on Eq. ( 55), it is evident that any violation of factorization is tied to M h −M q ̸ = 0; hence, is small for modest values of [M h −M q ]2/[M h +M q ].Similar conclusions were reached in Refs.[182,183].Consequently, one may reliably proceed with factorized representations for the LFWFs of light pseudoscalar mesons or, if desiring to express factorization violations, a simple wave function of the following form:

Generalized parton distributions
Generalized parton distributions (GPDs) are useful because they connect an array of hadron properties, such as those relating to DFs, EFFs, GFFs, and transverse spatial distributions -see, e.g., Refs.[117,177].In the case of pseudoscalar mesons at the hadron scale, it is especially convenient to work with the valence-quark GPD obtained via the overlap representation [184]: where: 2P = p ′ +p, p ′ , p are the final, initial meson momenta in the defining scattering process; ∆ = p ′ −p, P •∆ = 0; the skewness ξ The overlap representation is only valid within the so-called DGLAP kinematic domain (|x| ≥ ξ).This is sufficient for many purposes and the associated pion and kaon valence-quark GPDs are shown in Fig. 15.
To obtain a complete GPD, one must also have knowledge of the defining matrix element on the complementary ERBL domain (|x| < ξ).Methods have been developed to extend a DGLAP-defined GPD onto the ERBL domain [178,[185][186][187]. Considering Eq. (21a), it is apparent that a meson's hadron-scale DF and forwardlimit GPD are identical: Moreover, the following identities have been established: A meson's complete FF is obtained by properly summing the valence quark and antiquark contributions.In principle, both the direct amplitude approach to form factor calculations, discussed in Sec. 4, and the GPD representations, Eqs.(64), are entirely equivalent.
Novelty is provided by the ξ = 0 impact parameter space (IPD) GPD: This quantity is a true density, which relates to the number of partons within the lightfront at a transverse distance |b ⊥ | from the meson's centre of transverse momentum.
[right] K + and K 0 charge distributions.The charged case has been multiplied by 0.2 so that it shares the same range as the neutral one.
Herein, we limit ourselves to a discussion of π, K charge and mass distributions derived from the IPD GPD.For a q-in-P quark: so, the total meson distributions are: The results are displayed in Fig. 16.The left panel contrasts π + , K + , and π + c , where this last system is a pion-like state with quark and antiquark possessing c-quark masses, all calculated using a LFWF as in Eq. (60a).Consistent with Fig. 12, mass distributions are more compact than charge distributions.Moreover, with increasing meson mass, both distributions are squeezed toward |b ⊥ | = 0 -recall, the area under each curve is the same.This can be quantified by recording the associated radii ratios, r ⊥ M /r ⊥ C : π + = 0.77 (6); K + = 0.80 (7); π + c = 0.97 (3).(The uncertainties were obtained as explained in Ref. [75].) Plainly, pion and kaon ratios are equal, within uncertainties.Notwithstanding that, they do hint at a trend, which is highlighted by the π + c ratio being nearly unity.
Figure 16 -right depicts the K +,0 transverse-plane charge distributions.The K + profile is similar to that of the π + .Regarding the K 0 , destructive interference between d and s distributions leads to a zero at |b ⊥ | ≈ 0.3 fm.Below this value, the distribution is positive, which indicates that the positively-charged s valence-quark is more likely to be found nearer to the centre of transverse momentum than the lighter d quark.

Empirical determination of the pion GPD
It will be very difficult to extract sufficient precise data from deeply virtual Compton scattering experiments to reconstruct meson GPDs [117,177,186].Considering Eqs. ( 63), (64a), it is plain that a different approach is feasible, however.Namely, exploiting the AO evolution scheme and the fact that a factorized representation of the pion LFWF is reliable for integrated quantities, a pion GPD can be recovered from independent data on the pion EFF and valence-quark DF [75].
These ideas were implemented in Ref. [75] as follows.Beginning with existing analyses of pion + nucleus Drell-Yan and electron + pion scattering data [70,139,140,[151][152][153], ensembles of model-independent representations of the pion GPD were developed.They are illustrated in Fig. 17.Using those GPD ensembles and Eq.(64b), data-driven predictions for the pion mass distribution form factor, θ 2 , were obtained.Compared with the pion elastic EFF obtained simultaneously from Eq. (64a), θ 2 is harder: the ratio of the radii derived from these two form factors is r θ2 π /r π = 0.79 (3), which is in accord with Eq. ( 49).The Ref. [75] data-driven predictions for the pion GPD, related form factors and distributions should prove valuable as constraints on theories of pion structure.Improvements to those results are possible if new, precise data relating to pion DFs are obtained and EFF data is secured at larger momentum transfers.An extension to the kaon is currently being explored, but in this case serious impediments are presented by the lack of data and the imprecision of that which is available.

Perspective
Emergent phenomena in QCD are responsible for a diverse array of measurable outcomes.Of primary importance is the fact that the parton degrees-of-freedom used to express the QCD Lagrangian are not directly observable.Only color-neutral systems, seeded by valence partons, can be captured in detectors.This is an empirical definition of confinement.A mathematical definition remains elusive and those being discussed can be contentious.Nevertheless, a growing body of evidence suggests that the source of confinement can be traced to the dynamical generation of running masses in the gauge and matter sectors.Moreover, that the associated masses explain the emergence of ∼ 99% of the visible mass in the Universe.Finally, that these masses, too, eliminate the Landau pole in QCD and the Gribov ambiguity, thereby enabling the calculation of a unique process-independent running coupling, which is everywhere finite and practically momentum-independent in the infrared, so that QCD is effectively a conformal theory at long range.This being the case, then QCD is potentially the first mathematically well defined quantum field theory ever formulated in four-dimensions [54].
The properties of pseudoscalar mesons provide ideal means by which to elucidate these facets of QCD; and we exploited this in providing many explanations and illustrations.Indeed, as both bound-states of valence-quark and valence-antiquark and Nambu-Goldstone bosons in the chiral limit, pseudoscalar mesons provide a clean link to emergent hadron mass (EHM) and a route to exposing its interference with the other known mechanism of mass generation in the Standard model, viz. the Higgs boson.The progress achieved and standpoint of theory today, continuum and lattice, were, to many, unimaginable a decade ago.Nevertheless, controversies remain.
Many of the things we discussed are or will be the focus of experiments at modern and planned facilities.Confirmation of the predictions will go far toward resolving the disputes and validating the EHM paradigm discussed herein.Their complements are also critical, e.g., studies of nucleon structure, of course; but importantly, too, development of an understanding of meson and nucleon excited states.EHM is expressed in each such system; yet, the manifestations need not be everywhere identical.With the advent of high-energy, high-luminosity facilities, it will become possible to map and link the expressions of EHM across the entire field of strong interaction phenomena.Finally, thereby, will science test QCD and decide whether it is one part of an effective field theory of Nature or a theory in its own right, with the potential to guide extensions of the Standard Model.
Declaration of Competing Interest.One of the authors (C.D. Roberts) is a Field Editor for Few Body Systames in the area of Elementary Particle Physics.Apart from this, the authors declare that they have no other known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Funding Information.This work was supported by Spanish MICINN grant PID2022-140440NB-C22; regional Andalusian project P18-FR-5057; National Natural Science Foundation of China grant no.12135007; and STRONG-2020 "The strong interaction at the frontier of knowledge: fundamental research and applications" which received funding from the European Union's Horizon 2020 research and innovation programme (grant agreement no.824093).

Fig. 4 :
Fig.4:[left] π, K DAs in comparison with the asymptotic profile φ asy .[right] Analogous predictions for π, η s , η c .In both panels the uncertainty bands reflect a ±5% variation in the value of the low-order moments ⟨ξ 1 ⟩ P and ⟨ξ 2 ⟩ P , ξ = 1 − 2x; and each thin vertical line identifies the location of the maximum of a given DA -only the K DA maximum is displaced from x = 1/2.The predictions are taken from Refs.[66,122], which used a renormalisation scale ζ 2 = 2 GeV for all Schwinger functions.

Fig. 5 :
Fig.5: π, K, D, B DAs.The uncertainty bands reflect a ±5% variation in the value of the low-order moments ⟨ξ 1 ⟩ P and ⟨ξ 2 ⟩ P , ξ = 1 − 2x; and each thin vertical line identifies the location of the maximum of a given DA.The predictions are taken from Refs.[66,101].