A&A 393, L37-L40 (2002)
F. Aharonian1 - A. Akhperjanian7 - M. Beilicke4 - K. Bernlöhr1 - H. Börst5 - H. Bojahr6 - O. Bolz1 - T. Coarasa2 - J. Contreras3 - J. Cortina2 - S. Denninghoff2 - V. Fonseca3 - M. Girma1 - N. Götting4 - G. Heinzelmann4 - G. Hermann1 - A. Heusler1 - W. Hofmann1 - D. Horns1 - I. Jung1 - R. Kankanyan1 - M. Kestel2 - J. Kettler1 - A. Kohnle1 - A. Konopelko1 - H. Kornmeyer2 - D. Kranich2 - H. Krawczynski9 - H. Lampeitl4 - M. Lopez3 - E. Lorenz2 - F. Lucarelli3 - N. Magnussen10 - O. Mang5 - H. Meyer6 - M. Milite4 - R. Mirzoyan2 - A. Moralejo3 - E. Ona3 - M. Panter1 - A. Plyasheshnikov1,8 - J. Prahl4 - G. Pühlhofer1 - G. Rauterberg5 - R. Reyes2 - W. Rhode6 - J. Ripken4 - A. Röhring4 - G. P. Rowell1 - V. Sahakian7 - M. Samorski5 - M. Schilling5 - F. Schröder6 - M. Siems5 - D. Sobzynska2 - W. Stamm5 - M. Tluczykont4 - H. J. Völk1 - C. A. Wiedner1 - W. Wittek2 (HEGRA Collaboration) - Y. Uchiyama11 - T. Takahashi11
1 - Max-Planck-Institut für Kernphysik, Postfach 103980, 69029 Heidelberg, Germany
2 - Max-Planck-Institut für Physik, Föhringer Ring 6, 80805 München, Germany
3 - Universidad Complutense, Facultad de Ciencias Físicas, Ciudad Universitaria, 28040 Madrid, Spain
4 - Universität Hamburg, Institut für Experimentalphysik, Luruper Chaussee 149, 22761 Hamburg, Germany
5 - Universität Kiel, Institut für Experimentelle und Angewandte Physik, Leibnizstraße 15-19, 24118 Kiel, Germany
6 - Universität Wuppertal, Fachbereich Physik, Gaußstr. 20, 42097 Wuppertal, Germany
7 - Yerevan Physics Institute, Alikhanian Br. 2, 375036 Yerevan, Armenia
8 - On leave from Altai State University, Dimitrov Street 66, 656099 Barnaul, Russia
9 - Now at Washington University, St. Louis, MO 63130, USA
10 - Now at IFAE, Unversitat Autònoma de Barcelona, Spain
11 - Institute of Space and Astronautical Science, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 229-8510, Japan
Received 24 July 2002 / Accepted 9 August 2002
Deep observation (113 hrs) of the Cygnus region at TeV energies using the HEGRA stereoscopic system of air Cerenkov telescopes has serendipitously revealed a signal positionally inside the core of the OB association Cygnus OB2, at the edge of the 95% error circle of the EGRET source 3EG J2033+4118, and north of Cyg X-3. The source centre of gravity is RA : , Dec : . The source is steady, has a post-trial significance of +4.6, indication for extension with radius at the level, and has a differential power-law flux with hard photon index of . The integral flux above 1 TeV amounts 3% that of the Crab. No counterpart for the TeV source at other wavelengths is presently identified, and its extension would disfavour an exclusive pulsar or AGN origin. If associated with Cygnus OB2, this dense concentration of young, massive stars provides an environment conducive to multi-TeV particle acceleration and likely subsequent interaction with a nearby gas cloud. Alternatively, one could envisage -ray production via a jet-driven termination shock.
Key words: gamma rays: observations - stars: early-type - Galaxy: open clusters and associations: individual: Cygnus OB2
The current generation of ground-based imaging atmospheric Cerenkov telescopes offer coverage of the multi GeV to TeV -ray sky at centi-Crab sensitivity and arc-minute resolution. Stereoscopy employed by the HEGRA CT-System at La Palma (Daum et al. 1997) offers highly accurate reconstruction of event directions at angles up to off-axis. Results here are taken from data originally devoted to Cyg X-3, and the EGRET source GeV J2035+4214 (Lamb & Macomb 1997). The separation between these objects () permits a combined analysis given the overlap in their CT-System fields of view (FOV). This letter presents analysis details and observational properties of a serendipitously discovered TeV source in these data. A brief discussion concerning astrophysical origin and location of this new source is also presented.
The HEGRA system of imaging atmospheric Cerenkov telescopes (IACT-System), consists of 5 identical telescopes operating in coincidence for the stereoscopic detection of air showers induced by primary -rays in the atmosphere. In data dedicated to Cyg X-3, alternate 20 min runs targeting the Cyg X-3 position in declination were taken during moonless nights of Aug.-Sept. 1999, Sept.-Oct. 2000 and Jun.-Oct. 2001. Likewise in data dedicated to GeV J2035+4124, 20 min runs were obtained tracking directly the GeV source during Jul.-Aug. 2001. In total, three tracking positions are present in combined data. After quality checks, a total of 112.9 hours data are available for analysis. Preferential selection of -ray-like events (against the cosmic-ray background) is achieved by using the difference between the reconstructed and assumed event direction, , and the mean-scaled-width parameter, (Konopelko 1995). In searching for weak point-like and marginally extended sources, so-called tight cuts are considered optimal given the angular resolution of the CT-System (<): and , where we use algorithm "3'' as described by Hofmann et al. (1999) for the event direction reconstruction. The number of images per event, , used for calculating and was also a priori chosen at . Monte Carlo simulations (Konopelko et al. 1999) and tests on real sources have shown that events contribute little to the overall sensitivity.
In searching for new TeV sources, skymaps of event
direction excesses over the RA and Dec plane are generated after having estimated the
background over the FOV. A new empirically-based template background model has
been developed with the goal of simple generation of skymaps.
The template background comprises events normally rejected according
to the criterion. We define the number of events in the -ray regime s from
and for the template background b from
A necessary correction applied
to the template background accounts for differences in radial
profile between the two regimes.
A normalisation ,
to derive excess events
some position in the FOV, accounts for differences
in the total number of events in the two regimes.
A full description of the template model appears in Rowell (2002).
Figure 1 presents the resulting excess skymap.
The template model was used in discovering the TeV source which is evident
north of Cyg X-3. An event-by-event
centre of gravity (COG) calculation (Table 1a),
weighting events with 1 from
the s and regimes
respectively is performed. The COG accuracy is limited by a systematic pointing error of
(Pühlhofer et al. 1997).
A pre-trial significance at the COG position of +5.9
is obtained, summing
(Table 1b). Statistical trial factors arise from the initial "discovery''
skymap (different to that in Fig. 1) in which event directions are independently
summed in 1100 bins of size 0.1
Assuming 1100 trials are accrued in locating the COG, the
for P the pre-trial
or +5.9), is then calculated as
This gives a post-trial significance of +4.6.
1100 is actually a slightly conservative trial estimate since oversampling of the
-ray point spread function (PSF) by a factor 1.5occurs in the discovery skymap.
To verify results using the template model, we make use of
a conventional type of background model employing background regions
displaced from the on-source region spatially in the FOV but derived
from the same
regime. Background events are taken from ring-segments with
matching trigger characteristics to that of the source region.
according to the solid angle ratio between background and on-source regions is then applied.
Results using this so-called ring model (Table 1b) are consistent with those from the template model.
|Figure 1: Skymap (1.5 view at 0.05 binning) of excess events , using the template background model. At each bin, the excess is estimated from events within a radius . Included are 95% error ellipses of various EGRET sources, the core of Cygnus OB2 (Knödlseder 2000), the TeV COG (star), its 2 error circle, and the location of Cyg X-3. ASCA GIS contours (2-10 keV) are overlayed.|
|Open with DEXTER|
|Figure 2: Distribution of for calculated from the COG (solid histogram and filled dots), against a background derived from the template model (dashed hist). The convolved radial Gaussian fit is indicated by the solid line with the intrinsic source size. The PSF width is estimated from Crab data.|
|Open with DEXTER|
Splitting data firstly according to
their three tracking positions reveals commensurate source contributions
(Table 1c). The source is also found
to develop linearly with the cumulative number of background events.
Such tests suggest consistency with a steady source during the three years of data
collection. We have also verified that after cuts a constant background acceptance
throughout the dataset is observed and that the event excess in -space appears consistent
with that of a true -ray population.
To determine source size, we fit a radial Gaussian convolved with the point spread
function (determined from Crab data) to the
excess events as a function of ,
using a subset of events with
the best angular resolution (
for which errors
are minimised (Fig. 2).
The intrinsic size of the TeV source is estimated at
Correlations between the fit parameters suggest that the significance for a non-zero source size
is at the level rather than the 3.5 level indicated above.
A breakdown of the excess with
also shows that the
contributes strongly to the excess
(Table 1d). Such behaviour is
suggestive of a generally hard spectral index given that higher trigger
multiplicities are favoured by higher energy events.
For the energy spectrum calculation and selection cuts, we follow
the method of Aharonian et al. (1999) using
effective collection areas appropriate for on-axis and
off-axis sources as per the
exposure efficiency for the TeV source in these data.
A tight cut
as opposed to the
is also used. For energies below 0.8 TeV the effective collecting area
decreases markedly sharper at positions beyond
off-axis compared to positions nearer on-axis.
Limiting our fit therefore to energies >0.8 TeV reduces systematic
errors. A so-called loose cut
using the ring background model is used in deriving the energy bin-by-bin excess since a loose cut
improves the -ray selection efficiency according
to the moderately-extended nature of the source.
Results are shown in Table 1e and Fig. 3, with the spectrum being well fit
by a pure power law with generally hard photon
index. Systematic errors are estimated from changes in bin centres and uncertainties in Monte Carlo-derived
|Figure 3: Differential energy fluxes of the TeV source and other results. "H-A'' is the AIROBICC 90% confidence level upper limit (Prahl 1999) at the TeV COG converted to differential form at 20.8 TeV assuming a spectral photon index of -2.0. We interpret the 3EG J2033+4118 flux as an upper limit. The ASCA GIS 99% upper limit assumes a photon index of 2, and cm-2. Data fits (solid lines) and their 1 statistical errors in photon index (dashed lines) are shown. TeV data are well fit ( ) by a power law (Eq. (1)).|
|Open with DEXTER|
The OB association Cygnus OB2 is unique for its compact nature and extreme number of member OB and O stars (e.g. Knödlseder 2000), and both theoretical and observational grounds for non-thermal particle acceleration have long-been discussed (e.g. Montmerle 1979; Cassé & Paul 1980; Völk & Forman 1982; White & Chen 1992). The TeV source is positioned inside the core of Cygnus OB2 as defined by Knödlseder (2000). Assuming the TeV source is as distant as Cygnus OB2 (1.7 kpc), a luminosity 1032 erg s-1 above 1 TeV is implied, well within the kinetic energy (KE) budget of Cygnus OB2 estimated recently by Lozinskaya et al. (2002) at a erg s-1, and also within the KE budget of a number of notable member stars (e.g. Massey & Thompson 1991; Manchanda et al. 1996; Benaglia et al. 2001). So far no counterparts at other wavelengths are identified. No massive or luminous Cygnus OB2 star of note discussed recently (e.g. Massey & Thompson 1991; Romero et al. 1999; Herrero et al. 2001; Benaglia et al. 2001) is positioned within the 1 TeV error circle. No catalogued X-ray source from the ROSAT all-sky and pointed survey lies within the 2 TeV error circle. Our analysis of archival ASCA GIS data yields a 99% upper limit (2-10 keV) of erg cm-2 s-1 (Fig. 3). Such results may imply that the energy source for particle acceleration is not co-located with the TeV source, arising instead from the winds of the young/massive stars of Cygnus OB2, either individually or collectively, or from an alternative source. The former scenario would generally favour accelerated hadrons interacting with a local, dense gas cloud, giving rise to -decay TeV emission. The likely hard TeV spectrum can be explained by (energy-dependent) diffusion effects, accelerator age, and accelerator to target distance (see e.g. Aharonian 2001). There is however at present no strong indication from CO and HII surveys (Leung & Thaddeus 1992; Dame et al. 2001; Comeron & Torra 2001) for any co-located dense gas cloud, although see Martí et al. (2000) who discuss a nearby HII cloud in the context that follows immediately below. A suggested alternative scenario involves a jet-driven termination shock at which accelerated electrons produce synchrotron and TeV inverse-Compton (IC) emission (Aharonian & Atoyan 1998). Such a jet could emanate from a nearby microquasar, possibly a class of high energy -ray source (see e.g. Paredes et al. 2000). In fact two nearby sources, 3EG J2033+4118 and also the EGRET source possibly associated with Cyg X-3 (Mori et al. 1997) could be GeV indicators of such a microquasar. Remarkably, Cyg X-3 appears to have a bi-lobal jet (Martí et al. 2000,2001) well-aligned with the TeV source, the latter which would be 70 pc from Cyg X-3 in absolute terms if it is at the same distance (>8.5 kpc). Future X-ray observations will be a crucial constraint on the IC emission in this context. We interpret the flux from 3EG J2033+4118 presently as a MeV/GeV upper limit at the TeV COG (Fig. 3), and note that the directly extrapolated energy flux from the TeV source lies about two orders of magnitude below the 3EG J2033+4118 flux at overlapping energies. Finally we note that earlier claims for a TeV source (Neshpor et al. 1995, at a 1 Crab flux level) and flaring episodes coincident with a Cyg X-3 radio flare at energies >40 TeV (Merck 1993; Krawczynski 1995) positionally consistent with our TeV COG have been reported. These results are however in conflict with our estimates of the flux level and steady nature of the TeV source assuming they all have the same origin. Further observations with the HEGRA CT-System aimed at confirmation and improving our spectral and source morphology studies are now underway.
The support of the German ministry for Research and technology BMBF and of the Spanish Research Council CICYT is gratefully acknowledged. We thank the Instituto de Astrofísica de Canarias for the use of the site and for supplying excellent working conditions at La Palma. We gratefully acknowledge the technical support staff of the Heidelberg, Kiel, Munich, and Yerevan Institutes. GPR acknowledges receipt of a von Humboldt fellowship.