Issue |
A&A
Volume 567, July 2014
|
|
---|---|---|
Article Number | L8 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201424261 | |
Published online | 25 July 2014 |
Discovery of TeV γ-ray emission from the pulsar wind nebula 3C 58 by MAGIC
1 IFAE, Campus UAB, 08193 Bellaterra, Spain
2 Università di Udine and INFN Trieste, 33100 Udine, Italy
3 INAF National Institute for Astrophysics, 00136 Rome, Italy
4 Università di Siena and INFN Pisa, 53100 Siena, Italy
5 Croatian MAGIC Consortium, Rudjer Boskovic Institute, University of Rijeka and University of Split, 10000 Zagreb, Croatia
6 Max-Planck-Institut für Physik, 80805 München, Germany
7 Universidad Complutense, 28040 Madrid, Spain
8 Inst. de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain
9 University of Łódź, 90236 Lodz, Poland
10 Deutsches Elektronen-Synchrotron (DESY), 15738 Zeuthen, Germany
11 ETH Zurich, 8093 Zurich, Switzerland
12 Universität Würzburg, 97074 Würzburg, Germany
13 Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas, 28040 Madrid, Spain
14 Institute of Space Sciences, 08193 Barcelona, Spain
15 Università di Padova and INFN, 35131 Padova, Italy
16 Technische Universität Dortmund, 44221 Dortmund, Germany
17 Unitat de Física de les Radiacions, Departament de Física, and CERES-IEEC, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
18 Universitat de Barcelona, ICC, IEEC-UB, 08028 Barcelona, Spain
19 Japanese MAGIC Consortium, Division of Physics and Astronomy, Kyoto University, Japan
20 Finnish MAGIC Consortium, Tuorla Observatory, University of Turku and Department of Physics, University of Oulu, Finland
21 Inst. for Nucl. Research and Nucl. Energy, 1784 Sofia, Bulgaria
22 Università di Pisa and INFN Pisa, 56126 Pisa, Italy
23 ICREA and Institute of Space Sciences, 08193 Barcelona, Spain
24 Università dell’Insubria and INFN Milano Bicocca, Como, 22100 Como, Italy
25 Now at: NASA Goddard Space Flight Center, Greenbelt MD 20771, USA; and Department of Physics and Department of Astronomy, University of Maryland, College Park MD 20742, USA
26 Now at École polytechnique fédérale de Lausanne (EPFL), Lausanne, Switzerland
27 Now at Institut für Astro- und Teilchenphysik, Leopold-Franzens- Universität Innsbruck, 6020 Innsbruck, Austria
28 Now at Finnish Centre for Astronomy with ESO (FINCA), Turku, Finland
29 Now at Astrophysics Science Division, Bhabha Atomic Research Centre, 400085 Mumbai, India
30 Also at INAF-Trieste, Italy
31 Now at School of Chemistry & Physics, University of Adelaide, 5005 Adelaide, Australia
32 Inst. de Astrofísica de Andalucía (CSIC), 18080 Granada, Spain
33 Also at Depto. de Física Teórica, Facultad de Ciencias de la Universidad de Zaragoza, Spain
34 Now at Centro de Física del Cosmos de Aragón, Teruel, Spain
Received: 23 May 2014
Accepted: 19 June 2014
Context. The pulsar wind nebula (PWN) 3C 58 is one of the historical very high-energy (VHE; E> 100 GeV) γ-ray source candidates. It is energized by one of the highest spin-down power pulsars known (5% of Crab pulsar) and it has been compared with the Crab nebula because of their morphological similarities. This object was previously observed by imaging atmospheric Cherenkov telescopes (Whipple, VERITAS and MAGIC), although it was not detected, with an upper limit of 2.3% Crab unit (C.U.) at VHE. It was detected by the Fermi Large Area Telescope (LAT) with a spectrum extending beyond 100 GeV.
Aims. We aim to extend the spectrum of 3C 58 beyond the energies reported by the Fermi Collaboration and probe acceleration of particles in the PWN up to energies of a few tens of TeV.
Methods. We analyzed 81 h of 3C 58 data taken in the period between August 2013 and January 2014 with the MAGIC telescopes.
Results. We detected VHE γ-ray emission from 3C 58 with a significance of 5.7σ and an integral flux of 0.65% C.U. above 1 TeV. According to our results, 3C 58 is the least luminous VHE γ-ray PWN ever detected at VHE and has the lowest flux at VHE to date. The differential energy spectrum between 400 GeV and 10 TeV is well described by a power-law function dφ/dE = f0(E/1 TeV)−Γ with f0 = (2.0 ± 0.4stat ± 0.6sys) × 10-13 cm-2 s-1 TeV-1 and Γ = 2.4 ± 0.2stat ± 0.2sys. The skymap is compatible with an unresolved source.
Conclusions. We report the first significant detection of PWN 3C 58 at TeV energies. We compare our results with the expectations of time-dependent models in which electrons upscatter photon fields. The best representation favors a distance to the PWN of 2 kpc and far-infrared (FIR) values similar to cosmic microwave background photon fields. If we consider an unexpectedly high FIR density, the data can also be reproduced by models assuming a 3.2 kpc distance. A low magnetic field, far from equipartition, is required to explain the VHE data. Hadronic contribution from the hosting supernova remnant (SNR) requires an unrealistic energy budget given the density of the medium, disfavoring cosmic-ray acceleration in the SNR as origin of the VHE γ-ray emission.
Key words: gamma rays: general / pulsars: general / ISM: supernova remnants
© ESO, 2014
1. Introduction
The supernova remnant 3C 58 (SNR G130.7+3.1) has a flat radio spectrum and is brightest near the center, therefore it was classified as a pulsar wind nebula (PWN; Weiler & Panagia 1978). It is centered on PSR J0205+6449, a pulsar discovered in 2002 with the Chandra X-ray observatory (Murray et al. 2002). It is widely assumed that 3C 58 is located at a distance of 3.2 kpc (Roberts et al. 1993), but a new H I measurement suggests a distance of 2 kpc (Kothes 2013). The age of the system is estimated to be ~2.5 kyr (Chevalier 2005) from the PWN evolution and energetics. However, the historical association of the PWN with the supernova of 1181 (Stephenson & Green 2002) and different measurements such as the velocity of the optical knots (Fesen et al. 2008), neutron star cooling models (proposing that an exotic cooling mechanism must operate in this pulsar; Slane et al. 2002), the proper motion of the pulsar (Gotthelf et al. 2007), and radio expansion of the nebula (Bietenholz 2006) derive ages ranging from 0.8 kyr up to 7 kyr. The pulsar has one of the highest spin-down powers known (Ė = 2.7 × 1037 erg s-1). The PWN has a size of 9′ × 6′ in radio, infrared (IR), and X-rays (Bietenholz et al. 2001; Bocchino et al. 2001; Slane et al. 2004, 2008). Its luminosity is L0.5−10 keV = 2.4 × 1034 erg s-1 in the X-ray band, which is more than 3 orders of magnitude lower than that of the Crab nebula (Torii et al. 2000). 3C 58 has been compared with the Crab because the jet-torus structure is similar (Slane et al. 2004). Because of these morphological similarities with the Crab nebula and its high spin-down power (5% of Crab), 3C 58 has historically been considered one of the PWNe most likely to emit γ rays.
The pulsar J0205+6449 has a period P = 65.68 ms, a spin-down rate Ṗ = 1.93 × 10-13 s s-1, and a characteristic age of 5.38 kyr (Murray et al. 2002). It was discovered by the Fermi-LAT in pulsed γ rays. The measured energy flux is Fγ> 0.1 GeV = (5.4 ± 0.2) × 10-11 erg cm-2 s-1 with a luminosity of Lγ> 0.1 GeV = (2.4 ± 0.1) × 1034 erg s-1, assuming a distance for the pulsar of 1.95 kpc (Xu et al. 2006). The spectrum is well described by a power-law with an exponential cutoff at Ecutoff = 1.6 GeV (Abdo et al. 2013). No pulsed emission was detected at energies above 10 GeV (Ackermann et al. 2013). In the off-peak region, defined as the region between the two γ-ray pulsed peaks (off-peak phase interval φ = 0.64–0.99), the Fermi Collaboration reported the detection of emission from 3C 58 (Abdo et al. 2013). The reported energy flux is (1.75 ± 0.68) × 10-11 erg cm-2 s-1 and the differential energy spectrum between 100 MeV and 316 GeV is well described by a power-law with photon index Γ = 1.61 ± 0.21. No hint of spatial extension was reported at those energies. The association of the high-energy unpulsed steady emission with the PWN is favored, although an hadronic origin related to the associated SNR can not be ruled out. 3C 58 was tagged as a potential TeV γ-ray source by the Fermi Collaboration (Ackermann et al. 2013).
The PWN 3C 58 was previously observed in the very high-energy (VHE; E> 100 GeV) range by several imaging atmospheric Cherenkov telescopes. The Whipple telescope reported an integral flux upper limit of 1.31 × 10-11 cm-2 s-1 (~19% C.U. at an energy threshold of 500 GeV; Hall et al. 2001), and VERITAS established upper limits at the level of 2.3% C.U. above an energy of 300 GeV (Aliu 2008). MAGIC-I observed the source in 2005 and established integral upper limits above 110 GeV at the level of 7.7×10-12 cm-2 s-1 (~4% C.U.; Anderhub et al. 2010). The improved sensitivity of the MAGIC telescopes with respect to previous observations and the Fermi-LAT results motivated us to perform deep VHE observations of the source.
2. Observations
MAGIC is a stereoscopic system of two imaging atmospheric Cherenkov telescopes situated in the Canary island of La Palma, Spain (28.8°N, 17.9° W at 2225 m above sea level). During 2011 and 2012 it underwent a major upgrade of the digital trigger, readout systems, and one of the cameras (Mazin et al. 2013). The system achieves a sensitivity of (0.71 ± 0.02)% of the Crab nebula flux above 250 GeV in 50 h at low zenith angles (Sitarek et al. 2013).
MAGIC observed 3C 58 in the period between 4 August 2013 to 5 January 2014 for 99 h, and after quality cuts, 81 h of the data were used for the analysis. The source was observed at zenith angles between 36° and 52°. The data were taken in wobble-mode pointing at four different positions situated 0.4° away from the source to evaluate the background simultaneously with 3C 58 observations (Fomin et al. 1994).
3. Analysis and results
The data were analyzed using the MARS analysis framework (Zanin et al. 2013). For every event, we calculated the size of the image (number of photoelectrons (phe) in the image of the shower) and the angular distance θ between the reconstructed arrival direction of the γ-ray and the source position. We determined the arrival direction combining of the individual telescope information using the Disp method (Zanin et al. 2013). A random-forest algorithm was used to estimate a global variable, the hadronness, which was used for background rejection (Albert et al. 2008). The energy of a given event was estimated from the image sizes of the shower, its reconstructed impact parameter, and the zenith angle with Monte Carlo (MC) filled look-up tables. We optimized the cut parameters for detecting a 1% C.U. point-like source on an independent Crab nebula data sample by maximizing the Li & Ma significance (Li & Ma 1983, Eq. (17), hereafter LiMa). We selected events with a θ2 angle <0.01 deg2, hadronness <0.18, and size in both telescopes >300 phe. For each pointing we used five off-source regions to estimate the background. We used MC γ rays to calculate the energy threshold of the analysis, defined here as the peak of the true energy distribution. To calculate the flux and spectral energy distribution (SED), the effective area of the telescopes was estimated using the MC simulation of γ-rays. The SED was finally unfolded using the Schmelling method to account for the energy resolution and energy bias of the instrument (Albert et al. 2007).
The applied cuts and the zenith angle of the observations yield an energy threshold of 420 GeV. The significance of the signal, calculated with the LiMa formula, is 5.7σ, which establishes 3C 58 as a γ-ray source. The θ2 distribution is shown in Fig. 1. As the five OFF positions were taken for each of the wobble positions, the OFF histograms were re-weighted depending on the time taken on each wobble position.
![]() |
Fig. 1 Distribution of squared angular distance, θ2, between the reconstructed arrival directions of γ-ray candidate events and the position of PSR 0205+6449 (red points). The distribution of θ2 for the OFF positions is also shown (gray filled histogram). The vertical dashed line defines the signal region ( |
We show in Fig. 2 the relative flux (excess/background) skymap, produced using the same cuts as for the θ2 calculation. The test statistics (TS) significance, which is the LiMa significance applied on a smoothed and modeled background estimate, is higher than 6 at the position of the pulsar PSR J0205+6449. The excess of the VHE skymap was fit with a Gaussian function. The best-fit position is RA(J2000) = 2h05m31(09)stat(11)sys s; Dec (J2000) = 64°51′(1)stat(1)sys. This position is statistically deviant by 2σ from the position of the pulsar, but is compatible with it if systematic errors are taken into account. In the bottom left of the image we show the point spread function (PSF) of the smeared map at the corresponding energies, which is the result of the sum in quadrature of the instrumental angular resolution and the applied smearing (4.7′ radius, at the analysis energy threshold). The extension of the VHE source is compatible with the instrument PSF. The VLA contours are coincident with the detected γ-ray excess.
![]() |
Fig. 2 Relative flux (excess/background) map for MAGIC observations. The cyan circle indicates the position of PSR J0205+6449 and the black cross shows the fitted centroid of the MAGIC image with its statistical uncertainty. In green we plot the contour levels for the TS starting at 4 and increasing in steps of 1. The magenta contours represent the VLA flux at 1.4 GHz (Condon et al. 1998), starting at 0.25 Jy and increasing in steps of 0.25 Jy. |
Figure 3 shows the energy spectrum for the MAGIC data, together with published predictions for the γ-ray emission from several authors, and two spectra obtained with three years of Fermi-LAT data, which were retrieved from the Fermi-LAT second pulsar-catalog (2PC, Abdo et al. 2013) and the Fermi high-energy LAT catalog (1FHL, Ackermann et al. 2013). The 1FHL catalog used events from the Pass 7 Clean class, which provides a substantial reduction of residual cosmic-ray background above 10 GeV, at the expense of a slightly smaller collection area, compared with the Pass 7 Source class that was adopted for 2PC (Ackermann et al. 2012). The two γ-ray spectra from 3C 58 reported in the 2PC and 1FHL catalogs agree within statistical uncertainties. The differential energy spectrum of the source is well fit by a single power-law function dφ/dE = f0(E/1 TeV)−Γ with f0 = (2.0 ± 0.4stat ± 0.6sys) × 10-13 cm-2 s-1 TeV-1, Γ = 2.4 ± 0.2stat ± 0.2sys and χ2 = 0.04/2. The systematic errors were estimated from the MAGIC performance paper (Aleksić et al. 2012) including the upgraded telescope performances. The integral flux above 1 TeV is FE> 1 TeV = 1.4 × 10-13 cm-2 s-1. Taking into account a distance of 2 kpc, the luminosity of the source above 1 TeV is erg s-1, where d2 is the distance normalized to 2 kpc.
![]() |
Fig. 3 3C 58 spectral energy distribution in the range between 0.1 GeV and 20 TeV. Red circles are the VHE points reported in this work. The best-fit function is drawn in red and the systematic uncertainty is represented by the yellow shaded area. Black squares and black arrows are taken from the Fermi-LAT second pulsar-catalog results (Abdo et al. 2013). Blue squares are taken from the Fermi high-energy LAT catalog (Ackermann et al. 2013). The magenta line is the SED prediction for 3C 58 taken from Fig. 10 of Bednarek & Bartosik (2003). The clear green dashed-dotted line is the SED predicted by Tanaka & Takahara (2013), assuming an age of 1 kyr, and the dark green dotted line is the prediction from the same paper, assuming an age of 2.5 kyr. The blue dashed line represents the SED predicted by Torres et al. (2013) assuming that the Galactic far-infrared (FIR) background is high enough to reach a flux detectable by the MAGIC sensitivity in 50 h. |
4. Discussion
Several models have been proposed that predict the VHE γ-ray emission of PWN 3C 58. Bucciantini et al. (2011) presented a one zone model of the spectral evolution of PWNe and applied it to 3C 58 using a distance of 3.2 kpc. The VHE emission from this model consists of inverse Compton (IC) scattering of cosmic microwave background (CMB) photons and optical-to-IR photons, and also of pion decay. The flux of γ rays above 400 GeV predicted by this model is about an order of magnitude lower than the observation.
Bednarek & Bartosik (2003) proposed a time-dependent model in which positrons gain energy in the process of resonant scattering by heavy nuclei. The VHE emission is produced by IC scattering of leptons off CMB, IR, and synchrotron photons and by the decay of pions due to the interaction of nuclei with the matter of the nebula. The age of 3C 58 is assumed to be 5 kyr, using a distance of 3.2 kpc and an expansion velocity of 1000 km s-1. According to this model, the predicted integral flux above 400 GeV is ~10-13 cm-2 s-1, while the integral flux above 420 GeV measured here is 5 × 10-13 cm-2 s-1. Calculations by Bednarek & Bartosik (2005), using the same model with an initial expansion velocity of 2000 km s-1 and considering IC scattering only from the CMB, are consistent with the observed spectrum. However, the magnetic field derived in this case is B ~ 14μG and it underestimates the radio emission of the nebula, although a more complex spectral shape might account for the radio nebula emission.
Tanaka & Takahara (2010) developed a time-dependent model of the spectral evolution of PWNe including synchrotron emission, synchrotron self-Compton, and IC. They evolved the electron energy distribution using an advective differential equation. To calculate the observability of 3C 58 at TeV energies they assumed a distance of 2 kpc and two different ages: 2.5 kyr and 1 kyr (Tanaka & Takahara 2013). For the 2.5 kyr age, they obtained a magnetic field B ≥ 17 μG, while for an age of 1 kyr, the magnetic field obtained is B = 40μG. The emission predicted by this model is closer to the Fermi result for an age of 2.5 kyr.
Martín et al. (2012) presented a different time-dependent leptonic diffusion-loss equation model without approximations, including synchrotron emission, synchrotron self-Compton, IC, and bremsstrahlung. They assumed a distance of 3.2 kpc and an age of 2.5 kyr to calculate the observability of 3C 58 at high energies (Torres et al. 2013). The predicted emission, without considering any additional photon source other than the CMB, is more than an order of magnitude lower than the flux reported here. It predicts VHE emission detectable by MAGIC in 50 h for an FIR-dominated photon background with an energy density of 5 eV/cm3. This would be more than one order of magnitude higher than the local IR density in the Galactic background radiation model used in GALPROP (~0.2 eV cm-3; Porter et al. 2006). The magnetic field derived from this model is 35 μG. To reproduce the observations, a large FIR background or a revised distance to the PWN of 2 kpc are required. In the first case, a nearby star or the SNR itself might provide the necessary FIR targets, although no detection of an enhancement has been found in the direction of the PWN. As we mentioned in Sect. 1, a distance of 2 kpc has recently been proposed by Kothes (2013) based on the recent H I measurements of the Canadian Galactic Plane Survey. At this distance, a lower photon density is required to fit the VHE data (Torres, priv. comm.).
We have shown different time-dependent models in this section that predict the VHE emission of 3C 58. The SEDs predicted by them are shown in Fig. 3. They use different assumptions for the evolution of the PWN and its emission. Bucciantini et al. (2011) divided the evolution of the SNR into phases and modeled the PWN evolution inside it. In Bednarek & Bartosik (2003) model, nuclei play an important role in accelerating particles inside the PWN. Torres et al. (2013) and Tanaka & Takahara (2010) modeled the evolution of the particle distribution by solving the diffusion-loss equation. Torres et al. (2013) fully solved the diffusion-loss equation, while Tanaka & Takahara (2010) neglected an escape term in the equation as an approximation. Another difference between these latter two models is that Tanaka & Takahara (2010) took synchrotron emission, synchrotron self-Compton and IC into account, while Torres et al. (2013) also consider the bremsstrahlung. The models that fit the γ-ray data derived a low magnetic field, far from equipartition, very low for a young PWN, but comparable with the value derived by Slane et al. (2008) using other data.
We also evaluated the possibility that the VHE γ rays might be produced by hadronic emission in the SNR that is associated with 3C 58. For the calculations, we assumed a distance of 2 kpc, a density of the medium of 0.38 cm-3 (Slane et al. 2004), and an initial energy of the supernova explosion of 1051 erg. We calculated the expected flux above 1 TeV following the approach reported in Drury et al. (1994). We find that the cosmic-ray nuclei have to be accelerated with an efficiency above 100% to account for our observed flux, thus we conclude that cosmic rays from the SNR are unlikely to be the origin of the TeV emission. Note, however, that dense clumps embedded within the remnant may change these conclusions (Fukui et al. 2012).
5. Conclusions
We have for the first time detected VHE γ rays up to TeV energies from the PWN 3C 58. The measured luminosity and flux make 3C 58 the least-luminous VHE γ-ray PWN known and the object with the lowest flux at VHE to date. Only a closer distance of 2 kpc or a high local FIR photon density can qualitatively reproduce the multiwavelength data of this object in the published models. Since the high FIR density is unexpected, the closer distance with FIR photon density comparable with the averaged value in the Galaxy is favored. The models that fit the γ-ray data derived magnetic fields which are very far from equipartition. Following the assumptions in Drury et al. (1994), it is highly unlikely that the measured flux comes from hadronic emission of the SNR. The VHE spectral results we presented will help to explain the broadband emission of this source within the available theoretical scenarios.
Acknowledgments
We would like to thank the Instituto de Astrofísica de Canarias for the excellent working conditions at the Observatorio del Roque de los Muchachos in La Palma. The support of the German BMBF and MPG, the Italian INFN, the Swiss National Fund SNF, and the Spanish MINECO is gratefully acknowledged. This work was also supported by the CPAN CSD2007-00042 and MultiDark CSD2009-00064 projects of the Spanish Consolider-Ingenio 2010 programme, by grant 127740 of the Academy of Finland, by the DFG Cluster of Excellence “Origin and Structure of the Universe”, by the Croatian Science Foundation (HrZZ) Project 09/176, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3, and by the Polish MNiSzW grant 745/N-HESS-MAGIC/2010/0. We would like to thank S. J. Tanaka for providing us useful information about his model.
References
- Abdo, A. A., Ajello, M., Allafort, A., et al. 2013, ApJS, 208, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Ackermann, M., Ajello, M., Albert, A., et al. 2012, ApJS, 203, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Ackermann, M., Ajello, M., Allafort, A., et al. 2013, ApJS, 209, 34 [NASA ADS] [CrossRef] [Google Scholar]
- Albert, J., Aliu, E., Anderhub, H., et al. 2007, Nucl. Instr. Meth. Phys. Res. A, 583, 494 [Google Scholar]
- Albert, J., Aliu, E., Anderhub, H., et al. 2008, Nucl. Instr. Meth. Phys. Res. A, 588, 424 [Google Scholar]
- Aleksić, J., Alvarez, E. A., Antonelli, L. A., et al. 2012, Astropart. Phys., 35, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Aliu, E. 2008, in AIP Conf. Ser. 1085, eds. F. A. Aharonian, W. Hofmann, & F. Rieger, 324 [Google Scholar]
- Anderhub, H., Antonelli, L. A., Antoranz, P., et al. 2010, ApJ, 710, 828 [NASA ADS] [CrossRef] [Google Scholar]
- Bednarek, W., & Bartosik, M. 2003, A&A, 405, 689 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bednarek, W., & Bartosik, M. 2005, J. Phys. G Nucl. Phys., 31, 1465 [NASA ADS] [CrossRef] [Google Scholar]
- Bietenholz, M. F. 2006, ApJ, 645, 1180 [NASA ADS] [CrossRef] [Google Scholar]
- Bietenholz, M. F., Kassim, N. E., & Weiler, K. W. 2001, ApJ, 560, 772 [NASA ADS] [CrossRef] [Google Scholar]
- Bocchino, F., Warwick, R. S., Marty, P., et al. 2001, A&A, 369, 1078 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bucciantini, N., Arons, J., & Amato, E. 2011, MNRAS, 410, 381 [Google Scholar]
- Chevalier, R. A. 2005, ApJ, 619, 839 [NASA ADS] [CrossRef] [Google Scholar]
- Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 [NASA ADS] [CrossRef] [Google Scholar]
- Drury, L. O., Aharonian, F. A., & Voelk, H. J. 1994, A&A, 287, 959 [NASA ADS] [Google Scholar]
- Fesen, R., Rudie, G., Hurford, A., & Soto, A. 2008, ApJS, 174, 379 [NASA ADS] [CrossRef] [Google Scholar]
- Fomin, V. P., Stepanian, A. A., Lamb, R. C., et al. 1994, Astropart. Phys., 2, 137 [NASA ADS] [CrossRef] [Google Scholar]
- Fukui, Y., Sano, H., Sato, J., et al. 2012, ApJ, 746, 82 [NASA ADS] [CrossRef] [Google Scholar]
- Gotthelf, E. V., Helfand, D. J., & Newburgh, L. 2007, ApJ, 654, 267 [NASA ADS] [CrossRef] [Google Scholar]
- Hall, T. A., Wakely, S. P., & VERITAS Collaboration. 2001, in Int. Cosmic Ray Conf., 6, 2485 [Google Scholar]
- Kothes, R. 2013, A&A, 560, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Li, T.-P., & Ma, Y.-Q. 1983, ApJ, 272, 317 [NASA ADS] [CrossRef] [Google Scholar]
- Martín, J., Torres, D. F., & Rea, N. 2012, MNRAS, 427, 415 [NASA ADS] [CrossRef] [Google Scholar]
- Mazin, D., Tescaro, D., Garczarczyk, M., et al. 2013, Proc. ICRC, 1071 [Google Scholar]
- Murray, S. S., Slane, P. O., Seward, F. D., Ransom, S. M., & Gaensler, B. M. 2002, ApJ, 568, 226 [NASA ADS] [CrossRef] [Google Scholar]
- Porter, T. A., Moskalenko, I. V., & Strong, A. W. 2006, ApJ, 648, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Roberts, D. A., Goss, W. M., Kalberla, P. M. W., Herbstmeier, U., & Schwarz, U. J. 1993, A&A, 274, 427 [NASA ADS] [Google Scholar]
- Sitarek, J., Carmona, E., Colin, P., et al. 2013, Proc. of 33rd ICRC, Rio de Janeiro, Brazil [arXiv:1308.0141] [Google Scholar]
- Slane, P. O., Helfand, D. J., & Murray, S. S. 2002, ApJ, 571, L45 [NASA ADS] [CrossRef] [Google Scholar]
- Slane, P., Helfand, D. J., van der Swaluw, E., & Murray, S. S. 2004, ApJ, 616, 403 [NASA ADS] [CrossRef] [Google Scholar]
- Slane, P., Helfand, D. J., Reynolds, S. P., et al. 2008, ApJ, 676, L33 [NASA ADS] [CrossRef] [Google Scholar]
- Stephenson, F. R., & Green, D. A. 2002, Historical supernovae and their remnants, eds. F. R. Stephenson, & D. A. Green, Int. Ser. Astron. Astrophys., Vol. 5 (Oxford: Clarendon Press), 5 [Google Scholar]
- Tanaka, S. J., & Takahara, F. 2010, ApJ, 715, 1248 [NASA ADS] [CrossRef] [Google Scholar]
- Tanaka, S. J., & Takahara, F. 2013, MNRAS, 429, 2945 [NASA ADS] [CrossRef] [Google Scholar]
- Torii, K., Slane, P. O., Kinugasa, K., Hashimotodani, K., & Tsunemi, H. 2000, PASJ, 52, 875 [NASA ADS] [CrossRef] [Google Scholar]
- Torres, D. F., Cillis, A. N., & Martín, J. 2013, ApJ, 763, L4 [NASA ADS] [CrossRef] [Google Scholar]
- Weiler, K. W., & Panagia, N. 1978, A&A, 70, 419 [NASA ADS] [Google Scholar]
- Xu, Y., Reid, M. J., Zheng, X. W., & Menten, K. M. 2006, Science, 311, 54 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Zanin, R., Carmona, E., Sitarek, J., et al. 2013, Proc. ICRC 2013, 773 [Google Scholar]
All Figures
![]() |
Fig. 1 Distribution of squared angular distance, θ2, between the reconstructed arrival directions of γ-ray candidate events and the position of PSR 0205+6449 (red points). The distribution of θ2 for the OFF positions is also shown (gray filled histogram). The vertical dashed line defines the signal region ( |
In the text |
![]() |
Fig. 2 Relative flux (excess/background) map for MAGIC observations. The cyan circle indicates the position of PSR J0205+6449 and the black cross shows the fitted centroid of the MAGIC image with its statistical uncertainty. In green we plot the contour levels for the TS starting at 4 and increasing in steps of 1. The magenta contours represent the VLA flux at 1.4 GHz (Condon et al. 1998), starting at 0.25 Jy and increasing in steps of 0.25 Jy. |
In the text |
![]() |
Fig. 3 3C 58 spectral energy distribution in the range between 0.1 GeV and 20 TeV. Red circles are the VHE points reported in this work. The best-fit function is drawn in red and the systematic uncertainty is represented by the yellow shaded area. Black squares and black arrows are taken from the Fermi-LAT second pulsar-catalog results (Abdo et al. 2013). Blue squares are taken from the Fermi high-energy LAT catalog (Ackermann et al. 2013). The magenta line is the SED prediction for 3C 58 taken from Fig. 10 of Bednarek & Bartosik (2003). The clear green dashed-dotted line is the SED predicted by Tanaka & Takahara (2013), assuming an age of 1 kyr, and the dark green dotted line is the prediction from the same paper, assuming an age of 2.5 kyr. The blue dashed line represents the SED predicted by Torres et al. (2013) assuming that the Galactic far-infrared (FIR) background is high enough to reach a flux detectable by the MAGIC sensitivity in 50 h. |
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.