Bayesian multiscale deconvolution applied to gamma-ray spectroscopy

A common task in gamma-ray astronomy is to extract spectral information, such as model constraints and incident photon spectrum estimates, given the measured energy deposited in a detector and the detector response. This is the classic problem of spectral “deconvolution” or spectral inversion. The methods of forward folding (i.e., parameter fitting) and maximum entropy “deconvolution” (i.e., estimating independent input photon rates for each individual energy bin) have been used successfully for gamma-ray solar flares (e.g., Rank, 1997; Share and Murphy, 1995). These methods have worked well under certain conditions but there are situations were they don’t apply. These are: 1) when no reasonable model (e.g., fewer parameters than data bins) is yet known, for forward folding; 2) when one expects a mixture of broad and narrow features (e.g., solar flares), for the maximum entropy method; and 3) low count rates and low signal-to-noise, for both. Low count rates are a problem because these methods (as they have been implemented) assume Gaussian statistics but Poisson are applicable. Background subtraction techniques often lead to negative count rates. For Poisson data the Maximum Likelihood Estimator (MLE) with a Poisson likelihood is appropriate. Without a regularization term, trying to estimate the “true” individual input photon rates per bin can be an ill-posed problem, even without including both broad and narrow features in the spectrum (i.e., amultiscale approach). One way to implement this regularization is through the use of a suitable Bayesian prior. Nowak and Kolaczyk (1999) have developed a fast, robust, technique using a Bayesian multiscale framework that addresses these problems with added algorithmic advantages. We outline this new approach and demonstrate its use with time resolved solar flare gamma-ray spectroscopy

Extended γ‐ray emission in solar flares

During the solar flare events on 11 and 15 June 1991, COMPTEL measured extended emission in the neutron capture line for about 5 hours after the impulsive phase. The time profiles can be described by a double exponential decay with decay constants on the order of 10 min for the fast and 200 min for the slow component. Within the statistical uncertainty both flares show the same long‐term behaviour. The spectrum during the extended phase is significantly harder than during the impulsive phase and pions are not produced in significant numbers before the beginning of the extended emission. Our results with the measurements of others allow us to rule out long‐term trapping of particles in non‐turbulent loops to explain the extended emission of these two flares and our data favour models based on continued acceleration

COMPTEL upper limits for Seyfert galaxies

The gamma‐ray emission of Seyfert galaxies has fallen far short of pre‐GRO expectations. No single object of this class has been detected by either COMPTEL or EGRET, and OSSE has detected only a fraction of the Seyferts expected. To derive a more stringent upper limit to the emission from these objects in the energy ranges 0.75 to 1 and 1 to 3 MeV, we have summed a large number of COMPTEL observations acquired during Phase 1 of the GRO mission. From a total of 47 observations of 23 individual X‐ray selected Seyfert galaxies, we derive preliminary upper limits of 8×10−8 photons/(cm2 s keV) in the 0.75‐1 MeV band and 1×10−8 photons/(cm2 s keV) in the 1‐3 MeV band

Energetic proton spectra in the 11 June 1991 solar flare

We have studied a subset of the 11 June 1991 solar flare γ-ray data that we believe arise from soft proton or ion spectra. Using data from the COMPTEL instrument on the Compton Observatory we discuss the gamma-ray intensities at 2.223 MeV, 4–7 MeV, and 8–30 MeV in terms of the parent proton spectrum responsible for the emission

COMPTEL upper limits for the 56Co γ-rays from SN1998bu

The type Ia supernova SN 1998bu in M96 was observed by COMPTEL for a total of 88 days starting 17 days after the detection of the SN. A special mode improving the low-energy sensitivity was invoked. We obtained images in the 847 keV and 1238 keV lines of 56Co using an improved point-spread function for the low-energies. We do not detect SN1998bu. Sensitive upper limits at both energies constrain the standard supernova model for this event

Activation in the COMPTEL double-scattering gamma-ray telescope

Abstract-The COMPTEL gamma-ray telescope has been operating in low Earth orbit for six years, since the launch of the Compton Gamma-Ray Observatory in April 1991. Comparisons of data for different orbits and epochs show evidence of activation on time scales from minutes (27Mg, q,2=9.5 min) to years C2Na, q&.58 yr). The activation is correlated with both the orbital altitude and solar cosmic-ray modulation. Because it requires coincident measurements in two different detectors, COMPTEL is most susceptible to instrumental background events in which two or more photons are produced simultaneously

Energetic proton spectra in the 11 June 1991 solar flare

The June 11, 1991 gamma-ray flare seen by the Compton Gamma-ray Observatory (CGRO) displays several features that make it a dynamic and rich event. It is a member of a class of long duration gamma-ray events with both 2.223 MeV and greater than 8 MeV emission for hours after the impulsive phase. It also contains an inter-phase between the impulsive and extended phases that presents a challenge to the standard gamma-ray line (GRL) flare picture. This phase has strong 2.223 MeV emission and relatively weak 4.44 MeV emission indicative of a very hard parent proton spectrum. However, this would indicate emission greater than 8 MeV, which is absent from this period. We present the application of new spectroscopy techniques to this phase of the flare in order to present a reasonable explanation for this seemly inconsistent picture

Gamma ray measurements of the 1991 November 15 solar flare

The 1991 November 15 X1.5 flare was a well observed solar event. Comprehensive data from ground-based observatories and spacecraft provide the basis for a contextual interpretation of gamma-ray spectra from the Compton Gamma Ray Observatory (CGRO). In particular, spectral, spatial, and temporal data at several energies are necessary to understand the particle dynamics and the acceleration mechanism(s) within this flare. X-ray images, radio, Ca XIX data and magnetograms provide morphological information on the acceleration region [4,5], while gamma-ray spectral data provide information on the parent ion spectrum. Furthermore, time profiles in hard X-rays and gamma-rays provide valuable information on temporal characteristics of the energetic particles. We report the results of our analysis of the evolution of this flare as a function of energy (∼25 keV–2.5 MeV) and time. These results, together with other high energy data (e.g. from experiments on Yohkoh, Ulysses, and PVO) may assist in identifying and understanding the acceleration mechanism(s) taking place in this event