Single particle calculations for a Woods–Saxon potential with triaxial deformations, and large Cartesian oscillator basis

We present a computer program which solves the Schrodinger equation of the stationary states for an average nuclear potential of Woods–Saxon type. In this work, we take specifically into account triaxial (i.e. ellipsoidal) nuclear surfaces. The deformation is specified by the usual Bohr parameters. The calculations are carried out in two stages. In the first, one calculates the representative matrix of the Hamiltonian in the Cartesian oscillator basis. In the second stage one diagonalizes this matrix with the help of subroutines of the Eispack library. If it is wished, one can calculate all eigenvalues, or only the part of the eigenvalues that are contained in a fixed interval defined in advance. In this latter case the eigenvectors are given conjointly. The program is very rapid, and the run-time is mainly used for the diagonalization. Thus, it is possible to use a significant number of the basis states in order to insure a best convergence of the results. Program obtainable from:CPC Program Library, Queen's University of Belfast, N. Ireland Title of program:Triaxial Catalogue number:ADSK Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADSK Licensing provisions:None Computer:PC. AMD Athlon 1000 MHz Hard disk:40 Go Ram:256 Mo Swap file:4 Go Operating system:WINDOWS XP Software used:Microsoft Visual Fortran 5.0A (with full optimizations in the settings project options) Programming language:Fortran 77/90 (double precision) Number of bits in a word:32 Number of lines:7662 No. of bytes in distributed program, including test data, etc.:174 601 Distribution format:tar gzip file Nature of the 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 ( β, γ). Method of solution: 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. Restrictions: There are two restrictions for the code: The number of the major shells of the basis does not have to excee