COMPTEL detection of pulsed gamma-ray emission from PSR B1509-58 up to at least 10 MeV

We report on the first firm detection of pulsed gamma-ray emission from PSR B1509-58 in the 0.75-30 MeV energy range in CGRO COMPTEL data collected ov er more than 6 years. The modulation significance in the 0.75-30 MeV pulse- phase distribution is 5.40 and the lightcurve is similar to the lightcurves found earlier between 0.7 and 700 keV: a single broad asymmetric pulse rea ching its maximum 0.38 +/- 0.03 in phase after the radio peak, compared to the offset of 0.30 found in the CGRO BATSE soft gamma-ray data, and 0.27 +/ - 0.01 For RXTE (2-16 keV), compatible with ASCA (0.7-2,2 keV). Analysis in narrower energy windows shows that the single broad pulse is si gnificantly detected up to similar to 10 MeV. Above 10 MeV we do detect mar ginally significant (2.1 sigma) modulation with an indication for the broad pulse. However, imaging analysis shows the presence of a strong 5.6 sigma source at the position of the pulsar. To investigate this further, we have also analysed contemporaneous CGRO EGRET data (> 30MeV)collected over a nea rly 4 year period. In the 30-100 MeV energy window, adjacent to the COMPTEL 10-30MeV range, a 4.4 sigma source can be attributed to PSR B1509-58. Timi ng analysis in this energy window yields an insignificant signal of 1.1 sig ma, but with a shape somewhat similar to that of the COMPTEL 10-30 MeV ligh tcurve. Combining the two pulse-phase distributions results in a suggestive double-peaked pulsed signal above the background level estimated in the sp atial analyses, with one broad peak near phase 0.38 (aligned with the pulse observed at lower energies) and a second narrower peak near phase 0.85, wh ich is absent for energies below 10 MeV. The modulation significance is, ho wever, only 2.3 sigma and needs confirmation. Spectral analysis based on the excess counts in the broad pulse of the ligh tcurve shows that extrapolation of the OSSE power-law spectral fit with ind ex -1.68 describes our data well up to 10MeV. Above 10MeV the spectrum brea ks abruptly. The precise location of the break/bend between 10 and 30 MeV d epends on the interpretation of the structure in the lightcurve measured by COMPTEL and EGRET above 10 MeV. Such a break in the spectrum of PSR B1509-58 has recently been interpreted in the framework of polar cap models for the explanation of gamma-ray pulsa rs, as a signature of the exotic photon splitting process in the strong mag netic field of PSR B1509-58. For that interpretation our new spectrum const rains the co-latitude to similar to 2 degrees, close to the "classical" rad ius of the polar cap. In the case of an outer-gap scenario, our spectrum re quires a dominant synchrotron component.