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] abmq_v1_0.gz(8 Kbytes)
Manuscript Title: Quasi-bound state wavefunctions.
Authors: R.M. DeVries
Catalogue identifier: ABMQ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 11(1976)249
Programming language: Fortran.
Computer: IBM 360/65.
Operating system: OS MVT RELEASE 21.8.
RAM: 175K words
Word size: 64
Keywords: Nuclear physics, Unbound, Woods-saxon, Wavefunction, Schrodinger equation, Quasi-bound, Resonant state, Direct reaction, Shell model.
Classification: 17.11, 17.19.

Nature of problem:
A program has been written which calculates the solution to the time independent Schrodinger equation for particles trapped by a Coulomb barrier but whose energy is greater than zero.

Solution method:
For fixed Woods-Saxon radius and diffuseness, the depth of the nuclear potential is iteratively found starting from a radius where the total potential energy exceeds the particle energy. Using this nuclear potential the quasi-bound wavefunction is generated from the boundary conditions at r = infinity of G(r)/r where G(r) is the irregular Coulomb wavefunction.

The program is dimensioned for a maximum of 400 mesh points and a maximum angular momentum of 10. Both these limits are easily increased if desired.

Unusual features:
Because the solution is exactly on resonance no searching on the phase shift is required. Therefore this technique works even if the resonance widths are extremely narrow (10**-10 eV for example).

Running time:
For the test case (212 Po=208Pb + alpha), 7 s using the WATFIV compiler - this includes compilation time.