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.

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.

Present experimental data do not exclude fast oscillation of the neutron $n$ to its degenerate twin from a hypothetical parallel sector, the so called mirror neutron $n'$. We show that this effect brings to a remarkable modifications of the ultrahigh-energy cosmic ray spectrum testable by the present Pierre Auger Observatory (PAO) and Telescope Array (TA) detector, and the future JEM-EUSO experiment. In particular, the baryon non-conservation during UHECR propagation at large cosmological distances shifts the beginning of the GZK cutoff to lower energies, while in presence of mirror sources it may enhance the spectrum at $E \gtrsim 100$ EeV. As a consequence, a significant reduction of the expected diffuse cosmogenic neutrino flux is predicted.

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.

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.

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.

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.

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

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 present the high-energy neutrino Monte Carlo event generator ANIS (All Neutrino Interaction Simulation). The program provides a detailed and flexible neutrino event simulation for high-energy neutrino detectors, such as AMANDA, ANTARES or ICECUBE. It generates neutrinos of any flavor according to a specified flux and propagates them through the Earth. In a final step neutrino interactions are simulated within a specified volume. All relevant standard model processes are implemented. We discuss strengths and limitations of the program.

We present the high energy neutrino Monte Carlo event generator ANIS (All Neutrino Interaction Simulation). The aim of the program is to provide a detailed and flexible neutrino event simulation for high energy neutrino detectors, such as AMANDA and ICECUBE. It generates neutrinos of any flavor according to a specific flux, propagates them through the Earth and in a final step simulates neutrino interactions within a specified volume. All relevant standard model processes are implemented. We discuss strength and limitations of the program, and provide as an example event rates for atmospheric and E^-2 neutrino spectra.

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.

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.

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.

Various implications of new, non-perturbative pomeron inspired enhancement of small-x neutrino-nucleon structure functions for high-energy neutrino astrophysics are discussed. At x larger than 10^{-5} these functions are given by perturbative QCD, while at lower x they are determined by a specific generalization of F_2^{ep}(x,Q^2) description, proposed by A. Donnachie and P. V. Landshoff (their two-component model comprises hard and soft pomerons), to neutrino-nucleon scattering case. We found that i) such enhancement causes the most rapid growth of neutrino-nucleon cross-sections at high energies, ii) pomeron effects may be perceptible in the rates of neutrino induced events in future giant detectors and iii) the rate of high-energy neutrino flux evolution (due to absorption (CC+NC) and regeneration (NC)) on its pass through a large column depth of matter may be subjected to additional influence of hard pomeron. Solving transport equations for the initially power-law decreasing neutrino spectra, we have evaluated shadow factors for several column depths and spectrum indices. The results are compared with analogous calculations, performed within a trivial small-x extrapolation of structure functions. Hard pomeron enhanced high-energy shadow factors are found to be many orders of magnitude lower than those obtained within ordinary perturbative QCD.

Rapid growth of neutrino-nucleon cross-sections at high energies due to hypothetical hard pomeron enhancement of $\nu N$-structure functions is discussed. Differential and integral hadron moments, which, together with cross-sections, define rates of hadron-electromagnetic cascades in a neutrino detector are calculated for different power-law decreasing neutrino spectra. For comparison two small-x extrapolation schemes are discussed. First includes Regge theory inspired hard pomeron enhancement of $\nu N$-structure functions. The second is obtained with the help of trivial extrapolation of perturbative QCD structure functions from the large x region to the small x one. Implications of hard pomeron effects for cross-sections and hadron moments are demonstrated. The most pronounced manifestations are found in integral hadron moments for the case of charged current electron (anti)neutrinos scattering off nucleons.

An unknown small-x behavior of nucleon structure functions gives appreciable uncertainties to high-energy neutrino-nucleon cross-sections. We construct structure functions using at small x Regge inspired description by A. Donnachie and P. V. Landshoff with soft and hard pomerons, and employing at larger x the perturbative QCD expressions. The smooth interpolation between two regimes for each Q^2 is provided with the help of simple polynomial functions. To obtain low-x neutrino-nucleon structure functions $F_2^{\nu N, \bar \nu N}(x,Q^2)$ and singlet part of $F_{3}^{\nu N,\bar \nu N}(x,Q^2)$ from Donnachie-Landshoff function $F_2^{ep}(x,Q^2)$, we use the Q^2-dependent ratios R_2(Q^2) and R_3(Q^2) derived from perturbative QCD calculations. Non-singlet part of F_3 at low x, which is very small, is taken as power-law extrapolation of perturbative function at larger x. This procedure gives a full set of smooth neutrino-nucleon structure functions in the whole range of x and Q^2 at interest. Using these structure functions, we have calculated the neutrino-nucleon cross-sections and compared them with some other cross-sections known in literature. Our cross-sections turn out to be the highest among them at the highest energies, which is explained by contribution of the hard pomeron.