Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

Manuscript Title: Single particle calculations for a Woods-Saxon potential with triaxial deformations, and large Cartesian oscillator basis (new version code)
Authors: B. Mohammed-Azizi, D.E. Medjadi
Program title: Triaxial2007
Catalogue identifier: ADSK_v2_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 176(2007)634
Programming language: Fortran 77/90 double precision.
Computer: PC. Pentium 4, 2600Hz.
Operating system: Windows XP.
RAM: 256 Mb. Hard Disk - 40 Gb. Swap file - 4 Gb
Word size: 32 bits
Keywords: Nuclear physics, Energy levels, Wave functions, Schrödinger equation, Woods-Saxon potential.
PACS: 07.05.Tp, 21.60.-n, 21.60.-cs.
Classification: 17.17.

Does the new version supersede the previous version?: Yes

Nature of problem:
The Single particle energies and the single particle wave functions are calculated from one-body Hamiltonian including a central field of Woods-Saxon type, a spin-orbit interaction, and the Coulomb potential for the protons. We consider only ellipsoidal (triaxial) shapes. The deformation of the nuclear shape is fixed by the usual Bohr parameters (β, γ) .

Solution method:
The representative matrix of the Hamiltonian is built by means of the Cartesian basis of the anisotropic harmonic oscillator, and then diagonalized by a set of subroutines of the EISPACK library. Two quadrature methods of Gauss are employed to calculate respectively the integrals of the matrix elements of the Hamiltonian, and the integral defining the Coulomb potential

Reasons for new version:
More convenient handling of the eigenvectors

Summary of revisions:
One input file instead two. Reduced number of input parameters. Storage of eigenvalues and eigenvectors in memory in a very simple way which makes the code very convenient to the user.

Restrictions:
There are two restrictions for the code; The number of the major shells of the basis does not have to exceed Nmax=26. For the largest values of Nmax (~23-26), the diagonalization takes the major part of the running time, but the global run-time remains reasonable.

Running time:
(With full optimization in the project settings of the Compaq Visual Fortran on Windows XP ) With NMAX=23, for the neutrons case, and for both parities, if we need all eigenenergies and all eigenfunctions of the bound states, the running time is about 40 sec on the P4 computer at 2.6GHz. In this case, the calculation of the matrix elements takes only about 17 sec. If all unbound states are required, the runtime becomes larger.