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] abpf_v1_0.gz(13 Kbytes)
Manuscript Title: Deformed quasiparticle states in a Woods-Saxon potential and coupled to rotational states of the core.
Authors: B. Hird, K.H. Huang
Catalogue identifier: ABPF_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 10(1975)293
Programming language: Fortran.
Computer: IBM 360/65.
Operating system: HASP.
RAM: 23K words
Word size: 32
Keywords: Nuclear physics, Collective model, Schrodinger equation, Woods-saxon, Deformation, Quasiparticle state, Hexadecapole, Energies, Expansion coefficients.
Classification: 17.20.

Subprograms used:
Cat Id Title Reference

Nature of problem:
The energies are calculated for the eigenstates of an odd nucleon outside a deformed core which is represented by a Woods-Saxon potential with quadrupole and hexadecapole deformations. Account is taken of pair excitations from the core to states with energies above that occupied by the odd nucleon. Collective excitations of the core are included, and all the rotational-particle couplings between the particle and the rotational excitations are considered, so that the complete band mixed collective spectrum of an axially symmetric permanently deformed nucleus, apart from the vibrational excitations, is generated.

Solution method:
The single particle bound eigenstates of a spherically symmetric Woods- Saxon potential form the basis for subsequent interaction matrix elements. This basis is obtained by radial Runge-Kutta integration. The two diagonalisations which follow add to the Hamiltonian the potential corresponding to a quadrupole and hexadecapole deformation with a derivative radial dependence, and the rotational excitations of the core. Between the two diagonalisations, a transformation to deformed quasi-particle states is made.

The deformations are restricted to Y02and Y04 shapes, with no provision for triaxial deformations or vibrational degrees of freedom. The maximum spin value which can be used in the programme is 27/2. Dimensioned space is provided for 32 states in the basis. Should the space be exceeded a warning statement is printed out. If more states in the basis or higher spin values are required, then the only modifications needed are to the array sizes within the program.

Running time:
The search for each Woods-Saxon eigenstate takes about 1 s, and the number of these states varies from about 20 near mass 100 to a maximum of 32. The diagonalisations and the quasiparticle trasnformation each take about 1/2 s when there are 10 states in the basis, and this time is approximately proportional to the square of the number of states.