Knot a Bad Idea: Testing BLISS Mapping for Spitzer Space Telescope Photometry

Abstract

Much of transiting exoplanet science relies on high-precision photometry. The current generation of instruments can exhibit sensitivity variations greater than the astrophysical signals. For the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope, a popular way to handle this is BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping. As part of a Markov Chain Monte Carlo (MCMC), BLISS mapping estimates the sensitivity at many locations (knots) on the pixel, then interpolates to the target star's centroids. We show that such embedded optimization schemes can misfit or bias parameters. Thus, we construct a model of $\text{Spitzer}$ eclipse light curves to test the accuracy and precision of BLISS mapping. We compare standard BLISS mapping to a variant where the knots are fit during the MCMC, as well as to a polynomial model. Both types of BLISS mapping give similar eclipse depths, and we find that standard knots behave like real parameters. Standard BLISS mapping is therefore a reasonable shortcut to fitting for knots in an MCMC. BLISS maps become inaccurate when the photon noise is low, but typically approximate the real sensitivity well. We also find there is no perfect method for choosing the ideal number of BLISS knots to use on given data. BLISS mapping gives fits that are usually more accurate than precise (i.e. they are overly conservative), and the routine is more precise than polynomial models for significant eclipses or pixels with more varied sensitivities. BLISS mapping has better predictive power for most of these particular synthetic data, depending on how one treats time-correlated residuals. Overall, we conclude that BLISS mapping can be a reasonable sensitivity model for IRAC photometry.

arXiv comment32 pages, 13 figures; Accepted by PASP 9/26/16. Major updates: higher amplitude for main detector signals, estimating good number of BLISS knots for given data, more/better synthetic light curves, trying projected sensitivity along pixel axes, uncertainties in Figures 10-13. Most edits in Abstract, Sections 4-6. New findings but main conclusion same (i.e. BLISS Mapping can be acceptable)

References in corpus (40)

• emcee: The MCMC Hammer
• The Infrared Array Camera (IRAC) for the Spitzer Space Telescope
• The Spitzer Space Telescope Mission
• A map of the day-night contrast of the extrasolar planet HD 189733b
• Detection of Thermal Emission from an Extrasolar Planet
• Absolute Calibration of the Infrared Array Camera on the Spitzer Space Telescope
• The effect of red noise on planetary transit detection
• The 3.6-8.0 Micron Broadband Emission Spectrum of HD 209458b: Evidence for an Atmospheric Temperature Inversion
• The climate of HD 189733b from fourteen transits and eclipses measured by Spitzer
• 3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b
• A Gaussian process framework for modelling instrumental systematics: application to transmission spectroscopy
• A map of the large day-night temperature gradient of a super-Earth exoplanet
• Spitzer Secondary Eclipses of the Dense, Modestly-irradiated, Giant Exoplanet HAT-P-20b Using Pixel-Level Decorrelation
• Thermal Phase Variations of WASP-12b: Defying Predictions
• Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b
• Transmission Spectroscopy of the Hot-Jupiter WASP-12b from 0.7 to 5 microns
• Transit and Eclipse Analyses of Exoplanet HD 149026b Using BLISS Mapping
• Variability in the super-Earth 55 Cnc e
• New Analysis Indicates No Thermal Inversion in the Atmosphere of HD 209458b
• The HARPS-N Rocky Planet Search I. HD219134b: A transiting rocky planet in a multi-planet system at 6.5 pc from the Sun
• The Kepler-19 System: A Transiting 2.2 R_Earth Planet and a Second Planet Detected via Transit Timing Variations
• 3.6 and 4.5 $\mu$m $Spitzer$ Phase Curves of the Highly-Irradiated Hot Jupiters WASP-19b and HAT-P-7b
• A global analysis of Spitzer and new HARPS data confirms the loneliness and metal-richness of GJ 436 b
• Deciphering the Atmospheric Composition of WASP-12b: A Comprehensive Analysis of its Dayside Emission
• Spitzer observations of the thermal emission from WASP-43b
• WASP-80b has a dayside within the T-dwarf range
• 3.6 and 4.5 $\mu$m Phase Curves of the Highly-Irradiated Eccentric Hot Jupiter WASP-14b
• Spitzer Secondary Eclipses of WASP-18b
• A Search for a Sub-Earth Sized Companion to GJ 436 and a Novel Method to Calibrate Warm Spitzer IRAC Observations
• Features in the Broadband Eclipse Spectra of Exoplanets: Signal or Noise?
• WASP-8b: Characterization of a Cool and Eccentric Exoplanet with Spitzer
• A Search for Water in the Atmosphere of HAT-P-26b Using LDSS-3C
• Spitzer/MIPS 24 micron Observations of HD 209458b: 3 eclipses, 2.5 transits, and a Phase Curve Corrupted by Instrumental Sensitivity Variations
• Thermal Emission of WASP-14b Revealed with Three Spitzer Eclipses
• Of `Cocktail Parties' and Exoplanets
• Two nearby sub-Earth-sized exoplanet candidates in the GJ 436 system
• Search for a habitable terrestrial planet transiting the nearby red dwarf GJ 1214
• Tests of the planetary hypothesis for PTFO 8-8695b
• A new look at Spitzer primary transit observations of the exoplanet HD189733b
• Warm Spitzer Occultation Photometry of WASP-26b at 3.6μm and 4.5μm

Cited by in corpus (7)

• Detection of a Westward Hotspot Offset in the Atmosphere of a Hot Gas Giant CoRoT-2b
• Mass Loss from the Exoplanet WASP-12b Inferred from $\text{Spitzer}$ Phase Curves
• A Comprehensive Reanalysis of $\text{Spitzer}$'s 4.5 $\mu$m Phase Curves, and the Phase Variations of the Ultra-hot Jupiters MASCARA-1b and KELT-16b
• Zodiacal Exoplanets in Time (ZEIT) IX: a flat transmission spectrum and a highly eccentric orbit for the young Neptune K2-25b as revealed by Spitzer
• Investigating Trends in Atmospheric Compositions of Cool Gas Giant Planets Using Spitzer Secondary Eclipses
• Smaller than expected bright-spot offsets in Spitzer phase curves of the hot Jupiter Qatar-1b
• A High Geometric Albedo for LTT9779b Points Towards a Metal-rich Atmosphere and Silicate Clouds