# High energy protons from PKS 1333-33

\[

## Abstract

In this letter we give an account of the possible acceleration of protons in the outer radio lobes of the active galaxy PKS 1333-33. We also make estimates of the arrival energy spectrum.

PACS number(s): 96.40 - 98.70.S - 95.85.R - 13.85.T

\]

Radio galaxies considered as sources of cosmic rays (CR) beyond the Greisen-Zatsepin-Kuz’min cutoff  must be quite close to our own galaxy ($`z0.03`$). There are just a few objects within this range and, consequently, the expected ultra high energy CR distribution on the sky should be highly anisotropic. Assuming that the intergalactic magnetic field is near 1 nG, it should be expected an excess of CR detections at energies larger than $`5\times 10^{19}`$ eV within regions of angular radius $`\theta 15^{}`$ and centered at the positions of the nearest active galaxies . According to this picture, the southern CR-sky should be dominated by three outstanding sources: Cen A (the nearest radio galaxy) which would provide the most energetic particles detectable on Earth , Pictor A (a strong source with a flat radio spectrum) which would contribute with the larger CR flux , and PKS 1333-33 which might be a source of events similar to those recently detected in the Northern Hemisphere. In addition to these sources, there are other two southern candidates, Fornax A ($`z=0.057`$) and PKS 2152-69 ($`z=0.027`$), which could provide contributions to the CR flux above the cutoff. An approximate theoretical picture of the up to now unexplored ultra high energy CR southern sky is consequently at our disposal. The present letter is devoted to CR production in the outer region of the nearby southern active galaxy PKS 1333-33 and the spectral modifications arisen from their propagation through intergalactic space.

The radio galaxy PKS 1333-33 is made up of a core, two symmetric jets, and two extended radio lobes. Fig. 1 shows the entire source at 20 cm according to the VLA map obtained by Killeen et al. . The core has been identified with the E1 galaxy IC 4296, which has a redshift $`z`$=0.013 . The distance to the source is 35.2 $`h^1`$ Mpc, ($`h=H_0/100`$ km<sup>-1</sup> s Mpc, hereafter we shall adopt $`h=0.65`$ ). The large-scale structure of PKS 1333-33 has been studied in detail at 1.3, 2, 6, and 20 cm with an angular resolution of 3.2” . The jets are slightly bent, presumably as a consequence of the motion of the core with respect to the intergalactic medium. The total flux density of the source at 20 cm is $``$ 14.5 Jy. The integrated radio luminosity is $`2.5\times 10^{41}`$ erg s<sup>-1</sup> assuming a spectral index of $`\alpha =0.6`$ and frequency cutoffs at $`10^7`$ and $`10^{11}`$ Hz. Actually, the value $`\alpha =0.6`$ is correct only for the extended radio lobes. The spectral index steepens in the jets, reaching $`\alpha =1`$, while the core has flatter values: $`\alpha 0.3`$.

The hot spots of extended radio sources are regions of strong synchrotron emission . These regions are produced when the bulk kinetic energy of the jets ejected by a central active source (supermassive black hole plus accretion disk) is reconverted into relativistic particles and turbulent fields at a “working surface” in the head of the jets . The speed $`v_\mathrm{h}`$ with which the head of a jet advances into the intergalactic medium of particle density $`n_\mathrm{e}`$ can be obtained by balancing the momentum flux in the jet against the momentum flux of the surrounding medium. Measured in the frame comoving with the advancing head, $`v_\mathrm{h}v_\mathrm{j}[1+(n_\mathrm{e}/n_\mathrm{j})^{1/2}]^1`$, where $`n_\mathrm{j}`$ and $`v_\mathrm{j}`$ are the particle density and the velocity of the jet flow, respectively. Clearly, $`v_\mathrm{j}>v_\mathrm{h}`$ for $`n_\mathrm{e}n_\mathrm{j}`$, in such a way that the jet will decelerate. The result is the formation of a strong collisionless shock, which is responsible for particle reacceleration and magnetic field amplification. The acceleration of particles up to ultrarelativistic energies in the hot spots is the result of repeated scattering back and forth across the shock front. The particle deflection in this mechanism is produced by Alfvén waves in the turbulent magnetic field. Biermann and Strittmatter  have studied this process assuming that the energy density per unit of wave number of MHD turbulence is of Kolmogorov type. According to their calculations, the highest energy of protons injected in the intergalactic medium from the hot spot can be obtained by balancing the gains and losses in the diffusive shock acceleration process:

$`E_{\mathrm{max}}`$ $`=`$ $`7.8\times 10^5\beta _{_{\mathrm{jet}}}^{\mathrm{\hspace{0.33em}3}/2}u^{3/4}R_{_{\mathrm{hs}}}^{1/2}`$ (1)   
$`\times `$ $`B_5^{5/4}(1+Aa)^{3/4}\mathrm{EeV}`$ (2)   
where $`\beta _{_{\mathrm{jet}}}`$ stands for the jet velocity in units of $`c`$, $`u`$ is the ratio of turbulent to ambient magnetic energy density in the hot spot (of radius $`R_{\mathrm{hs}}`$ measured in kpc), $`B_5`$ is the magnetic field in units of $`10^5`$ G, $`a`$ is the ratio of photon to magnetic energy density, and $`A`$ gives a relative strength of $`\gamma p`$ interactions against the synchrotron emission.

An interesting feature of PKS 1333-33 is an intense region of synchrotron emission localized at the outer edge of the eastern lobe. This region can be considered as a “working surface” formed by the deceleration of the jet. This interpretation is supported by VLA polarimetric observations, which show a change in the field orientation from parallel to perpendicular to the jet axis . This change is probably due to the rearrangement of the field lines in the post-shock region. The degree of linear polarization in the eastern lobe is in the range 20%-40%. The synchrotron parameters estimated for this region are: the minimum energy density ($`ϵ_{\mathrm{min}}=2.1\times 10^{13}`$ erg cm<sup>-3</sup>), the minimum magnetic field ($`B_{\mathrm{min}}=1.5\times 10^6`$ G), the minimum total pressure ($`P_{\mathrm{min}}=1.1\times 10^{13}`$ dyne cm<sup>-2</sup>), and the degree of linear polarization ($`m`$ = 20-40%).

The leptonic component of the CRs will produce synchrotron emission with a spectrum given by $`S_\nu \nu ^\alpha `$, where $`\alpha =(\mathrm{\Gamma }1)/2`$. Since $`\alpha =0.6`$, we get $`\mathrm{\Gamma }=2.2`$. In what follows, we shall assume that electrons and protons in the source obey the same power law energy spectrum $`E^\mathrm{\Gamma }`$.

The degree of linear polarization expected for the synchrotron radiation when the magnetic field is homogeneous is

$$m_0(\mathrm{\Gamma })=\frac{3\mathrm{\Gamma }+3}{3\mathrm{\Gamma }+7}71\%$$
(3)   
However, the observed degree of polarization has a mean value of $`m(\mathrm{\Gamma })30\%`$. This fact can be explained by the presence of a turbulent component $`B_\mathrm{r}`$ in the field, in such a way that

$$m(\mathrm{\Gamma })=m_0(\mathrm{\Gamma })\frac{B_0^2}{B_0^2+B_\mathrm{r}^2}$$
(4)   
where $`B_0`$ stands for the homogeneous field. From Eq. (3) and (4) we get

$$B_\mathrm{r}1.2B_0$$
(5)   
and consequently $`u=B_\mathrm{r}^2/B_{\mathrm{total}}^20.6`$.

The radius of the acceleration region can be directly measured by means of a Gaussian fitting from the detailed VLA maps obtained by Killeen et al. , resulting $`R_{\mathrm{hs}}2.5h^1`$ kpc. The velocity of the jet is not well established. If the source is $`10^7`$ yr old, a velocity $`\beta _{\mathrm{jet}}0.2`$ can be estimated from an analysis of the energy budget . The value of $`A`$ in Eq. (2) has been computed by Biermann and Strittmatter  almost independently of the source parameters. They obtain $`A200`$. A value $`a0.01`$ seems to be reasonable for a source with the luminosity of PKS 1333-33 . Taking the above considerations into account (with a typical value for the total magnetic field in the spot of $`B_510`$ ), we obtain from Eq. (2) the maximun injection energy for protons, $`E_{\mathrm{max}}6\times 10^{20}`$ eV.

The diffusive shock acceleration process in the working surface of the eastern lobe leads to a power law particle spectrum ,

$$Q(E)=\kappa E^\mathrm{\Gamma }(E_0<E<E_{\mathrm{max}})$$
(6)   
where $`\kappa =(n_0/\mathrm{\Gamma }1)E_0^{1\mathrm{\Gamma }}`$ for $`E_{\mathrm{max}}E_0`$, and $`n_0`$ is the particle density in the source. Thus, we have the following proton injection spectrum

$$Q(E)E^{2.2},E_0<E<6\times 10^{20}\mathrm{eV}.$$
(7)   
Let us now estimate the losses that the high energy protons suffer during their journey through the intergalactic medium, and the consequent modification of the injection spectrum. This energy degradation originates in the interaction of the protons with the very low energy photons of the cosmic microwave background (CMB). For nearby sources (distances $`<`$ 100 Mpc) photomeson production is the dominant interaction at $`E>10^{19}`$ eV (if one considers the contribution coming from the tail of the Planck distribution) . The propagation of the protons through the relic gas of photons may be modelled by a kinetic transport equation. Although it is clear that the interactions of the ultra high energy protons with the CMB photons lead to a step-by-step energy loss (which should be included in the transport equation as a collision integral), hereafter we shall use the continuous energy loss approximation assuming straight line propagation which is expected to be in agreement with the physics of the problem at hand.

Using the formalism presented in  it is straightforward to compute the main characteristic of the evolved spectrum. It is well known that the conservation of the differential particle number implies that,

$$\frac{N}{t}=\frac{(b(E)N)}{E}+D^2N+Q,$$
(8)   
where, in the first term on the right, $`b(E)`$ is the mean rate at which particles lose energy ($`3.66\times 10^8E\mathrm{exp}\{287/E\}`$ EeV/yr ); the second term, the diffusion in the CMB, is found to be extremely small due to the low density of relic photons and is neglected in the following. The third term corresponds to the particle injection rate into the intergalactic medium given by Eq. (7). The solution of Eq. (8) reads,

$$N(E,t)=\kappa E_g^\mathrm{\Gamma }\frac{b(E_g)}{b(E)}$$
(9)   
with the injection energy $`E_g`$ fixed by the constraint,

$$3.66\times 10^8t\mathrm{yr}^1\mathrm{Ei}(287/E)+\mathrm{Ei}(287/E_g)=0$$
(10)   
being Ei$`(x)`$ the exponential-integral function. The evolution of the injection spectrum can be then conveniently represented by the modification factor $`\eta `$, defined as the ratio between the modified and the unmodified spectrum. In Fig. 2 we have plotted the relation obtained for $`\eta `$ in the case of PKS 1333-33. The energy loss by photomeson production creates the expected cutoff, and the resulting ultra high energy CRs (protons and neutrons) pile up just below the threshold energy of photopion production, forming a bump. However, the energies at which the cutoff and the bump appear seem to indicate that CRs above 100 EeV might be expected from the location in the sky of PKS 1333-33 ($`l313.7^{}`$, $`b27.7^{}`$).

The question whether this scenario is the correct one or not should be answered in a few years by the Pierre Auger Southern Observatory  (fluorescence detector plus ground array) as well as by the future eyes of the OWL  that will deeply watch into the CR-sky.

###### Acknowledgements.

We are indebted to N. E. B. Killeen, G. V. Bicknell, and R. D. Ekers for permission to reproduce the VLA image of PKS 1333-33. G. E. R. especially thanks the help by J. Blumina Romero. This work has been partially supported by FOMEC, CONICET, CLAF-CNPq, and the agency ANPCT.