Difference image photometry with bright variable backgrounds

• Published in: Monthly notices of the Royal Astronomical Society 2010-11, Vol.409 (1), p.247-258
• Main Authors: Kerins, E., Darnley, M. J., Duke, J. P., Gould, A., Han, C., Newsam, A., Park, B. G., Street, R.
• Format: Article
• Language: English
• Subjects: Algorithms; Astronomy; Density; Earth, ocean, space; Exact sciences and technology; galaxies: individual: M31; Imaging; Matching; Photometry; Stars & galaxies; Surface brightness; techniques: image processing; techniques: photometric; Variable stars
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1111/j.1365-2966.2010.17290.x

Over the last two decades the Andromeda galaxy (M31) has been something of a test-bed for methods aimed at obtaining accurate time-domain relative photometry within highly crowded fields. Difference imaging methods, originally pioneered towards M31, have evolved into sophisticated methods, such as the optimal image subtraction (OIS) method of Alard & Lupton, that today are most widely used to survey variable stars, transients and microlensing events in our own Galaxy. We show that modern difference image analysis (DIA) algorithms such as OIS, whilst spectacularly successful towards the Milky Way bulge, may perform badly towards high surface brightness targets such as the M31 bulge. Poor results can occur in the presence of common systematics which add spurious flux contributions to images, such as internal reflections, scattered light or fringing. Using data from the Angstrom Project microlensing survey of the M31 bulge, we show that very good results are usually obtainable by first performing careful photometric alignment prior to using OIS to perform point spread function (PSF) matching. This separation of background matching and PSF matching, a common feature of earlier M31 photometry techniques, allows us to take full advantage of the powerful PSF matching flexibility offered by OIS towards high surface brightness targets. We find that difference images produced this way have noise distributions close to Gaussian, showing significant improvement upon results achieved using OIS alone. We show that with this correction light curves of variable stars and transients can be recovered to within ∼10 arcsec of the M31 nucleus. Our method is simple to implement and is quick enough to be incorporated within real-time DIA pipelines. We also demonstrate that OIS is remarkably robust even when, as in the case of the central regions of the M31 bulge, the sky density of variable sources approaches the confusion limit.