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] adds_v5_0.tar.gz(127 Kbytes)
Manuscript Title: ISICS2011, an updated version of ISICS: A program for calculation K-, L-, and M-shell cross sections from PWBA and ECPSSR theories using a personal computer.
Authors: Sam J. Cipolla
Program title: ISICS2011
Catalogue identifier: ADDS_v5_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 182(2011)2439
Programming language: C.
Computer: 80486 or higher-level PCs.
Operating system: WINDOWS XP and all earlier operating systems.
Keywords: Atomic K-, L-, M-Shell ionization, PWBA and ECPSSR cross sections, Ion-Atom collisions, C package, Gauss-Legendre quadrature integration techniques.
PACS: 34.50.Fa, 31.15.-p.
Classification: 16.7.

Does the new version supersede the previous version?: Yes

Nature 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:
General need for higher precision in output format for projectile energies; some built-in binding energies needed correcting; some anomalous results occur due to faulty read-in data or calculated parameters becoming un-physical; erroneous calculations could result for the L and M shells when restricted K-shell options are inadvertently chosen; to achieve general compatibility with ISICSoo, a companion C++ version that is portable to Linux and MacOS platforms, has been submitted for publication in the CPC Program Library approximately at the same time as this present new standalone version of ISICS [1].

Summary of revisions:
The format field for projectile energies in the output has been expanded from two to four decimal places in order to distinguish between closely spaced energy values.
There were a few entries in the executable binding energy file that needed correcting; K-shell of Eu, Mshells of Zn, M1 shell of Kr. The corrected values were also entered in the ENERGY.DAT file. In addition, an alternate data file of binding energies is included, called ENERGY_GW.DAT, which is more up-to-date [2]. Likewise, an alternate atomic parameters data file is now included, called FLOURE_JC.DAT, which is more up-to-date [3] fluorescence yields for the K and L shells and Coster-Kronig parameters for the L shell. Both data files can be read in using the -f usage option. To do this, the original energy file should be renamed and saved (e.g., ENERGY_BB.DAT) and the new file (ENERGY_GW.DAT ) should be duplicated as ENERGY.DAT to be read in using the -f option. Similarly for reading in an alternate FLOURE.DAT file. As with previous versions, the user can also simply input different values of any input quantity by invoking the "specify your own parameters" option from the main menu. You can also use this option to simply check the values of the built-in values of the parameters.
If it still happens that a zero binding energy for a particular sub-shell is read in, the program will not completely abort, but will calculate results for the other sub-shells while setting the affected sub-shell output to zero. In calculating the Coulomb deflection factor, if the quantity inside the radical sign of the parameter zs
       zs = √(1 - (4/Mζsθs) (ζss)2)
becomes zero or negative, to prevent the program from aborting, the PWBA cross sections are still calculated while the ECPSSR cross sections are set to zero. This situation can happen for very low energy collisions, such as were noticed for helium ions on copper at energies of E ≤ 11.2 keV.
It was observed during the engineering of ISICSoo [1] that erroneous calculations could result for the L and M shell cases when restricted K-shell R or hsr scaling options were inappropriately chosen. The program has now been fixed so that these inappropriate options are ignored for the L and M shells.
In the previous versions, the usage for inputting a batch data file was incorrectly stated in the Users Manual as -Bxxx; the correct designation is -Fxxx, or alternatively, -Ixxx, as indicated on the usage screen in running the program.
A revised Users Manual is also available.

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. The running time is not significantly changed from the previous version.

[1] M. Batic, M. G. Pia, S. Cipolla, submitted to Comput. Phys. Comm. (2011).
[2] www.jlab.org/~gwyn/ebindene.html
[3] J. Campbell, At. Data and Nucl. Data Tables 85 (2003) 291.