The physics of the FLUKA code: Recent developments

FLUKA is a Monte-Carlo code able to simulate interaction and transport of hadrons, heavy ions and electromagnetic particles from few keV (or thermal neutron) to cosmic ray energies in whichever material. The highest priority in the design and development of the code has always been the implementation and improvement of sound and modern physical models. A summary of the FLUKA physical models is given, while recent developments are described in detail: among the others, extensions of the intermediate energy hadronic interaction generator, refinements in photon cross sections and interaction models, analytical on-line evolution of radio-activation and remnant dose. In particular, new developments in the nucleus–nucleus interaction models are discussed. Comparisons with experimental data and examples of applications of relevance for space radiation are also provided.     

Normalized ionization yield function for various nuclei obtained with full Monte Carlo simulations

Several recent results important for production of ion pairs in the Earth atmosphere by various primary cosmic ray nuclei are presented. The direct ionization by various primary cosmic ray nuclei is explicitly obtained. The longitudinal profile of atmospheric cascades is sensitive to the energy and mass (charge) of the primary particle. In this study different cosmic ray nuclei are considered as primaries, namely Helium, Oxygen and Iron nuclei. The cosmic ray induced ionization is obtained on the basis of CORSIKA 6.52 code simulations using FLUKA 2006 and QGSJET II hadronic interaction models. The energy of the primary particles is normalized to GeV per nucleon. In addition, the ionization yield function Y is normalized as ion pair production per nucleon. The obtained ionization yield functions Y for various primaries are compared. The presented results and their application are discussed.     

Ionization effect of solar protons in the Earth atmosphere – Case study of the 20 January 2005 SEP event

The cosmic ray ground level enhancement on January 20, 2005 is among the largest recorded events in the history of cosmic ray measurements. The solar protons of MeV energies cause an excess of ionization in the atmosphere, specifically over polar caps following major solar disturbances. The ionization effect in the Earth atmosphere is obtained for various latitudes on the basis of solar proton energy spectra, reconstructed from GOES 11 measurements and subsequent full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The estimation of ionization rates is based on a numerical model for cosmic ray induced ionization. The evolution of atmospheric cascade is performed with the CORSIKA 6.52 code using FLUKA 2006b and QGSJET II hadron interaction models. The atmospheric ion rate ionization is explicitly obtained for 40°N, 60°N and 80°N latitudes. The time evolution of obtained ion rates is presented. It is demonstrated that ionization effect is negative for 40°N and small for 60°N, because of accompanying Forbush decrease. The ionization effect is significant only in sub-polar and polar atmosphere during the major ground level enhancement of 20 January 2005.     

空间辐射场蒙特卡罗模拟计算方法研究

对利用蒙特卡罗方法对由银河宇宙射线引起的空间辐射场各成分进行计算的方法进行了调研,对计算模型的建立以及计算过程中通常使用的方差减小技术进行分析,给出了美国的Roesler等人利用FLUKA程序以及加拿大Anid等人利用MCNPX程序计算得到的由银河宇宙射线引起的空间辐射场各量值及其与实验结果的比较,验证了计算方法与计算模型的可靠性。对任意航线空间辐射场剂量分布预评估方法进行分析,给出了由银河宇宙射线引起的空间辐射场的基本特征。     

Resolving electrons from protons in ATIC

The Advanced Thin Ionization Calorimeter (ATIC) experiment is designed for high energy cosmic ray ion detection. The possibility to identify high energy primary cosmic ray electrons in the presence of the ‘background’ of cosmic ray protons has been studied by simulating nuclear-electromagnetic cascade showers using the FLUKA Monte Carlo simulation code. The ATIC design, consisting of a graphite target and an energy detection device, a totally active calorimeter built up of 2.5 cm × 2.5 cm × 25.0 cm BGO scintillator bars, gives sufficient information to distinguish electrons from protons. While identifying about 80% of electrons as such, only about 2 in 10,000 protons (@ 150 GeV) will mimic electrons. In September of 1999 ATIC was exposed to high-energy electron and proton beams at the CERN H2 beam line, and this data confirmed the electron detection capabilities of ATIC. From 2000-12-28 to 2001-01-13 ATIC was flown as a long duration balloon test flight from McMurdo, Antarctica, recording over 360 h of data and allowing electron separation to be confirmed in the flight data. In addition, ATIC electron detection capabilities can be checked by atmospheric gamma-ray observations.     

Extension of INCL4 between 2 and 15 GeV

The intranuclear cascade model INCL4 has been shown to be very successful for describing, without adjustable parameters, a whole set of data for p-induced reactions in the 40 MeV–2 GeV energy range. In view of its possible application to cosmic ray interactions, the INCL4 code has been extended to the 2–15 GeV energy range, so covering a large part of the spectrum of the incident energy of the cosmic rays.     

Depth dependency of neutron density produced by cosmic rays in the lunar subsurface

Depth dependency of neutrons produced by cosmic rays (CRs) in the lunar subsurface was estimated using the three-dimensional Monte Carlo particle and heavy ion transport simulation code, PHITS, incorporating the latest high energy nuclear data, JENDL/HE-2007. The PHITS simulations of equilibrium neutron density profiles in the lunar subsurface were compared with the measurement by Apollo 17 Lunar Neutron Probe Experiment (LNPE). Our calculations reproduced the LNPE data except for the 350–400 mg/cm² region under the improved condition using the CR spectra model based on the latest observations, well-tested nuclear interaction models with systematic cross section data, and JENDL/HE-2007.     

Comparison of measured cosmic ray LET spectra with models and predictions.

Radiation effects of cosmic ray nuclei are generally described as a function of the particle LET. For a large number of space missions LET spectra have been measured and models have been developed to calculate these spectra that include the effects of geomagnetic shielding and shielding provided by material. In this paper we compare measured and calculated LET spectra. For low earth orbits events with high local energy deposition, i.e., short range secondaries, contribute significantly to the measured spectra. These events are produced by nuclear interactions, mainly induced by protons from the south atlantic anomaly. The technique to include these contributions in the models depends on the size of radiation sensitive volumes. For sizes comparable to or larger than the range of target secondaries it is essential to separate contributions by target interactions from those of cosmic rays. This separation is possible in experiments which use stacks of plastic nuclear track detectors. The yield of short range events generated by protons and measured in the detector can be calibrated from accelerator experimental data. We present first results for CR-39 detectors.     

Assessment and requirements of nuclear reaction databases for GCR transport in the atmosphere and structures

The transport properties of galactic cosmic rays (GCR) in the atmosphere, material structures, and human body (self-shielding) are of interest in risk assessment for supersonic and subsonic aircraft and for space travel in low-Earth orbit and on interplanetary missions. Nuclear reactions, such as knockout and fragmentation, present large modifications of particle type and energies of the galactic cosmic rays in penetrating materials. We make an assessment of the current nuclear reaction models and improvements in these model for developing required transport code data bases. A new fragmentation data base (QMSFRG) based on microscopic models is compared to the NUCFRG2 model and implications for shield assessment made using the HZETRN radiation transport code. For deep penetration problems, the build-up of light particles, such as nucleons, light clusters and mesons from nuclear reactions in conjunction with the absorption of the heavy ions, leads to the dominance of the charge Z = 0, 1, and 2 hadrons in the exposures at large penetration depths. Light particles are produced through nuclear or cluster knockout and in evaporation events with characteristically distinct spectra which play unique roles in the build-up of secondary radiation's in shielding. We describe models of light particle production in nucleon and heavy ion induced reactions and make an assessment of the importance of light particle multiplicity and spectral parameters in these exposures.     

On ultra-high energy cosmic rays: Origin in AGN jets and transport in expanding universe

The cosmic ray source spectrum produced by AGN (active galactic nucleus) jets is calculated. A distinctive feature of these calculations is the account for the jet distribution on kinetic energy. The expected cosmic ray spectrum at the Earth is determined with the use of a simple numerical code which takes into account interactions of ultra-high energy protons and nuclei with the background radiation in an expanding universe.     

Galactic cosmic ray composition and energy spectra.

Galactic cosmic ray nuclei represent a significant risk to long-duration spaceflight outside the magnetosphere. We review briefly existing measurements of the composition and energy spectra of heavy cosmic ray nuclei, pointing out which species and energy ranges are most critical to assessing cosmic ray risks for spaceflight. Key data sets are identified and a table of cosmic ray abundances is presented for elements from H to Ni (Z = 1 to 28). Because of the 22-year nature of the solar modulation cycle, data from the approaching 1998 solar minimum is especially important to reducing uncertainties in the cosmic ray radiation hazard. It is recommended that efforts to model this hazard take advantage of approaches that have been developed to model the astrophysical aspects of cosmic rays.     

The high-LET radiation component measured during the EUROMIR-94 mission.

Stacks of CR-39 plastic nuclear track detectors were mounted inside the MIR-station during the EUROMIR-94-mission. We present LET-spectra determined separately for long range cosmic ray heavy ions and for short range target fragments produced in nuclear interactions of cosmic rays and measured charge distributions for relativistic and stopping particles.     

Heavy particle fluxes in Salyut space station orbit

Results of a detailed analysis of heavy ion fluences measured on the Salyut and MIR orbital stations from 1978 to 1990 are presented. The analysis has made use of new models that describe the cosmic ray fluxes and their transition through the magnetosphere. The penetration of solar cosmic ray particles to the orbit, the increase of the (ScCr)/Fe flux ratio in the orbit, and the occurrences of 200–500 MeV/nucl heavy nuclei in <30° latitudes have been analyzed.     

Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

电化学过程中的局域核反应及动态Casimir效应

通过CR-39固体径迹探测器、钯阴极的辐射自照相等方法观察到了电解过程中的局域核反应.综合国内外相关研究数据,认为电解过程中,电极尖端处出现的动态Casimir效应,及涡漩运动与零点能相干效应而提取真空能,是电化学异常放热的来源,少量具有滞后和积累效应的、高度定向的极化核反应,可以用具有轴向高能宇宙射线的类星体涡旋模型加以解释.讨论了真空零点电磁涨落和挠场理论,及有关实验现象.利用推出的类星体涡旋模型理论,可以对自然界和实验室中出现的许多异常现象进行解释.      

Observation of enhanced sub-iron (Sc---Cr) to iron ration in low energy cosmic rays of 50–100 MeV/N in spacelab-3

Relative abundances of sub-iron (Sc-Cr) to iron nuclei in low energy (50–100 MeV/N) galactic cosmic rays have been determined from an analysis of about 100 events of heavy ions (Z = 10−28) recorded in a detector assembly flown in the Anuradha cosmic ray experiment in the Spacelab-3 on a six day mission in April–May 1985. The measured abundance ratio of (Sc-Cr)/Fe nuclei in 50–100 MeV/N energy range is 1.1 ± 0.3, and the present result of enhanced ratio of sub-iron to iron nuclei is in agreement with other experimental results in 200–800 MeV/N range. The over-abundance of iron secondaries at these low energies cannot be explained in the conventional models for propagation of cosmic rays. Available experimental data indicate a very different time history for the low energy iron-group, as compared to those of lighter nuclei in galactic cosmic rays.     

Calculated cosmic rays dynamics before the geomagnetic storms onset caused by the interplanetary shocks

Dynamics of the galactic cosmic ray intensity caused by their interactions with a shock front in the September 9, 1992 event has been determined. Corresponding variations of the cosmic ray intensity have been calculated for different stations of the world network of neutron monitors and muon telescopes of stations Nagoya and Sakashita. Comparison calculated results with observational data shows, in general, satisfactory consensus both on amplitude and in time. The developed method can be used for investigation dynamics of the solar wind disturbances precursors in the cosmic rays.