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_v4_0.tar.gz(104 Kbytes)
Manuscript Title: ISICS2008: An expanded version of ISICS for calculating K-, L- and M-shell cross sections from PWBA and ECPSSR theory.
Authors: Sam J. Cipolla
Program title: ISICS2008
Catalogue identifier: ADDS_v4_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 180(2009)1716
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:
addition of relativistic treatment of both projectile and K-shell electrons.

Summary of revisions:
A new addition to ISICS is the option (R) to calculate ECPSSR cross sections that account for the relativistic treatment of both projectile and K-shell electron, as proposed recently by Lapicki [1], accordingly as

σKRECPSSR = CK • (1 + 0.07(θkςk)8(√(mRKυR1) ⁄ZzK)2 • σPWBA(√(mRKυR1) ⁄ ZzK , ςKθK)

where υR1 is the relativistic projectile velocity. The option can also be invoked in calculating ECPSShsR, where hsR stands for the Hartree-Slater description of the K-shell electron, which was already incorporated into ISICS2006 [2,3], and is now expressed in this option as,

σKRECPSShsR = CKhsR((2υR1) ⁄ ZzK ςKθK) , Zz⁄137) • (1 + 0.07(θkςk)8R1ZzK)2) • σPWBA(υR1) ⁄ ZzK , ςKθK)

using the function hsR that is already incorporated into ISICS2006. It should be noted that these expressions are corrected versions [4] from the ones published in ref[1].
In this new version, ISICS2008, the option line in the main menu that read "Use Relativisitic Proj. velocity" has been replaced by "R option for K shell....Uses Rel. Proj. vel.". As before, various combinations of options can be utilized and each is denoted in the output.

The consumed CPU time increases with the atomic shell (K,L,M), but execution is still very fast.

Additional comments:
A revised User Manual is included in the distribution file.

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. As before, 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.

[1] G. Lapicki, J. Phys. B: At. Mol. Op. Phys. 41 (2008) 115201.
[2] S. Cipolla, Comput. Phys. Comm. 176 (2007) 157.
[3] S. Cipolla, Nucl. Instrum. Meth Phys. Res. B 261 (2007) 142.
[4] G. Lapicki, private communication.