A systematic study of mass spectra and strong decay of strange mesons

The mass spectrum of the kaon family is analyzed by the {modified} Godfrey-Isgur model with a color screening effect approximating the kaon as a heavy-light meson system. This analysis gives us the structure and possible assignments of the observed kaon candidates, which can be tested by comparing the theoretical results of their two-body strong decays with the experimental data. Additionally, prediction of some partial decay widths is made on the kaons still missing in experiment. This study is crucial to establishing the kaon family and searching for their higher excitations in the future.


I. INTRODUCTION
As an important part of the meson family, the kaon subfamily has become more and more abundant with experimental progress on the observations of kaons in the past decades. Until now, Particle Data Group (PDG) has collected dozens of kaons [1]. When facing so abundant kaons, it is one of the main tasks of the present study of light hadron spectroscopy how to categorize them into the family and another task is to investigate its higher radial and orbital excitations.
Before the present work, there were some theoretical papers related to kaons. For example, thirty years ago, Godfrey and Isgur [2] developed a relativistic quark model, the socalled Godfrey-Isgur (GI) model, by which they studied the mass spectrum of hadrons including kaons. In 2002, Barnes et al. [3] further investigated the strong decays of the observed kaons, which have masses less than 2.2 GeV, where the 3 P 0 quark model associated with a simple harmonic oscillator (SHO) wave function was adopted in their calculation. In 2009, Ebert et al. [4] analyzed the mass spectrum and Regge trajectories of kaons by their relativistic quark model.
Due to the present experimental progress on kaons, it is a suitable time to systematically carry out phenomenological study of kaons. In this work, we first calculate the mass spectra of the kaon family by applying the modified GI model [5,6], where the screening effect is taken into account. Fitting some well established kaon states, we fix the parameters in the model, which are adopted when calculating the masses of other kaon states. Comparing theoretical results with experimental data, we obtain the structure information of the discussed kaons. Especially, we predict some radial ground states of kaon which are still missing in experiments, e.g. K 4 (2310)(1 1 G 4 ). Using our potential model approach, the spatial wave functions of the kaons studied can be numerically calculated, which we take as input when studying their * Corresponding author † Electronic address: xiangliu@lzu.edu.cn ‡ Electronic address: matsuki@tokyo-kasei.ac.jp Okubo-Zweig-Iizuka (OZI)-allowed two-body strong decays. For further testing the properties of the kaons, we study their OZI-allowed two-body strong decays, which provide valuable information of their partial and total decay widths, where one uses the quark pair creation (QPC) model which was proposed in Ref. [7] and extensively applied to studies of other hadrons in Refs. . Analyzing mass spectra and calculating strong decay behaviors, we finally identify their n 2S +1 L J quantum numbers, which reflect the inner structure of the kaons under discussion. Here, we predict the strong decay behaviors of some kaon states, e.g. K 4 (2310)(1 1 G 4 ) has a wide width about 710-880 MeV, and mainly decays into K * 4 (2045)π, K * 3 (1780)π, Kρ 3 (1690) and Ka 2 . The study presented in this work is helpful for establishing the kaon family by including more higher radial and orbital excitations. This paper is organized as follows. After Introduction, in Sect. II we explain the modified Godfrey-Isgur model and the QPC model. In Sect. III, we adopt the modified Godfrey-Isgur model by including the screening effect to study the mass spectra of the kaon family. Making a comparison between theoretical and experimental results, we further obtain the structure information of the observed kaons. In Sect. IV, we present the detailed study of the OZI-allowed two-body strong decays of the discussed kaons. The paper ends with conclusions and discussion.

II. PHENOMENOLOGICAL QUARK MODELS ADOPTED IN THIS WORK
In our calculation, two phenomenological quark models are adopted, i.e., the modified GI model with the color screening effect 1 , and the QPC model. The modified GI quark model 1 When studying the mass spectrum of mesons, there are approaches like Dyson-Schwinger and Bethe-Salpeter equations, which are directly related to QCD. However, such theory-based or theory-linked approaches still have some limitations for describing higher excitations of mesons. Instead of theoretical approaches, one may apply phenomenological models to deal is applied to calculate the mass spectrum of the kaon family, by which we obtain the structure information of the observed kaon candidates. Then, we further test the possible assignments by comparing the theoretical results of their twobody OZI-allowed decays with the experimental data, where the QPC model is used to calculate their strong decays. In the following, we will introduce these two models.
A. The modified GI model First, we introduce the Godfrey-Isgur (GI) relativized quark model and discuss how the GI model is modified by including the color screening effect. Below we describe the detailed procedure and equations actually done by us because those are necessary in our work but are not familiar to general readers. Some are common to Godfrey and Isgur. The interaction between quark and antiquark in the GI model [2] is described by the Hamiltoniañ where m u/d and m s are the masses of u/d and s quarks, respectively, i.e., m u = m d = 220 MeV, m s = 419 MeV. V eff (p, r) =H conf +H hyp +H SO is the effective potential of the qq interaction which can be obtained from on-shell qq scattering amplitudes in the center-of-mass (CM) frame [2] and relativistic effect corrections. The quantities with tilde will be defined later. On the other handṼ eff (p, r) also consists of two main parts. The first one is a γ µ ⊗ γ µ short-distance interaction of one-gluon-exchange and the second part is a 1 ⊗ 1 longdistance color confining linear interaction which is suggested by the lattice QCD [35][36][37][38][39].
In the nonrelativistic limit, V eff (p, r) without tilde is transformed into the familiar nonrelativistic potential V eff (r) [2,40] V eff (r) = H conf + H hyp + H so (2) with where H conf is the spin-independent potential which contains a linear confining potential S (r) = br + c and the one-gluon exchange potential G(r) = −4α s (r)/3r, H hyp and H SO are the with such subjects. Here, the modified GI model is adopted to calculate the mass spectrum of pseudoscalar mesons.
In Eq. (5), the spin-orbit interaction can be divided into two types in which H so(cm) is the color-magnetic term and H so(tp) is the Thomas-precession term. Their expression can be written as where L is the orbital momentum between quark and antiquark. Noting that the above interaction potentials are obtained in the nonrelativistic limit, and they can optimized by introducing the phenomenological relativistic effects. In the GI model, the relativistic effects are imposed into the model mainly by two ways. Firstly, a smearing function ρ(r − r ) is introduced to incorporate the effects of an internal motion inside a meson and nonlocality of interactions between quark and antiquark. A smearing transformation is given bỹ with where σ 0 =1.80 GeV and s=1.55, are the universal parameters in the GI model, f (r) is a arbitrary function and notation tilde stands for that the expression has been performed smearing transformation. By smearing transformation, the one-gluon exchange potential G(r) = −4α s (r)/(3r) and linear confined potential S (r) = br + c are changed as (20) where E 1 = m 2 u/d + p 2 and E 2 = m 2 s + p 2 are the energies of the quark and antiquark in the meson, and m 1 = m u/d , m 2 = m s , i is parameter for a different type of hyperfine and spin-orbit interactions, which include the contact, tensor, vector spin-orbit and scalar spin-orbit potentials. Here, vector spin-orbit and scalar spin-orbit potentials correspond to the Eq. (5) related to one gluon exchange and confinement term, respectively. So the total Hamiltonian can be written aŝ whereH wherẽ For solving Schrödinger equationĤΨ = EΨ withĤ shown in Eq. (21), the simple harmonic oscillators (SHO) wave function will be employed. In the configuration space, SHO wave function has the form with (28) where Y LM L (Ω) is spherical harmonic function with orbital angular momentum quantum number L, and L L+1/2 n−1 (x) is an associated Laguerre polynomial, and β is a parameter of oscillator radial wave function. A series of SHO wave function with different radial quantum number n can be regarded as a complete basis to expand the exact radial wave function of meson state, in this case, the meson mass spectrum can be obtained by diagonalizing the Hamiltonian matrix of Eq. (21) based the above SHO basis. The total wave function of meson is composed by color, flavor, spin, space wave function, and the spin wave functions χ are The space-spin wave function R nL (r, β)φ LS JM with total angular quantum number J can be constructed by coupling L ⊗ S and has form After calculating each matrix element, the mass and wave function of meson could be obtained and they also are available to the following strong decay process. Although the GI model has achieved great success in describing the meson spectrum, there still exists a discrepancy between the predictions given by the GI model and recent experimental observations. The previous work [5] presents a modified GI model with a screening potential whose predictions can be well consistent with the experiment data for the charm-strange mesons. For higher excitation states, the authors of Ref. [5] believe that a screening effect plays a very important role which could be introduced by the transformation br + c → b(1−e −µr ) µ + c, where µ is a screening parameter whose particular value is need to be fixed by the comparisons between theory and experiment. Modified confinement potential also need to make similar relativistic correction which has been mentioned in the GI model. Then, we further write V scr (r) as the way given in Eq. (13), By inserting the form of ρ(r − r ) in Eq. (11) into the above expression and finishing this integration, the concrete expression forṼ scr (r) is given bỹ It is worth mentioning that after the confinement potential is replaced with a screening potential, other treatments are similar to the original GI model including the calculation of matrix elements of the Hamiltonian.
For a decay process A → B + C, we can write where P B(C) is a three-momentum of a meson B(C) in the rest frame of a meson A. A superscript M J i (i = A, B, C) denotes an orbital magnetic momentum. The transition operator T is introduced to describe a quark-antiquark pair creation from vacuum, which has the quantum number J PC = 0 ++ , i.e., T can be expressed as which is constructed in a completely phenomenological way to reflect the creation of a quark-antiquark pair from vacuum, where the quark and antiquark are denoted by indices 3 and 4, respectively. A dimensionless parameter γ depicts the strength of the creation of qq from vacuum, where the concrete values of the parameter R which will be discussed in the later section. Y m (p) = |p| Y m (p) are the solid harmonics. χ, φ, and ω denote the spin, flavor, and color wave functions respectively, which can be treated separately. Subindices i and j denote the color of a qq pair. By the Jacob-Wick formula [50], the decay amplitude is expressed as and the general decay width reads where m A is the mass of an initial state A. In our calculation, we need the spatial wave functions of the discussed kaons and iso-scalar and iso-vector light mesons. which can be numerically obtained by the modified GI model.

III. MASS SPECTRUM ANALYSIS
Although the GI model has succeeded in describing the ground states of the kaon family, it does not well describe the excited states. Since unquenched effects are important for a heavy-light system, it is better to adopt the modified GI model (MGI) [5,6] which uses a screening potential with a new parameter µ. The parameter µ describes inverse of the size of screening. To use the MGI model to calculate the kaon family spectra, it is better to determine the value of a new parameter µ considering two features: The first is when we use the same parameter set as in Ref. [2] and add a new parameter µ, the mass of the ground state of the kaon family will be lower than the experiments. The second one is the value of µ may be not so small like the one in Refs. [5,6]. In fact, in bottomonium and charmonium states, Refs. [51,52] give the µ value about 0.1 GeV which is larger than the one in Refs. [5,6]. Since we do not know the real value of µ in the kaon family beforehand, we need to adjust the parameters by fitting with the experiments data. At first, the quark masses should be the same for all meson families. Secondly, we do not adjust the values of Λ and α s for the same reason. Since σ 0 and s are universal parameters which are resolved by the QQ system, we do not vary them in our fit. The confining term br + c will be replaced by the screening potential, so their parameters should be fitted again. The relativistic effects should be adapted to a different system with the different quark masses. So we fix the following seven parameters listed in Table I by fitting eleven experimental data which is listed in Table II.   In Table II, we select eleven experimental data of kaons listed in PDG and optimize these kaon masses to determine seven parameters in Table I. This optimization has χ 2 /n = 12.6 which is smaller than 90.2 for the GI model as shown in Table II. Another reasons why we choose these kaons to fix the parameters in our model is that there does not exist mixture between n 1 L L and n 3 L L states for these kaons. In order to obtain the optimum values of parameters and global and good fit of eleven data, we set "Error in fitting" in Table II so that the first two experimental data, corresponding to K and K * , have artificial larger error values instead of the real errors in the brackets in the fourth column. The results listed in Table II show that the MGI model is better than the GI model since the value of χ 2 /n of the MGI model is about 7 times smaller than that of the GI model and hence it is safely applied to describe the masses of the selected eleven kaons.
Although the MGI model is better than the GI model to depict eleven experimental data, we need to indicate that there may exist ∼ O(100 MeV) deviation between experimental and fitting results for several kaons, which is shown in Table II. Such a difference of experimental and theoretical results may be due to precision of experiment. For example, there is only one experiment [1] for K * 5 (2380) and K * 0 (1950). The confirmation to K * 5 (2380) and K * 0 (1950) is still absent. Thus, further experimental measurement of the resonance parameters of these kaons will be helpful to clarify this difference of experimental and theoretical results.
By using the parameters shown in Table I as input, we further calculate the masses of other kaons, which are collected in Table III, where we do not consider the mixing of states with n 1 L L and n 3 L L . Usually, there exists mixture of the n 1 L L and n 3 L L states, i.e., [53], where |nL and |nL are two mixing physical states and θ nL is the corresponding mixing angle. Introducing such mixing states, we find two mass relations m(nL) < m(n 1 L L ) and m(nL ) > m(n 3 L L ), which can be applied to identify these observed kaons with the same J P quantum number. Thus, we need to combine the mass relations and mass spectrum of kaons listed in Table III with the experimental data to further shed light on the properties of other observed kaons. We conclude that have J P = 1 + quantum number. K 1 (1270) and K 1 (1400) are the mixture of 1 1 P 1 and 1 3 P 1 states, i.e., K 1 (1270) and K 1 (1400) correspond to 1P and 1P states, respectively.
3. K 1 (1650) has J P = 1 + and M = (1650 ± 50) MeV [1]. Since the mass of K 1 (1650) is smaller than that of the 4. We suggest that K(1830) is a 3 1 S 0 state. Later, we will test this assignment by studying its decay behavior.
The possible assignment of K 3 (2320) is the 2F state, which is the mixture of states K(2 1 F 3 ) and K(2 3 F 3 ). As the partner of K 3 (2320), K(2F ) is till absent in experiments. In addition, we should mention that the 1F and 1F in the kaon family are still missing.
Surely, the above conclusions of possible quantum states are only from the point of mass spectra view. If we want to clearly study particle properties further, we also need to investigate decay behaviors, especially strong decays, and detailed study will be given in the next section.

IV. OZI-ALLOWED TWO-BODY STRONG DECAYS
In the previous section, calculating the spectra of the kaon family, we obtain kaon wave functions, too, at the same time, which can be used in the QPC model to study the strong decay of the kaon family. The parameter γ in the QPC model is determined by fitting with the experiment data [47]. Thus, there is no free parameter in the QPC model. We obtain γ =10.5 as shown in Table IV.

Channels
Experimental data Numerical result In the following, we mainly focus on the OZI-allowed twobody strong decay behaviors of these discussed kaons, by which we not only test these possible assignments to the observed kaons, but also provide more abundant predictions of higher radial and orbital excitations in the kaon family.

A. S -wave kaons
Since K(498) and K * (892) were established to be the 1 1 S 0 and 1 3 S 1 states in the kaon family, respectively, in this work we do not discuss them, but present the phenomenological analysis of the 2S and 3S states.

2S states
As the candidate of the 2 1 S 0 state, the K(1460) was listed in PDG. If further checking the experimental data, we find that the K(1460) was only reported in Refs. [55,56]. However, in the past thirty years, further experiment of the K(1460) was missing, which is the reason why the K(1460) was removed from the its summary table of PDG.
In Table V, we give the information on the partial and total decay widths of the K(1460) as an 2 1 S 0 state, in which one can find the K(1460) mainly decays into K * π, Kρ, and Kω. Here, our results are larger than the experimental data for K * π, Kρ, and the total width. If establishing the K(1460) to be an 2 1 S 0 state, we need to clarify these different between our calculation and experimental data. We expect an independent experiment to confirm the observation of the K(1460). Especially, we suggest precise measurement of the resonance parameters and partial decay widths of the K(1460).
In PDG, the K * (1410) is possible candidate of the 2 3 S 1 state. However, we must face the following puzzling facts: (1) the mass of the K * (1410) is smaller than that of the K(1460). Usually, an 2 3 S 1 state has mass higher than that of an 2 1 S 0 state. In addition, we also notice the theoretical results of the mass of an 2 3 S 1 state, i.e., papers of Refs. [2,4] and this work give the mass of an 2 3 S 1 state to be 1579, 1675, 1548, and 1580 MeV, respectively, all of which larger than the experimental data, 1414 MeV, if K * (1410) is an 2 3 S 1 state. Thus, we need to understand why there exists such puzzling mass relation for the K * (1410) and K(1460). (2) If the K * (1410) together with the ρ(1450), ω(1420), and φ(1680) forms an 2 3 S 1 nonet, one can notice that the mass of the K * (1410) as an 2 3 S 1 state is a bit small which was also indicated in Ref. [3].
In Table V, the obtained partial and total decay widths of the K * (1410) as an 2 3 S 1 are given, where we also compared our result with the experimental data. The main decay modes of the K * (1410) include the Kπ, K * π, Kρ, and Kη channels. The obtained total decay width of the K * (1410) is consistent with the experiment result. We also notice that the ratio Γ Kπ /Γ T otal obtained in this work is a little bit larger than the experimental value (Γ Kπ /Γ T otal = (6.6 ± 1 ± 0.8)%). The above result is gotten by assuming the K * (1410) as an 2 3 S 1 pure state. In fact, the K * (1410) could be a mixture of 2 3 S 1 and 1 3 D 1 states. Thus, in the following, we further discuss such an S-D mixing effect on the ratio Γ Kπ /Γ T otal of the K * (1410). The K * (1410) and K * (1680) as the mixture of 2 3 S 1 and 1 3 D 1 states can be expressed as where θ sd denotes the mixing angle. Under this scenario, we present the decay behavior of the K * (1410) dependent on θ sd as shown in Fig. 1. The result shows that the experimental total width [57] of the K * (1410) can be described when θ sd is taken as ∼ 90 • or ∼ −90 • , which supports the K * (1410) as a pure 2 3 S 1 state. We need to emphasize that the branching ratio Γ Kπ /Γ T otal becomes larger when |θ sd | becomes smaller. Thus, the S-D mixing effect on K * (1410) state is not obvious if describing the experimental data. Of course, we must admit that there still exists a small difference between theoretical and experimental results for the ratio Γ Kπ /Γ T otal . Finally, we give a conclusion for the K * (1410). The mass of K * (1410) as a 2 3 S 1 state is relatively small, and there exists some disagreement in the branching ratios with experiments. Obviously, confirmation of this state assignment needs more experimental information and theoretical study.

3S states
Although the K(1830) is not listed in the summary table of PDG, we still select the K(1830) as a possible candidate of the 3 1 S 0 state and study its decay behavior.
In Table V, the partial and total decay widths of the K(1830) as a 3 1 S 0 state are shown. Our results show that the largest decay width of K(1830) is given by the channel K * (1410)π instead of K * ρ given by Ref. [3]. The other main decay channels contain K * π, K * 0 (1430)π, Kρ, K * 2 (1430)π, and K * ρ. The total width of the theory agrees with the experimental data. Our prediction of the decay information on this state will be helpful for the future experimental study, since there exists only two experimental studies on K(1830) until now.
At present, the 3 3 S 1 state in the kaon family is still absent. Thus, in this work we predict its decay property, where we take the predicted mass of the 3 3 S 1 state by the MGI model as the input. The results shown in Table V indicate that its important decay modes are πK * (1410), K * π, Kρ, and K * 2 π, Kπ, πK(1460), Kη(1295). Additionally, Ka 2 also have sizable contribution to the total width. This predicted decay information is useful to the future experimental search for this missing In Table VI, we show the allowed decay channels of the K * 0 (1430), and the corresponding partial and total decay widths. Here, its dominant decay channel of the K * 0 (1430) is Kπ, which has decay width 314 MeV, which is comparable with the experiment data (267 ± 36) MeV listed in PDG [1]. Besides, the Kη decay channel also has sizable contribution to the total decay width of K * 0 (1430). In addition, the obtained total decay width is consistent with the experimental measurement just shown in Table VI. The above study indicates that the 1 3 P 0 assignment to the K * 0 (1430) is suitable. The K * 2 (1430) together with a 2 (1320), f 2 (1270), and f 2 (1525) may form a 1 3 P 2 nonet. In Table VI, we give the partial decay widths of the K * 2 (1430). It dominantly decays into Kπ and K * π, while the Kρ, Kω, and Kη modes also have sizable contributions in which Kη was already observed in experiment [1]. According to Table VI, we can find that our results are consistent with experimental data. Thus, the K * 2 (1430) as a 1 3 P 2 state in the kaon family can be supported by our study of its decays. The K 1 (1270) and K 1 (1400) as the 1P and 1P states re-spectively satisfy where θ 1P denotes the mixing angle, which makes us discuss the θ 1P dependence of the partial and total decay widths of the K 1 (1270) and K 1 (1400). According to Fig. 2 which describes mixing angle θ 1P dependence of the K 1 (1270) decay width, we find that θ 1P should be taken as either 22.5 • ∼ 29 • or 41.5 • ∼ 48 • by fitting the CNTR data of Γ K * π [55], which is fortunately in the same range when fitted with the ratio Γ(K * π) D−wave /Γ(K * π) S −wave = 1 ± 0.7 [55]. Here, the central value of this mixing angle is We further investigate the decays of the K 1 (1400). CNTR [55] also gave Γ K * π /Γ T otal = (94 ± 6)% [60] for the K 1 (1400), by which we obtain 38 • < θ 1P < 68 • with central value θ 1P = 45 • , where the details can be found in Fig.  3 2 Hence, the above analysis shows that the mixing angle θ 1P favors 45 • which agrees with the conclusion made in Refs. [3,47] but disagrees with Refs. [61][62][63], in which they obtained θ 1P = 34 • and ∼ 60 • , respectively.

2P states
As shown in Table III, papers of Refs. [2][3][4] and this work give the mass of the 2 3 P 0 state 1.890, 1.791, 1.850, and 1.829 GeV, respectively, which are all smaller than the experimental value 1945 MeV if the K * 0 (1950) is assumed to be a 2 3 P 0 state. Under the assignment of the 2 3 P 0 state to the K * 0 (1950), we study the strong decay behavior of the K * 0 (1950), which is presented in Table VI. Our results show that the K * ρ mode is its dominant decay channel. Its total decay width can reach up to 1000 MeV which is 5 times larger than the experimental value 200 MeV. We also notice the result of the total decay width of a 2 3 P 0 state given by Ref. [3], which is two times larger than the experimental value, where they use a smaller phase space (their mass of a 2 3 P 0 state is 1850 MeV). We also obtain the branching ratio Γ Kπ /Γ T otal = 6.4%, which is close to 10.5% calculated by Ref. [3], but are smaller than the experimental value 52%. Besides, we also confirm that K 1 (1270)π has sizable contribution to the width of the K * 0 (1950) [3]. It is obvious that there exists a difference between the present theoretical and experimental results. Until now, the K * 0 (1950) has not been established in experiment since this state was omitted from the summary table of PDG [1]. For clarifying it, we suggest further experimental study of the K * 0 (1950), where its resonance parameter and partial decay widths are crucial information.
Then, we discuss possibility of two different assignments to the K * 2 (1980) from two aspects, mass and decay information. In 1987, LASS reported a structure in the reaction K − p → K 0 π + π − n [64], and the obtained resonance parameter are M = (1973 ± 8 ± 25) MeV and Γ = (373 ± 33 ± 60) MeV. This is the particle called K * 2 (1980) listed in PDG [1]. Barnes et al. [3] have the viewpoint that the K * 2 (1980) is a 1 3 F 2 state, and give the total width 300 MeV. However, our results show that the mass of a 1 3 F 2 state is about 2093 MeV. Thus, the mass of the K * 2 (1980) is a bit small if K * 2 (1980) is a 1 3 F 2 state, which can be supported by another fact, i.e., as an iso-vector 1 3 F 2 state, a 2 (2030) is well established in Ref. [32]. In the same 1 3 F 2 nonet, the meson which contains one s quark is heavier than the mesons which only contain u/d quarks. Along this line, the mass of the 1 3 F 2 state in the kaon family should be heavier than 2030 MeV.
Assuming the 1 3 F 2 state assignment to K * 2 (1980), we illustrate its decay behavior. The present work (see Table VII) shows that the K 1 (1270)π is the dominant decay channel when we treat the K * 2 (1980) as a 1 3 F 2 state, even though the channel is not observed in experiments. K 2 (1770)π, Kb 1 , Ka 1 , Kπ, Kρ, and K * π modes, among which Kρ and K * π have been reported in the experiment [1], also have sizable contributions, where we take θ 1D = −39 • . Our prediction for the channels K 1 (1270)π, K 2 (1770)π, Kb 1 , Ka 1 , and Kπ will be helpful for the experimental test of this assignment.
Besides the assignment of the 1 3 F 2 state to the K * 2 (1980), there exists another possibility, the K * 2 (1980) as a 2 3 P 2 state. The analysis of mass spectra in Refs. [2][3][4] and this work shows that the mass of a 2 3 P 2 state is 1938, 1896, 1850, and 1870 MeV, respectively. Thus, the experimental mass value of K * 2 (1980) is a bit larger as a 2 3 P 2 state. If the K * 2 (1980) is a 2 3 P 2 state, its main decay modes are K * ρ, Kπ, K * π, Kρ, K * η, Kη , and K * ω. Besides the Kρ and K * π modes, one notices that K f 2 has been observed in experiments which has a sizable contribution in theory. Hence, the K * 2 (1980) as a 2 3 P 2 state is also a possible assignment.
Just presented above, we discuss two assignments to the K * 2 (1980), where the decay behaviors of the K * 2 (1980) under two assignments are different. Thus, we should combine further experimental decay information of the K * 2 (1980) with our results to determine which possibility of its assignments we should take.
The K 1 (1650) and its partner K 1 (2P ) satisfy In Fig. 4, we show the partial and total decay widths of the K 1 (1650) depending on the mixing angle θ 2P if the mass of the K 1 (1650) is adored to be M = 1650 ± 50 MeV [1]. Since the decays of K 1 (1650) into Kππ and Kφ were observed in experiments [55,[65][66][67], we can roughly conclude that θ 2P is probably less than zero as seen from Fig. 4, where Kρ, Kφ and K * π have sizable contributions to the total decay width of the K 1 (1650) in our calculation.
In experiment, the K 1 (1650) is also not well established since this state is omitted from the summary table of PDG [1]. More experimental and theoretical efforts are necessary to establish the K 1 (1650). We notice new experimental information of the K 1 (1650) from LHCb [66], where the measured mass of the K 1 (1650) is (1793 ± 59 +153 −101 ) MeV which is about 150 MeV larger than experimental data given by Ref. [1]. Taking the LHCb mass result as an input, we investigate the strong decay behaviors of the K 1 (1650) again, which are shown in Fig. 5. Here, K * π, K * ρ, Kρ and K * 2 (1430)π are dominant decay channels. However, we cannot give further constraint on the mixing angle θ 2P by comparing an experimental width with our theoretical result due to a large experimental error of the LHCb experimental data.
In the following, we discuss the partner of the K 1 (1650). Ref. [68] gives the following equation about the mass relation between the pure states and physical states Substituting m K 1 (2 1 P 1 ) = 1840 MeV and m K 1 (2 3 P 1 ) = 1861 MeV given in Table III (1980), where the values in brackets and without brackets in the third and fourth columns represent those for the K * 2 (1980) as the 1 3 F 2 and 2 3 P 2 states, respectively. The unit of the width is MeV.
C. D-wave kaons 1. 1D states K * (1680) together with ρ(1700) and ω(1650) forms a 1 3 D 1 nonet. Barnes et al. [3] predicted that this state should have the mass 1850 MeV, but we obtain 1.766 GeV which is closer to the experimental value 1.717 GeV. The mass spectrum anal-ysis supports K * (1680) as a 1 3 D 1 state.
As shown in Table IX, K * (1680) as a pure 1 3 D 1 state mainly decays into final states K 1 (1270)π, Kh 1 , Kπ, and Kη, while the K * π and Kρ modes also have sizable contributions. We notice that the obtained ratios of partial decay widths of Kπ, K * π, and Kρ in this work are comparable with experimental data given in PDG. Since the branching ratios of the Kπ, K * π, and Kρ decay channels given by PDG are 38.7%, 29.9%, and 31.4%, we conclude that the remaining K 1 (1270)π decay channel has a very small width. However, our calculation shows that K 1 (1270)π is a main contribution to the total width which is consistent with conclusion from the former analysis in [3] but contradicts with the present experimental data.
Here we and the authors of Ref. [3] adopted the mixing angle θ 1P = 45 • [4,47] in the corresponding calculations. It is obvious that we need to face this puzzle in this channel. More experimental and theoretical efforts are needed to clarify this point. The main strong decay widths of the K 1 (2P ) state which as the partner of the K 1 (1650). Here, the, the mass of the K 1 (2P ) is taken as 1968 MeV [66]. c = cos θ 2P and s = sin θ 2P . The unit of the width is MeV.
K * 3 (1780) together with ρ 3 (1690), ω 3 (1670), and φ 3 (1850) forms a 1 3 D 3 nonet. We give its mass 1.781 GeV by the MGI model, which is consistent with experiment data 1.776 GeV. As shown in Table IX, even though K * ρ is the dominant decay mode of K * 3 (1780), it is not observed in experiments so far. The channel K * ω has sizable contribution to its total decay width, which is still missing in experiment. A final state Kπ largely contributes to the total width and theory and experiments are consistent to each other. The branching ratio Γ K * π /Γ Kπ agrees with experimental data [70]. Next, let us focus on K 2 (1820). According to Fig 7, one notices that K 2 (1820) probably decays to Kππ, K * 2 (1430)π, K f 2 (1270), K * π, and Kω, in which Kππ comes from Kρ channel. As seen from θ 1D dependence of the widths of K 2 (1820) in Fig. 7, we notice that contributions of Kρ and K f 2 (1270) are large when θ 1D < 0, which indicates that it is very likely that θ 1D is smaller than zero that is consistent with the previous analysis for K 2 (1770). Because of absence of the experimental information, we cannot confirm the angle θ 1D , while our results will be helpful for the future experiments to study this state.

2D states
As one of the 2D states is missing, using the familiar program with 2P states, we obtain where K 2 (1990) is obtained from an equation similar to Eq. (42).
According to Table X, one finds that K 2 (2250) as a 2D state mainly decays into Ka 2 (1700), K * 2 (1980)π, K * 3 (1780)π, K * (1410)π, and K * 2 π. K * 2 π and K f 2 (1270) have been observed in experiments which have sizable contributions to the total width. Besides, Kρ, which is an important decay channel in our result, can decay into Kππ that is observed in experiment. On the other hand, the theoretical total width is larger than the experimental value 180 MeV given in PDG. We need more experimental information to study this 2D state, to test our results, and to have more detailed decay widths to ascertain the value of θ 2D .
We use 1994 MeV as the mass of the partner of K 2 (2250) with θ 2D ≈ ±44 • , and calculate the strong decay of this state as shown in Table XI. According to this Table, we can find that its main decay channels are K * π, Kρ, K * 2 π, and K * 3 (1780)π.
D. F-wave kaons

1F states
In this subsection, we discuss possibility of different assignments of K * 2 (1980) from two aspects, mass and decay information. In 1987, LASS reported a structure in the reaction K − p →K 0 π + π − n [64], and they obtained resonance parameters M = (1973 ± 8 ± 25) MeV and Γ = (373 ± 33 ± 60) MeV. This is the particle called K * 2 (1980) listed in PDG [1]. Barnes et al. [3] have the viewpoint that K * 2 (1980) is a 1 3 F 2 state, and give the total width 300 MeV. On the other hand, our results show that the mass of a 1 3 F 2 state is about 2093 MeV. Ebert et al. [4] predict 1 3 F 2 state with the mass 1964 MeV. As the partner of an iso-vector of 1 3 F 2 , a 2 (2030) is well established in Ref. [32]. In the same nonet, the meson which contains one s quark is probably heavier than the mesons which only contain u/d quarks. Along this line, the mass of 1 3 F 2 state of the kaon should be lager than 2030 MeV, so that the mass K * 2 (1980) is a bit small as a 1 3 F 2 state. The papers of Refs. [2][3][4] and this work give the mass for a 2 3 P 2 state 1938, 1896, 1850, and 1870 MeV, respectively, and for its iso-vector partner is a 2 (1700) [32], the mass of K * 2 (1980) is a bit larger as a 2 3 P 2 state. We should, of course, combine the decay information of K * 2 (1980) to determine which possibility of its assignment we should take.
Both Ref. [3] and our results show that K * 2 (1980) is 2 3 P 2 state, where Ref. [3] takes the mass 1850 MeV and we take the experimental value 1973 MeV. For this reason, the results between theirs and this work have some difference. The main decay modes are K * ρ, Kπ, K * π, Kρ, K * η, Kη , and K * ω. Besides the Kρ and K * π, one notices that K f 2 has been observed in experiments which has sizable contributions in theory. Hence K * 2 (1980) assigned to 2 3 P 2 is also reasonable. Finally, let us draw a rough conclusion for K * 2 (1980). According to the mass analysis, the mass of K * 2 (1980) is a bit small when assigned to a 1 3 F 2 state and a bit large when assigned to a 2 3 P 2 state. According to the decay information, K * 2 (1980) is in favor of a 2 3 P 2 state. We still, however, need more experimental information to test our assignment for K * 2 (1980). What is more important is that we give the prediction that the partial widths of K * 2 (1425)π, K * (1410)π, K * ω, K * ρ, and K * η treating K * 2 (1980) as 1 3 F 2 will be much larger than those of the case of a 2 3 P 2 state. The experimental study of these decay modes combined with our prediction will help us determine assignment of K * 2 (1980). Besides the above, our prediction can help the future experiments find the missing1 3 F 2 or 2 3 P 2 state. According to Table III, one can notice that our spectral results are consistent with the mass of K * 4 (2045) given by PDG when we treat it as 1 3 F 4 state. As for the strong decay of K * 4 (2045), one can notice that both the results of Ref. [3] and this work (Table XII) show that K * ρ, Kπ, K * π, K * ω, and Kρ are the main decay channels. The PDG gives two partial width ratios: one is Γ Kπ /Γ T otal = (9.9±1.2)% and our result is Γ Kπ /Γ T otal = 8.4% which is consistent with the experiment. Another is Γ K * φ /Γ T otal = (1.4 ± 0.7)% and we obtain Γ K * φ /Γ T otal = 1.54%, which is consistent with the experiment as well. On the other hand, Ref. [3] obtained 21% and 3.1% for these two partial width ratios, which are different from the experiment. these results, of course, prove the superiority of the accurate meson wave functions we have ob-  Papers of Refs. [2][3][4] and this work give the mass of a 1F state 2131, 2009, 2050 , and 2075 MeV (which we call K 3 (2075), and strictly speaking, this state is a pure 1 1 F 3 state, here we assume the physical state 1F has this mass), respectively, among which one notices that the last two results are almost identical. K 3 (2075) is assigned to the missing 1F state. We present the θ nF dependence of widths for these two cases in Tables XIII. The total width of a 1F state with the mass 2075 MeV is about (400 − 600) MeV, which means that the predicted K 3 (2075) is a broad state and it is not easy to identify K 3 (2075) in experiments. Its main decay channels are K * 3 (1780)π, K * ρ, K * π, Ka 2 , Kρ and K * 2 π.
E. G-wave kaons 1. 1G states K 4 (1G) and K 4 (1G ) mixing satisfies We assume the mass of a 1G state is about 2309 MeV, which we call K 4 (2310). The GI model [2] and Ref. [4] give this mass 2422 MeV and 2255 MeV, respectively, while the mass 3.1% 1.4 ± 0.7% [74] of the K 4 state in PDG is 2490 MeV. According to Table III, K 4 (2500) may be a 2G state. We predict the strong decay information of these two G wave states in Tables XV. As shown in Table XV, main decay modes of K 4 (2310) are K * 4 (2045)π, K * 3 (1780)π, Kρ 3 (1690), Ka 2 and K * 2 π when K 4 (2310) is assigned to a 1G state. Its total width will be (710 − 880) MeV, which is not easy to observe in experiments.

V. CONCLUSIONS AND DISCUSSION
In this paper, we have given the analysis of mass spectra of the kaon family via the modified Godfrey-Isgur quark model that includes a color screening effect, and have obtained the structure information of the observed kaon candidates. Then, we have further tested the possible assignments by comparing the theoretical results of their two-body strong decays with the experimental data. Additionally, we have also predicted the behaviors of some partial decay widths of the kaons, which are still missing in experiments. In Table XVII, we summarize the mass and main decay modes of these states, by which experiment may carry out the search for them.
This study is crucial to establish the kaon family and future search for their higher excitations. We have discussed the possible assignments to the kaons listed in PDG. The main task of the present work has been a calculation of the spectra and OZI-allowed two-body strong decays of the kaon family, which can test the possible assignments to the kaons. In Sections II and III, we have discussed these points in detail. The predicted decay behaviors of the discussed kaons can provide valuable information for further experimental study in the future.
At present, experimental information on the kaons is not abundant. Thus, we suggest to do more experimental measurements of the resonance parameters and to search for the missing main decay channels. Such an effort will be not only helpful to establish the kaon family in experiments but is also valuable to study the production of hidden-charm pentaquarks P c (4380) and P c (4450) by analyzing Λ b → J/ψpK [77] which has a close relation to the understanding the kaon family. With experimental progress, the exploration of the kaons will become a major issue in hadron physics, provided by good platforms in the BESIII, BelleII, and COMPASS experiments. We hope that, inspired by this work, more experimental and theoretical studies of high-spin states are conducted in the future.   The widths of (predicted) K 4 (2310) and K 4 (2500) depending on their mixing angle, where s and c represent sine and cosine