Program MSTAR

Description, and instructions for the user

MSTAR is a program that calculates electronic stopping powers for heavy ions ( 3Li to 18Ar). The basic description of its first version has been published [1]. For version 2.00 , the fitting functions and various details have been improved [2]. From version 2.11 on, the program can also calculate nickel targets. For version 3.00 (November 2002), data from recent publications (collected until September, 2002), and data for nickel targets have been included in the data base. All the fits for solid data have been redone. This has not changed mode 'c' much. But for mode 'd', it has become possible to determine a small residual Z 2-dependence of the normalized stopping power ratio Srel , well hidden in the statistical scatter of the data, and to produce a corresponding small correction to mode 'c'. Thus, mode 'd' is now available for all ions in solids, and it represents the data for most ions slightly better than mode 'c' (see statistical analysis).
Version 3.11 (April 2003) will also calculate B, Zr, Gd and Ta targets, but the data base has not been changed since version 3.00. A table based on this version has been published [3]. The present version 3.12 (May 2004) is mathematically identical to version 3.11, except for  an improved error function subroutine.

The program has been written in ANSI FORTRAN 77 in the form of a subroutine package with top-level routine MSTAR1 that can be linked to any other FORTRAN program requiring electronic stopping powers. Links to the particular routines are contained in the file Mssubs.f. A sample main program, MSTAR.FOR, has been written to demonstrate the use of subprogram MSTAR1; it can be used to calculate single values or tables of stopping powers.

MSTAR is mainly based upon the alpha stopping powers contained in ICRU Report 49 [4]. It will calculate the stopping of all the elements, mixtures and compounds contained in that report (and also of B, Ni, Zr, Gd, and Ta), for ions from 3Li to 18Ar. As in the report, targets are identified by an identification number (ID) which is identical to the atomic number for elemental targets. The ID numbers can be found in Tables 1.1 and 1.2 of the report [4], and in the (slightly augmented) Table 1 below. Table 1.2 of the report also gives the elemental compositions of mixtures and compounds.

MSTAR can also produce the alpha stopping powers that serve as input; except for B, Ni, Zr, Gd, and Ta, these are identical to those contained in ICRU Report 49.

Table 1. Identification (ID) Numbers of Targets

ID

Target Description

1

HYDROGEN

2

HELIUM

4

BERYLLIUM

5

BORON

6

AMORPHOUS CARBON (density 2.0 g/cm3)

7

NITROGEN

8

OXYGEN

10

NEON

13

ALUMINUM

14

SILICON

18

ARGON

22

TITANIUM

26

IRON

28

NICKEL

29

COPPER

32

GERMANIUM

36

KRYPTON

40

ZIRCONIUM

42

MOLYBDENUM

47

SILVER

50

TIN

54

XENON

64

GADOLINIUM

73

TANTALUM

74

TUNGSTEN

78

PLATINUM

79

GOLD

82

LEAD

92

URANIUM

99

A-150 TISSUE-EQUIVALENT PLASTIC

101

ACETYLENE

103

ADIPOSE TISSUE (ICRP)

104

AIR, DRY (NEAR SEA LEVEL)

106

ALUMINUM OXIDE

111

B-100 BONE-EQUIVALENT PLASTIC

119

BONE, COMPACT (ICRU)

120

BONE, CORTICAL (ICRP)

126

C-552 AIR-EQUIVALENT PLASTIC

130

CALCIUM FLUORIDE

134

CARBON DIOXIDE

138

CELLULOSE NITRATE

139

CERIC SULFATE DOSIMETER SOLUTION

141

CESIUM IODIDE

155

ETHYLENE

160

FERROUS SULFATE DOSIMETER SOLUTION

169

GLASS, BOROSILICATE (PYREX)

179

KAPTON POLYIMIDE FILM

185

LITHIUM FLUORIDE

189

LITHIUM TETRABORATE

191

M3 WAX

197

METHANE

200

MS20 TISSUE SUBSTITUTE

201

MUSCLE, SKELETAL (ICRP)

202

MUSCLE, STRIATED (ICRU)

203

MUSCLE-EQUIVALENT LIQUID, WITH SUCROSE

204

MUSCLE-EQUIVALENT LIQUID, WITHOUT SUCROSE

209

NYLON, TYPE 6 AND TYPE 616

213

PARAFFIN WAX

215

PHOTOGRAPHIC EMULSION

216

PLASTIC SCINTILLATOR (VINYLTOLUENE BASED)

219

POLYCARBONATE (MAKROLON, LEXAN)

221

POLYETHYLENE

222

POLYETHYLENE TEREPHTHALATE (MYLAR)

223

POLYMETHYL METHACRALATE (LUCITE, PERSPEX, PLEXIGLAS)

225

POLYPROPYLENE

226

POLYSTYRENE

227

POLYTETRAFLUOROETHYLENE (TEFLON)

232

POLYVINYL CHLORIDE

238

PROPANE

245

SILICON DIOXIDE

252

SODIUM IODIDE

255

STILBENE

263

TISSUE-EQUIVALENT GAS (METHANE BASED)

264

TISSUE-EQUIVALENT GAS (PROPANE BASED)

266

TOLUENE

276

WATER, LIQUID

277

WATER VAPOR

906

GRAPHITE (density 1.7 g/cm3)

For downloading, we have prepared an archive file MStar312.zip with three subdirectories. With the source codes we include a DOS executable of MSTAR that can be run, e.g., in the DOS window of Windows 9x/ME/NT (at the DOS prompt or started by double-clicking Mstar.exe in the Explorer). Then follow the instructions of the program.

The program MSTAR or any application calling MSTAR1 needs the database file Msdbs1.d in the same directory as the executable program. Following the first call to MSTAR1 this file is read into memory, which can take a few seconds on low performance systems. The user should keep this in mind, especially if the application is time critical. Additional calls to MSTAR1 do not read the database file again. It is not intended that the end-user should call any routines directly, except MSTAR1 and MSEMSG. For details on how to call these subroutines see the comment header of Mssubs.f and/or the headers of the particular routine.

In the top level subroutine MSTAR1, six modes of operation are available, which are distinguished by the use of different built-in fitting coefficients af, b f, cf :

Table 2. Modes of operation

Mode

Target [5]

Ion [5,6]

Remarks

Coefficients

'c'

any condensed

any

based on all the data for solids
(all ions taken together)

coefficients af,bf,cf Z2 -independent

'd'

any condensed [7]

any

based on the data for the particular ion only

coefficients af,bf,cf
Z2-dependent

'g'

any gaseous

any

based on data for all gases (all ions taken together), but H 2 and He treated separately.

coefficients af,bf,cf
Z2-independent

'h'

any gaseous,
except H, He

any except Mg,Al,Si,P

based on the data for all gases (except H2, He), but only one particular ion

coefficients af,bf,cf
Z2-independent [8]

'a'

any; select normal state of aggregation

any

always use c|g, even if d|h is available (MSTAR)

'b'

any; select normal state of aggregation

any

use d|h if available, otherwise use c|g (MSTAR)

In the sample main program MSTAR, only the modes 'a' and 'b' are available. The program prints out the actual mode ('c'|'d'|'g'|'h') used [9].

Concerning solids: mode 'd' should be superior to mode 'c', since it is specific to the ion, although mode 'c' is based on more data [2]. Indeed, mode 'd' represents the data better than mode 'c' (see table 3 below).
Concerning gases: mode 'h', though based on less data, appears to represent the measurements slightly better than mode 'g' (see table 4 below).

We suggest to use mode 'b' which should be good for most purposes.

REFERENCES AND FOOTNOTES.

[1] H.Paul and A.Schinner, "An empirical approach to the stopping power of solids and gases for ions from 3Li to 18Ar, Nucl. Instr. Meth. Phys. Res. B 179 (2001) 299

[2] H.Paul and A.Schinner, "An empirical approach to the stopping power of solids and gases for ions from 3Li to 18Ar, Part II, Nucl. Instr. Meth. Phys. Res. B 195 (2002) 166

[3] H. Paul and A. Schinner, “Empirical stopping power tables for ions from 3Li to 18Ar and from 0.001 to 1000 MeVnucleon in solids and gases”, Atomic Data Nucl. Data Tables 85 (2003) 377

[4] ICRU Report 49, "Stopping Powers and Ranges for Protons and Alpha Particles", Intern. Commission on Radiation Units and Measurements, Bethesda, Maryland (1993), and program ASTAR, version 2

[5] 'any target' or 'any ion' in this context means the supported range of targets (cf. table 1) or ions (3Li to 18Ar).

[6] MSTAR will also produce results for He ions; these are identical to those given in ref. [3].

[7] For targets where not enough data are available, or where the deviation from mode 'c' is negligible, MSTAR will choose to calculate mode 'c' even if mode 'd' was requested by the user.

[8]Not enough data are available for gases to determine the Z 2 -dependence in mode 'h'.

[9] In Ref. [1], modes 'c', 'd', 'g', and 'h' were called 1, 2, 3, and 4, respectively.

LICENSE AGREEMENT

This program and/or the corresponding subroutine package is hereby released as freeware. This means you can freely use, copy, modify, and/or distribute the source code. However, please do not remove the authors' names from the source code comments.

The authors cannot be made responsible for any damage caused by the use of their program code. YOU ARE USING THE ENTIRE PROGRAM SYSTEM OR ANY PART OF IT ENTIRELY AT YOUR OWN RISK!

On the other hand, we welcome bug reports and/or other feedback; please send emails to: nds.contact-point@iaea.org

By downloading this program you are accepting the above license agreement!

CLICK HERE TO DOWNLOAD MSTAR V3.12