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.     

Forbush decreases and solar events seen in the 10–20 GeV energy range by the Karlsruhe Muon Telescope

Since 1993, a muon telescope located at Forschungszentrum Karlsruhe (Karlsruhe Muon Telescope) has been recording the flux of single muons mostly originating from primary cosmic-ray protons with dominant energies in the 10–20 GeV range. The data are used to investigate the influence of solar effects on the flux of cosmic rays measured at Earth. Non-periodic events like Forbush decreases and ground level enhancements are detected in the registered muon flux. A selection of recent events will be presented and compared to data from the Jungfraujoch neutron monitor. The data of the Karlsruhe Muon Telescope help to extend the knowledge about Forbush decreases and ground level enhancements to energies beyond the neutron monitor regime.     

Fluxes and nuclear abundances of cosmic rays inside the magnetosphere using a transmission function approach

At 1 AU and outside the Earth’s magnetosphere, the relative abundances to protons for He (He/p), C (C/p) and Fe (Fe/p) nuclei were calculated using the observation data of AMS-01 (for p and He) and HEAO-3 (for C and Fe) above 0.8 GeV/nucleon. In addition, the transmission function (TF) for the GCR propagation inside the magnetosphere was evaluated using the IGRF and T96 (introduced by Tsyganenko and Stern) models to obtain permitted and forbidden trajectories inside the magnetosphere. The TF allowed one to derive the primary He-nuclei fluxes in the same geomagnetic regions of AMS-01 observations. These fluxes were found in good agreement with the observation data. Furthermore inside the magnetosphere in addition to the flux of helium, it allowed one to obtain those of the primary p, C, and Fe nuclei and the relative abundances of He, C and Fe nuclei to protons from the same observation data of AMS-01 and HEAO-3 above ≈0.8 GeV/nucleon. Up to a geomagnetic latitude of ≈45.84°, the relative isotopic abundances were found to depend on the mass number I_isot and, on average, range from a factor ≈2.31 up to ≈3.35 larger than those outside the magnetosphere at 1 AU. Thus, the magnetospheric isotopic/nuclear relative abundances differ from those inside the solar cavity and those in the interstellar space. The usage of the TF approach can allow one to determine the nuclear abundances in the magnetosphere at any geomagnetic latitude and, thus, any orbit, provided that the CR spectra are determined at 1 AU.     

An update about recent developments of the PHITS code

PHITS (Particle and Heavy-Ion Transport code System) is a general-purpose three-dimensional Monte Carlo code, developed and maintained by RIST, JAEA and KEK in Japan together with Sihver et al. at Chalmers in Sweden. PHITS can deal with the transports of all varieties of hadrons and heavy ions with energies up to around 100 GeV/nucleon, and in this paper the current status of PHITS is presented. We introduce a relativistically covariant version of JQMD, called R-JQMD, that features an improved ground state initialization algorithm, and we will present the introduction of electron and photon transport in PHITS using EGS5, which have increased the energy region for the photon and energy transport from up to around 3 GeV to up to several hundred GeV depending on the atomic number of the target. We show how the accuracy in dose and fluence calculations can be improved by using tabulated cross sections. Benchmarking of shielding and irradiation effects of high energy protons in different materials relevant for shielding of accelerator facilities is also presented. In particular, we show that PHITS can be used for estimating the dose received by aircrews and personnel in space. In recent years, many countries have issued regulations or recommendations to set annual dose limitations for aircrews. Since estimation of cosmic-ray spectra in the atmosphere is an essential issue for the evaluation of aviation doses, we have calculated these spectra using PHITS. The accuracy of the atmospheric propagation simulation of cosmic-ray performed by PHITS has been well verified by experimental cosmic-ray spectra taken under various conditions. Based on a comprehensive analysis of the simulation results, an analytical model called “PARMA” has been proposed for instantaneously estimating the atmospheric cosmic-ray spectra below the altitude of 20 km. We have also performed preliminary simulations of long-term dose distribution measurements at the ISS performed with the joint ESA-FSA experiment MATROSHKA-R (MTR-R) led by the Russian Federation Institute of Biomedical Problems (IMBP) and the ESA supported experiment MATROSHKA (MTR), led by the German Aerospace Center (DLR). For the purpose of examining the applicability of PHITS to the shielding design in space, the absorbed doses in a tissue equivalent water phantom inside an imaginary space vessel has been estimated for different shielding materials of different thicknesses. The results confirm previous results which indicate that PHITS is a suitable tool when performing shielding design studies of spacecrafts.     

Space observations of comets during solar flares: A possible explanation for comet brightness outbursts

Problems connected with mechanisms for comet brightness outbursts as well as for gamma-ray bursts remain open. Meantime, calculations show that irradiation of a certain class of comet nuclei, having high specific electric resistance, by intense fluxes of energetic protons and positively charged ions with kinetic energies more than 1 MeV/nucleon, ejected from the Sun during strong solar flares, can produce a macroscopic high-voltage electric double layer with positive charge in the subsurface zone of the nucleus, during irradiation times of the order of 10–100 h at heliocentric distances around 1–10 AU. The maximum electric energy accumulated in such layer will be restricted by the electric discharge potential of the layer material. For comet nuclei with typical radii of the order of 1–10 km the accumulated energy of such natural electric capacitor is comparable to the energy of large comet outbursts that are estimated on the basis of ground based optical observations. The impulse gamma and X-ray radiation together with optical burst from the comet nucleus during solar flares, anticipated due to high-voltage electric discharge, may serve as an indicator of realization of the processes above considered. Multi-wavelength observations of comets and pseudo-asteroids of cometary origin, having brightness correlation with solar activity, using ground based optical telescopes as well as space gamma and X-ray observatories, during strong solar flares, are very interesting for the physics of comets as well as for high energy astrophysics.     

Statistical study of the detection of solar protons of highest energies at 20 January 2005

On January 20, 2005, 7:02–7:04 UT the Aragats Multichannel Muon Monitor (AMMM) registered enhancement of the high energy secondary muon flux (energy threshold ∼5 GeV). The enhancement, lasting 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES satellite and the ground level enhancement detected by the world-wide network of neutron monitors and by muon detectors. The most probable proton energy corresponding to the measured 5 GeV muon flux is within 23–30 GeV. Due to upmost importance of the detection of solar particles of highest energies in presented paper we perform detailed statistical analysis of the detected peak. The statistical technique introduced in the paper is also appropriate for the searches of sources of ultra-high energy cosmic rays.     

Modulation of galactic cosmic rays in a north–south asymmetrical heliosphere

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock. Voyager 1 crossed this termination shock at ∼94 AU in 2004, while Voyager 2 crossed it in 2007 at a different heliolatitude, about 10 AU closer to the Sun. These different positions of the termination shock confirm the dynamic and cyclic nature of the shock’s position. Observations from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that apart from the dynamic nature caused by changing solar activity there also may exist a global asymmetry in the north–south (polar) dimensions of the heliosphere, in addition to the expected nose–tail asymmetry. This relates to the direction in which the heliosphere is moving in interstellar space and its orientation with respect to the interstellar magnetic field. In this paper we focus on illustrating the effects of this north–south asymmetry on modulation of galactic cosmic ray Carbon, between polar angles of 55° and 125°, using a numerical model which includes all four major modulation processes, the termination shock and the heliosheath. This asymmetry is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. When comparing the computed spectra between the two polar angles, we find that at energies E < ∼1.0 GeV the effects of the assumed asymmetry on the modulated spectra are insignificant up to 60 AU from the Sun but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > ∼1.0 GeV, these effects remain insignificant throughout the heliosphere even very close to the heliopause. Furthermore, we find that a higher local interstellar spectrum for Carbon enhances the effects of asymmetric modulation between the two polar angles at lower energies (E < ∼300 MeV). In conclusion, it is found that north–south asymmetrical effects on the modulation of cosmic ray Carbon depend strongly on the extent of the geometrical asymmetry of the heliosheath together with the assumed value of the local interstellar spectrum.     

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.     

PHITS simulations of the Matroshka experiment

The radiation environment in space is very different from the one encountered on Earth. In addition to the sparsely ionizing radiation, there are particles of different Z with energies ranging from keV up to hundreds of GeV which can cause severe damage to both electronics and humans. It is therefore important to understand the interactions of these highly ionizing particles with different materials such as the hull of space vehicles, human organs and electronics. We have used the Particle and Heavy-Ion Transport code System (PHITS), which is a three-dimensional Monte Carlo code able to calculate interactions and transport of particles and heavy ions with energies up to 100 GeV/nucleon in most matter. PHITS is developed and maintained by a collaboration between RIST (Research Organization for Information Science & Technology), JAEA (Japan Atomic Energy Agency), KEK (High Energy Accelerator Research Organization), Japan and Chalmers University of Technology, Sweden. For the purpose of examining the applicability of PHITS to the shielding design we have simulated the ESA facility Matroshka (MTR) designed and lead by the German Aerospace Center (DLR). Preliminary results are presented and discussed in this paper.     

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.     

R3D-B2 – Measurement of ionizing and solar radiation in open space in the BIOPAN 5 facility outside the FOTON M2 satellite

Solar and space radiation have been monitored using the R3D-B2 radiation risks radiometer-dosimeter on board a recent space flight on the Russian satellite Foton M2 within the ESA Biopan 5 facility mounted on the outside of the satellite exposed to space conditions. The solar radiation has been assayed in four wavelength bands (UV-C, 170–280 nm, UV-B, 280–315 nm), UV-A (315–400 nm) and PAR (photosynthetic active radiation, 400–700 nm). The data show an increasing tumbling rotation of the satellite during the mission. The photodiodes do not show a cosine response to the incident light which has been corrected. After calibration of the signals using the extraterrestrial spectrum, doses have been calculated for each orbit, for each day and for the total mission as basic data for the biological material which has been exposed in parallel in the Biopan facility. Cosmic ionizing radiation has been monitored and separated in 256 deposited energy spectra, which were further used for determination of the absorbed dose rate and flux. Basic data tables were prepared to be used by other Biopan 5 experiments. The paper summarizes the results for the Earth radiation environment at the altitude (262–304 km) of the Foton M2 spacecraft. Comparisons with the predictions of NASA Earth radiation environment experimental models AE-8 and AP-8, and the PSB97 model are also presented, which calculate the fluxes of ionizing radiation from a simulation. AP-8 is a model for trapped radiation.     

Relative ion Fe,C and O abundances in quiet time particle fluxes in the 23 SC

Using ACE and SOHO data the origin of quiet-time low-energy particle fluxes at 1 AU is studied in the 23rd solar cycle. One of the selection criteria of quiet-time periods is to demand that H/He < 10 provided that periods with noticeable contribution of remnants of gradual events have been excluded from consideration. Our results suggest different origin of 0.03–3 MeV/nucleon particles – different seed populations accelerated and different acceleration processes. During the ascending, maximum and descending phases of solar activity quiet-time ions consist of coronal particles accelerated to suprathermal energies in about a half of the quiet periods, the rest of quiet-time fluxes originates from particle acceleration in processes similar to those in small impulsive solar flares rich in Fe. At solar minimum the bulk solar wind particles serve as seed population.     

Status of the GAMMA-400 project

The preliminary design of the new space gamma-ray telescope GAMMA-400 for the energy range 100 MeV–3 TeV is presented. The angular resolution of the instrument, 1–2° at E_γ ∼ 100 MeV and ∼0.01° at E_γ > 100 GeV, its energy resolution ∼1% at E_γ > 100 GeV, and the proton rejection factor ∼10⁶ are optimized to address a broad range of science topics, such as search for signatures of dark matter, studies of Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission, gamma-ray bursts, as well as high-precision measurements of spectra of cosmic-ray electrons, positrons, and nuclei.     

High energy charged particle experiment – A payload proposal to KuaFu-B mission

High Energy Charged Particle Experiment (HECPE) is to measure the fluxes of MeV electrons and tens of MeV protons. The two satellites of KuaFu-B are in the same polar orbit with apogee 7.0R_E, perigee 1.8R_E. They can sweep large L values and pass through the inner and outer radiation belts. The high energy electrons and protons in the radiation belts are principal sources for failures of satellites and spacecrafts in the Earth orbits. The enhancements of the high energy electrons and protons, so-called energetic particle events, are important phenomena of the Space Weather. The energy ranges monitored by HECPE are 0.3–0.5 MeV, 0.5–1.0 MeV, 1.0–2.0 MeV, and E > 2.0 MeV for electrons, 5–10 MeV, 10–20 MeV, 20–40 MeV, and 40–80 MeV for protons.     

Simulation of depth–dose distributions for various ions in polyethylene medium

Study of depth–dose distributions for intermediate energy ion beams in tissue-like media such as polyethylene (CH₂)_n provides a good platform for further improvements in the fields of hadrontherapy and space radiation shielding. The depth–dose distributions for ¹²C ions at various energies and for light and intermediate ion beams (³He, ¹⁶O, ²⁰Ne and ²⁸Si) as well as for heavy ions ⁵⁶Fe in polyethylene were estimated by using simulation toolkit: Geant4. Calculations were performed mainly by considering two different combinations of standard electromagnetic (EM), binary cascade (BIC), statistical multifragmentation (SMF) and Fermi breakup (FB) models. The energies of the ion beams were selected to achieve the Bragg peaks at predefined position (∼60 mm) and as per their availability. Variations of peak-to-entrance ratio (from 7.44 ± 0.05 to 8.87 ± 0.05), entrance dose (from 2.89 ± 0.01 to 203.71 ± 0.63 MeV/mm) and entrance stopping power (from 3.608 to 208.858 MeV/mm, calculated by SRIM) with atomic number (Z) were presented in a systematic manner. The better peak-to-entrance ratio and less entrance dose in the region Z = 2 to 8 (i.e. ³He to ¹⁶O) may provide the suitability of the ion beams for hadrontherapy.     

A new modeling method of solar energetic proton events for ISO specification

Solar energetic protons degrade performance and reliability of spacecraft systems due to single-event effects, total dose effects and displacement damage in electronics components including solar cells. On designing a solar cell panel, a total fluence of solar energetic protons (SEPs) which cause solar cell damage is needed to estimate power loss of the solar cells over mission life. Nowadays a solar panel area of spacecraft is increasing as spacecraft mission life becomes longer (15–18 years). Thus an accurate SEP model is strongly required for the cost-minimum design from the aerospace industry. The SEP model, JPL-91 proposed by Feynman et al., is currently used widely for solar cell designing. However, it is known that the JPL-91 model predicts higher fluences of protons than values actually experienced in space, especially after 7 years on orbit. In addition, the model is based on several assumptions, and also needs Monte-Carlo simulations for calculating fluences. In this study, we propose a new method for modeling SEPs especially focused on solar cell degradation. The newly-proposed method is empirical, which constructs a model based directly upon proton flux measurement data taken by instruments onboard spacecraft. This method has neither assumptions nor dependence on SEP event selection, both of which are needed in JPL-91. The model fluences derived from this method show lower fluences in longer missions compared to JPL-91. This method has been proposed to ISO (International Organization for Standardization) and has been discussed for a new standard SEP model.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Correlation between the cosmic noise absorption calculated from the SARINET data and the energetic particles measured by MEPED: Simultaneous observations over SAMA region

The cosmic noise absorption is presented in terms of two-dimensional images obtained from the imaging riometers operated at the Southern Space Observatory (geographic coordinate: 29.4° S, 53.1° W), in São Martinho da Serra, Brazil, Concepcion (geographic coordinate: 36.5° S, 73.0° W) and Punta Arenas (geographic coordinate: 53.0° S, 70.5° W) in Chile, which belong to the South American Riometer Network and are located at the central and periphery regions of the South American Magnetic Anomaly. Correlations are performed between the maximum cosmic noise absorption observed at these stations and the energetic electron flux in two energy channels (>30 and >300 keV) and the proton flux in three energy channels (80–240, 800–2500 and >6900 keV) as measured by the Medium Energy Proton and Electron Detector, during a moderate geomagnetic storm that occurred on September 3, 2008. The results show high correlations between the cosmic noise absorption detected at São Martinho da Serra and the flux of protons with energy between 80 and 240 keV, and the flux of electrons with energies higher than 300 keV, while an additional ionization at Concepcion was correlated with electrons of energies higher than 30 keV. The cosmic noise absorption detected at Punta Arenas was probably caused by the increase of the protons flux with energy between 80 and 240 keV.     

Space Shuttle drops down the SAA doses on ISS

Long-term analysis of data from two radiation detection instruments on the International Space Station (ISS) shows that the docking of the Space Shuttle drops down the measured dose rates in the region of the South Atlantic Anomaly (SAA) by a factor of 1.5–3. Measurements either by the R3DE detector, which is outside the ISS at the EuTEF facility on the Columbus module behind a shielding of less than 0.45 g cm⁻², and by the three detectors of the Liulin-5 particle telescope, which is inside the Russian PEARS module in the spherical tissue equivalent phantom behind much heavier shielding demonstrate that effect. Simultaneously the estimated averaged incident energies of the incoming protons rise up from about 30 to 45 MeV. The effect is explained by the additional shielding against the SAA 30–150 MeV protons, provided by the 78 tons Shuttle to the instruments inside and outside of the ISS. An additional reason is the ISS attitude change (performed for the Shuttle docking) leading to decreasing of dose rates in two of Liulin-5 detectors because of the East–West proton fluxes asymmetry in SAA. The Galactic Cosmic Rays dose rates are practically not affected.     

High energy radiation from the direction of the galactic black hole Sgr A

X-ray observations indicate that the Galactic black hole Sgr A^∗ is inactive now, however, we suggest that Sgr A^∗ can become active when a captured star is tidally disrupted and matter is accreted into the black hole. Consequently the Galactic black hole could be a powerful source of relativistic protons with a characteristic energy ∼10⁵² erg per capture. The diffuse GeV and TeV γ-rays emitted in the direction of the Galactic Center (GC) are the direct consequences of p–p collisions of such relativistic protons ejected by very recent capture events occurred ?10⁵ yr ago. On the other hand, the extended electron-positron annihilation line emission observed from GC is a phenomenon related to a large population of thermalized positrons, which are produced, cooled down and accumulated through hundreds of past capture events during a period of ∼10⁷ yr. In addition to explaining GeV, TeV and 511 keV annihilation emissions we also estimate the photon flux of several MeV resulting from in-flight annihilation process.