A new study towards PSR J1826-1334 and PSR J1826-1256in the region of HESS J1825-137 and HESS J1826-130

# A new study towards PSR J1826−1334 and PSR J1826−1256 in the region of HESS J1825−137 and HESS J1826−130

L. Duvidovich CONICET - Universidad de Buenos Aires, Instituto de Astronomía y Física del Espacio (IAFE), Buenos Aires, Argentina.
Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina.
E. Giacani CONICET - Universidad de Buenos Aires, Instituto de Astronomía y Física del Espacio (IAFE), Buenos Aires, Argentina.
Facultad de Arquitectura Diseño y Urbanismo, Universidad de Buenos Aires, Buenos Aires, Argentina.
G. Castelletti CONICET - Universidad de Buenos Aires, Instituto de Astronomía y Física del Espacio (IAFE), Buenos Aires, Argentina.
Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina.
A. Petriella CONICET - Universidad de Buenos Aires, Instituto de Astronomía y Física del Espacio (IAFE), Buenos Aires, Argentina.
Ciclo BÃ¡sico ComÃºn, Universidad de Buenos Aires, Buenos Aires, Argentina.
L. Supan CONICET - Universidad de Buenos Aires, Instituto de Astronomía y Física del Espacio (IAFE), Buenos Aires, Argentina.
Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina.
Received 7 November 2018 / Accepted 21 January 2019
###### Key Words.:
ISM: individual objects: HESS J1825$-$137, HESS J1826$-$130-rays: ISM – X-rays: ISM – Radio continuum: ISM
offprints: L. Duvidovich
###### Abstract

Context:

Aims:The goal of this paper is to detect synchrotron emission from the relic electrons of the crushed pulsar wind nebula (PWN) HESS J1825137 and to investigate the origin of the -ray emission from HESS J1826130.

Methods:The study of HESS J1825137 was carried out on the basis of new radio observations centred at the position of PSR J18261334 performed with the Karl G. Jansky Very Large Array at 1.4 GHz in configurations B and C. To investigate the nature of HESS J1826130, we reprocessed unpublished archival data obtained with XMM-Newton.

Results:The new radio continuum image towards PSR J18261334 reveals a bright radio source, with the pulsar located in its centre, which suggests that this feature could be the radio counterpart of the compact component of the PWN detected at high energy. The new 1.4 GHz radio data do not reveal emission with an extension comparable with that observed in -rays for the HESS J1825137 source. On the other hand, the XMM-Newton study of the region including PSR J18261256 reveals an elongated non-thermal X-ray emitting nebula with the pulsar located in the northern border and a tail towards the peak of the very high energy source. The spectrum is characterized by a power law with a photon index going from 1.6 around the pulsar to 2.7 in the borders of the nebula, a behaviour consistent with synchrotron cooling of electrons. From our X-ray analysis we propose that HESS J1826130 is likely produced by the PWN powered by PSR J18261256 via the inverse Compton mechanism.

Conclusions:

## 1 Introduction

The H.E.S.S. Galactic Plane Survey (HGPS) has recently completed nine years of continuum observations in the 250-65 longitude range for latitudes b 3 (H. E. S. S. Collaboration et al., 2018). Observations towards the region containing the known pulsar wind nebula (PWN) HESS J1825137 show a new source named HESS J1826130, which emerges towards the north of HESS J1825137 and whose origin is still unknown (see Fig. 2).

HESS J1825137 is one of the brightest and most extensive PWN detected at very high energies (VHEs) (H. E. S. S. Collaboration et al., 2019). It has been proposed that this TeV source is powered by the high spin-down luminosity ( = 2.8  erg s) pulsar PSR~J1826$-$1334 (Manchester et al., 2005). This pulsar is located at a distance of 4 kpc based on dispersion measurements (Cordes & Lazio, 2002) and its spin-down age is 20 kyr. The pulsar is about 10 offset from the peak of HESS J1825137. This TeV source was the first to exhibit an energy dependent morphology characterized by a spectral steeping with increasing distance from the pulsar. In the GeV range, observations made with Fermi-LAT showed both morphological similarity and spectrum continuity with the H.E.S.S. data (Grondin et al., 2011).

XMM-Newton observations towards PSR J18261334 show a bright and elongated core that is 30 in size with a hard photon index (  ), which is embedded in a fainter diffuse component of emission that extends with a softer photon index (  ) and is located on the southern side of the pulsar (Gaensler et al., 2003). This spectral behaviour is remarkably similar to that observed in -rays, except that the spatial scales are different (few arcmin in X-rays, fraction of a degree in -rays). Inside the compact nebula, Chandra data reveal a small component (7) elongated in the NE-SW direction. The extended X-ray emission is detected in these observations up to at least south of the pulsar (Pavlov et al., 2008a).

The offset and asymmetrical nature of the emission in the X-ray and -ray ranges have been explained by Gaensler et al. (2003) and Aharonian et al. (2006) using the scenario known as “crushed PWN” (Blondin et al., 2001). In this framework, the supernova remnant (SNR) blastwave expands into an inhomogeneous interstellar medium (ISM). Consequently, an asymmetric reverse shock is generated, which returns from higher density regions first, reaching one side of the PWN sooner than the other side, then pushing the X-ray and TeV emission mainly towards regions of lower ambient density.

In the radio band, the region around PSR J18261334 has been explored at 8.3 GHz by Frail & Scharringhausen (1997) and at 1.4 GHz by Gaensler et al. (2000) using the Very Large Array (VLA). However, in both cases the observations failed to detect nebular emission down to the respective sensitivity limits (rms noise levels of 0.2 and 0.5 mJy beam at 8.3 and 1.4 GHz, respectively).

Regarding HESS J1826130, it is a new unidentified extended TeV source previously hidden in the emission from the bright nearby PWN HESS J1825137 (H. E. S. S. Collaboration et al., 2018). Fitting the excess map with a 2D symmetric Gaussian, Angüner et al. (2017) estimated the best fit-position of HESS J1826130 to be R.A.=18260.2, Dec.=1302 (J2000) with an extension of 0.17 .

Based on their spatial coincidence, HESS J1826130 has been associated with the PWN~G18.5$-$0.4 (also known as ‘Eel’ PWN), which is a diffuse and weak non-thermal X-ray nebula detected with Chandra (Roberts et al., 2007). This PWN is probably powered by the pulsar PSR J18261256, which lies within the VHE emission region offset from its centroid, and now is notable for being one of the brightest radio-quiet -ray pulsars (Acero et al., 2015). Nevertheless, within the extension of HESS J1826130 there are other potential sources of energetic particles capable to power the TeV emission, such as radio SNRs and HII regions (Paron et al., 2013), and this fact hampers a decisive conclusion on the origin of the -rays.

In this paper we report on new high-resolution and high-sensitivity radio observations carried out with the Karl G. Jansky Very Large Array (JVLA)333The Very Large Array of the National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. towards PSR J18261334 with the aim to detect the radio counterpart of the PWN detected at high and VHEs. In addition, we present a new X-ray study carried out with archival data acquired with the XMM-Newton satellite of the PWN G18.50.4 to investigate its connection with HESS J1826130 and shed light on the origin of this TeV source.

## 2 Observations and data reduction

We carried out observations centred on PSR J18261334 using the B and C configurations of the JVLA.444Further information on the JVLA array configurations is reported at https://science.nrao.edu/facilities/vla/docs/manuals/propvla/array_configs. The data were taken at L band using the wide-band 1 GHz receiver system centred at 1.4 GHz, which consists in 16 spectral windows with a bandwidth of 64 MHz each, spread into 64 channels. These observations are under the project 12A-166, which was intended to carry out deep full-synthesis imaging of intensity of the nebular radio emission from the vicinity of the scientific target source.

The observing details are summarized in Table 1. All the data were reduced using the capabilities within the CASA software infrastructure.

After an initial data inspection and editing, the calibration was done according a standard procedure separately for each of the four observing dates. This includes determinations of antenna position and delay corrections, as well as complex antenna gains and bandpass solutions. For the four datasets observations of the primary flux density calibrator 3C~286 were used, in reference to a model of the flux distribution of the calibrator source, to calculate the antenna-based gain (amplitude and phase) corrections and to derive the flux of the phase calibrator J1834$-$1237. The source 3C 286 was also used to correct the bandpass.

After the calibration was applied, the data of the scientific target PSR J18261334 was split out from all of the observations listed in Table 1 and the individual files were imaged together using the capability incorporated in the CASA task tclean. Two rounds of phase-only self-calibration were performed to improve the dynamic range. It was found that a further amplitude calibration iteration was not necessary as it showed very little additional improvements on the image results. The self-cal process was made with a multi-scale technique with three different scale sizes and a robust parameter of 0.5. A wide-field imaging technique based on the W-projection algorithm was also applied during the imaging process to take account the non-coplanar nature of the JVLA baselines as a function of the distance from the phase centre.

The synthesized beam of the final L-band image after combining the four datasets described in Table 1 and including primary beam corrections is 9.246.43 with a position angle PA=0.05. The sensitivity achieved in the L-band image, without primary beam corrections applied, is 0.24 mJy beam.

### 2.2 X-Ray data

We reprocessed unpublished archival XMM-Newton data (obs. ID #0744420101) towards PSR~J1826$-$1256 to analyse the properties of the extended X-ray emission. The EPIC-pn camera was operated in the small-window mode, which partially covers the extent of the PWN around PSR J18261256. On the other hand, EPIC-MOS1 and MOS2 cameras were set in the full-frame mode and hence mapped the full extension of the nebula. We employed the two EPIC-MOS cameras for the analysis of the X-ray emission of the PWN and the three cameras for the spectral study of the pulsar.

We reduced the observation and extracted the scientific products using the software packages SAS 16.1.0 and Heasoft 6.22.1, following standard procedures described in the SAS Sciences Threads. After obtaining calibrated event files using the latest current calibration files (CCF), we constructed a light curve in the high-energy band ( keV) to filter periods of high count-rate. We applied FLAG = 0 and PATTERN and 4 for MOS and pn cameras, respectively. We obtained filtered event files with 124, 128, and 91 ks for the MOS1, MOS2, and pn cameras, respectively. We constructed individual images from MOS1 and MOS2 cameras in the 1-7 keV band, the energy range where we detect emission, and combined these into a mosaicked and exposure-corrected image with a spatial scale of 3/pix. The final image was smoothed by convolving it with a 2D Gaussian function with pixel.

## 3 Results and discussion

### 3.1 Hess J1825−137

The main panel of Fig. 2 shows the new JVLA radio image at 1.4 GHz centred on the position of PSR J18261334. As the primary beam of the interferometer at the mentioned frequency is about 30, the data only cover the northern half of the source HESS J1825137. Especially prominent is the emission from two objects in the surveyed field. One of these sources, with its brightness peak at R.A.=18261.7, Dec.=133814.6 (J2000), is in positional coincidence with the HII region G017.9280.677 and most likely represents thermal emission associated with ionized gas (Anderson et al., 2011). The catalogued kinematical distance to this thermal source is 13 kpc, based on recombination-line emission detected at 39 km s with the NRAO Green Bank 140-foot telescope Anderson et al. 2011, see also Roman-Duval et al. 2009. The emission centred at R.A.=182537.5, Dec.=133641 (J2000) is catalogued as NVSS  J182538133637 in the 1.4 GHz NRAO VLA Sky Survey (Condon et al., 1998) and referred to as source “C” in Pavlov et al. (2008a). To the best of our knowledge, the literature does not contain a counterpart that helps to discern its nature.

A zoomed view of the region surrounding PSR J18261334 is shown in the upper left corner of Fig 2. This image shows, for the first time, emission from a bright region around the pulsar with a size of about 20, in good agreement with the size of the X-ray compact component of the PWN reported by Gaensler et al. (2003). We estimated the total flux density of this structure, detected at 4 above the noise level, in 4 mJy. The brightest central portion of the compact radio component, which is comparable in size to the angular resolution of our observations, harbours the pulsar. Its contribution to the flux calculated over all the compact component is about 3 mJy. We highlight that this result is largely consistent with the flux density determination for PSR J18261334 reported by Kijak et al. (2007), on the basis of the total on-pulse energy measured over the pulse period around 1 GHz using the Giant Metrewave Radio Telescope see also the compatible estimates presented in Jankowski et al. (2018) and references therein.

The enlargement included in the new image at 1.4 GHz (Fig. 2) also shows, at 2 above the noise, a feature extending to the north of the pulsar position. Even though the detection occurs at a low level, the internal structure of this faint feature suggests that it could be real. If it were, we are likely observing part of the X-ray PWN. To confirm this, further observations with higher sensitivity are required. In addition, further polarization and spectrum studies are needed since a high degree of polarization () and a flat radio spectral index (, S , where is the flux density, the frequency and the radio spectral index) are two unmistakable properties of PWNe in the radio band (Gaensler & Slane, 2006).

The non-detection in the radio band of a PWN detected at high energies is not an unusual result. Indeed, there are a large number of pulsars for which the search of radio PWNe yielded negative results (e.g. Gaensler et al., 2000; Giacani et al., 2014; Castelletti et al., 2016). In most of these cases, the failure to detect the nebular emission has been explained as a consequence of either a high magnetic field of the pulsar that inhibits the production of synchrotron radiation at longer wavelengths, or severe adiabatic losses that occur in young and energetic pulsars producing under luminous radio PWNe if they evolve in a very low ambient density (0.003 cm). Nevertheless, neither of the two conditions are satisfied for PSR J18261334. On one hand, taking into account the parameters of the pulsar ( = 101 ms and = 7.52 10 s s, Clifton et al. (1992)), we estimated its magnetic field strength in about 2.8  10 G, a value typical of young pulsars in SNRs. On the other side, the requirement of an ambient of very low density is not consistent with the presence of several dense molecular regions in the vicinity of the pulsar at its kinematic distance (Voisin et al., 2016).

Even though the two conditions invoked to account for the lack of an associated radio nebula fail to explain the case for PSR J18261334, we can quantify our non-detection. To do this, we characterized the luminosity integrated between 10 Hz and 10 Hz of the putative radio PWN by L erg s, where is the associated spin-down luminosity of the pulsar and is the efficiency of conversion of in radio emission. Assuming that typical values for the efficiency (10) and the spectral index (=0.3) measured in radio PWNe (Gaensler & Slane, 2006) are also valid for the case of PSR J18261334, we obtained an upper limit for the flux density at 1.4 GHz of 400 mJy; this upper limit is calculated using the relation  mJy, where is the distance to the pulsar in kpc and is in units of 10 erg. If the extension of the supposed PWN were at least the same as that of the X-ray nebula observed with XMM-Newton satellite (ellipsoid 5 4 in size), then the limit of the detection would be 0.3 mJy beam, i.e. of the order of the rms noise level of our radio observations. Morever, in the case in which the size of the radio PWN were comparable to that observed in rays, the corresponding limit for its detection in radio would be even one order of magnitude lower.

### 3.2 Hess J1826−130

#### 3.2.1 X-ray imaging and spectral analysis

Figure 3a shows the X-ray emission towards PSR J18261256 in the 1-7 keV energy band with overlaid contours of the TeV emission from HESS J1826130, for reference, while panel b of Fig. 3 indicates a brightness profile. This profile was constructed from a mosaicked image using data from the MOS1 and MOS2 cameras with a spatial scale of . No exposure correction and smoothing was applied. We defined rectangular boxes of placed along the most intense X-ray emission. For the background, we used rectangular regions of orientated in the same way as the source boxes and placed outside the diffuse emission. We defined the box containing the pulsar J18261256 as the reference position to measure boxes distances.

The large-exposure XMM-Newton observation allows us to identify some morphological features not detected in previous imaging studies performed with Chandra. Indeed, the asymmetric and non-uniform distribution of the X-ray emitting gas around PSR J18261256 is shown in Fig. 3. Bright emission is observed ahead of the pulsar up to a distance of (Fig. 3b). We notice, however, that this emission may be due to projection effects. An enlargement showing the detailed X-ray morphology is presented in Fig. 4a. The bulk of the emission comes from an elongated feature with an elliptical shape of (major minor-axis), which appears brighter around the pulsar and extends towards the southwest in the direction of the TeV source. We refer to this feature as the “bar” and we associate it with the “tail” identified by Roberts et al. (2007) in the Chandra image. These authors suggested that the shape of the X-ray PWN could be caused by the motion of the pulsar which leaves behind a high-energy trailing nebula. The bar is surrounded by faint and diffuse emission that is more prominent towards the southeast. PSR J18261256 appears slightly displaced from the apparent axis of the bar (indicated with a green line in Fig.4a).

Several point-source candidates are seen in the field in the Chandra Source Catalog (CSC; Evans et al. 2010) and in the XMM-Newton Serendipitous Source Catalog (3XMM DR8; Rosen et al. 2016). We identified the brightest sources that appear over the bar (indicated with red circles in Fig.4a) and excluded these from the spectral analysis. To obtain the overall properties of the PWN, we extracted the spectrum of the bar from the ellipse shown in Fig. 4a. The background was chosen from a circular region free of diffuse emission and excluding the point sources over it. We obtained individual spectra for MOS1 and MOS2 with minimum of 30 cts/bin. Inspection and fitting of the spectra was done using Xspec 12.9.1. The spectra are hard and show a continuum with no hint of emission lines (see Fig. 4b). We simultaneously fitted the MOS1 and MOS2 spectra in the 1-7 keV energy band with an absorbed power-law model (). The results are shown in Table 2. The reported flux was calculated in the 0.5-8.0 keV energy band and this flux is used in Sect. 3.2.2 to analyse the possible connection between the X-ray PWN and HESS J1826130.

Taking advantage of the high number of counts of the XMM-Newton observations, we performed a spatially resolved spectral analysis to look for variations of the photon index along the diffuse emission. We extracted individual spectra from nine rectangular boxes located along the bar (see Fig. 4a) from MOS1 and MOS2 cameras excluding the point sources. We defined the distance among boxes as the separation between their centres and the centre of box 1, which contains the pulsar J18261256. This is indicated in Fig. 4a with a green arrow. Spectra were binned to a minimum of 30 cts/bin. For each box, the two spectra are simultaneously fitted with a power-law model in the 1-7 keV energy band, keeping the hydrogen column density frozen to the best-fit value obtained for the whole diffuse emission, namely  cm. The results are summarized in Table 3 and plotted in Fig. 5. The top panel in Fig. 5 shows the variation of the photon index variation with respect to . The softening of the spectrum is remarkable, where increases from 1.6 around the pulsar (boxes 0-2) up to 2.7 in the outer regions (boxes 7 and 8). This behaviour has been observed in several Galactic PWNe, such as G0.9+0.1 (Porquet et al., 2003), 3C 58 (Bocchino et al., 2001), and G21.50.9 (Safi-Harb et al., 2001), among others. The spectral softening is a consequence of the underlying synchrotron radiation mechanism responsible for the non-thermal emission in the keV band. Highly energetic electrons cool faster than low energetic electrons after having travelled a short distance from the powering pulsar while they diffuse throughout the pulsar wind. On the other hand, the surface brightness SB, which is defined as the ratio of the absorption corrected flux in the 1-7 keV band and the area of the corresponding region, decreases with increasing distance to the powering pulsar, as shown in the bottom panel of Fig. 5. This is a consequence of the synchrotron radiation mechanism responsible for the X-ray emission from PWNe (Holler et al., 2012).

For the spectral study of the pulsar PSR J18261256 we used the three cameras (MOS1, MOS2, and pn) and extracted the corresponding spectra from a circular region of radius . Spectra were binned to a minimum of 15 cts/bin. For the background, we defined an annular region with inner and outer radii of and , respectively. X-ray emission from pulsars may correspond to black-body (BB) emission from the neutron star surface and/or to non-thermal synchrotron emission from the magnetosphere. To look for non-thermal emission from the pulsar, we fitted the spectra with an absorbed power-law model. Letting to vary freely or freezing it to the best-fit value of the diffuse emission ( cm) gives statistically acceptable fits () and model parameters that agree within the uncertainties (see Table 2). Fitting the spectra with an absorbed BB () and fixed to cm yields a statistically poorer fit (), which fails to model the low-energy part of the spectra ( keV). Letting to vary freely results in a better fit () but a hydrogen column density of 0.2 cm is obtained in this case, which is about an order of magnitude lower than the column density for the diffuse emission. If PSR J18261256 is powering the PWN, we would expect a similar ISM absorption. Moreover, the obtained BB temperature (1.6 keV) is higher than expected for young pulsars, which usually have BB temperatures  keV (Mereghetti, 2011). We therefore conclude that the model which better fits the X-ray emission from PSR J18261256 is a power law with the hydrogen column density similar to that obtained towards the PWN, as expected if they are associated sources.

#### 3.2.2 Origin of TeV emission

In a search for counterparts to HESS J1826130, the X-ray PWN driven by the PSR J18261256 appears to be the most plausible candidate. In a leptonic scenario the TeV emission is expected to arise from inverse Compton (IC) scattering between the ambient low energetic photons and the same population of electron producing synchrotron radiation in the keV band. In our analysis of the spectral behaviour (see Sect. 3.2.1), we found that the photon index of the extended X-ray PWN is , close to the value of the TeV spectrum of HESS J1826130 (H. E. S. S. Collaboration et al., 2018). This fact suggests a scenario in which the X-ray and TeV emissions arise from the same population of electrons (Pavlov et al., 2008b). Based on our X-ray results and assuming that the softening of the spectrum with the distance from the pulsar is due to cooling effects, we estimated the magnetic field in the nebula from the relation  kyr, where is the magnetic field expressed in and being the energy of photons produced in the synchrotron radiation. As a consequence of the radiative cooling, the synchrotron cooling time should be smaller than the age of the PWN, which we roughly approximated to the pulsar’s characteristic age  kyr. The condition , gives for  keV (the peak of the X-ray spectrum shown in Fig. 4a). The corresponding energy of the electrons is  TeV. For  keV and  G, we constrained the energy of the synchrotron emitting electrons to be  TeV.

The population of relativistic electrons analysed in the paragraph above, with energy , can also interact with background photons of energy via IC scattering producing -ray photons in the TeV band, revealed as the VHE counterpart to the X-ray PWN. Hence, we estimated the energy of the TeV photons through the relation  TeV, where . If the cosmic microwave background is the main source of background photons with  K, then  eV. Using  TeV, we found that the energy of the photons produced by the IC scattering is  TeV. This analysis shows that the VHE counterpart to the X-ray nebula around PSR J18261256 is expected to radiate IC photons with energies up to a few tens of TeV. This result is largely compatible with the detection of rays in the 0.5-40 TeV range from HESS J1826130 (Angüner et al., 2017).

Another approach to analyse the possible connection between PSR J18261256, its PWN, and HESS J1826130, is to compare the high-energy properties of the system with those observed in others PWNe. To do this, we used the luminosities and to estimate the efficiency of pulsar’s energy conversion in the keV and TeV bands as and , where and are the luminosities of the PWN in the 0.5-8.0 keV and 1-10 TeV energy bands, respectively. From Table 2, the flux of the PWN in the 0.5-8.0 keV band is  erg s cm. This translates into a luminosity  erg s, where is the distance in kpc. Using the -ray flux reported by Angüner et al. (2017), the luminosity in the 1-10 TeV band is  erg s. From the obtained luminosities, together with the pulsar spin-down power erg s, we calculated the efficiency of pulsar’s energy conversion and as a function of the distance. In the case of PSR J18261256 its distance is not well known. Wang (2011) suggested a distance to the pulsar of about 1.2 kpc, based on the correlation found between the efficiency in the GeV band and other pulsar parameters in a sample of -ray pulsars detected by Fermi-LAT. However, from a molecular study of a wide region containing the source HESS J1826 130, Voisin et al. (2016) reported the presence of molecular gas in spatial coincidence with it at kinematical velocities between 60 and 80 km s that correspond to intermediate near/far distances of about 4.6/11.4 kpc. Besides, the authors pointed out that a distance of 1.2 kpc to the pulsar implies a TeV -ray efficiency lower than the typical values. Thus, by adopting in our calculations the distances reported by Voisin et al. (2016) for the PSR-PWN-HESS source system, we compared their properties with the population of Galactic PWNe and listed them in Table 4. Photon indexes , in both keV and TeV bands, are in good agreement with those measured for other PWNe, and the obtained efficiencies and show that PSR J18261256 has a spin-down energy sufficient to power the observed nebulae in X- and -rays (Kargaltsev et al., 2013; H. E. S. S. Collaboration et al., 2019).

We also mention that a number of HII regions, OB stars, and high-mass star-forming regions are located in the vicinity of the TeV source at a distance of 4 kpc (Paron et al., 2013). These objects are able to accelerate particles to multi-TeV energy via hadronic interactions with dense material (e.g, Araudo et al., 2008; Bosch-Ramon et al., 2010). At the distance of 4 kpc, Voisin et al. (2016) reported the presence of molecular gas overlapping the TeV source, thus the contribution to the -ray emission via a hadronic process cannot be completely ruled out.

## 4 Summary

In this paper we presented the highest angular resolution and sensitivity radio continuum emission image from JVLA observatory at 1.4 GHz towards the pulsar PSR J18261334. We sought to search for the radio counterpart to the PWN HESS J825137, which has been also detected in the X-ray band. The new radio data allowed us to detect an extended bright radio source with the pulsar located in its centre, in positional coincidence, and with a size similar to that of the compact component of the X-ray PWN. Although with the available radio data it is not possible to go deeper in the analysis of the origin of the observed emission, we suggest that it could be the radio counterpart of the PWN detected in the X- and -ray spectral regimes. Regarding the diffuse emission of the X-ray PWN, its radio counterpart was not detected in the 1-2 GHz band with an rms noise of 0.24 mJy beam.

We also reported the analysis of unpublished archival X-ray observations performed with XMM-Newton towards PSR J18261256. The new image revealed an elongated nebula, brighter around the pulsar with the emission in the direction of the peak of the TeV source HESS J1826130. The spectral analysis demonstrated a non-thermal origin for the X-ray emission and a photon index softening with the distance to the pulsar. This study shows that the most plausible origin for HESS J1826130 is due to the IC mechanism within the PWN powered by PSR J18261256.

###### Acknowledgements.
The authors wish to thank the anonymous referee since his/her comments greatly improved our manuscript. This work is based on observations done with the XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and the US (NASA). This research was partially funded by Argentina Grants awarded by ANPCYT (PICT 0571/11) and University of Buenos Aires (UBACYT 20020150100098BA). G. C, E. G. and A. P. are Members of the Carrera del Investigador Científico of CONICET, Argentina. L. D. and L. S. are fellows of CONICET, Argentina.

## References

• Abdo et al. (2010) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJS, 187, 460
• Acero et al. (2015) Acero, F., Ackermann, M., Ajello, M., et al. 2015, ApJS, 218, 23
• Aharonian et al. (2006) Aharonian, F., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. 2006, A&A, 460, 365
• Anderson et al. (2011) Anderson, L. D., Bania, T. M., Balser, D. S., & Rood, R. T. 2011, The Astrophysical Journal Supplement Series, 194, 32
• Angüner et al. (2017) Angüner, E. O., Casanova, S., Oya, I., et al. 2017, ArXiv e-prints [\eprint[arXiv]1708.04844]
• Araudo et al. (2008) Araudo, A. T., Romero, G. E., Bosch-Ramon, V., & Paredes, J. M. 2008, International Journal of Modern Physics D, 17, 1889
• Blondin et al. (2001) Blondin, J. M., Chevalier, R. A., & Frierson, D. M. 2001, ApJ, 563, 806
• Bocchino et al. (2001) Bocchino, F., Warwick, R. S., Marty, P., et al. 2001, A&A, 369, 1078
• Bosch-Ramon et al. (2010) Bosch-Ramon, V., Romero, G. E., Araudo, A. T., & Paredes, J. M. 2010, A&A, 511, A8
• Castelletti et al. (2016) Castelletti, G., Giacani, E., & Petriella, A. 2016, A&A, 587, A71
• Clifton et al. (1992) Clifton, T. R., Lyne, A. G., Jones, A. W., McKenna, J., & Ashworth, M. 1992, MNRAS, 254, 177
• Condon et al. (1998) Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693
• Cordes & Lazio (2002) Cordes, J. M. & Lazio, T. J. W. 2002, ArXiv Astrophysics e-prints [\eprintarXiv:astro-ph/0207156]
• Evans et al. (2010) Evans, I. N., Primini, F. A., Glotfelty, K. J., et al. 2010, ApJS, 189, 37
• Frail & Scharringhausen (1997) Frail, D. A. & Scharringhausen, B. R. 1997, ApJ, 480, 364
• Gaensler et al. (2003) Gaensler, B. M., Schulz, N. S., Kaspi, V. M., Pivovaroff, M. J., & Becker, W. E. 2003, ApJ, 588, 441
• Gaensler & Slane (2006) Gaensler, B. M. & Slane, P. O. 2006, ARA&A, 44, 17
• Gaensler et al. (2000) Gaensler, B. M., Stappers, B. W., Frail, D. A., et al. 2000, MNRAS, 318, 58
• Giacani et al. (2014) Giacani, E., Rovero, A. C., Cillis, A., Pichel, A., & Dubner, G. 2014, ArXiv e-prints [\eprint[arXiv]1412.1673]
• Grondin et al. (2011) Grondin, M.-H., Funk, S., Lemoine-Goumard, M., et al. 2011, ApJ, 738, 42
• H. E. S. S. Collaboration et al. (2018) H. E. S. S. Collaboration, Abdalla, H., Abramowski, A., et al. 2018, A&A, 612, A1
• H. E. S. S. Collaboration et al. (2019) H. E. S. S. Collaboration, Abdalla, H., Aharonian, F., et al. 2019, A&A, 621, 116
• Holler et al. (2012) Holler, M., Schöck, F. M., Eger, P., et al. 2012, A&A, 539, A24
• Jankowski et al. (2018) Jankowski, F., van Straten, W., Keane, E. F., et al. 2018, MNRAS, 473, 4436
• Kargaltsev et al. (2013) Kargaltsev, O., Rangelov, B., & Pavlov, G. G. 2013, ArXiv e-prints [\eprint[arXiv]1305.2552]
• Kijak et al. (2007) Kijak, J., Gupta, Y., & Krzeszowski, K. 2007, A&A, 462, 699
• Manchester et al. (2005) Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993
• Mereghetti (2011) Mereghetti, S. 2011, Astrophysics and Space Science Proceedings, 21, 345
• Paron et al. (2013) Paron, S., Weidmann, W., Ortega, M. E., Albacete Colombo, J. F., & Pichel, A. 2013, MNRAS, 433, 1619
• Pavlov et al. (2008a) Pavlov, G. G., Kargaltsev, O., & Brisken, W. F. 2008a, ApJ, 675, 683
• Pavlov et al. (2008b) Pavlov, G. G., Kargaltsev, O., & Brisken, W. F. 2008b, ApJ, 675, 683
• Porquet et al. (2003) Porquet, D., Decourchelle, A., & Warwick, R. S. 2003, A&A, 401, 197
• Roberts et al. (2007) Roberts, M., Brogan, C., & Lyutikov, M. 2007, in Bulletin of the American Astronomical Society, Vol. 39, American Astronomical Society Meeting Abstracts, 997
• Roman-Duval et al. (2009) Roman-Duval, J., Jackson, J. M., Heyer, M., et al. 2009, ApJ, 699, 1153
• Rosen et al. (2016) Rosen, S. R., Webb, N. A., Watson, M. G., et al. 2016, A&A, 590, A1
• Safi-Harb et al. (2001) Safi-Harb, S., Harrus, I. M., Petre, R., et al. 2001, ApJ, 561, 308
• Voisin et al. (2016) Voisin, F., Rowell, G., Burton, M. G., et al. 2016, MNRAS, 458, 2813
• Wang (2011) Wang, W. 2011, Research in Astronomy and Astrophysics, 11, 824
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