Comparison Between Two Numerical Schemes to Study the Spectra of Charmed Quarkonium

TL;DR

This paper compares two numerical schemes for solving the radial Schrödinger equation to study the spectra of charmed quarkonium, evaluating their accuracy, convergence, and agreement with experimental data.

Contribution

It introduces and compares two novel numerical methods for analyzing heavy quarkonium spectra within nonrelativistic and QCD-inspired models.

Findings

One scheme shows higher precision in spectral calculations.

Both methods reliably model charmonium bound states.

Convergence behavior varies between the two techniques.

Abstract

Two numerical methods are developed to reduce the solution of the radial Schr\"odinger equation for proposed heavy quark-antiquark interactions, into the solution of the eigenvalue problem for the infinite system of tridiagonal matrices. Our perspective is a numerical approach relies on finding the proper numerical method to investigate the static properties of heavy quarkonia-mesons, such as spectrum, radius ... etc., with implantation of both the nonrelativistic quark model and the ingredients of the quantum chromodynamics (QCD) theory. The application of these proposed schemes resulted in mass spectra of charmed-quarkonium (charmonium) multiplets, which are compared with the experimental published profiles of Particle Data Group (PDG). Besides, the normalized radial wave-functions of the charmonium various bound states are represented. The convergence of each numerical recipe versus the iteration number N and the radial distance is investigated through this work. Although it was observed that our numerical treatments are reliable to study charmed Quarkonium bound states profile, we found that one of these proposed techniques is favored over the other in terms of high precision comparisons with experiments and convergence analysis.