In this thesis, I describe in great detail the physics of the decay of any Super-Heavy X particle (with masses up to the grand unification scale ~ 10^16 GeV and possibly beyond), and the computer code I developed to model this process - which currently is the most complete available one. The general framework for this work is the Minimal Supersymmetric Standard Model (MSSM). The results are presented in the form of fragmentation functions of any (s)particle of the MSSM into any final stable particle (proton, photon, electron, three types of neutrino, lightest superparticle LSP) at a virtuality Q = M_X, over a scaled energy range x = 2E/M_X in [10^{-13}, 1]. At very low x values, color coherence effects have been taken into account through the Modified Leading Log Approximation (MLLA). The whole process is explicitely shown to conserve energy with a numerical accuracy up to a few part per mille, which allows to make quantitative predictions for any N-body decay mode of any X particle. I then apply the results to the old - and yet unsolved - problem of Ultra High Energy Cosmic Rays (UHECRs). In particular, I provide quantitative predictions of generic ``top-down'' models for the neutrino and neutralino fluxes which could be observed in the next generation of detectors.

Decays of superheavy X particles with mass M_X ~ 10^12 - 10^16 GeV have been proposed as origin of the observed ultra high energy cosmic rays (UHECR). We describe in detail the physics involved in the different steps of the decay of such a particle. In particular, we give for the first time the complete set of splitting functions needed to model a parton shower in the minimal supersymmetric extension of the Standard Model (MSSM). We present our results in the form of fragmentation functions of any (s)particle of the MSSM into any final stable particle (proton, photon, electron, three types of neutrino, lightest superparticle LSP) at a virtuality Q = M_X, over a scaled energy range x = 2E/M_X in [10^{-13}, 1]. Extending the coverage to such small fractional energies is necessary since the energy region around 10^18 eV and below could be of considerable interest in testing this kind of model for generating UHECR. We explicitly demonstrate that our treatment conserves energy, and discuss the dependence of the final result on SUSY parameters. We also show that our results are essentially independent of the necessary extrapolation of the input fragmentation functions, which are known only for x >= 0.1, towards small x. Finally, we added a new treatment of the color coherence effects at very small x, using the analytic ``MLLA'' solution. Our computer code will soon be made available.

Top-down models of cosmic rays produce more neutrinos than photons and more photons than protons. In these models, we reevaluate the fluxes of neutrinos associated with the highest energy cosmic rays in light of mounting evidence that they are protons and not gamma rays. While proton dominance at EeV energies can possibly be achieved by efficient absorption of the dominant high-energy photon flux on universal and galactic photon and magnetic background fields, we show that the associated neutrino flux is inevitably increased to a level where it should be within reach of operating experiments such as AMANDA II, RICE and AGASA. In future neutrino telescopes, tens to a hundred, rather than a few neutrinos per kilometer squared per year, may be detected above 1 PeV.

We present a new and complete numerical analysis of the decay of super-heavy X particles, assumed to be the origin of cosmic rays with energy beyond the GZK cut-off. The decay of X initiates a ``parton shower'', where we include all degrees of freedom contained in the Minimal Supersymmetric Standard Model (MSSM). Since at energies near $M_X$ all gauge couplings are of similar magnitude, we include all of them, as well as third generation Yukawa couplings. Technically the shower development is described through the DGLAP evolution of the relevant fragmentation functions (FFs). We also carefully treat the decay of the superparticles as well as heavy SM particles created in the shower. Nonperturbative physics is parameterized through the input values of the FFs, which we take from the literature. The final result is the the complete spectrum of all stable particles at the very end of the shower : protons, electrons, neutrinos, photons and neutralinos, for an energy range from $10^{-7} M_X$ to $M_X$. In particular, the flux of high-energy neutralinos is sizable ; it might serve as ``smoking gun'' signature for this kind of scenario.