skip to main content

Title: Photometry on Structured Backgrounds: Local Pixel-wise Infilling by Regression
Abstract

Photometric pipelines struggle to estimate both the flux and flux uncertainty for stars in the presence of structured backgrounds such as filaments or clouds. However, it is exactly stars in these complex regions that are critical to understanding star formation and the structure of the interstellar medium. We develop a method, similar to Gaussian process regression, which we term local pixel-wise infilling (LPI). Using a local covariance estimate, we predict the background behind each star and the uncertainty of that prediction in order to improve estimates of flux and flux uncertainty. We show the validity of our model on synthetic data and real dust fields. We further demonstrate that the method is stable even in the crowded field limit. While we focus on optical-IR photometry, this method is not restricted to those wavelengths. We apply this technique to the 34 billion detections in the second data release of the Dark Energy Camera Plane Survey. In addition to removing many >3σoutliers and improving uncertainty estimates by a factor of ∼2–3 on nebulous fields, we also show that our method is well behaved on uncrowded fields. The entirely post-processing nature of our implementation of LPI photometry allows it to easily improve more » the flux and flux uncertainty estimates of past as well as future surveys.

« less
Authors:
;
Award ID(s):
2019786
Publication Date:
NSF-PAR ID:
10369014
Journal Name:
The Astrophysical Journal
Volume:
933
Issue:
2
Page Range or eLocation-ID:
Article No. 155
ISSN:
0004-637X
Publisher:
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Calcium-rich (Ca-rich) transients are a new class of supernovae (SNe) that are known for their comparatively rapid evolution, modest peak luminosities, and strong nebular calcium emission lines. Currently, the progenitor systems of Ca-rich transients remain unknown. Although they exhibit spectroscopic properties not unlike core-collapse Type Ib/c SNe, nearly half are found in the outskirts of their host galaxies, which are predominantly elliptical, suggesting a closer connection to the older stellar populations of SNe Ia. In this paper, we present a compilation of publicly available multiwavelength observations of all known and/or suspected host galaxies of Ca-rich transients ranging from far-UV to IR, and use these data to characterize their stellar populations withprospector. We estimate several galaxy parameters including integrated star formation rate, stellar mass, metallicity, and age. For nine host galaxies, the observations are sensitive enough to obtain nonparametric star formation histories, from which we recover SN rates and estimate probabilities that the Ca-rich transients in each of these host galaxies originated from a core-collapse versus Type Ia-like explosion. Our work supports the notion that the population of Ca-rich transients do not come exclusively from core-collapse explosions, and must either be only from white dwarf stars or a mixed populationmore »of white dwarf stars with other channels, potentially including massive star explosions. Additional photometry and explosion site spectroscopy of larger samples of Ca-rich host galaxies will improve these estimates and better constrain the ratio of white dwarf versus massive star progenitors of Ca-rich transients.

    « less
  2. Abstract

    We explore the fascinating eclipses and dynamics of the compact hierarchical triple-star system KOI-126 (KIC 5897826). This system is composed of a pair of M-dwarf stars (KOI-126 B and C) in a 1.74 day orbit that revolve around an F star (KOI-126 A) every 34 days. Complex eclipse shapes are created as the M stars transit the F star, due to two effects: (1) the duration of the eclipse is a significant fraction of the M-star orbital period, so the prograde or retrograde motion of the M stars in their orbit lead to unusually short or long duration eclipses; (2) due to 3-body dynamics, the M-star orbit precesses with an astonishingly quick timescale of 1.74 yr for the periastron (apsidal) precession, and 2.73 yr for the inclination and nodal angle precession. Using the full Kepler data set, supplemented with ground-based photometry, plus 29 radial velocity measurements that span 6 yr, our photodynamical modeling yields masses ofMA= 1.2713 ± 0.0047M(0.37%),MB= 0.23529 ± 0.00062M(0.26%), andMC= 0.20739 ± 0.00055M(0.27%) and radii ofRA= 1.9984 ± 0.0027R(0.14%),RB= 0.25504 ± 0.00076R(0.3%), andRC= 0.23196 ± 0.00069R(0.3%). We also estimate the apsidal motion constant of the M dwarfs, a parameter that characterizes the internal mass distribution. Althoughmore »it is not particularly precise, we measure a mean apsidal motion constant,k2¯, of0.0460.028+0.046, which is approximately 2σlower than the theoretical model prediction of 0.150. We explore possible causes for this discrepancy.

    « less
  3. ABSTRACT

    We present SCUBA-2 850 $\mathrm{ \mu}$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850 $\mathrm{ \mu}$m sources in 9 of 13 of these candidate protoclusters show significant overdensities compared to the field, with the probability <10−2 assuming the sources are randomly distributed in the sky. Using the 250, 350, 500, and 850 $\mathrm{ \mu}$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution templates with an MCMC method. The photometric redshift distribution, peaking at 2 < z < 3, is consistent with that of known z > 2 protoclusters and the peak of the cosmic star formation rate density (SFRD). We find that the 850 $\mathrm{ \mu}$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim 10^{12}\, \mathrm{L}_{\odot }$ and star formation rates of SFR  = (500–1500) M⊙ yr−1. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorized into four groups: six of them being high-redshift starbursting protoclusters, one being a lower redshift cluster or protocluster, three being protoclusters that contain lensed dusty star-forming galaxies or are rich in 850 $\mathrm{ \mu}$m sources, and three regionsmore »without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850 $\mathrm{ \mu}$m sources.

    « less
  4. Abstract

    We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the Hubble Space Telescope–based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300–9400 Å as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all presubmission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures ≳5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12 mag range. Overall, the mean of our system is calibrated to themore »mean of CALSPEC at the level of ∼3 mmag. With our large number of observations, careful cross-checks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys.

    « less
  5. Abstract

    Massive, star-forming clumps are a common feature of high-redshift star-forming galaxies. How they formed, and why they are so rare at low redshift, remains unclear. In this paper we identify the largest sample yet of clumpy galaxies (7050) at low redshift using data from the citizen science project Galaxy Zoo: Clump Scout, in which volunteers classified 58,550 Sloan Digital Sky Survey (SDSS) galaxies spanning redshift 0.02 <z< 0.15. We apply a robust completeness correction by comparing with simulated clumps identified by the same method. Requiring that the ratio of clump to galaxy flux in the SDSSuband be greater than 8% (similar to clump definitions used by other works), we estimate the fraction of local star-forming galaxies hosting at least one clump (fclumpy) to be3.220.34+0.38%. We also compute the same fraction with a less stringent relative flux cut of 3% (12.680.88+1.38%), as the higher number count and lower statistical noise of this fraction permit finer comparison with future low-redshift clumpy galaxy studies. Our results reveal a sharp decline infclumpyover 0 <z< 0.5. The minor merger rate remains roughly constant over the same span, so we suggest that minor mergers are unlikely to be themore »primary driver of clump formation. Instead, the rate of galaxy turbulence is a better tracer forfclumpyover 0 <z< 1.5 for galaxies of all masses, which supports the idea that clump formation is primarily driven by violent disk instability for all galaxy populations during this period.

    « less