Magnetic Deflections of Ultra-High Energy Cosmic Raysfrom Centaurus A

Magnetic Deflections of Ultra-High Energy Cosmic Rays
from Centaurus A

Azadeh Keivani, Glennys R. Farrar, and Michael Sutherland Department of Physics & Astronomy, Louisiana State University, Baton Rouge, LA
Center for Cosmology and Particle Physics, Department of Physics, New York University, New York, NY
Department of Physics & Center for Cosmology and AstroParticle Physics, The Ohio State University, Columbus, OH

We present the results of a study that simulates trajectories of ultra-high energy cosmic rays from Centaurus A to Earth, for particle rigidities from EV to 100 EV, i.e., covering the possibility of primary particles as heavy as Fe nuclei with energies exceeding 50 EeV. The Galactic magnetic field is modeled using the recent work of Jansson and Farrar (JF12) which fitted its parameters to match extragalactic Faraday rotation measures and WMAP7 synchrotron emission maps. We include the random component of the GMF using the JF12 3D model for and explore the impact of different random realizations, coherence length and other features on cosmic ray deflections. Gross aspects of the arrival direction distribution such as mean deflection and the RMS dispersion depend mainly on rigidity and differ relatively little from one realization to another. However different realizations exhibit non-trivial substructure whose specific features vary considerably from one realization to another, especially for lower rigidities. At the lowest rigidity of 2 EV, the distribution is broad enough that it might be compatible with a scenario in which Cen A is the principle source of all UHECRs. No attempt is made here to formulate a robust test of this possibility, although some challenges to such a scenario are noted.

Ultra-High Energy Cosmic Rays, Galactic Magnetic Fields, Centaurus A
journal: Astroparticle Physics

1 Introduction

One of the biggest mysteries in particle astrophysics is the origin of ultra-high energy cosmic rays (UHECRs). There are few source types that appear capable of accelerating the particles to energies as high as has been observed eV). Moreover, UHECRs with energies above eV (50 EeV) lose energy through interactions with cosmic microwave background photons so that any that are observed must have originated within about 200 Mpc zk (); g (). Centaurus A (Cen A) is the closest active galactic nucleus (AGN) to Earth and is thus one of the prime source candidates. The Pierre Auger Observatory has reported a modest excess of cosmic ray arrival directions from within about 20 of the celestial location of Cen A in its UHECR data Auger2010 (). This further increases the interest to study the propagation of UHECRs through the Galactic magnetic fields (GMF) from Cen A.

UHECRs are mostly or entirely charged particles Auger2009 () so they are deflected in the magnetic fields that they traverse on their way to Earth. In this study, we simulate trajectories of cosmic rays from Cen A to Earth to examine the deflections that occur due to the GMF and turbulent field realizations. We construct the GMF according to the recent model of Jansson and Farrar JF12model (); JF12random () (JF12). The functional form of the field used in JF12 model is more general than previously considered and is tuned to fit the observed rotation measure and polarized synchrotron data. An overview of the model is given in Section 2.

This study is an extension of the work of Farrar et al. FJFG () which simulated UHECRs from the center of Cen A considering only the coherent field and four rigidities (energy divided by the particle’s charge: ) of 160, 80, 40 and 20 EV. Here we simulate additional rigidities extending down to 2 EV from the center of Cen A as well as consider the northern and southern radio lobes and explore the impact of adding several realizations of a Kolmogorov random field to the coherent JF12 GMF model.

Cen A was proposed to be the possible source of most cosmic rays by Farrar and Piran FP () assuming constraints on the values of turbulent extragalactic magnetic fields; see also Ginzburg (). Since the Pierre Auger Observatory released its first set of UHECR events, there have been several studies on the correlation of UHECRs with Cen A. For example, Gorbunov et al. discuss the possibility that Cen A could be the source of a significant portion of the observed events Gorbunov (). Wibig and Wolfendale assume that all observed UHECRs originate from three sources including Cen A Wibig (). Hardcastle et al. argue the possibility that the giant radio lobes of Cen A can accelerate UHECRs Hardcastle (). Related works include Kachelriess (); Rachen (); Fargion1 (); Fargion2 (); Rieger (); Moskalenko (); Nagar (). Ref. Giacinti () also examined the effects of turbulence on UHECRs, however with an older model of the GMF and turbulent realizations that were not constrained by data.

The propagation method is briefly explained in Section 3, some properties of Cen A are given in Section 4, and finally the results of the simulation are discussed in Sections 5 and 6.

2 JF12 GMF model

JF12 is a recent GMF model that was fit to the WMAP7 Galactic synchrotron emission maps (120,000 pixels of data) wmap () and more than forty thousand extragalactic rotation measures to constrain its parameters. This model includes a large scale regular component, a so-called “striated” random component that is found to follow the larger geometric features of the regular component JF12model (), and a turbulent random field component JF12random (). The random component of the JF12 model is constrained by the WMAP total Galactic synchrotron emission map JF12random ().

In this study, we first consider the large scale regular component of the JF12 model. The component contains a disk field component, a toroidal halo component and an out-of-plane component as suggested by observations of external galaxies. Then, for a subset of rigidities, we add to this several realizations of a turbulent random field whose r.m.s value is that of the JF12 random field model JF12random (), each characterized by a Kolmogorov power spectrum with different maximum coherence scales. The striated field is not included. Its effect is analogous to smearing the cosmic ray rigidity, i.e., it smears arrival directions along the arc predicted as a function of rigidity for the coherent field alone. Using the striation parameter determined in JF12random (), the smearing from the striated field can be estimated to be of the deflection in the coherent field, where is the distance propagated through the region of significant and is the typical coherence length of the striated random field.

We ignore the extragalactic magnetic field in this study because Cen A is relatively close to us and the strength of the extragalactic magnetic field is expected to be much smaller than of the GMF dolag04 (); dar (). The typical cosmic ray deflection magnitude in a turbulent extragalactic magnetic field () can be estimated FJFG () to be,


where is the source distance, and are the r.m.s value and coherence length of the magnetic field, is the charge of the UHECR in units of the proton charge, and is the UHECR energy in units of 100 EeV. A 100 EeV proton originating from Cen A would thus be deflected by only in such a turbulent extragalactic magnetic field, which is small compared to that from the regular and random Galactic fields of JF12 determined below.

3 Tracking

We use the latest release (v3.0.4) of CRT sutherland () to inject and propagate UHECRs in the JF12 model. CRT uses adaptive Runge-Kutta integration methods to determine the trajectory of a charged particle through a magnetic field according to the relativistic Lorentz force.

The method we use to track the cosmic rays through the GMF is “backtracking” particles from the Earth and recording their locations and directions when they exit the Galactic field region. This is accomplished by following the antiparticle, which simply means changing the sign of the intended particle’s charge and reversing its incoming velocity vector. The propagation of a particle stops when it hits the pre-defined boundary of the Galaxy since the extragalactic magnetic field is assumed to have negligible magnitude.

Figure 1: High resolution radio image of Cen A at 1.4 GHz showing the center of Cen A and its radio lobes in equatorial coordinates. The northern lobe is at the top of the image and the southern lobe is at the bottom. Image is Figure 3 of Feain-etal ().

We backtrack cosmic rays from isotropically distributed initial locations using HEALPix healpix (). In this study, we pick HEALPix resolution index 11 (or 9) which tiles the sky with more than (or ) pixels, each covering deg (or deg). We backtrack one particle with a fixed rigidity from the center of each pixel. The Larmor radius of a relativistic charged particle in a magnetic field depends only on its rigidity, so the use of rigidity allows a single simulation to represent a variety of primary particle compositions with appropriately scaled total primary energies. An event with rigidity could be a proton with energy or an iron nucleus with energy . Rigidities considered start with log()=2 (or EV) and end at log()=0.3 (or EV); for the coherent field, steps of 0.02-0.05 in log() are used. When the random field is included the simulations are significantly more time-consuming therefore only four values of rigidity have been used: log()=2, 1.5, 1.0, and 0.5 (or EV), respectively. For each rigidity, we identify particles that backtrack to regions of a specified solid angle within a radius of 1 or 3 around a chosen source location; in the first case this is the center of Cen A. We repeat the process for source locations centered on the tips of the northern and southern radio lobe.

4 Centaurus A

Centaurus A (NGC 5128) is a Fanaroff-Riley class I FanaroffRiley () massive elliptical radio galaxy in the constellation of Centaurus. It is located at a distance of Mpc Harris () and is the closest radio galaxy to Earth. It has a radio luminosity of Cooper () and is a uniquely well-studied source that has been characterized over the entire electromagnetic spectrum. Cen A has an active galactic nucleus in its center which ejects matter in relativistic jets: a pair of asymmetric nuclear jets and a pair of inner jets Feain-etal (). The inner jets extend about 1.4 kpc from the center Burns () whereas the inner radio lobes extend a further 5 kpc from the nucleus Clarke (). Cen A has also middle and outer radio lobes that extend to about 14 and 500 kpc, respectively.

Figure 1 shows a high resolution radio image of Cen A Feain-etal (). The image has angular dimensions of translating to a physical size of 500 kpc 200 kpc  Israel () at the distance of Cen A. The center of Cen A is taken to be in Galactic coordinates. We define and to be the tip of the northern and southern radio lobes, respectively. In equatorial coordinates these directions correspond to , , and .

5 UHECR Deflections in JF12 Coherent Field

(a) 2 - 100 EV
(b) EV
Figure 2: Arrival directions of the simulated events originating from within 3 around the center of Cen A, in the JF12 coherent field, for (a) EV to 100 EV and (b) focussing on EV, for which the centroid and standard deviations of the direction distribution are shown by the dot and error bars.
Figure 3: Centroids and standard deviations of the arrival direction distributions of events with the same rigidity. Their source directions lie within 3 of the center of Cen A. The map is projected in Galactic coordinates.

In the following we consider UHECRs originating within or within of the center of Cen A, for different cases. For the lower resolution HEALPix level 9 simulations we use the locus. The HEALPix level 11 simulations give us a similar event statistics with a radius. A angular window encompasses a sizeable portion of the radio galaxy (see Figure 1) which may be sites of acceleration and emission of cosmic rays. Since Cen A is about 3.8 Mpc away from the Earth, the () radius circle encompasses a physical distance of about 200 (66) kpc in radius from the center. As discussed earlier, the limits of the outer radio lobes lie about 500 kpc from the center so a value of 200 (66) kpc is a reasonable estimate for including acceleration regions associated with the central engine.

Figure 2(a) shows the arrival direction distribution of simulated events from within of the center of Cen A, for different rigidities. We calculate the longitude and latitude centroids of the arrival direction distributions as well as standard deviations in each direction separately; these are shown as black points and lines respectively for EV in Figure 2(b). The defined source region is a circular patch on the sky and the arrival distribution generally retains this shape.

The process is repeated for 35 different values of rigidity for the JF12 coherent field. Figure 3 shows the centroids and standard deviations of the arrival direction distributions from the center of Cen A for all simulated rigidities. Each color represents a different rigidity as indicated in the legend. These points illustrate how the arrival directions of UHECRs from the center of Cen A deviate from the source direction as a function of the rigidity of the particle. Particles with large values deviate much less than those with smaller values. Note that for EV, the arrival directions branch into two distinct regions on the plot; the centroid and standard deviations of such split distributions are computed and plotted separately. The distributions resemble the circular nature of the defined source region down to low simulated rigidities ( EV) below which the distributions deform into clearly noncircular regions. None of the arrival distributions has a significant number of outliers that would skew the appearance of Figure 3; the bars accurately reflect the distributions. Note that the regions for EV and 3.55 EV at the bottom-right side of Figure 2(a) appear to be cut but this is a real effect and the events with slightly smaller latitudes backtrack to a region to the northwest of Cen A and lie outside of the radius circle.

We separately study the observed directions of the simulated events coming from within 3 near the tips of the northern and southern radio lobes. We find that the general stream of the events are the same, but the appearance of separated regions split away from the main thread begins at lower rigidities ( EV) for the northern tip source and at higher rigidities ( EV) for the southern tip source. The main thread is shifted towards greater latitude values for the source at the northern tip and is shifted towards smaller latitude values for the source at the southern tip compared to the source at the center of Cen A. These shifts are only a few degrees for the lower rigidities. The shape and orientation of the threads remain the same as for the Cen A center source case in the sense that progressively lower rigidities trail towards lower longitude and latitude values. The multiple images observed for lower rigidities shift primarily in longitude by about ten degrees but maintain the same curved rigidity ordering and structure.

6 UHECR Deflections in JF12 GMF Model including Random Fields

(a) EV
(b) EV
(c) EV
(d) EV
Figure 4: Rigidity dependence of deflections of cosmic rays within of the center of Cen A with regular (grey) and regular plus KRF1 realization (red). The black asterisk shows the location of the center of Cen A.

In this section we discuss the effects of adding the detailed Kolmogorov random field component to the JF12 regular component. UHECR propagation in the full GMF including the random component requires an explicit realization of the field, so we have developed a method for propagation by implementing a box containing a grid of field magnitudes that follow a Kolmogorov spectrum AK-thesis (). This box is repeated throughout the GMF volume and is locally scaled by the JF12 r.m.s random field strength at the point of interest, i.e., the particle position. The total magnetic field is the vector sum of the regular and random field vectors. Independent random field realizations can be constructed by generating boxes with different random number seeds and different maximum coherence lengths.

Since the simulations are computationally expensive we first consider the regular field plus one random field realization (called KRF1), backtracking one event per HEALPix resolution index 9 pixel ( events over the sky), for four different rigidities: and 3.16 EV. The arrival direction distributions of events originating within of the center of Cen A in the regular field only and regular plus KRF1 field are compared in Figure 4. We see that the results of propagation only begin to change significantly at lower rigidity, such as 31.62 EV, between the regular-only and the regular-plus-random models. Even at lower rigidities, the arrival distribution centroids with and without the random field are only separated by a few degrees, but a large broadening of the distribution is evident when the random field is included.

Next, we consider a suite of random field realizations, to examine the effects of random field characteristics on the arrival direction distributions. We do this for a single rigidity which we choose to be EV. At EV the effects of the random field are clearly important while the dispersion remains small enough that the study is computationally feasible. For these simulations we increase the HEALPix resolution index to 11 ( events all-sky); this enables us to reduce the angular size of the source region to 1 instead of 3 around the center of Cen A and maintain an adequate number of events in the sample.

Realization (pc) Orient. Method Res. (deg) (EV) of Evts (deg) STDV (deg) STDV (deg)
KRF1 100 flipped 9 3.16 540 44.6 (19.8)
10 1349 24.7 (8.8)
31.6 1543 8.1 (2.4)
100 1421 2.5 (0.8)
KRF1 100 flipped 11 10 2292 24.7 (7.9)
KRF2 100 flipped 11 10 1877 20.7 (7.2)
KRF2 100 fixed 11 10 1827 21.6 (6.4)
KRF3 512 fixed 11 2 679 58.4 (32.5)
3.16 723 42.2 (24.6)
10 800 27.7 (13.9)
31.6 1592 14.8 (3.2)
100 4719 1.8 (0.4)
KRF4 512 fixed 11 10 1037 33.4 (12.4)
GARF 512 fixed 11 10 2067 19 (9.1)
KRF3-disk/3 512 fixed 11 10 1329 25.3 (5.0) 4.7
Table 1: Properties of each random field realization used in this paper. Column 2 shows the used in each field realization; Col. 3 indicates whether the realization box is fixed or randomly flipped throughout the Galactic volume; Cols. 4 and 5 are the HEALPix resolution index and the angular radius around the center of Cen A from which we select event source directions. Col. 6 shows the rigidity for the given simulation. Col. 7 gives the total number of received events from the given source locus; in the absence of magnetic lensing this would be 2155 and 3832 for res 9 in 3 and res 11 in 1  respectively. Col. 8 gives the mean deflection with respect to the source direction, and the standard deviation of the angular separation of the events relative to that mean. Columns 9 and 10 give the centroid and width of the arrival direction distribution in Galactic coordinates.
Figure 5: Histogram of field strength in the unit box for KRF2, KRF3, KRF4 and GARF.

Four realizations of the JF12 random field are generated using CRT code (KRF1, KRF2, KRF3, KRF4). Two of these realizations (KRF2 and KRF3) are implemented multiple times as will be discussed below. A single realization (GARF) is produced using the commonly used GARFIELD method garfield () in order to compare different field generation implementations. A final realization (KRF-disk/3) reduces by a factor of 3 the strength of the disk component of the JF12 random field, motivated by the estimates in GF-CRAS () of the potential impact of replacing the GALPROP cosmic ray electron model used in JF12 by a more structured such as predicted in some recent source models such as benyamin+13 (). A summary of the realizations and simulations is given in Table 1. The fields and results are explained in greater detail below.

Comparing Realizations

We can characterize the different random field realizations by histograms of the magnetic field strength at each grid point in the unit box, before scaling with the JF12 r.m.s field strengths. The field strengths of the unit box are normalized so that their r.m.s value is one. Figure 5 shows the field strength histogram for KRF2, KRF4 and GARF. It should be noted that KRF3 and KRF4 have the same field parameters, KRF2 has a lower than the other realizations, and GARF has the same parameters as KRF3 and KRF4 but is generated by a different method.

Table 1 lists specific details of and differences between the realizations, as described in the caption. The number of events arriving at Earth that originate within the source locus is an interesting diagnostic. In the absence of magnetic lensing, this number would simply be the number of events originating within the source locus: 2155 for 3 and res 9 and 3832 for 1 and res 11. However due to magnetic lensing, the number varies with the field model and rigidity, as recorded in Column 7 of Table 1. The mean angular deflections and RMS angular size of the arrival direction distribution, and the mean and standard deviations of the longitude and latitude of the events, are also listed in Table 1 for each simulation.

Taken together, the suite of realizations allow us to investigate how the following impact the arrival direction distributions: A) re-orientation or not of a realization within the simulation volume, B) variations between independent random realizations of fields with the same parameters, C) maximum coherence length, and D) method of producing the Kolmogorov spectrum:

Figure 6: Arrival direction distribution of events with EV, for source locations lying within of the center of Cen A; the black regions give the direction distribution in the regular only field and the black star shows the location of the center of Cen A. Left: KRF2 used in fixed orientation throughout the entire volume versus re-orienting different instances when filling the volume; KRF2-flipped in blue and KRF2-fixed in red. Right: The effect of different random number seeds with the same procedure in other respects; KRF1 in red and KRF2 in blue.

Independent Realizations with Same Parameters

KRF2 is used in a pair of simulations where the orientation of the realization box is either held fixed throughout the entire volume or is randomly re-oriented (flipped) from one location to another. The left panel of Fig. 6 illustrates the effect of different orientation methods implemented for a single realization. The gross features of size and location of the arrival direction distribution is relatively insensitive to the orientation method (c.f.,Table 1).

Figure 7: Sensitivity to coherence length with pc in red and pc in blue. The left panel shows KRF1&3 and the right panel KRF2&4; otherwise as in Fig. 6.

The right panel of Fig. 6 compares the arrival directions for KRF1 and KRF2 to see the effect of two independent realizations (random number seeds) with the same random field parameters. The last two columns of Table 1 show that these two distributions are only different up to statistical fluctuations.

Sensitivity to Maximum Coherence Length and Random Field Generation

Maximum coherence length is related to the maximum wavelength of the Fourier modes of the random field realization, . According to Equation 2.6 of HMRS2002 (), for a Kolmogorov spectrum and ample range between and , the coherence length is so that our choices of allow us to span the range of 20-100 pc coherence length. KRF1 and 2 are generated using pc and KRF3 and 4 with 512 pc. Figure 7 show comparisons between KRF1 and KRF3, and KRF2 and KRF4, respectively, illustrating the effect of varying the maximum coherence length on the arrival direction distribution for 10 EV rigidity. As expected, the differences from one realization to another are smaller for the smaller coherence length.

We can study the sensitivity to the tails of the distribution of values, by comparing KRF3 and KRF4 to GARF, since the coherence length (512 pc) is the same for each of these; the box orientation is also fixed throughout the simulation volume for all three realizations. The left panel of Fig. 8 shows arrival direction distributions for KRF3, KRF4 and GARF. We see that GARF produces less dispersion, which is not surprising given its narrower tails in values (Fig. 5).

A larger coherence length, or broader distribution, seems to increase the “cosmic variance” from one realization to another. It appears that for 10 EV UHECRs coming from the source direction of Cen A, differences from one realization to another can be sufficient to significantly impact the fraction of events whose trajectories come close to the Galactic plane, where the random field can be quite large in some arms JF12random (). One sees from Table 1 that the variance about the mean deflection is much larger for the realizations in which the mean arrival direction is below the plane than for those above. A subject of future work will be to see how this behavior changes for different source directions. (Note that KRF1 was implemented with a flipped orientation scheme whereas the others use a fixed scheme, but this was found to make only a minor difference in the distribution.)

Figure 8: Left panel: sensitivity to random field production method seen comparing use of GARF (red), KRF3 (blue) and KRF4 (green). Right panel: sensitivity to the RMS strength of the disk component of the random field; blue points are standard KRF3 and red points are with 3-times smaller random field in the disk (KRF3-disk/3).

Reduced random field in the Disk

Finally we see that reducing the disk random field strength by a factor of 3 relative to the JF12 random field model has little impact on the mean deflection but significantly decreases the arrival-direction spread. This study was motivated to examine sensitivity of Cen A UHECR deflections to uncertainties in the GMF model GF-CRAS ().

7 Summary and Conclusions

We have performed a detailed exploration of how the UHECR image of Cen A moves and fragments as a function of CR rigidity from 100 EV to 2 EV in the JF12 model of the GMF. This study shows that in the coherent field alone, even for deflections which can be large, the size and shape of the image changes only mildly until rigidities are below a few EV, as seen in Figs. 2(a) and  3. However the picture changes dramatically when random fields are included. To explore this, we have performed high angular resolution simulations for a series of fixed realizations of the random field, for a variety of random-field realizations, to learn how much the arrival direction distribution changes from one realization to another. (In many earlier studies of magnetic deflections, the random field was generated “on-the-fly” during propagation, such that each CR experiences a different field, leading to unrealistic additional smearing and smoothing of the average deflection map.) These simulations lead us to several conclusions, applicable for rigidities such that the UHECRs come close enough to the Galactic plane to feel the relatively strong random field in the disk ( EV for Cen A as the source). We used the JF12 model which is surely imperfect, but can be a guide to general features which can be expected:

  • At rigidities of EV and below, individual deflection magnitudes can be very large and the arrival direction distributions span a large fraction of the sky.

  • The mean deflection and the RMS spread in arrival directions for a given rigidity varies a few degrees from one realization to another. These properties depend much more on the rigidity of the CR than on the specific realization of the random field.

  • The distribution of arrival directions has considerable structure: it has hot-spots and edges, and would not typically be well-described by a 2-D gaussian.

  • The location and features of the structures vary significantly from one realization to another, even though the gross properties (mean, RMS) are similar.

  • The “transmission efficiency” of UHECRs from a given source – the number of CRs from that source which reach Earth, relative to the number which would reach us in the absence of the GMF – varies from one realization to another as well as with rigidity.

Figure 9: Sky map in Galactic coordinates of the arrival direction distributions, for UHECRs from within of the center of Cen A, propagated through the regular JF12 GMF plus KRF3 realization of the random JF12 model, for rigidities from 2-100 EV. The observed Pierre Auger events with energies above 55 EeV (from Auger2010 ()) are shown with black asterisks; note that the arrival locus for 100 EV is somewhat obscured by 2 Auger event symbols.

We close by displaying (Fig. 9) the simulated arrival direction distributions for a particular random field realization, for rigidities , 3.2, 10, 32 and 100 EV, in a sky map using Mollweide projection. Observed events with EeV published by the Pierre Auger Observatory are shown for comparison. The data are from Auger2010 (). A striking but difficult to interpret feature of Fig. 9, is that the simulated arrival direction distribution corresponding to EV, is so broad as to encompass most of the Auger events. A detailed analysis taking into account the Auger exposure, including the recently-released TA events TAaniso14 () and studying more realizations of the field model would be required to determine whether a single-source (Cen A) scenario could conceivably be viable. Since photo-disintegration during propagation is unimportant because Cen A is so close, the UHECR spectral cut-off would reflect the maximum-rigidity that Cen A can produce. At lower energies, one would expect a mixed composition reflecting the relative abundances at the source, with each component having a maximum energy . A potential Achilles-heel of such a model is that – no matter what the composition mix is, and even with a pure-Fe composition – the events above 55 EeV must cover a range of rigidities extending at a minimum between and EV. As can be seen from the distributions for and EV in Fig. 9, events with these higher rigidities would be expected to be much more clustered than the EV events, which may be incompatible with the observed approximate isotropy.

In concluding, we stress that the JF12 random field model implemented here was only a first attempt to estimate the random component of the GMF. As such, it is subject to significant revision as improved modeling and constraining data becomes available. In particular, the possible very broad distribution of arrival directions of iron nuclei from Cen A our studies reveal, may not apply to the true GMF. This underlines the importance of the observational challenge: to find a robust way to test such scenarios, given that the specifics of the structure due to the random magnetic fields cannot be predicted in detail. Information from individual event charge-assignments and from extending the energy range at which the and distributions give composition information, for example as envisaged in the Auger Upgrade proposal, will give valuable complementary constraints to arrival-direction studies.

Acknowledgments: We thank Jonathan Roberts and Ronnie Jansson for their collaboration in the initial phase of this work as well as Jim Matthews for his useful comments. We also thank our many colleagues in the Pierre Auger Collaboration for their various interactions and input. A.K. and M.S. acknowledge support from the Department of Energy under grants DE-0009926 and DE-FG02-91-EF0617. G.F. acknowledges support from the National Science Foundation and NASA under grants NSF-PHY-1212538 and NNX10AC96G; resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center, consisting of time on the Pleiades supercomputing cluster awarded to G. Farrar.


  • (1) G. T. Zatsepin, V. A. Kuz’min, Upper limit of the spectrum of cosmic rays, Journal of Experimental and Theoretical Physics Letters 4 (1966) 78–80.
  • (2) K. Greisen, End to the cosmic-ray spectrum?, Physical Review Letters 16 (17) (1966) 748–750. doi:10.1103/PhysRevLett.16.748.
  • (3) P. Abreu, et al. (Pierre Auger Collaboration), Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter, Astroparticle Physics 34 (5) (2010) 314–326. doi:10.1016/j.astropartphys.2010.08.010.
  • (4) J. Abraham, et al. (Pierre Auger Collaboration), Upper limit on the cosmic-ray photon fraction at EeV energies from the Pierre Auger Observatory, Astroparticle Physics 31 (6) (2009) 399–406. doi:10.1016/j.astropartphys.2009.04.003.
  • (5) R. Jansson, G. R. Farrar, A new model of the galactic magnetic field, The Astrophysical Journal 757 (1) (2012) 14. doi:10.1088/0004-637X/757/1/14.
  • (6) R. Jansson, G. R. Farrar, The Galactic magnetic field, The Astrophysical Journal 761 (1) (2012) L11. doi:10.1088/2041-8205/761/1/L11.
  • (7) G. R. Farrar, R. Jansson, I. J. Feain, B. M. Gaensler, Galactic magnetic deflections and Centaurus A as a UHECR source, Journal of Cosmology and Astroparticle Physics 1 (2013) 23. doi:10.1088/1475-7516/2013/01/023.
  • (8) G. R. Farrar, T. Piran, Deducing the source of ultrahigh energy cosmic rays (2000). arXiv:0010370.
  • (9) V. L. Ginzburg, S. I. Syrovatskii, Cosmic rays in metagalactic space, Soviet Astronomy 7 (1963) 357.
  • (10) D. S. Gorbunov, P. G. Tinyakov, I. I. Tkachev, S. V. Troitsky, On the interpretation of the cosmic-ray anisotropy at ultra-high energies (2008). arXiv:0804.1088.
  • (11) T. Wibig, A. W. Wolfendale, Heavy cosmic ray nuclei from extragalactic sources above ‘the ankle’, The Open Astronomy Journal 2 (1) (2009) 95–101. doi:10.2174/1874381100902010095.
  • (12) M. J. Hardcastle, C. C. Cheung, I. J. Feain, Ł. Stawarz, High-energy particle acceleration and production of ultra-high-energy cosmic rays in the giant lobes of Centaurus A, Monthly Notices of the Royal Astronomical Society 393 (3) (2009) 1041–1053. doi:10.1111/j.1365-2966.2008.14265.x.
  • (13) M. Kachelrieß, S. Ostapchenko, R. Tomàs, High energy radiation from Centaurus A, New Journal of Physics 11 (6) (2009) 065017. doi:10.1088/1367-2630/11/6/065017.
  • (14) J. P. Rachen, Ultra-high energy cosmic rays from radio galaxies revisited (2008). arXiv:0808.0349.
  • (15) D. Fargion, Light nuclei solving the AUGER puzzles: the Cen-A imprint, Physica Scripta 78 (4) (2008) 045901. doi:10.1088/0031-8949/78/04/045901.
  • (16) D. Fargion, Coherent and random UHECR spectroscopy of lightest nuclei along Cen A: Shadows on GZK Tau neutrinos spread in a near sky and time, Nuclear Instruments and Methods in Physics Research Section A 630 (1) (2011) 111–114. doi:10.1016/j.nima.2010.06.040.
  • (17) F. M. Rieger, F. A. Aharonian, Centaurus A as TeV -ray and possible UHE cosmic-ray source, Astronomy and Astrophysics 506 (3) (2009) L41–L44. doi:10.1051/0004-6361/200912562.
  • (18) I. V. Moskalenko, L. Stawarz, T. A. Porter, C. C. Cheung, On the possible association of ultra high energy cosmic rays with nearby active galaxies, The Astrophysical Journal 693 (2) (2009) 1261–1274. doi:10.1088/0004-637X/693/2/1261.
  • (19) N. M. Nagar, J. Matulich, Ultra-high energy cosmic rays detected by Auger and AGASA, Astronomy & Astrophysics 523 (2010) A49. doi:10.1051/0004-6361/200912702.
  • (20) G. Giacinti, M. Kachelrieß, D. V. Semikoz, G. Sigl, Ultrahigh energy nuclei in the turbulent galactic magnetic field, Astroparticle Physics 35 (2011) 192–200. doi:10.1016/j.astropartphys.2011.07.006.
  • (21) J. L. Weiland, N. Odegard, R. S. Hill, E. Wollack, G. Hinshaw, M. R. Greason, N. Jarosik, L. Page, C. L. Bennett, J. Dunkley, B. Gold, M. Halpern, A. Kogut, E. Komatsu, D. Larson, M. Limon, S. S. Meyer, M. R. Nolta, K. M. Smith, D. N. Spergel, G. S. Tucker, E. L. Wright, Seven-year Wilkinson microwave anisotropy probe (WMAP) observations: Planets and celestial calibration sources, The Astrophysical Journal Supplement Series 192 (2) (2011) 19. doi:10.1088/0067-0049/192/2/19.
  • (22) K. Dolag, D. Grasso, V. Springel, I. Tkachev, Mapping deflections of extragalactic ultrahigh-energy cosmic rays in magnetohydrodynamic simulations of the local universe, Journal of Experimental and Theoretical Physics Letters 79 (12) (2004) 583–587. doi:10.1134/1.1790011.
  • (23) A. Dar, A. De Rújula, Magnetic field in galaxies, galaxy clusters, and intergalactic space, Physical Review D 72 (12) (2005) 123002. doi:10.1103/PhysRevD.72.123002.
  • (24) M. Sutherland, B. Baughman, J. Beatty, CRT: A numerical tool for propagating ultra-high energy cosmic rays through Galactic magnetic field models, Astroparticle Physics 34 (4) (2010) 198–204. doi:10.1016/j.astropartphys.2010.07.002.
  • (25) I. J. Feain, T. J. Cornwell, R. D. Ekers, M. R. Calabretta, R. P. Norris, M. Johnston-Hollitt, J. Ott, E. Lindley, B. M. Gaensler, T. Murphy, E. Middelberg, S. Jiraskova, S. O’Sullivan, N. M. McClure-Griffiths, J. Bland-Hawthorn, The radio continuum structure of Centaurus A at 1.4 GHz, The Astrophysical Journal 740 (1) (2011) 17. doi:10.1088/0004-637X/740/1/17.
  • (26) K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, M. Bartelmann, HEALPix: A framework for high-resolution discretization and fast analysis of data distributed on the sphere, The Astrophysical Journal 622 (2) (2005) 759–771. doi:10.1086/427976.
  • (27) B. L. Fanaroff, J. M. Riley, The morphology of extragalactic radio sources of high and low luminosity, Monthly Notices of the Royal Astronomical Society 167 (1974) 31p–36p.
  • (28) G. L. H. Harris, M. Rejkuba, W. E. Harris, The distance to NGC 5128 (Centaurus A), Publications of the Astronomical Society of Australia 27 (04) (2013) 457–462. doi:10.1071/AS09061.
  • (29) B. F. C. Cooper, R. M. Price, D. J. Cole, A study of the decimetric emission and polarization of Centaurus A, Australian Journal of Physics 18 (1965) 589.
  • (30) J. O. Burns, E. D. Feigelson, E. J. Schreier, The inner radio structure of Centaurus A - clues to the origin of the jet X-ray emission, The Astrophysical Journal 273 (1983) 128. doi:10.1086/161353.
  • (31) D. A. Clarke, J. O. Burns, M. L. Norman, VLA observations of the inner lobes of Centaurus A, The Astrophysical Journal 395 (1992) 444. doi:10.1086/171663.
  • (32) F. Israel, Centaurus A - NGC 5128, Astronomy and Astrophysics Review 8 (4) (1998) 237–278. doi:10.1007/s001590050011.
  • (33) A. Keivani, Magnetic deflections of ultra-high energy cosmic rays from Centaurus A, Ph.D. thesis, Louisiana State University, (2013).
  • (34) F. S. Kitaura, T. A. Ensslin, Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization, Monthly Notices of the Royal Astronomical Society 2 (2008) 497–544. doi:10.1111/j.1365-2966.2008.13341.x.
  • (35) G. R. Farrar, The galactic magnetic field and ultrahigh-energy cosmic ray deflections (2014). arXiv:1405.3680.
  • (36) D. Benyamin, E. Nakar, T. Piran, N. J. Shaviv, Recovering the Observed B/C Ratio in a Dynamic Spiral-armed Cosmic Ray Model, The Astrophysical Journal 782 (2014) 34. arXiv:1308.1727, doi:10.1088/0004-637X/782/1/34.
  • (37) D. Harari, S. Mollerach, E. Roulet, F. Sánchez, Lensing of ultra-high energy cosmic rays in turbulent magnetic fields, Journal of High Energy Physics 03 (2002) 045. doi:10.1088/1126-6708/2002/03/045.
  • (38) Telescope Array Collaboration, Indications of intermediate-scale anisotropy of cosmic rays with energy greater than 57 EeV in the northern sky measured with the surface detector of the Telescope Array experiment (2014). arXiv:1404.5890.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description