Masses of positive- and negative-parity hadron ground-states, including those with heavy quarks

A symmetry-preserving treatment of a vector ×\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times $$\end{document} vector contact interaction is used to compute spectra of ground-state JP=0±,1±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P = 0^\pm , 1^\pm $$\end{document}(fg¯)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(f{\bar{g}})$$\end{document} mesons, their partner diquark correlations, and JP=1/2±,3/2±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P=1/2^\pm , 3/2^\pm $$\end{document} (fgh) baryons, where f,g,h∈{u,d,s,c,b}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f,g,h \in \{u,d,s,c,b\}$$\end{document}. Results for the leptonic decay constants of all mesons are also obtained, including scalar and pseudovector states involving heavy quarks. The spectrum of baryons produced by this chiefly algebraic approach reproduces the 64 masses known empirically or computed using lattice-regularised quantum chromodynamics with an accuracy of 1.4(1.2)%. It also has the richness of states typical of constituent-quark models and predicts many baryon states that have not yet been observed. The study indicates that dynamical, nonpointlike diquark correlations play an important role in all baryons; and, typically, the lightest allowed diquark is the most important component of a baryon’s Faddeev amplitude.


Introduction
The empirical hadron spectrum is rich [1], even without including the array of recently discovered exotic states [2][3][4][5][6]. Yet, contemporary theory still predicts more states than have been observed. This is especially true when one goes beyond systems comprised of light {u, d}-quarks and includes hadrons seeded by all flavours of valence quarks and/or antiquarks whose lifetime is long enough to produce measurable bound states, i.e. also includes systems containing {s, c, b} quarks and/or related antiquarks. The empiria e-mail: yinpl@njupt.edu.cn b e-mail: phycui@nju.edu.cn c e-mail: cdroberts@nju.edu.cn (corresponding author) d e-mail: jsegovia@upo.es cally "missing" states are being sought in a worldwide effort [7][8][9][10][11].
Explaining the mass and structure of parity partners in hadron spectra is crucial to understanding strong interactions because they would be degenerate if chiral symmetry were not dynamically broken, as discussed for the a 1 -and ρmesons in Ref. [23]. This dynamical chiral symmetry breaking (DCSB) is a corollary of emergent hadronic mass (EHM), the mechanism responsible for the m p ∼ 1 GeV scale characteristic of visible matter [24][25][26][27]. Current-quark masses enter into quantum chromodynamics (QCD) via the Higgs mechanism of explicit chiral symmetry breaking; so studying the spectra of states built from {u, d, s, c, b} quarks and/or related antiquarks opens the way to exploring constructive interference between Nature's two distinct mass generating mechanisms.
Quark models have also been widely employed in the calculation of hadron spectra; see, e.g. Refs. [28][29][30][31][32]. Owing to complications introduced by the need to preserve chiral symmetry and its pattern of breaking in QCD, i.e. to faithfully express DCSB, quark models are most reliable for baryons and mesons containing heavy quarks. In such systems, Higgsinduced current-quark masses play the dominant role and differences between current-and constituent-quarks are least noticeable; see, e.g. Ref. [27,Fig. 2.5].
During the past decade, continuum Schwinger function methods, particularly Dyson-Schwinger equations, have increasingly been used in spectrum calculations following improvements in both [5,[33][34][35][36][37]: (i) understanding the capacities and limitations of the approach; and (ii) the quality and range of the description of hadron properties. Recently, predictions for meson and baryon spectra in some of the low-lying flavour-SU(N f = 5) multiplets were delivered [38,39]. They were made using the rainbow-ladder (RL) truncation of the continuum bound-state equations. This is the leading-order in a systematic scheme [40][41][42][43][44][45][46][47], which enables both: (i) a unified, symmetry-preserving description of mesons and baryons, accounting properly for DCSB; and (ii) a traceable connection to quantum chromodynamics (QCD). A next challenge here is to advance beyond the leading-order truncation. That effort will likely benefit from the use of high-performance computing.
Complementing such studies, Refs. [48,49] followed a different, largely algebraic path to calculating the spectra of ground-state pseudoscalar and vector ( fḡ) mesons and J P = 1/2 + , 3/2 + ( f gh) baryons, where f, g, h ∈ {u, d, s, c, b}. Exploiting the fact that a hadron's mass is a volumeintegrated (long-wavelength) quantity, so not very sensitive to details of the system's wave function, Refs. [48,49] used a symmetry-preserving treatment of a vector × vector contact interaction (CI) [50] to deliver insights into features of these systems that can be obscured in approaches that rely heavily on computer resources. For instance, the analysis found diquark correlations to be important in all baryons studied; and owing to the dynamical character of the diquarks, the lightest allowed diquark correlation is typically that which defines the dominant component of a baryon's Faddeev amplitude.
Herein, we extend Refs. [48,49] and explore the capacity of the contact interaction to also explain the spectra of scalar and pseudovector mesons and J P = 1/2 − , 3/2 − baryons. The meson spectrum is canvassed in Sect. 2, along with a discussion of leptonic decay constants, including those of scalar and pseudovector states. The masses and decay constants of π, K , η c , η b mesons are used to determine the current-masses of the {u = d, s, c, b} quarks and basic CI parameters. Predictions are then made for another sixty-eight quantities.
Connections with diquark masses and correlation amplitudes are drawn in Sect. 3. They are important because soft diquarks seemingly play an important role in hadron structure [5]; hence, serve usefully in developing a quark + diquark approximation to the baryon three-body problem [51][52][53][54].
Results for the spectra of ground-state J = 1/2 ± , 3/2 ± baryons are discussed in Sect. 4. Following Refs. [55][56][57], eight parameters are introduced to complete the CI definition in the baryon sector, with four serving to correct for omission of resonant (meson cloud) contributions in the quark + diquark Faddeev equation. Masses for 88 baryons subsequently emerge as predictions along with their CI spinflavour Faddeev amplitudes.
A summary and perspective is provided in Sect. 5. Appendices are included to sketch the CI, Appendix A; and detail the structure of flavour-SU(5) baryon Faddeev amplitudes, Appendix B.

Meson spectrum
Reference [48] calculated the masses and leptonic decay constants of ground-state pseudoscalar and vector ( fḡ) mesons, f, g ∈ {u, d, s, c, b}. Herein, we expand the coverage to include kindred scalar and pseudovector mesons, viz. their parity partners. All CI bound-states satisfy a Bethe-Salpeter equation of the form in Eq. (A.10). The pseudoscalar solution amplitude is given in Eq. (A.9) and the other systems are described by the following amplitudes: where Q μ γ ⊥ μ = 0 and I D is the 4×4 identity matrix in spinor space. When a momentum-dependent interaction is used, the amplitudes have a richer structure [58,59]; but the terms in Eqs. (1) are dominant [60][61][62].
It has long been known that RL truncation has some defects in the 0 + , 1 + channels because these systems possess significant orbital angular momentum and RL fails to generate sufficient spin-orbit repulsion [47,61,71,72]. Nonperturbatively generated corrections to RL truncation, such as those originating in EHM-induced dressed-quark anomalous chromomagnetic moments [73], can remedy this defect and generate the empirical splitting between parity partners [61] and radial excitations [47]. (Such corrections largely cancel in 0 − and 1 − channels.) An ameliorating expedient was introduced in Ref. [55], and refined in Refs. [56,57]. Namely, one includes a multiplicative factor in the Bethe-Salpeter kernel for each of these channels, g 0 + SO , g 1 + SO , with g 0 + SO chosen to produce a mass difference of approximately 0.3 GeV between the quark-core of the 0 + ud and that of the ρ-meson (as obtained with beyond-RL kernels [61]) and g 1 + SO fixed to obtain the empirical size of the a 1 -ρ mass-splitting: To assist with understanding the comparisons in Table 1, we also represent them in Fig. 1. Considering Fig. 1a, one sees that the CI delivers good estimates for the masses of both positive and negative parity ground-state flavour-SU(5) mesons. Aspects of its symmetry-preserving formulation skew results toward overestimates in most cases [74], but the mean relative-difference between theory and experiment/lQCD is only 5(6)%. (This value is obtained by comparing columns 1 and 2 in Table 1, omitting the underlined entries in column 1.) Furthermore, as found with pseudoscalar and vector mesons in Ref. [48], the computed masses of scalar and pseudovector mesons neatly follow a pattern prescribed by equal spacing rules (ESRs) [38,39,[75][76][77], e.g.
These relations expose the scales that may be identified with the splittings between the spectrum-generating constituentlike quark masses, e.g. Eqs. (4d) reveal a b − c spectrum mass difference d bc M = 3.18(4) GeV. This is a fair match with M b − M c = 3.23 GeV from Table 7. As evident in this table, the Higgs mechanism of mass generation is the origin of such splittings.
Interesting patterns are also found by comparing parity partners, e.g.
These results reveal the scale of DCSB in meson spectra.
Overall, parity partners are split by 0.35(10) GeV. However, comparing the impact of Higgs couplings into QCD, the expression of DCSB becomes weaker with increasing current-quark mass: Eq. (5a) cf. (5b) cf. (5d); and Eq. (5c) cf. (5e). Moreover, DCSB generates larger scalar-pseudoscalar splittings than pseudovector-vector: Eq. (5e) cf. (5d); and Eq. (5c) cf. (5b). CI results for meson leptonic decay constants are also listed in Table 1 and drawn in Fig. 1b. These quantities describe quark + antiquark annihilation at a single spacetime point, so they are sensitive to ultraviolet physics. In QCD, this is expressed through the appearance of a logarithmic ultraviolet divergence in integrals like those in Eqs. (1), which is compensated by the dressed-quark wave-function renormalisation constant [78,79]. It is therefore unsurprising that the cutoff-regularised CI supplies a poorer description of meson decay constants than it does of their masses. It would be worse if Eqs. (A.18) and (A.21) were not implemented. Nevertheless, in comparison with known empirical or lQCD values, the picture is fair: trends are generally reproduced and the mean-absolute-relative-difference between the entries in columns 5 and 6 of Table 1 is 15(11)%. (Given that very little is known about the decay constants of scalar and pseudovector mesons, they are excluded from this comparison.) Peculiarities of the CI's formulation mean that the description is better for pseudoscalars than it is for vector mesons [74]. Further improvement can be achieved by additional tuning of Eq. (A.18), including its extension into the light-quark sector and consideration of the particular properties of heavy + light mesons. Table 1 Mesons: computed masses, Bethe-Salpeter amplitudes, and decay constants. Empirical masses from Ref. [1]; entries marked by " * ", from Ref. [21]. Empirically unknown decay constants quoted from lQCD [20,[64][65][66][67][68][69]. (Dimensioned quantities in GeV. Underlined entries from Table 7. "-" indicates no empirical/lQCD results available for comparison)  [21]. Lower panelb CI predictions for meson leptonic decay constants compared with experiment [1], where known, and lQCD otherwise [64][65][66][67][68][69]. In both panels, contact-interaction predictions are depicted as (blue) circles and comparison values by (green) bars. (Pictorial representation of results in Table 1) Notwithstanding these observations, as with pseudoscalar and vector mesons, analogues of Eqs. (4) are applicable, e.g.
Such relations are not surprising when one recalls that such decay constants are order parameters for chiral symmetry breaking [80]. Regarding Fig. 1b, one finds that using the CI the decay constants of vector mesons and those of pseudovector mesons have the opposite tendency with increasing current-quark mass, e.g.
It is currently unclear whether this is a CI artefact. However, it does follow the trend set by the meson Bethe-Salpeter amplitudes in Table 1 and match qualitatively with the behaviour of scalar-meson decay constants. The question could be answered by using a momentum-dependent quark+antiquark scattering kernel [62].

Diquark mass spectrum
In approaching the problem of calculating the spectrum of flavour-SU(5) ground-state baryons with J P = 1/2 ± , 3/2 ± , we use a quark+diquark approximation to the Faddeev equation introduced elsewhere [51 -54]. As apparent in Fig. 2, here the diquarks are fully dynamical, appearing in a Faddeev kernel that requires their continual breakup and reformation. In consequence, they are vastly different from the static, pointlike degrees-of-freedom considered in early models of baryons [81]. To proceed, therefore, one must calculate the masses and amplitudes for all diquark correlations that can exist in J P = 1/2 ± , 3/2 ± baryons. Using the CI, this means the following colour-3 correlations: flavour-10 scalar; flavour-15 pseudovector; flavour-10 pseudoscalar; and flavour-10 vector. 1 This is a straightforward exercise because the RL Bethe-Salpeter equation for a J P diquark is generated from that for a J −P meson by simply multiplying the meson kernel by a factor of 1/2 [83]. For instance, working from Eq. (A.10), the mass and amplitude for a flavour-10 scalar diquark is obtained from the following equation: where the correlation amplitude is Γ [ f g] (Q) and with C = γ 2 γ 4 being the charge-conjugation matrix. The canonical normalisation conditions are likewise amended, with the multiplicative factor being 2/3 in this case. In a study of flavour-SU(3) baryon spectra [57], an additional step was found to be necessary in order to arrive at realistic results. Namely, considering that valence-quarks within a diquark are less tightly correlated than the valence-quark and -antiquark in a bound-state meson, then spin-orbit repulsion in diquarks should be weaker than it is in mesons. This was effected by writing so that the modification factor in the RL-like diquark Bethe-Salpeter equations is nearer unity; hence, generates less repulsion. The value of s SO was chosen in concert with g π Ψ π d D B , Eq. (13), in order to ensure that the dressed-quark core mass of the nucleon's parity partner is not lower than that of its first radial excitation.
Following this approach, using parameters determined as described here, in Sect. 2 and Appendix A, one obtains the diquark masses and amplitudes listed in Table 2. Consistent with RL studies using realistic interactions [84,85], in the Bethe-Salpeter equation for a given J P diquark correlation we use the values of α IR , Λ uv associated with its J −P meson partner. A ±10% change in s SO alters the 0 − , 1 − values by ∓2%. Evidently, the antisymmetric Dirac-scalar combination of any two quark flavours (scalar diquark) is always lighter than the symmetric J P = 1 + combination (pseudovector diquark, denoted { f g} 1 + ), followed by the pseudoscalar and then the vector diquarks. Just like mesons, the pattern of masses can be understood in terms of ESRs.
It is worth reiterating that diquark correlations are coloured and it is only in connection with partnering coloured objects (quark, another diquark, etc.) that a colour singlet system is obtained. Hence, diquarks are confined. That is not true in RL truncation [84]; but corrections to this leading-order analysis have been studied using an infrared-dominant interaction [86]. In fully self-consistent symmetry-preserving treatments, such corrections eliminate diquark bound-state poles from the quark + quark scattering matrix, whilst preserving the strong correlations [87]. In these cases, one still has a physical interpretation of the RL masses; to wit, the quantity ( f g) J P := 1/m ( f g) J P is a spacetime distance that the diquark correlation can propagate before fragmentation. Figure 3 compares calculated diquark masses with those of their partner mesons. The level ordering of diquark correlations is precisely the same as that for mesons. Scalar and pseudovector diquarks are heavier than their partner mesons. Owing to the introduction of s SO in Eq. (10), this ordering is reversed for pseudoscalar and vector diquark correlations. Notwithstanding this, omitting the [ud] 0 + and [us] 0 + correlations, the mass of a diquark's partner meson is a reasonable guide to the diquark's mass: the mean difference in absolute value is 0.08 (7) GeV.
Owing to DCSB, the light-quark scalar diquark channels are atypical because their partner mesons are the almost-Nambu-Goldstone modes, π and K . In two-colour chiral-limit QCD, scalar diquarks are also Nambu-Goldstone modes [88]. A hint of this remains in the three-colour theory, expressed in the relatively low masses of [ud] 0 + and [us] 0 + , even though they are split widely from the π and K .
It is the canonically normalised diquark Bethe-Salpeter amplitudes that appear in baryon Faddeev equations. Listed in Table 2, they exhibit a clear ordering: in any given flavour sector, scalar diquarks have the largest amplitudes ↔ strongest couplings. This means that the lighter flavour-10 scalar diquark correlations are typically favoured in all  J = 1/2 baryon amplitudes because the Faddeev equations involve the diquark-amplitude-squared. As will become apparent, under certain circumstances, e.g. in baryons whose valence-quarks have widely different masses, this preference can be overcome. Notably, the flavour-exchange symmetries of J P = 3/2 ± baryons mean that such CI states may only contain flavour-15 pseudovector diquark correlations, viz. positive-parity diquarks. Using momentum-dependent interactions, flavour-15 vector diquarks are possible. However, these correlations are weak; hence, even then, are unlikely to play a significant role. Nevertheless, this expectation should be checked.

Baryon spectrum
Following Ref. [48], we use the Faddeev equation drawn in Fig. 2 to calculate the spectrum of ground-state flavour-SU(5) J P = 1/2 ± , 3/2 ± baryons. In the isospin-symmetry limit, 88 distinct states are supported. Details are presented in Appendix B.

Completing the Faddeev kernels
As highlighted above, the kernel in Fig. 2 introduces binding through diquark breakup and reformation via exchange of a dressed-quark. Again following Ref. [48], we exploit an oft used simplification, viz. in the Faddeev equation for a baryon of type B, the quark exchanged between the diquarks, Eq. (B.1), is represented as where g = l, s, c, b is the quark's flavour and g B is a coupling constant. This is a variation on the "static approximation" introduced in Ref. [89]. It makes the Faddeev amplitudes momentum-independent, just like the diquark Bethe-Salpeter amplitudes. Calculations reveal that it has little impact on the calculated masses [90]. The couplings g N , g Δ are treated as parameters, with values chosen to obtain desired masses for the nucleon and Δ-baryon, m N , m Δ . The Faddeev equation, Fig. 2, generates what may be called a baryon's dressed-quark core [91]: it omits what are typically described as meson-cloud contributions to the baryon masses, which work to reduce m N , m Δ . The sizes of such corrections has been estimated: for the nucleon, the reduction is roughly 0.2 GeV and for the Δ it is 0.16 GeV. The choices [55] g N = 1.18, g Δ = 1.56, produce m N = 1.14 GeV and m Δ = 1.39 GeV, i.e. inflated masses that leave room for correction by meson cloud effects.
(N.B. The values in Eq. (12) are deliberately different from those used in Ref. [48].) Using the framework outlined above, Ref. [57] studied flavour-SU(3) baryons, finding that the ground-state positive-parity octet baryons are primarily constituted from like-parity diquarks, with negligible contributions from negative-parity correlations. Somewhat surprisingly, it found that their parity partners are also dominated by positive-parity diquark correlations; hence, too light. As with mesons and diquarks, the missing element was identified as too little spinorbit repulsion generated by RL-like kernels. This was remedied by inserting an additional parameter into the Faddeev equation for J P = (1/2) P baryons, i.e. g π Ψ π d D B , a linear multiplicative factor, expressing EHM-induced enhancement of spin-orbit repulsion, attached to each diquark amplitude in the baryon's Faddeev equation kernel, Eq. (B.1): i.e. nothing is done when the parity of the diquark correlation, π d , matches that of the host baryon, π Ψ , but suppression is introduced when the parity is different. The magnitude of the effect was chosen so that, in concert with Eqs. (10), the dressed-quark core mass of the nucleon's parity partner is not lower than that of its first radial excitation [57]. Every element in the Faddeev kernel is now specified. It remains only to select the channel, derive the explicit form for the associated algebraic equation, then solve the resulting matrix equation to obtain the baryon mass and amplitude. Specific examples of the Faddeev equations for the various light-quark systems, along with their derivations, are presented in Ref. [70]; and equations for Ξ + c , Ξ + c , Ξ * + c are given in Ref. [48]. The equations solved herein are not any more complex, although they do involve larger kernel matrices because we also allow all baryons to contain pseudoscalar and vector diquarks. As in Ref. [48], the ultraviolet cutoff in each baryon channel is identified with that of the lightest diquark in the system. This is always the smallest value; hence, the dominant regularising influence.

J P = 1/2 +
Our results for the masses and amplitudes of flavour-SU(5) J P = 1/2 + baryons are listed in Table 3. As remarked in association with Eq. (12), the Faddeev kernel employed omits resonant contributions that serve to reduce baryon masses. From each of the values in column 1 we have therefore subtracted in order to express an empirically informed estimate of such corrections. This value is the average difference between the unmodified CI result and experiment in the first four rows.     Table 3 compared with either experiment (first 15) [1] or lQCD (last 9) [18,21]. Lower panelb Sim-ilar comparison for ground-state flavour-SU(5) J P = 3/2 + baryons in Table 4: experiment (first 9) [1]; or lQCD (last 11) [18,21]. Both panels: CI results -(blue) circles; reference values -(green) bars It matches the expectations described in connection with Eq. (12). As indicated after Eq. (12), Ref. [48] achieved a similar effect by inflating the values of g N ,Δ , but we find that can produce unexpected (perhaps unphysical) modifications of the Faddeev amplitudes.
The computed masses in Table 3 are also compared with empirical or lQCD values in Fig. 4a: the mean absoluterelative-difference is 1.
Compared with the results in column 4, drawn from Ref. [39], the analogous difference is 5.2(2.8)%. In that study, the Faddeev equations were solved in a fully-consistent RL truncation, eschewing a quark + diquark approximation. Moreover, the results are ab initio predictions, whereas we introduced g N , z N + -Eqs. (12) and (14) -to adjust the overall scale of the J P = 1/2 + spectrum. Notwithstanding these things, the level of agreement confirms the validity of the ESRs used to complete the spectrum calculations in Ref. [39]; and the fact that our results exceed the accuracy of Ref. [39] indicates that we have implemented a phenomenologically efficacious CI formulation.
The 1/2 + Faddeev amplitudes are represented in Table 3 Table 3 confirm the conclusions in Ref. [48]. Namely, concerning ground-state flavour-SU (5) J P = 1/2 + baryons: (a) the lightest participating diquark correlation usually defines the most important component of a baryon's Faddeev amplitude and this is true even if a pseudovector diquark is the lightest channel available; and (b) light-diquark dominance may be overcome in flavour channels for which the bound-state's spin-flavour structure and the quark-exchange character of the kernel in Fig. 2 lead dynamically to a preference for mixed-flavour correlations.
In corroborating these conclusions, we reaffirm the arguments against treatments of the baryon problem which assume they can be described as effectively two-body in nature, e.g. as being built from a constituent-quark and static constituent-diquark. The dynamical character of diquark correlations is paramount because their breakup and reformation are crucial in defining baryon structure. This is confirmed for light-quark baryons by the fact that lQCD calculations produce a spectrum whose richness cannot be explained by a two-body model [17]. The implications extend to baryons involving one or more heavy quarks, challenging both (i) the treatment of singly-heavy baryons as two-body lightdiquark + heavy-quark bound-states and (ii) analyses of doubly-heavy baryons which assume such systems can be considered as two-body light-quark + heavy-diquark boundstates. These observations may also be relevant to few-body studies of tetra-and penta-quark problems; in particular, those involving {u, d, s} quarks.

J P = 3/2 +
Our results for the masses and amplitudes of flavour-SU(5) J P = 3/2 + baryons are listed in Table 4. Akin to Eq. (14), from each of the values in column 1 we have subtracted so as to express an empirically informed estimate of mesoncloud corrections: this value is the average difference between the unmodified CI result and experiment in the first four rows. Again, it matches expectations described in connection with Eq. (12).
The computed J P = 3/2 + masses are compared with empirical/lQCD values in Fig. 4b: the mean-absoluterelative-difference is 1.0(0.8)%. Furthermore, the results compare well with the three-body calculation described in Ref. [39], for which the analogous difference is 2.6(1.6)%. Notwithstanding the fact that we used g Δ , z Δ + -Eqs. (12), (15) -to adjust the overall scale of the J P = 3/2 + spectrum, this improvement over the results in Ref. [39] again highlights the utility of our CI formulation.
The CI Faddeev amplitudes of J = 3/2 + baryons are also listed in Table 4. Once more, the fully dynamical nature of the diquarks and the character of the Faddeev kernel work together to ensure a continual reshuffling of the dressedquarks within the diquark correlations. Consequently, in all cases involving more than one quark flavour, the diquark combination with maximal flavour shuffling is favoured because it is fed by twice as many exchange processes as the less-mixed correlation. This differs from the result for the Σ * b -baryon in Ref. [48]. Therein, however, g Δ was 26% larger than here, thereby providing additional enhancement for the single u-quark exchange process over the doubly active, but b−u mass-splitting-suppressed, b-quark exchange contribution which enables the shuffling. Table 5 Computed mass (in GeV) and Faddeev amplitude for each ground-state flavour-SU  Table 5 compared with experiment [1] (green bars), lQCD [22] (gold triangles), or three-body Faddeev equation results [39] (orange asterisks). Lower panelb Similar com-parison for ground-state flavour-SU(5) J P = 3/2 − baryons, compiled using Table 6 Table 5. From each of the values in column 1, we subtracted z N − = 0.27 GeV (16) in order to express an empirically informed estimate of meson cloud corrections: this value is the average difference between the unmodified CI result and experiment in the first three rows. The computed masses in Table 5 are also depicted in Fig. 5a. Compared with empirical values, or lQCD results when empirical values are unavailable, the mean absoluterelative-difference is 1.4(1.1)%; and compared with the results in column 4, drawn from Ref. [39], the analogous difference is 3.7(3.5)%.
The CI Faddeev amplitudes of J P = 1/2 − baryons are also listed in Table 5. Their characters are readily understood because they typically follow the pattern of the J P = 1/2 + amplitude: the lightest like-parity diquark correlation usually dominates. The exceptions to this simple rule are chiefly driven by the absence of flavour-symmetric vector diquark correlations, e.g. whereas Σ c is predominantly c{uu} 1 + , Σ P=− c is dominated by u[uc] 1 − because this is the lightest J = 1 correlation available and J = 1 is favoured in Σ c .

J P = 3/2 −
These states are somewhat unusual. Given that they must be completely symmetric under permutations of flavour labels but the only negative-parity diquarks supported by the CI are flavour-10 (antisymmetric) states, then the J = 3/2 − baryons can only contain positive parity pseudovector diquarks, Eq. (B.16). Owing to this peculiarity, we judged that opposite-parity diquark correlations should not receive additional suppression in these systems; so we solved the associated Faddeev equations with Table 6 Computed mass (in GeV) and Faddeev amplitude for each ground-state flavour-SU (5) J P = 3/2 − baryon: the last column identifies the baryon's dominant spin-flavour correlation. The possibilities are given in Eqs. (B.10), (B.12), (B.14), (B.16). Empirical mass values, M e , are taken from Ref. [1]; and lQCD results, M l , from Ref. [22]. The results under heading M 3 are the three-body Faddeev equation predictions in Ref. [ Working with the results thus obtained, the average difference between the CI result and experiment in the first three rows of Table 6 is Consequently, in column 1 of Table 6 we list the result obtained by subtracting z Δ − from the directly calculated CI result. Comparing column 1 with experiment, where known, and available lQCD results otherwise, the mean absoluterelative-difference is 2.2(1.4)%, and with the results in column 4, drawn from Ref. [39], the analogous difference is 1.4(1.9)%. The computed masses in Table 6 are also depicted in Fig. 5b alongside, as available, in this order: experiment, lQCD, or three-body Faddeev equation results.
The CI Faddeev amplitudes of J = 3/2 − baryons are also listed in Table 6. They are similar to those of the J = 3/2 + states in Table 4, with the dominant diquark correlation being the same in all cases except Σ * , Σ * c , Σ * b . Considering Eq. (B.12a), one sees that in these cases there is a competition between the lightest diquark correlation and the heavylight diquark, which is doubly-fed by the exchange kernel. The former wins in the negative-parity systems, whereas the latter is dominant in the positive-parity systems. This ordering reverses in the 3/2 − states as g is reduced from unity. The Σ * , Σ * c , Σ * b Faddeev amplitudes should therefore be considered as somewhat uncertain.

Summary and outlook
A confining, symmetry-preserving treatment of a vector× vector contact interaction (CI) was used to compute spectra of ground-state J P = 0 ± , 1 ± ( fḡ) mesons and J P = 1/2 ± , 3/2 ± ( f gh) baryons, where f, g, h ∈ {u, d, s, c, b}. The calculated meson masses agree well with experiment (Sect. 2): the mean-relative-difference for 33 states is 5(6)%. Expressing effects tied to the emergence of hadronic mass (EHM) was crucial to achieving this level of agreement. Regarding meson leptonic decay constants, empirical (or lattice QCD -lQCD) values are reproduced with an accuracy of 15(11)%; and predictions were made for the currently unknown decay constants of 19 0 + , 1 + states (Table 1).
A quark+diquark approximation to the baryon Faddeev equation was used herein. Its formulation required the calculation of masses and correlation strengths for all 38 distinct participating diquarks (Sect. 3). As in all studies to date, the level ordering of the J P diquarks matches that of their J −P meson partners. Scalar and pseudovector diquarks are heavier than their partner mesons; but in our CI formulation, this ordering is reversed for pseudoscalar and vector diquark correlations. Nevertheless, omitting the diquarks partnered with would-be Nambu-Goldstone mode mesons, the mass of a diquark's partner meson is a reasonable guide to the diquark's mass: the mean difference in absolute value is 0.08(7) GeV.
A static approximation to the quark exchange kernel was used to solve the Faddeev equation. It produces momentumindependent Faddeev amplitudes, thereby ensuring a level of consistency with CI two-body bound-state amplitudes, and introduces four parameters. Correcting also for the omission of meson cloud contributions to the quark-core Faddeev kernels, four additional parameters were necessary to complete the CI definition in the baryon sector. In total, the CI predicts 88 distinct baryons, i.e. every possible three-quark 1/2 ± , 3/2 ± ground-state baryon is realised. Of this number, 34 states are already known empirically and lQCD-computed masses are available for another 30. For this collection of 64 states, the mean absolute-relative-difference between CI prediction and experiment/lQCD mass is ard = 1.4(1.2)% (Figs. 4, 5). Implementation of EHM-induced effects associated with spin-orbit repulsion in 1/2 − baryons was important to achieving this outcome. The same 88 ground-states are also predicted by a three-body Faddeev equation [39]; and in comparison with those results, the analogous difference is 3.4(3.0)%.
A primary merit of the framework employed herein is its simplicity, enabling all analyses and calculations to be completed algebraically. In total, there are twelve parameters: four used to define the interaction and its scale dependence via π , K , η c , η b properties; and eight introduced to complete the baryon Faddeev equations. From this foundation, the CI delivers predictions for 164 distinct quantities. Thus far, 114 of these observables have either been measured or computed using lQCD; and a comparison on this subset yields ard = 4.5(7.1)%. This level of quantitative success suggests that some credibility be given to the qualitative conclusions that follow from our CI analysis. (I) Nonpointlike, dynamical diquark correlations play an important role in all baryons; and, typically, the lightest allowed diquark is the most important component of a baryon's Faddeev amplitude. (II) Positive-parity diquarks dominate in positive-parity baryons, with J = 1 + diquarks being prominent in all of them. (III) Negativeparity diquarks can be neglected when studying positiveparity baryons; but owing to EHM, they are significant, even dominant, in J = 1/2 − baryons. (IV) On the other hand, J = 3/2 − baryons are built (almost) exclusively from J = 1 + diquark correlations. Naturally, these conclusions should be checked using more sophisticated Faddeev equa-tions with momentum-dependent exchange interactions, e.g. extending Refs. [77,82].
An extension of this analysis to radial excitations of the states considered herein is possible. However, that requires additional intervention because one must, by-hand, force zeros into what would otherwise be momentum-independent Faddeev amplitudes [55]. A potentially more worthwhile direction would be to adapt the framework to the challenge of understanding tetra-and penta-quark states [5,6].
A key longer-term goal is extension of the rainbow-ladder truncation three-body Faddeev equation study in Ref. [39] so that it can: (a) directly handle non-degenerate valence-quark systems (reducing/eliminating the need to use equal spacing rules); and (b) include beyond-rainbow-ladder contributions and thus both search for signals indicating the appearance of diquark correlations and deliver improved ab initio predictions. In such an effort, the use of high-performance computing resources will be necessary. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Contact interaction
We employ the symmetry preserving treatment of the vector × vector contact interaction (CI) described in Ref. [50].
The key element is the quark+antiquark scattering kernel, which in rainbow-ladder (RL) truncation can be written (k = p 1 − p 1 = p 2 − p 2 ): where k 2 T μν (k) = k 2 δ μν − k μ k ν . (Our Euclidean metric and Dirac-matrix conventions are given in Ref. [70,Appendix A].) The defining quantity isG. Following two decades of study, much has been learnt about its pointwise behaviour. The qualitative conclusion is that owing to the emergence of a gluon mass-scale in QCD [92,93],G saturates at infrared momenta. Hence, one may writẽ In QCD, m G ≈ 0.5 GeV, α IR ≈ π [94]. Following Ref. [48], we retain this value of m G , but reduce α I R to a parameter. This is needed because the integrals appearing in CI bound-state equations require ultraviolet regularisation and this undermines the link between infrared and ultraviolet scales that is distinctive of QCD. Furthermore, since a CI cannot support relative momentum between bound-state constituents, it is sensible to simplify the tensor structure in Eqs. (A.1) such that, in operation: As remarked already, any use of Eq. (A.3) in an equation related to the continuum bound-state problem will require ultraviolet regularisation of any integrals involved. Additionally, the theory is not renormalisable, so the associated massscales, Λ uv , are additional physical parameters. They may usefully be interpreted as upper bounds on the momentum domains within which the properties of the associated systems are effectively momentum-independent. For instance, the π -meson is larger in size than the η b ; hence, one should expect to use 1/Λ π uv > 1/Λ η b uv . For the purpose of expressing confinement, we also introduce an infrared regularisation scale, Λ ir , when defining bound-state equations [95]. This has the effect of excising momenta less than Λ ir , thereby eliminating quark+antiquark production thresholds [96]. A typical choice for this scale is Λ ir = 0.24 GeV [50].
In implementing Eq. (A.3), Ref. [48] fixed the parameters via the masses and leptonic decay constants of the π , K , η c , η b , i.e. 0 −+ ff mesons, f ∈ {l = u = d, s, c, b}, where isospin symmetry is assumed. The starting point is the dressed-quark gap equation, which in RL truncation, using Eq. (A.3), takes the following form: where m f is the quark's current-mass. Regularising the integral in a Poincaré-invariant manner, the solution is with M f obtained as the solution of: Γ (α, y) is the incomplete gamma-function.
The following functions arise in solving the bound-state equations considered herein (τ 2 uv = 1/Λ 2 uv , τ 2 ir = 1/Λ 2 ir ): iu n (σ ), n = 0, 1, 2, . . . . The contact interaction Bethe-Salpeter amplitude for a 0 −+ fḡ meson has the following form [50,97]: Q is the bound-state's total momentum, Q 2 = −m 2 0 − , m 0 − is the meson's mass; and The amplitude is obtained by solving (t + = t + Q): (A.10) Using the symmetry-preserving regularisation scheme introduced in Refs. [50,97], which requires 12) one arrives at the following Bethe-Salpeter equation: is an eigenvalue problem. It has a solution for Q 2 = −m 2 0 − , at which point the eigenvector is the meson's Bethe-Salpeter amplitude. In the calculation of observables, one must use the canonically normalised amplitude, viz. the amplitude rescaled such that (A. 16) The pseudoscalar meson's leptonic decay constant is given by: Following Ref. [48], we use the light-quark results from Refs. [55,97], listed in Table 7. The fitted value of m s /m l = 24 is compatible with estimates in QCD [1], even though the individual current-masses are too large by a factor of 2 because of the CI's deficiencies in connection with ultraviolet quantities. The result M s /M l = 1.4 is commensurate with the value obtained in efficacious RL studies with momentumdependent interactions [24]: M s /M l = 1.25 (9).
In heavy-quark systems, Λ 0 − uv is allowed to vary with the meson's mass and the associated coupling is fixed by requiring (A.18) This procedure serves to limit the number of parameters, so that in fitting the η c,b quantities in Table 7 there are only two parameters for each case: m c,b , Λ η c,b uv . Considering the η b , a lQCD calculation reports f η b = 0.472(4) [67], but this is larger than the result for f Υ = 0.459 (22) [69]; hence, it is contrary to the experimental pattern: f π < f ρ , f η c < f J/ψ . Therefore, Ref. [48] chose to constrain m b , Λ η b uv via known experimental results [1]: The The physical origin of these outcomes is clear: the f 0 − integral diverges logarithmically with increasing Λ uv ; and the flow of α IR compensates for analogous behaviour in the Bethe-Salpeter kernel, thereby maintaining the given meson's mass.

B.2 Flavour structure of the amplitudes
Working in the isospin-symmetry limit, the Faddeev equation in Fig. 2 supports 88 distinct baryon states, which we list below. The flavour structure of the amplitudes reflects the quark+diquark approximation.