IAEA1382 SRNA-2K5.
SRNA-2K5, Proton Transport Using 3-D by Monte Carlo Techniques
NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHODS, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES, RELATED OR AUXILIARY PROGRAMS, STATUS, REFERENCES, HARDWARE REQUIREMENTS, LANGUAGE, SOFTWARE REQUIREMENTS, OTHER RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES
[ top ]
[ top ]
To submit a request, click below on the link of the version you wish to order.
Only liaison officers are authorised to submit online requests. Rules for requesters are
available here.
| Program name | Package id | Status | Status date |
|---|---|---|---|
| SRNA-2K5 | IAEA1382/03 | Tested | 18-MAY-2011 |
Machines used:
| Package ID | Orig. computer | Test computer |
|---|---|---|
| IAEA1382/03 | Linux-based PC | Linux-based PC,PC Windows |
[ top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION
SRNA-2K5 performs Monte Carlo transport simulation of proton in 3D source and 3D geometry of arbitrary materials. The proton transport based on condensed history model, and on model of compound nuclei decays that creates in nonelastic nuclear interaction by proton absorption.
SRNA-2K5 performs Monte Carlo transport simulation of proton in 3D source and 3D geometry of arbitrary materials. The proton transport based on condensed history model, and on model of compound nuclei decays that creates in nonelastic nuclear interaction by proton absorption.
[ top ]
4. METHODS
The SRNA-2K5 package is developed for time independent simulation of proton transport by Monte Carlo techniques for numerical experiments in complex geometry, using PENGEOM from PENELOPE with different material compositions, and arbitrary spectrum of proton generated from the 3D source. This package developed for 3D proton dose distribution in proton therapy and dosimetry, and it was based on the theory of multiple scattering. The compound nuclei decay was simulated by our and Russian MSDM models using ICRU 49 and ICRU 63 data. If protons trajectory is divided on great number of steps, protons passage can be simulated according to Berger's Condensed Random Walk model. Conditions of angular distribution and fluctuation of energy loss determinate step length.
Physical picture of these processes is described by stopping power, Moliere's angular distribution , Vavilov's distribution with Sulek's correction per all electron orbits, and Chadwick's cross sections for nonelastic nuclear interactions, obtained by his GNASH code. According to physical picture of protons passage and with probabilities of protons transition from previous to next stage, which is prepared by SRNADAT program, simulation of protons transport in all SRNA codes runs according to usual Monte Carlo scheme:
(i) proton from the spectrum prepared for random choice of energy, position and space angle is emitted from the source;
(ii) proton is loosing average energy on the step;
(iii) on that step, proton experience a great number of collisions, and it changes direction of movement randomly chosen from angular distribution;
(iv) random fluctuation is added to average energy loss; (v) protons step is corrected with data about protons position before and after scattering;
(vi) there is final probability on step for nonelastic nuclear interaction to happen, and for proton to be absorbed. Compound nucleus decays with emission of protons, neutrons, deuterons, tritons, alpha particles or photons. Particular decay particle is sampled from Poisson's distribution with appropriate average values of multiplication factor of each particle. Energy and angle of particle emission and factors of multiplication are determined from the cross section that obtained by the integration of differential cross section for nonelastic nuclear interaction. Energy and angle of secondary neutron are sampled from emission spectrum. Neutron and photon transport are not included in the current model. They are registered in data file and can be used by other code to simulate their transport. Emitted deuteron, triton and alpha particles are absorbed at the place their creation.
The SRNA-2K5 package is developed for time independent simulation of proton transport by Monte Carlo techniques for numerical experiments in complex geometry, using PENGEOM from PENELOPE with different material compositions, and arbitrary spectrum of proton generated from the 3D source. This package developed for 3D proton dose distribution in proton therapy and dosimetry, and it was based on the theory of multiple scattering. The compound nuclei decay was simulated by our and Russian MSDM models using ICRU 49 and ICRU 63 data. If protons trajectory is divided on great number of steps, protons passage can be simulated according to Berger's Condensed Random Walk model. Conditions of angular distribution and fluctuation of energy loss determinate step length.
Physical picture of these processes is described by stopping power, Moliere's angular distribution , Vavilov's distribution with Sulek's correction per all electron orbits, and Chadwick's cross sections for nonelastic nuclear interactions, obtained by his GNASH code. According to physical picture of protons passage and with probabilities of protons transition from previous to next stage, which is prepared by SRNADAT program, simulation of protons transport in all SRNA codes runs according to usual Monte Carlo scheme:
(i) proton from the spectrum prepared for random choice of energy, position and space angle is emitted from the source;
(ii) proton is loosing average energy on the step;
(iii) on that step, proton experience a great number of collisions, and it changes direction of movement randomly chosen from angular distribution;
(iv) random fluctuation is added to average energy loss; (v) protons step is corrected with data about protons position before and after scattering;
(vi) there is final probability on step for nonelastic nuclear interaction to happen, and for proton to be absorbed. Compound nucleus decays with emission of protons, neutrons, deuterons, tritons, alpha particles or photons. Particular decay particle is sampled from Poisson's distribution with appropriate average values of multiplication factor of each particle. Energy and angle of particle emission and factors of multiplication are determined from the cross section that obtained by the integration of differential cross section for nonelastic nuclear interaction. Energy and angle of secondary neutron are sampled from emission spectrum. Neutron and photon transport are not included in the current model. They are registered in data file and can be used by other code to simulate their transport. Emitted deuteron, triton and alpha particles are absorbed at the place their creation.
[ top ]
[ top ]
6. TYPICAL RUNNING TIME
The adopted parameters (energy cutoffs, geometry zones, etc.) have an influence on the computing time. As an example, a pencil beam depth-dose distribution of 250 MeV protons incident on a water phantom, obtained by simulating 100.000 histories, can be obtained with a running time of 4.8 minutes on a Pentium IV 1.67 GHz 512 MB RAM.
The adopted parameters (energy cutoffs, geometry zones, etc.) have an influence on the computing time. As an example, a pencil beam depth-dose distribution of 250 MeV protons incident on a water phantom, obtained by simulating 100.000 histories, can be obtained with a running time of 4.8 minutes on a Pentium IV 1.67 GHz 512 MB RAM.
[ top ]
[ top ]
[ top ]
IAEA1382/03, included references:
- Radovan D. Ilic:SRNA-2K5 - Protons Transport Simulation by Monte Carlo Techniques
Version 2K5, March, 2005 (User's Guide)
- Radovan D. Ilic et al.:
SRNA - Monte Carlo Codes for Proton Transport Simulation in Combined and
Voxelized Geometries
NT&RP 2002 XVII 1-2 p. 27-36 http://ntrp.vin.bg.ac.yu/
- Radovan D. Ilic et al.:
The Monte Carlo SRNA-VOX Code for 3D Proton Dose Distribution in Voxelized
Geometry Using CT Data
Phys. Med. Biol. 50 (2005) 1011-1017 (Feb. 2005)
[ top ]
IAEA1382/03
This package was tested at the NEA Databank on:- COMPUTER : apple iMac i7 2800MHz
IBM T61p Intel Core2Duo T7700 @2400MHz
- OPERATING SYSTEM : MS Windows XP Professional
ubuntu Linux 10.4 Kernel 2.6
- COMPILER : Compac Visual Fortran 6.1
GNU gfortran Version 4.4.3 (based on gcc 4.4.3)
[ top ]
[ top ]
[ top ]
[ top ]
IAEA1382/03
SRNADOC\ Documentation directory including manual and references in PDFSRLINPEN\ SRNA-2K5 files on the working folder
GEOMPEN\ PENGEOM examples and viewer
SRNATEST\ Input and output SRNA-2K5 test case files
Keywords: Monte Carlo method, geometry, protons, three-dimensional, transport.