Evaluation of negative energy Coulomb (Whittaker) functions

This paper describes a code for evaluating exponentially decaying negative energy Coulomb functions and their first derivatives with respect to the radial variable. The functions, which correspond to Whittaker functions of the second kind, are obtained to high accuracy for a wide range of parameters using recurrence techniques.: Program summaryTitle of program: whittaker_wCatalog identifier: ADSZProgram summary URL: http://cpc.cs.qub.ac.uk/summaries/ADSZProgram obtainable from: CPC Program Library, Queen''s University of Belfast, N. IrelandComputer: Cray T3E, Sun Ultra-5_10 sparc, Origin2000, Compaq EV67, IBM SP3, Toshiba 460CDTOperating systems under which the program has been tested: Windows NT4, Redhat Linux, SunOS 5.8Programming language used: Fortran 95Memory required to run with typical data: 500 KBNumber of bytes in distributed program, including test data, etc.: 39728Number of lines in distributed program, including test data, etc.: 2900Distribution format: tar gzip fileNature of physical problem: The closed-channel components of the asymptotic radial wave function corresponding to electron or positron scattering by atomic or molecular ions may be expressed in terms of negative energy Coulomb functions. The scattering observables are obtained from <f>S</f> or <f>T</f> matrices which in turn are obtained by matching the radial and asymptotic wavefunctions at a finite radial point. Recent large scale scattering calculations have required accurate values of the Coulomb functions at smaller <f>ρ</f> values and larger negative <f>η</f> values than previous work. The present program is designed to extend the range of parameters for which the function may be calculated.Method of solution: Recurrence relations, power series expansion, numerical quadrature.Restrictions on the complexity of the problem: The program has been tested for the parameter ranges: <f>0<ρ&les;1000</f>, <f>&z.sfnc;η&z.sfnc;&les;120</f> and <f>0&les;l&les;100</f>. These ranges may, with appropriate scaling to avoid underflow and overflow, be extended.References: A. Sunderland, C.J. Noble, P.G. Burke, V.M. Burke, Comp. Phys. Commun. 145 (2002) 311. [Copyright &y& Elsevier]