Gamma-ray astronomy and cosmic ray origin problem

We analyse the results of the nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) in order to describe their relevant properties: the remnant dynamics and the characteristics of nonthermal emission produced by CRs. It is shown that the theory fits the existing data in a satisfactory way and that the magnetic field in SNRs is significantly amplified due to efficient acceleration of the nuclear CR component. From the fact that magnetic field amplification occurs in all the young SNRs for which relevant data exist, and given the strong theoretical connection between magnetic field amplification and efficiently accelerating the nuclear CRs, we tentatively conclude that the Galactic SNRs are the source population of the Galactic CRs. Due to high interior magnetic fields in young SNRs the π⁰-decay γ-rays generated by the nuclear CR component as a rule dominate over γ-rays generated by the electron CR component, and the calculated γ-ray flux fits existing data.     

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.     

Evolution of magnetic fields and cosmic ray acceleration in supernova remnants

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh–Taylor (R–T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an Adaptive Mesh Refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions.     

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.     

New evidence for strong non-thermal effects in Tycho’s supernova remnant

For the case of Tycho’s supernova remnant (SNR), we present the relation between the blast wave and contact discontinuity radii calculated within the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is demonstrated that these radii are confirmed by recently published Chandra measurements which show that the observed contact discontinuity radius is very close to the shock radius. Therefore a consistent explanation of these observations can be given in terms of efficient CR acceleration which makes the medium more compressible.     

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.     

Contributions from nearby pulsars to the local cosmic ray electron spectrum

PSR J0437-4715 is one of the closest millisecond pulsars (MSPs) to Earth, lying at a distance of ∼140 pc. This pulsar has a characteristic age of 4.9 Gyr and a relatively low spindown power of ∼10³⁴ ergs/s. During its rather long lifetime, a large fraction of the energy output has been in the form of multi-TeV electrons. In this paper, we investigate the possible contribution of this nearby MSP to the local interstellar electron spectrum (LIS). The old age of the system justifies a steady-state evaluation of the contribution from this pulsar to the LIS. We calculate the electron spectrum at the light cylinder in the framework of a General Relativistic polar cap (PC) model, and use this as an injection spectrum in a diffusion model. The younger Geminga pulsar is also very close to Earth and warrants investigation. A steady-state approach is however no longer justified, so we use an impulsive injection model. We will present results of a study of the contribution from these pulsars to the cosmic ray (CR) LIS. Our calculations show that pulsars like Geminga can make a non-negligible contribution to the LIS.     

The N4,5 - OO Auger and “N3” N4,5O2,3 Coster-Kronig spectra of xenon induced by electron impact

The N_4,5OO Auger electron spectrum in the kinetic energy region (14–37) eV has been measured with high resolution at electron incident energies between 71 eV and 2019 eV. The Auger and satellite lines are assigned by comparison with previous literature data. Several states which correspond to the resonant Auger transitions from the 4d_{3/2, 5/2} 6p states are observed. The PCI effects on the N₅-O_2,3O_2,3(¹S₀) Auger peak when the incident electron energy approaches to the N₅ edge (67.55 eV) is investigated and it is found that the dependence of the energy shift on the excess energy, E_exc, is well represented by the function E_exc^β with β = (-1.40 ± 0.05). The ^“N₃^” N_4,5O_2,3 Coster-Kronig spectrum is measured at various electron incident energies. The assignments of the features are made in comparison with a similar spectrum from synchrotron radiation measured by Kivimäki et al. [J. Electron Spectrosc. Relat. Phenom. 1999, 101–103, 43–47]. The Auger spectrum in the kinetic energy regions (68–78) eV and (90–138) eV is observed for the first time in an electron impact experiment and a comparison is made with the photoionization experiment by Hikosaka et al. [Phys. Rev. A 2007, 76, 032708].     

Propagation and origin of ultra high-energy cosmic rays

Propagation of UHE protons through CMB radiation leaves the imprint on energy spectrum in the form of Greisen–Zatsepin–Kuzmin (GZK) cutoff, bump (pile-up protons) and dip. The dip is a feature in energy range 1 × 10¹⁸–4 × 10¹⁹ eV, caused by electron–positron pair production on CMB photons. Calculated for power-law generation spectrum with index γ_g = 2.7, the shape of the dip is confirmed with high accuracy by data of Akeno-AGASA, HiRes, Yakutsk and Fly’s Eye detectors. The predicted shape of the dip is robust: it is valid for the rectilinear and diffusive propagation, for different discreteness in the source distribution, for local source overdensity, deficit, etc. This property of the dip allows us to use it for energy calibration of the detectors. The energy shift λ for each detector is determined by minimum χ² in comparison of observed and calculated dip. After this energy calibration the absolute fluxes, measured by AGASA, HiRes and Yakutsk detectors remarkably coincide in energy region 1 × 10¹⁸–1 × 10²⁰ eV. Below the characteristic energy E_c ≈ 1 × 10¹⁸ eV the spectrum of the dip flattens for both diffusive and rectilinear propagation, and more steep galactic spectrum becomes dominant at E < E_c. The energy of transition E_tr < E_c approximately coincides with the position of the second knee E_2kn, observed in the cosmic ray spectrum. The dip-induced transition from galactic to extragalactic cosmic rays at the second knee is compared with traditional model of transition at ankle, the feature observed at energy 1 × 10¹⁹ eV.     

Antiparticle content in the magnetosphere

In the present work we assess the stable and transient antiparticle content of planetary magnetospheres, and subsequently we consider their capture and application to high delta-v space propulsion. We estimate the total antiparticle mass contained within the Earth’s magnetosphere to assess the expediency of such usage. Using Earth’s magnetic field region as an example, we have considered the various source mechanisms that are applicable to a planetary magnetosphere, the confinement duration versus transport processes, and the antiparticle loss mechanisms. We have estimated the content of the trapped population of antiparticles magnetically confined following production in the exosphere due to nuclear interactions between high energy cosmic rays (CR) and constituents of the residual planetary upper atmosphere.The galactic antiprotons that directly penetrate into the Earth’s magnetosphere are themselves secondary by its nature, i.e. produced in nuclear reactions of the cosmic rays passing through the interstellar matter. These antiproton fluxes are modified, dependent on energy, when penetrating into the heliosphere and subsequently into planetary magnetospheres. During its lifetime in the Galaxy, CR pass through the small grammage of the interstellar matter where they produce secondary antiprotons. In contrast to this, antiprotons generated by the same CR in magnetosphere are locally produced at a path length of several tens g/cm² of matter in the ambient planetary upper atmosphere. Due to the latter process, the resulting magnetically confined fluxes significantly exceed the fluxes of the galactic antiprotons in the Earth’s vicinity by up to two orders of magnitude at some energies.The radiation belt antiparticles can possibly be extracted with an electromagnetic-based “scoop” device. The antiparticles could be concentrated by and then stored within the superimposed magnetic field structure of such a device. In future developments, it is anticipated that the energy of the captured antiparticles (both rest energy and kinetic energy) can be adapted for use as a fuel for propelling spacecraft to high velocities for remote solar system missions.     

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.     

New insights into SNR evolution revealed by the discovery of recombining plasmas

We report the discovery of recombining plasmas in three supernova remnants (SNRs) with the Suzaku X-ray astronomy satellite. During SNR’s evolution, the expanding supernova ejecta and the ambient matter are compressed and heated by the reverse and forward shocks to form an X-ray emitting hot plasma. Since ionization proceeds slowly compared to shock heating, most young or middle-aged SNRs have ionizing (underionized) plasmas. Owing to high sensitivity of Suzaku, however, we have detected radiative recombination continua (RRCs) from the SNRs IC 443, W49B, and G359.1–0.5. The presence of the strong RRC is the definitive evidence that the plasma is recombining (overionized). As a possible origin of the overionization, an interaction between the ejecta and dense circumstellar matter is proposed; the highly ionized gas was made at the initial phase of the SNR evolution in dense regions, and subsequent rapid adiabatic expansion caused sudden cooling of the electrons. The analysis on the full X-ray band spectrum of IC 443, which is newly presented in this paper, provides a consistent picture with this scenario. We also comment on the implications from the fact that all the SNRs having recombining plasmas are correlated with the mixed-morphology class.     

Peculiarities of galactic cosmic ray Forbush effects during October–November 2003

Features of two successive Forbush effects of the galactic cosmic ray intensity in October–November 2003 have been studied based on the neutron monitors data. The rigidity spectrum of the galactic cosmic ray intensity in the course of the first Forbush effect (22–27 October) is gradually hardening, while the rigidity spectrum of the second Forbush effect (28 October–10 November) from the starting moment is very hard. As far, the energy range of the turbulence of the interplanetary magnetic field is in general responsible for the diffusion of galactic cosmic ray particles of the energy 5–50 GeV (to which neutron monitors are sensitive), we postulate that the gradually hardening (from day to day) of the rigidity spectrum of the first Forbush effect is associated with the enhancement of the power spectral density in the energy range of the interplanetary magnetic field turbulence caused by the large scale irregularities generated due to the interaction of the extending high speed disturbances with the background solar wind. The very hard rigidity spectrum (from the starting moment) of the second Forbush effect is generally associated with the well established new structure of the energy range of the interplanetary magnetic field turbulence enriched by the already created large scale irregularities. The gradually softening of the rigidity spectrum during the recovery phase of the second Forbush effect confirms that the disturbed interplanetary magnetic field turbulence step by step returns to the initial state.     

AVS-F observations of γ-ray emission during January 20, 2005 solar flare up to 140 MeV

The solar flare of January 20, 2005 (X7.1, 06:36–07:26 UT, maximum at 07:01 UT by the GOES soft X-ray data) was the most powerful one in January 2005 series. The AVS-F apparatus onboard CORONAS-F registered γ-emission during soft X-ray rising phase of this flare in two energy ranges of 0.1–20 MeV and 2–140 MeV. The highest γ-ray energy registered during this flare was ∼140 MeV. Six spectral features were registered in energy spectrum of this solar flare: annihilation + αα (0.4–0.6 MeV), ²⁴Mg + ²⁰Ne + ²⁸Si + neutron capture (1.7–2.3 MeV), ²¹Ne + ²²Ne + ¹⁶O + ¹²С (3.2–5.0 MeV), ¹⁶O (5.3–6.9 MeV), one from neutral pions decay (25–110 MeV) and one in energy band 15–21 MeV. Four of them contain typical for solar flares lines – annihilation, nuclear de-excitation and neutron capture at ¹H. Spectral feature caused by neutral pions decay was registered during several flares too. Some spectral peculiarities in the region of 15–21 MeV were first observed in solar flare energy spectrum.     

Possible influence of cosmic rays on climate through thunderstorm clouds

With decreasing of cosmic ray (CR) intensity caused by increasing of solar activity (SA) or in some short periods of Forbush-decreases, the intensity of secondary CR relativistic electrons decreases and the probability of formation of thunderstorm clouds and discharges between clouds or between clouds and ground is also expected to decrease. This will influence on weather and climate. In this case is very important to have more detail information on the atmospheric electric field distribution in the atmosphere, additional to information what gave now electric field sensors (EFS) only in about one point near the ground. We show that CR not only influenced on atmospheric electric field phenomenon, but can give practically continuous information on the atmospheric electric field distribution in the atmosphere. We extend our theory of CR atmospheric electric field effect on electron–photon, muon and neutron component including different multiplicities. We take into account that about 0.07 of neutron monitor counting rate caused by negative soft muons captured by lead nucleons and formed mesoatoms with generation of several MeV energy neutrons from lead. In this case the neutron monitor or neutron super-monitor works as analyzer that detects muons of only one, negative sign. It is very important because the atmospheric electric field effect have opposite signs for positive and negative muons that main part of this effect in the muon telescope or in ionization chamber is compensated and we can observe only small part of total effect of one sign muons. On the basis of our general theory of CR atmospheric electric field effects with taking into account of negative soft muon acceleration and deceleration in the Earth atmosphere (in dependence of direction and intensity of electric field) we discuss the possibility of existing this effect in CR neutron monitor counting rate and in different multiplicities and calculate the expected effects in dependence of atmospheric electric field distribution in the atmosphere. We show that the comparison of observed effects with theoretically expected will give important information on the value of atmospheric electric field and its distribution in the atmosphere. We consider also the possible influence of secondary relativistic electrons of CR and relativistic electrons precipitated from the Earth’s radiation belts on thunderstorms and lightnings, and through this – on climate change.     

Modeling and data analysis of a Forbush decrease

We study the Forbush decrease of the galactic cosmic ray intensity observed in 9–25 September 2005 using the experimental data and a newly developed time-dependent three dimensional modeling. We analyze neutron monitors and muon telescopes, and the interplanetary magnetic field data. We demonstrate a clear relationship between the rigidity (R) spectrum exponent (γ) of the Forbush decrease and the exponent (ν) of the power spectral density of the components of the interplanetary magnetic field in the frequency range of ∼ 10⁻⁶–10 ⁻⁵ Hz. We confirm that an inclusion of the time-dependent changes of the exponent ν makes the newly developed nonstationary three dimensional model of the Forbush decrease compatible with the experimental data. Also, we show that the changes of the rigidity spectrum exponent γ does not depend on the level of convection of the galactic cosmic rays stream by solar wind; depending on the changes of the exponent ν, i.e. on the state of the turbulence of the interplanetary magnetic field.     

Propagation of secondary antiprotons and cosmic rays in the Galaxy

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper, we discuss one possibility to overcome these difficulties. Using the measured antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local “unprocessed” component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z ⩽ 28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three “radioactive clock” isotopes in CR, ¹⁰Be, ²⁶Al, ³⁶Cl, are all consistent and indicate a halo size z_h ∼ 4 kpc based on the most accurate data taken by the ACE spacecraft.     

Measurement of cosmic ray elemental composition from the CAKE balloon experiment

CAKE (Cosmic Abundances below Knee Energies) was a prototype balloon experiment for the determination of the charge spectra and abundances of the primary cosmic rays (CR) with Z > 10. It was a passive instrument made of layers of CR39® and Lexan®/Makrofol® nuclear track detectors; it had a geometric acceptance of ∼0.7 m² sr for Fe nuclei. Here, the scanning and analysis strategies, the algorithms used for the off-line filtering and for the tracking in automated mode of the primary cosmic rays are presented, together with the resulting CR charge distribution and their abundances.     

The impact of coronal mass ejection on the horizontal geomagnetic fields and the induced geoelectric fields

This work investigates the influence of coronal mass ejection (CME) on the time derivatives of horizontal geomagnetic and geoelectric fields, proxy parameters for identifying GICs. 16 events were identified for the year 2003 from the CORONAS-PHOTON spacecraft. Five of the events (May 29, June 9, October 28, October 29, and November 4) were extensively discussed over four magnetic observatories, were analyzed using the time derivatives of the horizontal geomagnetic (dH/dt) and geoelectric (E_H) fields obtained from data of the INTERMAGNET network. It was observed that energy distributions of the wavelet power spectrum of the horizontal geoelectric field are noticed at the nighttime on both 29 May and 9 June 2003 across the stations. Daytime and nighttime intensification of energy distribution of the wavelet power spectrum of the horizontal geoelectric field are observed on both 28 and 29 October 2003 due to strong westward electrojet. The 4 November 2003 event depicts daytime amplification of energy distributions of the wavelet power spectrum across the stations. The highest correlation magnitude is obtained in the event of 4 November 2003 between dH/dt and E_H relationships during the intense solar flare of class X 17.4. We observed that the correlation magnitude between dH/dt and E_H increases with increase in CME activity. We concluded that the response of the surface impedance model for different stations plays a key role in determining the surface electric field strength, due to large electric field changes at different stations.     

Alpha Particle X-Ray Spectrometer (APXS) on-board Chandrayaan-2 rover

Alpha Particle X-ray Spectrometer (APXS) payload configuration for Chandrayaan-2 rover has been completed recently and fabrication of mechanical assembly, PCB layout design and fabrication are in progress. Here we present the design and performance evaluation of various subsystems developed for APXS payload. The low energy threshold of <1 keV and the energy resolution of ∼150 eV at 5.9 keV, for the Silicon Drift Detector (SDD), as measured from the developed APXS electronics is comparable to the standard spectrometers available off-the-shelf. We have also carried out experiments for measuring fluorescent X-ray spectrum from various standard samples from the USGS catalog irradiated by the laboratory X-ray source ²⁴¹Am with 1 mCi activity. It is shown that intensities of various characteristic X-ray lines are well correlated with the respective elemental concentrations.