An improved version of ISICS: a program for calculating K-, L- and M-shell cross sections from PWBA and ECPSSR theory using a personal computer

The C program, ISICS [Z. Liu, S.J. Cipolla, Comput. Phys. Comm. 97 (1996) 315-330], which calculates ionization and X-ray production cross-sections using PWBA and ECPSSR theory, has been enhanced to include new options, correct some minor flaws, and to make the program more versatile.

Program summary

Title of program: ISICSCatalog identifier: ADDS_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADDS_v2_0Program available from: CPC Program Library, Queen's University of Belfast, N. IrelandOperating system under which the program has been tested: WINDOWS XPProgram language used: CComputer: 80486 or higher-level PCsNo. of lines in distributed program, including test data, etc.: 5343No. of bytes in distributed program, including test data, etc.: 151 838Distribution format: tar.gzCatalogue identifier of previous version: ADDSJournal reference of previous version: Comput. Phys. Comm. 97 (1996) 315-330Does the new version supersede the previous version: YesNature of the physical problem: Ionization and X-ray production cross-section calculations for ion-atom collisions.Reasons for new version: Increased functionality and new options.Summary of revisions: Option for the united-atom approximation for binding-energy correction; easier inputting of updated atomic parameters; extension of projectile energy down to eV range; accounting for DHS wave function in K-shell ionization; other miscellaneous changes.Method of solution: Numerical integration of form factor using a logarithmic transform and Gaussian quadrature, plus exact integration limits.Restrictions on the complexity of the problem: The consumed CPU time increases with the atomic shell (K, L, M), but execution is still very fast.Typical running time: No change from previous version.Unusual features of the program: No     

An improved version of ISICS: A program for calculating K-, L- and M-shell cross sections from PWBA and ECPSSR theory using a personal computer

The formatting of the M-shell atomic parameters imbedded in file XCSC.H in ISICS has been corrected. The problem only affected cross section calculations for Uranium and heavier elements. The corrected version of ISICS has been re-compiled and is now available.

New version program summary

Program title: ISICSCatalogue identifier: ADDS_v3_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADDS_v3_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 4645No. of bytes in distributed program, including test data, etc.: 106 731Distribution format: tar.gzProgramming language: C++Computer: 80486 or higher-level PCsOperating system: WINDOWS 98 through WINDOWS XPClassification: 16.7Does the new version supersede the previous version?: YesNature of problem: Ionization and X-ray production cross section calculations for ion-atom collisions.Solution method: Numerical integration of form factor using a logarithmic transform and Gaussian quadrature, plus exact integration limits.Reasons for new version: The formatting of the M-shell atomic parameters involving cross section calculations for Uranium and heavier elements needed to be corrected.Summary of revisions: The affected file XCSC.H in ISICS has been corrected and ISICS has been recompiled.Restrictions: The consumed CPU time increases with the atomic shell (K, L, M), but execution is still very fast.Running time: This depends on which shell and the number of different energies to be used in the calculation. For example, to calculate K-shell cross sections for protons striking carbon for 19 different proton energies it took less than 10 s; to calculate M-shell cross sections for protons on gold for 21 proton energies it took 4.2 min.     

JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices

A new software code for computing selected eigenvalues and associated eigenvectors of a real symmetric matrix is described. The eigenvalues are either the smallest or those closest to some specified target, which may be in the interior of the spectrum. The underlying algorithm combines the Jacobi-Davidson method with efficient multilevel incomplete LU (ILU) preconditioning. Key features are modest memory requirements and robust convergence to accurate solutions. Parameters needed for incomplete LU preconditioning are automatically computed and may be updated at run time depending on the convergence pattern. The software is easy to use by non-experts and its top level routines are written in FORTRAN 77. Its potentialities are demonstrated on a few applications taken from computational physics.

Program summary

Program title: JADAMILUCatalogue identifier: ADZT_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZT_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 101 359No. of bytes in distributed program, including test data, etc.: 7 493 144Distribution format: tar.gzProgramming language: Fortran 77Computer: Intel or AMD with g77 and pgf; Intel EM64T or Itanium with ifort; AMD Opteron with g77, pgf and ifort; Power (IBM) with xlf90.Operating system: Linux, AIXRAM: problem dependentWord size: real:8; integer: 4 or 8, according to user's choiceClassification: 4.8Nature of problem: Any physical problem requiring the computation of a few eigenvalues of a symmetric matrix.Solution method: Jacobi-Davidson combined with multilevel ILU preconditioning.Additional comments: We supply binaries rather than source code because JADAMILU uses the following external packages:

•

MC64. This software is copyrighted software and not freely available. COPYRIGHT (c) 1999 Council for the Central Laboratory of the Research Councils.

•

AMD. Copyright (c) 2004-2006 by Timothy A. Davis, Patrick R. Amestoy, and Iain S. Duff. All Rights Reserved. Source code is distributed by the authors under the GNU LGPL licence.

•

BLAS. The reference BLAS is a freely-available software package. It is available from netlib via anonymous ftp and the World Wide Web.

•

LAPACK. The complete LAPACK package or individual routines from LAPACK are freely available on netlib and can be obtained via the World Wide Web or anonymous ftp.

•

For maximal benefit to the community, we added the sources we are proprietary of to the tar.gz file submitted for inclusion in the CPC library. However, as explained in the README file, users willing to compile the code instead of using binaries should first obtain the sources for the external packages mentioned above (email and/or web addresses are provided).

Running time: Problem dependent; the test examples provided with the code only take a few seconds to run; timing results for large scale problems are given in Section 5.     

KANTBP 2.0: New version of a program for computing energy levels, reaction matrix and radial wave functions in the coupled-channel hyperspherical adiabatic approach

A FORTRAN 77 program for calculating energy values, reaction matrix and corresponding radial wave functions in a coupled-channel approximation of the hyperspherical adiabatic approach is presented. In this approach, a multi-dimensional Schrödinger equation is reduced to a system of the coupled second-order ordinary differential equations on a finite interval with homogeneous boundary conditions: (i) the Dirichlet, Neumann and third type at the left and right boundary points for continuous spectrum problem, (ii) the Dirichlet and Neumann type conditions at left boundary point and Dirichlet, Neumann and third type at the right boundary point for the discrete spectrum problem. The resulting system of radial equations containing the potential matrix elements and first-derivative coupling terms is solved using high-order accuracy approximations of the finite element method. As a test desk, the program is applied to the calculation of the reaction matrix and radial wave functions for 3D-model of a hydrogen-like atom in a homogeneous magnetic field. This version extends the previous version 1.0 of the KANTBP program [O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Commun. 177 (2007) 649-675].

Program summary

Program title: KANTBPCatalogue identifier: ADZH_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZH_v2_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 20 403No. of bytes in distributed program, including test data, etc.: 147 563Distribution format: tar.gzProgramming language: FORTRAN 77Computer: Intel Xeon EM64T, Alpha 21264A, AMD Athlon MP, Pentium IV Xeon, Opteron 248, Intel Pentium IVOperating system: OC Linux, Unix AIX 5.3, SunOS 5.8, Solaris, Windows XPRAM: This depends on

1.

the number of differential equations;

2.

the number and order of finite elements;

3.

the number of hyperradial points; and

4.

the number of eigensolutions required.

The test run requires 2 MBClassification: 2.1, 2.4External routines: GAULEG and GAUSSJ [2]Nature of problem: In the hyperspherical adiabatic approach [3-5], a multidimensional Schrödinger equation for a two-electron system [6] or a hydrogen atom in magnetic field [7-9] is reduced by separating radial coordinate ρ from the angular variables to a system of the second-order ordinary differential equations containing the potential matrix elements and first-derivative coupling terms. The purpose of this paper is to present the finite element method procedure based on the use of high-order accuracy approximations for calculating approximate eigensolutions of the continuum spectrum for such systems of coupled differential equations on finite intervals of the radial variable ρ∈[ρmin,ρmax]. This approach can be used in the calculations of effects of electron screening on low-energy fusion cross sections [10-12].Solution method: The boundary problems for the coupled second-order differential equations are solved by the finite element method using high-order accuracy approximations [13]. The generalized algebraic eigenvalue problem AF=EBF with respect to pair unknowns (E,F) arising after the replacement of the differential problem by the finite-element approximation is solved by the subspace iteration method using the SSPACE program [14]. The generalized algebraic eigenvalue problem (A−EB)F=λDF with respect to pair unknowns (λ,F) arising after the corresponding replacement of the scattering boundary problem in open channels at fixed energy value, E, is solved by the LDLT factorization of symmetric matrix and back-substitution methods using the DECOMP and REDBAK programs, respectively [14]. As a test desk, the program is applied to the calculation of the reaction matrix and corresponding radial wave functions for 3D-model of a hydrogen-like atom in a homogeneous magnetic field described in [9] on finite intervals of the radial variable ρ∈[ρmin,ρmax]. For this benchmark model the required analytical expressions for asymptotics of the potential matrix elements and first-derivative coupling terms, and also asymptotics of radial solutions of the boundary problems for coupled differential equations have been produced with help of a MAPLE computer algebra system.Restrictions: The computer memory requirements depend on:

1.

the number of differential equations;

2.

the number and order of finite elements;

3.

the total number of hyperradial points; and

4.

the number of eigensolutions required.

Restrictions due to dimension sizes may be easily alleviated by altering PARAMETER statements (see Section 3 and [1] for details). The user must also supply subroutine POTCAL for evaluating potential matrix elements. The user should also supply subroutines ASYMEV (when solving the eigenvalue problem) or ASYMS0 and ASYMSC (when solving the scattering problem) which evaluate asymptotics of the radial wave functions at left and right boundary points in case of a boundary condition of the third type for the above problems.Running time: The running time depends critically upon:

1.

the number of differential equations;

2.

the number and order of finite elements;

3.

the total number of hyperradial points on interval [ρmin,ρmax]; and

4.

the number of eigensolutions required.

The test run which accompanies this paper took 2 s without calculation of matrix potentials on the Intel Pentium IV 2.4 GHz.References:[1] O. Chuluunbaatar, A.A. Gusev, A.G. Abrashkevich, A. Amaya-Tapia, M.S. Kaschiev, S.Y. Larsen, S.I. Vinitsky, Comput. Phys. Commun. 177 (2007) 649-675; http://cpc.cs.qub.ac.uk/summaries/ADZHv10.html.[2] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes: The Art of Scientific Computing, Cambridge University Press, Cambridge, 1986.[3] J. Macek, J. Phys. B 1 (1968) 831-843.[4] U. Fano, Rep. Progr. Phys. 46 (1983) 97-165.[5] C.D. Lin, Adv. Atom. Mol. Phys. 22 (1986) 77-142.[6] A.G. Abrashkevich, D.G. Abrashkevich, M. Shapiro, Comput. Phys. Commun. 90 (1995) 311-339.[7] M.G. Dimova, M.S. Kaschiev, S.I. Vinitsky, J. Phys. B 38 (2005) 2337-2352.[8] O. Chuluunbaatar, A.A. Gusev, V.L. Derbov, M.S. Kaschiev, L.A. Melnikov, V.V. Serov, S.I. Vinitsky, J. Phys. A 40 (2007) 11485-11524.[9] O. Chuluunbaatar, A.A. Gusev, V.P. Gerdt, V.A. Rostovtsev, S.I. Vinitsky, A.G. Abrashkevich, M.S. Kaschiev, V.V. Serov, Comput. Phys. Commun. 178 (2007) 301 330; http://cpc.cs.qub.ac.uk/summaries/AEAAv10.html.[10] H.J. Assenbaum, K. Langanke, C. Rolfs, Z. Phys. A 327 (1987) 461-468.[11] V. Melezhik, Nucl. Phys. A 550 (1992) 223-234.[12] L. Bracci, G. Fiorentini, V.S. Melezhik, G. Mezzorani, P. Pasini, Phys. Lett. A 153 (1991) 456-460.[13] A.G. Abrashkevich, D.G. Abrashkevich, M.S. Kaschiev, I.V. Puzynin, Comput. Phys. Commun. 85 (1995) 40-64.[14] K.J. Bathe, Finite Element Procedures in Engineering Analysis, Englewood Cliffs, Prentice-Hall, New York, 1982.