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

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 measurements of MeV gamma-ray burst spectra

We present results from the on-going spectral analysis of gamma-ray bursts measured by the COMPTEL instrument in its main Compton “Telescope” observing mode (0.75–30 MeV). Thus far, 18 bursts have been analyzed from three years (April 1991–April 1994) of observations. The time-averaged spectra of these events above 1 MeV are all consistent with a simple power law model with spectral index in the range 1.5–3.5. Exponential, thermal bremsstrahlung and thermal synchrotron models are statistically inconsistent with the burst sample, although they can adequately describe some of the individual burst spectra. We find good agreement between burst spectra measured simultaneously by BATSE, COMPTEL and EGRET, which typically show a spectral transition or “break” in the BATSE energy range around a few hundred keV followed by simple power law emission extending to hundreds of MeV. However, the temporal relation between MeV and GeV (e.g., as measured by EGRET) burst emission is still unclear. Measurement of rapid variability at MeV energies in the stronger bursts provides evidence that either the sources are nearby (within the Galaxy) or the gamma-ray emission is relativistically beamed

COMPTEL 1.8 MeV all sky survey: The Cygnus region

We present an updated version of COMPTEL’s 1.809 MeV sky survey. Based on eight years of observations we compare results from different imaging techniques using background from adjacent energy bands. We confirm the previously reported characteristics of the galactic 1.809 MeV emission, specifically an extended galactic ridge emission, mainly concentrated towards the inner galaxy, a peculiar emission feature in the Cygnus region, and a low-intensity ridge extending towards Carina and Vela. Because this gamma ray line is due to the decay of radioactive 26Al, predominantly synthesized in massive stars, one anticipates flux enhancements aligned with regions of recent star formation. This is born out by the observations. In particular the Cygnus feature, first presented in 1996 based on three years of COMPTEL data, is confirmed. Based on the stellar population we distinguish three prominent areas in this region, for which we separately derive fluxes, and discuss interpretations

COMPTEL gamma ray and neutron measurements of solar flares

COMPTEL on the Compton Gamma Ray Observatory has measured the flux of x‐rays and neutrons from several solar flares. These data have also been used to image the Sun in both forms of radiation. Unusually intense flares occurred during June 1991 yielding data sets that offer some new insight into of how energetic protons and electrons are accelerated and behave in the solar environment. We summarize here some of the essential features in the solar flare data as obtained by COMPTEL during June 1991

Neutron and gamma‐ray measurements of the solar flare of 1991 June 9

The COMPTEL Imaging Compton Telescope on‐board the Compton Gamma Ray Observatory measured significant neutron and γ‐ray fluxes from the solar flare of 9 June 1991. The γ‐ray flux had an integrated intensity (≳1 MeV) of ∼30 cm−2, extending in time from 0136 UT to 0143 UT, while the time of energetic neutron emission extended approximately 10 minutes longer, indicating either extended proton acceleration to high energies or trapping and precipitation of energetic protons. The production of neutrons without accompanying γ‐rays in the proper proportion indicates a significant hardening of the precipitating proton spectrum through either the trapping or extended acceleration process

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

Improved COMPTEL maps of the milky way

In the course of the mission we have gradually developed an analysis method that separates in an iterative manner the celestial emission and the (a priori unknown) instrumental background. It has become our standard analysis tool for point sources. We illustrate here that this method is widely applicable now. It provides mutually consistent sets of model-fitting parameters (spectra) and sky maps, both for continuum and line studies. Because of the wide applicability, it has been possible to make various cross-checks while building up confidence in this procedure

The angular distribution of COMPTEL Gamma-Ray bursts

The superior burst location capability of the COMPTEL instrument aboard the Compton Gamma-Ray Observatory allows us to study the small-scale angular distribution of burst sources with good sensitivity even though the number of burst detections is small. We accumulate four years (April 1991–April 1995) of observations to form a catalog of 27 burst locations whose mean 1σ uncertainty is ∼1°. We find that the COMPTEL bursts are consistent with an isotropic distribution of sources, yet the spatial coincidence of two of the bursts within COMPTEL’s angular resolution indicates the possibility of repetition. This possibility is studied using the two-point angular correlation function and the nearest neighbor statistic. Model dependent upper limits on the fraction of repeating sources are derived

Recent results from COMPTEL observations of Cygnus X‐1

The COMPTEL experiment on the Compton Gamma‐Ray Observatory (CGRO) has now observed Cyg X‐1 on four separate occasions during phase 1 and phase 2 of its orbital mission (April, 1991 to August, 1993). Here we report on the results of the latest analysis of these data, which provide a spectrum extending to energies greater than 2 MeV. A spectral analysis of these data, in the context of a classical Comptonization model, indicates an electron temperature much higher than previous hard X‐ray measurements would suggest (200 keV vs 60–80 keV). This implies either some limitations in the standard Comptonization model and/or the need to incorporate a reflected component in the hard X‐ray spectrum. Although significant variability near 1 MeV has been observed, there is no evidence for any ‘MeV excess.