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] aama_v1_0.gz(49 Kbytes)
Manuscript Title: Muonic atom cascade program.
Authors: V.R. Akylas, P. Vogel
Program title: MUONIC ATOM CASCADE
Catalogue identifier: AAMA_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 15(1978)291
Programming language: Fortran.
Computer: CDC6600, CDC7600.
Operating system: BKY24B-V2.
RAM: 30K words
Word size: 60
Keywords: Atomic physics, Muon negative, Cascade, Auger, Radiative, Polarization, X-rays, Multipole, Angular momentum, Monopole, Dipole, Quadrupole, Octopole, Shells, Depolarization, Penetration, Radiation physics.
Classification: 2.8, 21.1.

Nature of problem:
The radiative and Auger transition rates of a negative muon in an atomic cascade, as well as the muon depolarization are calculated, the intensities of X-rays are calculated, sorted and printed. User selected transitions and their intensities are optionally punched.

Solution method:
Using the multipole radiative and Auger transition rates and an initial distribution of angular momentum of the negative muon, the population of the final states and the intensities of the resulting muonic X-rays are calculated in a systematic way and a catalogue of X-ray intensities is printed. This program is an upgraded version of an earlier work, with the following improvements: 1) Complete monopole, dipole, quadrupole, and octupole treatment, 2) K,L and M shells are included for all multipolarities, 3) Penetration is included in all cases, 4) The muon depolarization due to the spin-orbit coupling is calculated, 5) The input and output options are flexible. The underlying principles are discussed in the long write-up and in the thesis of one of the authors.

Running time:
58.0s on CDC 7600 (including 8.1 s to compile and load) for one full cascade with mainly default options (multipoles up to gamma = 3, penetration included with cutoff at y=1) starting from principal quantum number 18. The running time will depend very strongly on the starting principal quantum number (proportional to a power between 4 and 5 thereof), on the length of the X-ray intensity catalogue, and on the other options selected (mainly on the treatment of penetration).