MAGIC Upper Limits G65.1+0.6

MAGIC Upper Limits for two Milagro-detected,
Bright Fermi Sources in the Region of SNR G65.1+0.6

J. Aleksićaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , L. A. Antonellibbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , P. Antoranzccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , M. Backesddaffiliation: Technische Universität Dortmund, D-44221 Dortmund, Germany , J. A. Barrioeeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , D. Bastieriffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , J. Becerra Gonzálezggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain affiliationmark: hhaffiliation: Depto. de Astrofisica, Universidad, E-38206 La Laguna, Tenerife, Spain , W. Bednarekiiaffiliation: University of Łódź, PL-90236 Lodz, Poland , A. Berdyuginjjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , K. Bergerggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain , E. Bernardinikkaffiliation: Deutsches Elektronen-Synchrotron (DESY), D-15738 Zeuthen, Germany , A. Bilandllaffiliation: ETH Zurich, CH-8093 Switzerland , O. Blanchaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , R. K. Bockmmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , A. Bollerllaffiliation: ETH Zurich, CH-8093 Switzerland , G. Bonnolibbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , P. Bordasnnaffiliation: Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain , D. Borla Tridonmmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , V. Bosch-Ramonnnaffiliation: Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain , D. Boseeeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , I. Braunllaffiliation: ETH Zurich, CH-8093 Switzerland , T. Bretzooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , M. Camaraeeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , E. Carmonammaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , A. Carosibbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , P. Colinmmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , J. L. Contreraseeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , J. Cortinaaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , S. Covinobbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , F. Dazzippaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy affiliationmark: **affiliation: supported by INFN Padova , A. De Angelisppaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , E. De Cea del Pozoqqaffiliation: Institut de Ciències de l’Espai (IEEC-CSIC), E-08193 Bellaterra, Spain , B. De Lottoppaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , M. De Mariappaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , F. De Sabatappaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , C. Delgado Mendezggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain affiliationmark: ****affiliation: now at: Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas , A. Diago Ortegaggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain affiliationmark: hhaffiliation: Depto. de Astrofisica, Universidad, E-38206 La Laguna, Tenerife, Spain , M. Doertddaffiliation: Technische Universität Dortmund, D-44221 Dortmund, Germany , A. Domínguezrraffiliation: Inst. de Astrofísica de Andalucía (CSIC), E-18080 Granada, Spain , D. Dominis Presterssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , D. Dornerllaffiliation: ETH Zurich, CH-8093 Switzerland , M. Doroffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , D. Elsaesserooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , M. Errandoaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , D. Ferencssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , M. V. Fonsecaeeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , L. Fontttaffiliation: Universitat Autònoma de Barcelona, E-08193 Bellaterra, Spain , R. J. García Lópezggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain affiliationmark: hhaffiliation: Depto. de Astrofisica, Universidad, E-38206 La Laguna, Tenerife, Spain , M. Garczarczykggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain , M. Gaugggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain , G. Giavittoaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , N. Godinovićssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , D. Hadaschqqaffiliation: Institut de Ciències de l’Espai (IEEC-CSIC), E-08193 Bellaterra, Spain , A. Herreroggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain affiliationmark: hhaffiliation: Depto. de Astrofisica, Universidad, E-38206 La Laguna, Tenerife, Spain , D. Hildebrandllaffiliation: ETH Zurich, CH-8093 Switzerland , D. Höhne-Mönchooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , J. Hosemmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , D. Hrupecssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , T. Joglermmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , S. Klepseraaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , T. Krähenbühlllaffiliation: ETH Zurich, CH-8093 Switzerland , D. Kranichllaffiliation: ETH Zurich, CH-8093 Switzerland , J. Krausemmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , A. La Barberabbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , E. Leonardoccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , E. Lindforsjjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , S. Lombardiffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , F. Longoppaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , M. Lópezffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , E. Lorenzllaffiliation: ETH Zurich, CH-8093 Switzerland affiliationmark: mmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , P. Majumdarkkaffiliation: Deutsches Elektronen-Synchrotron (DESY), D-15738 Zeuthen, Germany , G. Manevauuaffiliation: Inst. for Nucl. Research and Nucl. Energy, BG-1784 Sofia, Bulgaria , N. Mankuzhiyilppaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , K. Mannheimooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , L. Maraschibbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , M. Mariottiffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , M. Martínezaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , D. Mazinaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , M. Meucciccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , J. M. Mirandaccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , R. Mirzoyanmmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , H. Miyamotommaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , J. Moldónnnaffiliation: Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain , A. Moralejoaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , D. Nietoeeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , K. Nilssonjjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , R. Oritommaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , I. Oyaeeaffiliation: Universidad Complutense, E-28040 Madrid, Spain , R. Paoletticcaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , J. M. Paredesnnaffiliation: Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain , S. Partiniccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , M. Pasanenjjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , F. Paussllaffiliation: ETH Zurich, CH-8093 Switzerland , R. G. Pegnaccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , M. A. Perez-Torresrraffiliation: Inst. de Astrofísica de Andalucía (CSIC), E-18080 Granada, Spain , M. Persicppaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy affiliationmark: vvaffiliation: INAF/Osservatorio Astronomico and INFN, I-34143 Trieste, Italy , L. Peruzzoffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , J. Pochonggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain , F. Pradarraffiliation: Inst. de Astrofísica de Andalucía (CSIC), E-18080 Granada, Spain , P. G. Prada Moroniccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , E. Prandiniffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , N. Puchadesaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , I. Puljakssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , I. Reichardtaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , R. Reinthaljjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , W. Rhodeddaffiliation: Technische Universität Dortmund, D-44221 Dortmund, Germany , M. Ribónnaffiliation: Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain , J. Ricowwaffiliation: ICREA, E-08010 Barcelona, Spain affiliationmark: aaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , M. Rissillaffiliation: ETH Zurich, CH-8093 Switzerland , S. Rügamerooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , A. Saggionffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , K. Saitommaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , T. Y. Saitommaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , M. Salvatibbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , M. Sánchez-Condeggaffiliation: Inst. de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain affiliationmark: hhaffiliation: Depto. de Astrofisica, Universidad, E-38206 La Laguna, Tenerife, Spain , K. Sataleckakkaffiliation: Deutsches Elektronen-Synchrotron (DESY), D-15738 Zeuthen, Germany , V. Scalzottoffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , V. Scapinppaffiliation: Università di Udine, and INFN Trieste, I-33100 Udine, Italy , C. Schultzffaffiliation: Università di Padova and INFN, I-35131 Padova, Italy , T. Schweizermmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , M. Shaydukmmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , A. Sierpowska-Bartosikiiaffiliation: University of Łódź, PL-90236 Lodz, Poland , A. Sillanpääjjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , J. Sitarekmmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany affiliationmark: iiaffiliation: University of Łódź, PL-90236 Lodz, Poland , D. Sobczynskaiiaffiliation: University of Łódź, PL-90236 Lodz, Poland , F. Spanierooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , S. Spirobbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , A. Stamerraccaffiliation: Università di Siena, and INFN Pisa, I-53100 Siena, Italy , B. Steinkemmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , J. Storzooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , N. Strahddaffiliation: Technische Universität Dortmund, D-44221 Dortmund, Germany , J. C. Struebigooaffiliation: Universität Würzburg, D-97074 Würzburg, Germany , T. Suricssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , L. Takalojjaffiliation: Tuorla Observatory, University of Turku, FI-21500 Piikkiö, Finland , F. Tavecchiobbaffiliation: INAF National Institute for Astrophysics, I-00136 Rome, Italy , P. Temnikovuuaffiliation: Inst. for Nucl. Research and Nucl. Energy, BG-1784 Sofia, Bulgaria , T. Terzićssaffiliation: Croatian MAGIC Consortium, Institute R. Boskovic, University of Rijeka and University of Split, HR-10000 Zagreb, Croatia , D. Tescaroaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , M. Teshimammaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , D. F. Torreswwaffiliation: ICREA, E-08010 Barcelona, Spain affiliationmark: qqaffiliation: Institut de Ciències de l’Espai (IEEC-CSIC), E-08193 Bellaterra, Spain , H. Vankovuuaffiliation: Inst. for Nucl. Research and Nucl. Energy, BG-1784 Sofia, Bulgaria , R. M. Wagnermmaffiliation: Max-Planck-Institut für Physik, D-80805 München, Germany , Q. Weitzelllaffiliation: ETH Zurich, CH-8093 Switzerland , V. Zabalzannaffiliation: Universitat de Barcelona (ICC/IEEC), E-08028 Barcelona, Spain , F. Zandanelrraffiliation: Inst. de Astrofísica de Andalucía (CSIC), E-18080 Granada, Spain , R. Zaninaaaffiliation: IFAE, Edifici Cn., Campus UAB, E-08193 Bellaterra, Spain , klepser@ifae.es, decea@ieec.uab.es
Abstract

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around . We therefore report differential flux upper limits, assuming the emission to be point-like (°) or within a radius of 0.3°. In the point-like scenario, the flux limits around are at the level of 3% and 2% of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond , resulting in a photon index harder than 2.2 in the TeV band.

pulsars: individual(PSR J1957+2831, LAT PSR J1954+2836, LAT PSR J1958+2846) — ISM: supernova remnants

1 Introduction

In February 2009, the Fermi collaboration published a list of the most significant gamma-ray sources above , detected by the large area telescope (LAT) within 3 months of observation (Abdo et al., 2009b). The energy spectra of these sources often extend to several GeV, where at some point the steeply falling flux levels are too low to be detected by the limited detection area of a satellite instrument. Many of the LAT sources are hosted by our galaxy, and 34 of those are within the field of view of the Milagro gamma-ray observatory, which was located near Los Alamos, New Mexico (Abdo et al. (2009a) and references within). The Milagro collaboration, therefore, reinvestigated their previously gained skymap to look for counterparts to these GeV sources (Abdo et al., 2009a). The sensitivity of the instrument peaks in the energy range , although it ultimately depends on the energy spectrum and the declination of the source. Due to the reduced trial factor, 14 new Milagro sources could be claimed with confidence levels above .

Two of these Milagro-detected Fermi bright sources are in the vicinity of G65.1+0.6, a faint supernova remnant (SNR) first reported by Landecker et al. (1990). Tian & Leahy (2006) suggested its distance to be and a Sedov age of . An association to the radio pulsar PSR J1957+2831 was suggested by Lorimer et al. (1998).

The gamma-ray emission in the region of G65.1+0.6 was first detected by the COS-B satellite (Swanenburg et al., 1981) as 2CG065+00, and later confirmed by the EGRET satellite (3EG J1958+2909) in Hartman et al. (1999), where a possible extension or multiple sources were denoted. As of now, the two sources could be detected as two individual sources by the LAT as 1FGL J1954.3+2836 and 1FGL J1958.6+2845, as reported in the first year catalog of Fermi sources (Abdo et al., 2010b). They were analysed and reported as gamma-ray pulsars found through blind search in Saz Parkinson et al. (2010), Abdo et al. (2009c) and Abdo et al. (2010a). Their periods (, ), spin-down luminosities (, ) and energy cutoffs (, ), but also the characteristic magnetic fields and ages (, ) lie in the average range of all Fermi pulsars. In the following, we will use the names J1954 and J1958 for these 1-year catalog sources, and J1954/J1958 to specifically refer to the 3-month bright source list positions.

The detections by Milagro revealed significances of for J1954 and for J1958 (Abdo et al., 2009a). Flux values are stated for a characteristic median energy of . The angular resolution of Milagro is about , so these values can be expected to hold also for the 1-year catalog positions of the sources, which are offset by . Since gamma-ray pulsars typically have energy cutoffs as low as several GeV, the Milagro signals, if real, can be expected not to be caused directly by the pulsars, but possibly by associated objects, such as a pulsar wind nebula (PWN) or, in the case of J1954, the shell of the SNR, which surrounds the pulsar. Gamma rays at TeV energies may also be produced in an interaction of the shell with a possibly coincident molecular cloud. This cloud may be located at the position of the infra-red source IRAS 19520+2759, which has associated CO line, HO and OH maser emission at a similar distance as the SNR (Arquilla & Kwok, 1987).

The MAGIC telescopes use the Cherenkov imaging technique and are located on the Canary Island of La Palma (28.8°N, 17.8°W, 2220 m a.s.l.). It is the instrument with the lowest energy threshold among all Cherenkov telescopes. In single-telescope observations, as presented here, the nominal threshold in low zenith angle observations is (Albert et al., 2008a). It is therefore the instrument of choice to connect the upper ends of the Fermi spectra, which typically end at some tens of GeV, with the detections provided by Milagro at some tens of TeV. To bridge these spectra in the TeV range, and identify possible object associations for the Milagro signals were the main motivations for our investigation.

2 Observations

Motivated by the fact that J1954 was not marked as a pulsar in the initial Fermi Bright Source List, we observed J1954 as a main target in July and August 2009. The observations were carried out in false source tracking (wobble) mode (Fomin et al., 1994), which yielded two datasets with offsets of 0.4° in RA from this source, see Figure 1. The wobble position was altered every , and the data were taken at zenith angles between 0 and 43 degrees. At the time, the second MAGIC telescope was still under comissioning, so the analysis presented here used only the data from the stand-alone MAGIC-I telescope.

Quality selection cuts were applied to the event rate, the spread of hardware-sensitive shower image parameters, and few parameters that characterize the transparency of the atmosphere, such as the sky temperature and humidity. After all data selection cuts, 25.5 hours of high quality data were left for the analysis of J1954.

Besides this, we took advantage of the fact that J1958 is in the field of view of one of the two wobble positions. It could be analysed in a specific way described in the next section, yielding 13.8 hours of effective observation time.

3 Analysis and Results

The data were analysed in the MARS analysis framework (Moralejo et al., 2009), which is the standard software used in the analysis of MAGIC data. After the air shower images of the photomultiplier tube camera are calibrated, and times and charges of each pixel are extracted, a three-stage image cleaning is applied to filter out uncorrelated noise from the data acquisition electronics and the night sky background (Aliu et al., 2009). Shower image parameters are calculated, and the Random Forest method (Albert et al., 2008b) is used to derive estimators for the shower direction, its energy, and its likeliness to be of hadronic origin.

The observational setup described in the previous section requires a different analysis treatment of both sources. The main target of the observation, J1954, can be analysed using standard wobble analysis procedures. This means that the photon flux from the source is compared to the one on the opposite side of the camera (anti-source, solid squares in Figure 1). In this way, exposure inhomogeneities that can arise from imperfections in the photo multiplier tubes, trigger electronics or the signal transmission, cancel out, because both wobble positions were equally populated.

In the case of J1958, a wobble analysis, using only one of the two wobble positions, is not guaranteed to cancel out these inhomogeneities. Therefore, the analysis was done in ON/OFF manner, using the near wobble sample as ON-source data and the far wobble sample as OFF-source data (hollow squares in Figure 1). Having the OFF-source at the same position in relative camera coordinates as the source in the ON sample, the exposure inhomogeneities cancel out.

To test the presence of a gamma-ray signal, the distributions of squared angular distances () between photon directions and the source positions were used. In these -plots, a signal is expected to produce an excess where approaches zero. However, the integrals over the expected signal regions of the distributions agree, for both sources, well with the corresponding integrals done with respect to the OFF regions. Furthermore, the shapes of the ON and OFF distributions agree well with each other, and are sufficiently flat to exclude the unlikely case of an emission occuring by chance at both ON and OFF locations at a similar flux level, supressing the significance of the comparison. Such a coincidence was also excluded by cross-checking the analysis using several OFF regions, and by thoroughly investigating the skymaps with different background estimation algorithms. Consequently, we report the absence of a significant signal for both sources.

Besides this analysis, which takes advantage of a-priori defined source locations, also a 2.8°2.0°  skymap of the area was investigated at different energy ranges. The trial factors implied by these searches were typically between 120 and 260, depending on the PSF, which is smaller for high energies. Taking into account these trials, no significant signal was found at any energy.

We convert our data into three differential flux upper limits for each source. To do so, the data was divided into three bins of estimated energy, delimited by , , and . For each bin, an event number upper limit is calculated from the above-mentioned -plots, at confidence level (c.l.) after Rolke et al. (2005), and assuming an efficiency systematic error of (Albert et al., 2008a). For each limit, a power law energy spectrum with a photon index of is assumed both for the simulation of the effective area, and the conversion of event number upper limits to flux upper limits. The influence of this photon index is minor, since the energy ranges are sufficiently small.

Finally, since the data can only be selected by estimated energy, a Monte Carlo simulation (MC) is used to estimate the median true energy of the data that remains after all cuts.

This analysis assumes a source extension similar or smaller than the point spread function (PSF) of MAGIC. In an identically conducted analysis of contemporary Crab Nebula data the width of this PSF (defined as the Sigma of a two dimensional Gaussian function) was found to be about 0.08°. Since the Milagro source may be extended, a second analysis was done, assuming an extension of 0.3° instead. Since this is done simply by increasing the signal integration radius, more background events are included, which leads to higher upper limits.

4 Discussion

The derived 95% c.l. flux upper limits are summarized in Table 1. Figures 2 and 3 display them in the context of the published satellite data and the Milagro flux estimation. The differences between the limits of J1954 and J1958 are all compatible with the statistical fluctuations that can be expected for 95% upper limits. Around , where MAGIC is most sensitive, the flux is limited to 3% of the Crab Nebula flux for J1954 and 2% for J1958. Assuming that the Milagro emissions originate from objects spacially coinciding with the Fermi pulsars within our PSF, the photon index in the energy range of to must be harder than 2.2 for J1954, and 2.1 for J1958. The spectral energy distribution is thus likely to peak at energies in excess of .

If an extension up to 0.3° is assumed, the corresponding flux limits in Crab Nebula units are 14% for J1954 and 3% for J1958. In this extended case, the photon indices are limited to and , respectively. It shall be noted that the biggest TeV pulsar wind nebulae have sizes of few tens of parsecs, which at the distance of G65.1+0.6 would be within these 0.3°.

Bringing together the existing flux data and our upper limits, we conclude that the most likely scenario to explain the gamma ray production measured by Milagro might be the existence of two PWN, associated with the Fermi pulsars. With the ages of the pulsars being and , it is reasonable to expect an inverse compton component that dominates their energy outflows and may additionally be extended (de Jager & Djannati-Ataï, 2009; Tanaka & Takahara, 2009). Such old, extended PWN are common TeV gamma-ray sources and frequently have an emission spectrum that peaks at TeV energies or above.

5 Summary

We took 25.5 hours of good quality data in the area of the faint supernova remnant G65.1+0.6. In that region, the Milagro collaboration reported the emission of gamma rays with median energy of in the vicinity of the two GeV Fermi pulsars 1FGL J1954.3+2836 and 1FGL J1958.6+2845. Our observations, which were aimed to locate the source of the Milagro emission, yielded no significant gamma-ray signal for these two a-priori known source locations. Also, no post-trial significant signal from a skymap of the area could be established. We extracted three differential flux upper limits for each source, assuming two different extension radii. They are summarized in Figures 2 and 3.

Assuming that the two 4 and detections of Milagro are not statistical fluctuations, but real signals, the flux upper limits support the scenario in which the multi-TeV emission measured by Milagro is caused by a different mechanism or object than the Fermi emission. Given the ages of the pulsars and the SNR, the existence of two very old pulsar wind nebulae, powered by the two GeV pulsars, seems very likely.

We would like to thank the Instituto de Astrofisica 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 MICINN is gratefully acknowledged. This work was also supported by the Polish MNiSzW Grant N N203 390834, by the YIP of the Helmholtz Gemeinschaft, and by grant DO02-353 of the the Bulgarian National Science Fund.

References

  • Abdo et al. (2009a) Abdo, A. et al. 2009a, ApJ, 700, L127
  • Abdo et al. (2009b) Abdo, A. et al. 2009b, ApJS, 183, 46
  • Abdo et al. (2009c) Abdo, A. et al. 2009c, Science, 325, 840
  • Abdo et al. (2010a) Abdo, A. et al. 2010a, ApJS, 187, 460
  • Abdo et al. (2010b) Abdo, A. et al. 2010b, ApJS, 188, 405
  • Albert et al. (2008a) Albert, J. et al. 2008a, ApJ, 674, 1037
  • Albert et al. (2008b) Albert, J. et al. 2008b, Nucl. Instr. Meth. A, 588, 424
  • Aliu et al. (2009) Aliu, E. et al. 2009, Astropart. Phys.,30, 293
  • Arquilla & Kwok (1987) Arquilla, R., & Kwok, S. 1987, A&A, 173, 271
  • Fomin et al. (1994) Fomin, V. P., Stepanian, A. A., Lamb, R. C., Lewis, D. A., Punch, M., Weekes, T. C. 1994, Astropart. Phys., 2, 137
  • Hartman et al. (1999) Hartman, R.C. 1999, ApJS, 123, 79
  • Landecker et al. (1990) Landecker, T.L., Clutton-Brock, M., & Purton, C.R. 1990, A&A, 232, 207
  • de Jager & Djannati-Ataï (2009) de Jager, O.C., & Djannati-Ataï, A. 2009, Neutron Stars and Pulsars (Berlin, Heidelberg: Springer)
  • Lorimer et al. (1998) Lorimer, D.R., Lyne, A.G., & Camilo, F. 1998, 331, 1002
  • Moralejo et al. (2009) Moralejo, A. et al. 2009, in Proc. 31st ICRC (Łódź), preprint arXiv:0907.0943
  • Rolke et al. (2005) Rolke, W.A., Lopez, A.M., Conrad, J. 2005, Nucl. Instr. Meth. A, 551, 493
  • Saz Parkinson et al. (2010) Saz Parkinson, P.M. et al. 2010, submitted to ApJ, preprint arXiv:1006.2134v1
  • Seiradakis et al. (1985) Seiradakis, J. H., Reich, W., Sieber, W., Schlickeiser, R., & Kühr, H. 1985, A&A, 143, 478
  • Swanenburg et al. (1981) Swanenburg, B.N. et al. 1981, ApJ, 243, L69
  • Tanaka & Takahara (2009) Tanaka, S.J., Takahara, F. 2009, in Proc. 2009 Fermi Symposium (Washington, DC), preprint arXiv:1001.2360
  • Tian & Leahy (2006) Tian, W.W., & Leahy, D.A. 2006, A&A, 455, 1053
Figure 1: Observation setup for the two Fermi sources J1954 and J1958 in the context of SNR G65.1+0.6 and a Milagro significance contour. J1958 appears only in one wobble position (W1), so the OFF data is taken from the other wobble sample, using the same position relative to the pointing direction. The outline of the remnant is taken from the radio map in Landecker et al. (1990). The extension of the Milagro significance contour (Abdo et al., 2009a) is compatible with their point spread function.
Figure 2: Compilation of flux measurements and upper limits for 1FGL J1954.3+2836 from Fermi (Saz Parkinson et al., 2010; Abdo et al., 2010b), MAGIC and Milagro (Abdo et al., 2009a). The 3% fraction of the MAGIC Crab spectrum (Albert et al., 2008a) is shown for comparison.
Figure 3: Compilation of flux measurements and upper limits for 1FGL J1958.6+2845 from EGRET (Hartman et al., 1999), Fermi (Abdo et al., 2010a, b), MAGIC and Milagro (Abdo et al., 2009a). The 3% fraction of the MAGIC Crab spectrum (Albert et al., 2008a) is shown for comparison.
Source Name Assumed Extension Significance
(deg) (GeV) () (10TeVcms) (10TeVcms)
1FGL J1954.3+2836 228 -1.8 27 1.38
942 +1.1 1.22 1.08
5123 +1.9 0.055 1.45
234 -1.3 78 4.3
963 +2.6 5.0 4.66
4956 +2.0 0.157 3.85
1FGL J1958.6+2845 228 -1.5 50 2.60
966 -0.9 0.65 0.60
5123 +0.9 0.063 1.65
234 -0.3 220 12.0
963 -1.9 1.11 1.03
4956 +0.8 0.161 3.9
Table 1: Differential Upper limits
Comments 0
Request Comment
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
   
Add comment
Cancel
Loading ...
229595
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description