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] abpd_v1_0.gz(18 Kbytes)
Manuscript Title: PATIWEN - a code for Coulomb-nuclear interference calculations.
Authors: D.H. Feng, A.R. Barnett
Program title: PATIWEN
Catalogue identifier: ABPD_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 10(1975)401
Programming language: Fortran.
Computer: IBM 370/165.
Operating system: OSRLSE 21.
RAM: 344K words
Word size: 32
Keywords: Nuclear physics, Differential, Scattering, Elambda coulomb matrix, Electric multipole, Cross sections, Radial integrals, Continued fractions, Inelastic, Partial-wave expansion, Optical model.
Classification: 17.9.

Subprograms used:
Cat Id Title Reference
ABPC_v1_0 RCWFF CPC 8(1974)377
ABPC_v1_0 0001 RCWFF CPC 11(1976)141

Nature of problem:
The program PATIWEN calculates the differential and total cross sections for inelastic scattering of a structureless charged particle interacting via an electric multipole of order E2-E6 with nuclear system. It is also designed to accept nuclear matrix elements (up to 150 partial waves) from any conventional DWBA program and to complete a full Coulomb-nuclear interference calculation.

Solution method:
The numerical calculation closely follows the quantum-mechanical theory of nuclear and Coulomb excitations in the distorted-wave Born approximation. The projectile is assumed to be a point charge. The radial matrix elements are integrated by the gaussian quadrature technique. This represents an attractive procedure as the Coulomb functions needed for the integrand can be calculated independently at any radius using subroutine RCWFN. An adaptation of RCWFN which requires less array space can be used.

The calculation is non-relativistic and is dimensioned for 600 partial waves; more can be included if necessary. There is no inherent restriction to the scattering of like charges, although only this case has been tested.

Running time:
The running time depends directly on the integration range and the number of partial waves, i.e. on rho L. For a typical heavy-ion E2 case with a rho = 800 (r=250fm) and L=500 the running time is about 3 min, while for the E3 test case rho = 400 (r=250fm) and L=145 the time is 50 s. These times are central processor times for the IBM 370/165 computer; they could be drastically reduced by exploiting various recurrence relations among the Coulmob matrix elements.