A novel image reconstruction method applied to deep Hubble Space Telescope Images

We have developed a method for the linear reconstruction of an image from undersampled, dithered data, which has been used to create the distributed, combined Hubble Deep Field images -- the deepest optical images yet taken of the universe. The algorithm, known as Variable-Pixel Linear Reconstruction (or informally as "drizzling"), preserves photometry and resolution, can weight input images according to the statistical significance of each pixel, and removes the effects of geometric distortion both on image shape and photometry. In this paper, the algorithm and its implementation are described, and measurements of the photometric accuracy and image fidelity are presented. In addition, we describe the use of drizzling to combine dithered images in the presence of cosmic rays.Comment: Invited paper, to appear in Applications of Digital Image Processing XX, ed. A. Tescher, Proc. S.P.I.E. vol. 3164, in press; 6 pages, 4 included figures, SPIE LaTex style file include

Fast and Slow Rotators in the Densest Environments: a SWIFT IFS study of the Coma Cluster

We present integral-field spectroscopy of 27 galaxies in the Coma cluster observed with the Oxford SWIFT spectrograph, exploring the kinematic morphology-density relationship in a cluster environment richer and denser than any in the ATLAS3D survey. Our new data enables comparison of the kinematic morphology relation in three very different clusters (Virgo, Coma and Abell 1689) as well as to the field/group environment. The Coma sample was selected to match the parent luminosity and ellipticity distributions of the early-type population within a radius 15' (0.43 Mpc) of the cluster centre, and is limited to r' = 16 mag (equivalent to M_K = -21.5 mag), sampling one third of that population. From analysis of the lambda-ellipticity diagram, we find 15+-6% of early-type galaxies are slow rotators; this is identical to the fraction found in the field and the average fraction in the Virgo cluster, based on the ATLAS3D data. It is also identical to the average fraction found recently in Abell 1689 by D'Eugenio et al.. Thus it appears that the average slow rotator fraction of early type galaxies remains remarkably constant across many different environments, spanning five orders of magnitude in galaxy number density. However, within each cluster the slow rotators are generally found in regions of higher projected density, possibly as a result of mass segregation by dynamical friction. These results provide firm constraints on the mechanisms that produce early-type galaxies: they must maintain a fixed ratio between the number of fast rotators and slow rotators while also allowing the total early-type fraction to increase in clusters relative to the field. A complete survey of Coma, sampling hundreds rather than tens of galaxies, could probe a more representative volume of Coma and provide significantly stronger constraints, particularly on how the slow rotator fraction varies at larger radii.Comment: Accepted for publication in MNRA

M32+/-1

WFPC-2 images are used to study the central structure of M31, M32, and M33. The dimmer peak, P2, of the M31 double nucleus is centered on the bulge to 0.1", implying that it is the dynamical center of M31. P2 contains a compact source discovered by King et al. (1995) at 1700 A. This source is resolved, with r_{1/2} approx0.2 pc. It dominates the nucleus at 3000 A, and is consistent with late B-early A stars. This probable cluster may consist of young stars and be an older version of the cluster of hot stars at the center of the Milky Way, or it may consist of heavier stars built up from collisions in a possible cold disk of stars orbiting P2. In M32, the central cusp rises into the HST limit with gamma approx0.5, and the central density rho_0>10^7M_sol pc^-3. The V-I and U-V color profiles are flat, and there is no sign of an inner disk, dust, or any other structure. This total lack of features seems at variance with a nominal stellar collision time of 2 X 10^10 yr, which implies that a significant fraction of the light in the central pixel should come from blue stragglers. InM33, the nucleus has an extremely steep gamma=1.49 power-law profile for 0.05"<r<0.2" that becomes shallower as the HST resolution limit is approached. The profile for r<0.04" has either a gamma approx 0.8 cusp or a small core with r_c ~<0.13 pc. The central density is rho_0 > 2 10^6M_sol pc^-3, and the implied relaxation time is only ~3 X 10^6 yr, indicating that the nucleus is highly relaxed. The accompanying short collision time of 7 X 10^9 yr predicts a central blue straggler component quantitatively consistent with the strong V-I and B-R color gradients seen with HST and from the ground.Comment: 44 pages, 22 figures (7 as separate JPEG images), submitted to The Astronomical Journal. Full postscript image available at http://www.noao.edu/noao/staff/lauer/lauer_paper

Estimation from quantized Gaussian measurements: when and how to use dither

Subtractive dither is a powerful method for removing the signal dependence of quantization noise for coarsely quantized signals. However, estimation from dithered measurements often naively applies the sample mean or midrange, even when the total noise is not well described with a Gaussian or uniform distribution. We show that the generalized Gaussian distribution approximately describes subtractively dithered, quantized samples of a Gaussian signal. Furthermore, a generalized Gaussian fit leads to simple estimators based on order statistics that match the performance of more complicated maximum likelihood estimators requiring iterative solvers. The order statistics-based estimators outperform both the sample mean and midrange for nontrivial sums of Gaussian and uniform noise. Additional analysis of the generalized Gaussian approximation yields rules of thumb for determining when and how to apply dither to quantized measurements. Specifically, we find subtractive dither to be beneficial when the ratio between the Gaussian standard deviation and quantization interval length is roughly less than one-third. When that ratio is also greater than 0.822/K^0.930 for the number of measurements K > 20, estimators we present are more efficient than the midrange.https://arxiv.org/abs/1811.06856Accepted manuscrip

A near-IR variability study of the Galactic black hole: a red noise source with no detected periodicity

We present the results of near-infrared (2 and 3 microns) monitoring of Sgr A*-IR with 1 min time sampling using the natural and laser guide star adaptive optics (LGS AO) system at the Keck II telescope. Sgr A*-IR was observed continuously for up to three hours on each of seven nights, between 2005 July and 2007 August. Sgr A*-IR is detected at all times and is continuously variable, with a median observed 2 micron flux density of 0.192 mJy, corresponding to 16.3 magnitude at K'. These observations allow us to investigate Nyquist sampled periods ranging from about 2 minutes to an hour. Using Monte Carlo simulations, we find that the variability of Sgr A* in this data set is consistent with models based on correlated noise with power spectra having frequency dependent power law slopes between 2.0 to 3.0, consistent with those reported for AGN light curves. Of particular interest are periods of ~20 min, corresponding to a quasi-periodic signal claimed based upon previous near-infrared observations and interpreted as the orbit of a 'hot spot' at or near the last stable orbit of a spinning black hole. We find no significant periodicity at any time scale probed in these new observations for periodic signals. This study is sensitive to periodic signals with amplitudes greater than 20% of the maximum amplitude of the underlying red noise component for light curves with duration greater than ~2 hours at a 98% confidence limit.Comment: 37 pages, 2 tables, 17 figures, accepted by Ap