QQ¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q\bar{Q}$$\end{document} (Q∈{b,c}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q\in \{b, c\}$$\end{document}) spectroscopy using the Cornell potential

The mass spectra and decay properties of heavy quarkonia are computed in nonrelativistic quark-antiquark Cornell potential model. We have employed the numerical solution of Schrödinger equation to obtain their mass spectra using only four parameters namely quark mass (mc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_c$$\end{document}, mb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_b$$\end{document}) and confinement strength (Acc¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{c\bar{c}}$$\end{document}, Abb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{b\bar{b}}$$\end{document}). The spin hyperfine, spin-orbit and tensor components of the one gluon exchange interaction are computed perturbatively to determine the mass spectra of excited S, P, D and F states. Digamma, digluon and dilepton decays of these mesons are computed using the model parameters and numerical wave functions. The predicted spectroscopy and decay properties for quarkonia are found to be consistent with available data from experiments, lattice QCD and other theoretical approaches. We also compute mass spectra and life time of the Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} meson without additional parameters. The computed electromagnetic transition widths of heavy quarkonia and Bc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_c$$\end{document} mesons are in tune with available experimental data and other theoretical approaches.


Introduction
Mesonic bound states having both heavy quark and antiquark (cc, bb and cb) are among the best tools for understanding the quantum chromodynamics. Many experimental groups such as CLEO, LEP, CDF, D0 and NA50 have provided data and BABAR, Belle, CLEO-III, ATLAS, CMS and LHCb are producing and expected to produce more precise data in upcoming experiments. Comprehensive reviews on the status of experimental heavy quarkonium physics are found in literature [1][2][3][4][5][6].
In 90's, the nonrelativistic potential models predicted not only the ground state mass of the tightly bound state of c andb in the range of 6.2-6.3 GeV [81,82] but also predicted to have very rich spectroscopy. In 1998, CDF collaboration [83] reported B c mesons in pp collisions at √ s = 1.8 TeV and was later confirmed by D0 [84] and LHCb [85] collaborations. The LHCb collaboration has also made the most precise measurement of the life time of B c mesons [86]. The first excited state is also reported by ATLAS Collaborations [87] in pp collisions with significance of 5.2σ .
It is important to show that any given potential model should be able to compute mass spectra and decay properties of B c meson using parameters fitted for heavy quarkonia. Attempts in this direction have been made in relativistic quark model based on quasi-potential along with one loop radiative correction [27], quasistatic and confinement QCD potential with confinement parameters along with quark masses [88] and rainbow-ladder approximation of Dyson-Schwinger and Bethe-Salpeter equations [67].
The interaction potential for mesonic states is difficult to derive for full range of quark antiquark separation from first principles of QCD. So most forms of QCD inspired potential would result in uncertainties in the computation of spectroscopic properties particularly in the intermediate range. Different potential models may produce similar mass spectra matching with experimental observations but they may not be in mutual agreement when it comes to decay properties like decay constants, leptonic decays or radiative transitions. Moreover, the mesonic states are identified with masses along with certain decay channels, therefore the test for any successful theoretical model is to reproduce the mass spectrum along with decay properties. Relativistic as well as nonrelativistic potential models have successfully predicted the spectroscopy but they are found to differ in computation of the decay properties [22,[47][48][49][50][51]55,[78][79][80]. In this article, we employ nonrelativistic potential with one gluon exchange (essentially Coulomb like) plus linear confinement (Cornell potential) as this form of the potential is also supported by LQCD [89][90][91]. We solve the Schrödinger equation numerically for the potential to get the spectroscopy of the quarkonia. We first compute the mass spectra of charmonia and bottomonia states to determine quark masses and confinement strengths after fitting the spin-averaged ground state masses with experimental data of respective mesons. Using the potential parameters and numerical wave function, we compute the decay properties such as leptonic decay constants, digamma, dilepton, digluon decay width using the Van-Royen Weiskopf formula. These parameters are then used to compute the mass spectra and life-time of B c meson. We also compute the electromagnetic (E1 and M1) transition widths of heavy quarkonia and B c mesons.

Methodology
Bound state of two body system within relativistic quantum field is described in Bethe-Salpeter formalism. However, the Bethe-Salpeter equation is solved only in the ladder approximations. Also, Bethe-Salpeter approach in harmonic confinement is successful in low flavor sectors [92,93]. Therefore the alternative treatment for the heavy bound state is nonrelativistic. Significantly low momenta of quark and antiquark compared to mass of quark-antiquark system m Q,Q QC D ∼ |p| also constitutes the basis of the nonrelativistic treatment for the heavy quarkonium spectroscopy. Here, for the study of heavy bound state of mesons such as cc, cb and bb, the nonrelativistic Hamiltonian is given by where where m Q and mQ are the masses of quark and antiquark respectively, p is the relative momentum of the each quark and V Cornell (r ) is the quark-antiquark potential of the type coulomb plus linear confinement (Cornell potential) given by Here, 1/r term is analogous to the Coulomb type interaction corresponding to the potential induced between quark and antiquark through one gluon exchange that dominates at small distances. The second term is the confinement part of the potential with the confinement strength A as the model parameter. The confinement term becomes dominant at the large distances. α s is a strong running coupling constant and can be computed as where n f is the number of flavors, μ is renormalization scale related to the constituent quark masses as μ = 2m Q mQ/(m Q + mQ) and is a QCD scale which is taken as 0.15 GeV by fixing α s = 0.1185 [7] at the Z -boson mass. The confinement strengths with respective quark masses are fine tuned to reproduce the experimental spin averaged ground state masses of both cc and bb mesons and they are given in Table 1. We compute the masses of radially and orbitally excited states without any additional parameters. Similar work has been done by [47,51,52] and they have considered different values of confinement strengths for different potential indices. The Cornell potential has been shown to be independently successful in computing the spectroscopy of ψ and ϒ families. In this article, we compute the mass spectra of the ψ and ϒ families along with B c meson with minimum number of parameters.
Using the parameters defined in Table 1, we compute the spin averaged masses of quarkonia. In order to compute masses of different n m L J states according to different J PC values, we use the spin dependent part of one gluon exchange potential (OGEP) V S D (r ) perturbatively. The OGEP includes spin-spin, spin-orbit and tensor terms given by [20,22,59,68] The spin-spin interaction term gives the hyper-fine splitting while spin-orbit and tensor terms gives the fine structure of the quarkonium states. The coefficients of spin dependent terms of the Eq. (5) can be written as [20] Where V V (r ) and V S (r ) correspond to the vector and scalar part of the Cornell potential in Eq. (3) respectively. Using all the parameters defined above, the Schrödinger equation is numerically solved using Mathematica notebook utilizing the Runge-Kutta method [94]. It is generally believed that the charmonia need to be treated relativistically due to their lighter masses, but we note here that the computed wave functions of charmonia using relativistic as well as nonrelativistic approaches do not show significant difference [33]. So we choose to compute the charmonium mass spectra nonrelativistically in present study. The computed mass spectra of heavy quarkonia and B c mesons are listed in Tables 2, 3, 4, 5, 6 and 7.

Decay properties
The mass spectra of the hadronic states are experimentally determined through detection of energy and momenta of daughter particles in various decay channels. Generally, most phenomenological approaches obtain their model parameters like quark masses and confinement/Coulomb strength by fitting with the experimental ground states. So it becomes necessary for any phenomenological model to validate their fitted parameters through proper evaluation of various decay rates in general and annihilation rates in particular. In the nonrelativistic limit, the decay properties are dependent on the wave function. In this section, we test our parameters and wave functions to determine various annihilation widths and electromagnetic transitions.

Leptonic decay constants
The leptonic decay constants of heavy quarkonia play very important role in understanding the weak decays. The matrix elements for leptonic decay constants of pseudoscalar and vector mesons are given by where k is the momentum of pseudoscalar meson, * μ is the polarization vector of meson. In the nonrelativistic limit, the decay constants of pseudoscalar and vector mesons are given by Van Royen-Weiskopf formula [96] Here the QCD correction factorC 2 (α S ) [97,98] With δ P = 2 and δ V = 8/3. Using the above relations, we compute the leptonic decay constants f p and f v for charmonia, bottomonia and B c mesons. The results are listed in Tables  8, 9, 10, 11, 12 and 13 in comparison with other models including LQCD.

Annihilation widths of heavy quarkonia
Digamma, digluon and dilepton annihilation decay widths of heavy quarkonia are very important in understanding the dynamics of heavy quarks within the mesons. The measurement of digamma decay widths provides the information regarding the internal structure of meson. The decay η c → γ γ , χ c0,2 → γ γ was reported by CLEO-c [103], BABAR [104] and then BESIII [105] collaboration have reported high accuracy data. LQCD is found to underestimate the decay widths of η c → γ γ and χ c0 → γ γ when compared to experimental data [106,107]. Other approaches to attempt computation of annihilation rates of heavy quarkonia include NRQCD [108][109][110][111][112], relativistic quark model [31,32], effective Lagrangian [113,114] and next-to-next-to leading order QCD correction to χ c0,2 → γ γ in the framework of nonrelativistic QCD factorization [115]. The meson decaying into digamma suggests that the spin can never be one [116,117]. Corresponding digamma decay width of a pseudoscalar meson in nonrelativistic limit is given by Van Royen-Weiskopf formula [96,118] where the bracketed quantities are QCD next-to-leading order radiative corrections [118,119]. Digluon annihilation of quarkonia is not directly observed in detectors as digluonic state decays into various hadronic states making it a bit complex to compute digluon annihilation widths from nonrelativistic approximations derived from first principles. The digluon decay width of pseudoscalar meson along with the QCD leading order radiative correction is given by [113,[118][119][120] Here, the coefficients in the bracket have values of C Q = 4.8, C 0Q = 9.5, C 2Q = −2.2 for the charm quark and C Q = 4.4, C 0Q = 10.0, C 2Q = −0.1 for the bottom quark [118]. The vector mesons have quantum numbers 1 −− and can annihilate into dilepton. The dileptonic decay of vector meson along with one loop QCD radiative correction is given by [96,118] Here, α e is the electromagnetic coupling constant, α s is the strong running coupling constant in Eq. (4) and e Q is the charge of heavy quark in terms of electron charge. In above relations, |R ns P/V (0)| corresponds to the wave function of S-wave at origin for pseudoscalar and vector mesons while |R n P (0)| is the derivative of P-wave function at origin. The annihilation rates of heavy quarkonia are listed in Tables 14,  15, 16, 17, 18 and 19.

Electromagnetic transition widths
The electromagnetic transitions can be determined broadly in terms of electric and magnetic multipole expansions and their study can help in understanding the non-perturbative regime of QCD. We consider the leading order terms i.e. electric (E1) and magnetic (M1) dipoles with selection rules L = ±1 and S = 0 for the E1 transitions while L = 0 and S = ±1 for M1 transitions. We now employ the numerical where, mean charge content e Q of the QQ system, magnetic dipole moment μ and photon energy ω are given by and  respectively. Also the symmetric statistical factor is given by The matrix element |M i f | for E1 and M1 transition can be written as and The electromagnetic transition widths are listed in Tables 20,  21 In the spectator model [126], the total decay width of B c meson can be broadly classified into three classes. (i) Decay of b quark considering c quark as a spectator, (ii) Decay of c quark considering b quark as a spectator and (iii) Annihilation channel B c → + ν . The total width is given by (28) In the calculations of total width we have not considered the interference among them as all these decays lead to different channel. In the spectator approximation, the inclusive decay width of b and c quark is given by where C q = 3|V cs | for D s mesons and m q is the mass of heaviest fermions. V cs and V cb are the CKM matrices and we have taken the value of CKM matrices from the PDG. G f is the Fermi coupling constant. Here we have used the

Numerical results and discussion
Having determined the confinement strengths and quark masses, we are now in position to present our numerical results. We first compute the mass spectra of heavy quarkonia and B c meson. In most of the potential model computations, the confinement strength is fixed by experimental ground state masses for cc, bb and cb independently. We observe here that the confinement strength A for B c meson is arithmetic mean of those for cc and bb which discards the         need to introduce additional confinement strength parameter for computation of B c spectra. Similar approach has been used earlier within QCD potential model [88]. Using model parameters and numerical wave function we compute the various decay properties of heavy quarkonia and B c mesons namely leptonic decay constants, annihilation widths and electromagnetic transitions. In Tables 2 and 3, we present our result for charmonium mass spectra. We compare our results with PDG data [7], lattice QCD [17] data, relativistic quark model [27], nonrelativistic quark model [65,68], QCD relativistic functional approach [67], relativistic potential model [39] and nonrelativistic potential models [59,70,73,76]. Our results for S-wave are in excellent agreement with the experimental data [7]. We determine the mass difference for S-wave charmonia i.e. M J/ψ − M η c = 105 MeV and M ψ(2S) −M η c (2S) = 79 MeV while that from experimental data are 113 and 47 MeV respectively [7]. Our results for P-waves are also consistent with the PDG data [7] as well    as LQCD [17] with less than 2% deviation. Since experimental/LQCD results are not available for P-wave charmonia beyond n = 2 states, we compare our results with the relativistic quark model [27] and it is also observed to have 1-2 % deviation throughout the spectra. For charmonia, only 1P states are available and for 2P only one state is available namely χ c2 . Our results for 1P and 2P states are also satisfactory. We also list the mass spectra of D and F wave and find it to be consistent with the theoretical predictions. Overall, computed charmonium spectra is consistent with PDG and other theoretical models.
In Tables 4 and 5, we compare our results of bottomonium spectra with PDG data [7], relativistic quark model [27,64], nonrelativistic quark model [66], QCD relativistic functional approach [67], relativistic potential model [40], nonrelativistic potential models [73,77] and covariant constituent quark model [69]. Similarly for S-wave bottomonia, up to n = 3 vector states are known experimentally and for pseudoscalar states, only n = 1 and 2 are available. Our results for ϒ(1S) and ϒ(3S) are in good agreement with the PDG data while for ϒ(2S), ϒ(4S) and ϒ(5S), slight deviation (within 1%) is observed. Our results for η b (1S) and η b (3S) also match well with less than 0.5% deviation. We obtain M ϒ(1S) − M η b = 35 MeV and for M ϒ(2S) − M η b (2S) = 24 MeV against the PDG data of 62 and 24 MeV respectively. For P-wave, 1P and 2P states are reported and for 3P, only χ b1 is reported. Our results for 1P bottomonia deviate by 0.3% from the experimental results but for 2P, they are quite satisfactory and devi- ating by 0.2% only from the PDG data. Our result for ϒ(1D) also agrees well with the experimental data with 0.8% deviation. The F-wave mass spectra is also in good agreement with the theoretical predictions. Looking at the comparison with PDG data Ref. [7] and relativistic quark model Ref. [27], present quarkonium mass spectra deviate less than 2% for charmonia and less than 1% for bottomonia. We now employ the quark masses and confinement strengths used for computing mass spectra of quarkonia to predict the spectroscopy of B c mesons without introducing any additional parameter. Our results are tabulated in Tables  6 and 7. For B c mesons, only 0 −+ states are experimentally observed for n = 1 and 2 and our results are in very good agreement with the experimental results with less than 0.3 % error.
We note here that the masses of orbitally excited states (especially n = 1 states) of charmonia are systematically lower than the other models and experimental data. This tendency decreases as one moves to higher n states. Absence of similar trend in case of B c and bottomonia systems suggests that relativistic treatment might improve the results in lower energy regime of charmonia.
Using the mass spectra of heavy quarkonia and B c meson, we plot the Regge trajectories in (J, M 2 ) and (n r , M 2 ) planes where n r = n − 1. We use the following relations [27]   where α, β are slopes and α 0 , β 0 are the intercepts that can be computed using the methods given in Ref. [27]. In Figs. 1, 2 and 3, we plot the Regge trajectories. Regge trajectories from present approach and relativistic quark model [27] show similar trend i.e. for charmonium spectra, the computed mass square fits very well to a linear trajectory and found to be almost parallel and equidistant in both the planes. Also, for bottomonia and B c mesons, we observe the nonlinearity in the parent trajectories. The nonlinearity increases as we go from cb to bb mesons indicating increasing contribution from the inter-quark interaction over confinement.
According to the first principles of QCD, while the onegluon-exchange interaction gives rise to employment of Coulomb potential with a strength proportional to the strong   running coupling constant at very short distances, nonperturbative effect like confinement becomes prominent at larger distances. Charmonium belongs to neither purely nonrelativistic nor the relativistic regime where chiral symmetry breaking is more significant from physics point of view. Though Lattice QCD calculations in the quenched approximation have suggested a linearly increasing potential in the confinement range [8][9][10][11][12][13][14][15][16][17][18], a specific form of interaction potential in the full range is not yet known. At short distances relativistic effects are more important as they give rise to quark-antiquark pairs from the vacuum that in turn affect the nonrelativistic Coulomb interaction in the presence of sea quarks. The mass spectra of quarkonia is not sensitive to these relativistic effects at short distances. However, the decay properties show significant difference with inclusion of relativistic corrections. We have used the most accepted available correction terms for computation of decay properties [113,[118][119][120] that improves the results significantly in most cases. Using the potential parameters and numerical wave function, we compute the various decay properties of heavy quarkonia. We first compute the leptonic decay constants of pseudoscalar and vector mesons and our numerical results are tabulated in Tables 8, 9, 10, 11, 12 and 13. For the case of charmonia, our results are higher than those using LQCD and QCDSR [100]. In order to overcome this discrepancy, we include the QCD correction factors given in Ref. [97] and the results are tabulated in Tables 8 and 9. After introducing the correction factors our results match with PDG, LQCD and QCDSR [100] along with other theoretical models. We also compute the decay constants for excited S-wave charmonia and we found that our results are consistent with the other theoretical predictions. We also compute the decay constants of bottomonia and B c mesons. In this case, our results match with other theoretical predictions without incorporating the relativistic corrections. In the case of vector decay constants of bottomonia, our results are very close to experimental results as well as those obtained in LQCD Ref. [102]. For the decay constants of B c mesons, we compare our results with nonrelativistic potential models [52,95].
Next we compute the digamma, digluon and dilepton decay widths using the relations Eqs. (13)- (16). Where the bracketed quantities are the first order radiative corrections to the decay widths. We compare our results with the available experimental results. We also compare our results with the theoretical models such as screened potential model [76,77], Martin-like potential model [73], relativistic quark model (RQM) [31,32], heavy quark spin symmetry [114], relativistic Salpeter model [121] and other theoretical data. Tables 14 and 15 we present our results for digamma decay widths for charmonia and bottomonia. Our results for (η c → γ γ ) and (η c (2S) → γ γ ) are higher than the experimental results. Experimental observation of the two photon decays of pseudoscalar states are considered as an important probe for identification of flavour as well as internal structure of mesons. The first order radiative correction [bracketed terms in Eq. (13)] was utilized to incorporate the difference and it is observed that our results along with the correction match with the experimental results [7]. We also compute the digamma decay width of excited charmonia.
Our results for P-wave charmonia are higher than that of screened potential model [76] and relativistic quark model [32]. Our results for (η b → γ γ ) match quite well with the experimental data while computed (η b (2S) → γ γ ) value is overestimated when compared with the PDG data. For the excited state of S-wave bottomonia, our results fall in between those obtained in screened potential model [77] and relativistic quark model with linear confinement [64]. The scenario is similar with P-wave bottomonia and charmonia.
Di-gluon decay has substantial contribution to hadronic decay of quarkonia below cc and bb threshold. In Tables  16 and 17 we represent our results for digluon decay width of charmonia and bottomonia respectively. Our results for (η c → gg) match perfectly with the PDG data [7] but in the case of (η c (2S) → gg) our result is higher than the PDG data. We also compare the results obtained with that of the relativistic Salpeter method [121] and an approximate potential model [70]. It is seen from Table 16 that the relativistic corrections provide better results in case of P-wave charmonia where as that for bottomonia are underestimated in present calculations when compared to relativistic QCD potential model [122] and power potential model [47]. As the experimental data of digluon annihilation of bottomonia are not available, the validity of either of the approaches can be validated only after observations in forthcoming experiments.
We present the result of dilepton decay widths in the Tables  18 and 19 and it is observed that our results matches with the PDG data [7] upto n = 3 for both charmonia and bottomonia. The contribution of the correction factor is more significant in the excited states with compared to that in the ground states of the quarkonia, indicating different dynamics in the intermediate quark-antiquark distance. Our results are also in good accordance with the other theoretical models.
We present our results of E1 transitions in Tables 20, 21 and 22 in comparison with theoretical attempts such as relativistic potential model [39], quark model [30], nonrelativistic screened potential model [66,76,77]. We also compare our results of charmonia transitions with available experimental results. Our result for (ψ(2S) → χ cJ (1P) + γ ) is in good agreement with the experimental result for J = 0 but our results for J = 1, 2 are higher than the PDG data. Our results also agree well for the transition (χ c2 (1P) → J/ψ + γ ). We also satisfy the experimental constraints for the transition (1 3 D 1 → χ cJ + γ ) for J = 0, 1, 2. Our results share the same range with the results computed in other theoretical models. The E1 transitions of bottomonia agree fairly well except for the channel (ϒ(3S) → χ bJ (3P)), where our results are higher than the experimental results. The comparison of our results of E1 transitions in B c mesons with relativistic quark model [30,63] and power potential model [46] are found to be in good agreement. In Tables 23, 24 and 25, we present our results of M1 transitions and also compared with relativistic potential model [39], quark model [30,64], nonrelativistic screened potential model [65,66], power potential [46] as well as with available experimental results. Our results of (nψ → n η c + γ ) are in very good agreement with the PDG data as well with the other theoretical predictions. Computed M1 transitions in B c mesons are also within the results obtained from theoretical predictions. The computed M1 transition of bottomonia are found to be higher than the PDG data and also theoretical predictions.

Conclusion
In this article, we have reported a comprehensive study of heavy quarkonia in the framework of nonrelativistic potential model considering linear confinement with least number of free model parameters such as confinement strength and quark mass. They are fine tuned to obtain the corresponding spin averaged ground state masses of charmonia and bottomonia determined from experimental data. The parameters are then used to predict the masses of excited states. In order to compute mass spectra of orbitally excited states, we incorporate contributions from the spin dependent part of confined one gluon exchange potential perturbatively.
Our results are found to be consistent with available PDG data, LQCD, relativistic quark model and other theoretical potential models. We also compute the digamma, digluon and dilepton decay widths of heavy quarkonia using nonrelativistic Van-Royen Weiskopf formula. The first order radiative corrections in calculation of these decays provide satisfactory results for the charmonia while no such correction is needed in case of bottomonia for being purely nonrelativistic system. We employ our parameters in computation of B c spectroscopy employing the quark masses and mean value of confinement strength of charmonia and bottomonia and our results are also consistent with the PDG data. We also compute the weak decays of B c mesons and the computed life time is also consistent with the PDG data. It is interesting to note here that despite having a c quark, the nonrelativistic calculation of B c spectroscopy is in very good agreement with experimental and other theoretical models.