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.     

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.     

Progress in laboratory performance verification of the R-process Isotope Observer

The source of galactic cosmic rays (GCR) remains unknown, despite many decades of research. It is generally accepted among high-energy astrophysicists that GCR are accelerated by supernova (SN) shocks in the interstellar medium, but the evidence for this remains indirect. However, the fact that the cosmic-ray spectrum continues far past the limit of SN shock acceleration is a serious problem for this model. Most SNae occur in superbubbles – regions of the interstellar medium containing massive stars and SN remnants. Thus, SNae should be accelerating the ejecta of other SNae. The local medium of the superbubble will be enhanced in freshly synthesized R-process material. We are currently developing the R-process Isotope Observer (RIO), which could make the first measurements of the isotopic abundances of the “ultraheavy” GCR (those in the range 32 ⩽ Z ⩽ 42) in order to determine the fractional contribution of freshly synthesized R-process material in GCR through the measurement of several key isotopic ratios. We have recently performed a series of tests of the RIO detector with 400 A MeV ⁸⁴Kr at the HIMAC accelerator in Japan. The results of these tests indicate that a mass resolution of ∼0.3 u can be achieved with RIO. We also discuss the prospects for accommodation of the RIO mission.     

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).     

The energy spectra of protons and helium measured with the ATIC experiment

The Advanced Thin Ionization Calorimeter (ATIC) balloon experiment is designed to investigate the composition and energy spectra of cosmic rays at the highest energies currently accessible by direct measurements, i.e., the region up to 100 TeV. The instrument consists of a silicon matrix for charge measurement, a graphite target (0.75 nuclear interaction length) to induce hadronic interactions, three layers of scintillator strip hodoscopes for triggering and trajectory reconstruction, and a Bismuth Germanate (BGO) crystal calorimeter (18 radiation lengths) to measure particle energies. ATIC has had two successful Long Duration Balloon (LDB) flights from McMurdo, Antarctica: one from 12/28/00 to 01/13/01 and the other from 12/29/02 to 01/18/03. We present the energy spectra of protons and helium extracted from the first flight, over the energy range from 100 GeV to 100 TeV, and compare them with the results from other experiments at both the lower and higher energies. ATIC-1 results do not indicate significant differences in spectral shape between protons and helium over the investigated energy range.     

RIO: the R-process isotope observer

Our galaxy is filled with a “gas” of relativistic nuclei and electrons – the galactic cosmic rays (GCRs). The source of GCR nuclei is unknown, but there is a general consensus among high-energy astrophysicists that GCRs are accelerated by supernova (SN) shocks in the interstellar medium. The evidence in support of this picture is strong but indirect. However, the fact that the cosmic-ray spectrum extends continuously, without steps or peaks, to more than five orders of magnitude in energy beyond the limit of SN shock acceleration appears to be incompatible with this picture, at least without extreme fine-tuning of models. A “smoking gun” is needed to definitively establish that SN shocks are indeed the accelerator of GCRs. If GCRs are accelerated in SN shocks, they will be enhanced in freshly-synthesized r-process material. We are currently studying the R-process Isotope Observer (RIO) as a Mission of Opportunity for the International Space Station. RIO will make the first measurements of the isotopic abundances of the “ultraheavy” GCRs (those in the range 32 ⩽ Z ⩽ 42) and will determine the fractional contribution of freshly-synthesized r-process material in GCRs through the measurement of several key isotopic ratios.     

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.     

Scientific ballooning in Japan – An over view of recent activities

Current status of scientific ballooning in Japan is reviewed. First, I describe successful application of balloon technologies to construct a vessel of transparent plastic film, to contain about 1000 tons of liquid scintillator in Kamioka Liquid Scintillator Anti-Neutrino Detector (KamLAND). KamLAND is a project to study neutrino oscillation phenomena, by detecting anti-neutrinos produced in distant nuclear reactors. Next, I describe high altitude balloons developed by the ISAS balloon group. They developed balloons made from ultra-thin polyethylene film, producing a balloon of volume 60,000 m³ which successfully reached an altitude of 53 km in 2002. This is a world record, the greatest altitude that a balloon has ever achieved. ISAS is applying further effort to develop balloons with even thinner films, to achieve a higher altitude than 53 km. Other recent activities by the ISAS balloon group are briefly described.I also review scientific ballooning projects now operating in Japan, particularly focusing on the Balloon-Borne Experiment with a Superconducting Spectrometer (BESS) program. This is a US–Japan collaborative program that has carried out very precise measurements of antiprotons, protons and other components in primary cosmic rays, as well as measuring the fluxes of atmospheric muons and other components. The results of these observations give us important information to improve our understanding of the production mechanism of antiprotons observed in the primary cosmic rays. The data are also important for analysis of atmospheric neutrino events observed by Super-Kamiokande and other ground-based neutrino detectors. Future prospects of BESS and other balloon-borne cosmic-ray research programs are also presented.     

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.     

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.     

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.     

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.     

On the 27-day variations of the galactic cosmic ray anisotropy and intensity for different periods of solar magnetic cycle

The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy calculated based on the neutron monitors experimental data is larger in the qA > 0 period than in the qA < 0 period of solar magnetic cycle. The amplitudes of the 27-day variation of the galactic cosmic rays anisotropy do not depend on the tilt angles of the heliospheric neutral sheet for different the qA > 0 and the qA < 0 periods of solar magnetic cycle. A good correlation has been revealed between the changes of the amplitudes of the 27-day variations of the galactic cosmic ray anisotropy and intensity versus the qA > 0 and the qA < 0 periods of solar magnetic cycle.     

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.     

High energy emission from pulsars: Outer gap scenario

I will give a brief review of the recent development in the emission models of isolated, rapidly rotating neutron stars, focusing on the γ-ray radiation mechanism in their outer magnetospheres. By examining the Poisson equation for the electrostatic potential, I show that an active particle accelerator must extend from the vicinity of the neutron star surface to the vicinity of light cylinder. Furthermore, combining the Poisson equation with the Boltzmann equations for electrons/positrons and γ-rays, and assuming that the gap trans-field thickness is large compared to the longitudinal width, I demonstrate that the energy distribution of ultra-relativistic particles cannot be described by a power-law but by a quasi-monoenergetic distribution at the terminal Lorentz factor. The particles are accelerated in the gap and escape from it with large Lorentz factors. Is is shown that such energetic particles migrating outside of the gap contribute significantly to the γ-ray luminosity and reproduce the observed soft γ-ray spectrum between 100 MeV and 3 GeV for the Vela pulsar.     

The superbubble origin for galactic cosmic rays

The source of energy for cosmic-ray acceleration appears to be shock waves driven by supernova (SNe) ejecta. The great majority (80–90%) of SNe (SNII and SNIb) are formed by the core collapse of young, massive O and B stars. However, it has been known for more than forty years that the births of such massive stars in stellar clumps, termed OB associations, are correlated in space and time. The combined ejecta of core-collapse SNe, occurring at the deaths of these massive stars, create low-density (∼5.0 × 10⁻³ cm⁻³) superbubbles that reach dimensions of several hundred pc. The occurrence of correlated SNe in superbubbles affects not just the source of cosmic-ray energy, SNe shock waves, it impinges as well on the elemental and isotopic source abundances of cosmic-ray nuclei. We argue that the well-known anomalous cosmic-ray ²²Ne/²⁰Ne ratio, a factor of five times the Solar System ratio, results from a mixing of freshly synthesized nucleosynthetic material in supernova active cores of superbubbles. Although diluted by mixing with older, lower metallicity interstellar gas, the mean metallicities in the superbubble, SNe-dominated cores are high ∼3 times the Solar System value.     

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.     

Supernova remnant 1987A: The latest report from the Chandra X-ray Observatory

We continue monitoring supernova remnant (SNR) 1987A with the Chandra X-ray Observatory. As of 2004 January, bright X-ray spots in the northwest and the southwest are now evident in addition to the bright eastern ring. The overall X-ray spectrum, since 2002 December, can be described by a planar shock with an electron temperature of ∼2.1 keV. The soft X-ray flux is now 8 × 10⁻¹³ ergs cm⁻² s⁻¹, which is about five times higher than four years ago. This flux increase rate is consistent with our prediction based on an exponential density distribution along the radius of the SNR between the HII region and the inner ring. We still have no direct evidence of a central point source, and place an upper limit of L_X = 1.3 × 10³⁴ ergs s⁻¹ on the 3–10 keV band X-ray luminosity.     

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.     

A search for the signature of microquasars in the cosmic ray iron spectrum measured by TIGER

In a recent paper Heinz and Sunyaev suggest that relativistic jets observed in microquasars might result in narrow features in the energy spectra of heavy cosmic rays with ≈1 to ≈10 GeV/nuc. They further argue that such features might be observable if there has been one or more microquasars nearby within the last few million years. We report preliminary results of a search for evidence of such features using data from a 32-day balloon flight of the Trans-Iron Galactic Element Recorder (TIGER). Although this flight took place near solar maximum, calculations of the broadening effects of solar modulation indicate that a narrow feature of sufficient intensity should still be observable. An energy spectrum for iron with high statistical significance has been derived from ≈100,000 Fe events in the energy range from about 2.5 to 10 GeV/nuc. Although our preliminary results do not reveal any obvious features, we will discuss the possibility of observing such features with TIGER and other instruments.