Parameterization of Gamma, e^+/- and Neutrino Spectra Produced by p-p   Interaction in Astronomical Environment
  We present the yield and spectra of stable secondary particles (gamma, e^+/-, nu_e, nubar_e, nu_mu, and nubar_mu) of p-p interaction in parameterized formulae to facilitate calculations involving them in astronomical environments. The formulae are derived on the up-to-date p-p interaction model by [Kamae05] which incorporates the logarithmically rising inelastic cross-section, the diffraction dissociation process, and the Feynman scaling violation. To improve fidelity to experimental data in lower energies, two baryon resonance contributions have been added: one representing Delta(1232) and the other multiple resonances around 1600MeV/c^2. The parametrized formulae predict that all secondary particle spectra be harder by about 0.05 in power-law indices than that of the incident proton and their inclusive cross-sections be larger than those predicted by p-p interaction models based on the Feynman scaling. 
Tuneyoshi Kamae, Niklas Karlsson, Tsunefumi Mizuno, Toshinori Abe and   Tatsumi Koi