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

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

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