Programs in Physics & Physical Chemistry
|[Licence| Download | New Version Template] abjm_v1_0.gz(14 Kbytes)|
|Manuscript Title: Perturbation functions: PAC probe nuclei, I=2, 5/2, and 3.|
|Authors: G.L. Catchen|
|Program title: ACOEFF|
|Catalogue identifier: ABJM_v1_0|
Distribution format: gz
|Journal reference: Comput. Phys. Commun. 55(1989)85|
|Programming language: Fortran.|
|Computer: IBM 3090-400E.|
|Operating system: OS/VS2.|
|RAM: 2217K words|
|Word size: 32|
|Keywords: Nuclear physics, Spectroscopy, Perturbation function, Attenuation coefficient, Perturbed angular Correlation, Level scheme.|
Nature of problem:
In Perturbed Angular Correlation (PAC) Spectroscopy, static nuclear- electric-quadrupole interactions are described by perturbation functions (attenuation coefficients). For a quadrupole interaction with a non- axially-symmetric electric field gradient, the perturbation function calculation involves a summation over projection quantum numbers, Wigner 3-j coefficients, and the associated unitary matrix elements.
The calculation can be performed for spins, I=2, 5/2, and 3, for the case of a polycrystalline source of PAC probe nuclei. The Wigner 3-j coefficients that correspond to all of the m-to-m' transitions are generated and stored as a matrix. For each value of the assymetry parameter, eta, the program generates the unitary matrix and performs the summation. The program calculates and prints the associated m- state energy differences, a-coefficients, and perturbation functions, Gkk(t) for k=2 and 4 (I=2 and 5/2) and k=2, 4, and 6 (I=3).
Only sources of probe nuclei in randomly-oriented crystals, i.e., polycrystalline sources, are considered. The perturbation functions are calculated in 101 steps of (Omega0 *t/2*n*pi) in which n is usually chosen to be 1 or 2.
The program contains all of the Wigner 3-j coefficients needed for the calculation. From the corresponding closed-forms expressions, the program generates the m-state energy differences and the eigenvectors that compose the unitary transformation.
14s (4s, VS compiler) (for I=3, perturbation functions are calculated for 11 values of eta) and 40s (9s, VS compiler) (for I=5/2, S- coefficients are calculated for 101 values of eta and no perturbation functions are calculated).
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|