Excited mesons, baryons, glueballs and tetraquarks: predictions of the Holography Inspired Stringy Hadron model

In this note we collect and summarize the predictions of the Holography Inspired Stringy Hadron (HISH) model. Following a brief review of the model, we list the masses and widths of predicted excited states across the spectrum, based on placing the different hadrons on the non-linear Regge trajectories of a string with massive endpoints. Our predicted states include: (i) Light, heavy–light and heavy–heavy mesons. (ii) Baryons, including charmed, doubly charmed and bottom baryons. (iii) Glueballs, together with a method to disentangle them from flavorless mesons. (iv) Genuine tetraquarks, which are not “molecules” of hadrons, and are characterized by their decay into a baryon and an anti-baryon.


Introduction
The Holography Inspired Stringy Hadron (HISH) model [1] describes hadrons: mesons, baryons, glueballs and exotic hadrons in terms of bosonic strings in flat four dimensional space-time. The open strings have particles: "quarks", "antiquarks", or "diquarks" on their endpoints. The particles are massive and may also carry electric charges and/or spin. The basic open string configuration with massive endpoints describes a meson. The corresponding trajectories, which we refer to as HISH modified Regge trajectories (HMRTs), were determined in [2] and fitted to the observed meson spectrum. The mass of the quark as an endpoint particle m sep is neither the QCD current mass nor the constituent quark mass. A baryon corresponds to a string connecting a quark on one end and a diquark on the other, the latter being composed of a baryonic vertex with two short strings attached to it. The spectra of the baryons and their fits were worked out in [3]. Folded closed strings describe glueballs [4]. Exotic tetraquarks are formed of strings connecting baryonic vertices and anti-baryonic vertices [5]. On top of determining the spectra of hadrons the HISH model enables us to compute the strong decay width of the various hadrons, both the total decay width as well as some branching ratios of the different decay channels [6].
An important feature of the HISH model is that it provides a unified description of both light and heavy hadrons. The endpoint masses allow us to describe the HMRTs even of hadrons containing c and b quarks, with the same slope as for the corresponding light hadrons' trajectories. Only for the bottomonium we find there are deviations from the universal slope.
The HISH model may be viewed as a "renaissance" of the old Regge model 1 enhanced with several new features that follow from the string/gauge holographic duality. The model is built as a phenomenological model, and emphasizes comparison of its results with experimental data. Based on fits to data, we have extracted the best fit values of the basic physical parameters of the model which are the string tension (or the slope α ) and the masses of the endpoint particles m sep . In addition we fit all the intercepts a associated with the different trajectories. The intercepts however are not free parameters and in principle should be calculable in the model. With these parameters we can predict the existence of additional yet undiscovered hadrons.
In this note we summarize the predictions made in our previous papers on the HISH model, add additional predictions and further discussions, and point out certain issues which we believe are of special interest.
We list predictions for higher excited states for each of the leading meson trajectories, starting with the light mesons, the likes of the π and ρ, going through heavy-light mesons like the D and B and then all the way to heavy quarkonia including ϒ and η b . This is done for both angular momentum, for instance the states 7 −− and 8 ++ on the ρ/a trajectory, as well as the radial excitations, for instance n = 4 and n = 5 on the trajectory. Unexpectedly, according to the HISH model there are certain excited states of the kaon that have not been confirmed, for instance K with 3 + .
Similarly we predict the masses and widths of higher states for the different baryon trajectories. The light flavor baryons admit a peculiar behavior that the states with even and odd angular momentum are on separate parallel trajectories. The odd and the even states have different intercepts. We still do not possess a theoretical explanation for this effect. 2 As was discussed in [3] (Sect. 6.2) there are a priori ambiguities about the structure of the diquark of the HISH baryon. First one must determine the content of the diquark. For instance in a doubly charmed baryon with quark content of ucc it could be either (uc) or (cc). Second, the diquark masses as string endpoints have to be determined separately from the quark masses. For instance one may ask what is the mass of the (cc) diquark. Naturally it is expected to be close to 2m c . However, in holography this is far from obvious, and it can have any value starting from m c and up to 2m c or even higher. This depends on the location and tension of the baryonic vertex, which is a model dependent quantity [9,10]. We elaborate in this note about the various different predictions that follow from this situation. In particular we list the predictions of both options of the diquark in the c who is either (ss)c or s(sc), following its recently discovered first excited states. We also include predictions for excited bottom baryons: b , b , b , and b . Glueballs are described as closed strings, as described in works such as [11][12][13]. The lore about the glueballs is that there is no way to separate them from flavorless mesons. As was discussed in [4] since mesons are open strings and glueballs are closed string, there should be a way to distinguish between them. The latter have a slope which is half that of the former. The slope of glueball trajectories can also be measured on the lattice as we did in [4] or as in [14,15]. In addition the low lying states of closed string should have only even angular momentum and even parity and charge conjugation. To find the possible trajectories of glueballs we use four candidates for the scalar glueball: f 0 (980), f 0 (1370), f 0 (1500), and f 0 (1710). They are chosen as the ground states of the glueball trajectories, while the rest of the flavorless states are placed on nearly linear trajectories for light mesons and mass corrected trajectories of ss. For each of these cases we predict the masses of the n = 2 and n = 4 excited scalar glueball states. In case that these states will be found for a particular ground state one will be able to declare an identification of a glueball trajectory and hence also of the individual glueball states. We similarly discuss tensor glueballs and the spectrum of f 2 states. One supportive evidence for the existence of the glueball among these states is the fact that one cannot put all the flavorless states only on meson trajectories [4].
Another clue in identifying glueballs can come from their decays, as the closed string has to tear twice to produce two mesons unlike the open string which can tear only once. We work out the ratios of the widths for a decay of two vector mesons built from light quarks, for a decay into two K * and into two φ mesons. The presence of these three decay modes is a distinguishing feature of the closed string.
The fourth class of hadrons which we predict using the HISH model are tetraquarks. In [5] we stated that in principle there should be a big zoo of exotic hadrons built from strings connecting baryonic vertices, anti-baryonic vertices, with various combinations of quarks and antiquarks. The focus there was on the simplest such states, the tetraquarks. We argued that if a hadron decays predominantly to a baryon and anti-baryon, it is a signal of a tetraquark. In fact this idea was raised already in the early days of hadron physics [16] and has remained unresolved to this day. In recent years various experimental observations, specifically for hadrons including heavy quarks reignited this subject (see for instance [17][18][19]). We use the state denoted by Y (4630) that has such a decay, namely to c¯ c , as our prototype case. We have determined where its excited angular momentum and radial trajectory partners should be found. We predicted that a similar bottomonium-like state should have a decay process to b¯ b and its own trajectories of excited states. There may be also the hidden strange equivalent decaying to ¯ . The Y (4630) state and its bottom or strange analogues are all 1 −− states, and hence should be found in e + e − collisions. We collected in this note the data about the predicted masses and widths for such tetraquarks. We also suggest the existence of a trajectory as a way to distinguish between a genuine tetraquark and a bound state/molecule of two mesons. To find the trajectory one should naturally look above the corresponding baryon-antibaryon threshold, although narrower below threshold states can also be predicted.
The paper is organized as follows: After this introduction we bring a brief summary of the HISH model in Sect. 2. We then describe the fitting model used to compare the outcome of the theoretical model with the experimental observations in Sect. 3, where we write down all the formulas used to predict excited states. We discuss the parameters of our model, namely the Regge slope α , the various string endpoint masses m sep and the intercepts associated with the different trajectories. We also discuss the uncertainties of our predictions.
The following sections are where we list all our predictions. In Sect. 4 we discuss excited states on the meson trajectories. We have both orbital trajectories of states of high angular momentum and radial trajectories. For each trajectory we list the next two states. This is done for light mesons, light-heavy mesons and heavy-heavy charmonium and bottomonium. Similarly predictions about baryons are described in Sect. 5. This include predicted masses for the excited states of − baryons, charmed baryons ( c , c , c ), the doubly charmed baryon ++ cc , and the various bottom baryons: b , b , b , and b . Next we present predictions and a method of identifying glueballs in Sect. 6. Assuming that each of four candidates among the flavorless f 0 scalars, namely, f 0 (980), f 0 (1370), f 0 (1500), and f 0 (1710) are the ground state of the glueball trajectory, we predict the spectrum of excited glueballs that will follow. The glueball trajectories are characterized by a slope which is half that of the mesons. We also assign the rest of the flavorless states to normal meson trajectories. A similar procedure is applied also to tensor glueball states and the spectrum of the f 2 states. We also discuss the total width of the glueballs and also the ratio between various different decay channels. Section 7 is devoted to predictions about tetraquarks. We analyze the Y (4630) state that was observed to decay to c¯ c . We analyze the corresponding massive modified Regge trajectory and predict the states that reside on it. We perform a similar analysis for a bottomonium-like tetraquark, and tetraquarks in the light-quark sector. In the last part, Sect. 8, we summarize and outline several additional types of predictions that could be extracted from the HISH model.

A brief review of the HISH model
The holographic duality is an equivalence between certain bulk string theories and boundary field theories. The original duality is between the N = 4 SYM theory and string theory in Ad S 5 × S 5 . The N = 4 theory is obviously not the right framework to describe hadrons that resemble those found in nature. Instead we need a stringy dual of a four dimensional gauge dynamical system which is non-supersymmetric, nonconformal, and confining. The two main requirements on the desired string background is that it admits confining Wilson lines, and that it is dual to a boundary that includes a matter sector, which is invariant under a chiral flavor symmetry that is spontaneously broken. There are by now several ways to get a string background which is dual to a confining boundary field theory. For a review paper and a list of relevant references see for instance [20].
Practically most of the applications of holography of both conformal and non-conformal backgrounds are based on relating bulk fields (not strings) and operators on the dual boundary field theory. This is based on the usual limit of α → 0 with which we go, for instance, from a closed string theory to a theory of gravity.
However, to describe realistic hadrons we need strings rather than bulk fields since in nature the string tension, which is related to α via T = (2πα ) −1 , is not very large. In gauge dynamical terms the IR region is characterized by λ = g 2 N c of order one rather than very large.
It is well known that in holography there is a wide sector of gauge invariant physical observables which cannot be faithfully described by bulk fields but rather require dual stringy phenomena. This is the case for Wilson, 't Hooft, and Polyakov lines.
In the holography inspired stringy hadron (HISH) model [1] we argue that in fact also the description of the spectra, decays and other properties of hadrons -mesons, baryons, glueballs and exotics -should be recast as a description in terms of holographic stringy configurations only, and not fields. The major argument against describing the hadron spectra in terms of fluctuations of fields, like bulk fields or modes on probe flavor branes, is that they generically do not reproduce the Regge trajectories (or more precisely the HMRT) in the spectra, namely, the (almost) linear relation between M 2 and the angular momentum J . Moreover, for top-down models with the assignment of mesons as fluctuations of flavor branes one can get mesons only with J = 0 or J = 1. Higher J states will have to be related to strings, but then there is a big gap of order λ (or some fractional power of λ depending on the model) between the low and high J mesons. Similarly the attempts to get the observed linearity between M 2 and the excitation number n are problematic whereas for strings it is an almost trivial property.
The construction of the HISH model is based on the following steps. (i) Analyzing string configurations in holographic string models that correspond to hadrons. (ii) Devising a transition from the holographic regime of large N c and large λ to the real world that bypasses formally taking the limits of 1 N c and 1 λ expansions. (iii) Proposing a model of stringy hadrons in four flat dimensions that is inspired by the corresponding holographic strings. (iv) Confronting the outcome of the models with the experimental data (as was done in [2][3][4]) to measure the physical parameters and learn what additional features the stringy model should have.
Confining holographic models are characterized by a "wall" that truncates in one way or another the range of the holographic radial direction. A common feature to all the holographic stringy hadrons is that there is a segment of the string that stretches along a constant radial coordinate in the vicinity of the "wall", as in Fig. 1 [21,22]. For a stringy glueball it is the whole folded closed string that rotates there and for an open string it is part of the string, the horizontal segment, that connects with vertical segments either to the boundary for a Wilson line or to flavor branes for a meson or baryon. This fact that the classical solutions of the flatly rotating strings reside at fixed radial direction is a main rationale behind the map between rotating strings in curved spacetime and rotating strings in flat spacetime described in Fig. 1.
A key player of the map is the "string endpoint mass", m sep , that provides in the flat spacetime description the dual of the vertical string segments. It is important to note that this mass is neither the QCD mass parameter (the current quark mass) nor the constituent quark mass, and that the massive endpoint as a map of an exactly vertical segment is an approximation that is more accurate the longer the horizontal string is.
The stringy picture of mesons has been thoroughly investigated in the past (see [7] and references therein). In recent years there have also been attempts to describe hadrons in terms of strings. Papers on the subject that have certain overlap with our approach are [23][24][25][26][27][28][29][30][31][32][33][34]. A somewhat different approach to the stringy nature of QCD, at least conceptually, is the approach of low-energy effective theory on long strings reviewed in [35].
The HISH model describes hadrons in terms of the following basic ingredients: • Open strings which are characterized by a tension T , or equivalently a slope α = (2π T ) −1 ). The open string generically has a given energy/mass and angular momentum associated with its rotation. The latter gets contribu-tion from the classical configuration and in addition there is also a quantum contribution characterized the intercept a. When considering trajectories of the orbital angular momentum (not including the spin) as a function of M 2 it is a universal experimental fact that the corresponding intercepts are always negative a < 0. This property is responsible to the fact that the ground state of the stringy hadrons is not tachyonic. The negative intercept also provides a repulsive Casimir force that guarantees that even a non-rotating stringy hadron has a finite length. We do not have yet a theoretical explanation to this essential property of the hadronic spectra. The quantization of the HISH string including the determination of the intercept and the spectrum of excited states was performed in [36]. In holography the BV is built from a D p brane wrapping a p-cycle connected to flavor branes by N c strings. A priori there could be different metastable configurations of the N c strings. In the HISH we take that the preferable string configuration is that of a single long string and two very short ones. Since the endpoints of the two short strings are one next to the other we can consider them as forming a diquark. There are two arguments in favor of this string setup. It was shown that the Y-shape configuration, which is the most symmetric form when N c = 3, is in fact unstable and an introduction of a small perturbation to a rotating Y-shaped string solution will cause it to evolve into the string between quark and diquark. Secondly, had the Y-shape string been stable, the baryon trajectory slope α B should have been where α m is that of a meson. However as was shown in [3] the slope of the baryonic trajectory is within 5% the same as that of a meson trajectory. The setup of a holographic baryon and its map to HISH  Fig. 2. We emphasize that unlike other models which have quarks and diquarks as elementary particles (such as [37]), in the HISH model the diquark is always attached to a baryonic vertex. A BV can be also connected to any of more complex combinations of quarks and anti-baryonic vertices as discussed in [5]. • Closed strings which have an effective tension that is twice the tension of the open string (α closed = 1 2 α open ). They can have non-trivial angular momentum by taking a configuration of a folded rotating string. The excitation number of a closed string is necessarily even, and so is the angular momentum (on the leading Regge trajectory). The close string intercept, the quantum contribution to the angular momentum, is also twice that of an open string (a closed = 2a open ). These different properties of closed strings, the slope, the fact that the angular momentum has to be even, and the intercept, can serve as a way to disentangle between flavorless mesons and glueballs [4].
Hadrons, namely mesons, baryons and glueballs are being constructed in the simplest manner from the HISH building blocks.
• A single open string attached to two massive endpoint particles corresponds to a meson. • A single string that connects on one end to a quark and on the other hand a baryonic vertex with a diquark attached to it is the HISH description of a baryon. • A single closed string is a glueball.
• A single string connecting a baryonic vertex with a diquark on one end and an anti-baryonic vertex and antidiquark, is a stringy description of a tetraquark which may be realized in nature.

Fitting model and parameters
All the different hadrons described by the HISH model will be grouped into the HISH modified Regge trajectories of strings with massive endpoints (HMRTs). From the relation between angular momentum and energy of such a string we can derive the equations Here T is the string tension, β i is the velocity of the endpoint with the mass m i , and i the distance of the mass from the center of mass around which the endpoint particles rotate. The endpoint velocities β i are related to each other from the condition that the angular velocity is the same for both arms of the string, implying while the boundary conditions of the string imply The last two equations provide enough conditions that together with the equations for M and J they define the function J (M) which is the classical HMRT of the string with massive endpoints. To fit the mesons and baryons we add the quantum intercept to the classical expressions by taking J → J − a. To discuss also radial excitations we take J → J + n − a, where n = 0, 1, . . . is the radial excitation number. As we have shown in our more recent work [36], this is not the exact form of the quantum corrections, but it works well as a first approximation, especially considering we leave the intercept as a free parameter. Note also that we define the trajectory in terms of the orbital angular momentum and not the total one, and if there is additional spin in the hadron the intercept is shifted accordingly. In other words, all trajectories start from N = J orb + n = 0.
In the high energy or low mass limit where the endpoint masses are ultra-relativistic with β → 1, we can write an expansion in m E (we write here the symmetric case m 1 = m 2 only for simplicity): from which one can see that the linear Regge behavior is restored in the limit m → 0, and that the first correction is proportional to α m 3/2 E 1/2 . The Regge slope α is related to the string tension as usual by α = (2π T ) −1 . The opposing low energy limit, β → 0, holds when (E − 2m)/2m 1. Then the expansion is The HMRT describes then the interpolation from this highly non-linear behavior of J ∝ (E − 2m) 3/2 for high endpoint masses to the familiar linearity of J ∝ E 2 for light quark masses.
The width of a string state is given by the formula

7)
Here L is the string length, which is a function of the mass (or angular momentum) of the state and depends on string ten-sion T and endpoint masses. It also receives quantum corrections which are manifested as a dependence on a. In addition to it we have a phase space suppression factor , which we introduce to fit those states that are just above the threshold for decays. 4 The factor is given by where M 1 and M 2 are the masses of the outgoing states for a given channel. Note that for the predictions for higher states listed here we do not use this phase space suppression, taking it to be one, but it is used in some of the fits to the lower states from which the free parameter A is determined. This A is in fact the dimensionless asymptotic ratio of /M at high energies. For mesons it was found to be A ≈ 0.1 for most trajectories examined.
The statement that the decay width is proportional to the string length is true in a frame where the string is static. For the rotating string we have to adjust for the position dependent time dilation along the string, which affects the probability to tear and decay at each point: Therefore the length L of the string that enters in formula 3.7 is in fact adjusted for this time dilation and given by The effective L differs from the length of the string 1 + 2 by a factor which is a function of β i , and this factor goes from 1 for β i 1 to π/4 ≈ 0.79 as β i → 1. For the decays the intercept is also introduced by taking J → J − a. This will introduce a quantum "zero point length" of the string, and so will give a finite length even to those states with zero orbital angular momentum whose length vanishes classically. The rationale behind choosing this fitting formula for the decays was discussed at length in section 8.1 of [6]. Note that the decay widths in [6] were calculated for the leading (n = 0) orbital trajectories. In this note we assume the replacement J → J + n − a is sufficient to give estimates for the widths of states on radial trajectories, even though the dependence of the widths on n was not fully examined. In any case, we distinguish in the following the widths determined by our model as in [6] and more naive estimates where given.
For glueballs (closed strings), the situation is simpler since we do not have the endpoint masses. We use The slope of the closed strings is half of the slope of open strings α . Here J + n takes only even values since left and right moving modes must have equal total excitation number for a closed string. The length of the folded closed string is related to the mass via The effective length adjusted for time dilation that affects the decay width is simply so the decay width of a glueball is expected to be proportional to its mass. Note that by taking ∝ M we implicitly assume how quantum corrections affect the string length. For example, the ground state with J = n = 0 which has a classically vanishing length will have M 2 = −2a/α and thus a "quantum length" proportional to √ α |a|. Whether this is the correct way of defining the quantum length of the closed string is an open question.

Parameters
For the endpoint masses and quarks we take the values obtained from the fits in [2,3] and most recently in [6]. The quark masses we use are The heavy quark masses, m c and m b are rather strongly determined as their constituent quark masses. The mass of the c quark is half the mass of the lightest charmonium, and likewise for the b. This is what gives the optimal fits in those cases. In [2] the mass of the s quark was seen to be anywhere between 200 and 400 MeV, with some trajectories being fitted better by lower values of m s and other with higher masses. The additional measurements from the decay width fits of [6] are more consistent with the higher values in that range, so we use 400 MeV here. The light quark masses are not uniquely determined. They could be anywhere between zero and up to 100 MeV. The string tension is approximately T = (440 MeV) 2 so the ratio m 2 u /T is small in any case. The value 60 MeV used here is an average value that fits well both the spectrum and decay widths of the various mesons. The Regge slope is largely universal, with some variations. The most common value is common to all the meson trajectories in the (J, M 2 ) plane, with the bottomonium trajectories being a sole exception.
In the (n, M 2 ) plane the slope is typically somewhat lower, approaching 0.80 GeV −2 in some cases [2]. The slope shared by the baryon trajectories was seen [3] to be a little higher For each HMRT there is also the intercept, which we leave at this point as a free parameter, to be determined for each trajectory separately by fitting to formula 3.2.

Diquark masses
Our basic assumption is that the baryon can be described as a single string between a quark and a diquark. For the meson we only had to measure the quark masses and string tension, but for the baryons we also need to determine which of the quarks in the baryon form the diquark, and then also the mass of any possible diquark configuration. From the cases examined in [3] we only have measurements for diquarks containing at least one u or d quark. Our fits indicate that the mass of such a diquark is roughly the same as the mass of the heavier quark in the diquark, that is In other words, the baryonic vertex and the short string stretched from it to the u/d flavor brane do not appear to have a contribution to the endpoint mass that is sufficiently large for us to measure. To give an example of how that works, take the charmed baryon with the quark content udc. There are three a priori configurations for the quark-diquark in this baryon, u(dc), c(ud), or d(cu), but for all of them we will only see a string with the mass m c at one of the endpoints and m u/d at the other. This could also be an indication that the light quarks themselves should be nearly massless, though we do not use that conclusion where lone light quarks are concerned.
In this work we address also baryons which may have an (ss), (cs), or (cc) diquark. Unfortunately there are too few states from which we can measure the mass of these diquarks. However, our predictions are sensitive to the diquark masses, so any future measurements of the states we predict should help determining them. For these diquarks we assume, for the purpose of providing predictions and as a zeroth order approximation, that the quark masses sum: m (ss) = 2m s , m (cs) = m c + m s , etc. We note again that this seemingly  Table 14. In blue are the data points, and the red asterisks along the trajectory are the theoretical values trivial property is not fully apparent in the holographic models.

Uncertainties
The predictions presented below are not expected to be fully accurate. The great degree of precision to which some hadrons' masses are known is clearly not achievable in a model as simple as that of a string with massive endpoints. There is then a theoretical systematic error that is dominant when describing the spectrum through modified Regge trajectories alone. A certain deviation of the observed hadron masses from our predicted formulas was seen in our previous works, and is usually at most two to three percent in the mass, with a few exceptional points deviating more. This can also be seen in the plots of the fitted HMRTs in Figs. 3, 4, and 6. A similar deviation of a few percent can also be expected for the masses of the states that we predict below.
For estimating the decay widths we use less data and the deviations can be larger, in particular for the baryons where our model was less successful in describing the widths of the known states. For states on the orbital trajectories of mesons our calculated widths usually fall within the experimental range of error, but deviations are more common (see Fig. 5 for the fits of the meson widths). The predicted widths listed in Table 1 for higher excited mesons on the orbital trajectories are expected to be the most accurate. For states predicted outside that table, including the baryons and the radial trajectories of the mesons, the widths are provided as an estimate.

Predictions for meson states
In [2] we studied the spectrum of mesons, and showed that both light and heavy mesons can be placed on the modified Regge trajectories predicted by the HISH model. There are two types of trajectories, the orbital trajectories of increasing angular momentum, or trajectories in the (J, M 2 ) plane, and secondly the radial trajectories in the (n, M 2 ) plane, where n is the quantum radial excitation number. We list predictions for both types of trajectories. In [6] we calculated the decay widths of mesons on leading orbital trajectories, based on  Table 14. In blue are the data points, and the red asterisks along the trajectory are the theoretical values the linear relation between the decay width of a string state and its length. The fits done there are used to calculate the width of the mesons on these trajectories. In this note we add predictions for many states not found in [2], specifically including more mesons with heavy quarks. Table 1 collects all the predictions for the (J, M 2 ) trajectories, for higher J states. In Table 2 we list the predictions for the (n, M 2 ) trajectories, for highly excited states with fixed J PC . The orbital trajectories are plotted in Figs. 3 and 4. In Appendix A we include a list of all the observed states whose trajectories we use. Below we discuss some of the predicted states in further detail and point out specific open issues. The meson decay width fits of [6]. We plot the decay width (in some cases normalized by phase space factor as / ) as a function of the orbital angular momentum. In blue are the data points, and the red asterisks along the trajectory are the theoretical values. Some exceptional states were excluded from the fits and are marked in red

Light mesons
The trajectories of the light mesons, made of u, d, and s quarks are well established. The trajectories of mesons containing only u and d quarks are nearly linear, while the s quark adds a deviation from linearity thanks to its mass. In this sector our predictions are the surest, since the HISH model describes well both the masses and widths of light mesons, and they chiefly concern highly excited states.  Among the light mesons, the pion is notable in that it does not belong on a Regge trajectory. This can be explained as the pions can be thought of as the pseudo-Goldstone bosons associated with chiral symmetry breaking and as a result have an uncharacteristically light mass. On the other hand, the excited states of the pion, both orbital and radial fall neatly enough on HMRTs with the usual meson slope.
One point that remains unclear is the trajectory of the kaon. It is not obvious that the lightest kaons should lie on a trajectory, since they share the pseudo-Goldstone nature of the pions when considering three flavors. Given the first excited state K 1 (1270)1 + , the pseudoscalar kaon does appear to connect to it on a trajectory. The next state should be the K 2 (1770), but there is another reported state with J P = 2 − that is lower than that, K 2 (1580). Since the lowest mass states for a given J P should lie on the leading HMRT the existence of the K 2 (1580) is puzzling. Similarly, the PDG lists only the states K 3 (2320) for J P = 3 + and K 4 (2500) for 4 − where there should also be lighter states for these quantum numbers, as listed in our tables, which are yet unobserved. The natural parity K * mesons on the other hand fall very neatly on a trajectory, and there is no apparent reason for the K mesons to deviate too much from their own trajectory.
We have omitted from this section the discussion of the spectrum of the f 0 and f 2 resonances. The quantum numbers of these states, I = 0 with J PC = 0 ++ or 2 ++ , are the same as those expected for the lowest glueball states. For this reason, their spectrum is discussed in conjunction with the glueball spectrum in Sect. 6.

Heavy-light mesons
One of the strengths of the HISH model is that it provides a unified description of both light and heavy hadrons. The endpoint masses allow us to describe the massive modified Regge trajectories of hadrons containing b and c quarks, and the fits, whenever there are two or more mesons that can be placed on a trajectory, point to them having the same slope as the light mesons trajectory.
In previous work we have described the trajectories of the charmed mesons, including the decay widths of the D and D * s in [6]. For the mesons containing a b quark and a light or s quark we can add now several examples of pairs of states that still fall on trajectories with the usual meson slope, as can be seen in Fig. 4. For each of the B, B * , B s , and B * s mesons we find such a pair and list the next two orbitally excited states on the leading HMRT.
Based on our trajectories we can assign to some observed states their quantum numbers. The state D(2740) should be the 2 − state. Similarly B J (5970) is identified as the 2 − bottom meson.

Charmonium and bottomonium
The spectrum of heavy quarkonia has been studied extensively. For current hadron phenomenology it is specially important to know the full spectrum since many of the most compelling exotic hadron candidates are found among the charmonium and bottomonium states.
Here we start seeing some deviations from the universal Regge slope of the light and heavy-light mesons. The (J, M 2 ) trajectories of the cc still have the normal meson slope of approximately 0.9 GeV −2 but in the (n, M 2 ) plane we measured a slope of 0.60 GeV −2 [5]. For the bb the slopes of orbital and radial trajectories were 0.55 GeV −2 and 0.42 GeV −2 respectively [2]. Explaining the lower slopes of these trajectories is another important puzzle. From the point of view of the HISH model, these deviations could mean that we start to feel the departure from the approximate picture of the holographic string as a string with massive endpoints in flat spacetime. It is possible that for the heavy quarkonia we have to solve the full holographic string model. But one would also need to explain why the disparity between the orbital and radial slopes grows so large, especially for the charmonium trajectories.
Also interesting is the B c meson containing both charmed and bottom quarks. The two known states are the lightest of these, the B c meson with J P = 0 − at 6275 MeV as well as its first excited partner, also a pseudoscalar at 6842 MeV. The slope of a trajectory between them is 0.56 GeV −2 , which is between the cc and bb values. It would be interesting if more states are observed, and specifically the first J P = 1 + state which would allow us to measure the slope of the orbital trajectory.

Predictions for baryon states
In Tables 3 and 4 are the predictions for higher J and n baryons respectively, based on the results of our fits. The orbital trajectories have baryons across the spectrum, from the light ones to the heaviest charmed and bottom baryons. The radial trajectories concern only the N and states, where the trajectories were best established and the slope measured.
The fits done in [6] showed a larger deviation of the decay widths of baryons from our model's predictions when compared with the mesons, so the width predictions are less certain. In particular for some cases, like for the N and baryons, it looked like the width grows faster than what our model would predict based on linearity in the string length.

Light and strange baryons
The N and trajectories are well established with relatively many states on each trajectory [3]. One interesting effect that they exhibit, which is still not understood, is the splitting of the trajectories into two parallel lines for even and odd orbital angular momentum states. The strange baryons do not seem to exhibit this phenomenon. We list predictions for higher spin light and strange baryons in Table 18. There are also many N baryons which fall neatly on nearly linear radial trajectories, and predictions for higher states for them are in Table 4.
Moving on to the spectrum of the sss, the excited baryons could have helped us determine the mass of the ss diquark directly, but there is not much data on the excited states. Our model predicts the first orbitally excited state of the to have a mass of around 2050 MeV, but there is no confirmed state there. The states that are known -(2250), (2380), and (2470) -could be higher excitations. The  (2380) is at the right mass for the 7 2 + state. One possibility is that the trajectories exhibit the same even-odd effect as the N and . Then, one might claim that there is no missing state, and the (2250) as the 5 2 − state belongs on a separate parallel trajectory of the odd states. On the other hand, this would make the even-odd effect larger for the than it is for the N and . Spin-orbit interactions which would be related to this splitting are suppressed for higher quark masses, so this possibility is less likely.

Charmed baryons
The c and c trajectories were examined in [3]. We have no other measurement of the mass of either diquark, so this adds further ambiguity, but from fitting the two chosen pairs of states, the s(sc) configuration is preferred. We list in the table the predictions of either option.

Doubly charmed baryons
For the doubly charmed baryons, only the ground states + cc and ++ cc are listed by the PDG. The status of the former is very uncertain, so we only use the ++ cc whose measured mass is 3621.4 ± 0.8 MeV.
A priori the diquark in a doubly charmed baryon can be either (uc)/(dc), or it can be (cc). In the former configuration we will see masses of m c = 1.5 GeV on each side, while in the later we will have one light endpoint and one heavy one. The mass of the (cc) diquark is naturally expected to be close to 2m c , although that is not obvious in the holographic picture where the mass of the diquark is determined by the vertical segments of the strings and the tension of the baryonic vertex. The predicted mass for the next state depends strongly on the choice of configuration, and we list both options in Table 3. The configuration with m c on either side gives lower masses of the excited states. For the orbital trajectory we take the usual slope of 0.95 GeV −2 .
We may estimate the width by comparison with the single charm baryons. The first excited state of the cc will be just heavy enough to decay to the ground state and one or two pions, so we can expect a width of a few MeV (for most direct comparison see the charmed-strange baryon c (2815) 3 2 − whose width is slightly less than 3 MeV). The second excited state will be above the threshold to decay to a charmed baryon plus charmed meson, c D, which would be the preferred channel unless the diquark is (cc). A decay width of around 20 MeV would be a reasonable estimate for the first above threshold state.

Bottom baryons
The spectrum of excitations of the bottom baryons is largely unexplored. The search and identification of spectrum of bottom baryon excitations is ongoing [40][41][42][43]. The only confirmed excitation is a b state, 0 b (5920) 3 2 − , first observed in [44,45]. It and the ground state 0 b (5620) 1 2 + are again connected by a line whose slope is α = 0.95 GeV −2 , as for the other baryons, as can be seen in Fig. 6h. The other bottom baryons do not have yet their excited partners. We will assume they can also be accommodated on trajectories with the same slope, so we can predict where the excited states will be for the b , b , and b . For the latter, first excited ssb state could be narrow, like the equivalent c . Here the predicted mass does not depend very strongly on whether we have the (ss)b or s(sb) configuration, but the latter gives a lighter mass for the first excited state. Recently, a b state was observed for the first time, b (6097) [46]. It has a mass slightly higher than the 6060 MeV predicted by us for the first orbital excitation of the b , but can certainly be identified as the 1P state. Another recent observation from LHCb is an excited − b resonance [47] at a mass of 6226.9 ± 2.4 MeV. There it is speculated to be either the first orbital (1P) or first radial (2S) state. According to our model the first orbital excitation should be much lower than that, and in fact the measured mass is closer to our prediction for the next excitation (1D). If the b (6227) is the first orbital excitation, then the mass of the radial trajectory connecting it and the ground state is 0.55 GeV −2 , very similar to the slope we find for the radial trajectories of heavy quarkonia. If we assume that this is the characteristic slope of the bottom baryons' radial HMRTs, we can predict the 2S states for other bottom baryons. Our predictions for these states are in Table 4. Fig. 6 The leading orbital trajectories of the baryons fitted with the common slope α = 0.95 GeV −2 . All the states fitted are listed in the appendix in Table 18  Table 5 The different assignments of the f 0 into radial trajectories and predicted higher states. The slope α (in units of GeV −2 ) was fitted for each assignment separately, as was done in [4], but is common to all three types of trajectories. Widths are provided as estimates, based on proportionality of the width to the string length

Predictions for glueball states
In [4] we have analyzed the spectra of the isoscalar resonances f 0 (J PC = 0 ++ ) and f 2 (2 ++ ). In Appendix A we list all the known f 0 and f 2 states. It is not too unlikely that one of the known f 0 resonances is in fact the glueball scalar ground state. In our paper we have proposed to organize the different resonances into the Regge trajectories of mesons and those of glueballs. To identify the glueball, one should measure the excited states.

Scalar glueballs
In [4] four candidates for the glueball ground state were examined: f 0 (980), f 0 (1370), f 0 (1500), and f 0 (1710). We analyzed the spectrum of the f 0 resonances by assigning them to radial trajectories. There were three types of trajectory, nearly linear trajectories for light mesons, mass corrected trajectories of ss, and the linear trajectories with half the meson slope for the glueballs. For each glueball candidate, we assign the rest of the f 0 to meson trajectories and predict where the next states should be. Table 5 is a summary of the different possibilities. Each group of three rows is a different assignment, and we list the next two predicted states on each trajectory. Note that trajectories begin with n = 0 and that for the glueball, only even n are allowed. Recurring in the table is the state f 0 (1800), whose nature and status are still unclear. For most of these assignments to be complete this state is needed in one of the meson trajectories.
Among the conclusions of [4] were that an assignment with meson trajectories alone is not fully consistent with all data, and that the f 0 (980) can be a light glueball, or it can be excluded from the trajectories as a KK bound state, which many suppose it is. The assignments that were most consistent were those with either f 0 (1370) or f 0 (1500) as the glueball ground state.
We also include an estimate for the widths of the excited states. When calculating the width of a glueball, we assume the simplified relation /M = Const., 5

Tensor glueballs
The first excited state of the glueball is expected to have J PC = 2 ++ . Again, we have many f 2 resonances with the quantum number expected of the glueball. Sorting the f 2 states into trajectories, the situation is somewhat simpler than with the scalar, as here we have two states that belong on meson trajectories in the (J, M 2 ) plane. The f 2 (1270) belongs to the trajectory of the ω meson, and the f 2 (1525) is ss and sits on the φ trajectory. Their decay modes and other properties are also well known and there is no real doubt about their nature. We take these two states as the heads of trajectories, and sort the remaining states.
In [4] we organized the f 2 resonances into the trajectories shown in Table 10. There are the meson trajectories headed by f 2 (1270) and f 2 (1525), as well as another pair of states on a light quark trajectory. Among the remaining f 2 states listed by the PDG are the f 2 (2010), which was observed to decay to φφ (despite the very small phase space). It could perhaps be classified as a φφ bound state, in an analogous fashion to the f 0 (980) being KK . The most interesting states after that remain the f 2 (1430) and f J (2220). While the latter has been considered a candidate for the glueball and has been the object of some research (see many papers citing [49]), the former is rarely addressed, despite its curious placement in the spectrum between the lightest 2 ++ light and ss mesons. It seems a worthwhile experimental question to clarify its status -and its quantum numbers, as the most recent observation [50] can not confirm whether it is a 0 ++ or 2 ++ state, a fact which led to at least one suggestion [51] that the f 2 (1430) could be itself the scalar glueball.
The predictions for the 2 ++ glueball were given in Tables 6, 7, 8 and 9, assuming the tensor is an excited partner of the scalar ground state. Then it is expected in the range 2385-2770, depending on which is the scalar ground state. In Table 10, we give the predicted states of 2 ++ mesons in Table 10 The meson trajectories of f 2 as discussed in [4], with their excited states. The fitted slope here is 0.85 GeV −2 , and half this value is used for the glueball trajectories. Widths are provided as estimates, based on proportionality of the width to the string length. However, the glueball candidates have very small widths, 13 MeV for the f 2 (1430) and 23 MeV for the f J (2220). The predicted widths for the excited states based on our model are between 30 and 50 MeV, but it is unlikely that they remain so narrow that range. There are multiple meson states expected in that range as well as glueballs. We also list the predicted excited states for the glueball candidates f 2 (1430) and f J (2220), assuming they are the head of a new trajectory with half the meson slope.

Glueball decays
In [6], we have used several measurements from meson decays to measure the strangeness suppression factor λ s , which is defined as the relative probability to create a pair of ss compared with a pair of uū or dd. In holography pair creation is realized by the string fluctuating in the radial direction to reach a flavor brane. The probability is suppressed exponentially in the distance from the wall to the flavor brane, which itself is proportional to the mass of the quark. Namely there is a suppression factor of the form exp(−cm 2 q /T ), for some model dependent constant c [52]. This type of exponential suppression in m 2 q is familiar also from Schwinger pair production in the QCD flux tube [53].
The cleanest measurements of λ s in [6] were from the radiative decays of the cc and bb. They are also most relevant to the glueball decays since we can think of radiative decays as the string closing and annihilating the quarks, then tearing again to create lighter quarks. The measured result was λ s = 0.36 ± 0.12.
For the glueball, the closed string must tear twice to decay, and so we can create zero, one, or two pairs of s ands in the decay. The glueball could then decay to a pair of light quark mesons, a pair of a strange and an anti-strange meson, or a pair of ss, with each creation of an ss pair adding a factor of λ s . The presence of these three types of two body decays, with the given hierarchy between them, is unique to closed strings, and hence could be useful in identifying glueballs.
To give a numeric example of the expected ratios, consider a glueball candidate of a mass of 2.2 GeV, where we assume a specific value of the mass in order to evaluate the phase space factors explicitly. A glueball of this mass could be an excited scalar glueball, or a tensor [4]. It is above threshold for decays into two ss mesons. So we can predict the existence of several decay channels and the relative ratios between them. For instance, for decays into a pair of chargeless vector mesons, we predict the ratio between the decay widths Taking λ s = 0.36, then λ 2 s = 0.13 and the numbers are If the glueball is a higher mass than 2.2 GeV, then the phase space is less restricted. Even in the 2.2 GeV case, the ratio of 7% to create the φφ means all three decay channels should be observable. Other predictions for holographic glueball decays are in [54][55][56][57][58].

Predictions for tetraquark states
Our paper [5] discussed at length the stringy tetraquark, the system of a diquark and an anti-diquark joined by a string. Like all hadrons, these states should lie on HMRTs. We discuss in particular the charmonium sector, since there we find a candidate for exactly such a state. We predict more of these states, which should be identifiable by their decays to a baryon anti-baryon pair. The focus on exotic hadron phenomenology is among the charmonium and bottomonium states [18,59,60]. In Sect. 4.3 we discussed the spectrum of quarkonia.

Tetraquarks in the charmonium sector
Of particular interest to us is the state Y (4630). It was observed by Belle in 2008, in the process e + e − → γ + c − c , with a significance of 8σ [61]. The parameters measured there for the Y (4630) were Near the Y (4630), there is another reported resonance, the Y (4660). It was observed by both Belle [62] and BaBar [63], also has J PC = 1 −− , and its measured mass and width are The channel in which it was observed is e + e − → γ π + π − (2S). It may well be that the two resonances Y (4630) and Y (4660) are really one and the same. The PDG lists the Y (4630) and Y (4660) measurements under a single entry for the Y (4660). What we have then is a single state with both decay channels, c c and (2S)π + π − . If this is the case, then the c¯ c is still the dominant decay, so this does not interfere so much with the tetraquark picture [64,65]. where there are two different slopes, one for orbital trajectories in J and one for radial trajectories in n. By extrapolating the trajectory backwards, to lower masses, we can check if the Y (4630) is an excited state of an already known resonance. In the angular momentum plane, the preceding state would be a scalar J PC = 0 ++ . There is no currently known scalar resonance at the mass predicted by the HMRT, of about 4480 MeV. In the (n, M 2 ) plane we find the Y (4360), another vector resonance that is at the right mass for the Y (4630) to be its first excited state. Its mass and width are M Y (4360) = 4354 ± 10, Y (4360) = 78 ± 16. (7.4) It was observed in decays to (2S)π + π − like the Y (4660).
The problem with assigning the Y (4360) as the unexcited partner of the Y (4630) is in its width. The widths of the two states are roughly the same even though the lower state, being beneath the c¯ c threshold, should be significantly narrower than the above threshold state. Whether or not there is a lower, unexcited version of Y (4630), we still predict higher excitations of the stringy tetraquark. For the higher states, the c¯ c decays should be dominant, and their masses should be such that they fall on a HMRT. The predictions are in Table 11.

Tetraquarks in the bottomonium sector
In the spectroscopy of heavy quarkonia, one usually expects analogies to exist between the cc and bb spectra. Tetraquarks and other exotics should be no exception. Therefore, we propose to search for a state directly analogous to the Y (4630), which we will call the Y b , that decays primarily to b¯ b . We will assume that, like the Y (4630), the Y b state is located a little above the relevant baryon-antibaryon threshold of b b .
The mass of the 0 b is 5619.51±0.23 MeV. As an estimate -take a mass of The Y b would also have its own Regge trajectory, so, based on the mass we choose for it, we can predict the rest of the states that lie on the trajectory. The trajectories are calculated using the slopes of the ϒ trajectories, which are α ≈ 0.64 GeV MeV [2]. The resulting masses are listed in Table 12.
The masses of the states that would precede the above threshold Y b tetraquark are also listed in Table 12. The observed resonance Y b (10890) is a potential match to be the "n = −2" state. This resonance was observed to decay to π + π − ϒ(nS) and has a mass of 10888.4 ± 0.3 [66][67][68]. It seems analogous to the Y (4360) from the charmonium sector, which decayed to π + π − (2S) and was at the right mass to be an unexcited partner of the tetraquark candidate  Y (4630). If the Y (10890) is on the trajectory, then we also predict another resonance following it, which being similarly below the baryon-antibaryon threshold should have similar properties, at a mass of ≈ 11080 MeV.

Lighter tetraquarks
One of the puzzles in exotic hadron spectroscopy is the absence of exotics in the light quark sector. We propose to search for the hidden strange analogue of the Y (4630), a vector resonance decaying to ¯ , and the excited states on its trajectory. We predict again a trajectory of states beginning with a near-threshold one. The masses are in Table 13. We take a ground state mass of with the usual slope for light quark (u, d, s) hadrons of 0.9 GeV −2 in the (J, M 2 ) plane. In the (n, M 2 ) the slope is lower, and we take 0.8 GeV −2 . The mass of the s quark is taken to be 400 MeV. Again we repeat the exercise of extrapolating the trajectory backwards, listing two preceding states. The vector resonance ρ(1570) is at the right mass to be a match. We note that while the natural way to search for molecules is by looking in the vicinity of different two hadron thresholds [69], a good way to identify tetraquarks is to look above the two baryon threshold. It is there that we can look for particles that decay into a baryon anti-baryon pair, the stringy tetraquark's natural and distinctive mode of decay, but it is also a region where data is still quite scarce. And while the predictions above are for the Y (4630), and its immediate analogues Y b and Y s , it should be noted that there could be similar states with non-symmetric decay channels, like c c for other hidden charm tetraquarks, c being charmed and strange, etc.. For resonances found in such channels we will also expect trajectories, which can be easily delineated using the HISH model.

Summary
The Regge trajectories are a central feature of the hadron spectrum, and are still widely used today to characterize the spectrum of hadrons [70][71][72][73][74][75][76][77][78][79], although these descriptions usually deal only with light hadrons. Not only do we have the linear or near-linear trajectories of the light hadrons such as the ρ meson or the nucleons, but the picture also extends to the HISH mass modified Regge trajectories of hadrons including heavy quarks. Charmed and bottom mesons can be sorted into trajectories like the light mesons, with the endpoint masses of the quarks describing the deviation from linearity of these trajectories. The Regge slope (or equivalently string tension) measured for these trajectories is almost always the same for heavy hadrons as it is for the light ones, which is strong evidence for the validity of the stringy picture for all the different types of hadrons. Furthermore, the stringy picture appears to remain valid even when describing low spin, unexcited states, contrary to the usual lore that is should hold only for excited states long flux tubes form, which can then be given an effective string description.
Using the measured endpoint masses of the quarksm u/d ≈ 60 MeV, m s = 400 MeV, m c = 1490 MeV, and m b = 4700 MeV -and the (near) universal Regge slope of α = 0.88 GeV −2 for mesons or 0.95 GeV −2 for baryons, it is easy to place various hadrons on trajectories, find the intercept for each trajectory, and then predict higher excited states by extrapolating the line. Given the slope and endpoint masses, we can use any single state to predict the spectrum of its excited partners. In Tables 1 and 2 we collected all our predictions for excited meson states. Tables 3 and 4 collect the results for baryons.
Some exotic hadrons are also expected to behave as strings. The glueball can be described as a folded closed string and as such will have a trajectory with half the slope of the mesons. A tetraquark can be thought of as a string connecting a diquark to an anti-diquark, and will also have a trajectory like that of the baryons. In both these cases, it could be possible to identify the exotics by their trajectories. And this is done not by the properties of a single state, but by finding a whole spectrum of exotic states. For the glueballs the first excited states are expected in the range of 2.5-3.0 GeV, where not much is known. The excited glueballs are expected to have a characteristic decay pattern as described in Sect. 6.2. To identify tetraquarks we have suggested to look for states above the threshold for baryon-antibaryon decays of the tetraquark, for instance above the ≈ 4570 MeV threshold for c c decays to identify charmoniumlike tetraquark candidates. Our predictions determine where these states should be found. If one accepts diquarks as natural objects in QCD, as they appear in all baryons in our model, then the tetraquark states must surely also be part of the spectrum.
The observation and understanding of higher excited states is crucial if we want to gain more insight into the stringy nature of hadrons. In particular, measurements excitations of heavy and doubly heavy baryons will teach us more about diquarks and the diquark structure of the baryons. In this note we assumed that a naive approximation of the diquark mass (as a string endpoint, as always) is m (q 1 q 2 ) ≈ m q 1 +m q 2 , even though in the holographic picture there is no reason that this simple sum rule should hold. It would certainly be interesting if one would discover significant deviations from this.
A natural question related to the predictions of the HISH model is what are the priorities in attempting to verify the predictions. It seems to us that the highest priority should be given to detecting glueballs and disentangling them from flavorless mesons. The main point here is that to identify them it is not enough to detect one state but rather one has to find, preferably three or more states that reside on a glueball trajectory. Only then by identifying the correct slope (of around 0.45 GeV −2 ) we will be able to declare a glueball state. The second priority we would assign to verifying or disproving the HMRT associated with tetraquarks where certain of the states residing on it decay predominately to a baryon and an antibaryon. As mentioned above the Y (4630) is such a state and it is important to verify the existence of its HMRT partners. There are also potential narrow below thresholds state that may belong to the trajectory, which may be seen in the (nS)π π channels. In particular, the Belle II experiment is well suited to observed this type of resonances. The next priority is to detect the HMRT states that include charm and bottom flavors and in particular baryonic configurations that have not been observed so far with three charmed quarks and two and three bottom quarks.
Among the challenges still facing the HISH model are • As shown in Appendix B the measured intercept, defined by the dependence of the orbital angular momentum on M 2 , is always negative. The first determination of the intercept in the HISH model was carried out in [36].
In symmetric systems where the two endpoints masses are the same the outcome intercept was found out to be always positive. For asymmetric strings, for instance for one massless and one massive endpoint the intercept can be negative. However it is obvious that this cannot explain the negative intercept of mesons with identical quarks at their ends. The analysis of [36] did not include endpoint electric charges and spins. The latter change the boundary conditions of the string and hence the eigenfrequencies of fluctuations, and thus also the intercept. • The basic structure of the HISH model is independent on the particular holographic background that is implemented. However, there are certain properties of the model that do vary from model to model. One such a property is the location along the holographic radial direction of the baryonic vertex [3,9,10]. This affects the determination of the mass of the diquark. As was mentioned above, our current understanding of the HISH model does not allow us to determine unambiguously the mass of a diquark built from two identical quarks. The mass of the diquark can vary between the mass of the quark and twice this mass, and can be heavier if the baryonic vertex contributes more than the strings to the endpoint mass.
In addition the HISH model cannot yet determine the structure of the diquark in the cases that there are several options for the structure. These two ambiguities of the HISH model deserve further investigation. • Since for light-light and heavy-light hadrons the HISH model with a universal slope fits nicely the experimental data, it seemed that the exploration of the differences between the HISH model and the stringy holographic model is not so crucial. Not unexpectedly for heavy quarkonia there are significant deviations from a universal slope. This obviously calls for a quantitative study of holographic hadrons and the limitations of the HISH model. A related question is to try and understand the differences between the orbital and radial slopes and between the baryon slope and the one for mesons. • Here in this note we have included only predictions about tetraquark exotic hadrons. However, as was discussed in [5] there is in principle a much richer spectrum of exotic hadrons. For instance one can add to each baryonic vertex a pair (or several pairs) of strings with opposing orientations, namely connecting it to a quark and an antiquark. This will create a pentaquark (or higher multi-quark) states. These, if exist, are genuine exotic hadrons and not molecules of several hadrons. • A major challenge for the HISH model is to determine scattering amplitudes of hadrons from their stringy pictures, and to compare them to the measured ones. This subject is currently under investigation.

A The states used in the fits
In this appendix we provide tables with all states used in our fits, and which form the basis for the predictions presented in this note. All states are as listed in the Review of Particle Physics published by the Particle Data Group (PDG) [80]. f.

Table 16
All the f 0 states as listed by the PDG. The last few states, marked here by asterisk, are classified as "further states"     Table 19 The baryon states used in the (n, M 2 ) trajectory fits and their assignments. The fourth column indicates the status given to each state by the PDG, ranging from * for unconfirmed states to **** for those whose properties are certain  In Tables 14 and 15 we list all the meson states that we placed on HMRTs, for orbital and radial trajectories respectively. The selection and identification of states was discussed in detail in [2]. Separately, in Tables 16 and 17 are listed all the flavorless f 0 (J PC = 0 ++ ) and f 2 (2 ++ ) states, among which the glueball could be found. Tables 18 and 19 list all the baryon states used in the fits. The selection of states was discussed in [3].

B The intercept
One of the most important challenges for the HISH model is to have a full theoretical computation of the intercept of the modified Regge trajectory. In [36] we have determined the leading order quantum corrections to the trajectory of a string with massive endpoints for long rotating strings. However, this is a correction to the result a = 1 of the usual boson string theory. To describe real world hadrons we need the intercept to be negative (recall that we define the intercept w.r.t. the trajectory of the orbital angular momentum as a function of the mass, and therefore all trajectories start from J or bital = 0 and the intercept is always negative).
We also wish to understand the dependence of the phenomenological intercept as a function of the endpoint masses. Below, in Table 20 we list the fitted values for the intercepts of the various orbital HMRTs of mesons described in this work, arranged by the quark contents of the meson. The values listed do not give a very clear picture of the mass dependence. It should also be noted that the intercept is very sensitive to changes in the quark masses and slope, and much more so percentage-wise than any of the predictions for the observables, the hadron masses or widths, given throughout this note.