Imaging the Gamma-Ray Sky with SPI aboard INTEGRAL

The spectrometer SPI on INTEGRAL allows for the first time simultaneous imaging of diffuse and point-like emission in the hard X-ray and soft gamma-ray regime. To fully exploit the capabilities of the instrument, we implemented the MREM image deconvolution algorithm, initially developed for COMPTEL data analysis, to SPI data analysis. We present the performances of the algorithm by means of simulations and apply it to data accumulated during the first 2 mission years of INTEGRAL. Skymaps are presented for the 1809 keV gamma-ray line, attributed to the radioactive decay of 26Al, and for continuum energy bands, covering the range 20 keV - 3 MeV. The 1809 keV map indicates that emission is clearly detected by SPI from the inner Galactic radian and from the Cygnus region. The continuum maps reveal the transition between a point-source dominated hard X-ray sky to a diffuse emission dominated soft gamma-ray sky. From the skymaps, we extract a Galactic ridge emission spectrum that matches well SPI results obtained by model fitting. By comparing our spectrum with the cumulative flux measured by IBIS from point sources, we find indications for the existence of an unresolved or diffuse emission component above ~100 keV.Comment: 12 pages, 7 figure

Imfit: A Fast, Flexible New Program for Astronomical Image Fitting

I describe a new, open-source astronomical image-fitting program called Imfit, specialized for galaxies but potentially useful for other sources, which is fast, flexible, and highly extensible. A key characteristic of the program is an object-oriented design which allows new types of image components (2D surface-brightness functions) to be easily written and added to the program. Image functions provided with Imfit include the usual suspects for galaxy decompositions (Sersic, exponential, Gaussian), along with Core-Sersic and broken-exponential profiles, elliptical rings, and three components which perform line-of-sight integration through 3D luminosity-density models of disks and rings seen at arbitrary inclinations. Available minimization algorithms include Levenberg-Marquardt, Nelder-Mead simplex, and Differential Evolution, allowing trade-offs between speed and decreased sensitivity to local minima in the fit landscape. Minimization can be done using the standard chi^2 statistic (using either data or model values to estimate per-pixel Gaussian errors, or else user-supplied error images) or Poisson-based maximum-likelihood statistics; the latter approach is particularly appropriate for cases of Poisson data in the low-count regime. I show that fitting low-S/N galaxy images using chi^2 minimization and individual-pixel Gaussian uncertainties can lead to significant biases in fitted parameter values, which are avoided if a Poisson-based statistic is used; this is true even when Gaussian read noise is present.Comment: pdflatex, 27 pages, 19 figures. Revised version, accepted by ApJ. Programs, source code, and documentation available at: http://www.mpe.mpg.de/~erwin/code/imfit

Photometry of supernovae in an image series : methods and application to the Supernova Legacy Survey (SNLS)

We present a technique to measure lightcurves of time-variable point sources on a spatially structured background from imaging data. The technique was developed to measure light curves of SNLS supernovae in order to infer their distances. This photometry technique performs simultaneous PSF photometry at the same sky position on an image series. We describe two implementations of the method: one that resamples images before measuring fluxes, and one which does not. In both instances, we sketch the key algorithms involved and present the validation using semi-artificial sources introduced in real images in order to assess the accuracy of the supernova flux measurements relative to that of surrounding stars. We describe the methods required to anchor these PSF fluxes to calibrated aperture catalogs, in order to derive SN magnitudes. We find a marginally significant bias of 2 mmag of the after-resampling method, and no bias at the mmag accuracy for the non-resampling method. Given surrounding star magnitudes, we determine the systematic uncertainty of SN magnitudes to be less than 1.5 mmag, which represents about one third of the current photometric calibration uncertainty affecting SN measurements. The SN photometry delivers several by-products: bright star PSF flux mea- surements which have a repeatability of about 0.6%, as for aperture measurements; we measure relative astrometric positions with a noise floor of 2.4 mas for a single-image bright star measurement; we show that in all bands of the MegaCam instrument, stars exhibit a profile linearly broadening with flux by about 0.5% over the whole brightness range.Comment: Accepted for publication in A&A. 20 page

Foreground removal from CMB temperature maps using an MLP neural network

One of the main obstacles in extracting the Cosmic Microwave Background (CMB) signal from observations in the mm-submm range is the foreground contamination by emission from galactic components: mainly synchrotron, free-free and thermal dust emission. Due to the statistical nature of the intrinsic CMB signal it is essential to minimize the systematic errors in the CMB temperature determinations. Following the available knowledge of the spectral behavior of the galactic foregrounds simple, power law-like spectra have been assumed. The feasibility of using a simple neural network for extracting the CMB temperature signal from the combined CMB and foreground signals has been investigated. As a specific example, we have analysed simulated data, like that expected from the ESA Planck Surveyor mission. A simple multilayer perceptron neural network with 2 hidden layers can provide temperature estimates, over more than 80 percent of the sky, that are to a high degree uncorrelated with the foreground signals. A single network will be able to cover the dynamic range of the Planck noise level over the entire sky.Comment: Accepted for publication in Astrophysics and Space Scienc

The Sloan Digital Sky Survey Quasar Lens Search. I. Candidate Selection Algorithm

We present an algorithm for selecting an uniform sample of gravitationally lensed quasar candidates from low-redshift (0.6<z<2.2) quasars brighter than i=19.1 that have been spectroscopically identified in the SDSS. Our algorithm uses morphological and color selections that are intended to identify small- and large-separation lenses, respectively. Our selection algorithm only relies on parameters that the SDSS standard image processing pipeline generates, allowing easy and fast selection of lens candidates. The algorithm has been tested against simulated SDSS images, which adopt distributions of field and quasar parameters taken from the real SDSS data as input. Furthermore, we take differential reddening into account. We find that our selection algorithm is almost complete down to separations of 1'' and flux ratios of 10^-0.5. The algorithm selects both double and quadruple lenses. At a separation of 2'', doubles and quads are selected with similar completeness, and above (below) 2'' the selection of quads is better (worse) than for doubles. Our morphological selection identifies a non-negligible fraction of single quasars: To remove these we fit images of candidates with a model of two point sources and reject those with unusually small image separations and/or large magnitude differences between the two point sources. We estimate the efficiency of our selection algorithm to be at least 8% at image separations smaller than 2'', comparable to that of radio surveys. The efficiency declines as the image separation increases, because of larger contamination from stars. We also present the magnification factor of lensed images as a function of the image separation, which is needed for accurate computation of magnification bias.Comment: 15 pages, 17 figures, 4 tables, accepted for publication in A

The all-sky distribution of 511 keV electron-positron annihilation emission

We present a map of 511 keV electron-positron annihilation emission, based on data accumulated with the SPI spectrometer aboard ESA's INTEGRAL gamma-ray observatory, that covers approximately 95% of the celestial sphere. 511 keV line emission is significantly detected towards the galactic bulge region and, at a very low level, from the galactic disk. The bulge emission is highly symmetric and is centred on the galactic centre with an extension of 8 deg. The emission is equally well described by models that represent the stellar bulge or halo populations. The disk morphology is only weakly constrained by the present data, being compatible with both the distribution of young and old stellar populations. The 511 keV line flux from the bulge and disk components is 1.05e-3 ph cm-2 s-1 and 0.7e-3 ph cm-2 s-1, respectively, corresponding to a bulge-to-disk flux ratio in the range 1-3. Assuming a positronium fraction of 0.93 this translates into annihilation rates of 1.5e43 s-1 and 3e42 s-1, respectively. The ratio of the bulge luminosity to that of the disk is in the range 3-9. We find no evidence for a point-like source in addition to the diffuse emission, down to a typical flux limit of 1e-4 ph cm-2 s-1. We also find no evidence for the positive latitude enhancement that has been reported from OSSE measurements; our 3 sigma upper flux limit for this feature is 1.5e-4 ph cm-2 s-1. The disk emission can be attributed to the beta+ decay of the radioactive species 26Al and 44Ti. The bulge emission arises from a different source which has only a weak or no disk component. We suggest that Type Ia supernovae and/or low-mass X-ray binaries are the prime candidates for the source of the galactic bulge positrons. Light dark matter annihilation could also explain the observed 511 keV bulge emission characteristics.Comment: accepted for publication in A&

The QUEST Data Processing Software Pipeline

A program that we call the QUEST Data Processing Software Pipeline has been written to process the large volumes of data produced by the QUEST camera on the Samuel Oschin Schmidt Telescope at the Palomar Observatory. The program carries out both aperture and PSF photometry, combines data from different repeated observations of the same portion of sky, and produces a Master Object Catalog. A rough calibration of the data is carried out. This program, as well as the calibration procedures and quality checks on the output are described.Comment: 17 pages, 1 table, 8 figure