# Spin Splitting Spectroscopy of Heavy Quark and Antiquarks Systems

Authors: Hesham Mansour; Ahmed Gamal; M. Abolmahassen
Journal: Advances in High Energy Physics
 (2020)

---

## Abstract

Phenomenological potentials describe the quarkonium systems likecc¯,bb,¯andb¯c where they give a good accuracy for the mass spectra. In the present work, we extend one of our previous works in the central case by adding spin-dependent terms to allow for relativistic corrections. By using such terms, we get better accuracy than previous theoretical calculations. In the present work, the mass spectra of the bound states of heavy quarks cc¯, bb¯, and 𝐵𝑐 mesons are studied within the framework of the nonrelativistic Schrödinger equation. First, we solve Schrödinger’s equation by Nikiforov-Uvarov (NU) method. The energy eigenvalues are presented using our new potential. The results obtained are in good agreement with the experimental data and are better than the previous theoretical estimates.


---

## Body

## 1. Introduction

In the twentieth century, quarkonium systems have been discovered. Theorists have been trying to explain some aspects of those systems like mass spectra and decay mode properties. [1–5]. Some of them used lattice quantum chromodynamics view [6–12], effective field theory [13], relativistic potential models [14, 15], semirelativistic potential models [16], and nonrelativistic potential models [17–19] which have shared in common the Coulomb and linear potentials. There are other groups which use confinement power potential rn [20–22], the Bethe-Salpeter approach [23–25]. In the present work, we use mixed potential: nonrelativistic potential models (Coulomb+linear) and confinement power potentials plus spin-dependent splitting terms as a correction. Schrödinger’s equation is solved by the Nikiforov-Uvarov (NU) method [26–34], which gives asymptotic expressions for the eigenfunctions and eigenvalues of the Schrödinger equation.


## 2. Methodology

In the quarkonium system which deals with quark and antiquark interaction in the center of mass frame, the masses of the quark and antiquark are bigger than chromodynamics scaling, i.e.,Mq,q¯≫ΛQCD. So, this allows for nonrelativistic treatment and is considered as heavy bound systems. Using Schrödinger’s equation of the two-body system in a spherical symmetric potential one obtains the following:
(1)d2Qrdr2+2μℏ2E−Vtotr−ll+1r2Qr=0,where μ is reduced mass, E is energy eigenvalue, l is orbital quantum number, Vtotr is total potential of the system, Qr=rRr, and Rr is a radial wave function solution of Schrödinger’s equation.Our radial potential is taken as follows:
(2)Vr=−br+ar+dr2+pr4.Also, we use spin-dependent splitting terms: spin-spin interaction, spin-orbital interaction, and tensor interaction, respectively [35–43].
(3)Vs−l−Tr=VS−Sr+Vs−lr+VTr,where
(4)VS−Sr=23mqmq¯∇2VVrSq→∙Sq¯→,Vs−lr=12mqmq¯r3dVVrdr−dVsrdrL→∙S→,VTr=112mqmq¯1rdVVrdr−d2VVrdr26Sq→∙r→rSq¯→∙r→r−2Sq→∙Sq¯→.VV is a vector potential term, and Vs is a scalar potential term.So, the total potential becomes(5)Vtotr=Vr+Vs−l−Tr,Vtotr=4av3mqmq¯rvss+vlsmqmq¯3av−aS2r+3b2r3−d−2pr2+vT12mqmq¯avr+3br3−br+ar+dr2+pr4,where vss=Sq→∙Sq¯→,vls=L→∙S→,vT=−2Sq→∙Sq¯→−3Sq→∙r→/rSq¯→∙r→/raS+av=a.By substituting in equation (1), we get
(6)d2Qdr2+ε−Br3−ll+1r2−Cr−Ar−Fr2−Pr4Q=0,where in natural units,
(7)ε=2μE+2μdvlsmqmq¯,(8)B=μb6vls+vt2mqmq¯.(9)A=2μa,F=2μmqmq¯−2pvls+dmqmq¯,(10)C=μ16avvss+63av−aSvls+avvT−12mqmq¯b6mqmq¯,P=2μp.Let x=1/r, and by substituting in equation (6), we obtain
(11)d2Qdx2+2xdQdx+1x2εx2−Bx−Cx−ll+1−Ax3−Fx4−Px6Q=0.In equation (11), one can use the Nikiforov-Uvarov method (NU) to get eigenvalue and eigenfunction equations.Due to the singularity point in equation (11), put y+δ=x, and using the Taylor series to expand to second order terms, one obtains
(12)d2Qdx2+2xdQdx+1x2−q+wx−zx2Q=0,where
(13)−6εδ2+3Cδ+ll+1+10Aδ3+15Fδ4+28Pδ6=q,(14)−8εδ3−B+3Cδ2+15Aδ4+24Fδ5+48Pδ7=w,(15)−3εδ4+Cδ3+6Aδ5+10Fδ6+21Pδ8=z.We get
(16)z=w2n+1+2q+1/42.By substituting equations (13)–(15) in equation (16) and arrange it, we get
(17)ε=Cδ3+2Aδ+10F3δ2+7Pδ4−133C+15A/δ2+24F/δ3+48P/δ5−8ϵ/8ϵ−Bδ22n+1+2l+1/22+3C/δ−6ϵ/δ2+10A/δ3+15F/δ4+28P/δ62.We substitute equations (8)–(10) into equation (17) to obtain the energy eigenvalue equation.
(18)E=ℵ−μ6ξ2n+1+2l+1/22+μη2,where δ=1/r0,ℵ=16avvss+63av−aSvls+avvt−12mqmq¯b/36mqmq¯r0+2ar0+10/3mqmq¯−2pvls+dmqmq¯r02+7pr04−dvls/mqmq¯,(19)ξ=16avvss+63av−aSvls+avvt−12mqmq¯b2mqmq¯+30ar02+48mqmq¯−2pvls+dmqmq¯r03+96pr05−8r02E+2dvlsmqmq¯−b6vls+vt2mqmq¯r02,η=16avvss+63av−aSvls+avvt−12mqmq¯br02mqmq¯+20ar03+30mqmq¯−2pvls+dmqmq¯r04+56pr06−6r022E+2dvlsmqmq¯.Knowing that(20)Mqq¯=E+mq+mq¯→E=Mqq¯−mq+mq¯.So, the mass spectra equation becomes(21)Mqq¯=ℵ−μ6ξ2n+1+2l+1/22+μη2+mq+mq¯.The eigenfunction equation is(22)Qr=Nnlsr1/2−q+1/4e−z/rLn2q+1/42zr,whereLn2q+1/42z/r is the Rodrigues formula of the associated Laguerre polynomial and Nnls is a normalization constant.So, the radial wave function solution of Schrödinger’s equation is given by(23)Rr=Nnlsr−1/2−q+1/4e−z/rLn2q+1/42zr.The energy eigenvalue equation (18) has spin-orbital-tensor coefficients vss,vsl,vT, and those can be given from references [35–37]. Also, it has potential parameters aS,av,b,d,p and r0 due to the expansion, so we have six parameters of the eigenvalue equation which can be obtained from the experimental data [44] by best fitting as shown in Figure 1.Figure 1
This is a flow chart about how one can calculate the theoretical data. It consists of three stages to do that. We found that after some iteration processes in the loops, the previous suggested values of the coefficients are the same in all last processes, as there is no change in the data. So, we used another program which gives a small change in the parameter values which was Origin lab program. In the third stage, after we got the values of the coefficients of the potential with the smallest chi-square values from the Origin program, we will substitute them in the energy eigenvalue equation and calculate the theoretical data (present work) for the differentn and L states by using “fsolve” built in function in MATLAB like first flow chart.


## 3. Numerical Results and Discussions

In Table1, potential parameters are shown for each system. It is noticed that the values of these parameters are different for different systems, and this is due to the properties of those systems like energy scale and decay mode. We use spectroscopic notation for the levels n2S+1LJ.Table 1
Parameter values of each system.
Variablesmqmq¯r0asavbdpSystemsUnitsGeVGeVGeV-1GeV2GeV2--GeV3GeV5ccsystem1.3171.31712.82-0.0796-0.73490.0090.10686-0.000184bbsystem4.5844.5847.237950.05052.577115.02-0.07969-0.00102b¯csystem4.5841.31711.434-2.55491.53690.0390.11453-0.00018438S is the total spin of the system, L is the orbital quantum number, n is the principal quantum number, and J is the total (orbital+spin) quantum number.By using equation (21) and Table 1, we get the mass spectra of different quantum states as shown in Tables 2–7. Previously, we used the phenomenological potential in equation (2) without spin-dependent corrections (central-dependent potential) [45]. The results obtained were good in comparison with the experimental data.Table 2
Charmonia mass spectrum ofS and P waves in GeV.
LevelPresent work[24][46][37][47][48][15][49][40][50][51]Exp. [44]11S03.0332.932.9812.9842.9892.9792.9802.9802.9823.0882.9792.984 (49)13S13.1263.113.0963.0973.0943.0973.0973.0973.0903.1683.0963.097 (29)21S03.6663.683.6353.6373.6023.6233.5973.6333.6303.6693.6003.639 (27)23S13.7013.683.6853.6793.6813.6733.6853.6903.6723.7073.6803.686 (15)31S04.158--3.9894.0044.0583.9914.0143.9924.0434.0674.011--33S14.0553.804.0394.0304.1294.0224.0954.0304.0724.0944.0774.039 (16)41S04.415--4.4014.2644.4484.2504.4334.2444.3844.3984.397--43S14.415--4.4274.2814.5144.2734.4774.2734.4064.4204.4544.421 (6)51S04.607--4.8114.4594.7994.446--4.440--------53S14.585--4.8374.4724.8634.463--4.464--------61S04.754--5.155--5.1244.595--4.601--------63S14.733--5.167--5.1854.608--4.621--------13P03.4073.323.4133.4153.4283.4333.4163.3923.4243.4483.4883.415 (8)13P13.4873.493.5113.5213.4683.5103.5083.4913.5053.5203.5143.511 (24)11P13.5023.433.5253.5263.4703.5193.5273.5243.5163.5363.5363.525 (23)13P23.5223.553.5553.5533.4803.5563.5583.5703.5563.5643.5653.556 (33)23P03.8993.833.8703.8483.8973.8423.8443.8453.8523.8703.9473.918 (19)23P13.7863.673.9063.9143.9383.9013.9403.9023.9253.9343.972--21P13.8213.753.9263.9163.9433.9083.9603.9223.9343.9503.996--23P23.905--3.9493.9373.9553.9373.9943.9493.9723.9764.0213.927 (21)33P04.120--4.3014.1464.2964.131--4.1924.2024.214----33P14.1233.914.3194.1924.3384.178--4.1784.2714.275----31P14.164--4.3374.1934.3444.184--4.1374.2794.291----33P24.144--4.3544.2114.3584.208--4.2124.3174.316----43P04.362--4.698--4.653--------------43P14.373--4.728--4.696--------------41P14.42--4.744--4.704--------------43P24.411--4.763--4.718--------------53P04.543------4.983--------------53P14.561------5.026--------------51P14.611------5.034--------------Table 3
Charmonia mass spectrum ofD and F waves in GeV.
LevelPresent work[47][46][37][24][48][15][49][40][50][51]13D33.3073.7553.8133.8083.8693.7993.8313.8443.8063.8093.79811D23.3763.7653.8073.8053.7393.7963.8243.8023.7993.8033.79613D23.3483.7723.7953.8073.5503.7983.8243.7883.8003.8043.79413D13.3743.7753.7833.792--3.7873.8043.7293.7853.7893.79223D33.7974.1764.2204.1123.8064.1034.2024.1324.1674.1674.42521D23.8364.1824.1964.108--4.0994.1914.1054.1584.1584.22423D23.8014.1884.1904.109--4.1004.1894.0954.1584.1594.22323D13.8004.1884.1054.095--4.0894.1644.0574.1424.1434.22233D34.1634.5494.5744.340--4.331--4.351------31D24.1764.5533.5494.336--4.326--4.330------33D24.1354.5574.5444.337--4.327--4.322------33D14.1134.5554.5074.324--4.317--4.293------43D34.4364.8904.920--------4.526------41D24.4294.8924.898--------4.509------43D24.3834.8964.896--------4.504----43D14.3454.8914.857--------4.480------13F23.4033.9904.041------4.068--4.029----13F33.3754.0124.068--3.999--4.070--4.029----11F33.4034.0174.071--4.037--4.066--4.026----13F43.3154.0364.093------4.062--4.021----23F23.8124.3784.361----------4.351----23F33.8234.3964.400----------3.352----21F33.8584.4004.406----------4.350----23F43.8144.4154.434----------4.348----33F24.1114.730------------------33F34.1524.746------------------31F34.1944.749------------------33F44.1864.761------------------Table 4
Bottomonia mass spectrum ofS and P waves in GeV.
LevelPresent work[52][46][37][24][53][15][49][54]Exp. [44]11S09.4729.4029.3989.3909.4149.3899.3939.3929.4559.398 (74)13S19.5259.4659.4609.4609.4909.4609.4609.4609.5029.460 (65)21S010.0289.9769.9909.9909.9879.9879.9879.9919.9909.999 (28)23S110.04910.00310.02310.01510.08910.01610.02310.02410.01510.023 (26)31S010.3610.33610.32910.326--10.33010.34510.32310.330--33S110.37110.35410.35510.34310.32710.35110.36410.34610.34910.355 (16)41S010.59210.52310.57310.584--10.59510.62310.558----43S110.59810.63510.58610.597--10.61110.64310.57510.60710.579 (19)51S010.7910.86910.85110.800--10.817--10.741----53S110.8710.87810.86910.811--10.831--10.75510.81810.876 (6)61S010.96111.09711.06110.997--11.011--10.892----63S111.02211.10211.08810.988--10.988--10.90410.99511.019 (3)13P09.849.8479.8599.8649.8159.8659.8619.8629.8559.859 (19)13P19.8759.8769.8929.9039.8429.8979.8919.8889.8749.893 (18)11P19.8849.8829.9009.9099.8069.9039.9009.8969.8799.899 (15)13P29.9039.8979.9129.9219.9069.9189.9129.9089.8869.912 (9)23P010.20210.22610.23310.22010.25410.22610.23010.24110.22110.232 (30)23P110.22910.24610.25510.24910.12010.25110.25510.25610.23610.255 (26)21P110.23710.25010.26010.25410.15410.25610.26210.26110.24010.260 (23)23P210.25410.26110.26810.264--10.26910.27110.26810.24610.269 (15)33P010.29910.55210.52110.490--10.502--10.51110.500--33P110.33910.53810.54110.51510.30310.524--10.50710.513--31P110.36210.54110.54410.519--10.529--10.49710.516--33P210.40610.55010.55010.528--10.540--10.51610.521--43P010.53210.77510.781----10.732--------43P110.57110.78810.802----10.753--------41P110.59410.79010.804----10.757--------43P210.63710.79810.812----10.767--------53P010.73111.004------10.933--------53P110.76911.014------10.951--------51P110.79211.016------10.955--------Table 5
Bottomonia mass spectrum ofD and F waves in GeV.
LevelPresent work[52][46][37][24][53][15][49][54]Exp. [44]13D39.84910.11510.16610.15710.23210.15610.16310.17710.127--11D29.76710.14810.16310.15310.19410.15210.15810.16610.123--13D210.09610.14710.16110.15310.14510.15110.15710.16210.12210.163 (66)13D19.66610.13810.15410.146--10.14510.14910.14710.117--23D310.17510.45510.44910.436--10.44210.45610.44710.422--21D210.09310.45010.44510.432--10.43910.45210.44010.419--23D210.07110.44910.44310.432--10.43810.45010.43710.418--23D19.99610.44110.43510.425--10.43210.44310.42810.414--33D310.44610.71110.717----10.680--10.652----31D210.36810.70610.713----10.677--10.646-----33D210.34510.70510.711----10.676--10.645----33D110.27210.69810.704----10.670--10.637----43D310.67610.93910.963----10.886--10.817----41D210.59910.93510.959----10.883--10.813----43D210.57610.93410.957----10.882--10.811----43D110.50410.92810.949----10.877--10.811----13F29.64210.35010.34310.338----10.353--10.315--13F39.75410.35510.34610.34010.302--10.356--10.321--11F39.77810.35510.34710.33910.319--10.356--10.322--13F49.89610.35810.34910.340----10.357------23F29.97110.61510.610------10.610------23F310.08110.61910.614------10.613------21F310.10410.61910.647------10.613------23F410.21910.62210.617------10.615------33F210.24610.850----------------33F310.35310.853----------------31F310.37610.853----------------33F410.48910.856----------------Table 6
Bc meson mass spectrum of S and P waves in GeV.
LevelPresent work[47][55][46][56][57]Exp.[44]11S06.2766.2726.2786.2726.2716.2756.275 (0)13S16.3136.3216.3316.3336.3386.314--21S06.8416.8646.8636.8426.8556.8386.842 (1)23S16.8676.9006.8736.8826.8876.850--31S07.2817.3067.2447.2267.250----33S17.3087.3387.2497.2587.272----41S07.6347.6847.5647.585------43S17.667.7147.5687.609------51S07.9178.0257.8527.928------53S17.9418.0547.8557.947------61S08.1448.3408.120--------63S18.1688.3688.122--------13P06.2236.6866.7486.6996.7066.672--13P16.2816.7056.7676.7506.7416.766--11P16.2906.7066.7696.7436.7506.828--13P26.3666.7126.7756.7616.7686.776--23P06.7827.1467.1397.0947.1226.914--23P16.8367.1657.1557.1347.1457.259--21P16.8467.1687.1567.0947.1507.322--23P26.9177.1737.1627.1577.1647.232--33P07.2277.5367.4637.474------33P17.2787.5557.4797.510------31P17.2877.5597.4797.500------33P27.3557.5657.4857.524------43P07.5837.885--7.817------43P17.6317.905--7.853------41P17.6407.908--7.844------43P27.7047.915--7.867------53P07.8678.207----------53P17.9138.226----------51P17.9228.230----------Table 7
Bc meson mass spectrum of D and F waves in GeV.
LevelPresent work[47][55][46][56][57]13D36.4296.9907.0267.0297.0456.98011D26.3086.9947.0357.0267.0417.00913D26.2996.9977.0257.0257.0367.15413D16.2006.9987.0307.0217.0287.07823D36.9757.3997.3637.405----21D26.8617.4017.3707.400----23D26.8527.4037.3617,399----23D16.7597.4037.3657.392----33D37.4097.761--7.750----31D27.30177.762--7.743----33D27.2927.764--7.741----33D17.2057.762--7.732----43D37.7548.092--------41D27.6538.093--------43D27.6438.094--------43D17.568.091--------13F26.1837.234--7.2737.269--13F36.3267.242--7.2697.276--11F36.3357.241--7.2687.266--13F46.5017.244--7.2777.271--23F26.7417.607--7.618----23F36.8767.615--7.616----21F36.8857.614--7.615----23F47.0417.617--7.617----33F27.1867.946--------33F37.3147.954--------31F37.3237.953--------33F47.4707.956--------


## 4. Conclusions

The above tables show that spin-dependent terms are important factors to give a better accuracy and complete quantitative description of the quarkonium systems for the cases where experimental values are available. The theoretical work agrees with the experimental data. This shows also that the Nikiforov-Uvarov method is a good method to get the energy eigenvalues for the meson spectra. The results are even better than other previous works.



---

*Source: 2356436-2020-11-24.xml*