Advanced search for the extension of unresolved TeV sources with H.E.S.S. First measurement of the extension of the Crab nebula at TeV energies
The resolution power of current Imaging Atmospheric Cherenkov Telescopes is presently restricted to scales of a few arcminutes. In the very high-energy (VHE; GeV) gamma-ray regime, the measurement of source sizes that are comparable to or smaller than the resolution of the instrument is usually limited by statistics and in particular by the uncertainties in the characterisation of the instrument Point Spread Function (PSF). The PSF varies strongly with observation and instrument conditions demanding time dependent simulations of these conditions. Employing such simulations, we substantially improve our understanding of the H.E.S.S. PSF and are now able to probe source extensions well below one arcminute scale. We present the results of this new approach applied to known VHE gamma-ray sources and show how this enables us to reveal for the first time the extension of the Crab nebula at TeV gamma-ray energies, with a width of assuming a Gaussian source shape.
Advanced search for the extension of unresolved TeV sources with H.E.S.S.
\FullConference35th International Cosmic Ray Conference - ICRC2017
12-20 July, 2017
Bexco, Busan, Korea
1 Instrument and Dataset
H.E.S.S. is an array of five Imaging Atmospheric Cherenkov Telescopes (IACTs) located in the Khomas Highlands of Namibia. The Crab nebula has been observed with the H.E.S.S. telescopes since the beginning of operations in 2004. For the measurement presented here we only use observations with the four IACTs of the first phase of H.E.S.S. In addition to the standard run-quality selection criteria , we apply additional cuts e.g. on the maximum wind speed (), the observation wobble offset (), and the zenith angle of the observations () to define a high-quality dataset for extension measurements. The resulting Crab nebula dataset consists of observations performed between February 2004 and November 2011, amounting to a total deadtime-corrected live time of h. Since the Crab nebula is a Northern source, the typical observation zenith angles at the H.E.S.S. site in the Southern hemisphere are .
2 Analysis and Results
The data analysis is performed using semi-analytical air-shower
templates . To further improve the angular
resolution, we use a tight cut on the direction reconstruction
uncertainty for the event selection. Furthermore, we demand energies
of TeV to eliminate potential systematic effects near the
energy threshold of the dataset
We show the resulting distribution of events as a histogram in squared angular distance () in the top panel of Fig. 1. The value of is calculcated with respect to the centroid of the gamma-ray excess count distribution in the sky. The centroid position in equatorial coordinates (J2000) is , (systematic error from ), at a distance of from, and within uncertainties compatible with, the Crab pulsar location.
Dedicated run-wise Monte-Carlo (MC) simulations of the dataset, including the actual instrument and observation conditions at the time the data were recorded, are generated as described in . We re-weight the simulations to mimic the energy spectrum of the Crab nebula and analyse them with the same algorithms and analysis configurations as the actual data. The resulting histogram of this MC analysis serves as the PSF for this source and dataset and is shown in the upper panel of Fig. 1. The , , and containment radii of our PSF are , , and , respectively.
As shown in Fig. 1, the PSF is highly inconsistent with the
distribution of the gamma-ray excess counts. The residuals in the
lower panel of Fig. 1 indicate a clear broadening of the
data compared to the simulated PSF. To study this further, we perform
a 2D morphology fit with Sherpa , using the
sky images of gamma-ray-like events around the Crab nebula, of
gamma-ray-like background events estimated from a ring well outside
the source , as well as the simulated
PSF. To determine the best fit, we convolve the PSF with a 2D Gaussian
with different widths. For each width, we calculate a likelihood value
to assess the compatibility of the data and the convolved PSF. We find
the best-fit extension to be
, with a preference of an extension of the Crab
nebula over a point-source assumption of
The resulting best-fit convolution is also plotted in Fig. 1. It clearly provides a good description of the data both in the upper panel and the residuals in the lower panel.
To demonstrate the robustness of our results, we apply the same analysis using time-dependent simulations to two other bright and highly significant extragalactic gamma-ray sources, the active galactic nuclei PKS2155-304 and Markarian421.
As illustrated in Fig. 2, we find both sources compatible with being point-like and show upper limits on their extension. In both cases, these limits are well below the measured extension of the Crab nebula. We note that Markarian421 culminates at large zenith angles of at the H.E.S.S. site, making this source a particularly convincing test of our PSF understanding: even under such challenging observation conditions, Markarian421 appears to be point-like. As we also show in Fig. 2, we tested the Crab nebula dataset for a zenith angle dependence by splitting the observations in two datasets above and below . The measured extensions are compatible with each other.
Also shown in Fig. 2 is the extension of the Crab nebula cross-checked with an independent calibration, reconstruction, and analysis method . We find this second extension measurement slightly larger than our nominal value, and incorporate the difference as one of our systematic uncertainties mentioned above.
Our VHE gamma-ray extension of the Crab nebula is compared to the morphology found at UV wavelengths and X-ray energies in the left and right panel of Fig. 3, respectively. The H.E.S.S. extension covers a good fraction of the optical nebula. Comparing the TeV gamma-ray emission from Inverse Compton scattering to the keV synchrotron X-ray emission, we find that the nebula measured with H.E.S.S. is significantly larger than when measured with Chandra. This result is naturally explained by the radiation cooling of electrons. The Crab nebula size therefore decreases with increasing electron energy, and the energies of electrons producing the TeV gamma rays are well below those of electrons emitting the hard X-rays measured by Chandra. The measured size of the TeV gamma-ray nebula can be reproduced within the standard magnetohydrodynamic model of Kennel and Coroniti [12, 13] assuming a magnetisation parameter .
Here we document the ability of the H.E.S.S. IACT array to robustly measure VHE gamma-ray source extensions down to . The performance boost provided by using time-dependent simulations allows us to resolve for the first time the Inverse Compton component of the Crab nebula. The emission region size we find is well below the previously most constraining upper limit of  and is determined with a high accuracy of . Compared to the synchrotron Crab nebula seen in keV X-rays, the VHE gamma-ray emission region is clearly more extended.
We gratefully acknowledge the support from the agencies and institutions listed here.
- The H.E.S.S. energy threshold for the Crab nebula is relatively large because the source is only visible under large zenith angles.
- is the likelihood ratio test statistic and can be interpreted as statistical significance.
- F. A. Aharonian et al., Observations of the Crab nebula with HESS, Astronomy & Astrophysics 457 (Oct., 2006) 899–915, [astro-ph/0607333].
- M. de Naurois and L. Rolland, A high performance likelihood reconstruction of -rays for imaging atmospheric Cherenkov telescopes, Astroparticle Physics 32 (2009), no. 5 231 – 252, [arXiv:0907.2610].
- T.-P. Li and Y.-Q. Ma, Analysis methods for results in gamma-ray astronomy, The Astrophysical Journal 272 (Sept., 1983) 317–324.
- S. Gillesen, Sub-Bogenminuten-genaue Positionen von TeV-Quellen mit H.E.S.S. Dissertation, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany, 2004.
- M. Holler et al., Run-Wise Simulations for Imaging Atmospheric Cherenkov Telescope Arrays, \posPoS(ICRC2017)755, 2017.
- P. Freeman, S. Doe, and A. Siemiginowska, Sherpa: a mission-independent data analysis application, in Astronomical Data Analysis (J.-L. Starck and F. D. Murtagh, eds.), vol. 4477 of Proceedings of SPIE, pp. 76–87, Nov., 2001. astro-ph/0108426.
- D. Berge, S. Funk, and J. Hinton, Background modelling in very-high-energy -ray astronomy, Astronomy & Astrophysics 466 (May, 2007) 1219–1229, [astro-ph/0610959].
- R. D. Parsons and J. A. Hinton, A Monte Carlo template based analysis for air-Cherenkov arrays, Astroparticle Physics 56 (Apr., 2014) 26–34, [arXiv:1403.2993].
- G. Dubner et al., Morphological Properties of the Crab Nebula: A Detailed Multiwavelength Study Based on New VLA, HST, Chandra, and XMM-Newton Images, The Astrophysical Journal 840 (May, 2017) 82, [arXiv:1704.0296].
- J. Albert et al., VHE -Ray Observation of the Crab Nebula and its Pulsar with the MAGIC Telescope, The Astrophysical Journal 674 (Feb., 2008) 1037–1055, [arXiv:0705.3244].
- F. A. Aharonian et al., Optimizing the angular resolution of the HEGRA telescope system to study the emission region of VHE gamma rays in the Crab Nebula, Astronomy & Astrophysics 361 (Sept., 2000) 1073–1078.
- C. F. Kennel and F. V. Coroniti, Confinement of the Crab pulsar’s wind by its supernova remnant, The Astrophysical Journal 283 (Aug., 1984) 694–709.
- C. F. Kennel and F. V. Coroniti, Magnetohydrodynamic model of Crab nebula radiation, The Astrophysical Journal 283 (Aug., 1984) 710–730.