Relativistic central-field Green's functions for the RATIP package
Abstract
From perturbation theory, Green's functions are known for providing a simple and convenient access to the (complete) spectrum of atoms and ions. Having these functions available, they may help carry out perturbation expansions to any order beyond the first one. For most realistic potentials, however, the Green's functions need to be calculated numerically since an analytic form is known only for free electrons or for their motion in a pure Coulomb field. Therefore, in order to facilitate the use of Green's functions also for atoms and ions other than the hydrogen-like ions, here we provide an extension to the RATIP program which supports the computation of relativistic (one-electron) Green's functions in an—arbitrarily given—central-field potential V(r). Different computational modes have been implemented to define these effective potentials and to generate the radial Green's functions for all bound-state energies E<0. In addition, care has been taken to provide a user-friendly component of the RATIP package by utilizing features of the Fortran 90/95 standard such as data structures, allocatable arrays, or a module-oriented design. Program summaryTitle of program:XGREENS Catalogue number: ADWM Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWM Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Licensing provisions:None Computer for which the new version has been tested: PC Pentium II, III, IV, Athlon Installations: University of Kassel (Germany) Operating systems: SuSE Linux 8.2, SuSE Linux 9.0 Program language used in the new version: ANSI standard Fortran 90/95 Memory required to execute with typical data: On a standard grid (400 nodes), one central-field Green's function requires about 50 kBytes in RAM while approximately 3 MBytes are needed if saved as two-dimensional array on some external disc space No. of bits in a word: Real variables of double- and quad-precision are used Peripheral used: Disk for input/output CPU time required to execute test data: 2 min on a 450 MHz Pentium III processor No. of lines in distributed program, including test data etc.: 82 042 No. of bytes in distributed program, including test data etc.: 814 096 Distribution format: tar.gz Nature of the physical problem: In atomic perturbation theory, Green's functions may help carry out the summation over the complete spectrum of atom and ions, including the (summation over the) bound states as well as an integration over the continuum [R.A. Swainson, G.W.F. Drake, J. Phys. A 24 (1991) 95]. Analytically, however, these functions are known only for free electrons ( V(r)≡0) and for electrons in a pure Coulomb field ( V(r)=-Z/r). For all other choices of the potential, in contrast, the Green's functions must be determined numerically. Method of solution: Relativistic Green's functions are generated for an arbitrary central-field potential V(r)=-Z(r)/r by using a piecewise linear approximation of the effective nuclear charge function Z(r) on some grid r(i=1,…,N): Z(r)=Z0i+Z1ir. Then, following McGuire's algorithm [E.J. McGuire, Phys. Rev. A 23 (1981) 186], the radial Green's functions are constructed from the (two) linear-independent solutions of the homogeneous equation [P. Morse, H. Feshbach, Methods of Theoretical Physics, McGraw-Hill, New York 1953 (Part 1, p. 825)]. In the computation of these radial functions, the Kummer and Tricomi functions [J. Spanier, B. Keith, An Atlas of Functions, Springer, New York, 1987] are used extensively. Restrictions onto the complexity of the problem: The main restrictions of the program concern the shape of the effective nuclear charge Z(r)=-rV(r), i.e. the choice of the potential, and the allowed energies. Apart from obeying the proper boundary conditions for a point-like nucleus, namely, Z(r→0)=Z>0 and Z(r→∞)=Z-N⩾0, the first derivative of the charge function Z(r) must be smaller than the (absolute value of the) energy of the Green's function, {∂Z(r)}/{∂r}<|E|. Unusual features of the program:XGREENS has been designed as a part of the RATIP package [S. Fritzsche, J. Elec. Spec. Rel. Phen. 114-116 (2001) 1155] for the calculation of relativistic atomic transition and ionization properties. In a short dialog at the beginning of the execution, the user can specify the choice of the potential as well as the energies and the symmetries of the radial Green's functions to be calculated. Apart from central-field Green's functions, of course, the Coulomb Green's function [P. Koval, S. Fritzsche, Comput. Phys. Comm. 152 (2003) 191] can also be computed by selecting a constant nuclear charge Z(r)=Z. In order to test the generated Green's functions, moreover, we compare the two lowest bound-state orbitals which are calculated from the Green's functions with those as generated separately for the given potential. Like the other components of the RATIP package, XGREENS makes careful use of the Fortran 90/95 standard.
- Publication:
-
Computer Physics Communications
- Pub Date:
- November 2005
- DOI:
- 10.1016/j.cpc.2005.01.023
- arXiv:
- arXiv:quant-ph/0409040
- Bibcode:
- 2005CoPhC.172..187K
- Keywords:
-
- 32.10;
- 32.80;
- 32.80.Fb;
- 32.80.Wr;
- Photoionization of atoms and ions;
- Other multiphoton processes;
- Quantum Physics;
- Physics - Atomic Physics
- E-Print:
- 20 pages, 1 figure