Free Access
Issue
A&A
Volume 605, September 2017
Article Number A115
Number of page(s) 7
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201527871
Published online 19 September 2017

© ESO, 2017

1. Introduction

HESS J1303631 is one of the most prominent examples of the so-called very high energy (VHE; E> 100 GeV) γ-ray dark sources, which were detected in the VHE band but do not have counterparts at other energy bands. HESS J1303631 was discovered in 2005 (Aharonian et al. 2005) but the nature of this source was unclear until 2012, when a detailed study of the energy-dependent morphology provided evidence of its association with the pulsar PSR J13016305 (Abramowski et al. 2012). With the increase of the energy threshold, a very extended emission region (~ 0.4°× 0.3° at the [ 0.84 − 2 ] TeV band) of VHE γ-rays shrinks towards the position of the pulsar at E> 10 TeV. While at lower energies the peak position of the extended emission region is significantly offset from the position of the pulsar, at energies above 10 TeV the pulsar is coincident with the peak of the γ-ray emission region (Abramowski et al. 2012). Such an energy-dependent morphology is expected for ancient pulsar wind nebulae (PWNe; Abramowski et al. 2012, and references therein). Young electrons located close to the pulsar are not cooled yet and, thus, very energetic. These energetic electrons generate the VHE emission around the pulsar via inverse Compton (IC) scattering on the cosmic microwave background (CMB) photons. Older, cooled down, lower energy electrons might be spread farther away from the pulsar for several reasons (e.g. proper motion of the pulsar that causes older particles to be left behind and/or particle diffusion), but these energetic electrons can still produce γ-rays via IC scattering, however at lower energies than the young electrons.

Table 1

Details of the ATCA data of the HESS J1303631 region analysed in this paper.

The association of HESS J1303631 with the pulsar is further supported by the detection of its X-ray counterpart with XMM-Newton (Abramowski et al. 2012). The size of the X-ray PWN is much smaller than the size of the VHE source, extending 2′ − 3′ from the pulsar position towards the centre of the VHE γ-ray emission region. The much smaller size of the X-ray emitting region can be explained by an effective synchrotron cooling of older electrons to energies too low to generate synchrotron emission in the X-ray energy range and/or owing to the decreasing magnetic field strength in the PWN with time (see e.g. de Jager et al. 2009; Vorster et al. 2013b). The tail-like extension of the X-ray source might be an indication of the proper motion direction of the pulsar triggering speculations about its possible birth place (Abramowski et al. 2012).

An analysis of archival data from the Parkes-MIT-NRAO (PMN) survey at 4.85 GHz (Condon et al. 1993) also revealed a hint of radio emission at the pulsar position with a size comparable to the X-ray emission region (Abramowski et al. 2012). Data analysis showed a ~ 3σ feature with a peak flux density of 30 mJy/beam, which is at the detection limit of the survey and was considered as an upper limit. This hint of a radio counterpart of HESS J1303631 triggered new dedicated observations with the Australian Telescope Compact Array (ATCA), which were conducted in September 2013. Results of these observations are presented in this paper.

Recently, the counterpart of HESS J1303631 was finally detected at GeV energies with Fermi-LAT (Acero et al. 2013). The morphology of the source is consistent with a Gaussian of width 0.45°. The source is contaminated by the emission of the nearby supernova remnant (SNR) Kes 17 (G304.6 + 0.1), but it is clearly seen above 31 GeV. The emission region of the GeV counterpart of HESS J1303631 is, as expected, larger than the TeV source, but the morphology of the emission region is very similar and features an extension in the same direction as the TeV source. However, the size of the GeV source might be slightly overestimated because of the contamination from Kes 17. Kes 17 is the closest known SNR to the pulsar PSR J13016305 located at the angular distance of 37′ (Wu et al. 2011; Gelfand et al. 2013). Assuming a distance to the pulsar of 6.6 kpc (Abramowski et al. 2012), this corresponds to the projected distance between the pulsar and the SNR of 71 pc. This large distance makes the association of the SNR with the pulsar very unlikely as it would require an unrealistically high pulsar velocity of ~6000 km s-1 for the characteristic age of the pulsar of 11 kyr (Manchester et al. 2005).

2. Observations

The ATCA observations of the region surrounding PSR J13016305 were conducted on September 5, 2013. Observations were performed using the CABB receiver with the 1.5A configuration of the array (minimum and maximum baselines of 153 m and 3000 m, respectively) at 5.5 and 7.5 GHz frequencies and centred at α = 13h02m10.00s, δ = − 63°05′34.8′′ (J2000.0), at the angular distance of about 3′ from the pulsar position. The array configuration was chosen to match the resolution of the XMM-Newton observations of ~ 4′′, while at the same time to remain sensitive to structures comparable to the size of the X-ray PWN of ~ 2′. However, the maximum angular scales to which the observations are sensitive are slightly smaller than the size of the X-ray PWN, namely ~ 1.7′ at 5.5 GHz and ~ 1.3′ at 7.5 GHz, respectively.

The observations were carried out in two modes: CFB 1M (no zooms), a bandwidth of 2 GHz with 2048 1-MHz channels in each intermediate frequency (IF) band; and pulsar binning, the same but with the addition of pulsar binning according to the provided ephemerids. The on-source scan time was 656.5 min. Pulsar binning mode was used in observations in order to be able to correctly subtract the pulsar contribution to the total emission to determine the intrinsic emission from HESS J1303631. However, since no significant emission corresponding to HESSJ1303631 was detected (see Sect. 3), the subtraction of the data taken in pulsar binning mode was not performed and thus these data were not used in this study. Primary (flux density) and secondary (phase) calibrators were J1934 − 638 and J1352 − 63, respectively. The flux density of J1934 − 638 is 5.00 Jy at 5.5 GHz and 2.97 Jy at 7.5 GHz. The phase calibrator was observed every ~30 min. The observation recorded all four linear polarization modes. Details of the collected data are listed in Table 1.

In this paper we also considered archival ATCA data obtained during observations of HESS J1303631  centred at α = 13h03m0.400s, δ = − 63°11′11.55′′ (the position of the peak VHE γ-ray emission), and performed in the 1.384 GHz and 2.368 GHz bands. The archival data used in the analysis, taken as part of the Reinfrank et al. project C1557, are presented in Table 1. Only the archival data taken with all six antennas and with observational time longer than 100 min were used in the analysis. The source J1934 − 638 was used as a primary calibrator with flux densities of 14.95 Jy at 1.384 GHz and 11.59 Jy at 2.368 GHz. The sources J1421 − 490 and J1329 − 665 were used for phase calibration. The maximum angular scale to which observations are sensitive is ~ 18′ at 1.384 GHz and ~ 11′ at 2.368 GHz.

3. Data analysis and results

The data reduction and image analysis was performed using the miriad (Sault et al. 1995) and karma (Gooch 1995) packages. The resulting clean primary beam corrected (restricted to the area of primary beam response above 30%, which corresponds to the radial distance of 5.8′) image (Stokes I) at 5.5 GHz is shown in Fig. 1. The pulsar PSR J13016305 is detected at the position α = 13h01m45.678s ± 0.013s, δ = − 63°05′34.85′′ ± 0.20′′. No significant extended emission coincident with the pulsar position was detected. The fitted image root mean square (rms) noise is calculated using the imsad task at the level of 0.011 mJy/beam. The synthesised beam is an ellipse with the major and minor axes of 3.79′′ and 3.65′′, respectively, and the positional angle of − 7.6°.

There is no extended emission coincident with the pulsar position detected at 7.5 GHz as well. The fitted image rms noise, calculated using the imsad task, is at the level of 0.011 mJy/beam. The major and minor axes of the beam are 3.06′′ and 2.90′′, respectively, and the positional angle is − 6.2°.

The analysis of the archival data at 1.384 GHz and 2.368 GHz, which combine all the observations listed in Table 1, does not reveal any significant emission coincident with the pulsar either. The observations at 1.384 GHz, however, reveal a detection of a shell-like structure to the east of the pulsar position that might potentially be an SNR (see discussion in Sect. 4.1). Figure 2 shows the cleaned and primary beam corrected (restricted to the area of primary beam response above 20%, which corresponds to the radial distance of 25.3′) image at 1.384 GHz. The SNR candidate G304.40.2 is positioned within the black circle. The centre of the structure is at α = 13h04m31.1s, δ = − 63°02′08′′. The fitted image rms noise is estimated to be 0.697 mJy/beam. The major and minor axes of the beam are 51.3′′ and 43.8′′, respectively, and the positional angle is 20.9°. The brightest parts of the SNR candidate reach the significance of 13σ. It is difficult to draw any conclusions about a possible emission from the SNR candidate at 2.368 GHz, as these observations are less sensitive to large scale structures, and only a fraction of the SNR candidate is located within the primary beam and the image is distorted by artefacts produced by the strong source MGPS J130237 − 625718.

thumbnail Fig. 1

Radio map of the HESS J1303631 region at 5.5 GHz overlaid with contours of the X-ray PWN (red) as detected with XMM-Newton (Abramowski et al. 2012). X-ray contours show the peak of the emission close to the pulsar with a trail expanding in the eastern direction. The blue circle indicates the position of the pulsar PSR J13016305. The black box determines the region used for the flux density upper limit calculation. New point-like sources are indicated with index numbers. The synthesised beam is determined by an ellipse with the major and minor axes of 3.79′′ and 3.65′′, respectively, and the positional angle of − 7.6°.

Most of the sources detected in the field at 1.384 GHz (Fig. 2), both compact and extended, have counterparts at other radio frequencies (Murphy et al. 2007; Condon et al. 1993; Schinzel et al. 2015). However, most of these sources are not classified, except for MGPS J130422 − 624859 (G304.5 − 0.1), which is identified as a HII region (Misanovic et al. 2002) and MGPS J130037 − 625324, which is coincident with an infrared bubble (Simpson et al. 2012). Extended emission at the south-western edge of the field is most probably related to the very extended source PMN J1259 − 6337 (Wright et al. 1994).

Observations at 1.384 GHz and 2.368 GHz might shed some light on the nature of the unidentified source MGPS J130237 − 625718, which is also visible at other radio wavelengths (e.g. it is detected in the PMN survey). Both 1.384 GHz and 2.368 GHz data reveal a complex morphology with a central source and two lobes, strongly suggesting that MGPS J130237 − 625718 is an active galactic nucleus (AGN). Figure 3 shows the map at 2.368 GHz overlaid with contours indicating significant emission from the source at 1.384 GHz. Only the observations with the best angular resolution (2007-Apr.-24) were used in this analysis. The fitted image rms noise is 0.448 mJy/beam at 1.384 GHz (beam size: major and minor axes of 10.98′′ and 9.78′′, respectively, with the positional angle of 1.7°) and 0.322 mJy/beam at 2.368 GHz (beam size: major and minor axes of 7.15′′ and 6.65′′, respectively, with the positional angle of 3.7°). A fit to the observed emission at 1.384 GHz and 2.368 GHz was carried out using the imfit task with three Gaussian components. The integrated flux density from the central component is 136.8 ± 5.0 mJy and 80.8 ± 3.7 mJy at 1.384 GHz and 2.368 GHz, respectively. The flux densities from the southern and northern components are 260.1 ± 7.3 mJy and 192.5 ± 6.1 mJy at 1.384 GHz and 156.1 ± 7.4 mJy and 99.9 ± 5.5 mJy at 2.368 GHz. Unfortunately, the source is outside the primary beam for observations at 5.5 GHz and 7.5 GHz with much better angular resolution.

thumbnail Fig. 2

Radio map of the HESS J1303631 region at 1.384 GHz. Red contours represent the X-ray PWN as detected with XMM-Newton and blue contours represent HESS J1303631 as detected by H.E.S.S. (Abramowski et al. 2012). The black circle indicates the SNR candidate G304.40.2 and the black cross indicates its centre. MGPS-2 catalogue sources (Murphy et al. 2007) are shown as red circles; PMN survey sources (Condon et al. 1993), excluding those which are coincident with MGPS-2 sources, are shown as green crosses and two ATCA sources (Schinzel et al. 2015) are shown as blue squares. The two new point-like radio sources reported in this paper are indicated with black diamonds. The synthesised beam is determined by an ellipse with the major and minor axes of 51.3′′ and 43.8′′, respectively, and the positional angle of 20.9°.

Table 2

New point-like radio sources detected in these observations.

In the ATCA data presented here, five new point-like sources were detected in the region of HESS J1303631 (Figs. 1 and 2). Their locations and flux densities estimated using the imsad task are collected in Table 2. Only sources with significance above 10σ at both 5.5 GHz and at 1.384 GHz were considered. Each of these sources has one or more potential counterparts in infrared and/or optical catalogues (see e.g. Cutri et al. 2003; Wright et al. 2010; Benjamin et al. 2003; Monet et al. 2003). Two of these sources, J1301 − 6306 and J1300 − 6311, have X-ray counterparts in the XMM-Newton catalogue (Rosen et al. 2016). The source J1304 − 6258 is actually visible in the MGPS-2 (see Fig. 4 right) but is not listed in the catalogue; this is probably because it is very difficult to separate its emission from the extended emission coincident with G304.40.2.

4. Discussion

The new radio observations reported in this paper were triggered by a hint of a signal (~ 3σ) from a feature coincident with the X-ray PWN in the analysis of archival data from the PMN survey at 4.85 GHz (Abramowski et al. 2012). This feature is compatible with the resolution of the survey of ~ 5′ (Condon et al. 1993). These deeper observations at 5.5 GHz and 7.5 GHz with ATCA were optimised for the detection of the putative radio PWN with a size comparable to the size of the X-ray PWN (~ 2′), while providing a resolution that is comparable to that of XMM-Newton. The assumption of the size was motivated by the size of the possible radio feature and by the hypothesis that the sizes of both radio and X-ray PWNe are constrained by the region of high magnetic field around the pulsar. Indeed, three-dimensional magnetohydrodynamic simulations of the Crab nebula (Porth et al. 2014) show that the magnetic field strength close to the termination shock is an order of magnitude higher than the average magnetic field strength in the rest of the PWN. Moreover, in the left-behind relic nebula the magnetic field is expected to be relaxed with the magnetic field strength comparable to the interstellar medium (ISM) magnetic field of about 3 μG. However, neither new observations at 5.5/7.5 GHz nor the analysis of the archival data at 1.384/2.368 GHz show any evidence of extended emission coincident with PSR J13016305 and/or the X-ray PWN. It should be noted that the largest scales to which observations at 5.5/7.5 GHz are sensitive are slightly smaller than the size of X-ray PWN, 1.7′ at 5.5 GHz and 1.3′ at 7.5 GHz, and therefore we cannot rule out the detectability of the putative radio PWN of the size of the X-ray PWN at these frequencies. However, observations at 1.384/2.368 GHz allow us to detect structures with an extension up to ~ 18′ (1.384 GHz) and ~ 11′ (2.368 GHz). Assuming the size of the X-ray PWN as reported in Abramowski et al. (2012), the upper limits on the radio flux density at 1.384 GHz was estimated at the level of three times the rms noise in the region of the X-ray PWN defined by a box (black box in Figs. 1 and 2). The upper limit on the flux density at 1.384 GHz is 2.6 mJy.

thumbnail Fig. 3

Radio map of MGPS J130237 − 625718 at 2.368 GHz overlaid with contours indicating the significance of the emission at 1.384 GHz at the level of 10, 40, and 80 times the rms noise. The synthesised beam is shown with a black ellipse in the bottom left corner.

However, the size of the putative radio PWN might exceed the size of the largest angular structure resolved in the observations or even the size of the primary beam. GeV observations (Acero et al. 2013) exhibiting a large size of the PWN (~ 0.9°) indicate that a large amount of relatively low energy relativistic electrons is spread out to large distances from the pulsar. If the size of the putative radio PWN is constrained by the existence of relativistic electrons, i.e. the magnetic field is strong enough for efficient synchrotron emission across the whole PWN, the radio PWN would be at least as large as the GeV PWN. The same electrons that emit GeV γ-rays via inverse Compton scattering on ambient photon fields are also responsible for the generation of the radio emission via the synchrotron mechanism (see e.g. Vorster et al. 2013b, and references therein).

4.1. G304.40.2: an SNR or not

The extended structure, G304.40.2, detected to the east of the pulsar at 1.384 GHz exhibits a shell-like morphology. Such a morphology naturally suggests that this might be an SNR. G304.40.2 is coincident with PMN J13036259, an extended source detected in the PMN survey at 4.85 GHz (Condon et al. 1993; Wright et al. 1994). The emission from PMN J13036259 was fitted with a two-dimensional asymmetric Gaussian with major and minor axis widths of σx = 14.28′ and σy = 5.88′ with a position angle of 1.5° measured eastwards from the north direction (Griffith & Wright 1993). The total flux density from the source is 235 ± 14 mJy at 4.85 GHz (Wright et al. 1994). Although PMN J13036259 is coincident with the western part of the detected shell-like structure the map also reveals extended radio emission coincident with southern part of the shell (see Fig. 4 middle). Overall, the 4.85 GHz map show a good agreement with the 1.384 GHz morphology of the SNR candidate. To estimate a spectral index of the western part of the SNR candidate, we smoothed the 1.384 GHz map with a Gaussian with a width equal to the angular resolution of the PMN survey of 4.2′ and fitted the region of PMN J13036259 assuming the same extension of the source as obtained for PMN J13036259 using the miriad task imfit. The obtained flux density at 1.384 GHz is 398 ± 3 mJy resulting in a spectral index of α = 0.42 ± 0.051, which is in a perfect agreement with a range of values observed for SNRs, i.e. from ~ 0.2 to ~ 0.8 with a peak at around 0.5 (Reynolds 2011; Green 2014), and very close to a canonical α = 0.5 expected in the diffusive shock acceleration in the case of strong shocks with a compression ratio of 4 (see e.g. Dubner & Giacani 2015, and references therein).

Neither the compact source catalogue nor the SNR catalogue of the second epoch Molonglo Galactic Plane Survey (MGPS-2) at 843 MHz (Murphy et al. 2007; Green et al. 2014) provide any counterpart for the SNR candidate. However, the MGPS-2 image of the G304.40.2 region at 843 MHz reveals a faint extended radio emission coincident with the SNR candidate, roughly following the 1.384 GHz contours in the west and south regions but also exhibiting emission towards north-east (Fig. 4, right). The peak flux density reaches 17.5 mJy/beam, which is well above the sensitivity of the survey of 2 mJy/beam. The image is strongly distorted by artefacts produced by nearby bright sources, which makes it impossible to classify the morphology of the extended emission. This might be the reason why the source did not appear in the list of the SNR candidates detected in the MGPS-2 survey as the search criteria include a condition on a morphology that has to be shell-like or composite (Green et al. 2014).

Infrared and optical observations, for example the Two Micron All-Sky Survey (2MASS) in the H-band (1.65 μm; Skrutskie et al. 2006), do not show any extended emission from the region of G304.40.2. This further suggests that the radio emission from the SNR candidate is most probably non-thermal as expected for SNRs.

thumbnail Fig. 4

Left: radio map of the G304.40.2 region at 1.384 GHz. Middle: radio map of the G304.40.2 region at 4.85 GHz obtained in the PMN survey (Condon et al. 1993). The blue ellipse shows the extension of the PMN survey source PMN J13036259 as presented in the catalogue (Wright et al. 1994). The resolution of this map is ~ 5′ (Condon et al. 1993). Right: radio map of the G304.40.2 region at 843 MHz obtained in the MGPS-2 survey (Murphy et al. 2007). The size of the beam is ~ 45′′ × 50′′. The black circle in all the panels indicates the extension of the shell-like structure and the black cross shows its centre. Middle and right images were obtained using the NASA SkyView online tool (http://skyview.gsfc.nasa.gov).

Available X-ray observations2 of this region show no evidence of the large scale structures coincident with G304.40.2. This, however, does not contradict the SNR hypothesis. The absence of the non-thermal X-ray emission can be explained by a potentially old age of the SNR candidate (see below). Indeed, X-ray synchrotron emission from SNRs requires relatively high shock velocities of 2000 km s-1 to accelerate electrons to energies high enough and these shock velocities are believed to be associated with young SNRs of ~ 1000 yr (see e.g. Vink 2012). Old age might also be a reason for the lack of the X-ray thermal emission. Slowing down of the shock to <200 km s-1 in old SNRs of 10 kyr results in cooling of the post-shock region to temperatures lower than required for the X-ray emission (see e.g. Vink 2012). Some old SNRs, however, feature thermal X-ray emission from their interior that might be due to, for example, interaction with dense cloudlets, which survive the forward shock crossing and slowly evaporate inside the remnant because of saturated thermal conduction (see e.g. White & Long 1991; Sushch et al. 2011). Another possible reason for the lack of the thermal emission might be a low density of the ambient medium around the remnant. How low should this density be is, however, model dependent and depends on whether ionization quilibrium is sustained. Non-equilibirum ionization models coupled to efficient particle acceleration shows that X-ray thermal emission would dominate over synchrotron emission only for relatively high densities of about 1 cm-3 (Ellison et al. 2007). The estimate of the density also depends on the distance to the remnant, which is not clear in this case. Finally, the lack of detection of X-ray emission from the source, both thermal and non-thermal, might be due to insufficient sensitivity or exposure.

At γ-ray energies most of the field of view is covered by HESS J1303631 and it is impossible to distinguish any emission associated with the possible SNR.

However, other possible interpretations of the detected extended emission cannot be firmly ruled out. Owing to the faint emission, the image of G304.40.2 appears to be patchy and this makes it difficult to firmly identify the shell-like morphology. The emission could also be contaminated by unrelated point sources and the image might be distorted by artefacts created by bright sources in the field of view. Since the shell-like morphology of the source is the strongest argument in favour of its SNR nature, we have to treat our conclusions with caution. The observed large scale structure might still be a composition of individual background sources and/or HII regions.

4.2. Possible birth place of PSR J13016305

If G304.40.2 is indeed an SNR it could be a birth place of the pulsar PSR J13016305. This hypothesis is supported by the direction of the X-ray “tail” which roughly points to the position of the SNR candidate (see Fig. 2). The pulsar is located outside the shell-like structure, which means that in case the SNR candidate is indeed the birth place of PSR J13016305, the pulsar has already escaped the remnant and continues to propagate in the ambient medium. The larger size of the PWN does not contradict this hypothesis as the present time nebula could have formed after the pulsar escaped the SNR. While the pulsar is still inside the SNR, its nebula is strongly disrupted by the reverse shock of the remnant (see e.g. Blondin et al. 2001; Vorster et al. 2013a) and at the moment of interaction with the shell it is very small. Moving outside the remnant the pulsar builds up a new nebula that can become very large because of the proper motion of the pulsar, i.e. left behind electrons, and the diffusion of electrons in the medium. Escaping the SNR, the pulsar should also damage the shell of the remnant. Although the emission in the direction of the pulsar is slightly fainter, exhibiting a gap in the shell, there is no clear evidence of distortion. This can be naturally explained if the pulsar is not moving in the projected plane but its velocity has a considerable component that is perpendicular to the projected plane. In this case the distorted part of the shell is facing the observer and is thus not visible. The angular distance between the pulsar and the centre of the SNR candidate of about 19′ corresponds to a projected distance of 36 pc assuming the distance to the pulsar of 6.6 kpc. This corresponds to a projected pulsar velocity of Vp ≃ 3100 km s-1 for the characteristic age of the pulsar of 11 ky (Manchester et al. 2005). This velocity would make PSR J13016305 the fastest known pulsar. The highest pulsar velocity detected so far is ~ 1600 km s-1 (Cordes & Chernoff 1998). However, there are indications that some pulsars might be much faster with velocities comparable to the estimate presented here for PSR J13016305. The estimate of the kick velocity of the possible pulsar IGR J11014 − 6103 is 2400 − 2900 km s-1 (Tomsick et al. 2012), but the source nature cannot be unambiguously proven yet as no pulsations were detected so far. Also, the real age of PSR J13016305 might be higher than the characteristic age if the braking index is lower than 3, which is the case for eight out of nine pulsars for which the braking index has been reliably measured (Lyne et al. 2015; Archibald et al. 2016). In this case the estimate of the pulsar velocity would be lower.

Alternatively, the real age of the pulsar can be estimated as the age of its SNR. In case G304.40.2 is indeed an SNR and the birth place of the pulsar PSR J13016305, its size can provide an estimate for the pulsar age. Assuming a distance to the pulsar of 6.6 kpc (Abramowski et al. 2012), the angular size of the SNR candidate of about 16′ corresponds to 30 pc in diameter. The Sedov solution (Sedov 1959), which describes the hydrodynamical expansion of an SNR in the adiabatic stage of evolution into the homogeneous medium, provides an estimate of the SNR age for a given size of the remnant tage=16(E1051[erg])1/2(nISM1[cm-3])1/2(R15[pc])5/2[kyr],\begin{equation} t_{\mathrm{age}} = 16 \left(\frac{E}{10^{51}\,[\mathrm{erg}]}\right)^{-1/2}\left(\frac{n_{\mathrm{ISM}}}{1\,[\mathrm{cm}^{-3}]}\right)^{1/2}\left(\frac{R}{15\,[\mathrm{pc}]}\right)^{5/2}\,[\mathrm{kyr}], \end{equation}(1)where E is the explosion energy, nISM is the number density of the interstellar medium, and R is the radius of the remnant. This value is somewhat higher than the characteristic age of the pulsar of 11 kyr (Manchester et al. 2005), but does not contradict it if the braking index is lower than 3 (see above). If we assume that the real age of the pulsar is 16 kyr then the projected velocity would be 2100 km s-1, which is still very high. The SNR age estimate, however, is dependent on the ambient medium density, which is often considerably lower than 1 cm-3 even by orders of magnitude, and it is also very sensitive to the estimated physical size of the SNR. Therefore, this estimate of age should be taken with caution.

A pulsar with this high velocity moves supersonically in the interstellar medium creating a bow shock. A bow shock driven through the neutral gas can generate optical emission in the Balmer lines at the forward shock (Bucciantini & Bandiera 2001; Bucciantini 2002) and such emission is already discovered for a few pulsars in Hα. There is no evidence for such emission around PSR J13016305, which however does not necessarily contradict our hypothesis, since a lot of pulsars are believed to be moving supersonically and bow-shock structures are only detected at optical wavelengths for a few of these objects. It has been noted that Hα bowshocks and X-ray tails are rarely seen together (Kargaltsev et al. 2015). Inside the bow-shock pulsar wind particles are accelerated at the termination shock subsequently generating non-thermal synchrotron emission that can be detected in radio and X-rays (Gaensler & Slane 2006). Although no extended radio emission associated with PSR J13016305 was detected, which can simply be because of the potentially large size of the radio nebula, the X-ray nebula exhibits a bow-shock morphology with a tail pointing in the direction of G304.40.2; this strongly suggests that the pulsar is moving with high velocity.

5. Summary

The ATCA observations of the HESS J1303631 region at 5.5 GHz and 7.5 GHz do not reveal any significant extended emission associated with PSR J13016305. Archival 1.384 GHz and 2.368 GHz data also do not show any evidence of a radio counterpart of HESS J1303631. Non-detection of this evolved PWN at radio wavelengths suggests that either the putative radio PWN is comparable to the size of the GeV counterpart and, thus, larger than the largest reliably imaged structure and even the primary beam of radio observations or that the magnetic field is rather low, which is in agreement with the evolved PWN identification because the magnetic field is expected to decrease with time in PWNe. The comparison of the X-ray emission to the TeV emission implies the average magnetic field in the PWN of ~ 0.2 − 2 μG (Abramowski et al. 2012), which however does not exclude the possibility of enhanced magnetic field around the pulsar.

Archival 1.384 GHz observations reveal a detection of an extended structure centred at the angular distance of 19′ from the pulsar. This extended structure might be a SNR and a potential birth place of the pulsar. If this is the case then the projected velocity of the pulsar would be 3100 km s-1 assuming the characteristic age of the pulsar. This would make PSR J13016305 the fastest known pulsar. However, the uncertainty of the true age of the pulsar can significantly change this estimate.


1

Flux density, Sν, scales with frequency, ν, as Sννα.

2

All observations of the region of interest available at NASA’s HEASARC archive (https://heasarc.gsfc.nasa.gov/docs/archive.html) were examined including XMM-Newton, Swift, ROSAT, Chandra, and ASCA.

Acknowledgments

We would like to thank the anonymous referee for valuable comments that strongly improved the paper. We would like to acknowledge the help of Marek Jamrozy (Astronomical Observatory of the Jagiellonian University of Kraków, Poland) and Michael Bietenholz (Hartebeesthoek Radio Observatory, South Africa; York University, Toronto, Canada) on understanding certain aspects of the ATCA analysis performed in the paper. The Australia Telescope Compact Array is part of the Australia Telescope National Facility, which is funded by the Australian Government for operation as a National Facility managed by CSIRO. This paper includes archived data obtained through the Australia Telescope Online Archive (http://atoa.atnf.csiro.au).

References

  1. Abramowski, A., Acero, F., et al. (H.E.S.S. Collaboration) 2012, A&A, 548, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Acero, F., Ackermann, M., Ajello, M., et al. 2013, ApJ, 773, 77 [CrossRef] [Google Scholar]
  3. Aharonian, F., Akhperjanian, A. G., Aye, K.-M., et al. 2005, A&A, 439, 1013 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Archibald, R. F., Gotthelf, E. V., Ferdman, R. D., et al. 2016, ApJ, 819, L16 [NASA ADS] [CrossRef] [Google Scholar]
  5. Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953 [NASA ADS] [CrossRef] [Google Scholar]
  6. Blondin, J. M., Chevalier, R. A., & Frierson, D. M. 2001, ApJ, 563, 806 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bucciantini, N. 2002, A&A, 393, 629 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bucciantini, N., & Bandiera, R. 2001, A&A, 375, 1032 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Condon, J. J., Griffith, M. R., & Wright, A. E. 1993, AJ, 106, 1095 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cordes, J. M., & Chernoff, D. F. 1998, ApJ, 505, 315 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  12. de Jager, O. C., Ferreira, S. E. S., Djannati-Ataï, A., et al. 2009, ArXiv e-prints [arXiv:0906.2644] [Google Scholar]
  13. Dubner, G., & Giacani, E. 2015, A&ARv, 23, 3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Ellison, D. C., Patnaude, D. J., Slane, P., Blasi, P., & Gabici, S. 2007, ApJ, 661, 879 [NASA ADS] [CrossRef] [Google Scholar]
  15. Gaensler, B. M., & Slane, P. O. 2006, ARA&A, 44, 17 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gelfand, J. D., Castro, D., Slane, P. O., et al. 2013, ApJ, 777, 148 [NASA ADS] [CrossRef] [Google Scholar]
  17. Gooch, R. 1995, in Astronomical Data Analysis Software and Systems IV, eds. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, ASP Conf. Ser., 77, 144 [Google Scholar]
  18. Green, D. A. 2014, Bull. Astron. Soc. India, 42, 47 [NASA ADS] [Google Scholar]
  19. Green, A. J., Reeves, S. N., & Murphy, T. 2014, PASA, 31, e042 [NASA ADS] [CrossRef] [Google Scholar]
  20. Griffith, M. R., & Wright, A. E. 1993, AJ, 105, 1666 [NASA ADS] [CrossRef] [Google Scholar]
  21. Kargaltsev, O., Cerutti, B., Lyubarsky, Y., & Striani, E. 2015, Space Sci. Rev., 191, 391 [NASA ADS] [CrossRef] [Google Scholar]
  22. Lyne, A. G., Jordan, C. A., Graham-Smith, F., et al. 2015, MNRAS, 446, 857 [NASA ADS] [CrossRef] [Google Scholar]
  23. Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993 [NASA ADS] [CrossRef] [Google Scholar]
  24. Misanovic, Z., Cram, L., & Green, A. 2002, MNRAS, 335, 114 [NASA ADS] [CrossRef] [Google Scholar]
  25. Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 [NASA ADS] [CrossRef] [Google Scholar]
  26. Murphy, T., Mauch, T., Green, A., et al. 2007, MNRAS, 382, 382 [NASA ADS] [CrossRef] [Google Scholar]
  27. Porth, O., Komissarov, S. S., & Keppens, R. 2014, MNRAS, 438, 278 [NASA ADS] [CrossRef] [Google Scholar]
  28. Reynolds, S. P. 2011, Ap&SS, 336, 257 [NASA ADS] [CrossRef] [Google Scholar]
  29. Rosen, S. R., Webb, N. A., Watson, M. G., et al. 2016, A&A, 590, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data Analysis Software and Systems IV, eds. R. A. Shaw, H. E. Payne, & J. J. E. Hayes, ASP Conf. Ser., 77, 433 [Google Scholar]
  31. Schinzel, F. K., Petrov, L., Taylor, G. B., et al. 2015, ApJS, 217, 4 [NASA ADS] [CrossRef] [Google Scholar]
  32. Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics (New York: Academic Press) [Google Scholar]
  33. Simpson, R. J., Povich, M. S., Kendrew, S., et al. 2012, MNRAS, 424, 2442 [NASA ADS] [CrossRef] [Google Scholar]
  34. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  35. Sushch, I., Hnatyk, B., & Neronov, A. 2011, A&A, 525, A154 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Tomsick, J. A., Bodaghee, A., Rodriguez, J., et al. 2012, ApJ, 750, L39 [NASA ADS] [CrossRef] [Google Scholar]
  37. Vink, J. 2012, A&ARv, 20, 49 [Google Scholar]
  38. Vorster, M. J., Ferreira, S. E. S., de Jager, O. C., & Djannati-Ataï, A. 2013a, A&A, 551, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Vorster, M. J., Tibolla, O., Ferreira, S. E. S., & Kaufmann, S. 2013b, ApJ, 773, 139 [NASA ADS] [CrossRef] [Google Scholar]
  40. White, R. L., & Long, K. S. 1991, ApJ, 373, 543 [NASA ADS] [CrossRef] [Google Scholar]
  41. Wright, A. E., Griffith, M. R., Burke, B. F., & Ekers, R. D. 1994, ApJS, 91, 111 [NASA ADS] [CrossRef] [Google Scholar]
  42. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
  43. Wu, J. H. K., Wu, E. M. H., Hui, C. Y., et al. 2011, ApJ, 740, L12 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Details of the ATCA data of the HESS J1303631 region analysed in this paper.

Table 2

New point-like radio sources detected in these observations.

All Figures

thumbnail Fig. 1

Radio map of the HESS J1303631 region at 5.5 GHz overlaid with contours of the X-ray PWN (red) as detected with XMM-Newton (Abramowski et al. 2012). X-ray contours show the peak of the emission close to the pulsar with a trail expanding in the eastern direction. The blue circle indicates the position of the pulsar PSR J13016305. The black box determines the region used for the flux density upper limit calculation. New point-like sources are indicated with index numbers. The synthesised beam is determined by an ellipse with the major and minor axes of 3.79′′ and 3.65′′, respectively, and the positional angle of − 7.6°.

In the text
thumbnail Fig. 2

Radio map of the HESS J1303631 region at 1.384 GHz. Red contours represent the X-ray PWN as detected with XMM-Newton and blue contours represent HESS J1303631 as detected by H.E.S.S. (Abramowski et al. 2012). The black circle indicates the SNR candidate G304.40.2 and the black cross indicates its centre. MGPS-2 catalogue sources (Murphy et al. 2007) are shown as red circles; PMN survey sources (Condon et al. 1993), excluding those which are coincident with MGPS-2 sources, are shown as green crosses and two ATCA sources (Schinzel et al. 2015) are shown as blue squares. The two new point-like radio sources reported in this paper are indicated with black diamonds. The synthesised beam is determined by an ellipse with the major and minor axes of 51.3′′ and 43.8′′, respectively, and the positional angle of 20.9°.

In the text
thumbnail Fig. 3

Radio map of MGPS J130237 − 625718 at 2.368 GHz overlaid with contours indicating the significance of the emission at 1.384 GHz at the level of 10, 40, and 80 times the rms noise. The synthesised beam is shown with a black ellipse in the bottom left corner.

In the text
thumbnail Fig. 4

Left: radio map of the G304.40.2 region at 1.384 GHz. Middle: radio map of the G304.40.2 region at 4.85 GHz obtained in the PMN survey (Condon et al. 1993). The blue ellipse shows the extension of the PMN survey source PMN J13036259 as presented in the catalogue (Wright et al. 1994). The resolution of this map is ~ 5′ (Condon et al. 1993). Right: radio map of the G304.40.2 region at 843 MHz obtained in the MGPS-2 survey (Murphy et al. 2007). The size of the beam is ~ 45′′ × 50′′. The black circle in all the panels indicates the extension of the shell-like structure and the black cross shows its centre. Middle and right images were obtained using the NASA SkyView online tool (http://skyview.gsfc.nasa.gov).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.