We investigate acceleration of cosmic rays by shocks and accretion flows in galaxy clusters. Numerical results for spectra of accelerated particles and nonthermal emission are presented. It is shown that the acceleration of protons and nuclei in the nearby galaxy cluster Virgo can explain the observed spectra of ultra high energy cosmic rays.

We consider the diffusive shock acceleration in interstellar bubbles created by powerful stellar winds of supernova progenitors. Under the moderate stellar wind magnetization the bubbles are filled by the strongly magnetized low density gas. It is shown that the maximum energy of particles accelerated in this environment can exceed the "knee" energy in the observable cosmic ray spectrum.

Nonlinear model of diffusive shock acceleration is used for investigation of the particle acceleration in old supernova remnants W28 and W44. We modeled the hydrodynamical evolution of the remnants, shock modification and streaming instability produced by accelerated particles. Comparison with available radio and gamma-ray data is given.

The inverse problem of cosmic ray transport of ultra-high energy cosmic rays is considered. The analysis of Auger data on energy spectrum, energy dependence of mean logarithm of atomic mass number and its variance allows definite conclusions on the shape of the source spectrum in the frameworks of the inverse problem approach. The discussion on regularization procedure for considered ill-posed problem is presented.

The GAMMA-400 gamma-ray telescope with excellent angular and energy resolutions is designed to search for signatures of dark matter in the fluxes of gamma-ray emission and electrons + positrons. Precision investigations of gamma-ray emission from Galactic Center, Crab, Vela, Cygnus, Geminga, and other regions will be performed, as well as diffuse gamma-ray emission, along with measurements of high-energy electron + positron and nuclei fluxes. Furthermore, it will study gamma-ray bursts and gamma-ray emission from the Sun during periods of solar activity. The energy range of GAMMA-400 is expected to be from ~20 MeV up to TeV energies for gamma rays, up to 20 TeV for electrons + positrons, and up to 10E15 eV for cosmic-ray nuclei. For high-energy gamma rays with energy from 10 to 100 GeV, the GAMMA-400 angular resolution improves from 0.1{\deg} to ~0.01{\deg} and energy resolution from 3% to ~1%; the proton rejection factor is ~5x10E5. GAMMA-400 will be installed onboard the Russian space observatory.

It is shown that high-energy astrophysical neutrinos observed in the IceCube experiment can be produced by protons accelerated in extragalactic Type IIn supernova remnants by shocks propagating in the dense circumstellar medium. The nonlinear diffusive shock acceleration model is used for description of particle acceleration. We calculate the neutrino spectrum produced by an individual Type IIn supernova and the spectrum of neutrino background produced by IIn supernovae in the expanding Universe. We also found that the arrival direction of one Icecube neutrino candidate (track event 47) is at 1.35$^{\circ }$ from Type IIn supernova 2005bx.

It is shown that high-energy astrophysical neutrinos observed in the IceCube experiment can be produced by protons accelerated in extragalactic Type IIn supernova remnants by shocks propagating in the dense circumstellar medium. The nonlinear diffusive shock acceleration model is used for description of particle acceleration.

The GAMMA-400 gamma-ray telescope is intended to measure the fluxes of gamma rays and cosmic-ray electrons and positrons in the energy range from 100 MeV to several TeV. Such measurements concern with the following scientific goals: search for signatures of dark matter, investigation of gamma-ray point and extended sources, studies of the energy spectra of Galactic and extragalactic diffuse emission, studies of gamma-ray bursts and gamma-ray emission from the active Sun, as well as high-precision measurements of spectra of high-energy electrons and positrons, protons, and nuclei up to the knee. The main components of cosmic rays are protons and helium nuclei, whereas the part of lepton component in the total flux is ~10E-3 for high energies. In present paper, the capability of the GAMMA-400 gamma-ray telescope to distinguish electrons and positrons from protons in cosmic rays is investigated. The individual contribution to the proton rejection is studied for each detector system of the GAMMA-400 gamma-ray telescope. Using combined information from all detector systems allow us to provide the proton rejection from electrons with a factor of ~4x10E5 for vertical incident particles and ~3x10E5 for particles with initial inclination of 30 degrees. The calculations were performed for the electron energy range from 50 GeV to 1 TeV.

The energy spectra and composition of ultra-high energy cosmic rays are changing in a course of propagation in the expanding Universe filled with background radiation. We developed a numerical code for solution of inverse problem for cosmic-ray transport equations that allows the determination of average source spectra of different nuclei from the cosmic ray spectra observed at the Earth. Employing this approach, the injection spectra of protons and Iron nuclei in extragalactic sources are found assuming that only these species are accelerated at the source. The data from the Auger experiment and the combined data from the Telescope Array + HiRes experiments are used to illustrate the method.

GAMMA-400 is a new space mission, designed as a dual experiment, capable to study both high energy gamma rays (from $\sim$100 MeV to few TeV) and cosmic rays (electrons up to 20 TeV and nuclei up to $\sim$10$^{15}$ eV). The full simulation framework of GAMMA-400 is based on the Geant4 toolkit. The details of the gamma-ray reconstruction pipeline in the pre-shower and calorimeter will be outlined. The performance of GAMMA-400 (PSF, effective area) have been obtained using this framework. The most updated results on them will be shown.

GAMMA-400 is a new space mission which will be installed on board the Russian space platform Navigator. It is scheduled to be launched at the beginning of the next decade. GAMMA-400 is designed to study simultaneously gamma rays (up to 3 TeV) and cosmic rays (electrons and positrons from 1 GeV to 20 TeV, nuclei up to 10$^{15}$-10$^{16}$ eV). Being a dual-purpose mission, GAMMA-400 will be able to address some of the most impelling science topics, such as search for signatures of dark matter, cosmic-rays origin and propagation, and the nature of transients. GAMMA-400 will try to solve the unanswered questions on these topics by high-precision measurements of the Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission and the spectra of cosmic-ray electrons + positrons and nuclei, thanks to excellent energy and angular resolutions.

The present design of the new space observatory GAMMA-400 is presented in this paper. The instrument has been designed for the optimal detection of gamma rays in a broad energy range (from ~100 MeV up to 3 TeV), with excellent angular and energy resolution. The observatory will also allow precise and high statistic studies of the electron component in the cosmic rays up to the multi TeV region, as well as protons and nuclei spectra up to the knee region. The GAMMA-400 observatory will allow to address a broad range of science topics, like search for signatures of dark matter, studies of Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission, gamma-ray bursts and charged cosmic rays acceleration and diffusion mechanism up to the knee.

The measurements of gamma-ray fluxes and cosmic-ray electrons and positrons in the energy range from 100 MeV to several TeV, which will be implemented by the specially designed GAMMA-400 gamma-ray telescope, concern with the following broad range of science topics. Searching for signatures of dark matter, surveying the celestial sphere in order to study gamma-ray point and extended sources, measuring the energy spectra of Galactic and extragalactic diffuse gamma-ray emission, studying gamma-ray bursts and gamma-ray emission from the Sun, as well as high precision measuring spectra of high-energy electrons and positrons, protons and nuclei up to the knee. To clarify these scientific problems with the new experimental data the GAMMA-400 gamma-ray telescope possesses unique physical characteristics comparing with previous and present experiments. For gamma-ray energies more than 100 GeV GAMMA-400 provides the energy resolution of ~1% and angular resolution better than 0.02 deg. The methods developed to reconstruct the direction of incident gamma photon are presented in this paper, as well as, the capability of the GAMMA-400 gamma-ray telescope to distinguish electrons and positrons from protons in cosmic rays is investigated.

Cosmic ray propagation in the Galaxy is shortly reviewed. In particular we consider the self-consistent models of CR propagation. In these models CR streaming instability driven by CR anisotropy results in the Alfv\'enic turbulence which in turn determines the scattering and diffusion of particles.

The processes responsible for the broad-band radiation of the young supernova remnant Cas A are explored using a new code which is designed for a detailed treatment of the diffusive shock acceleration of particles in nonlinear regime. The model is based on spherically symmetric hydrodynamic equations complemented with transport equations for relativistic particles. Electrons, protons and the oxygen ions accelerated by forward and reverse shocks are included in the numerical calculations. We show that the available multi-wavelength observations in the radio, X-ray and gamma-ray bands can be best explained by invoking particle acceleration by both forward and reversed shocks. Although the TeV gamma-ray observations can be interpreted by interactions of both accelerated electrons and protons/ions, the measurements by Fermi LAT at energies below 1 GeV give a tentative preference to the hadronic origin of gamma-rays. Then, the acceleration efficiency in this source, despite the previous claims, should be very high; 25% of the explosion energy (or approximately $3\cdot 10^{50}$ erg) should already be converted to cosmic rays, mainly by the forward shock. At the same time, the model calculations do not provide extension of the maximum energy of accelerated protons beyond 100 TeV. In this model, the acceleration of electrons is dominated by the reverse shock; the required $10^{48}$ erg can be achieved under the assumption that the injection of electrons (positrons) is supported by the radioactive decay of $^{44}$Ti.

Data obtained in the ATIC-2 (Advanced Thin Ionization Calorimeter), CREAM (Cosmic Ray Energetics and Mass)) and PAMELA (Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics) experiments suggest that elemental interstellar spectra of cosmic rays below the knee at a few times $10^{6}$ GeV are not simple power laws, but they experience hardening at magnetic rigidity above about 240 GV. Another essential feature is the difference between proton and Helium energy spectra, so that the He/p ratio increases by more than 50% in the energy range from $10^{2}$ to $10^{4}$ GV. We consider the concavity of particle spectrum resulting from the nonlinear nature of diffusive shock acceleration in supernova remnants (SNR) as a possible reason for the observed spectrum hardening. Helium-to-proton ratio increasing with energy can be interpreted as a consequence of cosmic ray acceleration by forward and reverse shocks in SNRs. The contribution of particles accelerated by reverse shocks makes the concavity of the produced overall cosmic ray spectrum more pronounced. The spectra of protons and helium nuclei accelerated in SNRs and released into the interstellar medium are calculated. The derived steady state interstellar spectra are in reasonably good agreement with observations.

The GAMMA-400 gamma-ray telescope is designed to measure the fluxes of gamma rays and cosmic-ray electrons + positrons, which can be produced by annihilation or decay of the dark matter particles, as well as to survey the celestial sphere in order to study point and extended sources of gamma rays, measure energy spectra of Galactic and extragalactic diffuse gamma-ray emission, gamma-ray bursts, and gamma-ray emission from the Sun. The GAMMA-400 covers the energy range from 100 MeV to 3000 GeV. Its angular resolution is ~0.01 deg (E{\gamma} > 100 GeV), the energy resolution ~1% (E{\gamma} > 10 GeV), and the proton rejection factor ~10E6. GAMMA-400 will be installed on the Russian space platform Navigator. The beginning of observations is planned for 2018.

We argue that the decays of radioactive nuclei related to $^{44}$Ti and $^{56}$Ni ejected during supernova explosions can provide a vast pool of mildly relativistic positrons and electrons which are further accelerated to ultrarelativistic energies by reverse and forward shocks. This interesting link between two independent processes - the radioactivity and the particle acceleration - can be a clue for solution of the well known theoretical problem of electron injection in supernova remnants. In the case of the brightest radio source Cas A, we demonstrate that the radioactivity can supply adequate number of energetic electrons and positrons for interpretation of observational data provided that they are stochastically pre-accelerated in the upstream regions of the forward and reverse shocks.

A new numerical model of the nonlinear diffusive shock acceleration is presented. It is used for modeling of particle acceleration in supernova remnants. The model contains coupled spherically symmetric hydrodynamic equations and the transport equations for energetic protons, ions and electrons. The forward and reverse shocks are included in the consideration. The spectra of cosmic rays released into interstellar medium from a supernova remnant are determined. The role of the reverse shock in the production of CR ions and positrons is discussed.

The cosmic ray source spectrum produced by AGN (Active Galactic Nucleus) jets is calculated. A distinctive feature of these calculations is the account for the jet distribution on kinetic energy. The expected cosmic ray spectrum at the Earth is determined with the use of a simple numerical code which takes into account interactions of ultra-high energy protons and nuclei with the background radiation in an expanding universe.

The spectra of high-energy protons and nuclei accelerated by supernova remnant shocks are calculated taking into account magnetic field amplification and Alfvenic drift both upstream and downstream of the shock for different types of supernova remnants during their evolution. The maximum energy of accelerated particles may reach $5\cdot10^{18}$ eV for Fe ions in Type IIb SNRs. The calculated energy spectrum of cosmic rays after propagation through the Galaxy is in good agreement with the spectrum measured at the Earth.

The spectra of high-energy protons and nuclei accelerated by supernova remnant shocks are calculated taking into account magnetic field amplification and Alfvenic drift for different types of SNRs during their evolution. The overall energy spectrum and elemental composition of cosmic rays after propagation through the Galaxy are found.

A new numerical code, designed for the detailed numerical treatment of nonlinear diffusive shock acceleration, is used for modeling of particle acceleration and radiation in young supernova remnants. The model is based on spherically symmetric hydrodynamic equations complemented with transport equations for relativistic particles. For the first time, the acceleration of electrons and protons by both forward and reverse shocks is studied through detailed numerical calculations. We model the energy spectra and spatial distributions of nonthermal emission of the young supernova remnant RX J1713.7-3946 and compare the calculations with the spectral and morphological properties of this object obtained in broad energy band from radio to very high energy gamma-rays. We discuss the advantages and shortcomings of the so-called hadronic and leptonic models which assume that the observed TeV gamma-ray emission is produced by accelerated protons and electrons, respectively. We discuss also a "composite" scenario when the gamma-ray flux from the main parts of the shell has inverse Compton origin, but with a non-negligible contribution of hadronic origin from dense clouds interacting with the shell.

Cosmic ray acceleration in SNRs in the presence of the Alfv\'enic drift is considered. It is shown that spectra of accelerated particles may be considerably softer in the presence of amplified magnetic fields.

We have studied the non-resonant streaming instability of charged energetic particles moving through a background plasma, discovered by Bell (2004). We confirm his numerical results regarding a significant magnetic field amplification in the system. A detailed physical picture of the instability development and of the magnetic field evolution is given.

We investigate the diffusive shock acceleration in the presence of the non-resonant streaming instability introduced by Bell (2004). The numerical MHD simulations of the magnetic field amplification combined with the analytical treatment of cosmic ray acceleration permit us to calculate the maximum energy of particles accelerated by high-velocity supernova shocks. The estimates for Cas A, Kepler, SN1006, and Tycho historical supernova remnants are given. We also found that the amplified magnetic field is preferentially oriented perpendicular to the shock front downstream of the fast shock. This explains the origin of the radial magnetic fields observed in young supernova remnants.

%context {Recent observations of hard X-rays and very high energy gamma-rays from a number of young shell type supernova remnants indicate the importance of detailed quantitative studies of energy spectra of relativistic electrons formed via diffusive shock acceleration accompanied by intense nonthermal emission through synchrotron radiation and inverse Compton scattering.} %aim {The aim of this work was derivation of exact asymptotic solutions of the kinetic equation which describes the energy distribution of shock-accelerated electrons for an arbitrary energy-dependence of the diffusion coefficient.} %method {The asymptotic solutions at low and very high energy domains coupled with numerical calculations in the intermediate energy range allow analytical presentations of energy spectra of electrons for the entire energy region.} %results {Under the assumption that the energy losses of electrons are dominated by synchrotron cooling, we derived the exact asymptotic spectra of electrons without any restriction on the diffusion coefficient. We also obtained simple analytical approximations which describe, with accuracy better than ten percent, the energy spectra of nonthermal emission of shock-accelerated electrons due to the synchrotron radiation and inverse Compton scattering.} %conclusions {The results can be applied for interpretation of X-ray and gamma-ray observations of shell type supernova remnants, as well as other nonthermal high energy source populations like microquasars and large scale synchrotron jets of active galactic nuclei.

Using a so-called hemispherical model we derive a general transport equation for cosmic ray and thermal particles scattered in pitch angle by magnetic inhomogeneities in a moving collisionless plasma. The weak scattering through 90 degrees results in isotropic particle distributions in each hemisphere. The consideration is not limited by small anisotropies and by the condition that particle velocities are higher than characteristic flow velocity differences. For high velocities and small anisotropies the standard cosmic ray transport equation is recovered. We apply the equations derived for investigation of injection and acceleration of particles by collisionless shocks.

The physical processes involved in diffusion of Galactic cosmic rays in the interstellar medium are addressed. We study the possibility that the nonlinear MHD cascade sets the power-law spectrum of turbulence which scatters charged energetic particles. We find that the dissipation of waves due to the resonant interaction with cosmic ray particles may terminate the Kraichnan-type cascade below wavelengths 10^13 cm. The effect of this wave dissipation has been incorporated in the GALPROP numerical propagation code in order to asses the impact on measurable astrophysical data. The energy-dependence of the cosmic-ray diffusion coefficient found in the resulting self-consistent model may explain the peaks in the secondary to primary nuclei ratios observed at about 1 GeV/nucleon.

We present an analytical model for the hydrodynamic outflow from the disk of a starburst galaxy. The model is used to calculate the cosmic ray propagation and the radio intensity distribution in the nuclear starbust region of NGC 253. We find that the cosmic ray energy production rate of the central 600 pc of NGC 253 is about $3\cdot 10^{41}$ erg s$^{-1}$, that is about 15 percent of the total mechanical supernova power. For this inner region we estimate a terminal outflow velocity of 900 km s$^{-1}$ and a mass loss rate of $\dot{M}=2\div 4M_\odot$ per year.

The cosmic-ray streaming instability creates strong magnetohydrodynamic turbulence in the precursor of a SN shock. The level of turbulence determines the maximum energy of cosmic-ray particles accelerated by the diffusive shock acceleration mechanism. The present consideration continues our work Ptuskin & Zirakashvili (2003). 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 over the course of adiabatic SNR evolution (the Sedov stage) is approximately $Q(p)p^{2}\propto p^{-2}$ at energies larger than 10-30 Gev/nucleon and with the maximum particle energy that is close to the knee position in cosmic ray spectrum observed at $\sim4\times10^{15}$ eV. At earlier stage of SNR evolution - the ejecta dominated stage described by the Chevalier-Nadyozhin solution, the particles are accelerated to higher energies and have rather steep power-law distribution on momentum. 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.

Cosmic ray acceleration by shocks related with Slipping Interaction Regions (SIRs) in the Galactic Wind is considered. SIRs are similar to Solar Wind Corotating Interaction Regions. The spiral structure of our Galaxy results in a strong nonuniformity of the Galactic Wind flow and in SIR formation at distances of 50 to 100 kpc. SIRs are not corotating with the gas and magnetic field because the angular velocity of the spiral pattern differs from that of the Galactic rotation. It is shown that the collective reacceleration of the cosmic ray particles with charge $Ze$ in the resulting shock ensemble can explain the observable cosmic ray spectrum beyond the "knee" up to energies of the order of $10^{17}Z$ eV. For the reaccelerated particles the Galactic Wind termination shock acts as a reflecting boundary.

We study the production of cosmic rays (CRs) in supernova remnants (SNRs), including the reacceleration of background galactic cosmic rays (GCRs) - thus refining the early considerations by Blandford & Ostriker (1980) and Wandel et al. (1987) - and the effects of the nuclear spallation inside the sources (the SNRs). This combines for the first time nuclear spallation inside CR sources and in the diffuse interstellar medium, as well as reacceleration, with the injection and subsequent acceleration of suprathermal particles from the postshock thermal pool. Selfconsistent CR spectra are calculated on the basis of the nonlinear kinetic model. It is shown that GCR reacceleration and CR spallation produce a measurable effect at high energies, especially in the secondary to primary (s/p) ratio, making its energy-dependence substantially flatter than predicted by the standard model. Quantitatively, the effect depends strongly upon the density of the surrounding circumstellar matter. GCR reacceleration dominates secondary CR production for a low circumstellar density. It increases the expected s/p ratio substantially and flattens its spectrum to an almost energy-independent form for energies larger than 100 GeV/n if the supernovae explode on average into a hot dilute medium with hydrogen number density $N_H=0.003$cm$^{-3}$. The contribution of CR spallation inside SNRs to the s/p ratio increases with increasing circumstellar density and becomes dominant for $N_H\gsim 1$ cm$^{-3}$, leading at high energies to a flat s/p ratio which is only by a factor of three lower than in the case of the hot medium. Measurements of the boron to carbon ratio at energies above 100 GeV/n could be used in comparison with the values predicted here as a consistency test for the supernova origin of the GCRs.

The diffusion of galactic cosmic rays (CR) is considered. It is anticipated that the nonlinear MHD cascade sets the power-law spectrum of turbulence from the principle scale 100 pc to much smaller scales. We found that the dissipation of waves due to the resonant interaction with energetic particles may terminate the cascade at less than 1.e13 cm. The shape of CR diffusion coefficient that was found may explain the observed peaks of secondary to primary nuclei ratios at a few GeV/n.

The instability in the cosmic-ray (CR) precursor of a SN shock is studied. The level of turbulence in this region determines the maximum energy of accelerated CRs. The consideration is not limited by the case of weak turbulence. It is assumed that 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. It can be as high as $10^{17}Z$ eV in young SNR and falls down to about $10^{10}Z$ eV at the end of Sedov stage ($Z$ is the particle charge). This finding may explain why the SNRs with the age more than a few thousand years are not prominent sources of very high energy $\gamma $-rays. The averaged spectrum of ultrarelativistic CR injected in the interstellar medium is close or somewhat steeper than $E^{-2}$.

The instability in the cosmic-ray precursor of a supernova shock moving in interstellar medium is studied. The level of magnetohydrodynamic turbulence in this region determines the maximum energy of particles accelerated by the diffusive shock acceleration mechanism. The high efficiency of cosmic ray acceleration is accepted, and the consideration is not limited by the case of weak turbulence. It is assumed that 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 supernova remnant. The maximum energy can be as high as $\sim10^{17}Z$ eV in young supernova remnants and falls down to about $\sim10^{10}Z$ eV at the end of Sedov stage. Thus the standard estimate of maximum particle energy based on the Bohm limit calculated for the interstellar magnetic field strength is not justified in this case. This finding may explain why the supernova remnants with the age more than a few thousand years are not prominent sources of very high energy gamma-rays.

Cosmic ray acceleration in turbulent interstellar medium is considered. Turbulence is treated as ensemble of moving magnetic traps. We derive equations for particle momentum distribution function that describes acceleration of particles in this case. Rate of acceleration calculated is estimated for our Galaxy and compared with ones given by other acceleration mechanisms.

Nonlinear damping of parallel propagating Alfv\'en waves in high-$\beta$ plasma is considered. Trapping of thermal ions and Coulomb collisions are taken into account. Saturated damping rate is calculated. Applications are made for cosmic ray propagation in the Galaxy.

The equation for temporary evolution of spectral energy of collisionless Alfven waves is derived in framework of weak turbulence theory. The main nonlinear processes for such conditions are induced scattering and two quantum absorption of Alfven waves by thermal ions. The equation for velocity distribution of thermal particles is derived that describes diffusion in momentum space due to this nonlinear processes. Comparison is done with the results of another authors. Results obtained are qualitatively differ from the ones obtained for the case of Alfven waves propagation along mean magnetic field.