Overview of nuclear fragmentation models and needs.

It has been known for some time that adequate assessment of spacecraft shield requirements and concomitant estimates of astronauts radiation exposures from galactic cosmic radiation requires accurate, quantitative methods for characterizing these radiation fields as they pass through thick absorbers. The main nuclear interaction processes involved are (1) nuclear elastic and inelastic collisions, and (2) nuclear breakup (fragmentation) and electromagnetic dissociation (EMD). Nuclear fragmentation and EMD are important because they alter the elemental and isotopic composition of the transported radiation fields. At present, there is no suitably accurate theory for predicting nuclear fragmentation cross sections for all collision pairs and energies of interest in space radiation protection. Typical cross-section differences between theory and experiment range from about 25 percent to a factor of two. The resulting errors in transported flux, for high linear energy transfer (LET) particles, are comparble to these cross-section errors. In this overview, theoretical models of heavy ion fragmentation currently used to generate input data bases for cosmic-ray transport and shielding codes are reviewed. Their shortcomings are discussed. Further actions needed to improve their accuracy and generality are presented.     

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

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

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

The atmospheric cosmic- and solar energetic particle radiation environment at aircraft altitudes

Galactic cosmic rays interact with the solar wind, the earth's magnetic field and its atmosphere to produce hadron, lepton and photon fields at aircraft altitudes. In addition to cosmic rays, energetic particles generated by solar activity bombard the earth from time to time. These particles, while less energetic than cosmic rays, also produce radiation fields at aircraft altitudes which have qualitatively the same properties as atmospheric cosmic rays. We have used a code based on transport theory to calculate atmospheric cosmic-ray quantities and compared them with experimental data. Agreement with these data is seen to be good. We have then used this code to calculate equivalent doses to aircraft crews. We have also used the code to calculate radiation doses from several large solar energetic particle events which took place in 1989, including the very large event that occurred on September 29^th and 30^th of that year. The spectra incident on the atmosphere were determined assuming diffusive shock theory.     

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.     

太阳宇宙线源物质的高色球层模式和重离子丰度过量机制的探讨

本文利用太阳能量粒子事件中重离子平均丰度过量的资料,计算得到太阳能量粒子源物质的温度,提出了描述太阳宇宙线能量粒子源物质的新模式——高色球层模式;太阳耀斑观测确定,太阳宇宙线耀斑的加速区一般最可能出现在低日冕甚至高达几万公里的高度,从而,太阳宇宙线的源和加速区通常不位于同一区域;进而提出了描述太阳能量粒子事件中重离子丰度过量的可能机制——其源物质是通过太阳黑子的冻结型无力场从高色球层输送到活动区,形成耀斑前加速区内重离子丰度大和耀斑后宇宙线中重元素丰度的过量.      

The physics of cosmic-ray modulation

The theory of the modulation of galactic cosmic rays by the solar wind is reviewed. The basic transport equation is presented, interpreted and then applied to cosmic-ray transport in a model heliosphere immersed in a constant uniform bath of galactic cosmic rays. The results of numerical modelling are presented and the dominant physical effects analyzed. A variety of observational tests of the model which were reported over the last several years are summarized and shown, generally, to support a model in which particle drifts play an important role. Recent measurements which show that the latitudinal gradient of cosmic rays changed sign in the recent sunspot minimum (relative the last sunspot minimum) are shown to provide additional, strong, support for the model. A new picture of the interplanetary magnetic field is presented, which gives promise of improving considerably the agreement between the theory and observations in the few remaining problem areas.     

Transport of light ions in matter

A recent set of light ion experiments are analyzed using the Green's function method of solving the Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 code reasonably represents the fragmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model including shell structure and short range correlations appearing as tightly bound clusters in the light ion nucleus. The most recent NUCFRG2 code is augmented with a quasielastic alpha knockout model and semiempirical adjustments (up to 30 percent in charge removal) in the fragmentation process allowing reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross sections must await the full development of a coupled channel reaction model in which shell structure and clustering can be accurately evaluated.     

On the transport of cosmic rays in the distant heliosheath

This paper discusses the transport of energetic charged particles through a sectored magnetic field in distant regions of the inner heliosheath. As the plasma flow slows down on approach to the stagnation point on the heliopause, the distance between the folds of the current sheet decreases to the point where it becomes comparable to the cyclotron radius of a cosmic ray particle. Under these conditions a particle can effectively drift across the stack of magnetic sectors with a speed comparable with the particle’s velocity. For a random distribution of current sheet separation distances, a diffusive transport across the stack of sectors occurs instead. The proposed mechanism could have contributed to unusually high intensities of galactic cosmic rays measured by Voyager 1 in the heliosheath during 2009–2010.     

The FLUKA code for space applications: recent developments.

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions.     

Cosmic ray reacceleration on the galactic wind termination shock

Reacceleration of cosmic rays produced by galactic sources on the galactic wind termination shock is considered. The problem of the cosmic ray spectrum continuity is investigated. Numeric results are presented and discussed. We found that a smooth spectral transition from the galactic cosmic rays to the cosmic rays reaccelerated at the galactic wind termination shock is difficult to produce, if the maximum energy of accelerated particles is the same throughout the surface of the termination shock. The possible solution of this problem is the non-spherical termination shock with different maximum energies at different places of the shock.     

Secondary protons and neutrons generated by galactic and solar cosmic ray particles behind 1–100 g/cm aluminium shielding

Proton and neutron energy spectra behind AI shielding affected by galactic and solar cosmic ray particles are Monte-Carlo calculated using the multipurpose hadron transport code (SHIELD).     

Computation of radiation environment during ground level enhancements 65, 69 and 70 at equatorial region and flight altitudes

The radiation environment in the troposphere of the Earth is governed by cosmic rays of galactic and solar origin. During major solar energetic particles events the radiation environment changes dramatically. As a results the risk of biological effects due to exposure to ionizing radiation of aircrew increases. Here we present a numerical model for computation of absorbed dose in air due to cosmic rays of galactic and solar origin. It is applied for computation of radiation environment at flight altitude in the equatorial region during several major ground level enhancements, namely GLE65 on 28 October 2003, GLE69 on 20 January 2005 and GLE70 on 13 December 2006. The model is based on a full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The cascade simulation is carried out with CORSIKA 6.990 code with corresponding hadron generators FLUKA 2011 and QGSJET II. The contribution of different cascade components, namely electromagnetic, hadron and muon is explicitly obtained. The spectra of arriving solar energetic particles are calculated from ground level measurements with neutron monitors and satellite data from GOES. The obtained results are discussed.     

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.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

Cosmic ray hit frequencies in critical sites in the central nervous system.

One outstanding question to be addressed in assessing the risk of exposure to space travelers from galactic cosmic rays (GCR) outside the geomagnetosphere is to ascertain the effects of single heavy-ion hits on cells in critical regions of the central nervous system (CNS). As a first step toward this end, it is important to determine how many "hits" might be received by a neural cell in several critical CNS areas during an extended mission outside the confines of the earth's magnetic field. Critical sites in the CNS: the macula, and an interior brain point (typical of the genu, thalamus, hippocampus and nucleus basalis of Meynert) were chosen for the calculation of hit frequencies from galactic cosmic rays for a mission to Mars during solar minimum (i.e., at maximum cosmic-ray intensity). The shielding at a given position inside the body was obtained using the Computerized Anatomical Man (CAM) model, and a radiation transport code which includes nuclear fragmentation was used to calculate yearly fluences at the point of interest. Since the final Mars spacecraft shielding configuration has not yet been determined, we considered the minimum amount of aluminum required for pressure vessel-wall requirements in the living quarters of a spacecraft, and a typical duty area as a pressure vessel plus necessary equipment. The conclusions are: (1) variation of the position of the "target site" within the head plays only a small role in varying hit frequencies; (2) the average number of hits depends linearly on the cross section of the critical portion of the cell assumed in the calculation; (3) for a three-year mission to Mars at solar minimum (i.e., assuming the 1977 spectrum of galactic cosmic rays), 2% or 13% of the "critical sites" of cells in the CNS would be directly hit at least once by iron ions, depending on whether 60 micrometers2 or 471 micrometers2 is assumed as the critical cross sectional area; and (4) roughly 6 million out of some 43 million hippocampal cells and 55 thousand out of 1.8 million thalamus cell nuclei would be directly hit by iron ions at least once on such a mission for space travelers inside a simple pressure vessel. Also, roughly 20 million out of 43 million hippocampal cells and 230 thousand out of 1.8 million thalamus cell nuclei would be directly hit by one or more particles with z > or = 15 on such a mission.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.