Programs in Physics & Physical Chemistry
|[Licence| Download | New Version Template] aadx_v1_0.gz(16 Kbytes)|
|Manuscript Title: Generalized fractional parentage coefficients for shell-model calculations.|
|Authors: L.D. Skouras, S. Kossionides|
|Program title: GFPC1|
|Catalogue identifier: AADX_v1_0|
Distribution format: gz
|Journal reference: Comput. Phys. Commun. 39(1986)197|
|Programming language: Fortran.|
|Operating system: PRIMOS VERSION 18.3 UPWARDS.|
|RAM: 400K words|
|Keywords: Nuclear physics, Antisymmetry, Angular momentum, Isospin, Fractional parentage Coefficients.|
Nature of problem:
In a space defined by several single-particle orbitals anti-symmetric states for any n-fermion system can be expanded in terms of states which are antisymmetric only with respect to the first m and the last (n-m) particles. These (m, n-m) expansion coefficients, the Generalized Fractional Parentage Coefficients (gfpc), can then be used to determine easily the matrix elements of any m-body operator. The states are classified according to C which denotes the distribution of the n particles among the orbitals, their total angular momentum J and isospin T and an index mu that distinguishes the orthogonal states that have the same C, J, T values. Program GFPC1 constructs the coefficients for m=1, the (1,n-1) gfpc. It shares its long write-up with program GFPCM which calculates the (m, n-m) gfpc for 2 < m < n.
The (1, n-1) gfpc are calculated with recursive formulae. The Schmid method is used to select for each C, J, T the orthogonal states.
The complexity is restricted only by the available Disk space for storage of the gfpc.
The program is written in standard FORTRAN-77 and can easily be transported. All machine dependent characteristics (e.g. word length, file units, array dimensions) are set with PARAMETER statements and can easily be adapted. The program will produce the complete set of gfpc in one or several consecutive runs. Upon entry it will determine the point to which the gfpc-file has been completed and continue to the point determined by the input data. For all error conditions, error messages will be printed.
It depends on the complexity of the problem and averages to about 10 coefficients per second on the PRIME-750.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|