Using the Auger mass-composition analysis of ultra high energy cosmic rays, based on the shape-fitting of $X_{max}$ distributions, we demonstrate that mass composition and energy spectra measured by Auger, Telescope Array and HiRes can be brought into good agreement. The shape-fitting analysis of $X_{max}$ distributions shows that the measured sum of proton and Helium fractions, for some hadronic-interaction models, can saturate the total flux. Such p+He model, with small admixture of other light nuclei, naturally follows from cosmology with recombination and reheating phases. The most radical assumption of the presented model is the assumed unreliability of the experimental separation of Helium and protons, which allows to consider He/p ratio as a free parameter. The results presented here show that the models with dominant p+He composition explain well the energy spectrum of the dip in the latest (2015 - 2017) data of Auger and Telescope Array, but have some tension at the highest energies with the expected Greisen-Zatsepin-Kuzmin cutoff. The Auger-Prime upgrade experiment has a great potential to reject or confirm this model.

An isotropic component of high energy $\gamma$-ray spectrum measured by Fermi LAT constrains the proton component of UHECR. The strongest restriction comes from the highest, $(580-820)$ GeV, energy bin. One more constraint on the proton component is provided by the IceCube upper bound on ultrahigh energy cosmogenic neutrino flux. We study the influence of these restrictions on the source properties, such as evolution and distribution of sources, their energy spectrum and admixture of nuclei. We also study the sensitivity of restrictions to various Fermi LAT galactic foreground models (model B being less restrictive), to the choice of extragalactic background light model and to overall normalization of the energy spectrum. We claim that the $\gamma$-ray-cascade constraints are stronger than the neutrino ones, and that however many proton models are viable. The basic parameters of such models are relatively large $\gamma_g$ and not very large $z_{\max}$. The allowance for H$e^4$ admixture also relaxes the restrictions. However we foresee that future CTA measurements of $\gamma$-ray spectrum at $E_\gamma \simeq (600 - 800)$ GeV, as well as resolving of more individual $\gamma$-ray sources, may rule out the proton-dominated cosmic ray models.

Using the analytic modeling of the electromagnetic cascades compared with more precise numerical simulations we describe the physical properties of electromagnetic cascades developing in the universe on CMB and EBL background radiations. A cascade is initiated by very high energy photon or electron and the remnant photons at large distance have two-component energy spectrum, $\propto E^{-2}$ ($\propto E^{-1.9}$ in numerical simulations) produced at cascade multiplication stage, and $\propto E^{-3/2}$ from Inverse Compton electron cooling at low energies. The most noticeable property of the cascade spectrum in analytic modeling is 'strong universality', which includes the standard energy spectrum and the energy density of the cascade $\omega_{\rm cas}$ as its only numerical parameter. Using numerical simulations of the cascade spectrum and comparing it with recent Fermi LAT spectrum we obtained the upper limit on $\omega_{\rm cas}$ stronger than in previous works. The new feature of the analysis is "$E_{\max}$ rule". We investigate the dependence of $\omega_{\rm cas}$ on the distribution of sources, distinguishing two cases of universality: the strong and weak ones.

Neutralino annihilation in the Galactic halo is the most definite observational signature proposed for indirect registration of the SUSY Dark Matter (DM) candidate particles. The corresponding annihilation signal (in the form of gamma-rays, positrons and antiprotons) may be boosted for one or three orders of magnitude due to the clustering of cold DM particles into the small-scale and very dense self-gravitating clumps. We discuss the formation of these clumps from the initial density perturbations and their successive fate in the Galactic halo. Only a small fraction of these clumps, $\sim0.1$%, in each logarithmic mass interval $\Delta\log M\sim1$ survives the stage of hierarchical clustering. We calculate the probability of surviving the remnants of dark matter clumps in the Galaxy by modelling the tidal destruction of the small-scale clumps by the Galactic disk and stars. It is demonstrated that a substantial fraction of clump remnants may survive through the tidal destruction during the lifetime of the Galaxy. The resulting mass spectrum of survived clumps is extended down to the mass of the core of the cosmologically produced clumps with a minimal mass. The survived dense remnants of tidally destructed clumps provide an amplification (boosting) of the annihilation signal with respect to the diffuse DM in the Galactic halo. We describe the anisotropy of clump distribution caused by the tidal destruction of clumps in the Galactic disk.

We use a kinetic-equation approach to describe the propagation of ultra high energy cosmic ray protons and nuclei and calculate the expected spectra and mass composition at the Earth for different assumptions on the source injection spectra and chemical abundances. When compared with the spectrum, the elongation rate $X_{max}(E)$ and dispersion $\sigma(X_{max})$ as observed with the Pierre Auger Observatory, several important consequences can be drawn: a) the injection spectra of nuclei must be very hard, $\sim E^{-\gamma}$ with $\gamma\sim 1-1.6$; b) the maximum energy of nuclei of charge $Z$ in the sources must be $\sim 5Z\times 10^{18}$ eV, thereby not requiring acceleration to extremely high energies; c) the fit to the Auger spectrum can be obtained only at the price of adding an {\it ad hoc} light extragalactic component with a steep injection spectrum ($\sim E^{-2.7}$). In this sense, at the ankle ($E_{A}\approx 5\times 10^{18}$ eV) all the components are of extragalactic origin, thereby suggesting that the transition from Galactic to extragalactic cosmic rays occurs below the ankle. Interestingly, the additional light extragalactic component postulated above compares well, in terms of spectrum and normalization, with the one recently measured by KASCADE-Grande.

Small-scale clumps of dark matter are gravitationally bounded structures that have masses comparable to or lower than stellar masses and consist of noninteracting or weakly interacting dark matter particles. In this paper, the current knowledge about the formation and evolution of such structures is reviewed, various types of spectra of primordial cosmological perturbations are considered, and various dark matter models are discussed. Depending on the particular spectrum type, dark matter clumps may differ considerably in their formation processes and ultimate characteristics. The role of clumps in experiments on indirect detection of dark matter particles via their annihilation products is discussed. A number of astrophysical problems and phenomena that are related to dark matter clumps are examined.

We discuss the formation mechanisms and structure of the superdense dark matter clumps (SDMC) and ultracompact minihaloes (UCMH), outlining the differences between these types of DM objects. We define as SDMC the gravitationally bounded DM objects which have come into virial equilibrium at the radiation-dominated (RD) stage of the universe evolution. Such objects can be formed from the isocurvature (entropy) density perturbations or from the peaks in the spectrum of curvature (adiabatic) perturbation. The axion miniclusters (Kolb and Tkachev 1994) are the example of the former model. The system of central compact mass (e. g. in the form of SDMC or primordial black hole (PBH)) with the outer DM envelope formed in the process of secondary accretion we refer to as UCMH. Therefore, the SDMC can serve as the seed for the UCMH in some scenarios. Recently, the SDMC and UCMH were considered in the many works, and we try to systematize them here. We consider also the effect of asphericity of the initial density perturbation in the gravitational evolution, which decreases the SDMC amount and, as the result, suppresses the gamma-ray signal from DM annihilation.

The signatures of Ultra High Energy (E >1 EeV) proton propagation through CMB radiation are pair-production dip and GZK cutoff. The visible characteristics of these two spectral features are ankle, which is intrinsic part of the dip, beginning of GZK cutoff in the differential spectrum and E_{1/2} in integral spectrum. Measured by HiRes and Telescope Array (TA) these characteristics agree with theoretical predictions. However, directly measured mass composition remains a puzzle. While HiRes and TA detectors observe the proton dominated mass composition, the data of Auger detector strongly evidence for nuclei mass composition becoming progressively heavier at energy higher than 4 EeV and reaching Iron at energy about 35 EeV. The models based on the Auger and HiRes/TA data are considered independently and classified using the transition from galactic to extragalactic cosmic rays. The ankle cannot provide this transition. since data of all three detector at energy (1 - 3) EeV agree with pure proton composition (or at least not heavier than Helium). If produced in Galaxy these particles result in too high anisotropy. This argument excludes or strongly disfavours all ankle models with ankle energy E_a > 3 EeV. The calculation of elongation curves, X_{\max}(E), for different ankle models strengthens further this conclusion. Status of other models, the dip, mixed composition and Auger based models are discussed.

The signatures of UHE proton propagation through CMB are pair-production dip and GZK cutoff. The visible manifestations of these spectral features are ankle, beginning of GZK cutoff in the differential spectrum and E_{1/2} in integral spectrum. Observed in all experiments, the ankle is usually interpreted as transition from galactic to extragalactic cosmic rays. Using the mass composition measured by HiRes, Telescope Array (TA) and Auger detectors at energy (1-3) EeV, calculated anisotropy of galactic cosmic rays at these energies, and the elongation curves we strongly argue against the interpretation of the ankle given above. The transition must occur at lower energy, most probably at the second knee as the dip model predicts. The other prediction of this model, the shape of the dip, is well confirmed by HiRes, TA, AGASA and Yakutsk detectors, and, after recalibration of energies, by Auger detector. Predicted beginning of GZK cutoff and E_{1/2} agree well with HiRes and TA data. However, directly measured mass composition remains a puzzle. While HiRes and TA detectors observe the proton-dominated mass composition, as required by the dip model, the data of Auger detector strongly evidence for nuclei mass composition becoming steadily heavier at energy higher than 4 EeV and reaching Iron at energy about 35 EeV. The Auger-based scenario is consistent with another interpretation of the ankle at energy E_a=4 EeV as transition from extragalactic protons to extragalactic nuclei. The heavy- nuclei dominance at higher energies may be provided by low-energy of acceleration for protons E_{max} \sim 4 EeV and rigidity-dependent E_{max}^A =Z E_{max}$ for nuclei. The highest energy suppression may be explained as nuclei-destroying cutoff.

The study of the transition between galactic and extragalactic cosmic rays can shed more light on the end of the Galactic cosmic rays spectrum and the beginning of the extragalactic one. Three models of transition are discussed: ankle, dip and mixed composition models. All these models describe the transition as an intersection of a steep galactic component with a flat extragalactic one. Severe bounds on these models are provided by the Standard Model of Galactic Cosmic Rays according to which the maximum acceleration energy for Iron nuclei is of the order of $E_{\rm Fe}^{\rm max} \approx 1\times 10^{17}$ eV. In the ankle model the transition is assumed at the ankle, a flat feature in the all particle spectrum which observationally starts at energy $E_a \sim (3 - 4)\times 10^{18}$ eV. This model needs a new high energy galactic component with maximum energy about two orders of magnitude above that of the Standard Model. The origin of such component is discussed. As observations are concerned there are two signatures of the transition: change of energy spectra and mass composition. In all models a heavy galactic component is changed at the transition to a lighter or proton component.

We present a systematic study of different methods for the analytic calculation of ultra-high energy nuclei diffuse spectra. Nuclei propagating in the intergalactic space are photo-disintegrated and decrease their Lorentz factor due to the interaction with cosmic microwave background and extragalactic background light. We calculate the evolution trajectories in the backward time, that describe how atomic mass number $A$ and Lorentz factor $\Gamma$ change with redshift $z$. Three methods of spectra calculations are investigated and compared: {\it (i)} trajectory method, {\it(ii)} kinetic equation combined with trajectory calculations and {\it (iii)} coupled kinetic equations. We believe that these three methods exhaust at least the principal possibilities for any analytic solution of the problem. In the most straightforward method {\it(i)} only trajectory calculations are used to connect the observed nuclei flux with the production rate of primary (accelerated) nuclei $A_0$. In the second method {\it (ii)} the flux (space density) of primary nuclei, and secondary nuclei and protons are calculated with the help of kinetic equation and trajectories are used only to determine the generation rates of these nuclei. The third method {\it (iii)} consists in solving the complete set of coupled kinetic equations, written starting with primary nuclei $A_0$, then for $A_0-1$ etc down to the $A$ of interest. The solution of the preceding equation gives the generation rate for the one which follows. An important element of the calculations for all methods is the systematic use of Lorentz factor instead of energy. We consider here the interaction of nuclei only with the cosmic microwave background, this case is particularly suitable for understanding the physical results.

We discuss the problem of ultra high energy nuclei propagation in extragalactic background radiations. The present paper is the continuation of the accompanying paper I where we have presented three new analytic methods to calculate the fluxes and spectra of ultra high energy cosmic ray nuclei, both primary and secondary, and secondary protons. The computation scheme in this paper is based on the analytic solution of coupled kinetic equations, which takes into account the continuous energy losses due to the expansion of the universe and pair-production, together with photo-disintegration of nuclei. This method includes in the most natural way the production of secondary nuclei in the process of photo-disintegration of the primary nuclei during their propagation through extragalactic background radiations. In paper I, in order to present the suggested analytical schemes of calculations, we have considered only the case of the cosmic microwave background radiation, in the present paper we generalize this computation to all relevant background radiations, including infra-red and visible/ultra-violet radiations, collectively referred to as extragalactic background light. The analytic solutions allow transparent physical interpretation of the obtained spectra. Extragalactic background light plays an important role at intermediate energies of ultra high energy cosmic ray nuclei. The most noticeable effect of the extragalactic background light is the low-energy tail in the spectrum of secondary nuclei.

We develop a model of production of the very dense clumps of DM in RD epoch due to the accretion of DM on the loops of cosmic strings as the seeds. At some time the loops disappear, for example due to the gravitational radiation, and the remaining dense clumps produce the enhancement of the annihilation signal. We take into account the velocity distribution of the strings, and consider the two extreme regimes of DM annihilation: fast decay and continuous evaporation. The produced annihilation flux of gamma radiation is detectable, and for some parameters of DM particles and the strings can exceed the extragalactic flux of the gamma-radiation observed by Fermi. For the fixed parameters of DM particles (e.g. neutralino with fixed masses and cross-section of annihilation) one can obtain the limits on the basic string parameter, tension \mu, which is stronger than (more general) limits obtained from WMAP observations, cosmological nucleosynthesis and gravitational lensing. In particular for the neutralino with 100 GeV mass we exclude the interval 5\times10^{-10}<G{\mu}/c^2<5.1\times10^{-9}.

We reconsider the model of neutrino production during the 'bright phase', first suggested in 1977, in the light of modern understanding of the role of Pop III stars and acceleration of particles in supernova shocks. We concentrate on the production of cosmogenic UHE neutrinos in supernova explosions that accompany the death of Pop III stars. Accelerated protons produce neutrinos in collisions with CMB photons. We deliberately use simplified assumptions which make our results transparent. Pop III stars are assumed to be responsible for the reionization of the universe as observed by WMAP. Since the evolution of Pop III stars is much faster than the Hubble rate, we consider the burst of UHE proton production to occur at fixed redshift (z_b=10-20). We discuss the formation of collisionless shocks and particle acceleration in the early universe. The composition of accelerated particles is expected to be proton dominated. A simple calculation is presented to illustrate the fact that the diffuse neutrinos flux from the bright phase burst is concentrated in a relatively narrow range around $7.5 \times 10^{15}(20/z_b)^2$ eV. The $\nu_\mu$ flux may be detectable by IceCube without violating the cascade upper limit and the expected energetics of SNe associated with Pop III stars. A possible signature of the neutrino production from Pop III stars may be the detection of resonant neutrino events. For the burst at $z_b=20$ and $\bar{\nu}_e$-flux at the cascade upper limit, the number of resonant events in IceCube may be as high as 10 events in 5 years of observations. These events have equal energies, $E=6.3\times 10^{15}$ eV, in the form of e-m cascades. Given the large uncertainties in the existing predictions of UHE cosmogenic neutrino fluxes, we argue that neutrinos from the first stars might become one of the most reliable hopes for UHE neutrino astronomy.

Superstring theory and other supersymmetric theories predict the existence of relatively light, weakly interacting scalar particles, called moduli, with a universal form of coupling to matter. Such particles can be emitted from cusps of cosmic strings, where extremely large Lorentz factors are achieved momentarily. Highly boosted modulus bursts emanating from cusps subsequently decay into gluons, they generate parton cascades which in turn produce large numbers of pions and then neutrinos. Due to very large Lorentz factors, extremely high energy neutrinos, up to the Planck scale and above, are produced. For some model parameters, the predicted flux of neutrinos with energies $\gtrsim 10^{21}$ eV is observable by JEM-EUSO and by the future large radio detectors LOFAR and SKA.

Data of Pierre Auger Observatory show a proton-dominated chemical composition of ultrahigh-energy cosmic rays spectrum at (1 - 3) EeV and a steadily heavier composition with energy increasing. In order to explain this feature we assume that (1 - 3) EeV protons are extragalactic and derive their maximum acceleration energy, E_p^{max} \simeq 4 EeV, compatible with both the spectrum and the composition. We also assume the rigidity-dependent acceleration mechanism of heavier nuclei, E_A^{max} = Z x E_p^{max}. The proposed model has rather disappointing consequences: i) no pion photo-production on CMB photons in extragalactic space and hence ii) no high-energy cosmogenic neutrino fluxes; iii) no GZK-cutoff in the spectrum; iv) no correlation with nearby sources due to nuclei deflection in the galactic magnetic fields up to highest energies.

The short review of theoretical aspects of ultra high energy (UHE) neutrinos. The accelerator sources, such as Supernovae remnants, Gamma Ray Bursts, AGN etc are discussed. The top-down sources include Topological Defects (TDs), Superheavy Dark Matter (SHDM) and Mirror Matter. The diffuse fluxes are considered accordingly as that of cosmogenic and top-down neutrinos. Much attention is given to the cascade upper limit to the diffuse neutrino fluxes in the light of Fermi-LAT data on diffuse high energy gamma radiation. This is most general and rigorous upper limit, valid for both cosmogenic and top-down models. At present upper limits from many detectors are close to the cascade upper limit, and 5 yr IceCube upper limit will be well below it.

We develop a model for explaining the data of Pierre Auger Observatory (Auger) for Ultra High Energy Cosmic Rays (UHECR), in particular, the mass composition being steadily heavier with increasing energy from 3 EeV to 35 EeV. The model is based on the proton-dominated composition in the energy range (1 - 3) EeV observed in both Auger and HiRes experiments. Assuming extragalactic origin of this component, we argue that it must disappear at higher energies due to a low maximum energy of acceleration, E_p^{\max} \sim (4 - 10) EeV. Under an assumption of rigidity acceleration mechanism, the maximum acceleration energy for a nucleus with the charge number Z is ZE_p^{\max}, and the highest energy in the spectrum, reached by Iron, does not exceed (100 - 200) EeV. The growth of atomic weight with energy, observed in Auger, is provided by the rigidity mechanism of acceleration, since at each energy E=ZE_p^{\max} the contribution of nuclei with Z' < Z vanishes. The described model has disappointing consequences for future observations in UHECR: Since average energies per nucleon for all nuclei are less than (2 - 4) EeV, (i) pion photo-production on CMB photons in extragalactic space is absent; (ii) GZK cutoff in the spectrum does not exist; (iii) cosmogenic neutrinos produced on CMBR are absent; (iv) fluxes of cosmogenic neutrinos produced on infrared - optical background radiation are too low for registration by existing detectors and projects. Due to nuclei deflection in galactic magnetic fields, the correlation with nearby sources is absent even at highest energies.

Ultrahigh energy cosmic ray (UHECR) protons interacting with the cosmic microwave background (CMB) produce UHE electrons and gamma-rays that in turn initiate electromagnetic cascades on CMB and infrared photons. As a result, a background of diffuse isotropic gamma radiation is accumulated in the energy range $E\lsim 100$ GeV. The Fermi-LAT collaboration has recently reported a measurement of the extragalactic diffuse background finding it less intense and softer than previously measured by EGRET. We show that this new result constrains UHECR models and the flux of cosmogenic neutrinos. In particular, it excludes models with cosmogenic neutrino fluxes detectable by existing neutrino experiments, while next-generation detectors as e.g. JEM-EUSO can observe neutrinos only for extreme parameters.

Superconducting cosmic strings naturally emit highly boosted charge carriers from cusps. This occurs when a cosmic string or a loop moves through a magnetic field and develops an electric current. The charge carriers and the products of their decay, including protons, photons and neutrinos, are emitted as a narrow jets with opening angle $\theta \sim 1/\gamma_c$, where $\gamma_c$ is the Lorentz factor of the cusp. The excitation of electric currents in strings occurs mostly in clusters of galaxies, which are characterized by magnetic fields $B \sim 10^{-6}$ G and a filling factor $f_B \sim 10^{-3}$. Two string parameters determine the emission of the particles: the symmetry breaking scale $\eta$, which for successful applications should be of order $10^9$--$10^{12}$ GeV, and the dimensionless parameter $i_c$, which determines the maximum induced current as $J_{max} =i_c e \eta$ and the energy of emitted charge carriers as $\epsilon_x \sim i_c \gamma_c \eta$, where $e$ is the electric charge of a particle. For the parameters $\eta $ and $B$ mentioned above, the Lorentz factor reaches $\gamma_c \sim 10^{12}$ and the maximum particle energy can be as high as $\gamma_c\eta \sim 10^{22}$ GeV. The diffuse fluxes of UHE neutrinos are close to the cascade upper limit, and can be detected by future neutrino observatories. The signatures of this model are: very high energies of neutrinos, in excess of $10^{20}$ eV, correlation of neutrinos with clusters of galaxies, simultaneous appearance of several neutrino-produced showers in the field of view of very large detectors, such as JEM-EUSO, and 10 TeV gamma radiation from the Virgo cluster. The flux of UHE protons from cusps may account for a large fraction of the observed events at the highest energies.

We discuss the superluminal problem in the diffusion of ultra high energy protons with energy losses taken into account. The phenomenological solution of this problem is found with help of the generalized J\"uttner propagator, originally proposed for relativization of the Maxwellian gas distribution. It is demonstrated that the generalized J\"uttner propagator gives the correct expressions in the limits of diffusive and rectilinear propagation of the charged particles in the magnetic fields, together with the intermediate regime, in all cases without superluminal velocities. This solution, very general for the diffusion, is considered for two particular cases: diffusion inside the stationary objects, like e.g. galaxies, clusters of galaxies etc, and for expanding universe. The comparison with the previously obtained solutions for propagation of UHE protons in magnetic fields is performed.

Superheavy dark matter (SHDM) exchanges energy with its environment much slower than particles with masses close to the electroweak (EW) scale and has therefore different small-scale clustering properties. Using the neutralino as candidate for the SHDM, we find that free-streaming allows the formation of DM clumps of all masses down to $\sim 260 m_\chi$ in the case of bino. If small-scale clumps evolve from a non-standard, spiky spectrum of perturbations, DM clumps may form during the radiation dominated era. These clumps are not destroyed by tidal interactions and can be extremely dense. In the case of a bino, a "gravithermal catastrophe" can develop in the central part of the most dense clumps, increasing further the central density and thus the annihilation signal. In the case of a higgsino, the annihilation signal is enhanced by the Sommerfeld effect. As a result annihilations of superheavy neutralinos in dense clumps may lead to observable fluxes of annihilation products in the form of ultrahigh energy particles, for both cases, higgsinos and binos, as lightest supersymmetric particles.

The formation and evolution of superdense clumps (or subhalos) is studied. Such clumps of dark matter (DM) can be produced by many mechanisms, most notably by spiky features in the spectrum of inflationary perturbations and by cosmological phase transitions. Being produced very early during the radiation dominated epoch, superdense clumps evolve as isolated objects. They do not belong to hierarchical structures for a long time after production, and therefore they are not destroyed by tidal interactions during the formation of larger structures. For DM particles with masses close to the electroweak (EW) mass scale, superdense clumps evolve towards a power-law density profile $\rho(r) \propto r^{-1.8}$ with a central core. Superdense clumps cannot be composed of standard neutralinos, since their annihilations would overproduce the diffuse gamma radiation. If the clumps are constituted of superheavy DM particles and develop a sufficiently large central density, the evolution of their central part can lead to a 'gravithermal catastrophe.' In such a case, the initial density profile turns into an isothermal profile with $\rho \propto r^{-2}$ and a new, much smaller core in the center. Superdense clumps can be bserved by gamma radiation from DM annihilations and by gravitational wave detectors, while the production of primordial black holes and cascade nucleosynthesis constrain this scenario.

UHE neutrinos with $E>10^{17}$ eV can be produced by ultra-high energy cosmic rays (UHECR) interacting with CMB photons (cosmogenic neutrinos) and by top-down sources, such as topological defects (TD), superheavy dark matter (SHDM) and mirror matter. Cosmogenic neutrinos are reliably predicted and their fluxes can be numerically evaluated using the observed flux of UHECR. The lower limit for the flux is obtained for the case of pure proton composition of the observed UHECR. The rigorous upper limit for cosmogenic neutrino flux also exists. The maximum neutrino energy is determined by maximum energy of acceleration, which at least for the shock acceleration is expected not to exceed $10^{21} - 10^{22}$ eV. The top-down sources provide neutrino energies a few orders of magnitude higher, and this can be considered as a signature of these models. Oscillations play important role in UHE neutrino astronomy. At production of cosmogenic neutrinos $\tau$-neutrinos are absent and $\bar{\nu}_e$ neutrinos are suppressed. These species, important for detection, appear in the observed fluxes due to oscillation. Mirror neutrinos cannot be observed directly, but due to oscillations to ordinary neutrinos they can provide the largest neutrino flux at the highest energies.

The status of the Greisen-Zatsepin-Kuzmin (GZK) cutoff and pair-production dip in Ultra High Energy Cosmic Rays (UHECR) is discussed.They are the features in the spectrum of protons propagating through CMB radiation in extragalactic space, and discovery of these features implies that primary particles are mostly extragalactic protons. The spectra measured by AGASA, Yakutsk, HiRes and Auger detectors are in good agreement with the pair-production dip, and HiRes data have strong evidences for the GZK cutoff. The Auger spectrum,as presented at the 30th ICRC 2007, agrees with the GZK cutoff, too. The AGASA data agree well with the beginning of the GZK cutoff at E \leq 80 EeV, but show the excess of events at higher energies, the origin of which is not understood. The difference in the absolute fluxes measured by different detectors disappears after energy shift within the systematic errors of each experiment.

What happened to the central cores of tidally destructed dark matter clumps in the Galactic halo? We calculate the probability of surviving of the remnants of dark matter clumps in the Galaxy by modelling the tidal destruction of the small-scale clumps. It is demonstrated that a substantial fraction of clump remnants may survive through the tidal destruction during the lifetime of the Galaxy if the radius of a core is rather small. The resulting mass spectrum of survived clumps is extended down to the mass of the core of the cosmologically produced clumps with a minimal mass. Since the annihilation signal is dominated by the dense part of the core, destruction of the outer part of the clump affects the annihilation rate relatively weakly and the survived dense remnants of tidally destructed clumps provide a large contribution to the annihilation signal in the Galaxy. The uncertainties in minimal clump mass resulting from the uncertainties in neutralino models are discussed.

We discuss the production of ultra high energy secondary protons by cosmic ray primary nuclei propagating in the intergalactic space through Cosmic Microwave Background (CMB) and Infrared (IR) radiations. Under the assumption that only primary nuclei with a fixed atomic mass number $A_0$ are accelerated, the spectrum of secondary protons is calculated. It is found that for all $A_0$ the diffuse flux of secondary protons starts to dominate over that of primary nuclei at energy $E \sim (1 - 2)\times 10^{19}$ eV, and thus the standard Greisen-Zatsepin -Kuzmin (GZK) cutoff is produced.

Puzzles often give birth to the great discoveries, the false discoveries sometimes stimulate the exiting ideas in theoretical physics. The historical examples of both are described in Introduction and in section ``Cosmological Puzzles''. From existing puzzles most attention is given to Ultra High Energy Cosmic Ray (UHECR) puzzle and to cosmological constant problem. The 40-years old UHECR problem consisted in absence of the sharp steepening in spectrum of extragalactic cosmic rays caused by interaction with CMB radiation. This steepening is known as Greisen-Zatsepin-Kuzmin (GZK) cutoff. It is demonstrated here that the features of interaction of cosmic ray protons with CMB are seen now in the spectrum in the form of the dip and beginning of the GZK cutoff. The most serious cosmological problem is caused by large vacuum energy of the known elementary-particle fields which exceeds at least by 45 orders of magnitude the cosmological vacuum energy. The various ideas put forward to solve this problem during last 40 years, have weaknesses and cannot be accepted as the final solution of this puzzle. The anthropic approach is discussed.

We discuss the signatures of the transition from galactic to extragalactic cosmic rays in different scenarios, giving most attention to the dip scenario. The dip is a feature in the diffuse spectrum of ultra-high energy (UHE) protons in the energy range $1\times 10^{18} - 4\times 10^{19}$ eV, which is caused by electron-positron pair production on the cosmic microwave background (CMB) radiation. The dip scenario provides a simple physical description of the transition from galactic to extragalactic cosmic rays. Here we summarize the signatures of the pair production dip model for the transition, most notably the spectrum, the anisotropy and the chemical composition. The main focus of our work is however on the description of the features that arise in the elongation rate and in the distribution of the depths of shower maximum $X_{\rm max}$ in the dip scenario. We find that the curve for $X_{\max}(E)$ shows a sharp increase with energy, which reflects a sharp transition from an iron dominated flux at low energies to a proton dominated flux at $E\sim 10^{18}$ eV. We also discuss in detail the shape of the $X_{\max}$ distributions for cosmic rays of given energy and demonstrate that this represents a powerful tool to discriminate between the dip scenario and other possible models of the transition.

The transition from galactic to extragalactic cosmic rays is discussed. One of critical indications for transition is given by the Standard Model of Galactic cosmic rays, according to which the maximum energy of acceleration for iron nuclei is of order of $E_{\rm Fe}^{\rm max} \approx 1\times 10^{17}$ eV. At $E > E_{\rm Fe}^{\rm max}$ the spectrum is predicted to be very steep and thus the Standard Model favours the transition at energy not much higher than $E_{\rm Fe}^{\rm max}$. As observations are concerned there are two signatures of transition: change of energy spectra and elongation rate (depth of shower maximum in the atmosphere $X_{\rm max}$ as function of energy). Three models of transition are discussed: dip-based model, mixed composition model and ankle model. In the latter model the transition occurs at the observed spectral feature, ankle, which starts at $E_a \approx 1\times 10^{19}$ eV and is characterised by change of mass compostion from galactic iron to extragalactic protons. In the dip model the transition occures at the second knee observed at energy $(4 -8)\times 10^{17}$ eV and is characterised by change of mass composition from galactic iron to extragalactic protons. The mixed composition model describes transition at $E \sim 3\times 10^{18}$ eV with mass composition changing from galactic iron to extragactic mixed composition of different nuclei. These models are confronted with observational data on spectra and elongation rates from different experiments, including Auger.

The status of the pair-production dip as a spectral feature, produced by interaction of Ultra High Energy extragalactic protons with CMB is discussed.

We describe the anisotropy of dark matter clump distribution caused by tidal destruction of clumps in the Galactic disk. A tidal destruction of clumps with orbit planes near the disk plane occurs more efficiently as compared with destruction of clumps at near-polar orbits. A corresponding annihilation of dark matter particles in small-scale clumps produces the anisotropic gamma-ray signal with respect to the Galactic disk. This anisotropy is rather small, 9%, and superimposed on that due to off-centering position of the Sun in the Galaxy. The anisotropy of annihilation signal with respect to the Galactic disk provides a possibility to discriminate dark matter annihilation from the diffuse gamma-ray backgrounds of other origin.

We study the solution of the diffusion equation for Ultra-High Energy Cosmic Rays in the general case of an expanding universe, comparing it with the well known Syrovatsky solution obtained in the more restrictive case of a static universe. The formal comparison of the two solutions with all parameters being fixed identically reveals an appreciable discrepancy. This discrepancy is less important if in both models a different set of best-fit parameters is used.

We consider the astrophysical implications of the diffusion equation solution for Ultra-High Energy Cosmic Rays (UHECR) in the expanding universe, obtained in paper I (V.Berezinsky & A.Gazizov, ApJ 643 (2006) 8). The UHECR spectra are calculated in a model with sources located in vertices of the cubic grid with a linear constant (source separation) d. The calculations are performed for various magnetic field configurations (B_c,l_c), where l_c is the basic scale of the turbulence and B_c is the coherent magnetic field on this scale. The main purpose of these calculations is to demonstrate the validity of the solution obtained in paper I and to compare this solution with the Syrovatsky solution used in previous works. The Syrovatsky solution must be necessarily embedded in the static cosmological model. The formal comparison of the two solutions with all parameters being fixed identically reveals the appreciable discrepancies between two spectra. These discrepancies are less if in both models the different sets of the best-fit parameters are used.

We argue that an astrophysical solution to UHECR problem is viable. The pectral features of extragalactic protons interacting with CMB are calculated in model-independent way. Using the power-law generation spectrum $\propto E^{-\gamma_g}$ as the only assumption, we analyze four features of the proton spectrum: the GZK cutoff, dip, bump and the second dip. We found the dip, induced by electron-positron production on CMB, as the most robust feature, existing in energy range $1\times 10^{18} - 4\times 10^{19}$ eV. Its shape is stable relative to various phenomena included in calculations. The dip is well confirmed by observations of AGASA, HiRes, Fly's Eye and Yakutsk detectors. The best fit is reached at $\gamma_g =2.7$, with the allowed range 2.55 - 2.75. The dip is used for energy calibration of the detectors. After the energy calibration the fluxes and spectra of all three detectors agree perfectly, with discrepancy between AGASA and HiRes at $E> 1\times 10^{20}$ eV being not statistically significant. The agreement of the dip with observations should be considered as confirmation of UHE proton interaction with CMB. The dip has two flattenings. The high energy flattening at $E \approx 1\times 10^{19}$ eV automatically explains ankle. The low-energy flattening at $E \approx 1\times 10^{18}$ eV provides the transition to galactic cosmic rays. This transition is studied quantitatively. The UHECR sources, AGN and GRBs, are studied in a model-dependent way, and acceleration is discussed. Based on the agreement of the dip with existing data, we make the robust prediction for the spectrum at $1\times 10^{18} - 1\times 10^{20}$ eV to be measured in the nearest future by Auger detector.

We discuss as a new signature for the interaction of extragalactic ultrahigh energy protons with cosmic microwave background radiation a spectral feature located at $E=6.3\times 10^{19}$ eV in the form of a narrow and shallow dip. It is produced by the interference of $e^+e^-$-pair and pion production. We show that this dip and in particular its position are almost model-independent. Its observation by future ultrahigh energy cosmic ray detectors may give the conclusive confirmation that an observed steepening of the spectrum is caused by the Greisen-Zatsepin-Kuzmin effect.

The dip is a feature in the diffuse spectrum of ultra-high energy (UHE) protons caused by electron-positron pair production on the cosmic microwave background (CMB) radiation. For a power-law generation spectrum $E^{-2.7}$, the calculated position and shape of the dip is confirmed with high accuracy by the spectra observed by the Akeno-AGASA, HiRes, Yakutsk and Fly's Eye detectors. When the particle energies, measured in these detectors, are calibrated by the dip, their fluxes agree with a remarkable accuracy. The predicted shape of the dip is quite robust. The dip is only modified strongly when the fraction of nuclei heavier than protons is high at injection, which imposes some restrictions on the mechanisms of acceleration operating in UHECR sources. The existence of the dip, confirmed by observations, implies that the transition from galactic to extragalactic cosmic rays occurs at $E \lsim 1\times 10^{18}$ eV. We show that at energies lower than a characteristic value $E_{\rm cr}\approx 1\times 10^{18}$ eV, the spectrum of extragalactic cosmic rays flattens in all cases of interest, and it provides a natural transition to a steeper galactic cosmic ray spectrum. This transition occurs at some energy below $E_{\rm cr}$, corresponding to the position of the so-called second knee. We discuss extensively the constraints on this model imposed by current knowledge of acceleration processes and sources of UHECR and compare it with the traditional model of transition at the ankle.

We describe the formation of a seed massive black hole (MBH) inside a supermassive star (SMS) in a distant galactic nucleus. The short-lived SMS is naturally formed due to collision destructions of normal stars in the evolving galactic nucleus. The neutron stars (NSs) and stellar-mass black holes form a compact self-gravitating subsystem deep inside a SMS. This subsystem is short-lived in comparison with a host SMS and collapses finally into the MBH. Just before gravitational collapse of compact subsystem the frequent NS collisions are accompanied by the generation of numerous ultra-relativistic fireballs. A combined ram pressure of multiple coexisting fireballs produces a quasi-stationary rarefied cavity in the central part of SMS. The protons are accelerated in the fireballs and by relativistic shocks in the cavity. All secondary particles, produced in collisions, except the high-energy neutrinos are absorbed in the SMS interiors. An estimated high-energy neutrino signal from this hidden source can be detected by the neutrino telescope with an effective area 1 km^2 providing the evidence for MBH formation in a distant galactic nucleus. A corresponding lifetime of this high-energy hidden neutrino source is ~0.1-1 yr.

We present an analytic solution to diffusion equation for high energy cosmic rays in the expanding universe. The particles are assumed to be ultra-relativistic and they can have energy losses arbitrarily dependent on energy and time. The obtained solution generalizes the Syrovatsky solution, valid for the case when energy losses and diffusion coefficient are time-independent.

Superheavy particles are a natural candidate for the dark matter in the universe and our galaxy, because they are produced generically during inflation in cosmologically interesting amounts. The most attractive model for the origin of superheavy dark matter (SHDM) is gravitational production at the end of inflation. The observed cosmological density of dark matter determines the mass of the SHDM particle as $m_X=$(a few)$\:\times 10^{13}$ GeV, promoting it to a natural candidate for the source of the observed ultra-high energy cosmic rays (UHECR). After a review of the theoretical aspects of SHDM, we up-date its predictions for UHECR observations: no GZK cutoff, flat energy spectrum with $dN/dE\approx 1/E^{1.9}$, photon dominance and galactic anisotropy. We analyze the existing data and conclude that SDHM as explanation for the observed UHECRs is at present disfavored but not yet excluded. We calculate the anisotropy relevant for future Auger observations that should be the conclusive test for this model. Finally, we emphasize that negative results of searches for SHDM in UHECR do not disfavor SHDM as a dark matter candidate. Therefore, UHECRs produced by SHDM decays and with the signatures as described should be searched for in the future as subdominant effect.

A mass function of small-scale dark matter clumps is calculated in the standard cosmological scenario with an inflationary-produced primordial fluctuation spectrum and with a hierarchical clustering. We take into account the tidal destruction of clumps at early stages of structure formation starting from a time of clump detachment from the Universe expansion. Only a small fraction of these clumps, ~0.1-0.5%, in each logarithmic mass interval $\Delta\log M\sim1$ survives the stage of hierarchical clustering. The surviving clumps can be disrupted further in the galaxies by tidal interactions with stars. We performed the detailed calculations of the tidal destruction of clumps by stars in the Galactic bulge and halo and by the Galactic disk itself. It is shown that the Galactic disc provides the dominant contribution to the tidal destruction of small-scale clumps outside the bulge. The results obtained are crucial for calculations of the dark matter annihilation signal in the Galaxy.

The sources and fluxes of superGZK neutrinos, $E>10^{20}$ eV, are discussed. The fluxes of {\em cosmogenic neutrinos}, i.e. those produced by ultra-high energy cosmic rays (UHECR) interacting with CMB photons, are calculated in the models, which give the good fit to the observed flux of UHECR. The best fit given in no-evolutionary model with maximum acceleration energy $E_{\rm max}=1\times 10^{21}$ eV results in very low flux of superGZK neutrinos an order of magnitude lower than the observed flux of UHECR. The predicted neutrino flux becomes larger and observable by next generation detectors at energies $10^{20} - 10^{21}$ eV in the evolutionary models with $E_{\rm max}=1\times 10^{23}$ eV. The largest cosmogenic neutrino flux is given in models with very flat generation spectrum, e.g. $\propto E^{-2}$. The neutrino energies are naturally high in the models of {\em superheavy dark matter and topological defects}. Their fluxes can also be higher than those of cosmogenic neutrinos. The largest fluxes are given by {\em mirror neutrinos}, oscillating into ordinary neutrinos. Their fluxes obey some theoretical upper limit which is very weak, and in practice these fluxes are most efficiently limited now by observations of radio emission from neutrino-induced showers.

The dip is a feature in the diffuse spectrum of UHE protons in energy range $1\times 10^{18} - 4\times 10^{19}$ eV, which is caused by electron-positron pair production on CMB photons. Calculated for power-law generation spectrum with index $\gamma_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 discretenesses in the source distribution, for local source overdensity and deficit, for source inhomogeneities on scale $\ell \lsim 100$ Mpc etc. Below the characteristic energy $E_c \approx 1\times 10^{18}$ 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_{\rm tr} < E_c$ approximately coincides with the position of the second knee $E_{2kn}$ observed in the cosmic ray spectrum. The critical energy $E_c$ is determined by the energy $E_{\rm eq} = 2.3\times 10^{18}$ eV, where adiabatic and pair-production energy losses are equal. Thus, position of the second knee is explained by proton energy losses on CMB photons.

In this paper we discuss the anti-GZK effect that arises in the framework of the diffusive propagation of Ultra High Energy (UHE) protons. This effect consists in a jump-like increase of the maximum distance from which UHE protons can reach the observer. The position of the jump is independent of the Intergalactic Magnetic Field (IMF) strength and depends only on the energy losses of protons, namely on the transition energy from adiabatic and pair-production energy losses. The Ultra High Energy Cosmic Rays (UHECR) spectrum presents a low-energy steepening approximately at this energy, which is very close to the position of the observed second knee. The dip, seen in the universal spectrum as a signature of the proton interaction with the Cosmic Microwave Background (CMB) radiation, is also present in the case of diffusive propagation in magnetic fields.

We discuss the antiGZK effect in the diffusive propagation of ultra high energy protons in intergalactic magnetic fields, which consists in a jump-like increase of the maximum distance from which ultra high energy protons can reach an observer. The position of this jump, $E_j \approx 2\times 10^{18}$ eV, is determined exclusively by energy losses (transition from adiabatic to pair-production energy losses) and it is independent of the diffusion parameters. The diffuse spectrum presents a low-energy steepening approximately at this energy, which is very close to the position of the second knee observed in the cosmic ray spectrum. The dip, seen in the universal spectrum as a signature of the interaction with the cosmic microwave background radiation, is also present in the case of diffusive propagation in magnetic fields.

The short review of theoretical aspects of ultra high energy (UHE) neutrinos and superGZK neutrinos. The sources and diffuse fluxes of UHE neutrinos are discussed. Much attention is given to comparison of the cascade and cosmic ray upper bounds for diffuse neutrino fluxes. Cosmogenic neutrinos and neutrinos from the mirror mater are considered as superGZK neutrinos.

We study the extragalactic protons with universal spectrum, which is independent of mode of propagation, when distance between sources is less than the propagation lengths, such as energy attenuation length or diffusion length (for propagation in magnetic fields). The propagation features in this spectrum, the GZK cutoff, dip and bump, are studied with help of modification factor, which weakly depends on the generation spectrum index $\gamma_g$. We argue that from the above features the dip is the most model-independent one. For the power-law generation spectrum with $\gamma_g=2.7$ the dip is very well confirmed by the data of all existing detectors, which gives the strong evidence for extragalactic protons propagating through CMB. We develop the AGN model for origin of UHECR, which successfully explains the observed spectra up to $1\times 10^{20}$ eV and transition from galactic to extragalactic cosmic rays. The calculated spectrum has the GZK cutoff, and the AGASA excess of events at $E \gsim 1\times 10^{20}$ eV needs another component, e.g. from superheavy dark matter. In case of weak extragalactic magnetic fields this model is consistent with small-angle clustering and observed correlation with BL Lacs.

Ultrahigh energy (UHE) extragalactic protons propagating through cosmic microwave radiation (CMB) acquire the spectrum features in the form of the dip, bump (pile-up protons) and the Greisen-Zatsepin-Kuzmin (GZK) cutoff. We have performed the analysis of these features in terms of the modification factor. This analysis is weakly model-dependent, especially in case of the dip. The energy shape of the dip is confirmed by the Akeno-AGASA data with $\chi^2=19.06$ for d.o.f. =17 with two free parameters used for comparison. The agreement with HiRes data is also very good. This is the strong evidence that UHE cosmic rays observed at energies $1\times 10^{18}$ eV - $4\times 10^{19}$ eV are extragalactic protons propagating through CMB. The dip is also present in case of diffusive propagation in magnetic field.

We study the cosmological origin of small-scale DM clumps in the hierarchical scenario with the most conservative assumption of adiabatic Gaussian fluctuations. The mass spectrum of small-scale clumps with M<10^3Msun is calculated with tidal destruction of the clumps taken into account within the hierarchical model of clump structure. Only 0.1-0.5% of small clumps survive the stage of tidal destruction in each logarithmic mass interval. The mass distribution of clumps has a cutoff at Mmin due to diffusion of DM particles out of a fluctuation and free streaming at later stage. Mmin is a model dependent quantity. In the case the neutralino DM, considered as a pure bino, Mmin~10^-8 Msun. The evolution of density profile in a DM clump does not result in the singularity because of formation of the core under influence of tidal interaction. The radius of the core is ~0.1R, where R is radius of the clump. The applications for annihilation of DM particles in the Galactic halo are studied. The number density of clumps as a function of their mass, radius and distance to the Galactic center is presented. The enhancement of annihilation signal due to clumpiness, valid for arbitrary DM particles, is calculated. In spite of small survival probability, the global annihilation signal in most cases is dominated by clumps, with major contribution given by small clumps. The enhancement due to large clumps with M>10^6 Msun is very small.

We present a detailed analytical study of the propagation of ultra high energy (UHE) particles in extragalactic magnetic fields. The crucial parameter which affects the diffuse spectrum is the separation between sources. In the case of a uniform distribution of sources with a separation between them much smaller than all characteristic propagation lengths, the diffuse spectrum of UHE particles has a {\em universal} form, independent of the mode of propagation. This statement has a status of theorem. The proof is obtained using the particle number conservation during propagation, and also using the kinetic equation for the propagation of UHE particles. This theorem can be also proved with the help of the diffusion equation. In particular, it is shown numerically, how the diffuse fluxes converge to this universal spectrum, when the separation between sources diminishes. We study also the analytic solution of the diffusion equation in weak and strong magnetic fields with energy losses taken into account. In the case of strong magnetic fields and for a separation between sources large enough, the GZK cutoff can practically disappear, as it has been found early in numerical simulations. In practice, however, the source luminosities required are too large for this possibility.

In this work we present a detailed computation of the spectra of UHECR in the top-down scenario. We compare the spectra of hadrons obtained by two different methods in QCD and supersymmetric (SUSY) QCD with large primary energies $\sqrt{s}$ up to $10^{16}$ GeV. The two methods discussed are a Monte Carlo (MC) simulation and the evolution of the hadron fragmentation functions as described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. The hadron spectra obtained by the two methods agree fairly well in the interesting energy range $10^{-5}M_{X}<E<0.3M_{X}$ ($M_X$ is the energy scale of the process $M_{X}\ge 10^{12}$ GeV). We have also computed the spectra of photons, neutrinos and nucleons obtaining a good agreement with other published results. The consistency of the spectra computed by different methods allows us to consider the spectral shape as a signature of the production model for UHECR, such as the decay of super heavy relic particles or topological defects.

We investigate the problem of transition from Galactic cosmic rays to extragalactic ultra high energy cosmic rays. Using the model for extragalactic ultra high energy cosmic rays and observed all-particle cosmic ray spectrum, we calculate the galactic spectrum of iron nuclei in the energy range $10^8 - 10^9$ GeV. The flux and spectrum predicted at lower energies agree well with the KASCADE data. The transition from Galactic to extragalactic cosmic rays is distinctly seen in spectra of protons and iron nuclei, when they are measured separately. The shape of predicted iron spectrum agrees with the Hall diffusion.

We present results from two different methods for the calculation of hadron spectra in QCD and SUSY QCD with large primary energies $\sqrt{s}$ up to $10^{16}$ GeV. The two methods considered are a Monte Carlo (MC) simulation and the evolution of fragmentation functions described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. We find that the pion, nucleon and all-hadron spectra calculated with the two methods agree well. The MC simulation is performed with new hadronization functions (in comparison with our previous work), motivated by low-energy ($\sqrt s < M_Z$) data and DGLAP. The hadron spectra calculated with both sets of hadronization functions agree well, which indicates that our method for calculating the hadronization function works successfully. The small difference in the calculated hadron spectra characterizes the uncertainties of this method. We calculate also the spectra of photons, neutrinos and nucleons and compare them with other published results. The agreement is good for all $x$ from $\sim 10^{-5}$ up to $x\leq 0.3$. The consistency of the spectra calculated by different methods allows to consider the spectral shape as a signature of models with decays or annihilations of superheavy particles, such as topological defects or superheavy DM. The UHECR spectra from these sources are calculated.

Production of small-scale DM clumps is studied in the standard cosmological scenario with an inflation-produced primeval fluctuation spectrum. Special attention is given to three following problems: (i) The mass spectrum of small-scale clumps with $M \lesssim 10^3 M_{\odot}$ is calculated with tidal destruction of the clumps taken into account within the hierarchical model of clump structure. Only 0.1 - 0.5% of small clumps survive the stage of tidal destruction in each logarithmic mass interval $\Delta\ln M\sim1$. (ii) The mass distribution of clumps has a cutoff at $M_{\rm min}$ due to diffusion of DM particles out of a fluctuation and free streaming at later stage. $M_{\rm min}$ is a model dependent quantity. In the case the neutralino, considered as a pure bino, is a DM particle, $M_{\rm min} \sim 10^{-8} M_{\odot}$. (iii) The evolution of density profile in a DM clump does not result in the singularity because of formation of the core under influence of tidal interaction. The radius of the core is $R_c \sim 0.1 R$, where $R$ is radius of the clump. The applications for annihilation of DM particles in the Galactic halo are studied. The number density of clumps as a function of their mass, radius and distance to the Galactic center is presented. The enhancement of annihilation signal due to clumpiness, valid for arbitrary DM particles, is calculated. In spite of small survival probability, the annihilation signal in most cases is dominated by clumps. For observationally preferable value of index or primeval fluctuation spectrum $n_p \approx 1$, the enhancement of annihilation signal is described by factor 2 - 5 for different density profiles in a clump.

We demonstrate that the energy spectra of Ultra High Energy Cosmic Rays (UHECR) as observed by AGASA, Fly's Eye, HiRes and Yakutsk detectors, have the imprints of UHE proton interaction with the CMB radiation as the dip centered at $E\sim 1\times 10^{19}$ eV, beginning of the GZK cutoff, and very good agreement with calculated spectrum shape. This conclusion about proton composition agrees with recent HiRes data on elongation rate that support the proton composition at $E\geq 1\times 10^{18}$ eV. The visible bump in the spectrum at $E \sim 4\times 10^{19}$ eV is not caused by pile-up protons, but is an artifact of multiplying the spectrum by $E^3$. We argue that these data, combined with small-angle clustering and correlation with AGN (BL Lacs), point to the AGN model of UHECR origin at energies $E \leq 1\times 10^{20}$ eV. The events at higher energies and the excess of the events at $E \geq 1\times 10^{20}$ eV, which is observed by AGASA (but absent in the HiRes data) must be explained by another component of UHECR, e.g. by UHECR from superheavy dark matter.

The sources and fluxes of superGZK neutrinos, $E>10^{20}$ eV, are discussed. The most promising sources are reionization bright phase, topological defects, superheavy dark matter and mirror matter. The energy of neutrinos can be above the GUT scale ($\sim 10^{16}$ GeV). The predicted fluxes are observable by future space detectors EUSO and OWL.

We demonstrate that the energy spectra of Ultra High Energy Cosmic rays (UHECR) as observed by AGASA, Fly's Eye, HiRes and Yakutsk detectors, have the imprints of UHE proton interaction with the CMB radiation in the form of the dip at $E\sim 1\times 10^{19}$~ eV, of the beginning of the GZK cutoff, and of very good agreement with calculated spectrum shape. We argue that these data, combined with small-angle clustering and correlation with AGN (BL Lacs), point to the AGN model of UHECR origin at energies $E \lsim 1\times 10^{20}$ eV. The excess of the events at $E \gsim 1\times 10^{20}$ eV, which is observed by AGASA (but absent in HiRes data) can be explained by another component of UHECR, e.g. by UHECR from superheavy dark matter.

We study decays of superheavy particles $X$ into leptons. We show that they initiate cascades similar to QCD parton jets, if $m_X\gsim 10^6$ GeV. Electroweak cascading is studied and the energy spectra of the produced leptons are calculated in the framework of a broken SU(2) model of weak interactions. As application, important for the Z-burst model for ultrahigh energy cosmic rays, we consider decays of superheavy particles coupled on tree-level only to neutrinos and derive stringent limit for these decays from the observed diffuse extragalactic $\gamma$-ray flux.

The energy losses and spectra of Ultra High Energy Cosmic Rays (UHECR) are calculated for protons as primary particles. The attention is given to the energy losses due to electron-positron production in collisions with the microwave 2.73 K photons. The energy spectra are calculated for several models, which differ by production spectra and by source distribution, namely: (i) Uniform distribution of the sources with steep generation spectra with indices 2.4 - 2.7, with cosmological evolution and without it. In this case it is possible to fit the shape of the observed spectrum up to 8.10^{19} eV, but the required CR emissivity is too high and the GZK cutoff is present. (ii) Uniform distribution of the sources with flat generation spectrum dE/E^2 which is relevant to GRBs. The calculated spectrum is in disagreement with the observed one. The agreement at E\lesssim 8.10^{19} eV can be reached using a complex generation spectrum, but the required CR emissivity is three orders of magnitude higher than that of GRBs, and the predicted spectrum has the GZK cutoff. (iii) The case of local enhancement within region of size 10 - 30 Mpc with overdensity given by factor 3- 30. The overdensity larger than 10 is needed to eliminate the GZK cutoff.

We study the proposal that the cosmic ray primaries above the Greisen-Zatsepin-Kuzmin (GZK) cutoff are gluino-containing hadrons ($\tilde g$-hadrons). We describe the interaction of $\tilde g$-hadrons with nucleons in the framework of the Gribov-Regge approach using a modified version of the hadronic interaction model QGSJET for the generations of Extensive Air Showers (EAS). There are two mass windows marginally allowed for gluinos: $m_{\tilde g}\lsim 3$ GeV and $25\lsim m_{\tilde g}\lsim 35$ GeV. Gluino-containing hadrons corresponding to the second window produce EAS very different from the observed ones. Light $\tilde g$-hadrons corresponding to the first gluino window produce EAS similar to those initiated by protons, and only future detectors can marginally distinguish them. We propose a beam-dump accelerator experiment to search for $\tilde g$-hadrons in this mass window. We emphasize the importance of this experiment: it can discover (or exclude) the light gluino and its role as a cosmic ray primary at ultra high energies.

With SNO data on electron-neutrino flux from the sun, it is possible to derive the $\nu_e$ survival probability $P_{ee}(E)$ from existing experimental data of Super-Kamiokande, gallium experiments and Homestake. The combined data of SNO and Super-Kamiokande provide boron $\nu_e$ flux and the total flux of all active boron neutrinos, giving thus $P_{ee}(E)$ for boron neutrinos. The Homestake detector, after subtraction of the signal from boron neutrinos, gives the flux of Be+CNO neutrinos, and $P_{ee}$ for the corresponding energy interval, if the produced flux is taken from the Standard Solar Model (SSM). Gallium detectors, GALLEX, SAGE and GNO, detect additionally pp-neutrinos. The pp-flux can be calculated subtracting from the gallium signal the rate due to boron, beryllium and CNO neutrinos. The ratio of the measured $pp$-neutrino flux to that predicted by the SSM gives the survival probability for $pp$-neutrinos. Comparison with theoretical survival probabilities shows that the best (among known models) fit is given by LMA and LOW solutions.

Comparison of solar-neutrino signals in SNO and Super-Kamiokande (SK) detectors results in discovery of $\nu_e \to \nu_{\mu,\tau}$ oscillations at level $3.1 - 3.3 \sigma$ \cite{SNO}. This comparison involves the assumption of neutrino spectrum and a choice for the thresholds of detection in both experiments. In this Letter we obtain an exact formula for the comparison of the signals which is valid for arbitrary spectra and thresholds. We find that the no-oscillation hypothesis is excluded at 3.3$\sigma$. If the energy-dependent component of the survival probability for electron neutrinos is small as compared with the average value, i.e. in the case of small distortion of the observed spectrum, the oscillation hypothesis can also be tested to similar accuracy. The oscillation to sterile neutrino only, is excluded at $3.3 \sigma$ level, and oscillation to active neutrinos is confirmed at the same level, though with some reservations.

About 400 years have passed since the great discoveries by Galilei, Kepler and Newton, but astronomy still remains an important source of discoveries in physics. They start with puzzles, with phenomena difficult to explain, and which in fact need for explanation the new physics. Are such puzzles existing now? There are at least three candidates: absence of absorption of TeV gamma radiation in extragalactic space (violation of Lorentz invariance?), absence of GZK cutoff in the spectrum of Ultra High Energy Cosmic Rays (new particle physics?), tremendous energy (up to $10^{54}$ ergs) released in Gamma Ray Bursts during a time scale of a second (collapsing stars or sources of a new type?). Do these puzzles really exist? A critical review of these phenomena is given.

Cusps of superconducting strings can serve as GRB engines. A powerful beamed pulse of electromagnetic radiation from a cusp produces a jet of accelerated particles, whose propagation is terminated by the shock responsible for GRB. A single free parameter, the string scale of symmetry breaking $\eta \sim 10^{14} GeV$, together with reasonable assumptions about the magnitude of cosmic magnetic fields and the fraction of volume that they occupy, explains the GRB rate, duration and fluence, as well as the observed ranges of these quantities. The wiggles on the string can drive the short-time structures of GRB. This model predicts that GRBs are accompanied by strong bursts of gravitational radiation which should be detectable by LIGO, VIRGO and LISA detectors. Another prediction is the diffuse X- and gamma-ray radiation at 8 MeV - 100 GeV with a spectrum and flux comparable to the observed. The weakness of the model is the prediction of too low rate of GRBs from galaxies, as compared with observations. This suggests that either the capture rate of string loops by galaxies is underestimated in our model, or that GRBs from cusps are responsible for only a subset of the observed GRBs not associated with galaxies.

Cusps of superconducting strings can serve as GRB engines. A powerful beamed pulse of electromagnetic radiation from a cusp produces a jet of accelerated particles, whose propagation is terminated by the shock responsible for GRB. A single free parameter, the string scale of symmetry breaking $\eta \sim 10^{14} GeV$, together with reasonable assumptions about the magnitude of cosmic magnetic fields and the fraction of volume that they occupy, explains the GRB rate, duration and fluence, as well as the observed ranges of these quantities. The wiggles on the string can drive the short-time structures of GRB. This model predicts that GRBs are accompanied by strong bursts of gravitational radiation which should be detectable by LIGO, VIRGO and LISA detectors.

We present results from a new Monte Carlo simulation for jet fragmentation in QCD and SUSY QCD for large primary energies $\sqrt s$ up to $10^{16}$ GeV. In the case of SUSY QCD the simulation takes into account not only gluons and quarks as cascading particles, but also their supersymmetric partners. A new model-independent hadronization scheme is developed, in which the hadronization functions are found from LEP data. An interesting feature of SUSY QCD is the prediction of a sizeable flux of the lightest supersymmetric particles (LSPs), if R-parity is conserved. About 10% of the jet energy is transferred to LSPs which, owing to their harder spectra, constitute an important part of the spectra for large $x=E/E_{jet}$. Spectra of protons and of secondary particles, photons and neutrinos, are also calculated. These results have implications for the decay of superheavy particles with masses up to the GUT scale, which have been suggested as a source of ultrahigh energy cosmic rays.

We propose the model of a short-lived very powerful source of high energy neutrinos. It is formed as a result of the dynamical evolution of a galactic nucleus prior to its collapse into a massive black hole and formation of high-luminosity AGN. This stage can be referred to as ``pre-AGN''. A dense central stellar cluster in the galactic nucleus on the late stage of evolution consists of compact stars (neutron stars and stellar mass black holes). This cluster is sunk deep into massive gas envelope produced by destructive collisions of a primary stellar population. Frequent collisions of neutron stars in a central stellar cluster are accompanied by the generation of ultrarelativistic fireballs and shock waves. These repeating fireballs result in a formation of the expanding rarefied cavity inside the envelope. The charged particles are effectively accelerated in the cavity and, due to pp-collisions in the gas envelope, they produce high energy neutrinos. All high energy particles, except neutrinos, are absorbed in the thick envelope. Duration of this pre-AGN phase is about 10 yr, the number of the sources can be \~ 10 per cosmological horizon. High energy neutrino signal can be detected by underground neutrino telescope with effective area ~1 km^2.

The He3 abundance is not constrained by helioseismic data. Mixing of He3 inside the solar core by processes not included in the solar standard model (SSM) has been recently proposed as possible solution to the solar neutrino problem. We have performed a model independent analysis of solar neutrino fluxes using practically arbitrary He3 mixing. In addition, we have been simultaneously varying within very wide ranges the temperature in the neutrino production zone and the astrophysical factors, S17 and S34, of the p+Be7 and He3+He4 cross sections. Seismic data are used as constraints, but the solar- luminosity constraint is not imposed. It is demonstrated that even allowing He3 abundances higher by factors up to 16 than in the SSM, temperatures 5% (or more) lower, the astrophysical factor S17 up to 40% higher, and varying S34 in the range (-20% ,+40%), the best fit is still more than 5 sigmas away from the observed fluxes. We conclude that practically arbitrary He3 mixing combined with independent variations of temperature and cross-sections cannot explain the observed solar neutrino fluxes.

The harmonic analysis of anisotropy of Ultra High Energy Cosmic Rays is performed for the Dark Matter halo model. In this model the relic superheavy particles comprise part of the Dark Matter and are concentrated in the Galactic halo. The Ultra High Energy Cosmic Rays are produced by the decays of these particles. Anisotropy is caused by the non-central position of the Sun in the Galactic halo. The calculated anisotropy is in reasonable agreement with the AGASA data. For more precise test of the model a comparison of fluxes in the directions of the Galactic Center and Anticenter is needed.

The excess of solar-neutrino events above 13 MeV that has been recently observed by Superkamiokande can be explained by vacuum oscillations (VO). If the boron neutrino flux is 20% smaller than the standard solar model (SSM) prediction and the chlorine signal is assumed 30% (or 3.5 sigmas) higher than the measured one, there exists a VO solution that reproduces both the observed boron neutrino spectrum, including the high energy distortion, and the other measured neutrino rates. This solution might already be testable by the predicted anomalous seasonal variation of the gallium signal. Its most distinct signature, a large anomalous seasonal variation of Be7 neutrino flux, can be easily observed by the future detectors, BOREXINO and LENS.

Ultra High Energy Cosmic Rays (UHECR) can be a signal from very early (post-inflationary) Universe. At this cosmological epoch Topological Defects (TD) and long-lived suprheavy (SH) particles are expected to be naturally and effectively produced. Both of these relics can produce now the particles, such as protons and photons, with energies in a great excess of what is observed in UHECR, $E \sim 10^{10} - 10^{11} GeV$. The Topological Defects as the UHECR sources are critically reviewed and cosmic necklaces and monopolonia are identified as most plausible sources. The relic superheavy particles and monopolonia are shown to be clustering in the halo of our Galaxy and their decays produce UHECR without the GZK cutoff. The observational signature of both models are discussed.

The excess of highest energy solar-neutrino events recently observed by Superkamiokande can be in principle explained by anomalously high $hep$-neutrino flux $\Phi_{\nu}(hep)$. Without using SSM calculations, from the solar luminosity constraint we derive that $\Phi_\nu(hep)/S_{13}$ cannot exceed the SSM estimate by more than a factor three. If one makes the additional hypothesis that $hep$ neutrino production occurs where the $^3$He concentration is at equilibrium, helioseismology gives an upper bound which is (less then) two times the SSM prediction. We argue that the anomalous $hep$-neutrino flux of order of that observed by Superkamiokande cannot be explained by astrophysics, but rather by a large production cross-section.

We argue that due to various restrictions cosmic strings and monopole-string networks are not likely to produce the observed flux of ultra-high energy cosmic rays (UHECR). Among the topological defects studied so far, the most promising UHECR sources are necklaces and monopolonia. Other viable sources which are similar to topological defects are relic superheavy particles. All these sources have an excess of pions (and thus photons) over nucleons at production. We demonstrate that in the case of necklaces the diffuse proton flux can be larger than photon flux, due to absorption of the latter on radiobackground, while monopolonia and relic particles are concentrated in the Galactic halo, and the photon flux dominates. Another signature of the latter sources is anisotropy imposed by asymmetric position of the sun in the Galactic halo. In all cases considered so far, including necklaces, photons must be present in ultra-high energy radiation observed from topological defects, and experimental discrimination between photon-induced and proton-induced extensive air showers can give a clue to the origin of ultra-high energy cosmic rays.

The supersymmetric generalization of the limiting and Gaussian QCD spectra is obtained. These spectra are valid for $x \ll 1$, when the main contribution to the parton cascade is given by gluons and gluinos. The derived spectra are applied to decaying superheavy particles with masses up to the GUT scale. These particles can be relics from the Big Bang or produced by topological defects and could give rise to the observed ultrahigh energy cosmic rays. General formulae for the fluxes of protons, photons and neutrinos due to decays of superheavy particles are obtained.

The current status of Ultra High Energy Cosmic Rays (UHECR) is reviewed, with emphasis given to theoretical interpretation of the observed events. The galactic and extragalactic origin, in case of astrophysical sources of UHE particles, have the problems either with acceleration to the observed energies or with the fluxes and spectra. Topological defects can naturally produce particles with energies as observed and much higher, but in most cases fail to produce the observed fluxes. Cosmic necklaces and monopole-antimonopole pairs are identified as most plausible sources, which can provide the observed flux and spectrum. The relic superheavy particles are shown to be clustering in the Galactic halo, producing UHECR without Greisen-Zatsepin-Kuzmin cutoff. The Lightest Supersymmetric Particles are discussed as UHE carriers in the Universe.

Decaying superheavy particles can be produced by Topological Defects or, in case they are quasi-stable, as relics from the early Universe. The decays of these particles can be the sources of observed Ultra High Energy Cosmic Rays ($E \sim 10^{10} - 10^{12} GeV$). The Topological Defects as the UHE CR sources are critically reviewed and cosmic necklaces and monopole-antiminopole pairs are identified as most plausible sources. The relic superheavy particles are shown to be clustering in the halo and their decays produce UHE CR without GZK cutoff. The Lightest Supersymmetric Particles with Ultra High Energies are naturally produced in the cascades accompanying the decays of superheavy particles. These particles are discussed as UHE carriers in the Universe.

Solar neutrinos are discussed in the light of the new data and of recent progress in helioseismology. Most attention is given to the new status of Standard Solar Models due to seismically measured density and sound speed in the inner solar core. The elementary particle solutions to the Solar Neutrino Problem and their observational signatures are discussed.

We argue that the lightest supersymmetric particles (LSP) can be produced with extremely high energies $E\gsim 10^{10}$ GeV in the Universe at the present epoch. Their most probable sources are decaying superheavy particles produced by topological defects or as relic Big Bang particle. We discuss the mechanisms of production of LSP at ultra-high energies (UHE) and the interaction of the UHE LSP with matter. The most attention is given to the neutralino as LSP, although the gluino is also considered as a phenomenological possibility.

We study the decays of ultraheavy (m_X \geq 10^{13}GeV) and quasistable (lifetime \tau_X much larger than the age of the Universe t_0) particles as the source of Ultra High Energy Cosmic Rays (UHE CR). These particles are assumed to constitute a tiny fraction \xi_X of CDM in the Universe, with \xi_X being the same in the halo of our Galaxy and in the intergalactic space. The elementary-particle and cosmological scenarios for these particles are briefly outlined. The UHE CR fluxes produced at the decays of X- particles are calculated. The dominant contribution is given by fluxes of photons and nucleons from the halo of our Galaxy and thus they do not exibit the GZK cutoff. The extragalactic components of UHE CR are suppressed by the smaller extragalactic density of X-particles and hence the cascade limit is relaxed. We discuss the spectrum of produced Extensive Air Showers (EAS) and a signal from Virgo cluster as signatures of this model.

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.

The central solar temperature T and its uncertainties are calculated in helioseismologically-constrained solar models. From the best fit to the convective radius, density at the convective radius and seismically determined helium abundance the central temperature is found to be T=1.58x10^7 K, in excellent agreement with Standard Solar Models. Conservatively, we estimate that the accuracy of this determination is Delta T/T=1.4 %, better than that in SSM. Neutrino fluxes are calculated. The lower limit to the boron neutrino flux, obtained with maximum reduction factors from all sources of uncertainties, is 2 sigma higher than the flux measured recently by SuperKamiokande.

The recent observations of 50 TeV gamma radiation by HEGRA have the potential of determining the extragalactic flux of infrared radiation. The fact that radiation is observed in the range between 30 and 100 TeV sets an upper limit on the infrared flux, while a cutoff at $E_{\gamma} \approx 50$ TeV fixes this flux with a good accuracy. If the intrinsic radiation is produced due to interaction of high energy protons with gas or low-energy target photons, then an accompaning high-energy neutrino flux is unavoidable. We calculate this flux and underground muon flux produced by it. The muon flux is dominated by muons with energies about 1 TeV and can be marginally detected by a 1 km$^2$ detector like an expanded AMANDA.

Cosmic necklaces are hybrid topological defects consisting of monopoles and strings, with two strings attached to each monopole. We argue that the cosmological evolution of necklaces may significantly differ from that of cosmic strings. The typical velocity of necklaces can be much smaller than the speed of light, and the characteristic scale of the network much smaller than the horizon. We estimate the flux of high-energy protons produced by monopole annihilation in the decaying closed loops. For some reasonable values of the parameters it is comparable to the observed flux of ultrahigh-energy cosmic rays.

Monopole-antimonopole pairs connected by strings and monopole-string networks with $N>2$ strings attached to each monopole can be formed at phase transitions in the early universe. In such hybrid defects, monopoles accelerate under the string tension and can reach ultrarelativistic Lorentz factors, $\gamma\gg 1$. We study the radiation of gauge quanta by accelerating monopoles. For monopoles with a chromomagnetic charge, we also discuss the high-energy hadron production through emission of virtual gluons and their subsequent fragmentation into hadrons. The relevant parameter for gauge boson radiation is $M/a$, where $M$ is the boson mass and $a$ is the proper acceleration of the monopole. For $M\ll a$, the gauge bosons can be considered as massless and the typical energy of the emitted quanta is $E\sim\gamma a$. In the opposite limit, $M\gg a$, the radiation power is exponentially suppressed and gauge quanta are emitted with a typical energy $E\sim\gamma M$ in a narrow range $\Delta E/E\sim (a/M)^{1/2}$. Cosmological monopole-string networks can produce photons and hadrons of extremely high energies. For a wide range of parameters these energies can be much greater than the Planck scale.

These are the Proceedings of the 4th discussion workshop in the Gran Sasso Laboratory. The structure of this workshop was as before: most of the time was allocated to the review talks and to the discussions. There were discussion sessions and short original presentations for each topic as well as the general discussion session at the end of the workshop.

The baryonic and cold dark matter are reviewed in the context of cosmological models. The theoretical search for the particle candidates is limited by supersymmetric extension of the Standard Model. Generically in such models there are just two candidates associated with each other: generalized neutralino, which components are usual neutralino and axino, and axion which is a partner of axino in supermultiplet. The status of these particles as DM candidates is described.

The microlensing objects, Machos, recently observed in the halo of our Galaxy, can be interpreted as dense neutralino objects, neutralino stars, produced by the gravitational instability of neutralino gas. Taking the mass and radius of these objects from microlensing observations we calculated the diffuse gamma-ray flux produced in neutralino--neutralino annihilation inside the objects. The resulting flux is many orders of magnitude higher than the observed one.

Neutrino telescopes of large area offer the possibility of searching for indirect signals of relic neutralinos in the galactic halo, due to annihilations in the Sun or the Earth. Here we investigate the sensitivity, using a supergravity scheme where the soft scalar mass terms are not constrained by universality conditions at the grand unification scale. We first discuss in which regions of the supersymmetric parameter space the neutralino may be considered as a good candidate for cold dark matter. The discovery potential of the search using neutrino telescopes is then compared to that of the direct search for relic neutralinos.

The best particle candidates for non--baryonic cold dark matter are reviewed, namely, neutralino, axion, axino and Majoron. These particles are considered in the context of cosmological models with the restrictions given by the observed mass spectrum of large scale structures, data on clusters of galaxies, age of the Universe etc.

We propose a new mechanism for high energy gamma ray burst in Supernova type I (SNI) explosion. Presupernova is assumed to be a binary system comprised of a red giant and a white dwarf with a wind accretion. The accretion flow terminates by the accretion shock in the vicinity of the white dwarf at the distance of the order of the accretion radius. The gas inside the accretion radius constitutes the main body of a target for gamma ray production by $pp$ collisions.\par The supernova explosion and the shock propagation in the white dwarf result in the hydrodynamical acceleration of the outer layers of the star. The spectrum of accelerated particles is steep and the maximum energy does not exceed $1000$ $GeV$.\par The gamma ray burst is produced due to interaction of the accelerated particles with the gas in the binary system. Most of the photons have energies about $100$ $MeV$, the total number of the emitted photons is between $10^{46}$ and $10^{47}$. The gamma ray burst might have one or two precursors. The first one is produced during the shock breakout when the shock approaches the star surface and crosses it.\par The calculations are valid in the case of a weak magnetic field. The case of strong magnetic field will be considered in part II (in preparation).

We discuss what appears the last hope for an astrophysical solution to the solar neutrino problem: a correlated variation of the astrophysical factors for the helium burning cross sections ($S_{33}$ and $S_{34}$) and either $S_{17}$ or the central temperature $T_c$. In this context, we recognize the important role played by the CNO neutrinos. In fact, we can obtain a fair fit to the experimental data only if three conditions are met simultaneously: the astrophysical factor $S_{33}$ is about 200 times what is presently estimated, the astrophysical factor $S_{17}$ is about 3 times larger and the $^{13}$N and $^{15}$O neutrino fluxes are negligible compared to the ones predicted by standard solar models. These conditions are not supported by the present data and their correlated combination is improbable.

Neutralino dark matter is studied in the context of a supergravity scheme where the scalar mass terms are not constrained by universality conditions at the grand unification scale. We analyse in detail the consequences of the relaxation of this universality assumption on the supersymmetric parameter space, on the neutralino relic abundance and on the event rate for the direct detection of relic neutralinos.

We study pion curvature radiation by a proton, i.e. pion emission by a proton moving along a curved trajectory. We suggest an approximate semiclassical solution and the exact solution for which we assume that a proton moves in a fictitious magnetic field with the Larmor radius equal to the curvature radius of the real trajectory. As possible application we consider the pion radiation by ultrahigh energy protons moving along curved magnetic field lines. Such situation can occur in the magnetosphere of a young pulsar, in the magnetosphere of the accretion disk around a black hole, and in the vicinity of a superconducting cosmic string. The decay products of these pions, such as high energy photons or neutrinos, can give the observable consequences of the considered mechanism.

We consider the production of $^3$He and $^2$H by $^4$He photodisintegration initiated by non-thermal energy releases during early cosmic epochs. We find that this process cannot be the predominant source of primordial $^2$H since it would result in anomalously high $^3$He/D ratios in conflict with standard chemical evolution assumptions. We apply this fact to constrain topological defect models of highest energy cosmic ray (HECR) production. Such models have been proposed as possible sources of ultrahigh energy particles and gamma-rays with energies above $10^{20}$eV. The constraints on these models derived from $^4$He-photodisintegration are compared to corresponding limits from spectral distortions of the cosmic microwave background radiation (CMBR) and from the observed diffuse gamma-ray background. It is shown that for reasonable primary particle injection spectra superconducting cosmic strings, unlike ordinary strings or annihilating monopoles, cannot produce the HECR flux at the present epoch without violating at least the $^4$He-photodisintegration bound. The constraint from the diffuse gamma-ray background rules out the dominant production of HECR by the decay of Grand Unification particles in models with cosmological evolution assuming standard fragmentation functions. Constraints on massive black hole induced photodisintegration are also discussed.

We describe a calculation of electromagnetic cascading in radiation and matter in the early universe initiated by the decay of massive particles or by some other process. We have used a combination of Monte Carlo and numerical techniques which enables us to use exact cross sections, where known, for all the relevant processes. In cascades initiated after the epoch of big bang nucleosynthesis $\gamma$-rays in the cascades will photodisintegrate $^4$He, producing $^3$He and deuterium. Using the observed $^3$He and deuterium abundances we are able to place constraints on the cascade energy deposition as a function of cosmic time. In the case of the decay of massive primordial particles, we place limits on the density of massive primordial particles as a function of their mean decay time, and on the expected intensity of decay neutrinos.

The difficulties for non-standard solar models (NSSM) in resolving the solar neutrino problem are discussed stressing the incompatibility of the gallium--Kamiokande data, and of the gallium--chlorine data. We conclude that NSSM's cannot explain simultaneously the results of any two of the solar neutrino data (chlorine, Kamiokande and gallium). We address further the question whether the MSW solution exists for NSSM's (e.g. models with $^8$B neutrino flux much lower than the standard one and/or central temperature $T_c$ very different from $T_c^{\text{SSM}}$). We demonstrate that the MSW solution exists and is very stable relative to changes of $S_{17}$ ($S$-factor for $p$ + $^7$Be reaction) and $T_c$. In particular, $\Delta m^2$ is almost constant while $\sin^2 2\theta$ depends on the exact values of $S_{17}$ (or $^8$B-neutrino flux) and $T_c$.