RSICC CODE PACKAGE CCC-667

1. NAME AND TITLE

SHIELD: Monte Carlo Transport Code for Simulating Interaction of High Energy Hadrons with Complex Macroscopic Targets.

2. CONTRIBUTOR

Institute for Nuclear Research of the Russian Academy of Science (RAS), Moscow, Russia.

3. CODING LANGUAGE AND COMPUTER

Fortran 77; Sun, IBM PC, Vax (C00667SUN0500).

4. NATURE OF PROBLEM SOLVED

The SHIELD code considers interaction of high energy particles with condensed matter, including hadron-nucleus interactions inside the target, generation and transportation of secondary particles, deposition of energy and production of radionuclides in the target. The modern version of the SHIELD code allows simulation of the transfer of nucleons (including low energy neutrons), pions, kaons, antinucleons, and muons in energy range up to 1 TeV. Recently, the transfer of ions (arbitrary A,Z-nuclei) was added. The ionization loss and straggling (optionally) are taken into account as well as the main modes of the mesons decay. The transfer of neutrons (En<14.5 MeV) is simulated on the basis of the 28-group neutron data system BNAB. A special interface allows one to use the known EGS4 code for simulation of EM-showers initiated by products of mesons decay. Programs for simulation of nuclear reactions, included in the hadron nucleus generator of the SHIELD code (Many Stage Dynamical Model-MSDM), were elaborated by N.S.Amelin, K.K.Gudima and V.D.Toneev, Joint Institute for Nuclear Research, Dubna, and by A.S.Botvina, Institute for Nuclear Research RAS, Moscow.

5. METHOD OF SOLUTION

The Monte Carlo method is used for simulation of the interaction process. Direct simulating is employed as a rule; weighted simulating mode is foreseen in the code also. Exclusive approach to modeling of inelastic hadron-nucleus interactions is realized.

6. RESTRICTIONS OR LIMITATIONS

The SHIELD code allows one to simulate transfer of nucleons (including low energy neutrons above thermal energy), pions, kaons, antinucleons, and muons in energy range up to 1 TeV at arbitrary configuration and chemical composition of the target. Within this scope any limitations are practically absent. The only restriction concerns transport of neutrons below 14.5 MeV: BNAB neutron data library, used in SHIELD, includes near 40 chemical elements, but not the whole Mendeleev's Periodical Table.

7. TYPICAL RUNNING TIME

Running time depends essentially on the target dimensions, incident energy, and task set-up (efficiency of registration). Roughly typical running time is 1-10 hours for PC Pentium-200.

8. COMPUTER HARDWARE REQUIREMENTS

The SHIELD code was successfully implemented by the author on: Vax, DEC ALPHA under VAX/VMS; CONVEX under Unix; Sun under Unix; and IBM PC under MS DOS with NDP Fortran compiler and under WINDOWS 95 with MS Fortran Power Station compiler.

9. COMPUTER SOFTWARE REQUIREMENT

A standard Fortran compiler with standard functions (SIN etc.) and built-in uniform random number generator are required. The source code was compiled and executed at RSICC on a SUN SPARC 60 running SUN OS 5.6 (SOLARIS 2.6) with Fortran f77 version 4.2. No executables are included.

10. REFERENCES

a. included in documentation:

N. M. Sobolevsky, "The SHIELD Code (Version 1996.hadr.0) Short User's Manual," Informal report (October 20, 1998). A comprehensive list of references is included in the User's Manual.

b. background information:

A. V. Dementyev and N. M. Sobolevsky, "SHIELD - A Monte Carlo Hadron Transport Code. Proc. of a Specialists' Meeting "Intermediate Energy Nuclear Data: Models and Codes," Issy-les-Moulineaux, France, p 237 (30 May-1 June 1994).

11. CONTENTS OF CODE PACKAGE

Included are the referenced document in 10.a and a GNU compressed tar file transmitted on one DS/HD diskette in DOS format. The tar file contains the User's Manual, Fortran source, data, and test case input and output.

12. DATE OF ABSTRACT

July 2001.

KEYWORDS: CHARGED PARTICLES; HADRON; HIGH ENERGY; MONTE CARLO; SPALLATION