Free Access
Issue
A&A
Volume 654, October 2021
Article Number A4
Number of page(s) 5
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/202140999
Published online 01 October 2021

© ESO 2021

1. Introduction

Pulsar wind nebulae (PWNe) interacting with the surrounding interstellar medium (ISM) can give rise to distinct extended morphological features. Torus-like structures and bipolar jets are observed in pulsars displaying slow proper motions, as well as bow-shaped shocks and extended cometary tails in pulsars moving at high velocities (see e.g., Gaensler & Slane 2006). In the X-ray band, these structures have been resolved in detail in a number of cases (Kargaltsev et al. 2008). A small subset of supersonically moving PWN (sPWN from now on) have, in addition, been observed to display unusually long X-ray outflows. Some of the most prominent examples of such outflows are found in the Guitar Nebula (Hui et al. 2012) and in the Lighthouse Nebula (Pavan et al. 2014), but the number of large-scale filaments from sPWNe is gradually increasing, see for example the recent findings in PSR J1135−6055 (Bordas & Zhang 2020), J1809–1917 (Klingler et al. 2020), PSR J2030+4415 (de Vries & Romani 2020), PSR J2055+2539 (Marelli et al. 2019), or the filaments found in the Galactic Center region such as G0.13-0.11 (Zhang et al. 2020), likely associated to a sPWN (see also previous findings in Kargaltsev et al. 2017). The origin of these extended jet-like structures is unknown, but their length, almost rectilinear geometry, and misaligned orientation with respect to the pulsar proper motion are difficult to be accommodated for in the framework of canonical pulsar jet formation theories, where jets are produced by the magnetic collimation of the plasma downstream of the pulsar wind termination shock (see e.g., Lyubarsky 2002; Komissarov & Lyubarsky 2004; Porth et al. 2014 and references therein). This led to alternative interpretations, including a scenario in which particles accelerated at the pulsar wind/medium shock interface escape the PWN, and diffused into the ambient medium magnetic field, yielding synchrotron emission which is then observed in the X-ray band as quasi-rectilinear jet-like structures (Bandiera 2008). Recent numerical simulations seem to support such a scenario, explaining also the asymmetric morphologies observed at opposite sides of the pulsar (Barkov et al. 2019a; Bucciantini et al. 2020).

In this paper we model the radio emission produced by high-energy electrons escaping the sPWN IGR J11014–6103, the Lighthouse Nebula, into the ISM, where they radiate in the ISM magnetic field. After a brief description of the source, in Sect. 2 we outline the model used to describe the emitting electron population and their diffusion into the ISM. Section 3 provides particle and radiation distributions, as well as synthetic maps at radio wavelengths used to predict the expected signal for dedicated observations of the source by highly sensitive facilities. Section 4 discusses the potential of radio observations to constrain the origin of the jet-like structures observed from sPWNe.

Originally discovered by INTEGRAL as an unidentified hard X-ray source, IGR J11014–6103 (also known as the “Lighthouse Nebula”; Pavan et al. 2014) is located ∼11 arcmins southwest of its likely progenitor supernova remnant SNR MSH 11-61A, at a distance of ∼7 kpc from the Sun (see García et al. 2012). The Lighthouse Nebula is powered by the pulsar PSR J1101−6101, displaying a period of P = 62.8 ms and an estimated spin-down energy of Lsd = 1.36 × 1036 erg s−1 (Halpern et al. 2014). Observations in the X-ray band have revealed a complex system morphology, including a trailing PWN, shaped in a narrow cone elongated toward the parent SNR, and prominent X-ray jet-like structures oriented nearly perpendicular to the PWN axis (see Pavan et al. 2016). At radio energies, extended emission has been detected with ATCA at 2 GHz from the IGR J11014–6103 tail (Pavan et al. 2014). This emission is not spatially coincident with that observed at X-rays, with the peak of the radio emission being separated further away along the tail by about 22″ with respect to the X-ray maximum.

The jets observed in the Lighthouse Nebula display a rectilinear structure extending for more than 10 pc into the surrounding medium (see Fig. 1). Their X-ray luminosity is unusually high, at the level of that retrieved from the pulsar and the tail structures, altogether representing a few percent of the pulsar spin-down power, Lsd ∼ 1036 erg s−1. The jet X-ray spectrum is moderately hard, with photon indices ΓX, jets ≈ 1.6, whereas the tail region displays ΓX, tail ≈ 1.9. Both the jet spectral hardness and its high X-ray efficiency may be expected if these are magnetic jets carrying a hard electron distribution, and there is an amplified magnetic field along the jet structures (see e.g., Bykov et al. 2017). However, accounting for the high-speed velocity of the source, vpsr ≳ 1000 km s−1, the absence of any significant bending along their ≳10 pc is striking prompts for alternative mechanisms for the production of such jet-like structures. Other sPWNe have been interpreted in a diffusion scenario (Bandiera 2008), and rectilinear structures or “filaments” have also been recently reported from the Galactic Center region (MeerKAT Collaboration 2018), and discussed as “kinetic jets” of escaping particles diffusing into the ambient magnetic field (Barkov & Lyutikov 2019). It is currently uncertain whether the jet-like structures in the Lighthouse Nebula are ballistic jets carrying an amplified magnetic field or filaments of escaping high-energy electrons from the sPWN into the ISM.

thumbnail Fig. 1.

Chandra/ACIS exposure-corrected broadband (0.5–7.0 keV) flux image of IGR J11014–6103 (the “Lighthouse Nebula”). The image has been produced using all Chandra archival observations of the Lighthouse Nebula, for a total of ∼365 ks, merged using CIAO’s merge_obs tool, and smoothed using a 2D Gaussian kernel with radius on a pixel size of 0.492″ and has been smoothed with a 2D-Gaussian kernel with radius = 1.7″ and σ = 0.7″.

2. Modeling the radio emission in a diffusion scenario

The pulsar PSR J1101-6101 in the Lighthouse Nebula displays a spin-down power of 1.36 × 1036 erg s−1 (Halpern et al. 2014). Observations with the Chandra satellite report a 2–10 keV X-ray flux from the jets at the level of 6 × 10−13 erg cm−2 s−1, whereas spectral analysis suggests a synchrotron origin rather than a thermal emission scenario (Pavan et al. 2014). Particles accelerated at the Lighthouse Nebula that are responsible for this emission could also emit in the radio band. Under these premises, we consider the presence of a nonthermal population of relativistic electrons injected close to the sPWN bow shock, carrying a total kinetic luminosity Linj ≈ LX ×  tsync/tdyn. Here LX is the jet X-ray luminosity derived from Chandra observations (we assume a distance to the Lighthouse Nebula of dLN = 7 kpc), tsync is the synchrotron cooling time for electrons emitting at the characteristic frequency c = 2 keV in a given magnetic field B, yr, and tdyn is the dynamical timescale to travel along the jet-like feature, tdyn ∼ ljet/vdyn; with ljet ∼ 10 pc and vdyn ∼ c, tdyn ∼ 32 yr. A total power Linj = 3 × 1034 erg s−1 is therefore injected, which is distributed in energy following a power-law with index Γ = 2 and an exponential cutoff at both the low- and the high-energy ends of the particle spectrum at and , respectively. A single power-law model has been assumed for simplicity in our model, although we note that PWNe could inject a more complex particle distribution displaying a spectral break following the presence of separate electron populations, for example those injected within the pulsar’s light cylinder, and electrons accelerated at the shock between the pulsar wind and the external medium (see e.g., Atoyan & Aharonian 1996; Venter & de Jager 2007; Torres et al. 2014). The low-energy cutoff is an assumption in our model, and it is treated as a free parameter, with two cases explored in this report: and 105. Such a low-energy cutoff may be expected in case an energy selection mechanism is in place in IGR J11014–6103 for particles able to escape the PWN. We use γmax = 108 for the high-energy particle distribution cutoff, which should be attainable in the bow shock of a runaway pulsar such as IGR J11014–6103 (compatible with the jet X-rays; see e.g., Bykov et al. 2017). As for the spectral index, we note that the value Γ = 2 is on the other hand adopted for simplicity, but also because it predicts potentially detectable fluxes. If the electron distribution were much harder, then little emission would be produced in the radio band.

Particles injected close to the bow shock are assumed to escape the immediate sPWN region and diffuse into the surrounding medium, where they evolve with time due to radiation losses, both synchrotron and inverse Compton (IC) emission, the latter accounting for the surrounding cosmic microwave background (CMB) and infrared radiation fields (UCMB ≈ 0.25 eV cm−3, UIR ≈ 0.2 eV cm−3, Popescu et al. 2017). We evaluate the resulting cooled particle distribution at the present time, tLN = 20 kyr (Halpern et al. 2014). A sufficient number of particle injection events are considered, which are equally spaced along the source trajectory during this time, with Δt shorter than the typical synchrotron and IC cooling timescales for particle Lorentz factors γe ≲ γmax. After injection, we consider electron diffusion in the surrounding ISM. A diffusion coefficient is assumed (Gabici et al. 2007). For simplicity, we assume that diffusion is isotropic from the injection point up to distances ldiff ∼ (Ddifft )1/2 (the medium is considered uniform, and thus adiabatic losses are not relevant).

Freshly injected particles, with injection times tinj ≤ 1 kyr, are presumably responsible for the X-ray emission observed from the jet-like features in the Lighthouse Nebula. For this injection period we consider distinct conditions of the ISM, which we label as the “JET region”, characterized by a magnetic field BJET = 5–30 μG that can be significantly higher than BISM ∼ 5μG under which electrons evolve after tinj > 1 kyr. The emission from particles injected until t ∼ 1 kyr diffusing into this JET region, and the emission resulting from their later diffusion in the ISM region (for t ≳ 1 kyr), was computed separately. We emphasize however that the injected electron particle population is the same in the JET and ISM regions, only the value of the corresponding magnetic field for each region varies. The former are evolved until tinj = 1 kyr only, whereas the latter evolves from tinj = 1 kyr up to tinj = 20 kyr. The different conditions from the JET and ISM regions follow the assumption that the JET region corresponds to a temporary feature that relocates as the pulsar moves, and to which we assign a typical crossing timescale of ∼1 kyr based on the JET apparent width and the pulsar velocity. The origin of the JET feature could be related to the pulsar motion pulling the threaded ISM magnetic field and, as the pulsar propagates, new ISM magnetic lines become entangled at the front of the PWN whereas the “older” ones relax and turn back into normal ISM magnetic field conditions.

Synchrotron emission is computed by employing the formulae in Pacholczyk (1970) together with the analytical approximations in Aharonian et al. (2010). The cumulative radio emission is evaluated throughout the area surrounding the source. Synthetic synchrotron emission maps are computed, which are then convolved with the beam of a radio telescope at GHz frequencies in order to obtain concrete observational predictions.

In this study we aim at predicting the radio emission produced in the surroundings of the sPWN, and therefore we do not attempt to model in detail the observed X-ray emission from the jet, including its morphological properties. In this regard, we only consider isotropic diffusion of particles injected in the JET and ISM regions. The effects that may result from anisotropies generated for example by the structured geometry of the underlying magnetic field, are therefore disregarded. Our estimates are based on assumed values for the injected kinetic energy, the energy distribution of the emitting particles, and on the ISM and JET regions’ magnetic field values. We do not consider higher-energy radiation channels, although we note that gamma-ray halos around PWN have long been discussed (see e.g., Aharonian 1995). In this regard, extended TeV emission has been recently reported from the PWN in Geminga (Abeysekara et al. 2017), which is consistent with free particle propagation in the ISM (see discussion in Giacinti et al. 2020). Finally, we focus in this study we focus on the outcomes that may be expected in the case of IGR J11014–6103. A general model extendable to any other sPWN is not explored in this paper. We note also that the values adopted for the parameters describing the injected particle distribution can be seen as rather optimistic (e.g. Γ could be harder than the adopted value of 2 and the low-energy cutoff could be more restrictive than 103). These, however, are relatively unconstrained parameters for the case of the source under study. At the same time, the adopted values are such that they can provide estimates of the expected radio emission that can then be contrasted with dedicated observations.

3. Results

3.1. Injected particle distribution

The energy distribution of particles injected along IGR J11014–6103’s lifetime is shown in Fig. 2, which displays the two cases described in Sect. 2 for the value of γmin = 103 and 105. For the high-energy cutoff, a value of is assumed in the two cases. Particles diffusing through the JET and the ISM regions are shown separately in Fig. 2. For the ISM region, individual single-injection events every 1 kyr (dotted lines) are considered. The particles evolve under synchrotron and IC losses while accumulating, forming the final particle distribution (solid lines). The older the injections are, the lower their contribution is at the high-energy end of the final distribution. Freshly accelerated particles are considered only for the JET region (dashed lines in Fig. 2).

thumbnail Fig. 2.

Energy distribution of non-thermal particles injected into the JET and ISM regions for different assumed values of the low-energy cutoff . Particle injection is discretized as occurring every tinj = 1 kyr (dotted lines). The final cumulative distribution at present time accounts for radiation cooling due to IC on the CMB and synchrotron losses. Distributions are shown for particles injected recently (tinj ≤ 1 kyr) that diffuse through the JET region (blue lines), and particles injected at tinj > 1 kyr diffusing into the ISM region (in red lines). In both cases, a magnetic field of BISM = BJET = 5 μG is assumed. A low cutoff for the particle distribution at injection is considered, adopting two different values and (solid and dashed lines, respectively). A high energy cutoff is also assumed in all cases. The distribution normalization is computed accounting for the source observed X-ray luminosity (see text for details).

3.2. Estimated synchrotron fluxes

The JET and the ISM regions differ in the assumed value for the magnetic field. We used BISM = 5μG to compute the synchrotron emission by each injection event and their total final contribution in the ISM. For the JET region, the value of BJET ranges between that of the ISM and the value inferred to produce the observed X-ray emission in the Lighthouse Nebula (for electrons at a maximal energy γmax ≈ 108) of BISM = 30 μG (Pavan et al. 2014). The energy flux of the predicted synchrotron emission is displayed in Fig. 3 for the cases and 105. Individual contributions and total emission in the ISM are displayed in red with dashed and solid lines, respectively. The energy flux from the JET region is in blue, for five different values of BJET: 5, 10, 15, 20 and 30μG.

thumbnail Fig. 3.

Left: spectral energy distribution of the emission produced by electrons injected in the ISM (in red) and JET regions (in blue). A magnetic field BISM = 5 μG is assumed for the ISM, whereas a range in the value of BJET = 5, 10, 15, 20 and 30 μG is probed for the JET region. Particle injection is considered up to tinj = 103 yr for the JET region and from tinj = 103 yr until tinj = 20 kyr for the ISM region (red-dotted lines). Particles are assumed to follow a power-law distribution with a low and a high energy cutoff and , respectively. Green and orange vertical bands denote Chandra and MeerKAT frequency coverage. Right: same as in the left panel but for a particle spectrum with a low-energy cutoff .

The radio emission for a low-energy cutoff reaches values of a few ×1029 erg s−1, which should be at levels detectable by existing facilities sensitive in this band. Radio emission is strongly suppressed for , which instead provides a larger number of particles at higher energies. Current X-ray satellites can detect the JET region in X-rays even for modest magnetic fields ∼15 μG. The total radio emission from the ISM region exceeds that of the JET region by a factor 2–10, depending on the value of BJET and the adopted value of γmin. The JET region is however brighter than any single injection event in the ISM when taking a similar injection time lapse ∼1 kyr. We include in Fig. 3 the observation bands of the radio telescope MeerKAT and the X-ray satellite Chandra, noted by orange and green bands, respectively.

3.3. Synthetic radio maps

The cooled particle distribution injected in each event producing the synchrotron emission described above is distributed isotropically from its injection point. The synchrotron maps displayed in Fig. 4 reflect this isotropy in the estimated spatial extension of the emission. The injection locations along the pulsar trajectory produce a “trail” of radio emission that is brighter along this direction, whereas the sum of the single events display an overall elliptical geometry encompassing the sPWN position. For the JET region, a single injection produces a radially symmetric extended source in our simplified model. In Fig. 4 the maps display the angular distribution of the predicted synchrotron fluxes, which has been convolved with the beam of a radio telescope of 15″ × 15″, such as that of the recently operative MeerKAT telescope (Jonas & MeerKAT Team 2016). We report only the case in Fig. 4, for which radio emission can reach detectable levels, above ∼50 μJy beam−1. For radio emission would be undetectable. The spatial distribution in our model provides evidence that the JET region is the brightest spot around IGR J11014–6103, but it is concentrated in a relatively small angular region of a few arcmins. We did not attempt to introduce further complexity with a nonisotropic particle distribution, which may produce extended structures in the radio maps, for example rectilinear jet-like features such as the ones reported in X-rays.

thumbnail Fig. 4.

Synthetic synchrotron specific flux maps displaying the expected radio emission from IGR J11014–6103 in the 0.9–1.7 GHz band in which MeerKAT operates, convolved with an instrument beam of 15″ × 15″ (Jonas & MeerKAT Team 2016). The contribution from particles injected all along the IGR J11014–6103’s lifetime, which escaped the nebula and diffused into the ISM (BISM = 5μG), is displayed in the left panel, and that from the JET region (BISM = 30μG) is shown in the right panel. A white star indicates the position of the pulsar, which is assumed to move leftwards along the horizontal dashed-line.

4. Discussion

The origin of the extended jet-like features observed in the Lighthouse Nebula and other sPWNs is unclear. The high velocity of these systems and the apparent rectilinear geometry of their associated outflows, with no appreciable signature of bending as expected due to the ambient medium ram pressure, suggest alternative mechanisms responsible for such jet-like features. Bandiera (2008) proposed a scenario in which the jets could be produced by high-energy electrons escaping the parent PWN and diffusing into the ambient magnetic field. Recent 3D magnetohydrodynamic (MHD) numerical simulations seem to support such a possibility (Barkov et al. 2019a; Bucciantini et al. 2020). Here we provide further observational predictions to probe this hypothesis in the context of the Lighthouse Nebula, modeling the radio emission from the immediate surroundings of the source.

The results reported in Sect. 3 indicate that, for suitable parameters describing the injected particle energy distributions escaping the sPWN, substantial radio emission around IGR J11014–6103 could be produced by these particles when diffusing into the ISM. We note that the value assumed for these parameters (e.g. Γ = 2 and ) were chosen so as to account for cases in which detectable radio fluxes may be expected, whereas in others (e.g., ) such emission would be suppressed. For the former case, the JET region provides the larger contribution, since BJET is larger than BISM. The fact that BJET > BISM could follow the presence of streaming instabilities generated in the flow of relativistic electrons, generating a turbulent component of the magnetic field responsible for the field amplification, and for the confinement of the radiating electrons themselves (Bykov et al. 2017). The total cumulative radio emission from the more diffuse ISM region, on the other hand, can be a factor of a few to ∼20 times higher than that in the JET region. According to our model, observations should therefore reveal a bright concentration toward the jet region, and a broader, dimmer diffuse radio nebular emission surrounding the source. Compared to the radio emission from the PWN tail reported in Pavan et al. (2014), the emission from the ISM region in our model displays somewhat softer fluxes. The nebular emission in the ISM region reported here might be intrinsically softer than that within the PWN in IGR J11014–6103. Otherwise, the presence of several populations of accelerated electrons in the PWN particle injection spectrum (see e.g., Atoyan & Aharonian 1996), could also make the radio fluxes harder. For the JET region, although we do not aim to describe the morphological and spectral properties of the extended jet-like features observed in the Lighthouse Nebula in X-rays, our model predicts fluxes compatible with the ones observed at this band, at a level L2−10 keV ≳ 1033 erg s−1, when assuming previously derived values of BJET of a few ×10 μG (Pavan et al. 2014). At higher energies, in the gamma-ray band, it is worth noting that, given the value of UIR the produced IC gamma-ray luminosity by the electrons responsible for the JET X-rays should be a factor times that observed in X-rays, L2−10 keV (a detailed treatment is left for future work).

As stated above the radio emission predicted in our model strongly depends on the assumed limits in the energy distribution of the emitting particles. As displayed in Fig. 3, right panel, a strong low-energy cutoff of would render these structures too dim to be detected. On the contrary, the synchrotron radio flux maps displayed in Fig. 4 for indicate that extended structures should be resolvable by current radio telescopes capable of reaching sensitivities at the level of ≳50 μJy beam−1. The value of depends on the mechanisms through which accelerated particles are able to escape the PWN region. For this to occur, the diffusion distance within the bow-shock region should be larger than the shocked wind layer thickness. Given that this is a several times smaller than the bow-shock size, the escape energy is just a factor order of 10 smaller than that given by the Hillas limit, which for the Lighthouse Nebula has an associated Lorentz factor of ≳108. Nevertheless, particles with lower energies may still be able to escape the PWN along the expanding PWN tail where the magnetic field density decreases and/or the geometry of the shock interface with the ambient ISM becomes more favorable for electrons to diffuse away. In this regard, it is worth noting that for the case of the Lighthouse Nebula, the base of the main jet-like feature seems to be shifted to a position behind the pulsar with respect to the incoming medium (see e.g., Pavan et al. 2016; see also the inset in Fig. 9 in Kargaltsev et al. 2017), where the PWN conditions could imply lower values for the particles’ Lorentz factor to escape the nebula. Alternatively, particle escape can follow magnetic reconnection between the sPWN and the ISM magnetic fields (see e.g., Barkov et al. 2019b; Bucciantini et al. 2020), as suggested in the case of the rectilinear structures detected in J1809−1917, for which the gyroradius of the X-ray emitting electrons is by far shorter than the typical size of the nebula (Klingler et al. 2020).

Finally, although in this report we have focused on the Lighthouse Nebula, we note that other sPWNe may also display extended radio emission from electrons escaping their nebulae. If confirmed, the detection of these radio structures would represent a footprint for an efficient particle escape mechanism taking place in these systems, allowing in turn the formation of extended, rectilinear, and misaligned jet-like structures as observed in some of the most extreme sPWNe.

Acknowledgments

PB, VBR, JMP and XZ acknowledge the financial support by the Spanish Ministerio de Economía, Industria y Competitividad (MINEICO/FEDER, UE) under grant AYA2016-76012-C3-1-P, from the State Agency for Research of the Spanish Ministry of Science and Innovation under grant PID2019-105510GB-C31 and through the “Unit of Excellence María de Maeztu 2020-2023” award to the Institute of Cosmos Sciences (CEX2019-000918-M), and by the Catalan DEC grant 2017 SGR 643.

References

  1. Abeysekara, A. U., Albert, A., Alfaro, R., et al. 2017, Science, 358, 911 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aharonian, F. A. 1995, Nucl. Phys. B Proc. Suppl., 39, 193 [NASA ADS] [CrossRef] [Google Scholar]
  3. Aharonian, F. A., Kelner, S. R., & Prosekin, A. Y. 2010, Phys. Rev. D, 82, 043002 [Google Scholar]
  4. Atoyan, A. M., & Aharonian, F. A. 1996, MNRAS, 278, 525 [NASA ADS] [Google Scholar]
  5. Bandiera, R. 2008, A&A, 490, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Barkov, M. V., & Lyutikov, M. 2019, MNRAS, 489, L28 [NASA ADS] [CrossRef] [Google Scholar]
  7. Barkov, M. V., Lyutikov, M., Klingler, N., & Bordas, P. 2019a, MNRAS, 485, 2041 [NASA ADS] [CrossRef] [Google Scholar]
  8. Barkov, M. V., Lyutikov, M., & Khangulyan, D. 2019b, MNRAS, 484, 4760 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bordas, P., & Zhang, X. 2020, A&A, 644, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bucciantini, N., Olmi, B., & Del Zanna, L. 2020, J. Phys. Conf. Ser., 1623, 012002 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bykov, A. M., Amato, E., Petrov, A. E., Krassilchtchikov, A. M., & Levenfish, K. P. 2017, Space Sci. Rev., 207, 235 [NASA ADS] [CrossRef] [Google Scholar]
  12. de Vries, M., & Romani, R. W. 2020, ApJ, 896, L7 [NASA ADS] [CrossRef] [Google Scholar]
  13. Gabici, S., Aharonian, F. A., & Blasi, P. 2007, Ap&SS, 309, 365 [CrossRef] [Google Scholar]
  14. Gaensler, B. M., & Slane, P. O. 2006, ARA&A, 44, 17 [NASA ADS] [CrossRef] [Google Scholar]
  15. García, F., Combi, J. A., Albacete-Colombo, J. F., et al. 2012, A&A, 546, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Giacinti, G., Mitchell, A. M. W., López-Coto, R., et al. 2020, A&A, 636, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Halpern, J. P., Tomsick, J. A., Gotthelf, E. V., et al. 2014, ApJ, 795, L27 [NASA ADS] [CrossRef] [Google Scholar]
  18. Hui, C. Y., Huang, R. H. H., Trepl, L., et al. 2012, ApJ, 747, 74 [NASA ADS] [CrossRef] [Google Scholar]
  19. Jonas, J., & MeerKAT Team 2016, MeerKAT Science: On the Pathway to the SKA, 1 [Google Scholar]
  20. Kargaltsev, O., & Pavlov, G. G. 2008, in 40 Years of Pulsars: Millisecond Pulsars, Magnetars and More, eds. C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi, AIP Conf. Ser., 983, 171 [NASA ADS] [CrossRef] [Google Scholar]
  21. Kargaltsev, O., Pavlov, G. G., Klingler, N., & Rangelov, B. 2017, J. Plasma Phys., 83, 635830501 [CrossRef] [Google Scholar]
  22. Klingler, N., Yang, H., Hare, J., et al. 2020, ApJ, 901, 157 [NASA ADS] [CrossRef] [Google Scholar]
  23. Komissarov, S. S., & Lyubarsky, Y. E. 2004, MNRAS, 349, 779 [NASA ADS] [CrossRef] [Google Scholar]
  24. Lyubarsky, Y. E. 2002, MNRAS, 329, L34 [NASA ADS] [CrossRef] [Google Scholar]
  25. Marelli, M., Tiengo, A., De Luca, A., et al. 2019, A&A, 624, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. MeerKAT Collaboration 2018, Mon. Notes Astron. Soc. S. Afr., 77, 102 [Google Scholar]
  27. Pacholczyk, A. G. 1970, Radio Astrophysics. Nonthermal Processes in Galactic and Extragalactic Sources [Google Scholar]
  28. Pavan, L., Bordas, P., Pühlhofer, G., et al. 2014, A&A, 562, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Pavan, L., Pühlhofer, G., Bordas, P., et al. 2016, A&A, 591, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Popescu, C. C., Yang, R., Tuffs, R. J., et al. 2017, MNRAS, 470, 2539 [NASA ADS] [CrossRef] [Google Scholar]
  31. Porth, O., Komissarov, S. S., & Keppens, R. 2014, MNRAS, 438, 278 [NASA ADS] [CrossRef] [Google Scholar]
  32. Torres, D. F., Cillis, A., Martín, J., & de Oña Wilhelmi, E. 2014, J. High Energy Astrophys., 1, 31 [NASA ADS] [CrossRef] [Google Scholar]
  33. Venter, C., & de Jager, O. C. 2007, in WE-Heraeus Seminar on Neutron Stars and Pulsars 40 years after the Discovery, eds. W. Becker, & H. H. Huang, 40 [Google Scholar]
  34. Zhang, S., Zhu, Z., Li, H., et al. 2020, ApJ, 893, 3 [Google Scholar]

All Figures

thumbnail Fig. 1.

Chandra/ACIS exposure-corrected broadband (0.5–7.0 keV) flux image of IGR J11014–6103 (the “Lighthouse Nebula”). The image has been produced using all Chandra archival observations of the Lighthouse Nebula, for a total of ∼365 ks, merged using CIAO’s merge_obs tool, and smoothed using a 2D Gaussian kernel with radius on a pixel size of 0.492″ and has been smoothed with a 2D-Gaussian kernel with radius = 1.7″ and σ = 0.7″.

In the text
thumbnail Fig. 2.

Energy distribution of non-thermal particles injected into the JET and ISM regions for different assumed values of the low-energy cutoff . Particle injection is discretized as occurring every tinj = 1 kyr (dotted lines). The final cumulative distribution at present time accounts for radiation cooling due to IC on the CMB and synchrotron losses. Distributions are shown for particles injected recently (tinj ≤ 1 kyr) that diffuse through the JET region (blue lines), and particles injected at tinj > 1 kyr diffusing into the ISM region (in red lines). In both cases, a magnetic field of BISM = BJET = 5 μG is assumed. A low cutoff for the particle distribution at injection is considered, adopting two different values and (solid and dashed lines, respectively). A high energy cutoff is also assumed in all cases. The distribution normalization is computed accounting for the source observed X-ray luminosity (see text for details).

In the text
thumbnail Fig. 3.

Left: spectral energy distribution of the emission produced by electrons injected in the ISM (in red) and JET regions (in blue). A magnetic field BISM = 5 μG is assumed for the ISM, whereas a range in the value of BJET = 5, 10, 15, 20 and 30 μG is probed for the JET region. Particle injection is considered up to tinj = 103 yr for the JET region and from tinj = 103 yr until tinj = 20 kyr for the ISM region (red-dotted lines). Particles are assumed to follow a power-law distribution with a low and a high energy cutoff and , respectively. Green and orange vertical bands denote Chandra and MeerKAT frequency coverage. Right: same as in the left panel but for a particle spectrum with a low-energy cutoff .

In the text
thumbnail Fig. 4.

Synthetic synchrotron specific flux maps displaying the expected radio emission from IGR J11014–6103 in the 0.9–1.7 GHz band in which MeerKAT operates, convolved with an instrument beam of 15″ × 15″ (Jonas & MeerKAT Team 2016). The contribution from particles injected all along the IGR J11014–6103’s lifetime, which escaped the nebula and diffused into the ISM (BISM = 5μG), is displayed in the left panel, and that from the JET region (BISM = 30μG) is shown in the right panel. A white star indicates the position of the pulsar, which is assumed to move leftwards along the horizontal dashed-line.

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.