We computed the eigenvalues and eigenfunctions via the Cornell potential, approximately, using Numerov approach. The experimental data for the quantum states ³ S ₁ were utilized to fit the mass spectra. Results were compared with observed data and the results obtained from other potential functions in the existing literature for the heavy ( c c ¯ , b b ¯ and b c ¯ ) as well as heavy-light (B, Bs, D and Ds) meson states. The yielded meson spectra based on this potential are matching well with the available experimental data.

An analytical study of a source of Beauty mesons (B mesons) by the Strong decay of $$\Upsilon$$Υ (11020) is performed to obtain various modes of B mesons. The non-relativistic quark model framework is used to predict the masses of the bottomonium mesons. The strong decay of the $$\Upsilon$$Υ (11020) state is calculated utilizing Quark Pair Creation model. In addition, the branching ratio of the strong decay has been calculated for each decay width and using them to estimate the number of B-meson pairs for each strong decay width. The results show reasonable agreement with recent Particle Data Group data, in particular with the total strong decay width of $$\Upsilon$$Υ (11020). The $$\Upsilon$$Υ (11020) particle is expected to be a good candidate for the conventional $$\Upsilon$$Υ (6S) state and provides both types of strange and non-strange B mesons. We found new conclusions for the strong decay widths of the $$\Upsilon$$Υ (11020) in the relativistic phase space. The strong width of the $$B^{*}B^{*}$$B ∗ B ∗ channel is predicted to be the highest width by the ratio ( $$\cong 51\%$$≅ 51 % ) concerning higher vector $$\Upsilon$$Υ (11020) State. In addition, the specific $$B_{s}B_{s}^{*}$$B s B s ∗ channel is expected to be the highest width relative to higher $$\Upsilon$$Υ resonances.

A prediction study for the not yet discovered $$c\overline{s}$$c s ¯ and $$b\overline{s}$$b s ¯ meson states of the $$1^{3}D_1$$1 3 D 1 quantum numbers is conducted, denoted as $$D_s$$D s and $$B_s$$B s , respectively. By employing our method, both states of ground and excited vector mesons are considered, which predominantly correspond to the ground state vector mesons with quantum numbers $$n~ ^{2S+1}L_{J} = 1 ~ ^{3}S_1$$n 2 S + 1 L J = 1 3 S 1 and to the orbitally excited vector mesons with quantum numbers $$n~ ^{2\,S+1}L_{J} = 1~ ^{3}D_1$$n 2 S + 1 L J = 1 3 D 1 . The considered method depends on the difference in the masses for the relevant states from orbitally and ground state vector mesons. Theoretical calculations are in good agreements with newly published data.