Decay of the pseudoscalar glueball into vector, axial-vector, scalar and pseudoscalar mesons

We resume the investigation of the ground-state pseudoscalar glueball, JPC=0-+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{PC}=0^{-+}$$\end{document}, by computing its two- and three-body decays into vector and axial-vector quark-antiquark meson fields additional to scalar and pseudoscalar mesons through the construction of an interaction chiral Lagrangian that produces these decays. We evaluate the branching ratio, via a parameter-free calculation, by setting the mass of the pseudoscalar glueball to 2.6 GeV as predicted by lattice QCD simulations. We duplicate the computation for the branching ratios for a pseudoscalar glueball mass 2.37 GeV which matches to the measured mass of the resonance X(2370) in the BESIII experiment. The present channels and states are potentially reached and are interesting for the running BESIII and Belle-II experiments and the planned PANDA experiment at FAIR/GSI which will be able to detect the pseudoscalar glueball within the accessible energy range.


I. INTRODUCTION
Glueballs, the composite particles containing gluons without valence quarks, are predicted by the confinement properties of Quantum chromodynamics (QCD) [1] and the non-Abelian nature [2] of the SU (3) c colour symmetry by virtue of the gauge fields of QCD, the gluons, self-interaction and strong vacuum fluctuations.A scalar glueball, J P C = 0 ++ , is the ground state and its mass range is estimated to be from 1000 to 1800 MeV.That followed by a pseudoscalar glueball, J P C = 0 −+ , at higher mass.Up to now, glueballs remain experimentally undiscovered [3] because of their mixing with ordinary meson states, and no meson is listed by the Particle Data Group (PDG) [4] to be unambiguous of predominant glueball nature.Therefore, the search for glueballs witnessed extensive and intensive investigations by both theoretical and experimentally studies [5][6][7][8][9][10][11][12][13][14].Actually, ten scalars including the isoscalars f 0 (1370), f 0 (1500), and f 0 (1710) are observed.The last resonance is classified to be predominantly a scalar glueball, according to the result of the numerical approach of the lattice QCD [15,16] and effective approaches [9,17,18].Theoretical studies of glueballs are performed by nonperturbative approaches, the flux tube model [19], constituent models [20], the holographic approach [21], effective chiral models [7,9,11,12] and lattice QCD simulations [15,16] which reach the spectrum of glueball states below 5 GeV, meaning that these simulations play an important role in the investigation of the low energy strong interaction phenomena.In the quenched approximation, lattice QCD simulation computed the masses of glueballs [5,15].For example, it predicted the mass of the ground-state pseudoscalar glueball around 2.6 GeV and the mass of the first excited state of the pseudoscalar glueball around 3.7 GeV.On the other hand, the production rates of glueballs in the J/ψ radiative decays [22,23] represent an additional role in determining the glueballs owing to the gluon-rich environment.In the quenched approximation, the pure gauge glueballs are well-defined hadron states.Therefore, the electromagnetic form factors of J/ψ radiatively decaying into glueballs are directly extracted from the calculation of the matrix elements of the electromagnetic current between glueballs and J/ψ.The BESIII collaboration studied the process J/ψ → γη π + π − and observed in the η π + π − decay the resonance X(2370) with quantum number J P C = 0 −+ [24].This resonance has the same quantum number of a pseudoscalar glueball and lies in the mass range of the pseudoscalar glueball in lattice QCD prediction.That leads us to suggest the assignment of the pseudoscalar glueball also to the X(2370) in our work previously published in Ref. [7] and references therein.There are also several candidates for the pseudoscalar glueball as seen for example in Ref. [25], where the lower pseudoscalar glueball state is suggested to be the upper iota η(1490).Moreover, the resonance I(1440) is required to be a pseudoscalar glueball in Ref. [26] out of the η − η − I(1440) − η C mixing investigation.The phenomenology of the pseudoscalar glueball is studied in a family of finite-width Gaussian sum rules upon a correction from instanton-gluon interference to the correlation function in Ref. [27], where the authors concluded that the interference contribution is gauge-invariant.
The present study of the decay properties of the ground state pseudoscalar glueball (denoted as G) is based on the chiral symmetric model of low-energy QCD called the extended Linear Sigma Model (eLSM) [28].It contains all quark-antiquark mesons with (pseudo)scalar and (axial)vector as well as a scalar and a pseudoscalar glueball and implements the symmetries of the QCD and their breaking.The eLSM is interesting for the study the hadron phenomenology.One can see in Ref. [29] that the eLSM compatible with chiral perturbation theory for what concerns low-energy pions (most notably, pion-pion scattering).In detail, Ref. [28] achieved a good description of PDG data by a fit to various experimental quantities.That fit allowed to fix the parameters of the model, subsequently making other predictions/postdiction.Both conventional qq-states and various non-conventional gluonic mesons were already studied in the eLSM.The phenomenology of the light mesons [28] and excited mesons [30] have been nicely described, as well as the properties of the open and hidden charmed mesons [31,32].As a consequence of dilatation invariance and its anomalous, the scalar glueball appears naturally in the eLSM which allowed to study the vacuum properties of the scalar glueball [9,33].The eLSM has been also applied to the hadronic decays of the pseudoscalar glueball(s) [7,11,12] and the vector glueball [34].As an additional advantage, in Ref. [35][36][37][38], the eLSM has been successfully applied in the baryonic sector within the so-called mirror assignment which the predictions turn out to be in agreement with the experimental data for pion-nucleon scattering and baryonic decays.The eLSM was also employed at a finite temperature [39,40] and density [38,41], to describe the chiral phase transition in the medium.On the other hand, the inclusion of hybrid mesons was presented in Ref. [42].
Within the eLSM [7], we have analysed the decay of the ground-state pseudoscalar glueball into scalar and pseudoscalar mesons and found that the channel KKπ is dominant and the KKη and KKη decay modes are sizable.Moreover, the hadronic decays of the pseudoscalar glueball into nucleons were computed [8] and into (pseudo)scalar mesons and their first excited state Ref. [12].Furthermore, the decay properties of the first excited pseudoscalar glueball have been also studied in Refs.[11,12].These efforts on the pseudoscalar glueball and its first excitation properties are important in the comprehension of the non-perturbative behaviour of QCD and useful in searching for the pseudoscalar glueball in future experiments.Based on that, we are interested in continuing our investigations on the properties of glueballs.
In this paper, we use both masses M G = 2.37 GeV relevant to the BESIII experiment candidate and M G = 2.6 GeV predicted by the lattice QCD simulation, to calculate the decay widths of the pseudoscalar glueball in the framework of the constructed effective model so to connect the pseudoscalar glueball, gg, to qq vector and axial-vector mesons in addition to scalar and pseudoscalar mesons.This work is a further step in our investigations of the pseudoscalar glueball and predicts new decay channels including vector and axial-vector mesons.We can thus compute the widths for the decays G → P S, G → P V , G → P P P , G → P P A, G → P P V , and G → P SV , where P, S, V , and A stand for pseudoscalar, scalar, vector, and axial-vector quark-antiquark states, respectively.The pseudoscalar field P corresponds to the well-known light mesons {π, K, η, η } and the scalar S corresponds to the scalars above 1 GeV: {a 0 (1450), K * 0 (1430), f 0 (1370), f 0 (1500)}, while the vector state V refers to {ρ(770), K * (892), ω, φ} and the axialvector A refers to {a 1 (1260), f 0 (1285), f 1 (1420), K 1 (1200)}.The results are presented as branching ratios in order to disregard the unknown coupling constant.The present results confirm all channels which were already predicted earlier in Ref. [7]: the decays of the pseudoscalar glueball into scalar and pseudoscalar mesons.This is particularly interesting, as it may help the community understand the hadron spectrum and the search for glueballs in future experiments.The present investigation for the two-and three-body decays of the ground-state pseudoscalar is a useful guideline for both the running BESIII/(Beijing, China) and Belle-II/(Tsukuba, Japan) experiments and the planned PANDA experiment at the FAIR/(GSI, Germany) [43] which are able to measure the proposed channels.
The present paper is organized as follows.In Sec.II we present the chiral multiplets.Then, in Sec.III we introduce the constructed effective Lagrangian which describes the two-and three-body decays of the ground-state pseudoscalar glueball into vector, axial-vector, scalar and pseudoscalar quark-antiquark degrees of freedom, allowing for the branching ratios prediction for the decays.Finally, in Sec.IV we present the conclusions.

II. CHIRAL MULTIPLETS
In this section, we present the quark-antiquark fields which represent the pseudoscalar glueball decay products.On the basic ingredients of the eLSM, the (pseudo)scalar and (axial-)vector field mesons are presented below, where the chiral combinations were properly taken into the account.
The multiplet of the scalar S a , and the pseudoscalar P a , quark-antiquark states is introduced [28] as where t a are the generators of the group U (N f ).The multiplet Φ transforms under are U (3) L(R) matrices, and under the charge conjugation C as Φ → Φ T as well as under the parity The vector V a and axial-vector A a , degree of freedom, are presented [28] in the following left-and right-handed matrices, L µ and R µ , as and Under In the present investigation, we are interested in studying the hadronic decays of the pseudoscalar glueball field G which is chirally invariant and transforms under charge conjugation as G → G and under the parity Consequently, these transformation properties of the multiplets Φ, L µ , R µ and the pseudoscalar glueball G have been used to construct the below effective invariant Lagrangian ( 6) and the extended Linear Sigma Model (eLSM), see Appendix A and Ref. [28] as well.
The identification of the quark-antiquark fields in the present model (6) with the physical resonances presented in details in Ref. [28], is straightforward in the light (pseudo)scalar and (axial-)vector states with mass 2 GeV.The Pseudoscalar sector P a includes the pion isotriplet − → π , the kaon isodoublet K [4], and the isoscalar fields η N ≡ ūu + dd / √ 2 and η S ≡ |ss , which represent the non-strange and strangeness mixing components of the physical states η and η [4] with mixing angle ϕ −44.6 • [28].
The scalar sector S a contains the isotriplet a 0 which refers the physical resonance a 0 (1450) and the kaon field K * 0 which is assigned to the physical isodoublet state K 0 (1430).In the scalar-isoscalar sector, the non-strange bare field σ N , the bare strange field σ S and the scalar glueball G mix and generate the three physical resonances f 0 (1370) , f 0 (1500), and f 0 (1710) through the following mixing matrix as constructed in Ref. [9]: However, the scalar-isoscalr fields, f 0 (1370), f 0 (1500), and f 0 (1700) are predominantly described by the bare configuration≡ ūu + dd / √ 2 state, ss states and a glueball gg state, respectively.We now turn to the assignment of the (axial-)vector states.The vector sector V a contains the iso-triplet field − → ρ , the kaon states − → K * , and the isoscalar states ω N and ω S which are assigned to the ρ(770), K * (892), ω and φ mesons, respectively [28].Note that the mixing between the strange and nonstrange isoscalars vanishes in the extended linear sigma model eLSM [28], whereas this mixing is so small as obtained in Ref. [44].In the end, for the axial-vector sector S a , the isotriplet − → a 1 , the isoscalar fields f 1N and f 1S correspond to the resonances a 1 (1260), f 1 (1285) and f 1 (1420) respectively.However, the four kaon states K 1 refer predominantly to the resonance K 1 (1200) and could also refer to K 1 (1400) because of the mixing between the pseudovector states and axial-vector states [45].

III. DECAY OF THE PSEUDOSCALAR GLUEBALL INTO CONVENTIONAL MESONS
We consider a chirally invariant Lagrangian which couples the ground-state pseudoscalar glueball G ≡ |gg with quantum numbers J P C = 0 −+ to the quark-antiquark vector, axial-vector, scalar and pseudoscalar field mesons which is invariant under U (3) R ×U (3) L , C, and P transformations.The coupling constant α is an unknown coupling constant and has a dimension of Energy −3 .Based on Ref. [46], we constructed this model where the heterochiral Φ and Φ † , involves (pseudoscalar) mesons, coupled to the homochiral L µ and R µ , consisting of (axial-)vector mesons, through only structures which contain derivatives of Φ.This Lagrangian describes the two-and three-body decays of the pseudoscalar glueball.Notice that the two-body decays appear only through the condensate and the interesting thing is that it does not lead to the two-body decays for the nonet of chiral partners A a .
According to the validity of the joint model the pseudoscalar glueball in the present work has a mass of about 2.37 GeV coinciding to a claimed BESIII experiment candidate and of 2.6 GeV from the lattice QCD prediction, while the eLSM is a low-energy chiral model valid up to 1.7 GeV.Therefore, it should be accepted that this model is suited to calculate exclusively the decays of the field G.
For instance in Refs.[31,32], the eLSM has found to be applicable to study the phenomenology of the heavy charmed mesons, which concerns the calculation of masses and large-N c dominant decays although it could be far from the natural domain of chiral symmetry.Moreover, the employed effective models in Ref. [34] have been used in the decay modes of the vector glueball with a mass of about 3.8 GeV.The present approach is consequently expected to be reliable within a similar accuracy, even if it is proposed here to test the decays of an (as of yet) unidentified glueball.From the Refs.[11,12] and refs.therein, we can prove the validity of the employed effective chiral model in Eq. ( 6) to study the decay modes of the pseudoscalar glueball.The novelty of our approach is that the qualitative outcomes do not depend on the precise value of the input parameters.
The lagrangian (6) shows that the pseudoscalar glueball state of such high mass could decay, which is not the only and first case, that appears in this field of investigations.There were widely used models which couple one heavy field to light mesons as seen for instance in Refs.[11,12,32,34,[47][48][49] and refs.therein.Also the decays of the heavy scalar and pseudoscalar charmonium states, χ c0 and η C , are studied by using the eLSM, that gave results in reasonable agreement with experimental data where available, see the details in Ref. [32], in additional to the accepted investigation of the decays of the heavy vector charmonium state J/ψ into light pseudoscalar mesons in Ref. [47].It is axiomatic that this basic assumption of effective hadronic models would be tested through only advanced lattice simulations and/or the future experimental discovering of glueballs.Definitely, when we couple the heavy field such as glueballs to the eLSM, we took into account that glueballs are flavour blind and chirally blind.So chiral symmetry, with its spontaneous breaking, does not affect the determination of the pseudoscalar hadronic decays.While up to now there are no data for a direct comparison, the decay ratios can only be predicted with model-dependent, and by neglecting a mixing influence and symmetry breaking terms.The present outcome branching ratios could be useful in the future search of the ground state pseudoscalar glueball.
For computing the decays of the pseudoscalar glueball, one has to perform the shift the scalar-isoscalar fields by their vacuum expectation values φ N and φ S to implement the effectiveness of the spontaneous symmetry breaking which takes place (m 2 0 < 0) [28] σ In matrix form: In addition, we shift also the axial-vector fields to redefine the wave-function renormalization constants of the pseudoscalar fields whereas i = 1, 2, 3, 4 indicates the four kaonic fields and j refers to N and S. The numerical values of all the parameters and the renormalization constants of the corresponding wave-functions appearing in the present paper expressions have been fixed in Ref. [28].Their values are summarized in Table I.The equivalence between the wave-function renormalization constants Z π and Z η N comes out from the isospin symmetry, accordingly for ω a1 and ω f 1N .The corresponding chiral condensates φ N and φ S read where the decay constant of the pion is f π = 0.0922 GeV and the kaon is f K = 0.110 GeV [4].
After performing the shift operations ( 8) and (10) in the Lagrangian ( 6), we obtain the relevant tree-level vertices for the decay processes of pseudoscalar glueball G.
The branching ratios of the pseudoscalar glueball G for the decays into two and three-body, P S, P V , P P P , P P V , P P A and P SV are reported in the following tables for two possible choices of the pseudoscalar glueball masses.The choice of the value 2.6 GeV is a consequence of the central value of a given lattice calculation and 2.37 GeV according to the obtained candidate by the BESIII experiment.The results are presented relative to the decay width of the pseudoscalar glueball into ππη, Γ G→ππη , to eliminate the unknown coupling constant.In Table II, we predict the same channels of the decay of the pseudoscalar glueball into P S, which were presented earlier in Ref. [7] additional to the only new channel G → KK * which describe the validity of the pseudoscalar glueball to decay into the vector meson K * and the pseudoscalar meson K.That confirm all the previous two-body P S channels and the decay of the pseudoscalar glueball into scalar-isoscalar states f 0 (1370), f 0 (1500) and f 0 (1700) including the full mixing pattern above 1 GeV.(For details of the two-body decay calculation, see Appendix B).The two-body decay channels a 0 π and KK * are sizable.Moreover, the two-body decay channel G → KK S can proceed through a sequential instance, K * 0 (1430) → Kπ, leading to the three-body decay G → KKπ.In order to obtain the total three-body decay width for G → KKπ, the two-and three-body decay amplitudes of this channel have to be added coherently before taking the modulus square as seen in Ref. [8].
Branching ratios for the decay of the pseudoscalar glueball G into a scalar, a pseudoscalar, a vector and an axial-vector meson.
In Table III, the results of the branching ratios of the pseudoscalar glueball G for three-body decays P P P are presented.Note that, these decay channels are the same decay channels that are predicted in our previous work Ref. [7], which kind of proves the validity of these decay channels.In Table IV, we turn out to list new processes for the branching ratios of the three-body decays of pseudoscalar glueball G into vector and axial vector mesons additional to scalar and pseudoscalar for both choices of M G = 2.6 GeV and M G = 2.37 GeV.That is important to widen our concept for the nature of the pseudoscalar glueball and also help researchers to detect the listed decay channels in experiments.The three-body decay channels ππη, KKπ and ππη are sizable, which were also sizable in Ref. [7].Comparing Table II and III, one sees constancies such as BR(KKπ) > BR(K * K) where the K * K is obviously a part of the KKπ channel.The decay channel Γ G→πππ is suppressed.The results depend only slightly on the glueball mass, which explains the similarity of their two columns.(For details of the three-body decay calculation, see Sec.C of the Appendix.)Fig.I shows the total decay action line of the pseudoscalar glueball, Γ tot G = Γ G→P S + Γ G→P V + Γ G→P P P + Γ G→P P A + Γ G→P P V + Γ G→P SV , as function of the coupling constant c for both masses suggested in the present work, where the decay into baryons is negligible.The coupling constant c has small value because of asymptotic freedom, and QCD decays are approximately those of free quarks and gluons at high energies.

IV. CONCLUSION
The two-and three-body decays of the ground state of the pseudoscalar glueball into a vector, an axial-vector, a scalar and a pseudoscalar quark-antiquark fields have been studied.We have started with the chiral invariant effective Lagrangian describing the interaction of the pseudoscalar glueball with (axial-)vector and (pseudo)scalar mesons for the three-flavour case N f = 3.The size of the coupling constant intensity can not be determined.That leads to predict the branching ratios for the decay channels which are expected to dominate.According to the mass of the pseudoscalar glueball, we considered two options: (i) M G = 2.6 GeV which is chosen to be in agreement with lattice QCD in the quenched approximation.(ii) M G = 2.37 GeV which assumes that the measured resonance X(2370) in the BESIII experiment is a pseudoscalar glueball (predominantly) [24].These masses can be tested in the planned PANDA experiment at FAIR/GSI [43], since glueballs could be directly formed in proton-antiproton fusion processes.The two-body decay of the pseudoscalar glueball produce P S (includes the scalar-isoscalar states f 0 (1370), f 0 (1500) and f 0 (1700)) and P V while the three-body decay produces P P P, P P A, P P V and P SV .The only new two-body decay channel is G → K * (892)K, see Table II for all results.The G → ππη channel is predicted in the dominant decay channel, followed by G → KKπ then G → KKη .On the contrary, the decay of the pseudoscalar glueball into πππ is predicted to vanish.From the investigation of the decay of the pseudoscalar glueball in the present and our previous work in Ref. [7], we confirm the validity of the two-and three-body decay processes, G → P S and G → P P P , that are presented in Table II and Table III, respectively.This indicates that experimentalists should search for glueballs through these channels.Moreover, the new three-body decay of the pseudoscalar glueball into (axial-)vector and (pseudo)scalar mesons is reported in Table IV.The present results of this work could be helpful for understand glueball spectroscopy and be used as guidelines in the search for the pseudoscalar glueball in the ongoing BESIII and Belle-II experiments and the future PANDA experiment at the FAIR/GSI.The discovery and study of the glueball spectroscopy is a crucial test of QCD.Therefore, we plan, in the framework of a chiral model, to study the decay of the first excited pseudoscalar glueball into vector and axial-vector mesons in addition to scalar and pseudoscalar mesons, and into nucleons as well.Forthcoming developments of the present work will be based on new results for the pseudoscalar glueball and its excitations when lattice QCD will include the effect of dynamic fermions and working beyond the quenched approximation.That would be very beneficial to understand quark and gluon confinement in QCD. and Appendix B: Two-body decay The general formula of two-body decay width is written as seen in Ref. [54]: with decaying particle A and the decay products B and C.Where is the center-of-mass momentum of the two particles production in the decay, M A→BC is the corresponding tree-level decay amplitude, and S A→BC refers to a symmetrization factor (it equals 1 if B and C are different and it equals 1/2 for two identical particles in the final state).
As an example of a two-body decay channel for G → P S, let us consider the case G → KK * 0 .This process is given from Eq. ( 6) as where f G→KK * 0 is the isospin factor, M G is the pseudoscalar glueball mass, and (B4) where m K and m K * 0 are the masses of the kaon and K * 0 mesons, respectively, while k GKK * 0 is the center of mass momentum of kaon and K * 0 and reads Next we turn to the decay of G → KK * as an example for G → P V , which its process is obtained from Eq. ( 6) as where f G→KK * = 4 and k GKK * is the center of mass momentum of K and K * .In an analogous way, all the decay processes G → P V in table I are calculated the corresponding change of the isospin factors, the decay products mass and the constants entering in the amplitudes.The quantities m 1 , m 2 , m 3 are the masses of the three decay products P 1 , P 2 , and P 3 , M G→P1P2P3 is the corresponding tree-level decay amplitude, and f G→P1P2P3 is a symmetrization factor which equals 1 if all decay products are different, equals 2 for two identical decay products in the final state, and it equals 6 if P 1 , P 2 , and P 3 , are identical in the final state.
For example, the amplitude for the process G → K * 0 K * 0 0 π 0 is The other three-body decays are calculated in analogous way.
Notice that there are several decay channels of the pseudoscalar glueball, G that appear in Eq. ( 6) but they are not kinematically allowed because the mass of the decaying particle is lower than the sum of the mass of the decay products M <

FIG. 1 :
FIG.1: Solid (orange) line: Total decay width of the pseudoscalar glueball with the bare mass M G= 2.6 GeV as function of the coupling c.Dashed (blue) line: the same curve for M G= 2.37 GeV.

TABLE I :
Parameters and wave-function renormalization constants.

TABLE II :
Branching ratios for the decay of the pseudoscalar glueball G into P V and P S.

TABLE III :
Branching ratios for the decay of the pseudoscalar glueball G into three pseudoscalar mesons.