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] abnb_v1_0.gz(37 Kbytes)
Manuscript Title: Static and dynamic muonic-atom codes MUON and RURP.
Authors: G.A. Rinker
Program title: MUON
Catalogue identifier: ABNB_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 16(1979)221
Programming language: Fortran.
Computer: CDC 6400, CDC 6600, CDC 7600.
Operating system: ANY BATCH PROCESSOR.
RAM: 25K words
Word size: 60
Keywords: Nuclear physics, Atomic physics, Muonic atoms, Energy levels, Qed corrections.
Classification: 2.8.

Nature of problem:
Muonic-atom energy levels and wave functions are calculated. The results are corrected for nearly all important static effects, including finite nuclear size as produced by a user-specified nuclear charge distribution.

Solution method:
Fourth-order Runge-Kutta integration of the radial Dirac equations, with the secant method determining the eigenvalues. Corrections are either included as potentials in the Dirac equations or computed as first-order perturbations.

Restrictions:
The code assumes spherical symmetry, so that effects of nuclear deformation may be treated only as static first-order perturbations. Dynamic muon-electron and muon-nuclear effects are neglected (the latter are calculated by RURP).

Unusual features:
The main arrays are allocated semi-dynamically in order to keep core requirements to a minimum. Although the code as supplied is set up for CDC computers, previous versions have been run successfully on IBM 370/ 168 using the H-compiler with the AUTODOUBLE option. Use of the code on low-precision machines is recommended only if extensive double precision is implemented.

Running time:
The code and test case as supplied are compiled, loaded, and executed in 59 s on a CDC-6600, of which 21 s is compilation time (FTN compiler with optimization level 2).