# A Multiquark Approach to Excited Hadrons and Regge Trajectories

Authors: S. S. Afonin
Journal: Advances in High Energy Physics
 (2019)

---

## Abstract

We propose a novel approach to construction of hadron spectroscopy. The case of light nonstrange mesons is considered. By assumption, all such mesons above 1 GeV appear due to creation of constituent quark-antiquark pairs insideπ or ρ(ω) mesons. These spin-singlet or triplet pairs dictate the quantum numbers of formed resonance. The resulting classification of light mesons turns out to be in a better agreement with the experimental observations than the standard quark model classification. It is argued that the total energy of quark components should be proportional to the hadron mass squared rather than the linear mass. As a byproduct a certain relation expressing the constituent quark mass via the gluon and quark condensate is put forward. We show that our approach leads to an effective mass counting scheme for meson spectrum and results in the linear Regge and radial Regge trajectories by construction. An experimental observation of these trajectories might thus serve as evidence not for string but for multiquark structure of highly excited hadrons.


---

## Body

## 1. Introduction

The regularities in the masses of atomic nuclei tell us that nuclei certainly consists of some building blocks, the nucleons, and the mass of nuclei is just proportional to the number of nucleons. The masses squared of light mesons also reveal many regularities [1, 2] but they are usually attributed to the properties of strong interactions acting between valent quarks. In the present paper, we will try to develop an alternative point of view on the masses of light hadrons: within the presented scheme, they could be also viewed as originating from some effective “building blocks”.Historically the main basic framework for description of hadrons is the quark model. This model had a great success in making order in numerous hadronic zoo. In particular, all mesons within this picture are composed of a quark and antiquark having the relative angular momentumL which dictates the spatial and charge parities for the quark-antiquark systems, P=(-1)L+1 and C=(-1)L+S, where S is the total quark spin. By now we know that the ensuing classification of mesons has many phenomenological flaws (see, e.g., the corresponding reviews in the Particle Data [3]). For instance, it does not explain the scalar mesons below 1 GeV for which a tetraquark structure is often assumed [4]. A question arises why do not we see many light multiquark hadrons above 1 GeV? Among other questions one can mention the experimental absence of many predicted states (e.g., with JPC=2--) and, on the other hand, observation of many not predicted states (e.g., too rich scalar sector and exotic π1).There are some deep theoretical questions as well. The light hadrons represent highly relativistic quantum systems in whichL is not a conserved quantity. The quark angular momentum L is nevertheless a standard ingredient in constructing dynamical models for light hadrons. Another question concerns the observable quantities like hadron masses: they must be renorminvariant in the field-theoretical sense. Stated simply, they must represent some constants independent of energy scale. The quark masses, for instance, are not such constants since they have anomalous dimension described by QCD. A relation of calculated observables to the renormalization invariance is obscure in many known phenomenological models of hadrons.The quark model represents a concept rather than a model. A real model appears only when a definite interaction between quarks is postulated. Any such model should be able to explain general features of observed hadron spectrum. In the light hadrons, perhaps the most spectacular general feature is the observation of approximately linear Regge and radial trajectories of the kind:(1)mJ,n2=aJ+a′n+b,J,n=0,1,2,….Here J denotes the spin, n is the radial quantum number (enumerating the daughter Regge trajectories), and a, a′, and b are the slope and intercept parameters. There is no general agreement on the position of different states on the trajectories (1) (see, e.g., [1, 2, 5–20]) but the existence of these trajectories seems to be certain for many phenomenologists [12–27]. There are also some indications on a similar behavior in the heavy quark sector [28–32]. The nonrelativistic and partially relativistic potential models cannot explain the recurrences (1) in a natural way. Usually the observation of linear trajectories (1) is interpreted as an evidence for string picture of mesons. In spite of many attempts (see, e.g., [33–40]), however, a satisfactory quantized hadron string has not been constructed. Among typical flaws of this approach one can mention the absence of spontaneous chiral symmetry breaking, rapidly growing (with mass) size of meson excitations that is not supported experimentally, unclear role of higher Fock components in the hadron wavefunction.The purpose of the present work is to propose a novel realization of quark model concept, a realization that leads to a natural (and alternative to hadron strings) explanation of Regge recurrences (1) and that is potentially free of typical shortcomings inherent to string, potential, and some other approaches.


## 2. Masses and Classification of Light Mesons

In hadron physics, the pion is known to be the most important and best studied meson both experimentally and theoretically. Concerning the theoretical aspect, the pion is the only hadron (along withK and η) for which we know a model-independent mass formula: the Gell-Mann–Oakes–Renner (GOR) relation [41],(2)mπ2=-q¯qfπ2mu+md=Λ·2mq,where we set mu=md≐mq and Λ≐-q¯q/fπ2. We will use the standard values for the quark condensate and masses of current quarks at the scale of the pion mass from the QCD sum rules [42, 43] and Chiral Perturbation Theory (ChPT) [44], q¯q=-(250MeV)3, mu+md=11 MeV. Using fπ=92.4 MeV (the value of weak pion decay constant), relation (2) yields mπ=140 MeV and Λ=1830 MeV. The famous relation (2) was derived in various approaches assuming the spontaneous chiral symmetry breaking (CSB) in strong interactions.The renormalization invariance of pion mass follows from the renorminvariance of operatormqq¯q in QCD Lagrangian. Masses of other mesons must be also renorminvariant. On the other hand, the light nonstrange mesons can be viewed, in one way or another, as some quantum excitations of pion. We will assume that the masses of these mesons can be represented as(3)mh2=ΛEh+2mq=ΛEh+mπ2.Below we will motivate a proportionality of the energy parameter Eh to the renorminvariant gluon condensate, Eh~αsGμν2/q¯q. But first we would like to develop a phenomenological mass counting scheme based on the relation (3) and demonstrate how it describes the meson spectroscopy.Our second assumption is that the parameterEh can be interpreted as an effective energy of some constituents different from the current quarks. The problem is to propose a model for these constituents which should represent some excitations inside the pion. We postulate three basic excitations. The first one appears when one of quarks absorbs a gluon of certain energy Eh=Eρ. The spin of excited quark changes its direction to the opposite one, converting the original spin-singlet qq¯-pair (π-meson) to the spin-triplet one (ρ-meson). The second kind of excitation emerges due to formation of spin-singlet qq¯-pair with effective mass Eh=E0 (the lower index stays for total spin of S-wave qq¯-pair). We will call this pair Q0. The third basic excitation is the formation of spin-triplet qq¯-pair with effective mass Eh=E1. This pair will be called Q1. The formation of constituent pairs Q0 and Q1 may be a QCD analogue of formation of para- and orthopositronia from photons. We assume that any excitation inside pion leading to an observable resonance can be represented as a combination of these basic excitations so that the total effective energy Eh in (3) is just a sum (appropriate number of times) of Eρ, E0, and E1. Simultaneously this will dictate the quantum numbers of constructed resonance (Q0 and Q1 possess JPC=0-+ and 1--, respectively). In a sense, the introduction of Q0 and Q1 is our model for excitations of higher Fock components in the pion wavefunction.It is convenient to divide the nonstrange meson resonances above 1 GeV into excited pions and excitedρ(ω) mesons. The excited ρ mesons containing n pairs Q0 and l pairs Q1 will be labeled as ρJQ0nQ1l, where J means the total spin. According to our mass counting scheme, masses of these excitations are given by the relation(4)mρQ0nQ1l2=ΛEρ+nE0+lE1+mπ2.By assumption, the whole system is in the S-wave state (one can include into this scheme the higher L>0 waves; we will consider the simplest possibility); hence, at fixed n and l, relation (4) defines mass of a set of degenerate states with various spins up to J=l+1. P and C parities obey the standard multiplicative law: (P,C)=((-1)n+l+1,(-1)l+1). The total isospin of system of n+l constituent qq¯ pairs is zero as this system, by assumption, arises from an isosinglet combination of gluons.The excited pions follow the same principle. The statesπJQ0nQ1l have mass(5)mπQ0nQ1l2=ΛnE0+lE1+mπ2,maximal spin J=l, and parities (P,C)=((-1)n+l+1,(-1)l).Since a simultaneous excitation of two close pairsQ0Q0 does not change neither spin nor parities, it is natural to refer to a meson containing even number 2k of Q0 and no one Q1 as the k-th radial excitation. The radial excitations appear thus with a “period” 2E0Λ.Consider a formation of two close vector pairsQ1Q1. As they are formed from originally massless gluons we must combine their spins as if they were massless particles. This follows from conservation of total spin of original gluons. We recall that the spin of massless gluons has only two projections: two possible helicities. The total helicity of two gluons can be 0 or 2. This rule (in other words, the addition law for spins of spin-1 massless particles is the same as for massive spin -1/2 particles; just the final result must be multiplied by 2) extends to arbitrary number of close Q1 pairs. It is easy to see now that the excitations of the kind ρJQ12k, k=0,1,2,…, will give rise to a degenerate family of resonances with JPC=(1,3,…,2k+1)--. The states on the main ρ-meson Regge trajectory follow thus with a “period” 2E1Λ. These excitations generate also daughter trajectories. For instance, the additional ρ mesons (JPC=1--) arise which appear with the same “period” 2E1Λ. They are different from the “radial” ρ mesons that appeared with another period 2E0Λ.In order to demonstrate how our scheme works in practice let us consider some examples.The mass ofρ-meson is mρ=ΛEρ+mπ2. This relation fixes Eρ≈310 MeV from the averaged mass of charged ρ in the hadronic processes, mρ=766.5±1.1 MeV [3].ρ   1 Q 0 and ω1Q0 have quantum numbers JPC=1+- and, according to (4), mass m=Λ(Eρ+E0)+mπ2. Fitting to the masses of experimental states h1(1170) and b1(1230) [3], we obtain the estimate E0≈430-510 MeV.The stateπ0Q0 has the quantum numbers of scalar particle, JPC=0++. According to (5) its mass is m=ΛE0+mπ2≈900-980 MeV. This estimate is close to the mass of a0(980) [3].ρ   J Q 1 and ωJQ1 have JPC=(0,1,2)++. They should be the series of states a1(1260), f1(1285), f2(1270), a2(1320), and likely f0(1370) [3]. Substituting the experimental masses of well measured resonances f1(1285) and f2(1270) to our relation for their mass, m=Λ(Eρ+E1)+mπ2, we obtain the estimate E1≈570-580 MeV.ρ   1 Q 0 2 and ω1Q02 have mass m=Λ(Eρ+2E0)+mπ2≈1470-1590 MeV. They are the first radial excitations of ρ and ω mesons and describe the resonance regions ρ(1450) and ω(1420) [3].ρ   J Q 0 Q 1 and ωJQ0Q1 have quantum numbers JPC=(0,1,2)-+, i.e., 3 possible spins and pion parities. These states describe resonances π, π1, and π2 above 1 GeV. In contrast to the standard quark model, the π1-meson is not exotic in our scheme! The relation (5) predicts mass about m=Λ(Eρ+E0+E1)+mπ2≈1550-1610 MeV. The possible candidates are π1(1600) and π2(1670) [3].The statesρJQ12 and ωJQ12 have JPC=(1,3)-- as the total spin of the pair Q1Q1 can be 0 or 2. Natural candidates are ρ(1700) and ρ3(1690) (and ω(1650) with ω3(1670) for ω) [3]. Our fits yield m=Λ(Eρ+2E1)+mπ2≈1630-1650 MeV. As in the standard quark model, the exotic states with quantum numbers JPC=0-- are absent and ρ(1700) is not the second radial excitation (in the potential quark models, ρ(1700) (and ω(1650)) is aD-wave state, while the radial excitations of ρ(770) should beS-wave ones) of ρ(770). On the other hand, we do not obtain states with JPC=2-- which are predicted by the standard quark model but were not observed.


## 3. Regge Trajectories

Relations (4) and (5) lead to linear Regge, equidistant daughter Regge, and radial trajectories. Below we give examples for some of them.The statesπ0Q02n, n=0,1,2,…, form linear trajectory for the radial excitations of pion,(6)mπ2n=2ΛE0n+mπ2.The radial excitations of ρ(770), the resonances ρ1Q02n, lie on the first radial ρ-trajectory,(7)mρ2nI=2ΛE0n+Eρ2E0+mπ2.The second radial ρ-trajectory is composed of the states ρ1Q12Q02n,(8)mρ2nII=2ΛE0n+E1E0+Eρ2E0+mπ2.It is evident that the states ρ1Q12(k-1)Q02n formally give rise to the k-th radial ρ-trajectory. The resonances having structure ρ1Q1Q02n form the first radial a1-trajectory,(9)ma12nI=2ΛE0n+E12E0+Eρ2E0+mπ2.The expressions for further axial radial trajectories can be easily written. It should be remarked that the standard quark model predicts only two radial ρ-trajectories (S- and D-wave ones) and one a1-trajectory (a P-wave one). The scheme under discussion is much richer.The spinρ-trajectory is composed of the states ρJQ1J-1 with J=1,3,5,…. The corresponding masses are(10)mρJ2=ΛE1J-1+EρE1+mπ2.For the even spins, J=2,4,6,…, trajectory (10) describes aJ(fJ) mesons. The Regge trajectory (10) describes thus states with alternating parities and this agrees with the phenomenology.The principalρ-meson Regge trajectory (10) is accompanied by the daughter trajectories following with the step 2ΛE1. The spin-1 ρ-mesons of the kind ρ1Q1J-1 are the lowest states on the daughters. For example, the lowest state in (8) is the lowest state on the first daughter. The spectrum of ρ1Q1J-1 excitations reads(11)mρ1Q1J-12=2ΛE1k+Eρ2E1+mπ2,where k=0,1,2,… enumerates the daughters. Similar relations can be written for the axial and other mesons. An obvious consequence of the emerging spectrum is the degeneracy of spin and daughter radial excitations of the type m2(J,k)~J+k, which is typical for the Veneziano dual amplitudes [45] and the Nambu–Goto open strings. This kind of degeneracy was observed in the experimental spectrum of light nonstrange mesons [5–11].The fitted values ofEρ, E0, and E1 are rather close. This allows considering reasonable limits where some of them are equal. The notions of radial and daughter radial trajectories coincide in the limit E0=E1. In the limit Eρ=E1, the radial vector and axial trajectories are related by(12)ma12n=mρ2n+mρ2-mπ2.This relation holds both for radial and for daughter radial trajectories. In the chiral limit, mπ=0, relation (12) for the ground states, n=0, reduces to the old Weinberg relation, ma12=2mρ2 [46]. In the most symmetric limit, Eρ=E0=E1, the vector and axial radial spectrum in the chiral limit reduce to a very simple form,(13)mρ2n=2mρ2n+12,ma12n=2mρ2n+1.Relations (13) first appeared in the variants of Veneziano amplitude which incorporated the Adler self-consistency condition [45]. This condition (the amplitude of ππ scattering is zero at zero momentum) incorporates the CSB removing degeneracy between the ρ and a1 spectra. Within the QCD sum rules, relations (13) may be interpreted as a large-Nc generalization of the Weinberg relation [47].We see thus that in certain limits the Regge phenomenology of our approach reproduces various known relations. In addition, it is not excluded that the assumptionEρ=E0=E1 is compatible with the available data if a global fit is performed. According to the phenomenological analysis in review [2] the averaged slopes of spin and radial trajectories in light nonstrange mesons, a and a′ in (1), are equal; the reported value is a≈a′≈1.14GeV2. This remarkable observation leads to a large degeneracy in the spectrum [5–11, 21–27]. From the given value we have E0≈E1≈a/(2Λ)≈1.14/(2·1.83)≈0.31 GeV that coincides with our previous estimate for Eρ.Let us clarify further how the excited resonances with identical quantum numbers can have different origin in the proposed scheme. They may represent the radial states, states on daughter Regge trajectories and various “mixed” ones. For instance, the secondρ-meson excitation with the same parities appears in three forms: the second radial excitation ρ1Q04, the vector state on the second daughter trajectory ρ1Q14, and the mixed one ρ1Q02Q12. They are degenerate only in the limit E0=E1. It is likely difficult to detect such a splitting experimentally because of overlapping widths. In reality, one would observe rather a “broad resonance region”.We note also that placing of the observed “radial” states on a certain trajectory should be made with care: an incorrect interpretation of states leads to a false (or more precisely, introduced by hands) nonlinearity of the trajectory. Take again theρ meson as a typical example. The first three radial excitations of ρ are the states ρ1Q02, ρ1Q04, and ρ1Q06. They are accompanied by the following states with the quantum numbers of ρ: ρ1Q12, ρ1Q12Q02, and ρ1Q14. Since 3/2E0>E1>E0 in our fits, the sequence of first 7 ρ-mesons is ρ, ρ1Q02, ρ1Q12, ρ1Q04, ρ1Q12Q02, ρ1Q14, and ρ1Q06. They likely correspond to the vector resonances ρ(770), ρ(1450), ρ(1700), ρ(1900), ρ(2000), ρ(2150), and ρ(2270) [3].


## 4. A Theoretical Motivation

It is interesting to get a theoretical hint both on the rulemh2~Eρ,0,1 underlying the presented mass counting scheme and on the value of energy parameters Eρ,0,1. The latter point is intriguing by itself since the scale about 0.3 GeV is ubiquitous in the hadron physics; for instance, it often emerges in mass relations between hadrons containing heavy quarks [48].The mass of a hadron stateh can be related to the trace of energy momentum tensor Θμν in QCD (the sign convention is mostly positive),(14)2mh2=-hΘμμh,where Θμμ is given by the scale anomaly,(15)Θμμ=β2gsGμν2+∑q=u,d,…mqq¯q.Here β denotes the QCD beta-function and gs is the coupling constant. Relation (14) represents trace taken in the following Ward identity known in deep inelastic phenomenology [49]: 2pμpν=-hpμΘμνhpν. As follows from derivation of the identity (14), the squared mass (not the linear one!) in the l.h.s. appears due to the relativistic invariance. The r.h.s. of (14) is a renorminvariant quantity. One can build two substantially different renormalization invariant operators in QCD; both of these operators are present in (15). The l.h.s. of (14) defines the gravitational mass of a hadron. In reality, the vast majority of hadron resonances probably have no well-defined gravitational (and inertial) mass as they do not propagate in space: their typical lifetime of the order of 10-23–10-24s does not allow leaving the reaction area of the order of 1 fm which is comparable to their size. They show up only as some structures in the physical observables at certain energy intervals. The resonance mass is commonly associated with the real part of an S-matrix pole on the second (unphysical) sheet. How this definition is related to the gravitational mass is by far not obvious in a theory with confinement.In the sector of lightu and d quarks, there are only two hadrons with well-defined gravitational mass: the pion and nucleon. The pion mass is given by the GOR relation (2). This relation should somehow follow from the Ward identity (14). An explicit derivation of (2) from (14) would likely give an analytical proof for spontaneous CSB in QCD. Our present aim is to show, at least on a rather intuitive level, how the mass gap may emerge in nonperturbative gluon vacuum and get some quantitative insight. Consider the case of nucleon, h=N, for which (14) must yield the proton mass. Let us insert a complete set of vacuum states 00 from both sides of Θμμ in (14); the nucleon mass is then given by(16)mN2=-120∣N20β2gsGμν2+∑q=u,dmqq¯q0.This expression suggests that, in the chiral limit, mq=0, the nucleon mass is determined by the vacuum average 0Gμν20.From the one-loop QCD beta-function we haveβ/2gs=-β0/8αs/π, where β0=11-2/3nf, αs=gs2/4π. Relation (16) can be then rewritten as(17)mN2=120∣N2q¯q-∑qmq+β08αs/πGμν2q¯q.Relation (17) demonstrates that the masses of current quarks in nucleon acquire a contribution from nonperturbative gluon vacuum. Since nucleon at rest is interpreted in the quark model as a bound system of three constituent quarks, we can estimate from (17) an effective energy per quark, i.e., the value of constituent quark mass; in the chiral limit, it is just one-third of gluon contribution in (17):(18)Mc≃-β024αs/πGμν2q¯q.Substituting into (18) the standard value of gluon condensate from QCD sum rules, αs/πGμν2=0.012(3)GeV4 [42, 43], and q¯q=-(0.25GeV)3, nf=2, we obtain the estimate Mc≃310 MeV. It is a typical estimate for the value of constituent quark mass (strictly speaking, the value of constituent quark mass is very model dependent; as far as we know, it ranges from 220 to 450 MeV in various models. The value Mc(p2)≃310 MeV at small momentum p, however, proves to be seen in unquenched lattice simulations in the chiral limit (see, e.g., [50]). We can indicate a simple qualitative way leading to this estimate. When one incorporates the vector and axial mesons into the low-energy models describing the spontaneous CSB in QCD and related physics, e.g., into the Nambu–Jona-Lasinio model [51] or the linear sigma-model, a model-independent relation emerges: ma12=mρ2+6Mc2. On the other hand, the idea of CSB was also exploited in the derivation of famous Weinberg relation, ma12=2mρ2 [46]. Combining both relations, one has Mc=mρ/6. Substitution of the value mρ≈766 MeV used in our fits leads to Mc≃0.31 GeV), Mc≈mp/3, where mp is the proton mass.It is interesting to observe the numerical coincidenceMc≃Eρ. This suggests that the “spin flip” converting pion to ρ is essentially equivalent to creation of constituent quark. Thus the ρ meson in our approach may be interpreted as a system of one constituent quark and one current antiquark (or vice versa) which interact with the QCD vacuum. This picture is different from the usual quark model where ρ consists of two constituent quarks interacting with each other (in view of these speculations one can ask what is the exact relation for the nucleon mass within the given approach? We can propose the following guess. In the case of ρ meson, we motivated “a rule of transition” from the nonrelativistic quark model to our scheme: mh=∑Mc+interactions→mh2=Λ·1/2∑Mc+mπ2. Applying this “rule” to the ground state of nucleon (3 constituent quarks) we arrive at the relation mN2=3/2ΛMc+mπ2 that gives mN≃0.93 GeV. The factor 3/2 could be also qualitatively interpreted as follows: the creation of 3 constituent quarks and 3 constituent antiquarks (i.e., the creation of NN¯ pair in vacuum) is equivalent to creation of 3 constituent quark-antiquark pairs, so the relation of energies of constituents in nucleon and vector meson is 3 to 2.).The given interpretation can be further substantiated by the observation that the same “spin flip” should convert the nucleon toΔ baryon. As this flip “costs” ΛMc for the hadron mass squared, we should expect the relation mΔ2=mN2+ΛMc. It indeed yields mΔ≃1.2 GeV in a good agreement with experimental data for Δ(1232) [3].Coming back to motivation of our approach, the presented mass counting scheme is based on the assumption that the origin of hadron masses is similar to the case of nucleon mass in (17); the hadron mass squared (not the linear one as in many other approaches!) is proportional to effective energy of hadron constituents. This rule is conjectured from the Ward identity (14).The second assumption refers to the postulated form of meson constituents; the valentqq¯ pair plus constituent qq¯ pairs which effectively parametrize contributions to hadron mass from strong gluon field. The given choice is a model which works in the meson spectroscopy above and near 1 GeV. Besides the spectroscopy it could explain why the highly excited nonstrange mesons prefer to decay to more than two pions or, say, to have sometimes ρ-meson in the decay products instead of ππ-pair [3]; the corresponding possibilities are likely enciphered in the expression for the mass of a resonance as is seen in the examples of Section 2.It should be mentioned that hadron structure strongly depends on a reference frame. For instance, the proton structure experimentally looks very different if it probed for a proton at rest or for ultrarelativistic proton. It might be that the hadron model developed in this paper is more appropriate as a picture of light hadrons near the light cone while the traditional quark models refer to the rest frame. In this case there is no contradiction between different approaches.The constituentqq¯ pairs are not operative degrees of freedom noticeably below 1 GeV. It is interesting to note, however, that the spectroscopy in this sector can be constructed following essentially the same relation (3) if the pseudoscalar pair Q0 is replaced by the pseudo-goldstone bosons π, K, η [52].


## 5. Summary

We have proposed a novel approach to classification of mesons and description of meson spectroscopy. The approach represents a new realization of the quark model concept, a realization in which highly excited states appear due to multiquark components in hadron wavefunctions. We constructed a mass counting scheme that allows estimating meson masses by a trivial arithmetics with accuracy comparable to numerical calculations in complicated dynamical models. The given scheme is based on introduction of constituent scalar and vector quark-antiquark pairs. These pairs in excited hadrons, in some sense, bear a superficial resemblance to neutrons and protons in atomic nuclei.The inclusion of strange quarks into our approach is straightforward. It would be interesting to extend the presented ideas to light baryons and to hadrons with heavier quarks.The proposed approach is broader than “just another one model” as it provides a new framework for analysis of hadron resonances and can be used as a starting point for construction of essentially new dynamical models in hadron physics.



---

*Source: 1701939-2019-02-12.xml*