Effect of random nature of cosmic ray sources – Supernova remnants – on cosmic ray intensity fluctuations,anisotropy, and electron energy spectrum

The random nature of sources (the supernova remnants) leads to the fluctuations of cosmic ray intensity in space and time. We calculate the expected fluctuations in a flat-halo diffusion model for particles with energies from 0.1 to 10³ TeV. The data on energy spectra and anisotropy of very high energy protons, nuclei and electrons, and the astronomical data on supernova remnants, the potential sources of cosmic rays, are used to constrain the value of the cosmic-ray diffusion coefficient and its dependence on energy.     

On the production of highest energy solar protons at 20 January 2005

On January 20, 2005, 7:02–7:05 UT the Aragats Multidirectional Muon Monitor (AMMM) located at 3200 m a.s.l. registered enhancement of the high energy secondary muon flux (threshold ∼5 GeV). The enhancement, lasting for 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES, and very fast (>2500 km/s) CME seen by SOHO, and the Ground Level Enhancements (GLE) #69 detected by the world-wide network of neutron monitors and muon detectors. The energetic and temporal characteristics of the muon signal from the AMMM are compared with the characteristics of other monitors located at the Aragats Space-Environmental Center (ASEC) and with other neutron and muon detectors. Since secondary muons with energies >5 GeV are corresponding to solar proton primaries with energies 20–30 GeV we conclude that in the episode of the particle acceleration at 7:02–7:05 UT 20 January 2005 solar protons were accelerated up to energies in excess of 20 GeV.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.     

Flux upper limits of diffuse TeV gamma rays from the Galactic plane using the effective area of the Tibet-II and -III arrays

We obtained new upper limits on the diffuse gamma rays from the inner Galactic (IG) and outer Galactic (OG) planes in 3–10 TeV region, using the Tibet air shower data and new Monte Carlo simulation results. A difference of the effective area of the air-shower array for observing gamma rays and cosmic rays was carefully taken into account in this analysis, resulting in that the flux upper limits of the diffuse TeV gamma rays were reduced by factors of 4.0–3.7 for 3–10 TeV than those in our previous results (Amenomori, M., Ayabe, S., Cui, S.W., et al. Observation of multi-TeV diffuse gamma rays from the Galactic plane with the Tibet air shower array. Astrophys. J. 580, 887–895, 2002.). This new result suggests that the inverse power index of the energy spectrum of source electrons responsible for generating diffuse TeV gamma rays through inverse Compton effect should be steeper than 2.2 and 2.1 for IG and OG planes, respectively, with 99%C.L.     

Evaluation of solar proton spectra using balloon cosmic ray observations and Monte Carlo simulation results

We have developed a method to evaluate the spectrum of solar energetic protons at the top of the Earth’s atmosphere from the measurements of our balloon cosmic ray experiment. By using the Monte Carlo PLANETOCOSMICS code based on Geant4 we compute the interaction of solar protons [10 MeV–10 GeV] with the Earth’s atmosphere. We obtain the angular and energy distributions of secondary particles (p, e⁻, e⁺, photons, muons) at different atmospheric levels as a function of primary proton spectra. By comparing the calculated depth dependence of the particle flux with the data obtained by our balloon experiment we can deduce the parameters of the solar proton spectrum that best fit the observations. In this paper we discuss our solar proton spectrum estimation method, and present results of its application to selected solar proton events from 2001 to 2005.     

Propagation of cosmic-ray electrons in the Galaxy

We formulate the global propagation model of cosmic-ray electrons including the source region, which is currently considered to be supernova remnants (SNRs). The model is characterized by the escape rate of electrons from SNRs into the interstellar space. It becomes clear that the energy index of the escape rate influences the high energy side of the interstellar spectrum and makes it possible to explain the observed data up to 2 TeV in the case of source spectral index smaller than 2.2 that is expected from the radio spectrum in SNRs. The escape lifetime of electrons in SNRs is also discussed by using the ratio of the radio flux in two regions: SNRs and the Galaxy. The result shows the mean lifetime in SNRs of ∼10⁴ yr around 1 GeV, which corresponds to the SNR age in the Sedov phase.     

Prospects for observing dark-matter remnants with GLAST

The satellite-based experiment, GLAST (Gamma-ray Large Area Space Telescope), is under construction and is planned to measure the cosmic γ-ray flux in the energy range 20 MeV to >300 GeV, with supporting measurements for γ-ray bursts from 10 keV to 25 MeV. With its launch in 2007, GLAST will open a new and important window on a wide variety of high-energy phenomena, including exotic relics from the Big Bang. Among these may be the decay/annihilation products of the hypothesized super symmetric image of the known particles. Single-photon energy thresholds for channels leading to such final states have been excluded in a model-dependent manner by accelerator searches to energies greater than 50 GeV. The ability of GLAST to set limits on this important component of cosmological evolution is presented along with an update on the present status of this mission.     

On the quality of mutations in mammalian cells induced by high LET radiations

The deleterious effects of accelerated heavy ions as component of the space radiation environment on living cells are of increasing importance for long duration human space flight activities. The most important aspect of such densely ionizing particle radiation is attributed to the type and quality of biological damage induced by them. This issue is addressed by investigating cell inactivation and mutation induction at the Hprt locus (coding for hypoxanthine-guanine-phosphoribosyl-transferase) of cultured V79 Chinese hamster cells exposed to densely ionizing radiation (accelerated heavy ions with different LETs from oxygen to gold, specific energies ranging from 1.9 to 69.7 MeV/u, corresponding LET values range from 62 to 13,223 keV/μm) and to sparsely ionizing radiation (200 kV X-rays). 30 spontaneous, 40 X-ray induced and 196 heavy ion induced 6-thioguanine resistant Hprt mutant colonies were characterized by Southern technique using the restriction enzymes EcoRI, PstI and BglII and a full length Hprt cDNA probe isolated from the plasmid pHPT12. Restriction patterns of the spontaneous Hprt mutants were indistinguishable from the wild type pattern, as these mutants probably contain only small deletions or even point mutations in the Hprt locus. In contrast, the overall spectrum of heavy ion induced mutations revealed a majority of partial or total deletions of the Hprt gene. With constant particle fluence (3 × 10⁶ particles/cm²) the quality of heavy ion induced mutations in the Hprt locus depends on physical parameters of the beam (atomic number, specific energy, LET). This finding suggests a relationship between the type of DNA damage and track structure. The fraction of mutants with severe deletions in the Hprt locus after exposure to oxygen ions increases from 65% at 60 keV/μm up to a maximum (100%) at 300 keV/μm and declines with higher LET values to 75% at 750 keV/μm. With heavier ions (Ca- and Au-ions) and even higher LET-values this mutant fraction decreases to 58% at 13,200 keV/μm. Heavy ion induced DNA break points in the Hprt locus are not randomly distributed.     

Modelling of protons spectra encountered in space using medical accelerator and its microdosimetric characterization

Radiation environments in space are mainly composed of protons coming from the Galactic Cosmic Rays (GCRs) pervading the universe, the Solar Particle Events (SPEs) resulting from solar flares and coronal mass ejections, and the two Van Allen Belts surrounding the Earth due to the presence of the geomagnetic field trapping charged particles. Their wide spectra of energies up to hundreds of GeV imply diverse radiobiological effects to astronauts and radiation damage to electronics in the spacecraft. Even if lower in abundance, heavy ions such as He, C, O, Si, Fe are present in space and constitute an even bigger hazard due to their high penetrability and high linear energy transfer (LET). Most irradiation facilities available for research and testing worldwide provide usually only monoenergetic beams of high-energy protons or other heavier particles limiting studies of radiobiological effects and effects on electronics to a set of discrete energies.This paper introduces a procedure where a proton fluence spectra of interest for space radiation protection, previously generated by Monte Carlo simulations was delivered using a clinical proton therapy accelerator. Particularly, it reports the first results of modelling a proton radiation field in space in the energy range from 70 to 230 MeV during a single experimental session by programming a treatment planning system (TPS) to deliver required proton irradiation energies. Moreover, the angular distribution of the proton irradiation field has been varied to reproduce the isotropic exposure experienced by humans in space. The obtained proton radiation field was characterized using a 3D sensitive volume SOI microdosimeter developed by the Centre for Medical Radiation Physics (CMRP), University of Wollongong, Australia.     

Evidence for a solar coronal thick-target hard X-ray source observed by RHESSI

We study a solar flare hard X-ray (HXR) source observed by the Reuven Ramaty high energy solar spectroscopic imager (RHESSI) in which the HXR emission is almost entirely in a coronal loop so dense as to be collisionally thick at electron energies up to ∼45−60 keV. This contrasts with most events previously reported in which the HXR emission is primarily from the loop footpoints in the collisionally dense chromosphere. In particular, we show that the high loop column densities inferred from the GOES and RHESSI soft X-ray emission measure and the volume of the flare loop are consistent with the coronal thick-target interpretation of the HXR images and spectra. The high column densities observed already at the very beginning of the impulsive phase are explained by chromospheric evaporation during a preflare which, as Nobeyama 17 GHz radio images reveal, took place in the same set of nested loops as the main flare.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.     

Particle acceleration and transport at an oblique CME-driven shock

In gradual solar energetic particle (SEP) events, protons and heavy ions are often accelerated to >100 MeV/nucleon at a CME-driven shock. In this work, we study particle acceleration at an oblique shock by extending our earlier particle acceleration and transport in heliosphere (PATH) code to include shocks with arbitrary θ_BN, where θ_BN is the angle between the upstream magnetic field and the shock normal. Instantaneous particle spectra at the shock front are obtained by solving the transport equation using the total diffusion coefficient κ, which is a function of the parallel diffusion coefficient κ_∥ and the perpendicular diffusion coefficient κ_⊥. In computing κ_∥ and κ_⊥, we use analytic expressions derived previously. The particle maximum energy at the shock front as a function of time, the time intensity profiles and particle spectra at 1 AU for five θ_BN’s are calculated for an example shock.     

Measurement of turbulent kinetic energy dissipation rates in the mesosphere by a 3 MHz Doppler radar

A new narrow beam Doppler radar operating at 3.17 MHz has been installed close to the Andøya Rocket Range in Andenes, Norway in summer 2002 in order to improve the ground based capabilities for measurements of turbulence in the mesosphere. The main feature of the radar is a Mills Cross transmitting/receiving antenna consisting of 29 crossed half-wave dipoles. In combination with the modular transceiver system this provides high flexibility in beam forming and pointing. In general, vertical and oblique beams with a minimum one way half-power full-beam width (HPFW) of 6.6° are used. The observations are usually performed with a height resolution of 1 km and with off-zenith beams at 7.3° directed towards NW, NE, SE, and SW. Turbulence intensities have been estimated from the width of the observed signal spectra using an computationally intensive correction method which requires precise knowledge of the antenna radiation pattern. The program uses real-time measurements of the wind field in all determinations. Turbulent kinetic energy dissipation rates based on radar observations are presented and compared with corresponding climatological summer and winter profiles from rocket measurements, as well as with single profiles from model runs for selected periods from September 2003 to Summer 2004. The mean turbulent kinetic energy dissipation rates based on these radar measurements are about 5 mW/kg at 60 km altitude and about 20 mW/kg at 80 km, in reasonable agreement with mean turbulence intensities obtained from previous rocket soundings at Andenes.     

Detection of meteors and sub-relativistic dust grains by the fluorescence detectors of ultra high energy cosmic rays

Fluorescence detectors of ultra high energy cosmic rays (UHECR) allow to record not only the extensive air showers, initiated by the UHECR particles, but also to detect light, produced by meteors and by the fast dust grains. It is shown that the fluorescence detector operated at the mountain site can register signals from meteors with kinetic energy threshold of about 25 J (meteor mass ∼ 5 × 10⁻⁶ g, velocity ∼ 3 × 10⁶ cm/s). The same detector might be used for recording of the dust grains of smaller mass (of about 10⁻¹⁰ g) but with velocity 10⁹ cm/s, close to the light velocity (sub-relativistic dust grains). The light signal from a sub-relativistic dust grain is expected in much shorter time scale (∼0.001 s), in comparison with the meteor signal (∼0.1–1 s), and much longer than duration of the UHECR signals (tens of μs). The fluorescence detector capable to register various phenomena: from meteors to UHECR – should have a variable pixel and selecting system integration time. A study of the new phenomenon of sub-relativistic grains will help to understand the mechanism of particle and dust grain acceleration in astrophysical objects (in SN explosions, for example).     

Hard X-ray footpoint motions in solar flares: Comparing magnetic reconnection models with observations

Hard X-ray observations from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of the October 29, 2003 GOES X10 two-ribbon flare are used together with magnetic field observations from the Michelson Doppler Imager (MDI) onboard SoHO to compare footpoint motions with predictions from magnetic reconnection models. The temporal variations of the velocity v of the hard X-ray footpoint motions and the photospheric magnetic field strength B in footpoints are investigated. The underlying photospheric magnetic field strength is generally higher (B ∼ 700–1200 G) in the slower moving (v ∼ 20–50 km s⁻¹) western footpoint than in the faster (v ∼ 20–100 km s⁻¹) moving eastern source (∼100–600 G). Furthermore, a rough temporal correlation between the HXR flux and the product vB² is observed.     

Direct measurement of the composition and spectra of cosmic rays above 1 TeV/amu from Jacee

Direct measurements on cosmic ray protons through iron above about 1 TeV/amu have been performed in a series of balloon-borne experiments with emulsion chambers. The measured energy spectra of protons and helium are power laws with exponents of 2.77 ± 0.09 and 2.72 ± 0.11 in the energy range 5 to 500 TeV and 2 to 50 TeV/amu, respectively. The proton spectrum shows no evidence of the steepening near 2 TeV which was reported by other experiments. Helium has a slightly higher intensity compared to extrapolations from lower energy measurements. The heavier elements, carbon to sulfur, show a small tendency for intensity enhancement in the relative abundance above 10 TeV/amu.     

Giant radio pulses from the Crab pulsar

Individual giant radio pulses (GRPs) from the Crab pulsar last only a few microseconds. However, during that time they rank among the brightest objects in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio frequencies. Our observations show that GRPs can be found in all phases of ordinary radio emission including the two high frequency components (HFCs) visible only between 5 and 9 GHz [Moffett, D.A., Hankins, T.H. Multifrequency radio observations of the Crab pulsar. Astrophys. J. 468, 779–783, 1996]. This leads us to believe that there is no difference in the emission mechanism of the main pulse (MP), inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent observations of giant pulses with the Effelsberg telescope at a center frequency of 8.35 GHz show distinct spectral maxima within our observational bandwidth of 500 MHz for individual pulses. Their narrow band components appear to be brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, there is an evidence for spectral evolution within and between those structures. High frequency features occur earlier than low frequency ones. Strong plasma turbulence might be a feasible mechanism for the creation of the high energy densities of ∼6.7 × 10⁴ erg cm⁻³ and brightness temperatures of ∼10³¹ K.     

Storm time spatial variations in TEC during moderate geomagnetic storms in extremely low solar activity conditions (2007–2009) over Indian region

The total electron content (TEC) measurements from a network of GPS receivers were analyzed to investigate the storm time spatial response of ionosphere over the Indian longitude sector. The GPS receivers from the GPS Aided Geo Augmented Navigation (GAGAN) network which are uniquely located around the ∼77°E longitude are used in the present study so as to get the complete latitudinal coverage from the magnetic equator to low mid-latitude region. We have selected the most intense storms but of moderate intensity (−100 nT < Dst < −50 nT) which occurred during the unusually extremely low solar activity conditions in 2007–2009. Though the storms are of moderate intensity, their effects on equatorial to low mid-latitude ionosphere are found to be very severe as TEC deviations are more than 100% during all the storms studied. Interesting results in terms of spatial distribution of positive/negative effects during the main/early recovery phase of storms are noticed. The maximum effect was observed at crest region during two storms whereas another two storms had maximum effect near the low mid-latitude region. The associated mechanisms like equatorial electrodynamics and neutral dynamics are segregated and explained using the TIMED/GUVI and EEJ data during these storms. The TEC maps are generated to investigate the storm time development/inhibition of equatorial ionization anomaly (EIA).     

A detailed FLUKA-2005 Monte-Carlo simulation for the ATIC detector

We have performed a detailed Monte-Carlo (MC) simulation for the Advanced Thin Ionization Calorimeter (ATIC) detector using the MC code FLUKA-2005 which is capable of simulating particles up to 10 PeV. The ATIC detector has completed two successful balloon flights from McMurdo, Antarctica lasting a total of more than 35 days. ATIC is designed as a multiple, long duration balloon flight, investigation of the cosmic ray spectra from below 50 GeV to near 100 TeV total energy; using a fully active Bismuth Germanate (BGO) calorimeter. It is equipped with a large mosaic of silicon detector pixels capable of charge identification, and, for particle tracking, three projective layers of x–y scintillator hodoscopes, located above, in the middle and below a 0.75 nuclear interaction length graphite target. Our simulations are part of an analysis package of both nuclear (A) and energy dependences for different nuclei interacting in the ATIC detector. The MC simulates the response of different components of the detector such as the Si-matrix, the scintillator hodoscopes and the BGO calorimeter to various nuclei. We present comparisons of the FLUKA-2005 MC calculations with GEANT calculations and with the ATIC CERN data.     

Relationship of solar flare accelerated particles to solar energetic particles (SEPs) observed in the interplanetary medium

Observations of hard X-ray (HXR)/γ-ray continuum and γ-ray lines produced by energetic electrons and ions, respectively, colliding with the solar atmosphere, have shown that large solar flares can accelerate ions up to many GeV and electrons up to hundreds of MeV. Solar energetic particles (SEPs) are observed by spacecraft near 1 AU and by ground-based instrumentation to extend up to similar energies as in large SEP events, but it appears that a different acceleration process, one associated with fast coronal mass ejections is responsible. Much weaker SEP events are observed that are generally rich in electrons, ³He, and heavy elements. The energetic particles in these events appear to be similar to those accelerated in flares. The Ramaty high energy solar spectroscopic imager (RHESSI) mission provides high-resolution spectroscopy and imaging of flare HXRs and γ-rays. Such observations can provide information on the location, energy spectra, and composition of the flare accelerated energetic particles at the Sun. Here, preliminary comparisons of the RHESSI observations with observations of both energetic electron and ion near 1 AU are reviewed, and the implications for the particle acceleration and escape processes are discussed.