COMPTEL observations of the blazars 3C 454.3 and CTA 102

We have analyzed the two blazars of 3C 454.3 and CTA 102 using all available COMPTEL data from 1991 to 1999. In the 10–30 MeV band, emission from the general direction of the sources is found at the 4σ-level, being consistent with contributions from both sources. Below 10 MeV only 3C 454.3 is significantly detected, with the strongest evidence (5.6 σ) in the 3–10 MeV band. Significant flux variability is not observed for both sources, while a low emission is seen most of the years in the 3–10 MeV light curve for 3C 454.3. Its time-averaged MeV spectrum suggests a power maximum between 3 to 10 MeV

COMPTEL observations of the Virgo blazars 3C 273 and 3C 279

We report the main MeV properties (detections, light curves, spectra) of the Virgo blazars 3C 273 and 3C 279 which were derived from a consistent analysis of all COMPTEL Virgo observations between 1991 and 1997

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

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

Variable precision arithmetic: a Fortran 95 module

Abstract. This paper describes the design and development of a software package supporting variable precision arithmetic as a semantic extension to the Fortran 95 language. The working precision of the arithmetic supported by this package can be dynamically and arbitrarily varied. The facility exploits the data-abstraction capabilities of Fortran 95 and allows the operations to be used elementally with array operands as well as with scalars. The number system is defined in such a way as to be closed under all of the basic operations of normal arithmetic; no programterminating numerical exceptions can occur. Precision loss situations like underflow and overflow are handled by defining special value representations that preserve as much of the numeric information as is practical and the operation semantics are defined so that these exceptional values propagate as appropriate to reflect this loss of information. The number system uses an essentially conventional variable precision floating-point representation. When operations can be performed exactly within the currently-set working precision limit, the excess trailing zero digits are not stored, nor do they take part in future operations. This is both economical in storage and improves efficiency