Programs in Physics & Physical Chemistry
|[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.|
|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.
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.
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.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|