A deconvolution-based algorithm for crowded field photometry with unknown Point Spread Function

A new method is presented for determining the Point Spread Function (PSF) of images that lack bright and isolated stars. It is based on the same principles as the MCS (Magain, Courbin, Sohy, 1998) image deconvolution algorithm. It uses the information contained in all stellar images to achieve the double task of reconstructing the PSFs for single or multiple exposures of the same field and to extract the photometry of all point sources in the field of view. The use of the full information available allows to construct an accurate PSF. The possibility to simultaneously consider several exposures makes it very well suited to the measurement of the light curves of blended point sources from data that would be very difficult or even impossible to analyse with traditional PSF fitting techniques. The potential of the method for the analysis of ground-based and space-based data is tested on artificial images and illustrated by several examples, including HST/NICMOS images of a lensed quasar and VLT/ISAAC images of a faint blended Mira star in the halo of the giant elliptical galaxy NGC5128 (Cen A).Comment: Institutes: (1) Institut d'Astrophysique et de Geophysique, Universite de Liege, allee du 6 Aout 17, B-4000 Liege, Belgium; (2) Ecole Polytechnique Federale de Lausanne (EPFL), Laboratoire d'Astrophysique, Observatoire, CH-1290 Sauverny, Switzerland; (3) Observatoire de Geneve, 51 Chemin des Maillettes, CH-1290 Sauverny, Switzerland. 8 pages, 8 figures. Accepted for publication in A&

A method for spatial deconvolution of spectra

A method for spatial deconvolution of spectra is presented. It follows the same fundamental principles as the ``MCS image deconvolution algorithm'' (Magain, Courbin, Sohy, 1998) and uses information contained in the spectrum of a reference Point Spread Function (PSF) to spatially deconvolve spectra of very blended sources. An improved resolution rather than an infinite one is aimed at, overcoming the well known problem of ``deconvolution artefacts''. As in the MCS algorithm, the data are decomposed into a sum of analytical point sources and a numerically deconvolved background, so that the spectrum of extended sources in the immediate vicinity of bright point sources may be accurately extracted and sharpened. The algorithm has been tested on simulated data including seeing variation as a function of wavelength and atmospheric refraction. It is shown that the spectra of severely blended point sources can be resolved while fully preserving the spectrophotometric properties of the data. Extended objects ``hidden'' by bright point sources (up to 4-5 magnitudes brighter) can be accurately recovered as well, provided the data have a sufficiently high total signal-to-noise ratio (200-300 per spectral resolution element). Such spectra are relatively easy to obtain, even down to faint magnitudes, within a few hours of integration time with 10m class telescopes.Comment: 18 pages, 6 postscript figures, in press in Ap

A weak lensing comparability study of galaxy mergers that host AGNs

We compared the total mass density profiles of three different types of galaxies using weak gravitational lensing: (i) 29 galaxies that host quasars at z~0.32 that are in a post-starburst (PSQ) phase with high star formation indicating recent merger activity, (ii) 22 large elliptical galaxies from the SLACS sample that do not host a quasar at z~0.23, and (iii) 17 galaxies that host moderately luminous quasars at z~0.36 powered by disk instabilities, but with no intense star formation. On an initial test we found no evidence for a connection between the merger state of a galaxy and the profile of the halo, with the PSQ profile comparable to that of the other two samples and consistent with the Leauthaud et al. (2014) study of moderately luminous quasars in COSMOS. Given the compatibility of the two quasar samples, we combined these and found no evidence for any connection between black hole activity and the dark matter halo. All three mass profiles remained compatible with isothermality given the present data.Comment: 6 pages, 3 figures, 1 table, ACCEPTED MNRA

A deep, narrow J-band search for proto-galactic Lyman Alpha emission at redshifts z~9

We present a deep, narrow J-band search for proto--galactic Lyman Alpha emission at redshifts z~9. We combine an exceptionally deep image of the Hubble Deep Field South, obtained using a narrow band filter centred on the wavelength 1.187 microns using the VLT/ISAAC facility, with existing deep, broad band images covering optical to near infrared wavelengths. We search for z~9 Lyman Alpha emitting galaxies displaying a significant narrow band excess relative to the Js-band that are undetected at optical wavelengths. We detect no sources consistent with this criterion to the 90% point source flux limit of the NB image, F_{NB} = 3.28 x 10^{-18} ergs/s/cm2. The survey selection function indicates that we have sampled a volume of approximately 340 h^{-3} Mpc^3 to a Lyman Alpha emission luminosity of 10^{43} h^{-2} ergs/s. We conclude by considering the potential implications for the physics of the high--redshift universe.Comment: 11 pages, accepted MNRAS. Please note that the automatic postscript generation may result in a corrupted postscript file. Please use ftp://astroftp.phys.uvic.ca/pub/jwillis/jwillis_zen1.ps.gz to obtain an uncorrupted versio

Confirmation of two extended objects along the line of sight to PKS1830-211 with ESO-VLT adaptive optics imaging

We report on new high-resolution near-infrared images of the gravitationally lensed radio source PKS1830-211, a quasar at z=2.507. These adaptive optics observations, taken with the Very Large Telescope (VLT), are further improved through image deconvolution. They confirm the presence of a second object along the line of sight to the quasar, in addition to the previously known spiral galaxy. This additional object is clearly extended in our images. However, its faint luminosity does not allow to infer any photometric redshift. If this galaxy is located in the foreground of PKS1830-211, it complicates the modeling of this system and decreases the interest in using PKS1830-211 as a means to determine H0 via the time delay between the two lensed images of the quasar.Comment: Accepted in A&A Letter

Interpolating point spread function anisotropy

Planned wide-field weak lensing surveys are expected to reduce the statistical errors on the shear field to unprecedented levels. In contrast, systematic errors like those induced by the convolution with the point spread function (PSF) will not benefit from that scaling effect and will require very accurate modeling and correction. While numerous methods have been devised to carry out the PSF correction itself, modeling of the PSF shape and its spatial variations across the instrument field of view has, so far, attracted much less attention. This step is nevertheless crucial because the PSF is only known at star positions while the correction has to be performed at any position on the sky. A reliable interpolation scheme is therefore mandatory and a popular approach has been to use low-order bivariate polynomials. In the present paper, we evaluate four other classical spatial interpolation methods based on splines (B-splines), inverse distance weighting (IDW), radial basis functions (RBF) and ordinary Kriging (OK). These methods are tested on the Star-challenge part of the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) simulated data and are compared with the classical polynomial fitting (Polyfit). We also test all our interpolation methods independently of the way the PSF is modeled, by interpolating the GREAT10 star fields themselves (i.e., the PSF parameters are known exactly at star positions). We find in that case RBF to be the clear winner, closely followed by the other local methods, IDW and OK. The global methods, Polyfit and B-splines, are largely behind, especially in fields with (ground-based) turbulent PSFs. In fields with non-turbulent PSFs, all interpolators reach a variance on PSF systematics $\sigma_{sys}^2$ better than the $1\times10^{-7}$ upper bound expected by future space-based surveys, with the local interpolators performing better than the global ones

Slit and integral-field optical spectroscopy of the enigmatic quasar HE0450-2958

Interest in the quasar HE0450-2958 arose following the publication of the non-detection of its expected massive host, leading to various interpretations. This article investigates the gaseous and stellar contents of the system through additional VLT/FORS slit spectra and integral field spectroscopy from VLT/VIMOS. We apply our MCS deconvolution algorithm on slit spectra for the separation of the QSO and diffuse components, and develop a new method to remove the point sources in Integral Field Spectra, allowing extraction of velocity maps, narrow-line images, spatially resolved spectra or ionization diagrams of the surroundings of HE0450-2958. The whole system is embedded in gas, mostly ionized by the QSO radiation field and shocks associated with radio jets. The observed gas and star dynamics are unrelated, revealing a strongly perturbed system. Despite longer spectroscopic observations, the host galaxy remains undetected.Comment: 9 pages, 13 figures, in press in A&

Evaluating the effect of stellar multiplicity on the PSF of space-based weak lensing surveys

The next generation of space-based telescopes used for weak lensing surveys will require exquisite point spread function (PSF) determination. Previously negligible effects may become important in the reconstruction of the PSF, in part because of the improved spatial resolution. In this paper, we show that unresolved multiple star systems can affect the ellipticity and size of the PSF and that this effect is not cancelled even when using many stars in the reconstruction process. We estimate the error in the reconstruction of the PSF due to the binaries in the star sample both analytically and with image simulations for different PSFs and stellar populations. The simulations support our analytical finding that the error on the size of the PSF is a function of the multiple stars distribution and of the intrinsic value of the size of the PSF, i.e. if all stars were single. Similarly, the modification of each of the complex ellipticity components (e1,e2) depends on the distribution of multiple stars and on the intrinsic complex ellipticity. Using image simulations, we also show that the predicted error in the PSF shape is a theoretical limit that can be reached only if large number of stars (up to thousands) are used together to build the PSF at any desired spatial position. For a lower number of stars, the PSF reconstruction is worse. Finally, we compute the effect of binarity for different stellar magnitudes and show that bright stars alter the PSF size and ellipticity more than faint stars. This may affect the design of PSF calibration strategies and the choice of the related calibration fields.Comment: 10 pages, 6 figures, accepted in A&

Stellar classification from single-band imaging using machine learning

Information on the spectral types of stars is of great interest in view of the exploitation of space-based imaging surveys. In this article, we investigate the classification of stars into spectral types using only the shape of their diffraction pattern in a single broad-band image. We propose a supervised machine learning approach to this endeavour, based on principal component analysis (PCA) for dimensionality reduction, followed by artificial neural networks (ANNs) estimating the spectral type. Our analysis is performed with image simulations mimicking the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) in the F606W and F814W bands, as well as the Euclid VIS imager. We first demonstrate this classification in a simple context, assuming perfect knowledge of the point spread function (PSF) model and the possibility of accurately generating mock training data for the machine learning. We then analyse its performance in a fully data-driven situation, in which the training would be performed with a limited subset of bright stars from a survey, and an unknown PSF with spatial variations across the detector. We use simulations of main-sequence stars with flat distributions in spectral type and in signal-to-noise ratio, and classify these stars into 13 spectral subclasses, from O5 to M5. Under these conditions, the algorithm achieves a high success rate both for Euclid and HST images, with typical errors of half a spectral class. Although more detailed simulations would be needed to assess the performance of the algorithm on a specific survey, this shows that stellar classification from single-band images is well possible.Comment: 10 pages, 9 figures, 2 tables, accepted in A&