The mass spectra and decay properties of dimesonic states, using the Hellmann potential

Mass spectra of the dimesonic (meson–antimeson) molecular states are computed using the Hellmann potential in a variational approach, which consists of a relativistic correction to the kinetic energy term as well as to the potential energy term. For the study of a molecular bound state system, the Hellmann potential, of the form V(r)=-αsr+Be-Crr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V(r)=-\frac{\alpha _{s}}{r} + \frac{B e^{-Cr}}{r}$$\end{document}, is being used. The one pion exchange potential is also incorporated in the mass calculation. The digamma decay width and decay width of the dimesonic system are evaluated using the wave function. The experimental states such as f0(980)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{0}(980)$$\end{document}, b1(1235)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_{1}(1235)$$\end{document}, h1(1380)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_{1}(1380)$$\end{document}, a0(1450)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{0}(1450)$$\end{document}, f0(1500)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{0}(1500)$$\end{document}, f2′(1525)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{2}'(1525)$$\end{document}, f2(1565)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{2}(15 65)$$\end{document}, h1(1595)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_{1}(1595)$$\end{document}, a2(1700)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{2}(1700)$$\end{document}, f0(1710)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{0}(1710)$$\end{document} and f2(1810)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{2}(1810)$$\end{document} are compared with dimesonic states. Many of these states (masses and their decay properties) are close to our theoretical predictions.


Introduction
The spectroscopy of hadrons imparts information of a basic force of nature. The theory of Quantum Chromo Dynamics (QCD) is a tool to understand the strong force in which hadron spectroscopy plays a key role in the non-perturbative and the perturbative regime. The new experimental developments at Belle, BES, CLEO, CDF, LHC, BABAR brought out enormous data and came with large numbers of surprises [1][2][3][4][5]. A number of new states are observed which do not fit in the conventional qq scheme, mainly in the meson sector. All these exotic states which do not fit in qqq and qq schemes require extra theoretical attention [6][7][8][9][10][11][12][13].
Meson spectroscopy, mainly light meson spectroscopy, gives information of the non-perturbative regime of QCD. The basic difficulty to study the light meson spectrum is that mostly resonances do not come out as narrow, isolated peaks [6,7]. Moreover, their large decay widths make them difficult to identify experimentally as well as to differentiate them as threshold or as interference effects, rather than pure resonances. Especially for scalar mesons, they are very difficult to identify due to large decay widths, and overlaps between resonances and background. Even the scalar mesons have the same quantum numbers as the vacuum. Thus, the understanding of the properties of scalar mesons may help us to understand the mechanism of symmetry breaking [6][7][8]. Various mesons in the low-lying sector like f 0 (980), a 0 (980), f 0 (1500) etc. were observed experimentally, but their specific structure and properties are not understood theoretically, nor experimentally. For example, the structure of the scalar states f 0 (1500) and f 0 (1710) have been debated since a long time [22][23][24][25][26][27][28][29][30][31][32][33][34]. The authors in [22][23][24] indicated them as mixed or rather pure glue states, respectively. References [26][27][28][29][30][31][32] indicated f 0 (1710) as a vector-vector molecular candidate and f 0 (1500) as a glueball. The status of these scalar states is still an open question. The study of these exotic states may give the answer to the basic questions of hadron masses, quark confinement, relevant degrees of freedom and interaction between the constituents of the same [6][7][8] ( Table 1).
In the present study, we have investigated the molecular structures in the light flavour meson regime. Our dimesonic system consists of a meson and an antimeson (qq − qq). The masses of these states would be less than the sum of the masses of two constituent mesons [35][36][37]. These loosely bound dimesonic system are similar to the deuteron like (proton-neutron) bound state system. This approximation had been taken previously in Refs. [36,37] and introduced as deusons'. Furthermore, the molecular like picture has been studied by Silvestre et al. [38,39] as diquonia.
In our model for the study of dimesonic states, we use the Hellmann potential with One Pion Exchange Potential (OPEP). The Hellmann potential is a superposition of the Coulomb and the Yukawa potential. The superposed potential of Coulomb and Yukawa was studied by Hellmann a long time ago [68][69][70][71][72]. There are several authors who applied this potential for the calculation of the bound state systems for different values of their strength parameters [73][74][75][76][77]. All the phenomenological potential model approaches used by the authors to study the mesons or multiquark spectra have dealt with complex interactions like, coulombic, confinement, instanton induced interaction, color, flavour and spin dependency to accomplish a description of the bound state properties. At very short distances, the bound state system is produced due to a delicate attractive and repulsive interaction. The Hellmann potential simply accomplishes the complicated theoretical calculation for the cancellation of the attractive and repulsive interaction at short distances, causing the small binding energy of the bound state. The Hellmann potential with the OPEP used in the present model takes care of the short and the long range behaviour of the interactions. As we are aware that mesons are color neutral, we have not introduced the confinement potential for the dimesonic systems. It is well studied that the OPEP is mainly responsible for the long range part of the strong nuclear force. The concept of meson exchange interaction has been used since a long time ago for the N-N interaction [78,79] as well as for hadronic molecules [36,37,[80][81][82].
To test the potential and mechanism employed, we first apply it to study the state f 0 (980), which was found in the literature to be interpreted as a molecular state [11][12][13]44]. Moreover, the molecular like picture has been supported by lattice QCD [83,84]. We calculate the spectra of f 0 (980) and obtained its mass, binding energy and root mean square radius. The results show fairly good agreement with the experimental results. Thus we are convinced and apply the same methodology to other dimesonic molecular states.
The study of the light meson sector needs a relativistic treatment. We are incorporating a relativistic correction to the kinetic energy term as well as to the meson-antimeson potential. The hydrogen like trial wave function is being used for this study. We have calculated the digamma decay width and decay width of the dimesonic systems. The study of the digamma decay of scalars helps to distinguish among different scenarios for the scalar meson structure [43]. The digamma decays of the dimesonic systems are calculated using the wave function at the origin [43]. For the decay calculation, we adopt the same formulation as used by Ref. [44], which was used to study the molecular picture of 1 −+ exotic states.
The article is organised as follows. After a brief introduction, we present the theoretical framework for the study of dimesonic molecular states in semirelativistic approach with the Hellmann potential and OPE potential in Sect. 2. The calculation of digamma decays and decay widths are discussed in Sect. 3. Then we present the calculated results in Sect. 4, and the last section is dedicated to our conclusions.

Theoretical framework
In the variational approach, we have solved the Schroedinger equation by using a trial wave function. The variational parameter is obtained by using the virial theorem. The (dimesonic) molecular system is assumed to consist of a meson-antimeson bound state. The Hamiltonian of the dimesonic system is given by [85][86][87][88] where m a and m b are the masses of mesons, P is the relative momentum of the two mesons and V (r ) is the molecular interaction potential of the dimesonic system. We expand the kinetic energy term of the Hamiltonian up to O (P 6 ): V (r ) is the molecular (meson-antimeson) potential, where We consider the meson-antimeson interacting through the Hellmann potential. It is a superposition of the Coulomb and Yukawa potentials of the form [68][69][70][71][72][73][74][75][76] Here α s and B are the strengths of the potential, C is the screening parameter and r 12 is the relative coordinates between meson and antimeson; the strength of the Coulomb potential is used as running coupling constant. B is the Yukawa potential strength and it takes both positive and negative values. Reference [74] has carried out a detailed study for different values of B and the screening parameter C for low-lying energy eigenvalues. For our dimesonic study, we have taken B positive. The value of the α s running coupling constant is determined through the model, namely GeV, and n f is number of flavours [89]. The non-perturbative form of V (1) (r ) is not yet known, but leading order perturbation theory yields where C F = 4/3 and C A = 3 are the Casimir charges of the fundamental and adjoint representation, respectively [90]. The relativistic mass correction is found to be similar to the Coulombic term of the static potential when applied to the charmonium and to be one-fourth of the Coulombic term for bottomonium [90]. Recently we have used this correction successfully in a study of the B c meson [85,86]. The long range OPEP is used based on the assumption that the molecular like structure of the multiquark system is as deuterium, like the structure of the nucleon [35,80]. The OPEP for the NN-interaction takes the form [91,92] where g 8 is a pion-nucleon coupling constant. τ and σ are isospin and spin factors, respectively. Here T π (r ), Y π (r ) and S 12 are The exchange meson having its own internal structure, to accomplish the finite size effect [91,92], the usual form factor appears due to the dressing of quarks and is assumed to be proportional to the exchange meson mass [80]. We have The functions T (r ) and Y (r ) with the finite size effects take the form A more detailed derivation of the potential can be found in Refs. [91,92]. S 12 is the usual tensor operator, it plays a prominent role in the NN-interaction. The matrix element of the tensor operator for L = 0 vanishes [82]. If two hadrons are in an L = 0 (ground) state, the term with a tensor operator vanishes as in our dimesonic case. Thus, the OPEP takes the form where g 8 = 0.69 is a meson pion coupling constant, π = km π , whereas m a and m b are constituents masses, m π = 0.134 GeV and k = 2.2. For PV states the values of the spin-isospin factor have been taken as τ i · τ j σ i · σ j = −3, 1 for I = 0, 1. For VV states the values are taken as τ i · τ j σ i · σ j = −6, −3, 3 for isospin I = 0 and spin S = 0, 1, 2, respectively; whereas for isospin I = 1 and spin The values of the spin-isospin factor are from Refs. [36,37] and the other parameters used in pion exchange potential are from Refs. [80][81][82].
For the hyperfine splitting, the spin dependent interaction potential is added perturbatively and takes the form The spin factor S 1 · S 2 can be found by the general formula We have used the hydrogen like trial wave function where μ is the variational parameter and L 2l+1 n−l−1 (μr 12 ) is the Laguerre polynomial. The ground state mass of the low-lying dimesonic states is calculated by obtaining the expectation value of the Hamiltonian and the variational parameter(μ) is obtained for each state by using the virial theorem [85][86][87][88]. We have We have included the kinetic energy term up to P 6 Table 2 [1]. The angular momentum, parity, spin and isospin all are conserved in the strong interaction, being good quantum numbers for the dimesonic system. The orbital angular momentum of meson and antimeson are L 1 and L 2 , respectively. In same way spin denote by S 1 , S 2 and isospin I 1 , I 2 . Employing the coupling rules, one has the relative orbital momentum and total spin of system L 12 and S 12 whereas the total angular momentum is J = L 12 + S 12 ; I = I 1 + I 2 .
The OPEP for the dimesonic system in Eq. (8) is spin-isospin dependent. Törnqvist in [36,37] has discussed the spin-isospin dependency of OPEP for deusons. The spin-isospin factor gives influence to the potential and decides whether a channel becomes attractive or repulsive. For VV states we got a repulsive channel for (S, I ) = (0, 1)(1, 1)(2, 0), while an attractive channel was found for (S, I ) = (0, 0)(1, 0)(2, 1). In the case of PV states the channel becomes attractive for (S, I ) = (1, 0), while it becomes repulsive for (S, I ) = (1, 1). Due to parity violation, the two pseudoscalars could not be bound by a pseudoscalar [35][36][37]. Thus, we could not consider OPEP for PP states. For such PP states, the Hellmann potential given in Eq. (5) is being used with a relativistic correction. The calculated masses are close to experimental measurements (PDG), tabulated in Table 3.
We have used a pseudoscalar and a vector meson for the dimesonic system, having the possible combinations (i) pseudoscalar-pseudoscalar (PP state), (ii) pseudoscalarvector (PV state), (iii) vector-vector (VV state).

Digamma width and decay width
The digamma decay of the dimesonic molecules are estimated, using the wave function at the origin, in analogy to the two photon decay of parapositronium [43]. The digamma width for all mass combinations is given by α = e 2 /4π indicates the fine-structure constant, we have the factor ξ = 1 √ 2 and m a and m b are the masses of the constituent mesons.
The decay width is calculated using the wave function at the origin. The dimesonic molecular state (constituent mesons a and b) decays into mesons c and d. In the decay through meson exchange, the amplitude is proportional to the wave function squared at the origin [44]. The formula of |ψ(0)| 2 given by [44,[94][95][96] is with M the reduced mass of the system, ψ(0) includes the effects of all ranges in r . The decay width is where m m is the mass of a dimesonic molecule and l is the magnitude of the 3-momentum of the decay product, given by m q and q are the mass and 3-momentum of the exchange mesons (for present work, we consider only a pion). is the free parameter which takes care of the off-shell effects at the vertices because of the internal structure of the mesons. The digamma decay widths are tabulated in Table 3 along with the masses. The decay width of the compared states as well as our predicated states are tabulated in Table 4.

Results and discussion
By using the Hellmann potential, the OPE potential and a relativistic correction, we have solved the Schroedinger equation to extract the masses, binding energy, decay widths and digamma decay widths for dimesonic states. In various combination of dimesonic states, like the PP, PV and VV states, we have calculated the possible molecular like structure in the light flavoured sector. Our calculated results are tabulated in Tables 3 and 4.
To describe the molecular picture of dimesonic states, we have incorporated the short range and long range behaviour of the potential. As the constituent of the molecule (meson) itself is color neutral and comes with a small binding energy and large radius (r in f m) in comparison to the constituent size, we do not need to introduce the confinement potential and so we have not dealt with the fundamental quark gluon condense state in the present study. Thus, we have not incorporated any mixing scheme of gluonia and quarkonia. Still, the results obtained in the present study have reasonably good agreement with experimental measurements and explain the molecular picture well.
We have expanded the kinetic energy term up to O(P 6 ). In the series expansion of the kinetic energy, For v c, the effect of the higher order term of the momentum P 2n (n > 2) is negligible, even more, the higher order term has poor convergence. While the expansion term up to P 4 does not have a lower bound. So, the usable expansion to incorporate the relativistic effect is being up to P 6 [91,92]. We have dealt with the systems which have constituent masses below 2 GeV.
For the lighter constituents, the effect of the momentum is small. As v < c tends to v c, the effect of the higher order terms contributes less than 1 %. The effect of the higher order terms (up to P 6 ) is very small, still, it has variation with increasing masses and contributes (v → c) to the net kinetic energy of the system. Therefore, it is justified to incorporate it into the expansion up to P 6 .
The correction used in the potential has its dominant effect on the potential energy. In our calculation, the effect of the correction to the potential energy part is increasing as mass of the system decreasing, and, it is about to be 25-40 %. The effect of correction to potential is approximately 38 % in the lightest dimesonic system (K − K ) of this work, while the heaviest system φ − φ has a contribution of about 28 %. We have calculated the partial decay width as per Eq. (22). The parameter is assumed to be proportional to the mass of the dimesonic state, = Km m . We have chosen the Yukawa strength parameter and the color screening parameter B = 0.4 and C = 0.35, respectively, to find the mass spectra for the K − K molecule, which is believed to be f 0 (980). For the decay calculation of K − K , we took K = 1.262. Our calculated results for mass, digamma and decay width are close to the PDG values [1]. Thus, we fixed these parameters for the calculation of all PP and PV states, while for the VV state we took B = 1.2. For the dimesonic states having masses between 1 to 1.6 GeV, the constant K = 1.085, whereas for the states between 1.6 to 2 GeV, K = 0.931. In the PP states, we used the Hellmann potential with a relativistic correction. Due to parity violation, the OPEP cannot be applied to the calculation of the PP states. We have Our results for f 0 (980) k K are in good agreement with the suggested mass and decay width (0.982 ± 0.003 and 80 ±10 MeV) of Ref. [97]. The calculated digamma width of f 0 (980) is in close agreement with other theoretical predictions [48,[97][98][99][100][101]. The ππ decay mode of f 0 (980) is dominant and also agrees with Ref. [1]. For the molecular state, one should have a small binding energy and a large radius (compared to the constituent energy and radius). Our results satisfied these conditions. The mass and decay properties are close to the experimental measurements. Thus, we employ the same methodology to calculate the properties of all the meson-antimeson combinations. The computed results are tabulated in Tables 3 and 4.

Vector-vector (VV) states
In the VV combination, the dependency of one pion interaction potential on the spin-isospin factor gives possibilities of the numbers of combinations for dimesonic molecular like states [36,37,82]. We have already discussed in the previous section about the matrix element that we have attractive and repulsive channels of the spin-isospin factors. We have found bound states in repulsive spin-isospin channels, showing the dominance of the Hellmann potential and the mass correction over pion exchange. This can be clearly seen from the binding energy tabulated in Table 4 for a state of the same spin with different isospin. With a different spin-isospin combination, we have compared experimentally observed states with our dimesonic molecular states, with comparable mass spectra and quantum numbers. The states f 0 (1500), f 2 (1525), f 2 (1565), a 0 (1450), f 0 (1710), a 2 (1700) and f 2 (1810) compared with different vector-vector combinations are shown in Table 3. f 0 (1500) is compared as a ρ − ρ molecule with (S, I ) = (0, 0), falling in the 1500 GeV mass regime. The computed mass and decay widths are comparable with experimental values [102,103] and with Refs. [102,103], with a suggested mass of 1.522±0.005 GeV and a decay width of 108±8 MeV. In Ref. [97] is suggested a decay width of 131 ± 15 MeV, and also in Refs. [104][105][106][107][108]. We have found that the ππ decay mode is dominant. Experimentally two photon decay has not been seen yet. We have m m f 0 (1500) = 1.489 GeV; γ γ = 0.079 KeV; (h 1 (1500)) = 110 MeV.
The state f 2 (1565) is compared with the ω − ω molecule. The calculated mass of this state is underestimated by around 40 MeV with the PDG value [1], and the calculated digamma decay width is being underestimated by about one-ninth of the recent PDG value [1]. Even more, the calculated decay width for the observed decay modes are also far away from the PDG values [1]. We have For f 2 (1565) m m f 2 (1565) = 1.514 GeV; γ γ = 0.082 KeV; ( f 2 (1565)) = 81 MeV.
The mass of a 2 (1700) k * K * is in excellent agreement with the experimental PDG value [1]. The digamma decay is found at one-third of the PDG value (0.30 ± 0.05 KeV) [1], but the decay width is in good agreement with the experimental values for the observed decay modes [1]. We have m m a 2 (1700) = 1.736 GeV; γ γ = 0.098 KeV; (a 2 (1700)) = 181 MeV.
Anisovich et al. [117] made a remark on the status of the state f 2 (1810) such that the state may be confused with f 0 (1790).
Reference [117] has note that, if there was not any confusion between f 2 (1810) and f 0 (1790), there should be a missing 0 ++ state. As we observed throughout our calculated results, all states (except a few states) are underestimated by approximately 20-40 MeV. If we consider it as a limitation of our model, we can compare the state f 0 (1790) with φ − ω, with (S, I ) = (0, 0). So it may be possible that these two states f 2 (1810) and f 0 (1790) have the molecular structure of φ − ω with spin state S = 0, 2, respectively. But the experimental status of neither state is confirmed. Thus, it needs more attention theoretically as well experimentally.
Our results favoured the molecular picture. The status of the scalar states f 0 (1500) and f 0 (1710) as a molecule, guleball or glue mixed state have been debated since a long time ago [22][23][24][25][26][27][28][29][30][31][32][33][34]. The glueball and quarkonia with the same quantum number and a nearly equal mass may have interfered with each other and formed a new state. If the one bare glueball and one or more quarkonia interfere, one needs the mixing scheme to explain them. This interference or mixing of the state could be solved in a mass matrix formalism or in a linear combination of pure quarkonia and gluonia states by finding the linear coefficients, which mainly depends on the decay properties [23]. The authors of Ref. [22] have predicated f 0 (1500) as a glue mixture with ss component and f 0 (1710) as predominantly a pure glueball, as also in Refs. [23,24]. In Ref. [25] one has indicated the f 0 (1710) as a glueball candidate. References [26][27][28][29][30][31][32] have indicated f 0 (1710) as a vector-vector molecular candidate and f 0 (1500) as a glueball. So far, as per the present literature, the structure has remained unclear and debatable of the f 0 (1500) and f 0 (1710). Moreover, in the recent review of particle data [1] one has noted the status of the scalar states, in which f 0 (1710) together with f 2 (1525) were interpreted as bound systems of two vector mesons [1,30]. The molecular picture could be tested in radiative J/ψ decays as well as radiative decays of the states themselves [1,33,34]. The vector component of f 0 (1710) might have the origin of the enhancement seen in J/ψ −→ γ ϕω near threshold as observed at BES [1,31,32]. The scenario of these states will may become clear in future experiments like PANDA at FAIR.
Beside these, we have also calculated the decay width of the remaining dimesonic states and as listed in Table 4, these states may be identified in future experimental measurements.

Conclusion
In this paper, we are able to calculate the S-wave masses, digamma and decay widths of the dimesonic states in the light meson sector. Here, we would like to say that the dimesonic model becomes more accurate if the binding energy of the constituents is small compared to their masses. We compared many states which are experimentally observed and predicted theoretically and having a promising non-qq structure as the dimesonic states. The calculated results are in fairly good agreement with experimental measurements as well as theoretical predictions. On the basis of mass spectra, decay widths and in some cases digamma widths, we are able to identify some dimesonic (meson-antimeson) states. We have calculated the decay widths for experimentally observed decay modes for each of the compared states.
In our present study, we strongly recommend the states f 0 (980), b 1 (1235), f 2 (1525), f 0 (1710) as meson-antimeson molecules. Furthermore, the states h 1 (1380), a 0 (1450), f 2 (1565) and h 1 (1595) in the calculation have been found with deviated masses and decay widths from the experimental measurements, which rules out their candidature as dimesonic molecules. Additionally, some states like f 0 (1500), a 2 (1700) and f 2 (1810) as regards the calculated results have a little bit of variation in the decay widths, but in some decay modes, the decay widths are fairly close to some experimental results [104][105][106]120,129]. Some states like a 2 (1700), f 2 (1810) and f 0 (1790) have an unconfirmed experimental status, needing more attention from experimentalist as well as from theoreticians.
Finally, we have predicted the masses, decay widths and digamma widths of the S-wave dimesonic molecular states (light sector). In the future we would like to employ this model to calculate the P-state masses and decay properties for light sector as well as for the heavy sector dimesonic systems.