The fate of ultrahigh energy nuclei in the immediate environment of young fast-rotating pulsars

Young, fast-rotating neutron stars are promising candidate sources for the production of ultrahigh energy cosmic rays (UHECRs). The interest in this model has recently been boosted by the latest chemical composition measurements of cosmic rays, that seem to show the presence of a heavy nuclear component at the highest energies. Neutrons stars, with their metal-rich surfaces, are potentially interesting sources of such nuclei, but some open issues remain: 1) is it possible to extract these nuclei from the stars surface? 2) Do the nuclei survive the severe conditions present in the magnetosphere of the neutron star? 3) What happens to the surviving nuclei once they enter the wind that is launched outside the light cylinder? In this paper we address these issues in a quantitative way, proving that for the most reasonable range of neutron star surface temperatures ($T<10^7,$K), a large fraction of heavy nuclei survive photo-disintegration losses. These processes, together with curvature losses and acceleration in the stars electric potential, lead to injection of nuclei with a chemical composition that is mixed, even if only iron is extracted from the surface. We show that under certain conditions the chemical composition injected into the wind region is compatible with that required in previous work based on purely phenomenological arguments (typically $sim 50%$ protons, $sim 30%$ CNO and $sim 20%$ Fe), and provides a reasonable explanation of the mass abundance inferred from ultra high energy data.

On the cosmic ray spectrum from type II Supernovae expanding in their red giant presupernova wind

While from the energetic point of view SNRs are viable sources of Galactic CRs, the issue of whether they can accelerate protons up to PeV remains unsolved. Here we discuss particle acceleration at the forward shock of SN and discuss the possibility that the escaping particle current may excite a non-resonant instability that in turn leads to the formation of resonant modes confining particles close to the shock and increasing the maximum energy. This mechanism works throughout the expansion of the SN explosion, from the ejecta dominated (ED) to the Sedov-Taylor (ST) phase. Because of their higher explosion rate,we focus on type II SNae expanding in the slow, dense red supergiant wind. When the explosion occurs in such winds, the transition between the ED and the ST phase is likely to take place within a few tens of years. As a result, the spectrum of accelerated particles shows a break in the slope, at the maximum energy (Em) achieved at the beginning of the ST phase. Above this energy, the spectrum becomes steeper but remains a power law than developing an exponential cutoff. We show that for type II SNae typical parameters, proton Em can easily reach PeV energies, confirming that type II SNRs are the best candidate sources for CRs at the knee. We have tried to fit KASCADE-Grande, ARGO -YBJ and YAC1-Tibet Array data with our model but we could not find any parameter combination that could explain all data sets. Indeed the recent measurement of the proton and helium spectra in the knee region, with the ARGO-YBJ and YAC1-Tibet Array, has made the situation very confused. These measurements suggest that the knee in the light component is at 650 TeV, appreciably below the overall spectrum knee. This finding would resolve the problem of reaching very high energies in SNae, but, on the other hand, it would open a critical issue in the transition region between Galactic and extragalactic CRs.

Spectral breaks as a signature of cosmic ray induced turbulence in the Galaxy

We show that the complex shape of the cosmic ray (CR) spectrum, as recently measured by PAMELA and inferred from Fermi-LAT gamma-ray observations of molecular clouds in the Gould belt, can be naturally understood in terms of basic plasma astrophysics phenomena. A break from a harder to a softer spectrum at blue rigidity Rsimeq 10 GV follows from a transition from transport dominated by advection of particles with Alfven waves to a regime where diffusion in the turbulence generated by the same CRs is dominant. A second break at Rsimeq 200 GV happens when the diffusive propagation is no longer determined by the self-generated turbulence, but rather by the cascading of externally generated turbulence (for instance due to supernova (SN) bubbles) from large spatial scales to smaller scales where CRs can resonate. Implications of this scenario for the cosmic ray spectrum, grammage and anisotropy are discussed.

A kinetic approach to cosmic ray induced streaming instability at supernova shocks

We show that a purely kinetic approach to the excitation of waves by cosmic rays in the vicinity of a shock front leads to predict the appearance of a non-alfvenic fastly growing mode which has the same dispersion relation as that previously found by cite{bell04} by treating the plasma in the MHD approximation. The kinetic approach allows us to investigate the dependence of the dispersion relation of these waves on the microphysics of the current which compensates the cosmic ray flow. We also show that a resonant and a non-resonant mode may appear at the same time and one of the two may become dominant on the other depending on the conditions in the acceleration region. We discuss the role of the unstable modes for magnetic field amplification and particle acceleration in supernova remnants at different stages of the remnant evolution.

Simulated synchrotron and Inverse Compton emission from Pulsar Wind Nebulae

We present a complete set of diagnostic tools aimed at reproducing synthetic non-thermal (synchrotron and/or Inverse Compton, IC) emissivity, integrated flux energy, polarization and spectral index simulated maps in comparison to observations. The time dependent relativistic magnetohydrodynamic (RMHD) equations are solved with a shock capturing code together with the evolution of the maximum particles energy. Applications to Pulsar Wind Nebulae (PWNe) are shown.