F. Aharonian¹ - A. G. Akhperjanian² - A. R. Bazer-Bachi³ - M. Beilicke⁴ - W. Benbow¹ - D. Berge^1, - K. Bernlöhr^1,5 - C. Boisson⁶ - O. Bolz¹ - V. Borrel³ - I. Braun¹ - E. Brion⁷ - A. M. Brown⁸ - R. Bühler¹ - I. Büsching⁹ - T. Boutelier¹⁷ - S. Carrigan¹ - P. M. Chadwick⁸ - L.-M. Chounet¹⁰ - G. Coignet¹¹ - R. Cornils⁴ - L. Costamante^1,23 - B. Degrange¹⁰ - H. J. Dickinson⁸ - A. Djannati-Ataï¹² - L. O'C. Drury¹³ - G. Dubus¹⁰ - K. Egberts¹ - T. Eifert^5, - D. Emmanoulopoulos¹⁴ - P. Espigat¹² - C. Farnier¹⁵ - F. Feinstein¹⁵ - E. Ferrero¹⁴ - A. Fiasson¹⁵ - G. Fontaine¹⁰ - Seb. Funk⁵ - S. Funk¹ - M. Füßling⁵ - Y. A. Gallant¹⁵ - B. Giebels¹⁰ - J. F. Glicenstein⁷ - B. Glück¹⁶ - P. Goret⁷ - C. Hadjichristidis⁸ - D. Hauser¹ - M. Hauser¹⁴ - G. Heinzelmann⁴ - G. Henri¹⁷ - G. Hermann¹ - J. A. Hinton^1,14, - A. Hoffmann¹⁸ - W. Hofmann¹ - M. Holleran⁹ - S. Hoppe¹ - D. Horns¹⁸ - A. Jacholkowska¹⁵ - O. C. de Jager⁹ - E. Kendziorra¹⁸ - M. Kerschhaggl⁵ - B. Khélifi^10,1 - Nu. Komin¹⁵ - K. Kosack¹ - G. Lamanna¹¹ - I. J. Latham⁸ - R. Le Gallou⁸ - A. Lemière¹² - M. Lemoine-Goumard¹⁰ - T. Lohse⁵ - R. N. Manchester²⁴ - J. M. Martin⁶ - O. Martineau-Huynh¹⁹ - A. Marcowith^3,15 - C. Masterson^1,23 - G. Maurin¹² - T. J. L. McComb⁸ - E. Moulin^15,7 - M. de Naurois¹⁹ - D. Nedbal²⁰ - S. J. Nolan⁸ - A. Noutsos⁸ - J.-P. Olive³ - K. J. Orford⁸ - J. L. Osborne⁸ - M. Panter¹ - G. Pelletier¹⁷ - P.-O. Petrucci¹⁷ - S. Pita¹² - G. Pühlhofer¹⁴ - M. Punch¹² - S. Ranchon¹¹ - B. C. Raubenheimer⁹ - M. Raue⁴ - S. M. Rayner⁸ - O. Reimer^21, - J. Ripken⁴ - L. Rob²⁰ - L. Rolland⁷ - S. Rosier-Lees¹¹ - G. Rowell^1, - V. Sahakian² - A. Santangelo¹⁸ - L. Saugé¹⁷ - S. Schlenker⁵ $^\star$ - R. Schlickeiser²¹ - F. Schmidt^5, - R. Schröder²¹ - U. Schwanke⁵ - S. Schwarzburg¹⁸ - S. Schwemmer¹⁴ - A. Shalchi²¹ - H. Sol⁶ - D. Spangler⁸ - F. Spanier²¹ - R. Steenkamp²² - C. Stegmann¹⁶ - G. Superina¹⁰ - P. H. Tam¹⁴ - J.-P. Tavernet¹⁹ - R. Terrier¹² - M. Tluczykont^10,23, - C. van Eldik¹ - G. Vasileiadis¹⁵ - C. Venter⁹ - J. P. Vialle¹¹ - P. Vincent¹⁹ - H. J. Völk¹ - S. J. Wagner¹⁴ - M. Ward⁸

1 - Max-Planck-Institut für Kernphysik, PO Box 103980, 69029 Heidelberg, Germany
2 - Yerevan Physics Institute, 2 Alikhanian Brothers St., 375036 Yerevan, Armenia
3 - Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, 9 av. du Colonel Roche, BP 4346, 31029 Toulouse Cedex 4, France
4 - Universität Hamburg, Institut für Experimentalphysik, Luruper Chaussee 149, 22761 Hamburg, Germany
5 - Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
6 - LUTH, UMR 8102 du CNRS, Observatoire de Paris, Section de Meudon, 92195 Meudon Cedex, France
7 - DAPNIA/DSM/CEA, CE Saclay, 91191 Gif-sur-Yvette Cedex, France
8 - University of Durham, Department of Physics, South Road, Durham DH1 3LE, UK
9 - Unit for Space Physics, North-West University, Potchefstroom 2520, South Africa
10 - Laboratoire Leprince-Ringuet, IN2P3/CNRS, École Polytechnique, 91128 Palaiseau, France
11 - Laboratoire d'Annecy-le-Vieux de Physique des Particules, IN2P3/CNRS, 9 Chemin de Bellevue - BP 110, 74941 Annecy-le-Vieux Cedex, France
12 - APC, 11 Place Marcelin Berthelot, 75231 Paris Cedex 05, France
13 - Dublin Institute for Advanced Studies, 5 Merrion Square, Dublin 2, Ireland
14 - Landessternwarte, Universität Heidelberg, Königstuhl, 69117 Heidelberg, Germany
15 - Laboratoire de Physique Théorique et Astroparticules, IN2P3/CNRS, Université Montpellier II, CC 70, Place Eugène Bataillon, 34095 Montpellier Cedex 5, France
16 - Universität Erlangen-Nürnberg, Physikalisches Institut, Erwin-Rommel-Str. 1, 91058 Erlangen, Germany
17 - Laboratoire d'Astrophysique de Grenoble, INSU/CNRS, Université Joseph Fourier, BP 53, 38041 Grenoble Cedex 9, France
18 - Institut für Astronomie und Astrophysik, Universität Tübingen, Sand 1, 72076 Tübingen, Germany
19 - Laboratoire de Physique Nucléaire et de Hautes Énergies, IN2P3/CNRS, Universités Paris VI & VII, 4 place Jussieu, 75252 Paris Cedex 5, France
20 - Institute of Particle and Nuclear Physics, Charles University, V Holesovickach 2, 180 00 Prague 8, Czech Republic
21 - Institut für Theoretische Physik, Lehrstuhl IV: Weltraum und Astrophysik, Ruhr-Universität Bochum, 44780 Bochum, Germany
22 - University of Namibia, Private Bag 13301, Windhoek, Namibia
23 - European Associated Laboratory for Gamma-Ray Astronomy, jointly supported by CNRS and MPG
24 - Australia Telescope National Facility, CSIRO, PO Box 76, Epping NSW 1710, Australia

Abstract
Aims. We present the results of a search for pulsed very-high-energy (VHE) $\gamma$ -ray emission from young pulsars using data taken with the HESS imaging atmospheric Cherenkov telescope system.
Methods. Data on eleven pulsars, selected according to their spin-down luminosity relative to distance, are searched for $\gamma$ -ray signals with periodicity at the respective pulsar spin period. Special analysis efforts were made to improve the sensitivity in the 100 GeV $\gamma$ -ray energy domain in an attempt to reduce the gap between satellite and ground-based $\gamma$ -ray instruments.
Results. No significant evidence for pulsed emission is found in any data set. Differential upper limits on pulsed energy flux are determined for all selected pulsars in the approximate $\gamma$ -ray energy range between 100 GeV and 50 TeV, using different limit determination methods, testing a wide range of possible pulsar light curves and energy spectra.
Conclusions. The limits derived here imply that the magnetospheric VHE $\gamma$ -ray production efficiency in young pulsars is less than 10^-4 of the pulsar spin-down luminosity, requiring spectral turnovers for the high-energy emission of four established $\gamma$ -ray pulsars, and constrain the inverse Compton radiation component predicted by several outer gap models.

Key words: gamma rays: observations - pulsars: individual: Crab pulsar - pulsars: individual: Vela pulsar - pulsars: individual: PSR B1259-63 - pulsars: individual: PSR B1509-58 - pulsars: individual: PSR B1706-44

1 Introduction

Rotating neutron stars are known to convert a significant part of their rotational energy into radiation originating from within the magnetosphere. This emission is observable as a periodic signal at the neutron star rotation frequency (the pulsar phenomenon). For many of the known young and energetic pulsars, the emitted luminosity peaks at X-ray or $\gamma$ -ray energies (see Thompson et al. 1999, and references therein), usually attributed to curvature radiation of accelerated electrons in the strong magnetic fields pervading the pulsar magnetosphere. The luminosity of the pulsed high-energy emission was found to correlate significantly with the energy loss rate of the pulsar, i.e. its spin-down power $\dot{E}$ (Cheng et al. 2004; Dermer & Sturner 1994), which relates to the magnetic field at the neutron star surface, $B_{\rm surf}$ , approximately as $\dot{E}\propto \Omega^4 B^2_{\rm surf}$ , with $\Omega$ as the angular frequency of the pulsar rotation. For most of the pulsars with established $\gamma$ -ray emission (Kuiper et al. 1999; Fierro et al. 1998), there is evidence for a turnover in the pulsed spectrum at a critical energy $E_{\rm c}$ in the sub-GeV to 10 GeV range. Note that no such turnover has been found up to 10 GeV in the case of PSR B1706-44 (Thompson et al. 1996).

Currently, two different scenarios for magnetospheric $\gamma$ -ray emission are discussed, placing the emission regions either near the magnetic poles of the neutron star (polar cap model, Dyks & Rudak 2004; Sturrock 1971; Daugherty & Harding 1996; Ruderman & Sutherland 1975; Rudak & Dyks 1999) or near the null surface in the outer magnetosphere of the pulsar (outer gap model, Cheng et al. 2000; Chiang & Romani 1994; Cheng et al. 1986b; Zhang & Cheng 1997; Romani 1996; Chiang & Romani 1992; Cheng et al. 1986a). Both models predict a cutoff in the curvature radiation spectrum at $\gamma$ -ray energies of the order of GeV up to several tens of GeV. However, the cutoff is of a different physical nature in the two models: while in the case of the polar cap model, the maximum energy of curvature radiation photons escaping the magnetosphere is limited by pair production in the polar magnetic field, the maximum energy of the accelerated electrons determines the maximum photon energy in the outer gap scenario. This results in different predicted light-curves, spectral shapes, and cutoff energies $E_{\rm c}$ of the emitted high-energy radiation in the two models. Note that the predictions of outer gap models generally depend more strongly on model parameters, such as the inclination between the magnetic and rotational axes and the viewing angle to the observer, than is the case for polar cap models. Additionally, in some outer gap model calculations, a spectral component in the TeV range due to inverse Compton (IC) up-scattering of soft ambient seed photons by the accelerated electrons is predicted (e.g. Takata et al. 2006; Hirotani et al. 2003; Hirotani & Shibata 2001b). The boosted TeV photons, the flux of which is determined by the seed photon density, are more likely to escape due to the lower magnetic field strength in the outer magnetosphere compared to the polar cap regions.

Bridging the space between the distinct polar cap and outer gap models, particle acceleration within the slot gap (Arons 1983), a narrower region between the last open magnetic field line and the pair-plasma region of the polar cap, has been discussed in the context of high-energy radiation. However, so far slot gap calculations have been used primarily to model details of pulsed high-energy light-curves (Dyks & Rudak 2003; Muslimov & Harding 2003) without reproducing the corresponding energy spectra in the GeV range or predicting spectra for TeV energies.

Pulsars with established $\gamma$ -ray emission components have been subject to intensive searches for pulsed very-high-energy (VHE, energies above $\sim$ 100 GeV) $\gamma$ -ray emission by ground based instruments. Up to now, no evidence for pulsed emission has been found in these observations (Yoshikoshi et al. 1997; Chadwick et al. 2000; Aharonian et al. 2004a; de Naurois et al. 2002; Lessard et al. 2000), and upper limits on the pulsed VHE $\gamma$ -ray flux have been derived under various assumptions on the characteristics of the pulsed emission. However, the IC component predicted by outer gap models has not yet been significantly constrained.

The superior sensitivity of the HESS detector with respect to previous instruments (see e.g. Ong 2006), combined with its Southern Hemisphere location, puts the predicted pulsed IC component from outer gap models within reach of testability, motivating a new search for pulsed VHE emission from the known $\gamma$ -ray pulsars. A significant amount of HESS data were accumulated in a survey of the inner Galaxy (Aharonian et al. 2006d), resulting in the detection of numerous new sources of VHE $\gamma$ -rays, many of which are apparently associated with pulsars. Thus, in addition to the known $\gamma$ -ray pulsars, these data can be used to search for pulsed emission from promising young, energetic pulsars.

The paper is structured as follows: in Sect. 2 the selection criteria of the studied pulsar sample are explained, Sect. 3 introduces the HESS instrument and data analysis, in Sect. 4 the methods employed in the search for variability and the flux limit determination are explained, and finally, the results are presented and discussed in Sects. 5 and 6, respectively.

2 Pulsar candidate selection

The prime candidates for the search for VHE $\gamma$ -ray emission are the pulsars with established $\gamma$ -ray emission at energies below $\sim$ 10 GeV which have been detected by CGRO instruments. Four of them, PSR B0531+21 (Crab pulsar), PSR B0833-45 (Vela pulsar), PSR B1706-44, and PSR B1509-58 were the subject of pointed HESS observations and many hours of observations were obtained for each of them in the years 2002 to 2005. In addition, a significant number of known pulsars were observed as part of the HESS survey of the Galactic Plane in 2004 and 2005, or lie within the field of view of observations of other HESS targets. The location of the HESS site in the southern hemisphere limits the declination of observations to < $30\hbox{$^\circ$ }$ .

Apart from the known $\gamma$ -ray pulsars, other candidates for which HESS data were available were selected from the ATNF pulsar catalogue (see Manchester et al. 2005, and references therein) if their spin-down flux $\dot{E}/D^2$ was greater than $10^{35}~\rm erg~s^{-1}~kpc^{-2}$ . Table 1 lists all candidates chosen along with selected measured and derived characteristics collected from the literature.

Table 1: The characteristics of the selected pulsars taken from Manchester et al. (2005). Period, P, distance, D, spin-down age, spin-down luminosity, $\dot{E}$ , and the corresponding value for $\dot{E}/D^2$ , and calculated magnetic field strength at the neutron star surface, $B_{\rm surf}$ , and the light cylinder, $B_{\rm LC}$ , are listed. The last column shows the rank in $\dot{E}/D^2$ within the ATNF catalogue.

Some peculiarities of individual pulsars are pointed out in the following. PSR J1747-2958, situated in a bow-shock pulsar wind nebula (PWN) is detected as a compact X-ray source (Gaensler et al. 2004). The only binary pulsar in the present sample, PSR B1259-63 orbits a Be-star companion. Its magnetic field is interestingly low ( $\sim$ 3 $\times$ 10¹¹ G), suggesting a high cutoff energy in polar cap models. It has been discovered as an emitter of transient X-rays (Cominsky et al. 1994) and VHE $\gamma$ -rays (Aharonian et al. 2005), however, pulsed high-energy emission has not so far been detected. PSR J1420-6048 is embedded in a PWN and identified with a compact X-ray source (Ng et al. 2005) within the Kookaburra complex, from which VHE $\gamma$ -ray emission was recently detected (Aharonian et al. 2006a). PSR J1811-1925, for which pulsed emission has been detected up to hard X-rays (Gavriil et al. 2004), is located in the composite supernova remnant G11.2-0.3. All of the selected pulsars are visible at radio wavelengths with the exception of PSR J1811-1925 which is only seen in X-rays.

All available HESS data were searched for the respective pulsar to be in the field of view with a maximum angular distance of $\theta=2$ $^\circ$ between the pulsar position and the HESS pointing direction.

3 VHE gamma-ray detection and data analysis

3.1 The HESS detector

The High Energy Stereoscopic System (HESS, Hofmann 2003), an array of imaging atmospheric Cherenkov telescopes located in Namibia at $23\degr16\hbox{$^\prime$ }18\hbox{$^{\prime\prime}$ }$ S $16\degr30\hbox{$^\prime$ }00\hbox{$^{\prime\prime}$ }$ E and 1800 m a.s.l., detects cosmic VHE $\gamma$ -rays by imaging the Cherenkov emission of their air showers in the atmosphere using optical telescopes. Each telescope has a tessellated spherical mirror with 13 m diameter and $107~\rm m^2$ area (Cornils et al. 2003; Bernlöhr et al. 2003) and is equipped with a camera of 960 $0.16\hbox{$^\circ$ }$ -photomultiplier tubes providing a total field of view of $5\hbox{$^\circ$ }$ in diameter (Vincent et al. 2003). During the stereoscopic observations, an array trigger requires the simultaneous detection of air showers by several telescopes at the hardware level (Funk et al. 2004), allowing a suppression of background events and the assignment of GPS timestamps to each triggered air shower, without the need of an offline synchronisation of the individual telescopes. After digitisation, the event data are stored on disk by the data acquisition system (Borgmeier et al. 2003).

Following the calibration of the recorded air shower data (Aharonian et al. 2004b), each telescope image was parametrised by its centre of gravity and second moments (Hillas 1985) followed by the stereoscopic reconstruction of the shower geometry providing an angular resolution of $\sim$ $0.1\hbox{$^\circ$ }$ for individual $\gamma$ -rays. The $\gamma$ -ray energy was estimated from the image intensity and the shower geometry with a typical resolution of $\sim$ $15\%$ . In order to reject the vast background of cosmic-ray showers, $\gamma$ -ray candidates are selected using cuts on image shape scaled with their expectation values obtained from Monte Carlo simulations. Details of the analysis technique can be found in Aharonian et al. (2006c). The selection cuts used for this analysis were optimised on simulations of a $\gamma$ -ray point source with 10% of the flux of, and the same spectral shape as, the VHE $\gamma$ -ray emission from the Crab Nebula. These cuts, in the following called standard, allow the detection of sources with 1% Crab Nebula flux at energies above $\sim$ 200 GeV with a significance of 5 $\sigma$ within 25 h at a zenith angle of 20 $\hbox{$^\circ$ }$ . Note that the threshold energy above which $\gamma$ -ray sources can be detected increases with the zenith angle of observation.

Table 2: Selection cuts used for the analysis of air shower events. The image amplitude (Size) is given in photo-electrons (p.e.), $E_{\rm max}$ corresponds to the maximum reconstructed event energy, MRSW and MRSL represent the mean scaled reduced width and length parameters, and $\theta ^2_{\rm cut}$ is the maximum squared angular distance with respect to the target position. The tail cuts represent the image cleaning thresholds applied. For details see Aharonian et al. (2006c).

3.2 Low energy optimisation

Since observational data indicate steep cut-offs in high-energy (GeV) $\gamma$ -rays, special low-energy cuts have been applied to reduce the gap in observational coverage between satellite and ground based $\gamma$ -ray observations of young pulsars. Table 2 lists the values of the two samples of selection cuts used in this analysis. See Aharonian et al. (2006c) for a detailed description of the individual cuts. The low energy cuts were optimised on simulated $\gamma$ -rays to yield the best sensitivity for steep spectrum sources ( ${\rm d}N/{\rm d}E\propto E^{-\Gamma}, \Gamma=5$ ) with energies below 500 GeV. Figure 1 shows the rate of simulated $\gamma$ -ray events for a Crab-like source spectrum ( $\Gamma=2.6$ ) passing the standard (dark grey) and low energy (light grey) selection cuts as a function of $\gamma$ -ray energy for a typical zenith angle of 20 $^\circ$ . The peak of the distribution ${\rm d}N/{\rm d}E$ is commonly defined as the threshold energy. The low energy cuts result in an energy threshold lowered by a factor of 1.5 with respect to the standard cuts. Although the number of selected events is higher by 50%, the background rejection efficiency of the low energy cuts is worse by more than a factor of 3. Thus, the extension of the detectable $\gamma$ -ray energy range is achieved only at the expense of decreased sensitivity, especially at higher energies.

$\begin{figure} \par\includegraphics[width=8.35cm,clip]{6737f1.eps}\end{figure}$	Figure 1: Event rate of simulated $\gamma$ -rays as a function of $\gamma$ -ray energy for two different sets of selection cuts. The dark grey histogram corresponds to the standard and the light grey histogram to the low energy cuts. The vertical arrows roughly represent the corresponding threshold energy.
Open with DEXTER

In order to quantify possible systematic effects of the low energy cuts on the energy spectrum determination, data from observations of the Crab Nebula were used to reconstruct the differential energy spectrum for both sets of cuts. In Fig. 2 the open points show the spectrum of unpulsed $\gamma$ -rays determined with the standard cuts. The solid line indicates the best-fit power law with exponential cutoff taken from Aharonian et al. (2006c), Table 5. The full points result from the analysis with low energy cuts, showing minor systematic differences in the absolute flux normalisation (<10%).

$\begin{figure} \par\includegraphics[width=8.35cm,clip]{6737f2.eps}\end{figure}$

Figure 2: Differential energy spectrum of unpulsed VHE $\gamma$ -rays from the Crab Nebula obtained with the different sets of selection cuts. The open points correspond to the results of the standard and the full points to the low energy cut analysis (but with $E_{\rm max}=100$ TeV as for the standard cuts). The solid line represents the spectrum parametrisation taken from Aharonian et al. (2006c).

Open with DEXTER

3.3 Point source analysis

Observation runs are selected according to standard quality selection criteria, ensuring stable detector and atmospheric conditions as described in Aharonian et al. (2006c). The data were analysed using the standard method for the search for point sources without considering any time variability. The number of events in a circular region around the source position with radius $\theta_{\rm cut}$ is compared to the number of counts in several background regions of the same size, arranged on a circle around the centre of the field of view with radius $\theta_{\rm off}$ . For each run, the background normalisation is determined by the number of background regions used. For data sets comprised of several different offsets, as is the case for most of the selected pulsars, the background normalisation for each run was weighted by its relative acceptance for a given offset. This acceptance was derived from data of HESS observations without significant $\gamma$ -ray emission within the field of view (see also Berge et al. 2006).

4 Periodicity analysis

In order to test for pulsed $\gamma$ -ray emission at the pulsar position, the timestamps of each recorded shower passing selection cuts were transformed from the observer's frame into the pulsar frame and then folded with the pulsar spin period. The resulting unbinned distribution of pulsar phases corresponding to each shower event was tested for variability. In cases where no significant periodicity was found, flux limits were determined under various assumptions concerning pulse shape and position.

4.1 Timing corrections

The time of arrival (TOA) for each shower event recorded by the telescope array is determined with a central trigger GPS clock (Funk et al. 2004), achieving an accuracy of better than $2~\mu$ s. Each TOA $t_{\rm TOA}$ was transformed to the pulsar frame using the following corrections:

4.2 Tests for pulsed emission

The corrected TOA was subsequently folded with the timing parameters of the target pulsar (pulsar ephemeris) obtained from radio or X-ray observations. For each event TOA t, the pulsar phase $\phi$ was determined according to

Table 3: Pulsar ephemerides. The first and last reference times (t₀) of the timing solutions used in this analysis are given. For a complete list of the solutions and the corresponding pulsar parameters see the references: 1: Lyne et al. (2005), 2: Manchester et al. (2005), 3: Camilo et al. (2002).

For the X-ray pulsar PSR J1811-1925 and the radio pulsars PSR J1801-2451 and PSR J1747-2958 there were no contemporaneous observations available from which the pulsar timing parameters could have been deduced. For these pulsars, the most recently published timing solution was used, assuming that the pulsar phase determined with the known parameters is accurate enough to reveal the periodicity of a potential signal. In particular, this assumes that no glitches or other major frequency changes occurred within the time range between the HESS observations and the measurement of the timing parameters.

The unbinned phase distribution of events was tested for significant periodicity, i.e. for a significant deviation from a flat distribution, using several statistical tests. The Pearson $\chi ^2$ -test tests for a constant distribution of phases in a number of phase bins, which was set to 20 in order to give acceptable coverage for any of the considered data sets. The Z²_m-test (Buccheri et al. 1983) checks for a combination of fundamental sine and cosine harmonics up to the m-th order and is independent of any binning. Here, the Z²₁- and Z²₂-tests for single and double peak profiles are employed. Additionally the H-test (de Jager et al. 1989), a more powerful generalisation of Z²_m, was performed to increase the sensitivity to unknown pulse profiles. Finally, a rather different approach which checks for maximum deviations from the uniform phase distribution, the Kuiper-test (Jetsu & Pelt 1996), was used.

In order to test for spurious signals induced by possible systematic effects in the HESS event data, the tests were also applied to the phase series of background events, i.e. events from sky regions used in the background determination for the unpulsed $\gamma$ -ray excess (cf. Sect. 3.3).

As an example, the distribution of event phases from observations of the Vela Pulsar (PSR B0833-45) is shown in Fig. 3, obtained using the standard cuts. The difference between on and off results from the known $\gamma$ -ray excess from HESS J0835-456 at the position of the pulsar (see Table 5). No significant deviation from uniformity was found with any of the statistical tests for pulsations.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f3.eps}\end{figure}$

Figure 3: Upper plot: distribution of event phases for the Vela pulsar (PSR B0833-45). The points represent the events in the on-region at the pulsar position and the histogram the normalised off-region events. The solid and dashed lines represent the $\chi ^2$ fits to a constant for the on- and off-events, respectively. For both regions the probabilities for being consistent with a uniform distribution according to the statistical tests on pulsations (see text) are listed along with the pulsed fraction limit as introduced in Sect. 4.3. Lower plot: phase distribution for $\gamma$ -rays with energies between 2 and 10 GeV as measured by EGRET (Fierro et al. 1998). The dotted vertical lines denote the edges of the on-pulse regions.

Open with DEXTER

4.3 Flux limit determination

Several methods were applied to obtain limits on the $\gamma$ -ray flux from the selected pulsars. They differ in the assumptions made concerning the characteristics of the pulsed emission, in particular concerning the form of the light-curve and the phase position.

In the most conservative approach, without any assumptions on the temporal structure or spectrum of a possible signal, the energy flux $E^2{\rm d}N/{\rm d}E$ is calculated using the standard analysis as described in Sect. 3.3. This method considers all events in the direction of the pulsar, i.e. in the on-region, ignoring any possible variations and is referred to as the unpulsed method. The flux is calculated event-wise for the on- and off- (background) regions within the same field of view separately, accumulated in bins of energy, and subsequently subtracted to yield the excess flux (Aharonian et al. 2005c). Note that due to the finite width of the energy bins, the reconstructed flux is weakly dependent on the $\gamma$ -ray energy spectrum for which the detector efficiency was calculated. Finally, the flux obtained in this way is converted into a 99% confidence level upper limit using the method of Feldman & Cousins (1998). As some of the pulsars in this sample show significant unpulsed $\gamma$ -ray emission, the flux limits derived can be considerably higher than the nominal HESS $\gamma$ -ray flux sensitivity since the unpulsed $\gamma$ -ray emission is an additional background to the pulsed emission.

The second method takes advantage of the pulse position and shape as measured in other energy domains, assuming similar characteristics for the hypothetical VHE $\gamma$ -ray emission. In particular for the known $\gamma$ -ray pulsars this approach seems reasonable, since the $\gamma$ -ray emission mechanism is expected to be similar at VHE $\gamma$ -ray energies, especially for photons originating from curvature radiation. For each pulsar, a known pulse profile was chosen according to the energy band of the emission, taking the highest photon energies. For these pulse profiles, an on- and an off-phase domain were assigned (on-off-pulse method), representing the pulsar phases for high and low (or no) emission, respectively. Table 4 lists the used on-phase intervals and their corresponding energy bands for all pulsars. The phase intervals are given relative to the reference phase which is coincident with the highest peak of the respective radio profile. For example for the Vela pulsar, the on- and off-phase domains were chosen based on the light-curve measured using EGRET at GeV energies (Fierro et al. 1998), as indicated by the dotted lines in Fig. 3. Note that any possible contribution from inter-pulse domains, the phase range between the on-pulse domains, to the pulsed signal are neglected in this method. For the pulsars without established $\gamma$ -ray emission, the selected phase region always covers 20% of the profile around the reference phase. The VHE $\gamma$ -ray flux was then calculated as for the unpulsed method, but using the off-phase events in the direction of the pulsar for background estimation, instead of the events in the off-regions within the field of view. If there was no significant unpulsed signal visible within a data set and no known VHE $\gamma$ -ray source coincident with the target pulsar, i.e. for PSR J1524-5625, PSR B1706-44, and J1801-2451, the off-event statistics were increased by including those events which were in the off-regions of the field of view and additionally had phases falling in the assigned off-phase domain. In these cases the background normalisation was adjusted accordingly. The on-off-pulse method is less reasonable for those pulsars for which pulsed emission is known only at low energies, especially at radio wavelengths, since the $\gamma$ -ray production mechanism may be completely unrelated to that of the low energy emission.

Therefore, a third approach was also applied, that does not rely on knowledge of the phase position of the pulse, but only assumes a certain pulse shape. Here, a combination of first and second order harmonics was chosen, similar to the Z²₂-test. A 99% confidence level upper limit on the fraction, p $\ll$ 1, of pulsed emission present in the series of event phases for events with shower directions from within the on-region around the pulsar position was calculated for a given energy bin according to Brazier (1994). The resulting event fraction was converted into a flux using the effective collection area of the detector and assuming a spectral distribution in energy following a power law with photon index $\Gamma$ = 2.7, roughly matching the shape of the cosmic ray background spectrum. More details on the flux calculations can be found in Aharonian et al. (2006c). This method, in the following referred to as the pulsed fraction method, is less constraining than the on-off-pulse method, since no background subtraction can be performed and thus the post-cut cosmic ray background is treated as an unpulsed component in addition to any unpulsed VHE $\gamma$ -ray excess.

Altogether, the pulsed fraction method was considered to be the most robust approach since the structure of the VHE $\gamma$ -ray light curve remains unknown a priori. Nevertheless, the on-off-pulse method was still applied to ensure comparability with results from prior searches for pulsed VHE $\gamma$ -ray emission where very similar techniques have been used. For the three pulsars for which valid timing parameters were not available for the HESS data, PSR J1747-2958, PSR J1801-2451 and PSR J1811-1925, flux limits have been derived only using the unpulsed method.

5 Results

5.1 Point source analysis

Table 5 shows the main characteristics of the data sets for each pulsar along with the results of the point source analysis for unpulsed emission. Additionally, any HESS VHE $\gamma$ -ray source close to the corresponding pulsar position is listed whenever a significant $\gamma$ -ray excess was detected which can be attributed to this known source. Note that the excess given may differ significantly from that given in the respective reference, since in all cases, with the exception of the Crab nebula, the VHE emission region is both extended and offset from the pulsar position.

Table 5: Results of the analysis of unpulsed $\gamma$ -ray emission for each pulsar data set using the standard cuts and assuming a point source at the pulsar position. Dead-time corrected observation time, $t_{\rm obs}$ , mean offset of the pulsar direction with respect to the array pointing direction, $\langle\theta_{\rm off}\rangle$ , approximate mean energy threshold, $\langle E_{\rm thresh}\rangle$ , number of excess counts, $N_{\rm excess}$ , and corresponding significance, S, are given. In addition to the mean threshold for the standard cuts, the corresponding threshold for the low energy cuts is given in parentheses if applicable. If the pulsar is in the vicinity (angular distance <1 $^\circ$ ) of, or associated with, a known VHE $\gamma$ -ray source, the source name is given along with the corresponding reference. Note that the given numbers can significantly differ from the published results for the known sources which are located at different positions and might be extended. References: A: Aharonian et al. (2006c), B: Aharonian et al. (2006b), C: Aharonian et al. (2005), D: Aharonian et al. (2006a), E: Aharonian et al. (2005b), F: Aharonian et al. (2005d), G: Aharonian et al. (2006d), H: Aharonian et al. (2005a).

Table 6: Results of statistical tests for pulsed emission. Number of events, $N_{\rm on}$ , $\chi ^2$ probability for a fit of the phase distribution to a constant, and probabilities for the event series being compatible with a flat phase distribution for the H-, Z²₁-, Z²₂-, and Kuiper-test (P_K) are shown. All quantities are given for both the low energy and standard cut analyses (see text).

Three of the four pulsars not coincident with known VHE $\gamma$ -ray sources - PSR J1524-5625, PSR B1706-44, and PSR J1801-2451 - show no significant signal, but for PSR J1811-1925 an excess with 3.8 $\sigma$ significance is seen. After taking into account all four trials for these target positions, the significance is reduced to 3.5 $\sigma$ .

Since the HESS energy threshold increases with the zenith angle of observation, the low energy selection cuts introduced in Sect. 3.2 effectively reject all data taken at zenith angles higher than 50 $^\circ$ since the effective collection area of the instrument vanishes below 500 GeV $\gamma$ -ray energy for these elevations. Therefore the effective exposure time for the low energy analysis is reduced for some data sets compared to that of the standard cuts.

For PSR B1706-44, only archival data from the early commissioning phase of HESS from observations with two telescopes and without an array level trigger are available. Therefore, event timing was determined from the GPS timestamps of the individual cameras. Additionally, the low energy analysis was not performed.

In the special case of PSR B1259-63, the data set used here differs from that used in Aharonian et al. (2005). Newer data from observations in 2005 are included and data taken before MJD 53135 are excluded. Before this date reliable radio ephemerides could not be provided as the radio pulsar was eclipsed by its binary companion. Additionally, the unpulsed $\gamma$ -ray emission, presumably arising from the interaction of the pulsar wind and the stellar outflow, would result in a strong background for any pulsed emission component.

5.2 Search for pulsations

Table 6 shows the results of the statistical tests applied to the phase distributions of events in the on-region around each pulsar position. The probabilities for the hypothesis of a flat phase distribution, i.e. the probability that no pulsed signal exists, are given for the two sets of selection cuts applied. It can be seen that these probabilities lie well above the one per cent level with a single exception. For the low energy analysis of PSR J1420-6048, the test probabilities seem significantly low, especially for the more powerful H- and Kuiper-tests. No such effect is seen for the respective probability of the background distributions of any data set, which are all above 1%. In Fig. 4, the corresponding phase distribution is shown, again for both on and off-events. The off-event distribution is much more compatible with a flat distribution, with probabilities well above 10%. The probability $P_{\rm H}\vert _{\rm PSR~J1420-6048} = 3.1$ $\times$ 10^-3 of the H-test for the on-events corresponds to $3.0\sigma$ . However, this probability has to be corrected for the number of trials for this test. Considering 19 trials for all data sets analysed with the low-energy (9 trials) and standard (10 trials) selection cuts respectively, the corresponding significance for the applied tests lies between $2.0\sigma$ (H-test) and $0.4\sigma$ (Z²₁-test) with an average of $1.5\sigma$ . Note that the test results for the pulsars PSR J1747-2958 and PSR J1801-2451 for which no valid timing information was available do not allow firm conclusions about the presence of periodicity within the HESS data with respect to the intrinsic pulsar rotation period. For PSR J1811-1925, the tests were not performed since there was no valid timing information available for any period within several years of the HESS observation time.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f4.eps}\end{figure}$	Figure 4: Distribution of event phases for PSR J1420-6048 using the low energy cut analysis for on- and off-regions (see also Fig. 3). Note that there is significant unpulsed emission at the pulsar position resulting in a difference between the mean number of events for the on and off-regions shown by the solid and dashed lines respectively.
Open with DEXTER

As a cross-check for systematic effects possibly biasing the statistical tests, the ensemble distribution of the test statistic was checked against the hypothesis that no pulsed signal exists in any of the data sets. Figure 5 shows the distribution of the test statistic H of the H-test for the on- and off-regions and for the different sets of cuts. Also shown is the expected distribution $N_{\rm H}(H) = N_0 \exp(-\lambda H)$ for the case when no pulsed signal is present in any data set, with N₀ depending on the number of tests and $\lambda=0.4$ . In order to quantitatively evaluate the match between data and expectation, each of the distributions for on- and off-events and for the two sets of cuts was fitted by the respective expected distribution with a log-likelihood method, leaving N₀ and $\lambda$ as free parameters, and the probability for being compatible with the expected parameters was calculated. No significant deviation is found for any of the four distributions with all probabilities above 20%.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f5.eps}\end{figure}$

Figure 5: Ensemble distributions of the H-test statistic for the selected pulsars of Table 6 and their corresponding background control samples. The results for the pulsar direction are shown as open and closed circles for the low energy and standard cut analysis, respectively. The distributions for the off-regions are displayed as grey filled and outlined histograms, respectively. The solid curve shows the expected distribution if no pulsed signals are present.

Open with DEXTER

5.3 Flux upper limits for pulsed emission

Differential energy flux upper limits (i.e. in $E^2{\rm d}N/{\rm d}E$ ) were derived for each data set, for the energy range from threshold energies up to the highest energies for which significant event statistics were accumulated. In Figs. 6-11 the differential flux limits calculated are shown separately for each pulsar. Different points are shown for the different limit determination methods (see Sect. 4.3). The points above energies of 0.5 TeV correspond to the results for the standard selection cuts and the points below 0.5 TeV to these for the low energy cuts. The horizontal position of the points denotes the mean energy of the events in each energy bin while the bin edges are indicated by the horizontal error bars. In the case of the pulsed fraction method, the mean energy was calculated for the assumed power law energy spectrum.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f6.eps}\end{figure}$

Figure 6: HESS energy flux limits (99% c.l.) for pulsed emission of the Crab pulsar. The full circles and full squares correspond to the on-off-pulse and pulsed fraction limit determination methods, respectively. Below energies of 0.5 TeV the results were obtained with the low energy selection cuts, otherwise the standard cuts were used. Note that the threshold energy of the HESS data lies at 310 GeV due to the observation at high zenith angles imposed by the limited accessibility of the northern sky. The crosses, open squares, and open diamonds denote limits from CELESTE (99% c.l., de Naurois et al. 2002), Whipple (99.9% c.l., Lessard et al. 2000), and HEGRA (99.865% c.l., Aharonian et al. 2004a), respectively. The open triangles show the EGRET phase-averaged spectrum (Fierro et al. 1998). The indicated polar cap curve was generated according to Rudak & Dyks (1999) and the outer gap model curve taken from Hirotani & Shibata (2001a).

Open with DEXTER

Note that the overall $\gamma$ -ray energy scale for the standard cuts has a systematic error estimated to be 20% for a power law energy spectrum with $\Gamma$ $\sim$ 2.6 (Aharonian et al. 2006c), dominated by the absolute shower energy calibration of the instrument. The impact of the systematic uncertainty in the energy reconstruction method on the flux upper limit for the low-energy cuts was determined to be 30% for energies around the energy threshold of the respective data set, mainly affecting the given flux limit point for the lowest energy bin.

6 Discussion

For the four $\gamma$ -ray pulsars, these new results allow some general conclusions on high-energy emission mechanisms of pulsars to be drawn. The results are put in the context of model predictions from selected polar cap and outer gap scenarios. In the particular case of the outer gap scenario, only those model calculations were considered which make explicit predictions for the IC component at TeV energies. Where appropriate, simulated spectra are shown which are calculated using simulation codes from Rudak & Dyks (1999) and Hirotani et al. (2003). Finally, the results can be compared with findings of other experiments in the TeV energy domain.

Crab pulsar. The Crab pulsar has been studied by all major ground-based $\gamma$ -ray instruments at energies between 60 GeV (de Naurois et al. 2002) and 50 TeV (Aharonian et al. 2004a). Figure 6 shows the new HESS results along with $\gamma$ -ray measurements from EGRET and selected results in the VHE domain. The southern location of HESS allows only observations at rather high zenith angles, prohibiting a deep exposure especially at low energy thresholds. Therefore the limits obtained using the Northern Hemisphere instruments HEGRA, with more than 200 h of observations, and Whipple, with a lower energy threshold (Lessard et al. 2000), could not be improved upon. However, the HEGRA limits are derived only for a rather narrow pulse region between -0.06 and 0.04 of the main pulse, whereas the Whipple and HESS limits include both the main and secondary pulse regions. Taking the HESS results for the on-off-pulse method and assuming a power-law spectrum with exponential cutoff ${\rm d}N/{\rm d}E\propto E^{-\Gamma}\exp(-E/E_{\rm c})$ , with $\Gamma=2.05$ , the cutoff energy is constrained to be $E_{\rm c}$ < 80 GeV, less constraining than the flux limit at $\sim$ 60 GeV reported by CELESTE (de Naurois et al. 2002). Our new differential flux upper limits for the Crab pulsar constrain the IC component in the TeV range predicted by Hirotani & Shibata (2001a), reinforcing the constraints placed by previous $\gamma$ -ray experiments, both at GeV and TeV energies. No constraints can be placed on the polar cap prediction, which exhibits a sharp cutoff at a few GeV.

Vela pulsar. For Vela, the limits set using the Durham Mark 6 telescope (Chadwick et al. 2000) and the CANGAROO (Yoshikoshi et al. 1997) detector are improved by more than one order of magnitude in flux and reach down to $\sim$ 100 GeV (see Fig. 7), constraining the cutoff energy for the EGRET power law with $\Gamma=2.38$ to be $E_{\rm c}<14$ GeV. Additionally, the HESS limits significantly constrain the parameter space for the IC component of the outer gap model (Hirotani et al. 2003), which was computed with a slight change of the gap size $D_\perp=0.177$ (see Hirotani et al. 2003, for details) to better match the EGRET data. Note that in the model calculation of Hirotani et al. (2003), no contribution from the photons of the optical pulse to the outer gap seed photons was considered (Hirotani 2006) which would result in a further increased IC component. On the other hand, a more detailed treatment of the outer gap anisotropy would lead to a decreased IC contribution (Hirotani 2006). A newer, two-dimensional version of the same outer gap model, which includes the contribution from the optical and infrared pulse, predicts a lower IC flux in the TeV range (Takata et al. 2006). The polar cap model again predicts pulsed emission only at energies far below the energy threshold of current ground-based $\gamma$ -ray detectors.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f7.eps}\end{figure}$

Figure 7: HESS energy flux limits (99% c.l.) for pulsed emission of the Vela pulsar (see Fig. 6 for point descriptions). The cross and open square denote limits from the Durham Mark 6 telescope (99.865% c.l., Chadwick et al. 2000), and the CANGAROO detector (95% c.l., Yoshikoshi et al. 1997), respectively. The open triangles resemble the EGRET phase-averaged spectrum (Fierro et al. 1998). The polar cap and outer gap model curves are generated according to the model of Rudak & Dyks (1999, black dotted curve) and from Hirotani et al. (2003, dashed curve), respectively. The dotted grey outer gap model curve is taken from Takata et al. (2006).

Open with DEXTER

PSR B1706-44. For this pulsar, limits on pulsed $\gamma$ -ray emission were previously reported by Chadwick et al. (2000). Although HESS data were mostly limited to 2-telescope data from the commissioning phase to which only the on-off-pulse method has been applied, more restrictive limits could be set (see Fig. 8), constraining the exponential cutoff energy to be $E_{\rm c}<71$ GeV using the power law fit of the EGRET spectrum with $\Gamma=2.25$ . However, the IC component in the outer gap scenario of Hirotani et al. (2003) is left unconstrained by several orders of magnitude.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f8.eps}\end{figure}$

Figure 8: HESS energy flux limits (99% c.l.) for pulsed emission of PSR B1706-44 for the on-off-pulse method (full points). The cross denotes a flux limit from the Durham group (99.865% c.l., Chadwick et al. 2000). The open triangles show the EGRET phase-averaged $\gamma$ -ray spectrum (Thompson et al. 1996). Polar cap and outer gap model curves were calculated according to Rudak & Dyks (1999) and taken from Hirotani et al. (2003).

Open with DEXTER

PSR B1509-58. The rather strong unpulsed VHE $\gamma$ -ray emission of the PWN surrounding PSR B1509-58 decreases the flux sensitivity for pulsed emission in comparison with other pulsars. The limits shown in Fig. 9 represent the first sensitive search for sub-TeV $\gamma$ -rays from this pulsar. The IC component of an outer gap model prediction (Hirotani & Shibata 2001b, dotted curve, parameters j₁=0.6, j₂=0) can be ruled out by the HESS limits, although the curvature radiation component seems to be already marginally constrained by the EGRET limits at energies <200 MeV.

$\begin{figure} \par\includegraphics[width=8.4cm,clip]{6737f9.eps}\end{figure}$	Figure 9: HESS energy flux limits (99% c.l.) for pulsed emission of PSR B1509-58 (see Fig. 6 for point descriptions). The open triangles resemble the EGRET upper limits for pulsed emission (Thompson et al. 1999). The outer gap model curve (dotted) originates from Hirotani & Shibata (2001b), Fig. 6.
Open with DEXTER

$\begin{figure} \par\includegraphics[width=8.5cm,clip]{6737f10a.eps}\hspace{2mm}... ...c.eps}\hspace{2mm} \includegraphics[width=8.5cm,clip]{6737f10d.eps}\end{figure}$	Figure 10: Differential upper limits on the energy flux of pulsed VHE $\gamma$ -ray emission for non- $\gamma$ -ray pulsars as measured by HESS (at 99% confidence, point descriptions as in Fig. 6).
Open with DEXTER

$\begin{figure} \par\includegraphics[width=8.2cm,clip]{6737f11a.eps}\par\includeg... ...p]{6737f11b.eps}\par\includegraphics[width=8.2cm,clip]{6737f11c.eps}\end{figure}$	Figure 11: Differential upper limits at 99% confidence on the energy flux of unpulsed VHE $\gamma$ -ray emission for PSR J1747-2958, PSR J1801-2451 and PSR J1811-1925 ( from top to bottom) as measured by HESS.
Open with DEXTER

The remaining pulsars have not been detected in the MeV-GeV range and no detailed calculations based on any high-energy emission models are available for the TeV range.

Independent of any $\gamma$ -ray emission scenario, for all selected pulsars the flux limits obtained imply an upper limit on the ratio of the total pulsed VHE $\gamma$ -ray luminosity compared to the pulsars' spin-down power in the range 10^-6 to 10^-4. This limit on the $\gamma$ -ray efficiency is significantly lower than the corresponding EGRET sensitivity in the GeV range.

7 Conclusions

HESS data on young and energetic pulsars, including four of the eight known GeV $\gamma$ -ray pulsars, have been searched for evidence of pulsed $\gamma$ -ray emission at energies between 100 GeV and 100 TeV. No conclusive evidence for pulsed emission has been found and differential upper limits on the pulsed flux were derived using a variety of different methods, constraining the pulsed flux for a wide range of possible pulse shapes and spectra in the VHE $\gamma$ -ray range.

In the case of the well-studied $\gamma$ -ray pulsars Crab, Vela, and PSR B1509-58, the pulsed flux upper limits reported here significantly constrain the inverse Compton component of selected outer gap models for flux predictions in the TeV range which were available at the time of writing. However, the results cannot rule out the outer gap scenario in general, since not all of the different calculations published give predictions for the IC component. Also, free model parameters such as seed photon density, inclination, gap geometry or viewing angles can be adjusted for the model cases considered here. However, it remains an open challenge for theory to provide reasonable parameters consistent with the flux limits reported here.

For the complete sample of pulsars, containing a large fraction of the most energetic pulsars currently known, the absence of pulsed VHE $\gamma$ -ray emission already disfavours a significant contribution of the IC component to the energy loss mechanism of these pulsars. This sample exhibits great variety in different pulsar characteristics, providing valuable input to further studies of the outer gap IC component. For example, in the case of the binary pulsar PSR B1259-63, a large spin-down luminosity is accompanied by a rather weak magnetic field, i.e. reduced $\gamma$ -ray absorption by pair creation, and additionally the photon field originating from the companion star provides a high density target with known energy spectrum for a possible pulsed IC component, fixing an otherwise uncertain model parameter.

Although in several cases there is spatial coincidence with extended TeV $\gamma$ -ray emission, pulsed emission is not detected in VHE $\gamma$ -rays. In particular, the flux upper limits derived are of the order of 10^-4 to 10^-6 of the pulsar spin-down flux, underlining the non-magnetospheric origin of the TeV radiation component and supporting the widely accepted scenario of an effective energy transport mechanism to, and strong particle acceleration in, the pulsar wind nebula.

Finally, the upper limits given here imply a steep turnover of the pulsed high-energy spectrum at energies of a few tens of GeV for all pulsars studied which have established high-energy (>keV) emission, confirming a prediction common to all pulsar models. Moreover, for the remaining pulsars, a turnover at sub-TeV energies would be likely even if pulsed GeV $\gamma$ -ray emission with spectral properties similar to that of the established $\gamma$ -ray pulsars is detected with future sensitive instruments (Thompson 2004). However, the different scenarios for pulsar emission differ significantly in their predictions of the exact shape and energy of the turnover. Thus, the search for pulsed $\gamma$ -ray emission from pulsars provides interesting prospects for future satellite-based and low-threshold ground-based $\gamma$ -ray instruments.

Search for pulsed VHE gamma-ray emission from young pulsars with HESS

2 Pulsar candidate selection

3 VHE gamma-ray detection and data analysis

3.2 Low energy optimisation

3.3 Point source analysis

4 Periodicity analysis

4.2 Tests for pulsed emission

4.3 Flux limit determination

5 Results

6 Discussion