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.     

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.     

A study on emission geometry and the separation of peaks in the γ-ray pulse profiles of pulsar

In the inner annular gap (IAG) model the γ-ray radiation sources are suggested to be located close to the null charge surface (NCS). A method to explore two geometric parameters (κ, λ) of the IAG model is given in this paper. Pulsar population statistics are proposed to test the radiation locations of γ-ray pulsars. Within the IAG model, predictions for the GLAST observation of γ-ray pulsars are also made.     

Gamma-ray observations of explosive nucleosynthesis products

In this review I discuss the various γ-ray emission lines that can be expected and, in some cases have been observed, from radioactive explosive nucleosynthesis products. The most important γ-ray lines result from the decay chains of ⁵⁶Ni, ⁵⁷Ni, and ⁴⁴Ti. ⁵⁶Ni is the prime explosive nucleosynthesis product of Type Ia supernovae, and its decay determines to a large extent the Type Ia light curves. ⁵⁶Ni is also a product of core-collapse supernovae, and in fact, γ-ray line emission from its daughter product, ⁵⁶Co, has been detected from SN1987A by several instruments. The early occurrence of this emission was surprising and indicates that some fraction of ⁵⁶Ni, which is synthesized in the innermost supernova layers, must have mixed with the outermost supernova ejecta.Special attention is given to the γ-ray line emission of the decay chain of ⁴⁴Ti (⁴⁴Ti → ⁴⁴Sc → ⁴⁴Ca), which is accompanied by line emission at 68, 78, and 1157 keV. As the decay time of ⁴⁴Ti is ∼86 yr, one expects this line emission from young supernova remnants. Although the ⁴⁴Ti yield (typically 10⁻⁵–10⁻⁴M_⊙) is not very high, its production is very sensitive to the energetics and asymmetries of the supernova explosion, and to the mass cut, which defines the mass of the stellar remnant. This makes ⁴⁴Ti an ideal tool to study the inner layers of the supernova explosion. This is of particular interest in light of observational evidence for asymmetric supernova explosions.The γ-ray line emission from ⁴⁴Ti has so far only been detected from the supernova remnant Cas A. I discuss these detections, which were made by COMPTEL (the 1157 keV line) and BeppoSAX (the 68 and 78 keV lines), which, combined, give a flux of (2.6 ± 0.4 ± 0.5) × 10⁻⁵ ph cm⁻² s⁻¹ per line, suggesting a ⁴⁴Ti yield of (1.5 ± 1.0) × 10⁻⁴M_⊙. Moreover, I present some preliminary results of Cas A observations by INTEGRAL, which so far has yielded a 3σ detection of the 68 keV line with the ISGRI instrument with a flux that is consistent with the BeppoSAX detections. Future observations by INTEGRAL-ISGRI will be able to constrain the continuum flux above 90 keV, as the uncertainty about the continuum shape, is the main source of systematic error for the 68 and 78 keV line flux measurements. Moreover, with the INTEGRAL-SPI instrument it will be possible to measure or constrain the line broadening of the 1157 keV line. A preliminary analysis of the available data indicates that narrow line emission (i.e., Δv < 1000 km s⁻¹) can be almost excluded at the 2σ level, for an assumed line flux of 1.9 × 10⁻⁵ ph cm⁻² s⁻¹.     

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.     

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.     

What gamma-ray deexcitation lines reveal about solar-flares

Gamma-ray emission from solar flares reveals information about the nature of the accelerated particles and about the physical conditions of the medium through which the accelerated particles are transported. In this paper, we present the gamma-ray line-production and loop transport models used in our calculations of high-energy emission. We discuss the calculated interaction time history, the depth distribution, the interacting-particle angular distribution, and fluence ratios of the narrow gamma-ray lines. We show the relationship between the γ-ray observables and the parameters of the transport and line-production models. For illustration, we use calculations of 4.44 MeV ¹²C nuclear deexcitation line-production. Applications of the calculations to flare observations by both SMM and RHESSI are also presented.     

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.     

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.     

Study of >100 GeV electrons with BETS detector using a long duration balloon flight in Antarctica

We have observed cosmic-ray electrons from 10 to 1000 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates. The geometrical factor of detector is about 600 cm²sr and the total thickness of lead absorber is 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. We have collected 5.7 × 10³ events over 100 GeV including nearly 100 candidates of primary electrons.     

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.     

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.     

Energetic particle transport in anisotropic magnetic turbulence

We study energetic particle transport in a magnetic field configuration which models the solar wind magnetic turbulence plus the background field. A power-law Fourier amplitude is used for the fully 3D turbulence model, and in order to model anisotropic turbulence, the constant amplitude surfaces in k space are ellipsoids. The turbulence correlation lengths parallel (perpendicular) to the background magnetic field l_∥ (l_⊥) are varied in a wide range, and proton energies from 1 MeV to 10 GeV are assumed. Considering propagation on a distance corresponding to 1 AU, it is found that transport parallel and perpendicular to the background field heavily depends on the turbulence anisotropy, that is on the ratio l_∥/l_⊥. The spatial distribution of energetic particle follows the shape of magnetic flux tube up to about 10 MeV, while for larger energies the structure of the magnetic flux tube is progressively washed out. The scatterplots of particle distribution show intermittent, non Gaussian structures for l_∥ ≪ l_⊥ (quasi slab turbulence), while a more diffusive, Gaussian structure is obtained for l_∥ ≫ l_⊥ (quasi 2D turbulence). The long time behavior of transport shows that anomalous (subdiffusive perpendicular and superdiffusive parallel) transport regimes are obtained for l_∥ ≪ l_⊥, while Gaussian diffusive transport is obtained for both l_∥ ≫ l_⊥ and the isotropic turbulence case.     

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.     

VHE Gamma-ray supernova remnants

Increasing observational evidence gathered especially in X-rays and γ-rays during the course of the last few years support the notion that Supernova remnants (SNRs) are Galactic particle accelerators up to energies close to the “knee” in the energy spectrum of Cosmic rays. This review summarises the current status of γ-ray observations of SNRs. Shell-type as well as plerionic type SNRs are addressed and prospect for observations of these two source classes with the upcoming GLAST satellite in the energy regime above 100 MeV are given.     

Highly significant detection of solar neutrons on 2005 September 7

We have successfully detected solar neutrons at ground level in association with the X17.0 solar flare that occurred on 2005 September 7. Observations were made with the solar neutron telescopes and neutron monitors located in Bolivia and Mexico. In this flare, large fluxes of hard X-rays and γ-rays were observed by the GEOTAIL and the INTEGRAL satellites. The INTEGRAL observations include the 4.4 MeV line γ-rays of ¹²C. The data suggest that solar neutrons were produced at the same time as these hard electromagnetic radiations. We have however found an apparent discrepancy between the observed and the expected time profiles. This fact suggests a possible extended neutron emission.     

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

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.