Cosmic-ray composition and its relation to shock acceleration by supernova remnants

An overview is given on the present status of the understanding of the origin of galactic cosmic rays. Recent measurements of charged cosmic rays and photons are reviewed. Their impact on the contemporary knowledge about the sources and acceleration mechanisms of cosmic rays and their propagation through the Galaxy is discussed. Possible reasons for the knee in the energy spectrum and scenarios for the end of the galactic cosmic-ray component are described.     

The variability of the proton cosmic ray flux on the Sun’s way around the galactic center

The discovery of direct evidences for the acceleration of high energetic particles at the shell supernova remnant RXJ1713.7-3946 underlined the need to calculate the cosmic ray (CR) distribution in the Galaxy on a spatial grid fine enough to resolve the changes in the CR density due to these kind of objects. It was shown before by Büsching et al. [Büsching, I., Kopp, A., Pohl, M., Schlickeiser, R., Perrot, C., Grenier, I. Cosmic-ray propagation properties for an origin in supernova remnants. ApJ 619, 314–326, 2005] that the discrete nature (both in space and time) of super novae (SN) as sources of Galactic CR leads to CR spectra changing in space and time, resulting in a range of possible CR spectra at a given location in the Galaxy. As the most frequent SN types Ib and II are found within spiral arms, one can expect a significant difference of the range of possible spectra in and outside spiral arms. We investigate the variability of the local interstellar CR proton spectrum during the motion of the Sun in and out of spiral arms in its journey around the Galactic center. Using the code described by Büsching et al. [Büsching, I., Kopp, A., Pohl, M., Schlickeiser, R., Perrot, C., Grenier, I. Cosmic-ray propagation properties for an origin in supernova remnants. ApJ 619, 314–326, 2005], the proton CR density in the Galaxy is calculated with high spatial and temporal resolution (75 pc in galactocentric radius r and azimuth φ at the position of the Sun and 20 pc in z (perpendicular to the galactic plane), with a time step of 1 kyr), assuming stochastically distributed point sources with a probability distribution that resembles the spiral structure of our Galaxy. We find the averaged CR flux outside the spiral arm to be about 50% of that inside a spiral arm. We further find spatial and temporal variations of the CR flux inside spiral arms.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

A monople model for annihilation line emission from the galactic center

A model of the galactic center containing magnetic monopoles has been suggested. We have estimated the number of monopoles preserved in the formation process of the galactic nuclei. The catalytic reaction between magnetic monopoles and baryons can produce positrons and release the energy as well. According to our model, the galactic center cannot collapse into a black, even in the case where the contents of magnetic monopoles are much less than the Parker limit. This model may account for the intensity of the positron-electron pair annihilation line and the high energy photons continuum (>500 keV) observed at the galactic center direction.     

A model of galactic cosmic rays for use in calculating linear energy transfer spectra.

The galactic cosmic rays (GCR) contain fully stripped nuclei, from Hydrogen to beyond the Iron group, accelerated to high energies and are a major component of the background radiation encountered by satellites and interplanetary spacecraft. This paper presents a GCR model which is based upon our current understanding of the astrophysics of GCR transport through interstellar and interplanetary space. The model can be used to predict the energy spectra for all stable and long-lived radioactive species from H to Ni over an energy range from 50 to 50,000 MeV/nucleon as a function of a single parameter, the solar modulation level phi. The details of this model are summarized, phi is derived for the period 1974 to present, and results from this model during the 1990/1991 CRRES mission are presented.     

Modulation of cosmic rays: Perpendicular diffusion and drifts in a heliosphere with a solar wind termination shock

In this study the roles of polar perpendicular diffusion and drifts are illustrated in a model containing a heliosheath and diffusive shock acceleration as applied to the solar wind termination shock. Of particular interest is the relation of polar perpendicular diffusion to particle drifts and how the effectiveness of the termination shock acceleration of galactic and anomalous protons is influenced by this relation. We found that drifts have a more prominent effect than the polar enhancement of perpendicular diffusion so that its omission from termination shock models would produce unrealistically large shock acceleration and consequently also larger modulation effects throughout the heliosphere. The computed spectra at a heliolatitude of 35° are almost similar for the two polarity epochs indicating that the two Voyager spacecraft might not observe differences between the two cycles in future.     

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

Estimation of long-term cosmic ray intensity variation in near future and prediction of their contribution in expected global climate change

On the basis of results obtained in our paper [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], we determine: the dimension of the Heliosphere (modulation region), radial diffusion coefficient and other parameters of convection–diffusion; drift mechanisms of long-term variations of cosmic ray (CR) dependence on particle energy; level of solar activity (SA); and generally, the solar magnetic field. We obtain this important information on the basis of CR and SA data in the past, taking into account the theory of convection–diffusion and global drift modulation of galactic CR in the Heliosphere. By using these results and other regularly published predictions of expected SA variation in the near future, as well as predictions of the next SA cycle, we may make predictions of long-term cosmic ray intensity variation expected in the near future (up to 10–12 years). In [Dorman, L.I. Long-term cosmic ray intensity variation and part of global climate change, controlled by solar activity through cosmic rays, Paper D2.1/C2.2/E3.1-0097-04. Adv. Space Res., 2004 (accepted)], properties of connections between long-term variation in CR intensity and some part of a global climate change were estimated, controlled by solar activity through CR. We show that in this way we may make predictions of some part of a global climate change expected in the near future (up to 10–12 years and maybe more, depending upon the period during which definite predictions of SA can be made), controlled by solar activity through CR. In this case, estimations of expected long-term changes in the planetary distribution of cutoff rigidities, which also influence CR intensity, as well as CR-influenced effects on global climate variation, become important.     

X-ray investigations of the hot ISM

At energies less than one keV, the intensity of the galactic x-ray background dominates that of the extragalactic background in almost every direction on the sky. Below 1/4 keV, the galactic x-ray background has a galactic stellar component, but the dominant emitter seems to be hot interstellar matter. The origin of the general 3/4 keV x-ray background remains uncertain, but one component must also be the contribution from hot interstellar matter. An overview is given of recent x-ray investigations of the hot interstellar medium using data from the ROSAT XRT/PSPC instrument. Several prominent features in the low energy x-ray background that are interpreted as fossil supernova remnants are discussed.     

The composition of cosmic rays at high energies

Measurements of the composition of the cosmic rays at high energies, and of the energy spectra of the individual components provide the basis for the understanding of the sources, of the acceleration mechanism, and of the galactic containment of these particles. We briefly review the presently available information, and we describe a recent measurement on the Space Shuttle that we performed in order to substantially extend the range of energies in which the elemental composition is known. We present and discuss the results, and we also summarize and discuss recent data on the electron component of cosmic rays. The body of data now available contains several features that are difficult to explain within current models of galactic shock acceleration and “leaky box” containment. We emphasize the need for further measurements, and we briefly discuss possible opportunities for future work.     

High resolution X-ray spectroscopy and chemical evolution in elliptical galaxies

We discuss the elemental composition in the hot ISM of elliptical galaxies derived from new and recent X-ray spectral analysis in the context of new phenomenological models of their chemical evolution. Star formation histories, the IMF, the astrophysics of supernovae, and the nature of galactic winds impact the metal content and relative abundances in the hot ISM. We evaluate how X-ray spectroscopy may be utilized to deconstruct how elliptical galaxies, and the stellar populations that compose them, form and evolve, with an emphasis on present and future high resolution spectroscopic analysis.     

太阳宇宙线源物质的高色球层模式和重离子丰度过量机制的探讨

本文利用太阳能量粒子事件中重离子平均丰度过量的资料,计算得到太阳能量粒子源物质的温度,提出了描述太阳宇宙线能量粒子源物质的新模式——高色球层模式;太阳耀斑观测确定,太阳宇宙线耀斑的加速区一般最可能出现在低日冕甚至高达几万公里的高度,从而,太阳宇宙线的源和加速区通常不位于同一区域;进而提出了描述太阳能量粒子事件中重离子丰度过量的可能机制——其源物质是通过太阳黑子的冻结型无力场从高色球层输送到活动区,形成耀斑前加速区内重离子丰度大和耀斑后宇宙线中重元素丰度的过量.      

Computation of radiation environment during ground level enhancements 65, 69 and 70 at equatorial region and flight altitudes

The radiation environment in the troposphere of the Earth is governed by cosmic rays of galactic and solar origin. During major solar energetic particles events the radiation environment changes dramatically. As a results the risk of biological effects due to exposure to ionizing radiation of aircrew increases. Here we present a numerical model for computation of absorbed dose in air due to cosmic rays of galactic and solar origin. It is applied for computation of radiation environment at flight altitude in the equatorial region during several major ground level enhancements, namely GLE65 on 28 October 2003, GLE69 on 20 January 2005 and GLE70 on 13 December 2006. The model is based on a full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The cascade simulation is carried out with CORSIKA 6.990 code with corresponding hadron generators FLUKA 2011 and QGSJET II. The contribution of different cascade components, namely electromagnetic, hadron and muon is explicitly obtained. The spectra of arriving solar energetic particles are calculated from ground level measurements with neutron monitors and satellite data from GOES. The obtained results are discussed.     

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.     

The onset of sunspot cycle 24 and galactic cosmic ray modulation

The annual mean sunspot number (SSN) has a minimum value in 2008, while the monthly mean SSN has a value of zero in August 2009. The galactic cosmic ray modulation for cycle 24 began at earth orbit in January 2010. We study the onset characteristics of the new modulation cycle using data from the global network of neutron monitors. They respond to time variations in different segments of the galactic cosmic ray rigidity spectrum. The corresponding temporal variations in the interplanetary magnetic field intensity (B) and solar wind velocity (V) as well as the tilt angle of the heliospheric current sheet are also studied. There is a lag of 3 months between a large, sharp increase of the tilt angle of the heliospheric current sheet and the onset of modulation. Some neutron monitors are undergoing long-term drifts of unknown origin.     

Improving the angular resolution of EGRET and new limits on supersymmetric dark matter near the galactic center

Using the EGRET data and an improved point source analysis, including an energy-dependent point spread function and an unbinned maximum likelihood technique, we have been able to place considerably lower limits on the γ-ray flux from the galactic center region. We also test this method on known sources, the Crab and Vela pulsars. In both cases, we find that our method improves the angular precision of EGRET data over the 3EG catalog.This new limit on γ-rays from the galactic center can be used to test models of annihilating supersymmetric dark matter and galactic halo profiles. We find that the present EGRET data can limit many supersymmetric models if the density of the galactic dark matter halo is cuspy or spiked toward the galactic center. We also discuss the ability of GLAST to test these models.     

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.     

Constraints on a hadronic model for unidentified off-plane galactic gamma-ray sources

Recently the H.E.S.S. collaboration announced the detection of an unidentified gamma-ray source with an off-set from the galactic plane of 3.5°: HESS J1507-622. If the distance of the object is larger than about one kpc it would be physically located outside the galactic disk. The density profile of the ISM perpendicular to the galactic plane, which acts as target material for hadronic gamma-ray production, drops quite fast with increasing distance. This fact places distance dependent constraints on the energetics and properties of off-plane gamma-ray sources like HESS J1507-622 if a hadronic origin of the gamma-ray emission is assumed. For the case of this source it is found that there seems to be no simple way to link this object to the remnant of a stellar explosions.     

Nature of gamma-ray burst sources

Observational evidence suggests that gamma ray bursts have a local galactic origin involving neutron stars. In this light we make a critical review of physics of the thermonuclear runaway model placing emphasis on self-consistency. We further show that some of the proposed models can be observationally excluded in the light of existing data from the Einstein Observatory. The possibility of gamma bursts arising in low mass binaries is finally discussed in the light of evolutionary scenarios leading to low luminosity systems.     

Analytical approach to cosmic ray ionization by nuclei with charge Z in the middle atmosphere – Distribution of galactic CR effects

The effects of galactic and solar cosmic rays (CR) in the middle atmosphere are considered in this work. A new analytical approach for CR ionization by protons and nuclei with charge Z in the lower ionosphere and middle atmosphere is developed in this paper. For this purpose the ionization losses (dE/dh) according to the Bohr–Bethe–Bloch formula for the energetic charged particles are approximated in three different energy intervals. More accurate expressions for energy decrease E(h) and electron production rate profiles q(h) are derived. The obtained formulas allow comparatively easy computer programming. The integrand in q(h) gives the possibility for application of adequate numerical methods – such as Romberg method or Gauss quadrature, for the solution of the mathematical problem. On this way the process of interaction of cosmic ray particles with the upper, middle and lower atmosphere will be described much more realistically. Computations for cosmic ray ionization in the middle atmosphere are made. The full CR composition is taken into account: protons, Helium (α-particles), light L, medium M, heavy H and very heavy VH group of nuclei.