Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aecq_v1_0.tar.gz(2279 Kbytes)
Manuscript Title: QQ-onia package: a numerical solution to the Schrödinger radial equation for heavy quarkonium
Authors: Juan-Luis Domenech-Garret, Miguel-Angel Sanchis-Lozano
Program title: QQ-onia package
Catalogue identifier: AECQ_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 180(2009)768
Programming language: PAW, Physics Analysis Workstation (http://wwwasd.web.cern.ch/wwwasd/paw/).
Computer: PC/Workstation.
Operating system: Windows-XX and Unix (Linux).
Keywords: Heavy Quarkonium potential, Wave function at the origin.
PACS: 14.40.-n, 12.39.-x, 14.40.Gx.
Classification: 11.1, 11.6.

Nature of problem:
Software to solve the Schrödinger radial equation using a suitable potential V (r) for the heavy quarkonium system, allowing to perform spectroscopy. It also allows the analysis of relevant quantities of those resonances such as the square of the wave functions at the origin, their corresponding derivatives for l ≠ 0 states, typical heavy-quark velocities, and mean square radii. The package is a (userfriendly) multipurpose tool for dealing with different heavy quarkonium systems, providing a way to study the influence of a given potential on a series of relevant physical quantities, by either varying parameterized values of a well-known potential form, or by including new terms.

Solution method:
Based upon the Numerov 0(h6) Method, we perform a matching procedure to the reduced wave function at the cut point. We also perform a normalization technique for these wave functions taking into account the different domains when we use a Numerov backward-forward technique. In the case of l ≥ 2 we present a way to find the corresponding derivatives at the origin by only calculating the reduced radial wave function and first derivative. When estimating the heavy quark velocity, we introduce an additional way to compute this quantity from the virial theorem. The calculated reduced wave functions and radial wave functions at the origin are later used to obtain the heavy quarkonia nL splitting and E1 transition rates.

Additional comments:
Using Windows, to optimize the edition of the files, please, open it with MFC-WORDPAD.

Running time:
It depends on the choice of the r range, and the number of energy steps.