We examine in detail the detectability of a signal of diffuse high energy neutrinos produced in the Milky Way by the interaction of cosmic rays (CRs) with the interstellar medium (ISM). We show that highly inhomogeneous galactic CR densities arise naturally from large scale drift effects and the inhomogeneous distribution of CR sources. In particular, the CR density in the Galactic Center region (where the ISM density is maximal) could be an order of magnitude larger than the local CR density. Hence, the expected diffuse flux of neutrinos from the Galactic Center region becomes enhanced and could be detected at energies above $\sim 10^4$ GeV in a km$^3$-size neutrino telescope within 10 years of operation. The galactic anisotropy, which is the main signature of this flux, allows to discriminate this signal from other possible extraterrestrial diffuse fluxes of high energy neutrinos.

We explore several aspects related to the propagation of high energy cosmic rays (CRs) of galactic origin (i.e. in the range $E\simeq 10^{15}- 10^{18}$ eV). In particular, we study in detail the diffusion/drift scenario, a rigidity-dependent model that explains the origin of the knee and the second knee of the CR spectrum, as well as the most recent observations of composition and anisotropies. Moreover, the behavior of the diffusion tensor for charged particles propagating in highly turbulent magnetic fields is determined numerically. We also investigate the high energy neutrino fluxes induced by CRs, and we show that a rigidity-dependent scenario for the CR knee reduces the background of atmospheric neutrinos, hence encouraging the search for astrophysical neutrino sources. After exploring this effect in the detection channel of lepton tracks, we investigate the less studied ``shower'' channel, which turns out to be particularly effective for isolating the prompt atmospheric neutrino component. Finally, we consider in detail the nuclear photodisintegration model as an alternative explanation for the knee in the CR spectrum, in which CR nuclei emitted by a source are disintegrated during their propagation through a surrounding optical/UV photon field. In this case, the composition is predicted to turn lighter across the knee, in constrast to the results obtained in rigidity-dependent scenarios.

At high energies, the very steep decrease of the conventional atmospheric component of the neutrino spectrum should allow the emergence of even small and isotropic components of the total spectrum, indicative of new physics, provided that they are less steeply decreasing, as generically expected. One candidate is the prompt atmospheric neutrino flux, a probe of cosmic ray composition in the region of the knee as well as small-$x$ QCD, below the reach of collider experiments. A second is the diffuse extragalactic background due to distant and unresolved AGNs and GRBs, a key test of the nature of the highest-energy sources in the universe. Separating these new physics components from the conventional atmospheric neutrino flux, as well as from each other, will be very challenging. We show that the charged-current {\it electron} neutrino "shower" channel should be particularly effective for isolating the prompt atmospheric neutrino flux, and that it is more generally an important complement to the usually-considered charged-current {\it muon} neutrino "track" channel. These conclusions remain true even for the low prompt atmospheric neutrino flux predicted in a realistic cosmic ray scenario with heavy and varying composition across the knee (Candia and Roulet, 2003 JCAP {\bf 0309}, 005). We also improve the corresponding calculation of the neutrino flux induced by cosmic ray collisions with the interstellar medium.

We determine numerically the parallel, perpendicular, and antisymmetric diffusion coefficients for charged particles propagating in highly turbulent magnetic fields, by means of extensive Monte Carlo simulations. We propose simple expressions, given in terms of a small set of fitting parameters, to account for the diffusion coefficients as functions of magnetic rigidity and turbulence level, and corresponding to different kinds of turbulence spectra. The results obtained satisfy scaling relations, which make them useful for describing the cosmic ray origin and transport in a variety of different astrophysical environments.

Scenarios in which the knee of the cosmic ray spectrum depends on the particle rigidities usually predict that the cosmic ray composition becomes heavier above the knee and have associated a change in the spectral slope of each individual nuclear component which is steeper than the change ($\Delta\alpha\simeq 0.3$) observed in the total spectrum. We show that this implies that the very high energy ($E_\nu>10^{14}$ eV) diffuse neutrino fluxes produced by cosmic rays hitting the atmosphere or colliding with the interstellar medium in the Galaxy will be significantly suppressed, making their detection harder but also reducing the background for the search of other (more challenging) astrophysical neutrino sources.

We consider the scenario in which the knee in the cosmic ray spectrum is due to a change in the escape mechanism of cosmic rays from the Galaxy from one dominated by transverse diffusion to one dominated by drifts. We show that this scenario explains not only the changes in spectral slope at the knee and at the second knee, but can also account for the main characteristics of the observed energy dependent anisotropy amplitude and phase of first harmonic in the energy range between $10^{15}$ and $10^{18}$ eV. This provides a useful handle to distinguish this diffusion/drift model from other scenarios proposed to explain the knee in the spectrum.

We reconsider the scenario in which the knee in the cosmic ray spectrum is explained as due to a change in the escape mechanism of cosmic rays from the Galaxy from one dominated by transverse diffusion to one dominated by drifts. We solve the diffusion equations adopting realistic galactic field models and using diffusion coefficients appropriate for strong turbulence (with a Kolmogorov spectrum of fluctuations) and consistent with the assumed magnetic fields. We show that properly taking into account these effects leads to a natural explanation of the knee in the spectrum, and a transition towards a heavier composition above the knee is predicted.

We show that the second knee in the cosmic ray spectrum (i.e. the steepening occurring at $E\simeq 4\times 10^{17}$ eV) could be related to drift effects affecting the heaviest nuclear component, the iron group nuclei, in a scenario in which the knee at $3\times 10^{15}$ eV indicates the onset of drift effects in the lighter proton component. We also study the anisotropies resulting from diffusion and drift currents in the Galaxy, showing their potential relevance to account for the AGASA observations at $E\sim 10^{18}$ eV, before the extragalactic component becomes dominant.

We explore in some detail the scenario proposed to explain the observed knee of the cosmic ray (CR) spectrum as due to the effects of photodisintegration of the CR nuclei by interactions with optical and soft UV photons in the source region. We show that the photon column densities needed to explain the experimental data are significantly lower than those obtained in previous estimations which neglected multinucleon emission in the photodisintegration process. We also treat more accurately the photodisintegration thresholds, we discuss the effects of photopion production processes and the neutron escape mechanism, identifying the physical processes responsible for the qualitative features of the results. This scenario would require the CR nuclei to traverse column densities of $\sim 5 \times 10^{27}- 2 \times 10^{28}$ eV/cm$^2$ after being accelerated in order to reproduce the observed knee, and predicts that the CR composition should become lighter above $\sim 10^{16}$ eV.