We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.

Core-collapse supernovae produce fast shocks which expand into the dense circumstellar medium (CSM) of the stellar progenitor. Cosmic rays (CRs) accelerated at these shocks can induce the growth of electromagnetic fluctuations in the pre-shock medium. Using a self-similar description for the shock evolution, we calculate the growth time-scales of CR driven instabilities for SNe in general, and SN 1993J in particular. We find that extended SN shocks can trigger fast intra-day instabilities, strong magnetic field amplification, and CR acceleration. In particular, the non-resonant streaming instability can contribute to about 50 per cent of the magnetic field intensity deduced from radio data. This results in the acceleration of CR particles to energies of 1-10 PeV within a few days after the shock breakout.

We study electron and positron acceleration at the termination shock of a striped pulsar wind. Drift motion along the shock surface keeps either electrons or positrons -but not both, close to the equatorial plane of the pulsar, where they are accelerated by the first-order Fermi process. Their energy spectrum is a power law, and both the X-ray flux and photon index of the Crab Nebula, as measured by NuSTAR, can be reproduced for sufficiently large downstream turbulence levels. The implication that one sign of charge is preferentially accelerated in pulsar wind nebulae is potentially important for the interpretation of the positron fraction in cosmic-rays.

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

Two-fluid (electron-positron) plasma modelling has shown that inductive acceleration can convert Poynting flux directly into bulk kinetic energy in the relativistic flows driven by rotating magnetized neutron stars and black holes. Here, we generalize this approach by adding an ion fluid. Solutions are presented in which all particles are accelerated as the flow expands, with comparable power channeled into each of the plasma components. In an ion-dominated flow, each species reaches the limiting rigidity, according to Hillas' criterion, in a distance significantly shorter than in a lepton-dominated flow. These solutions support the hypothesis that newly born magnetars and pulsars are potential sources of ultra-high energy cosmic rays. The competing process of Poynting flux dissipation by magnetic reconnection is shown to be ineffective in low-density flows in which the conventionally defined electron multiplicity satisfies $\kappa_{\rm e}\lesssim 10^5\left(4\pi L_{38}/\Omega\right)^{1/4} /\textrm{Max}\left(\eta_{\rm ion}^{1/2},1\right)$, where $L_{38}\times 10^{38}\textrm{erg s}^{-1}$ is the power carried by the flow in a solid angle $\Omega$, and $\eta_{\rm ion}$ is the ratio of the ion to lepton power at launch.

Core-collapse supernovae exploding in dense winds are favorable sites for cosmic-ray (CR) acceleration to very high energies. We present our CR-radiation-hydrodynamics simulations of the explosion of a red supergiant. We study the evolution of the shock wave during the first day following core collapse, and estimate the time at which CR acceleration can start. We then calculate the maximum CR energy at the forward shock as a function of time, and show that it may already exceed 100 TeV only a few hours after shock breakout from the surface of the star.

The discovery of extended TeV emission around the Geminga and PSR B0656+14 pulsars, with properties consistent with free particle propagation in the interstellar medium (ISM), has sparked considerable discussion on the possible presence of such halos in other systems. Here we make an assessment of the current TeV source population associated with energetic pulsars, in terms of size and estimated energy density. Based on two alternative estimators we conclude that a large majority of the known TeV sources have emission originating in the zone energetically and dynamically dominated by the pulsar (i.e. the pulsar wind nebula), rather than from a halo of particles diffusing in to the ISM. Furthermore, whilst the number of established halos will surely increase in the future, we find that it is unlikely that such halos contribute significantly to the total TeV $\gamma$-ray luminosity from electrons accelerated in PWN.

Cosmic ray electrons and positrons are tracers of particle propagation in the interstellar medium (ISM). A recent measurement performed using H.E.S.S. extends the all-electron (electron+positron) spectrum up to 20TeV, probing very local sources and transport due to the $\sim$10~kyr cooling time of these particles. An additional key local measurement was the recent estimation of the ISM diffusion coefficient around Geminga performed using HAWC. The inferred diffusion coefficient is much lower than typically assumed values. It has been argued that if this diffusion coefficient is representative of the local ISM, pulsars would not be able to account for the all-electron spectrum measured at the Earth. Here we show that a low diffusion coefficient in the local ISM is compatible with a pulsar wind nebula origin of the highest energy electrons, if a so far undiscovered pulsar with spin-down power $\sim 10^{33-34}$ erg/s exists within 30 to 80~pc of the Earth. The existence of such a pulsar is broadly consistent with the known population and may be detected in near future survey observations.

We place constraints on the properties of the interstellar turbulence that surrounds Geminga pulsar, using the recent measurements from the HAWC Observatory in this region. We propagate very-high-energy electrons in realizations of 3D isotropic Kolmogorov or Kraichnan turbulence, calculate their gamma-ray emission, and compare with HAWC measurements. We show that the measurements can be well fitted for both models of the turbulence and for reasonable values of its strength, $B_{\rm rms}$, and coherence length, $L_{\rm c}$. Our best fits are obtained for $B_{\rm rms} \simeq 3 \mu$G and $L_{\rm c} \simeq 1$ pc. Furthermore, the absence of strong asymmetries in the observed emission favours $L_{\rm c} \leq 5$ pc.

We calculate the shape of the large-scale anisotropy of TeV-PeV cosmic-rays (CR) in different models of the interstellar turbulence. In general, the large-scale CR anisotropy (CRA) is not a dipole, and its shape can be used as a new probe of the turbulence. The 400 TeV and 2 PeV data sets of IceTop can be fitted with Goldreich-Sridhar turbulence and a broad resonance function, but other possibilities are not excluded. We then present our first numerical calculations of the CRA down to 3 TeV energies in 3D isotropic Kolmogorov turbulence. At these low energies, the large-scale CRA aligns well with the direction of local magnetic field lines around the observer. In this type of turbulence, the CR intensity is flat in a broad region perpendicular to field lines. Even though the CRA is quite gyrotropic, we show that the local configuration of the turbulence around the observer does result in the appearance of weak, "non-gyrotropic" small-scale anisotropies, which contain information on the local turbulence level.

We study particle acceleration at the termination shock of a striped pulsar wind by integrating trajectories in a prescribed model of the magnetic field and flow pattern. Drift motion on the shock surface maintains either electrons or positrons on "Speiser" orbits in a ring-shaped region close to the equatorial plane of the pulsar, enabling them to be accelerated to very high energy by the first-order Fermi mechanism. A power-law spectrum results: $dN_{\rm e}/d\gamma\propto\gamma^{\alpha_{\rm e}}$, where $\alpha_{\rm e}$ lies in the range $-1.8$ to $-2.4$ and depends on the downstream turbulence level. For sufficiently strong turbulence, we find $\alpha_{\rm e} \simeq -2.2$, and both the photon index and the flux of $1-100$ keV X-rays from the Crab Nebula, as measured by NuSTAR, can be reproduced. The particle spectrum hardens to $\alpha_{\rm e} \simeq -1.8$ at lower turbulence levels, which may explain the hard photon index observed by the Chandra X-ray Observatory in the central regions of the Nebula.

We calculate the diffusion coefficients of charged cosmic rays (CR) propagating in regular and turbulent magnetic fields. If the magnetic field is dominated by an isotropic turbulent component, we find that CRs reside too long in the Galactic disc. As a result, CRs overproduce secondary nuclei like boron for any reasonable values of the strength and the coherence length of an isotropic turbulent field. We conclude therefore that the propagation of Galactic CRs has to be strongly anisotropic because of a sufficiently strong regular field and/or of an anisotropy in the turbulent field. As a consequence, the number of sources contributing to the local CR flux is reduced by a factor ${\cal O}(100)$ compared to the case of isotropic CR diffusion.

Core collapse supernovae (CCSNe) produce fast shocks which pervade the dense circum-stellar medium (CSM) of the stellar progenitor. Cosmic rays (CRs) if accelerated at these shocks can induce the growth of electromagnetic fluctuations in the foreshock medium. In this study, using a self-similar description of the shock evolution, we calculate the growth timescales of CR-driven instabilities. We select a sample of nearby core collapse radio supernova of type II and Ib/Ic. From radio data we infer the parameters which enter in the calculation of the instability growth times. We find that extended IIb SNe shocks can trigger fast intra day instabilities, strong magnetic field amplification and CR acceleration. In particular, the non-resonant streaming instability can contribute to about 50\% of the magnetic field intensity deduced from radio data. This results in the acceleration of CRs in the range 1-10 PeV within a few days after the shock breakout. In order to produce strong magnetic field amplification and CR acceleration a fast shocks pervading a dense CSM is necessary. In that aspect IIn supernov\ae~are also good candidates. But a detailed modeling of the blast wave dynamics coupled with particle acceleration is mandatory for this class of object before providing any firm conclusions. Finally, we find that the trans-relativistic object SN 2009bb even if it produces more modest magnetic field amplification can accelerate CRs up to 2-3 PeV within 20 days after the outburst.

Observations of extended gamma-ray emission around Galactic cosmic-ray (CR) sources can be used as novel probes of interstellar magnetic fields. Using very-high-energy gamma-ray data from the HAWC Observatory, we place constraints on the properties of the magnetic turbulence within $\approx 25$ pc from Geminga. We inject and propagate individual CR electrons in 3D realizations of turbulent magnetic fields, calculate the resulting gamma-ray emission, and compare with HAWC measurements of this region. We find that HAWC data is compatible with expectations for Kolmogorov or Kraichnan turbulence, and can be well fitted for reasonable coherence lengths and strengths of the turbulence, despite implying a CR diffusion coefficient significantly smaller than those suggested by Galactic CR propagation codes. The best fit is found for a coherence length $L_{\rm c} \approx 1$ pc and a magnetic field strength $B_{\rm rms} \approx 3 \mu$G, and the preferred value for $L_{\rm c}$ increases with $B_{\rm rms}$. Moreover, the apparent lack of strong asymmetries in the observed emission allows us to constrain the coherence length to $L_{\rm c} \lesssim 5$ pc in this region.

The ~400 MeV flaring emission from the Crab Nebula is naturally explained as the result of an abrupt reduction in the mass-loading of the pulsar wind. Very few particles are then available to carry the current required to maintain wave activity, causing them to achieve high Lorentz factors. When they penetrate the Nebula, a tightly beamed, high luminosity burst of hard gamma-rays results, with characteristics similar to the observed flares. This mechanism may operate in other powerful pulsars, such as J0537-6910 (PWN N 157B), B0540-69, B1957+20 and J0205+6449 (3C 58).

We study a scenario in which the Fermi bubbles are formed through a Galactocentric outflow of gas and pre-accelerated cosmic-rays (CR). We take into account CR energy losses due to proton-proton interactions with the gas present in the bubbles, and calculate the associated gamma-ray emission. We find that CRs diffusing and advecting within a breeze outflow result in an approximately flat surface brightness profile of the gamma-ray emission, as observed by Fermi satellite. Finally, we apply similar outflow profiles to larger Galactocentric radii, and investigate their effects on the CR spectrum and boron-to-carbon ratio. Hardenings can appear in the spectrum, even in cases with equal CR diffusion coefficients in the disk and halo.

We show that the mysterious, rapidly variable emission at ~400 MeV observed from the Crab Nebula by the AGILE and Fermi experiments could be the result of a sudden drop in the mass-loading of the pulsar wind. The current required to maintain wave activity in the wind is then carried by very few particles of high Lorentz factor. On impacting the Nebula, these particles produce a tightly beamed, high luminosity burst of hard gamma-rays, similar to those observed. This implies (i) the emission is synchrotron radiation in the toroidal field of the Nebula, and, therefore, linearly polarized and (ii) this mechanism potentially contributes to the gamma-ray emission from other powerful pulsars, such as the Magellanic Cloud objects J0537-6910 and B0540-69.

We calculate the large-scale cosmic-ray (CR) anisotropies predicted for a range of Goldreich-Sridhar (GS) and isotropic models of interstellar turbulence, and compare them with IceTop data. In general, the predicted CR anisotropy is not a pure dipole; the cold spots reported at 400 TeV and 2 PeV are consistent with a GS model that contains a smooth deficit of parallel-propagating waves and a broad resonance function, though some other possibilities cannot, as yet, be ruled out. In particular, isotropic fast magnetosonic wave turbulence can match the observations at high energy, but cannot accommodate an energy dependence in the shape of the CR anisotropy. Our findings suggest that improved data on the large-scale CR anisotropy could provide a valuable probe of the properties - notably the power-spectrum - of the interstellar turbulence within a few tens of parsecs from Earth.

We study the role of local interstellar turbulence in shaping the large-scale anisotropy in the arrival directions of TeV-PeV cosmic-rays (CRs) on the sky. Assuming pitch-angle diffusion of CRs in a magnetic flux tube containing the Earth, we compute the CR anisotropy for Goldreich-Sridhar turbulence, and for isotropic fast modes. The narrow deficits in the 400 TeV and 2 PeV data sets of IceTop can be fitted for some parameters of the turbulence. The data also rule out a part of the parameter space. The shape of the CR anisotropy may be used as a local probe of the still poorly known properties of the interstellar turbulence and of CR transport.

Motivated by the discovery of the non-thermal Fermi bubble features both below and above the Galactic plane, we investigate a scenario in which these bubbles are formed through Galacto-centric outflow. Cosmic rays (CR) both diffusing and advecting within a Galactic breeze outflow, interacting with the ambient gas present, give rise to gamma-ray emission, providing an approximately flat surface brightness profile of this emission, as observed. Applying the same outflow profile further out within the disk, the resultant effects on the observable CR spectral properties are determined. A hardening in the spectra due to the competition of advective and diffusive propagation within a particular energy range is noted, even in the limiting case of equal CR diffusion coefficients in the disk and halo. It is postulated that this hardening effect may relate to the observed hardening feature in the CR spectrum at a rigidity of $\approx 200$ GV.

We investigate how the extragalactic proton component derived within the "escape model" can be explained by astrophysical sources. We consider as possible cosmic ray (CR) sources normal/starburst galaxies and radio-loud active galactic nuclei (AGN). We find that the contribution to the total extragalactic proton flux from normal and starburst galaxies is only subdominant and does not fit the spectral shape deduced in the escape model. In the case of radio-loud AGN, we show that the complete extragalactic proton spectrum can be explained by a single source population, BL Lac/FR I, for any of the potential acceleration sites in these sources. We calculate the diffuse neutrino and $\gamma$-ray fluxes produced by these CR protons interacting with gas inside their sources. For a spectral slope of CRs close to $\alpha=2.1-2.2$ as suggested by shock acceleration, we find that these UHECR sources contribute the dominant fraction of both the isotropic $\gamma$-ray background and of the extragalactic part of the astrophysical neutrino signal observed by IceCube.

We show that the cosmic ray (CR) knee can be entirely explained by energy-dependent CR leakage from the Milky Way, with an excellent fit to all existing data. We test this hypothesis calculating the trajectories of individual CRs in the Galactic magnetic field. We find that the CR escape time $\tau_{\rm esc}(E)$ exhibits a knee-like structure around $E/Z={\rm few}\times 10^{15}$ eV for small coherence lengths and strengths of the turbulent magnetic field. The resulting intensities for different groups of nuclei are consistent with the ones determined by KASCADE and KASCADE-Grande, using simple power-laws as injection spectra. The transition from Galactic to extragalactic CRs is terminated at $\approx 2\times 10^{18}$ eV, while extragalactic CRs contribute sizeable to the subdominant proton flux already for $\gtrsim 2\times 10^{16}$ eV. The natural source of extragalactic CRs in the intermediate energy region up to the ankle are in this model normal and starburst galaxies. The escape model provides a good fit to $\ln(A)$ data; it predicts that the phase of the CR dipole varies strongly in the energy range between $1\times 10^{17}$ and $3\times 10^{18}$ eV, while our estimate for the dipole magnitude is consistent with observations.

During a supernova explosion, a radiation-dominated shock (RDS) travels through its progenitor. A collisionless shock (CS) is usually assumed to replace it during shock breakout (SB). We demonstrate here that for some realistic progenitors enshrouded in optically thick winds, such as possibly SN 2008D, a CS forms deep inside the wind, soon after the RDS leaves the core, and therefore significantly before SB. The RDS does not survive the transition from the core to the thick wind when the wind close to the core is not sufficiently dense to compensate for the $r^{-2}$ dilution of photons due to shock curvature. This typically happens when the shock velocity is $\lesssim 0.1 {\rm c} \, (\frac{u_{\rm w}}{10\,{\rm km/s}}) (\frac{\dot{M}}{5 \cdot 10^{-4} \, {\rm M}_\odot {\rm /yr}})^{-1} (\frac{r_\ast}{10^{13}\,{\rm cm}})$, where $u_{\rm w}$, $\dot{M}$ and $r_\ast$ are respectively the wind velocity, mass-loss rate and radius of the progenitor star. The radiative CS results in a hard spectrum of the photon flash at breakout, which would produce an X-ray flash. Cosmic ray acceleration would start before SB, for such progenitors. A fraction of secondary TeV neutrinos can reach the observer up to more than ten hours before the first photons from breakout, providing information on the invisible layers of the progenitor.

The propagation of TeV-PeV cosmic rays (CR) in our Galaxy can be described as a diffusive process. We discuss here two effects, with important observational consequences, that cannot be predicted by the diffusion approximation in its usual form. First, we present an explanation for the CR anisotropies observed at small angular scales on the sky. We show that the local magnetic field configuration within a CR mean free path from Earth naturally results in CR flux anisotropies at small and medium scales. Second, we point out that TeV-PeV CRs should be expected to diffuse strongly anisotropically in the interstellar medium on scales smaller than the maximum scale of spatial fluctuations of the field, $\sim 100$ pc. This notably questions the usual assumptions on CR diffusion around sources.

We investigate the possibility that the cosmic ray (CR) knee is entirely explained by the energy-dependent CR leakage from the Milky Way. We test this hypothesis calculating the trajectories of individual CRs with energies between $E/Z=10^{14}$ eV and $10^{17}$ eV propagating them in the regular and turbulent Galactic magnetic field. We find a knee-like structure of the CR escape time $\tau_{\rm esc}(E)$ around $E/Z={\rm few}\times 10^{15}$ eV for a coherence length $l_{\rm c} \simeq 2$ pc of the turbulent field, while the decrease of $\tau_{\rm esc}(E)$ slows down around $E/Z\simeq 10^{16}$ eV in models with a weak turbulent magnetic field. Assuming that the injection spectra of CR nuclei are power-laws, the resulting CR intensities in such a turbulence are consistent with the energy spectra of CR nuclei determined by KASCADE and KASCADE-Grande. We calculate the resulting CR dipole anisotropy as well as the source rate in this model.

Analyses of TeV-PeV cosmic ray (CR) diffusion around their sources usually assume either isotropic diffusion or anisotropic diffusion due to the regular Galactic magnetic field. We show that none of them are adequate on distances smaller than the maximal scale Lmax ~ 100 pc of fluctuations in the turbulent interstellar magnetic field. As a result, we predict anisotropic gamma-ray emissions around CR proton and electron sources, even for uniform densities of target gas. The centers of extended emission regions may have non-negligible offsets from their sources, leading to risks of misidentification. Gamma-rays from CR filaments have steeper energy spectra than those from surrounding regions. We point out that gamma-ray telescopes can be used in the future as a new way to probe and deduce the parameters of the interstellar magnetic field.

Galactic cosmic ray (CR) acceleration to the knee in the spectrum at a few PeV is only possible if the magnetic field ahead of a supernova remnant (SNR) shock is strongly amplified by CR escaping the SNR. A model formulated in terms of the electric charge carried by escaping CR predicts the maximum CR energy and the energy spectrum of CR released into the surrounding medium. We find that historical SNR such as Cas A, Tycho and Kepler may be expanding too slowly to accelerate CR to the knee at the present time.

In the energy range from ~ 10^12 eV to ~ 10^15 eV, the Galactic cosmic ray flux has anisotropies both on large scales, with an amplitude of the order of 0.1%, and on scales between ~ 10 and ~ 30 degrees, with amplitudes smaller by a factor of a few. With a diffusion coefficient inferred from Galactic cosmic ray chemical abundances, the diffusion approximation predicts a dipolar anisotropy of comparable size, but does not explain the smaller scale anisotropies. We demonstrate here that energy dependent smaller scale anisotropies naturally arise from the local concrete realization of the turbulent magnetic field within the cosmic ray scattering length. We show how such anisotropies could be calculated if the magnetic field structure within a few tens of parsecs from Earth were known.

We constrain the energy at which the transition from Galactic to extragalactic cosmic rays occurs by computing the anisotropy at Earth of cosmic rays emitted by Galactic sources. Since the diffusion approximation starts to loose its validity for E/Z >~ 10^(16-17) eV, we propagate individual cosmic rays (CRs) using Galactic magnetic field models and taking into account both their regular and turbulent components. The turbulent field is generated on a nested grid which allows spatial resolution down to fractions of a parsec. Assuming sufficiently frequent Galactic CR sources, the dipole amplitude computed for a mostly light or intermediate primary composition exceeds the dipole bounds measured by the Auger collaboration around E ~ 10^18 eV. Therefore, a transition at the ankle or above would require a heavy composition or a rather extreme Galactic magnetic field with strength >~ 10 muG. Moreover, the fast rising proton contribution suggested by KASCADE-Grande data between 10^17 eV and 10^18 eV should be of extragalactic origin. In case heavy nuclei dominate the flux at E >~ 10^18 eV, the transition energy can be close to the ankle, if Galactic CRs are produced by sufficiently frequent transients as e.g. magnetars.

We investigate the diffusion of cosmic rays (CR) close to their sources. Propagating individual CRs in purely isotropic turbulent magnetic fields with maximal scale of spatial variations Lmax, we find that CRs diffuse anisotropically at distances r <~ Lmax from their sources. As a result, the CR densities around the sources are strongly irregular and show filamentary structures. We determine the transition time t* to standard diffusion as t* ~ 10^4 yr (Lmax/150 pc)^b (E/PeV)^(-g) (Brms/4 muG)^g, with b ~ 2 and g = 0.25-0.5 for a turbulent field with Kolmogorov power spectrum. We calculate the photon emission due to CR interactions with gas and the resulting irregular source images.

In this work, we present detailed simulations for propagation of ultra-high energy (UHE) heavy nuclei, with E > 60 EeV, within recent Galactic Magnetic Field (GMF) models. We investigate the impacts of the regular and turbulent components of the GMF. We show that with UHE heavy nuclei, there is no one-to-one correspondence between the arrival directions of cosmic rays (CR) measured at Earth and the direction of their extragalactic sources. Sources can have several distorted images on the sky. We compute images of galaxy clusters and of the supergalactic plane in recent GMF models and show the challenges, and possibilities, of "UHECR astronomy" with heavy nuclei. Finally, we present a quantitative study of the impact of the GMF on the (de-)magnification of source fluxes, due to magnetic lensing effects. We find that for 60 EeV iron nuclei, sources located in up to about one fifth of the sky would have their fluxes so strongly demagnified that they would not be detectable at Earth, even by the next generation of UHECR experiments.

In this work we study how the turbulent component of the Galactic magnetic field (GMF) affects the propagation of ultrahigh energy heavy nuclei. We investigate first how the images of individual sources and of the supergalactic plane depend on the properties of the turbulent GMF. Then we present a quantitative study of the impact of the turbulent field on (de-) magnification of source fluxes, due to magnetic lensing effects. We also show that it is impossible to explain the Pierre Auger data assuming that all ultrahigh energy nuclei are coming from Cen A, even in the most favorable case of a strong, extended turbulent field in the Galactic halo.

We present a new method to search for heavy nuclei sources, on top of background, in the Ultra-High Energy Cosmic Ray data. We apply it to the 69 events with energies E>55 EeV published by the Pierre Auger Collaboration. We find a set of events for which the method reconstructs the source near the Virgo galaxy cluster. The probability to have a comparable set of events in some background is ~ 0.7 %. The reconstructed source is located at ~ 8.5 degrees from the active galaxy M87. The probability to reconstruct the source at less than 10 degrees from M87 for data already containing a comparable set of events is ~ 0.4 %. This may be a hint at the Virgo galaxy cluster as an ultra-high energy heavy nuclei source. We discuss the capability of current and near future experiments to test this possibility. Such a scenario gives a self-consistent description of the Auger anisotropy and composition data at the highest energies.

We propose a new way to detect individual bright Ultra-High Energy Cosmic Ray (UHECR) sources above background if the Galactic Magnetic Field (GMF) gives the main contribution to UHECR deflections. This method can be directly applied to maps given by experiments. It consists in starting from at least two high energy events above 6x10^19 eV, and looking at lower energy tails. We test the efficiency of the method and investigate its dependence on different parameters. In case of detection, the source position and the local GMF deflection power are reconstructed. Both reconstructions are strongly affected by the turbulent GMF. With the parameters adopted in this study, for 68 % of reconstructed sources, the angular position is less than one degree from the real one. For typical turbulent field strengths of 4 \mu G at the Earth position and 1.5 kpc extension in the halo, one can reconstruct the deflection power with 25 % precision in 68 % of cases.

We propose a new method to search for heavy nuclei sources, on top of background, in the Ultra-High Energy Cosmic Ray data. We apply this method to the 69 events recently published by the Pierre Auger Collaboration and find a tail of events for which it reconstructs the source at a few degrees from the Virgo galaxy cluster. The reconstructed source is located at ~ 8.5 degrees from M87. The probability to have such a cluster of events in some random background and reconstruct the source position in any direction of the sky is about 7 x 10^(-3). The probability to reconstruct the source at less than 10 degrees from M87 in a data set already containing such a cluster of events is about 4 x 10^(-3). This may be a hint at the Virgo cluster as a bright ultra-high energy nuclei source. We investigate the ability of current and future experiments to validate or rule out this possibility, and discuss several alternative solutions which could explain the existing anisotropy in the Auger data.

Observations are consistent with a significant fraction of heavy nuclei in the cosmic ray flux above a few times 10^19 eV. Such nuclei can be deflected considerably in the Galactic magnetic field, with important implications for the search of their sources. We perform detailed simulations of heavy nuclei propagation within recent Galactic magnetic field models. While such models are not yet sufficiently constrained to predict deflection maps in detail, we find general features of the distribution of (de-) magnified flux from sources. Since in most theoretical models sources of heavy nuclei are located in the local large scale structure of galaxies, we show examples of images of several nearby galaxy clusters and of the supergalactic plane. Such general features may be useful to develop efficient methods for source reconstruction from observed ultrahigh energy cosmic ray arrival directions.

In this paper, we suggest a new way to identify single bright sources of Ultra High Energy Cosmic Rays (UHECR) on the sky, on top of background. We look for doublets of events at the highest energies, E > 6 x 10^19 eV, and identify low energy tails, which are deflected by the Galactic Magnetic Field (GMF). For the sources which are detected, we can recover their angular positions on the sky within one degree from the real ones in 68% of cases. The reconstruction of the deflection power of the regular GMF is strongly affected by the value of the turbulent GMF. For typical values of 4 microG near the Earth, one can reconstruct the deflection power with 25% precision in 68% of cases.