We aim to elucidate cosmic ray (CR) propagation in the weakly ionized environments of supernova remnants (SNRs) basing our analysis on the cosmic ray cloud (CRC) model developed by \citet{2013ApJ...768...73M} and \citet{2016MNRAS.461.3552N}. We solved two transport equations simultaneously: one for the CR pressure and one for the Alfv\'en wave energy density where CRs are initially confined in the SNR shock. Cosmic rays trigger a streaming instability and produce slab-type resonant Alfv\'en modes. The self-generated turbulence is damped by ion-neutral collisions and by noncorrelated interaction with Alfv\'en modes generated at large scales. We show that CRs leaking in cold dense phases such as those found in cold neutral medium (CNM) and diffuse molecular medium (DiM) can still be confined over distances of a few tens of parsecs from the CRC center for a few thousand years. At 10 TeV, CR diffusion can be suppressed by two or three orders of magnitude. This effect results from a reduced ion-neutral collision damping in the decoupled regime. We calculate the grammage of CRs in these environments. We find that in both single and multi-phase setups at 10 GeV, CNM and DiM media can produce grammage in the range 10-20 $\rm{g/cm^2}$ in the CNM and DiM phases. At 10 TeV, because of nonlinear propagation the grammage increases to values in the range 0.5-20 $\rm{g/cm^2}$ in these two phases. We also present preliminary calculations in inhomogeneous interstellar medium combining two or three different phases where we obtain the same trends.

The paper is a script of a lecture given at the ISAPP-Baikal summer school in 2018. The lecture gives an overview of the Tunka Advanced Instrument for cosmic rays and Gamma Astronomy (TAIGA) facility including historical introduction, description of existing and future setups, and outreach and open data activities.

We study the problem of the escape and transport of Cosmic-Rays (CR) from a source embedded in a fully ionised, hot phase of the interstellar medium (HIM). In particular, we model the CR escape and their propagation in the source vicinity taking into account excitation of Alfv\'enic turbulence by CR streaming and mechanisms damping the self-excited turbulence itself. Our estimates of escape radii and times result in large values (100 pc, $2\times10^5$ yr) for particle energies $\lesssim20$ GeV and smaller values for particles with increasing energies (35 pc and 14 kyr at 1 TeV). These escape times and radii, when used as initial conditions for the CR propagation outside the source, result in relevant suppression of the diffusion coefficient (by a factor 5-10) on time-scales comparable with their (energy dependent) escape time-scale. The damping mechanisms are fast enough that even on shorter time scales the Alfv\'enic turbulence is efficiently damped, and the ratio between random and ordered component of the magnetic field is $\delta B/B_0\ll 1$, justifying the use of quasi-linear theory. In spite of the suppressed diffusion coefficient, and then the increased residence time in the vicinity (<200 pc) of their source, the grammage accumulated by CRs after their escape is found to be negligible (at all energies) as compared to the one accumulated while diffusing in the whole Galaxy, due to the low density of the HIM.

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.

Supernova remnants are believed to be the main sources of galactic Cosmic Rays (CR). Within this framework, particles are accelerated at supernova remnant shocks and then released in the interstellar medium. The mechanism through which CRs are released and the way in which they propagate still remain open issues. The main difficulty is the high non-linearity of the problem: CRs themselves excite the magnetic turbulence that confines them close to their sources. We solve numerically the coupled differential equations describing the evolution in space and time of the escaping particles and of the waves generated through the CR streaming instability. The warm ionized and warm neutral phases of the interstellar medium are considered. These phases occupy the largest fraction of the disc volume, where most supernovae explode, and are characterised by the significant presence of neutral particles. The friction between those neutrals and ions results in a very effective wave damping mechanism. It is found that streaming instability affects the propagation of CRs even in the presence of ion-neutral friction. The diffusion coefficient can be suppressed by more than a factor of $\sim 2$ over a region of few tens of pc around the remnant. The suppression increases for smaller distances. The propagation of $\approx 10$ GeV particles is affected for several tens of kiloyears after escape, while $\approx 1$ TeV particles are affected for few kiloyears. This might have a great impact on the interpretation of gamma-ray observations of molecular clouds located in the vicinity of supernova remnants.

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.

Recent advances on gamma-ray observations from SuperNova Remnants and Molecular Clouds offer the possibility to study in detail the properties of the propagation of escaping Cosmic Rays (CR). However, a complete theory for CR transport outside the acceleration site has not been developed yet. Two physical processes are thought to be relevant to regulate the transport: the growth of waves caused by streaming instability, and possible wave damping mechanisms that reduce the growth of the turbulence. Only a few attempts have been made so far to incorporate these mechanisms in the theory of CR diffusion. In this work we present recent advances in this subject. In particular, we show results obtained by solving the coupled equations for the diffusion of CRs and the evolution of Alfven waves. We discuss the importance of streaming instabilities and wave damping in different ISM phases.

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

The role of nearby galactic sources, the supernova remnants, in formation of observed energy spectrum and anisotropy of high-energy cosmic rays is studied. The list of these sources is made up based on radio, X-ray and gamma-ray catalogues. The distant sources are treated statistically as ensemble of sources with random positions and ages. The source spectra are defined based on the modern theory of cosmic ray acceleration in supernova remnants while the propagation of cosmic rays in the interstellar medium is described in the frameworks of galactic diffusion model. Calculations of anisotropy are made to reproduce the experimental procedure of getting the two-dimensional anisotropy. The detailed comparison of calculations with cosmic ray data is made. We explained simultaneously the new cosmic ray data on a fine structure of all particle spectra around the knee and the amplitude and direction of the dipole component of anisotropy in the wide energy range 1 TeV - 1 EeV. Made assumptions do not look exotic, and they confirm the modern understanding of cosmic ray origin.

In this paper, a probable interpretation of a remarkable fine structure of all particle spectra between the knee and the ankle, as well as a high content of heavy nuclei around 10^17 eV measured recently in Tunka-133 and KASCADE Grande experiments, is presented in the model where production of cosmic rays in Type Ia SNRs provides observed cosmic rays (CR) flux around the knee. Subtracting the contribution of these sources from all particle spectrum we obtained the residual flux of CR in the transition region between Galactic and Extragalactic CR in the range 10^17-5 10^18 eV. The obtained spectrum also has a pronounced knee at energy (2-3)10^17 eV and slopes g1 ~ 1.8+-0.3 below it and g2 ~3.4 above it. We analyzed the possible contribution from known Galactic sources and found that formally the best candidate to contribute significantly to the transition region is Cas A. This source posses a number of unusual properties and considered as Type IIb SNR. However we showed that the hypothesis of extragalactic origin of CR starting from 10^17 eV seems more realistic.

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.

A new EAS Cherenkov light array, Tunka-133, with ~1 km^2 geometrical area has been installed at the Tunka Valley (50 km from Lake Baikal) in 2009. The array permits a detailed study of cosmic ray energy spectrum and mass composition in the energy range 10^16 - 10^18 eV with a uniform method. We describe the array construction, DAQ and methods of the array calibration.The method of energy reconstruction and absolute calibration of measurements are discussed. The analysis of spatial and time structure of EAS Cherenkov light allows to estimate the depth of the EAS maximum X_max. The results on the all particles energy spectrum and the mean depth of the EAS maximum X_max vs. primary energy derived from the data of two winter seasons (2009 -- 2011), are presented. Preliminary results of joint operation of the Cherenkov array with antennas for detection of EAS radio signals are shown. Plans for future upgrades -- deployment of remote clusters, radioantennas and a scintillator detector network and a prototype of the HiSCORE gamma-telescope -- are discussed.

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.

We briefly review sources of cosmic rays, their composition and spectra as well as their propagation in the galactic and extragalactic magnetic fields, both regular and fluctuating. A special attention is paid to the recent results of the X-ray and gamma-ray observations that shed light on the origin of the galactic cosmic rays and the challenging results of Pierre Auger Observatory on the ultra high energy cosmic rays. The perspectives of both high energy astrophysics and cosmic-ray astronomy to identify the sources of ultra high energy cosmic rays, the mechanisms of particle acceleration, to measure the intergalactic radiation fields and to reveal the structure of magnetic fields of very different scales are outlined.

The simple leaky-box model of propagation of cosmic rays in the Galaxy is quite suitable for handling data on cosmic ray nuclei energy spectra and composition at E ? 1 GeV [1,2]. In the leakybox model a full information about cosmic ray propagation in Galaxy is compressed to the single parameter - escape length, Xe, characterizing mean grams of a matter passed by cosmic rays from sources to the Earth. In this paper we analyze the world data on proton and iron cosmic ray spectra collected in the past (HEAO, CRN et al.) and in series of recent electronic experiments (ATIC, CREAM, MS, BESS, Tracer et al.) and obtain the rigidity dependence of escape length, Xe(R) = R^(-0.47\pm-0.03), from the measured rigidity dependence of the protons/iron ratio. It quite agrees with the one stimated in standard manner from the secondary/primary nuclei ratio. But at R > 300 GV the behavior of Xe(R)distinctly changes, that can (variant of explanation)point out to the change of proton/iron ratio in cosmic ray sources.

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 new EAS Cherenkov array Tunka-133 with about 1 km**2 geometric acceptance area is installed in the Tunka Valley (50 km from Lake Baikal). The array willpermit a detailed study of cosmic ray energy spectrum and mass composition in the energy range of 10**15 - 10**18 eV with uniform method. The array consistsof 19 clusters, each composed of 7 optical detectors with 20 cm PMTs. Since November 2008, the first 12 clusters are in operation, commissioning of the whole array is planned for September 2009 (At the time of submission of this paperto electronic arXiv(February 2010) the completed Tunka-133 array is already taking data). We describe the array construction and DAQ, preliminary results and plans for the future development: deployment of radio-antennas and muon detectors network.

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.

The final analysis of the Extensive Air Shower (EAS) maximum X_max depth distribution derived from the data of Tunka-25 atmospheric Cherenkov light array in the energy range 3.10^15 - 3.10^16 eV is presented. The perspectives of X_max studies with the new Cherenkov light array Tunka-133 of 1 km^2 area, extending the measurements up to 10^18 eV, are discussed.

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.

The new EAS Cherenkov array TUNKA-133, with about 1 km**2 sensitive area, is being installed in the Tunka Valley. The investigated energy range is 10**15-10**18 eV. It will consist of 133 optical detectors based on EMI9350 PMTs. Optical detectors are grouped into 19 clusters with 7 detectors each. The detectors are connected to the cluster box with RG-58 cables. Every PMT signal is digitized in the cluster box with 200 MHz FADC. The cluster boxes are connected to the data acquisition center with a 1 Gb/s optical link. A detailed description of the data acquisition system (DAQ) is presented.

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.

The new EAS Cherenkov array Tunka-133, with about 1 km**2 sensitive area, is being installed in the Tunka Valley since the end of 2005. This array will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range of 10**15-10**18 eV with a unique method. The array will consist of 19 clusters, each composed of 7 optical detectors. The first cluster started operation in October 2006. We describe the data acquisition system and present preliminary results from data taken with the first cluster.

We survey the theory and experimental tests for the propagation of cosmic rays in the Galaxy up to energies of 10^15 eV. A guide to the previous reviews and essential literature is given, followed by an exposition of basic principles. The basic ideas of cosmic-ray propagation are described, and the physical origin of its processes are explained. The various techniques for computing the observational consequences of the theory are described and contrasted. These include analytical and numerical techniques. We present the comparison of models with data including direct and indirect - especially gamma-ray - observations, and indicate what we can learn about cosmic-ray propagation. Some particular important topics including electrons and antiparticles are chosen for discussion.

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.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

The project of an EAS Cherenkov array in the Tunka valley/Siberia with an area of about 1 km^2 is presented. The new array will have a ten times bigger area than the existing Tunka-25 array and will permit a detailed study of the cosmic ray energy spectrum and the mass composition in the energy range from 10^15 to 10^18 eV.

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.

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.

In a recent paper (Moskalenko et al., 2002), it has been shown that the flux of secondary cosmic ray (CR) antiprotons appears to be contradictory to measurements of secondary to primary nuclei ratios in cosmic rays when calculated in the same Galactic propagation model. The contradiction appears as a value of the diffusion coefficient necessary to match the secondary ratios antiproton/proton and B/C. In particlular, it was shown that the reacceleration models designed to match secondary to primary nuclei ratios produce too few antiprotons. It is, however, clear that some reacceleration is unavoidable in the turbulent interstellar medium. Here we discuss an idea of how to improve reacceleration model by allowing for the damping of interstellar turbulences on the small scale by cosmic rays, mostly protons. This would lead to increase in the mean free path lenghts at low energies, the well-known phenomena empirically discovered in the Leaky Box models, thus producing less secondary nuclei. Antiprotons will remain almost non-affected due to their high energy threshold of production cross section.

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.

In an attempt to understand the source and propagation of galactic cosmic rays we have employed the Modified Weighted Slab technique along with recent values of the relevant cross sections to compute primary to secondary ratios including B/C and Sub-Fe/Fe for different galactic propagation models. The models that we have considered are the disk-halo diffusion model, the dynamical halo wind model, the turbulent diffusion model and a model with minimal reacceleration. The modified weighted slab technique will be briefly discussed and a more detailed description of the models will be given. We will also discuss the impact that the various models have on the problem of anisotropy at high energy and discuss what properties of a particular model bear on this issue.

It is demonstrated that clusters of galaxies are able to keep cosmic rays for a time exceeding the age of the Universe. This phenomenon reveals itself by the production of the diffuse flux of high energy gamma and neutrino radiation due to the interaction of the cosmic rays with the intracluster gas. The produced flux is determined by the cosmological density of baryons, $\Omega_b$, if a large part of this density is provided by the intracluster gas. The signal from relic cosmic rays has to be compared with the flux produced by the late sources, which can be considered as a background in the search for cosmic ray production in the past. We calculate this flux considering the normal galaxies and AGN in the clusters as the sources of cosmic rays. Another potential cosmic ray source is the shock in the gas accreting to a cluster. We found that this background is relatively high: the diffuse fluxes produced by relic cosmic rays are of the same order of magnitude which can be expected from AGN in the clusters. In all cases the predicted diffuse gamma-ray flux is smaller than the observed one and the diffuse neutrino flux can be seen as the small bump at $E\sim 10^6~GeV$ over the atmospheric neutrino flux. A bright phase in the galaxy evolution can be a source of the relic cosmic rays in clusters, revealing itself by diffuse gamma and neutrino radiations. We found that the observation of a signal from the bright phase is better for an individual cluster.