Alireza Khalaghi,Majid Monemzadeh,and Narges Tazimi

Department of Physics,University of Kashan,Ghotb Ravandi Boulevard,Kashan,Iran

Abstract Truly by looking through the analytical model of constituent quarks and further the importance of the effects of relativity on quark dynamics in studying the interior structure of mesons,in this research we have strived to have a much more precise modeling for quark interior structure.Certainly by observing the constituent model of quarks,at first we consider the mesons as two-body system,then we place the considered calculated Potential,which is a function of location and spin,in Schr?dinger’s equation.Next we will solve the mentioned equation in analytical method.Moving on this solution,we will import the spin and isospin interactions as perturbation in our problem,and finally by using these solutions we can obtain both binding energy and wave function for bound state and excited states of meson.Eventually,by applying these calculations in the next and last step we will compare our data about meson’s binding energy and masses with others results.

Key words: mesons,binding energy,spin and isospin interactions,Schr?dinger’s equation,exact solution

1 Introduction

In particle physics,the mesons are hadron subatomic particles,which consist of a quark and an anti-quark binding with strong interaction.Since the mesons are composed of quarks,their physical size is about one femtometre,about two thirds of the size of the proton or neutron.All the mesons are unstable,with the longest-lived lasting for only a few hundredths of a microsecond.The discharge mesons may collapse into the photon and both of these disintegrations indicate that the color is no longer a side product.[1?5]

The potential between quarks and antiquarks in mesons can be considered in different ways,nonetheless the best practical potential is the one with all properties of quark-antiquark pair; In addition,the power of color is created by replacing the gluons between the quark and the antiquark,as each quark is influenced by the force which comes from another one; therefore,the central potential should be considered as the following configuration to be closer to the laboratory evidence.[6]

In order to study the two-quark systems known to be two-particle systems,also by Schrodinger and Dirac equation solution,various potentials in strong interactions between quarks have been proposed; however,the net harmonic oscillator potential plus the square of the distance term is applied in our calculation.The considered potential is:[7]

To calculate the mass spectrum of meson,we are struggling to solve the non-relativistic model of Schr?dingers equation with a QCD potential; nonetheless,some potential models are characterized by the flavor independence of the selected potential and also the presence of a confining term.The quark masses used in our model are the constituent quark masses that are the mass parameters appearing in the QCD Lagrangian.[7]They differ from current quark masses.Constituent quark masses are larger than current quark masses.Thus could result from gluonic condensate effects.Constituent quark masses are free parameters used to fit in potential models and therefore,various values assigned to them are observed in different studies in the literature.

In this form,this added term to the harmonic potential,(squared distance),provides us a precise calculation in the presence of short-range interactions,and the constant W is manually used for proportional masses beforehand in order to obtain an accurate outcome in the spectrum computation.

Following the information on this potential we can mention the constants as below:[7]

Then the spin-spin,isospin-isospin and spin-isospin interactions potentials are introduced into the problem solving sentences and we obtain the contribution to the mass spectrum calculation of each meson as follows.The mass fragmentation in several layers of hydronium may be due to various factors; as a result,an ultrafine fission can be seen.

In the following terms,the spin-spin,isospinisospin and interaction potentials can be described as follows:[8?13]

whereσST=2/3 fm2,andAST=106/2 fm2.[14?16]

The perturbation potential,the sum of the high interaction potentials,is now defined,hence we solve the Schrodinger equation for the first potential then we consider the potential for the disturbance,and the meson spectrum will be procured.By comparing the results with the masses obtained for each meson,it can be determined that the quarks which are forming each meson are in single or triple state,and this procedure can be used to compute the states of each meson.Given the interaction and the superconducting interaction potential,we can consider the interaction term of a meson as follows:

2 The Calculation of Mesons Binding Energy

Using the radial coordination we will predict the Schr?dinger radial equation as:[17]

Now,considering the variable alternation,which is a function of(r)asψ(r)=(1/r)?(r)and the substitution of that in Eq.(6) we can obtain the following term:

which the following mentioned function?(x) is proposed to be as below:

Furthermoref(x) is like Hermitian polynomials and the coefficientαviis obtained in terms of potential coefficient.Moreover,we assume and by importing it in mentioned Eq.(8) we can have:

By maintaining Eq.(11),it is enough to set the potential interactions that we have here to be the same as mentioned pure oscillator in squared distance to achieve the relation below:

Now in the ground state we choose (f(r)=1),then we paste the term:g(r)=(1/2)βr2+δlnrin Eq.(12) and we will reach the underneath equation:

Considering the relation below based on the equation above:

And takingl=0 we will have the ground state energy for mesons as below:

Applying the calculated numbers in Eq.(2)in meson mass relation:Mqˉq=mq+mˉq+ξ,shall afford us with meson’s mass in ground state.[7,18?20]Our calculation is to be done for the s,c,b quarks by recognizing their fitted masses as the listed in Table 1.[7]

Table 1 Quark fitted mass (×103 MeV).[7]

At the beginning,we set theW-value to fit each meson with proportional masses,so that the empirical experiment has an appropriate proximity,then we will consider an averageWfor all the mentioned meson structures.We can witness an example for Charmonium below:

We choose the=?4356 (MeV) so the mass can be the exact same quantity as the following number=2983 (MeV).[21]We shall do the similar stages for other structures and we can acquire various values onWwhich have been transitioned in Table 2.[21]

Table 2 The W constant calculated for every meson manually.

The last one is named?=,which is not the intentions of our computation as a reason of the light mass it has in comparison with other structures; however,it is probably a conclusive proof on our theory,which does not work on light mesons class that we will find out later on this article.Finally the interval can be easily summed up with 6 more constants,and their division into 6 equals toWrefrence=Wref=?4.306×103(MeV),hence new masses are exposed within Table 3.

Table 3 Meson masses (×103 MeV) by the constant Wref.

By repeating the calculations for all the particles with this reference constant,we can observe that the obtained constant has a suitable value.Moving on this computation,we can manage to have the amounts of ground state energy for mesons either by the reference constant (W)which is going to be achieved by the same direction we applied for masses although this time with Eq.(12) for ground state energy.The percentage of their very small errors both in masses and ground state energy in the next section is displayed as a table.

3 Calculation of Meson’s Masses in the Excited State

In order to calculate the meson mass in the first excited state we should consider:

In this case:

which we can interpret the equation above as below:

As a result it is led to:

Henceforth,the energy relation will be written as below:

Now,by multiplying a negative in Eq.(20e),and making it a perfect square,δwould be written as below:

Given the above equations,ifl=0 thenδ=1/2 is derived.

Acknowledging the reduced mass which was mentioned earlier,we can achieve the exited state energy relation,in addition the exited mass formula which will be derived as:

in which quark quantities and masses can be used from the previous sections as follows:[7]

Considering a new reference constant (W) is an obligation,which has been applied in previous section,toward having precise quantities on both mass and energy.Now,as before,we obtain a manual value forWin order to receive the closest answer to the experimental measure.Thus for strangeBmeson we have:[22]

By assuming=?3981.72 (MeV) the mass is obtained=51415.4 (MeV).For the three remaining structures we follow the same path and we shall find the terms below:[23?25]

In this part,as it has been investigated in Sec.2,by averaging on the four obtained equations,we introduce a reference constant though in an excited state:

The overriding point in this calculation is the quark masses,which should be constructed of the similar fitted cited masses as we mentioned earlier in the first table.In addition,the meson mass of the applied structures in the excited states must exist.Therefore,we are only calculating the masses and energies for mentioned structures given in Eq.(25).Based on Eq.(23)for excited state energy we shall have the values represented in Table 4.[22?25]

Regarding the values obtained from calculations in Secs.2 and 3,the percentage of errors of listed meson’s masses in the ground state and the excited state are available to be shown in Tables 5–7.(the following tables.)

Table 4 The ground state meson masses according to Refs.[21–25].

Table 5 The excited state meson masses according to Refs.[21–25].

Table 6 Investigating the values of energy in the ground states.[22?25]

Table 7 Investigating the values of energy in the excited states.[22?25]

4 Investigating the Spin-Spin Interaction Potential Impact

In this regard,in order to refine the solutions,we will consider the interaction sentence as a perturbation in our calculations:

whereS1andS2are the quark and antiquark spins,respectively.This interaction potential can be written as:

whichcan be derived as below:

And so far we define Ψ(r) and?(r) as below:

By replacing the achieved functions in the following mentioned interaction

we can reach to the main formula:

In this part we need to solve the integrals which by default limitations,the answer of the numerator integral and the total fraction will be divergent;hence to unravel this problem,the potential must be set to zero with the reference constant for the meson,in order to obtain an approximate numerical value for meson radius.Because of the strong interaction among these particles,they do not interact with distance far from their radius; therefore,we consider this distance as physical infinity,(physically unlimited).

Consideringα=1/8μand computing this amount for charmonium.Given the numerical values forK,Wrefused in Sec.2 to obtain an approximate radius for an infinite physical value,we have set the potential for this calculated value to zero,which is followed as:r=0.012.Thus we can write the first perturbation as below:

The integral relation of the numerator in the fraction is solved in the variable change method,which ultimately leads us to the following answer as a perturbation given below:

By taking the first sentence given in these series to be calculated,the result is as follows:

Considering the amounts of given constants below and substituting them in Eq.(37)we can calculate the desired term:

5 Investigating the Isospin-Isospin Interaction Potential Impact

The same applied method in earlier section can be used here both by taking the amountsAsIandσIas below:

The Potential perturbation can be written as the following:

6 Investigating the Spin-Isospin Interaction Potential Impact

Potential perturbation of the interaction is to be calculated as follows:

In which the next given fraction is our perturbation potential:

Considering the amounts ofASIandσSIas below:

We shall have the results of entire applied interaction potentials as it is shown in Table 8.

Simply by adding the three perturbation interaction to the configuration of mentioned mass for mesons,we are able to have much more precise amount of masses using the summation of mentioned factors below:[26?27]

The results for the calculated mass in each structure and their error percentage with their experimental values are obtained in Table 9.

Table 8 Investigating the values of energy in the excited states.[22?25]

Table 9 Calculated masses for mesons,including the spin and isospin interaction potentials and their percentage error with references.

7 Results

The application of the mentioned potentials has led us to remark,that by considering the proportional fitted masses; in addition,by fetching the manual constants,we can approximate the outcomes to the empirical results obtained; furthermore,the values given as constants of this potential could be applied in physical computation,to calculate the masses of mesons with the main mentioned potential assistance.Howbeit,the significant matter in numerical calculations of interaction potentials such as spinspin,isospin-isospin,and isospin-spin is the applied perturbation category which is not practical for those masons that they have been considered as light weighted class.The same approach is principally confirmed,as the main potential utilization for the single light meson was not efficient and its error rate was far higher in analogy with the empirical references,nevertheless,in the direction of this article purpose,on calculating the mass spectrum of heavy mesons,the deployed approach was an adequate route toward reducing our numerical computational error,which by the usage of the theoretical applied trend,the mentioned error is now close to zero for even some heavy mesons.