We present a new method employing machine learning techniques for measuring astrophysical features by correcting systematics in IRAC high precision photometry using Random Forests. The main systematic in IRAC light curve data is position changes due to unavoidable telescope motions coupled with an intrapixel response function. We aim to use the large amount of publicly available calibration data for the single pixel used for this type of work (the sweet spot pixel) to make a fast, easy to use, accurate correction to science data. This correction on calibration data has the advantage of using an independent dataset instead of using the science data on itself, which has the disadvantage of including astrophysical variations. After focusing on feature engineering and hyperparameter optimization, we show that a boosted random forest model can reduce the data such that we measure the median of ten archival eclipse observations of XO-3b to be 1459 +- 200 parts per million. This is a comparable depth to the average of those in the literature done by seven different methods, however the spread in measurements is 30-100% larger than those literature values, depending on the reduction method. We also caution others attempting similar methods to check their results with the fiducial dataset of XO-3b as we were also able to find models providing initially great scores on their internal test datasets but whose results significantly underestimated the eclipse depth of that planet.

We present a transmission spectrum for the Neptune-size exoplanet HD 106315 c from optical to infrared wavelengths based on transit observations from the Hubble Space Telescope/Wide Field Camera 3, K2, and Spitzer. The spectrum shows tentative evidence for a water absorption feature in the $1.1 - 1.7\mu$m wavelength range with a small amplitude of 30 ppm (corresponding to just $0.8 \pm 0.04$ atmospheric scale heights). Based on an atmospheric retrieval analysis, the presence of water vapor is tentatively favored with a Bayes factor of 1.7 - 2.6 (depending on prior assumptions). The spectrum is most consistent with either enhanced metallicity, high altitude condensates, or both. Cloud-free solar composition atmospheres are ruled out at $>5\sigma$ confidence. We compare the spectrum to grids of cloudy and hazy forward models and find that the spectrum is fit well by models with moderate cloud lofting or haze formation efficiency, over a wide range of metallicities ($1 - 100\times$ solar). We combine the constraints on the envelope composition with an interior structure model and estimate that the core mass fraction is $\gtrsim0.3$. With a bulk composition reminiscent of that of Neptune and an orbital distance of 0.15 AU, HD 106315 c hints that planets may form out of broadly similar material and arrive at vastly different orbits later in their evolution.

We report improved masses, radii, and densities for four planets in two bright M-dwarf systems, K2-3 and GJ3470, derived from a combination of new radial velocity and transit observations. Supplementing K2 photometry with follow-up Spitzer transit observations refined the transit ephemerides of K2-3 b, c, and d by over a factor of 10. We analyze ground-based photometry from the Evryscope and Fairborn Observatory to determine the characteristic stellar activity timescales for our Gaussian Process fit, including the stellar rotation period and activity region decay timescale. The stellar rotation signals for both stars are evident in the radial velocity data and are included in our fit using a Gaussian process trained on the photometry. We find the masses of K2-3 b, K2-3 c and GJ3470 b to be 6.48$^{+0.99}_{-0.93}$, 2.14$^{+1.08}_{-1.04}$, and 12.58$^{+1.31}_{-1.28}$ M$_\oplus$ respectively. K2-3 d was not significantly detected and has a 3-$\sigma$ upper limit of 2.80 M$_\oplus$. These two systems are training cases for future TESS systems; due to the low planet densities ($\rho$ $<$ 3.7 g cm$^{-3}$) and bright host stars (K $<$ 9 mag), they are among the best candidates for transmission spectroscopy in order to characterize the atmospheric compositions of small planets.

We present precision 4.5 $\mu$m Spitzer transit photometry of eight planet candidates discovered by the K2 mission: K2-52 b, K2-53 b, EPIC 205084841.01, K2-289 b, K2-174 b, K2-87 b, K2-90 b, and K2-124 b. The sample includes four sub-Neptunes and two sub-Saturns, with radii between 2.6 and 18 $R_\oplus$, and equilibrium temperatures between 440 and 2000 K. In this paper we identify several targets of potential interest for future characterization studies, demonstrate the utility of transit follow-up observations for planet validation and ephemeris refinement, and present new imaging and spectroscopy data. Our simultaneous analysis of the K2 and Spitzer light curves yields improved estimates of the planet radii, and multi-wavelength information which help validate their planetary nature, including the previously un-validated candidate EPIC 205686202.01 (K2-289 b). Our Spitzer observations yield an order of magnitude increase in ephemeris precision, thus paving the way for efficient future study of these interesting systems by reducing the typical transit timing uncertainty in mid-2021 from several hours to a dozen or so minutes. K2-53 b, K2-289 b, K2-174 b, K2-87 b, and K2-90 b are promising radial velocity (RV) targets given the performance of spectrographs available today or in development, and the M3V star K2-124 hosts a temperate sub-Neptune that is potentially a good target for both RV and atmospheric characterization studies.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations.

K2-55b is a Neptune-sized planet orbiting a K7 dwarf with a radius of $0.715^{+0.043}_{-0.040}R_\odot$, a mass of $0.688\pm0.069 M_\odot$, and an effective temperature of $4300^{+107}_{-100}$K. Having characterized the host star using near-infrared spectra obtained at IRTF/SpeX, we observed a transit of K2-55b with Spitzer/IRAC and confirmed the accuracy of the original K2 ephemeris for future follow-up transit observations. Performing a joint fit to the Spitzer/IRAC and K2 photometry, we found a planet radius of $4.41^{+0.32}_{-0.28} R_\oplus$, an orbital period of $2.84927265_{-6.42\times10^{-6}}^{+6.87\times10^{-6}}$ days, and an equilibrium temperature of roughly 900K. We then measured the planet mass by acquiring twelve radial velocity (RV) measurements of the system using HIRES on the 10m Keck I Telescope. Our RV data set precisely constrains the mass of K2-55b to $43.13^{+5.98}_{-5.80} M_\oplus$, indicating that K2-55b has a bulk density of $2.8_{-0.6}^{+0.8}$ g cm$^{-3}$ and can be modeled as a rocky planet capped by a modest H/He envelope ($M_{\rm envelope} = 12\pm3\% M_p$). K2-55b is denser than most similarly sized planets, raising the question of whether the high planetary bulk density of K2-55b could be attributed to the high metallicity of K2-55. The absence of a substantial volatile envelope despite the large mass of K2-55b poses a challenge to current theories of gas giant formation. We posit that K2-55b may have escaped runaway accretion by migration, late formation, or inefficient core accretion or that K2-55b was stripped of its envelope by a late giant impact.

While planets between the size of Uranus and Saturn are absent within the Solar System, the star K2-24 hosts two such planets, K2-24b and c, with radii equal to $5.4~R_E$ and $7.5~R_E$, respectively. The two planets have orbital periods of 20.9 days and 42.4 days, residing only 1% outside the nominal 2:1 mean-motion resonance. In this work, we present results from a coordinated observing campaign to measure planet masses and eccentricities that combines radial velocity (RV) measurements from Keck/HIRES and transit-timing measurements from K2 and Spitzer. K2-24b and c have low, but non-zero, eccentricities of $e_1 \sim e_2 \sim 0.08$. The low observed eccentricities provide clues regarding the formation and dynamical evolution of K2-24b and K2-24c, suggesting that they could be the result of stochastic gravitational interactions with a turbulent protoplanetary disk, among other mechanisms. K2-24b and c are $19\pm2~M_E$ and $15\pm2~M_E$, respectively; K2-24c is 20% less massive than K2-24b, despite being 40% larger. Their large sizes and low masses imply large envelope fractions, which we estimate at $26^{+3}_{-3}\%$ and $52^{+5}_{-3}\%$. In particular, K2-24c's large envelope presents an intriguing challenge to the standard model of core nucleated accretion that predicts the onset of runaway accretion when $f_{env} \approx 50\%$.

We present a new {\em Spitzer} transit observation of K2$-$28b, a sub-Neptune ($R_{\rm p} = 2.45\pm0.28 R_{Earth}$) orbiting a relatively bright ($V_{\rm mag} = 16.06$, $K_{\rm mag} = 10.75$) metal-rich M4 dwarf (EPIC 206318379). This star is one of only seven with masses less than 0.2 M_{Sun} known to host transiting planets, and the planet appears to be a slightly smaller analogue of GJ 1214b ($2.85\pm0.20 R_{Earth}$; \citealt{2013A&A...549A..10H}). Our new {\em Spitzer} observations were taken two years after the original K2 discovery data and have a significantly higher cadence, allowing us to derive improved estimates for this planet's radius, semi-major axis, and orbital period, which greatly reduce the uncertainty in the prediction of near future transit times for the {\em James Webb Space Telescope} ({\em JWST}) observations. We also evaluate the system's suitability for atmospheric characterization with {\em JWST} and find that it is currently the only small (< 3 R_{Earth}) and cool ($< 600$ K) planet aside from GJ 1214b with a potentially detectable secondary eclipse. We also note that this system is a favorable target for near-infrared radial velocity instruments on larger telescopes (e.g., the Habitable Planet Finder on the Hobby-Eberly Telescope), making it one of only a handful of small, cool planets accessible with this technique. Finally, we compare our results with the simulated catalog of the Transiting Exoplanet Survey Satellite ({\em TESS}) and find K2$-$28b to be representative of the kind of mid-M systems that should be detectible in the {\em TESS} sample.

We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the prospector fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). We find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. In the future, our approach can be applied to a large and diverse sample of UDGs down to faint surface brightness limits, with the goal of constraining their stellar ages, stellar metallicities, and circumstellar and diffuse interstellar dust content.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.

We study the mass-richness relation of 116 spectroscopically-confirmed massive clusters at $0.4 < z < 2$ by mining the $Spitzer$ archive. We homogeneously measure the richness at 4.5$\mu$m for our cluster sample within a fixed aperture of $2^{\prime}$ radius and above a fixed brightness threshold, making appropriate corrections for both background galaxies and foreground stars. We have two subsamples, those which have a) literature X-ray luminosities and b) literature Sunyaev-Zeldovich effect masses. For the X-ray subsample we re-derive masses adopting the most recent calibrations. We then calibrate an empirical mass-richness relation for the combined sample spanning more than one decade in cluster mass and find the associated uncertainties in mass at fixed richness to be $\pm 0.25$ dex. We study the dependance of the scatter of this relation with galaxy concentration, defined as the ratio between richness measured within an aperture radius of 1 and 2 arcminutes. We find that at fixed aperture radius the scatter increases for clusters with higher concentrations. We study the dependance of our richness estimates with depth of the [4.5]$\mu$m imaging data and find that reaching a depth of at least [4.5]= 21 AB mag is sufficient to derive reasonable mass estimates. We discuss the possible extension of our method to the mid-infrared $WISE$ all-sky survey data, and the application of our results to the $Euclid$ mission. This technique makes richness-based cluster mass estimates available for large samples of clusters at very low observational cost.

We measure the 4.5 $\mu$m thermal emission of five transiting hot Jupiters, WASP-13b, WASP-15b, WASP-16b, WASP-62b and HAT-P-22b using channel 2 of the Infrared Array Camera (IRAC) on the {\sl Spitzer Space Telescope}. Significant intrapixel sensitivity variations in Spitzer IRAC data require careful correction in order to achieve precision on the order of several hundred parts per million (ppm) for the measurement of exoplanet secondary eclipses. We determine eclipse depths by first correcting the raw data using three independent data reduction methods. The Pixel Gain Map (PMAP), Nearest Neighbors (NNBR), and Pixel Level Decorrelation (PLD) each correct for the intrapixel sensitivity effect in Spitzer photometric time-series observations. The results from each methodology are compared against each other to establish if they reach a statistically equivalent result in every case and to evaluate their ability to minimize uncertainty in the measurement. We find that all three methods produce reliable results. For every planet examined here NNBR and PLD produce results that are in statistical agreement. However, the PMAP method appears to produce results in slight disagreement in cases where the stellar centroid is not kept consistently on the most well characterized area of the detector. We evaluate the ability of each method to reduce the scatter in the residuals as well as in the correlated noise in the corrected data. The NNBR and PLD methods consistently minimize both white and red noise levels and should be considered reliable and consistent. The planets in this study span equilibrium temperatures from 1100 to 2000~K and have brightness temperatures that require either high albedo or efficient recirculation. However, it is possible that other processes such as clouds or disequilibrium chemistry may also be responsible for producing these brightness temperatures.

The recent detections of two transit events attributed to the super-Earth candidate K2-18b have provided the unprecedented prospect of spectroscopically studying a habitable-zone planet outside the Solar System. Orbiting a nearby M2.5 dwarf and receiving virtually the same stellar insolation as Earth, K2-18b would be a prime candidate for the first detailed atmospheric characterization of a habitable-zone exoplanet using HST and JWST. Here, we report the detection of a third transit of K2-18b near the predicted transit time using the Spitzer Space Telescope. The Spitzer detection demonstrates the periodic nature of the two transit events discovered by K2, confirming that K2-18 is indeed orbited by a super-Earth in a 33-day orbit and ruling out the alternative scenario of two similarly-sized, long-period planets transiting only once within the 75-day K2 observation. We also find, however, that the transit event detected by Spitzer occurred 1.85 hours (7-sigma) before the predicted transit time. Our joint analysis of the Spitzer and K2 photometry reveals that this early occurrence of the transit is not caused by transit timing variations (TTVs), but the result of an inaccurate K2 ephemeris due to a previously undetected data anomaly in the K2 photometry likely caused by a cosmic ray hit. We refit the ephemeris and find that K2-18b would have been lost for future atmospheric characterizations with HST and JWST if we had not secured its ephemeris shortly after the discovery. We caution that immediate follow-up observations as presented here will also be critical in confirming and securing future planets discovered by TESS, in particular if only two transit events are covered by the relatively short 27-day TESS campaigns.

We have used the Spitzer Space Telescope in February 2016 to obtain high cadence, high signal-to-noise, 17-hour duration light curves of Neptune at 3.6 and 4.5 $\mu$m. The light curve duration was chosen to correspond to the rotation period of Neptune. Both light curves are slowly varying with time, with full amplitudes of 1.1 mag at 3.6 $\mu$m and 0.6 mag at 4.5 $\mu$m. We have also extracted sparsely sampled 18-hour light curves of Neptune at W1 (3.4 $\mu$m) and W2 (4.6 $\mu$m) from the WISE/NEOWISE archive at six epochs in 2010-2015. These light curves all show similar shapes and amplitudes compared to the Spitzer light curves but with considerable variation from epoch to epoch. These amplitudes are much larger than those observed with Kepler/K2 in the visible (amplitude $\sim$0.02 mag) or at 845 nm with the Hubble Space Telescope in 2015 and at 763 nm in 2016 (amplitude $\sim$ 0.2 mag). We interpret the Spitzer and WISE light curves as arising entirely from reflected solar photons, from higher levels in Neptune's atmosphere than for K2. Methane gas is the dominant opacity source in Neptune's atmosphere, and methane absorption bands are present in the HST 763, and 845 nm, WISE W1, and Spitzer 3.6 $\mu$m filters.

We examine the repeatability, reliability, and accuracy of differential exoplanet eclipse depth measurements made using the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope during the post-cryogenic mission. We have re-analyzed an existing 4.5 {\mu}m data set, consisting of 10 observations of the XO-3b system during secondary eclipse, using seven different techniques for removing correlated noise. We find that, on average, for a given technique, the eclipse depth estimate is repeatable from epoch to epoch to within 156 parts per million (ppm). Most techniques derive eclipse depths that do not vary by more than a factor 3 of the photon noise limit. All methods but one accurately assess their own errors: for these methods, the individual measurement uncertainties are comparable to the scatter in eclipse depths over the 10 epoch sample. To assess the accuracy of the techniques as well as to clarify the difference between instrumental and other sources of measurement error, we have also analyzed a simulated data set of 10 visits to XO-3b, for which the eclipse depth is known. We find that three of the methods (BLISS mapping, Pixel Level Decorrelation, and Independent Component Analysis) obtain results that are within three times the photon limit of the true eclipse depth. When averaged over the 10 epoch ensemble, 5 out of 7 techniques come within 60 ppm of the true value. Spitzer exoplanet data, if obtained following current best practices and reduced using methods such as those described here, can measure repeatable and accurate single eclipse depths, with close to photon-limited results.

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

Over the past decade, observations of giant exoplanets (Jupiter-size) have provided key insights into their atmospheres, but the properties of lower-mass exoplanets (sub-Neptune) remain largely unconstrained because of the challenges of observing small planets. Numerous efforts to observe the spectra of super-Earths (exoplanets with masses of one to ten times that of Earth) have so far revealed only featureless spectra. Here we report a longitudinal thermal brightness map of the nearby transiting super-Earth 55 Cancri e revealing highly asymmetric dayside thermal emission and a strong day-night temperature contrast. Dedicated space-based monitoring of the planet in the infrared revealed a modulation of the thermal flux as 55 Cancri e revolves around its star in a tidally locked configuration. These observations reveal a hot spot that is located 41 +- 12 degrees east of the substellar point (the point at which incident light from the star is perpendicular to the surface of the planet). From the orbital phase curve, we also constrain the nightside brightness temperature of the planet to 1380 +- 400 kelvin and the temperature of the warmest hemisphere (centred on the hot spot) to be about 1300 kelvin hotter (2700 +- 270 kelvin) at a wavelength of 4.5 microns, which indicates inefficient heat redistribution from the dayside to the nightside. Our observations are consistent with either an optically thick atmosphere with heat recirculation confined to the planetary dayside, or a planet devoid of atmosphere with low-viscosity magma flows at the surface.

Motivated by a high Spitzer IRAC oversubscription rate, we present a new technique of randomly and sparsely sampling phase curves of hot Jupiters. Snapshot phase curves are enabled by technical advances in precision pointing as well as careful characterization of a portion of the central pixel on the array. This method allows for observations which are a factor of roughly two more efficient than full phase curve observations, and are furthermore easier to insert into the Spitzer observing schedule. We present our pilot study from this program using the exoplanet WASP-14b. Data of this system were taken both as a sparsely sampled phase curve as well as a staring mode phase curve. Both datasets as well as snapshot style observations of a calibration star are used to validate this technique. By fitting our WASP-14b phase snapshot dataset, we successfully recover physical parameters for the transit and eclipse depths as well as amplitude and maximum and minimum of the phase curve shape of this slightly eccentric hot Jupiter. We place a limit on the potential phase to phase variation of these parameters since our data are taken over many phases over the course of a year. We see no evidence for eclipse depth variations compared to other published WASP-14b eclipse depths over a 3.5 year baseline.

We have used the {\it Spitzer Space Telescope} to observe two transiting planetary systems orbiting low mass stars discovered in the \Kepler \Ktwo mission. The system K2-3 (EPIC 201367065) hosts three planets while EPIC 202083828 (K2-26) hosts a single planet. Observations of all four objects in these two systems confirm and refine the orbital and physical parameters of the planets. The refined orbital information and more precise planet radii possible with \Spitzer will be critical for future observations of these and other \Ktwo targets. For K2-3b we find marginally significant evidence for a Transit Timing Variation between the \Ktwo and \Spitzer\ epochs.

This article summarizes a workshop held on March, 2014, on the potential of the James Webb Space Telescope (JWST) to revolutionize our knowledge of the physical properties of exoplanets through transit observations. JWST's unique combination of high sensitivity and broad wavelength coverage will enable the accurate measurement of transits with high signal-to-noise. Most importantly, JWST spectroscopy will investigate planetary atmospheres to determine atomic and molecular compositions, to probe vertical and horizontal structure, and to follow dynamical evolution, i.e. exoplanet weather. JWST will sample a diverse population of planets of varying masses and densities in a wide variety of environments characterized by a range of host star masses and metallicities, orbital semi-major axes and eccentricities. A broad program of exoplanet science could use a substantial fraction of the overall JWST mission.

The dominant non-instrumental background source for space-based infrared observatories is the zo- diacal light. We present Spitzer Infrared Array Camera (IRAC) measurements of the zodiacal light at 3.6, 4.5, 5.8, and 8.0 {\mu}m, taken as part of the instrument calibrations. We measure the changing surface brightness levels in approximately weekly IRAC observations near the north ecliptic pole (NEP) over the period of roughly 8.5 years. This long time baseline is crucial for measuring the annual sinusoidal variation in the signal levels due to the tilt of the dust disk with respect to the ecliptic, which is the true signal of the zodiacal light. This is compared to both Cosmic Background Explorer Diffuse Infrared Background Experiment (COBE DIRBE) data and a zodiacal light model based thereon. Our data show a few percent discrepancy from the Kelsall et al. (1998) model including a potential warping of the interplanetary dust disk and a previously detected overdensity in the dust cloud directly behind the Earth in its orbit. Accurate knowledge of the zodiacal light is important for both extragalactic and Galactic astronomy including measurements of the cosmic infrared background, absolute measures of extended sources, and comparison to extrasolar interplanetary dust models. IRAC data can be used to further inform and test future zodiacal light models.

We present a detailed strong lensing, weak lensing and X-ray analysis of Abell 2744 (z = 0.308), one of the most actively merging galaxy clusters known. It appears to have unleashed `dark', `ghost', `bullet' and `stripped' substructures, each ~10^14 solar masses. The phenomenology is complex and will present a challenge for numerical simulations to reproduce. With new, multiband HST imaging, we identify 34 strongly-lensed images of 11 galaxies around the massive Southern `core'. Combining this with weak lensing data from HST, VLT and Subaru, we produce the most detailed mass map of this cluster to date. We also perform an independent analysis of archival Chandra X-ray imaging. Our analyses support a recent claim that the Southern core and Northwestern substructure are post-merger and exhibit morphology similar to the Bullet Cluster viewed from an angle. From the separation between X-ray emitting gas and lensing mass in the Southern core, we derive a new and independent constraint on the self-interaction cross section of dark matter particles sigma/m <~ 3 \pm 1 cm^2 g^-1. In the Northwestern substructure, the gas, dark matter, and galaxy components have become separated by much larger distances. Most curiously, the `ghost' clump (primarily gas) leads the `dark' clump (primarily dark matter) by more than 150 kpc. We propose an enhanced `ram-pressure slingshot' scenario which may have yielded this reversal of components with such a large separation, but needs further confirmation by follow-up observations and numerical simulations. A secondary merger involves a second `bullet' clump in the North and an extremely `stripped' clump to the West. The latter appears to exhibit the largest separation between dark matter and X-ray emitting baryons detected to date in our sky.

We present 3.6 and 4.5 micron Spitzer IRAC imaging over 0.77 square degrees at the Virgo cluster core for the purpose of understanding the formation mechanisms of the low surface brightness intracluster light features. Instrumental and astrophysical backgrounds that are hundreds of times higher than the signal were carefully characterized and removed. We examine both intracluster light plumes as well as the outer halo of the giant elliptical M87. For two intracluster light plumes, we use optical colors to constrain their ages to be greater than 3 & 5 Gyr, respectively. Upper limits on the IRAC fluxes constrain the upper limits to the masses, and optical detections constrain the lower limits to the masses. In this first measurement of mass of intracluster light plumes we find masses in the range of 5.5 x 10^8 - 4.5 x 10^9 and 2.1 x 10^8 - 1.5 x 10^9 solar masses for the two plumes for which we have coverage. Given their expected short lifetimes, and a constant production rate for these types of streams, integrated over Virgo's lifetime, they can account for the total ICL content of the cluster implying that we do not need to invoke ICL formation mechanisms other than gravitational mechanisms leading to bright plumes. We also examined the outer halo of the giant elliptical M87. The color profile from the inner to outer halo of M87 (160 Kpc) is consistent with either a flat or optically blue gradient, where a blue gradient could be due to younger or lower metallicity stars at larger radii. The similarity of the age predicted by both the infrared and optical colors (> few Gyr) indicates that the optical measurements are not strongly affected by dust extinction.

We present 20 band photometry from the far-IR to X-ray in the Spitzer IRAC dark field. The bias for the near-IR camera on Spitzer is calibrated by observing a ~20 arcminute diameter "dark" field near the north ecliptic pole roughly every two-to-three weeks throughout the mission duration of Spitzer. The field is unique for its extreme depth, low background, high quality imaging, time-series information, and accompanying photometry including data taken with Akari, Palomar, MMT, KPNO, Hubble, and Chandra. This serendipitous survey contains the deepest mid-IR data taken to date. This dataset is well suited for studies of intermediate redshift galaxy clusters, high redshift galaxies, the first generation of stars, and the lowest mass brown dwarfs, among others. This paper provides a summary of the data characteristics and catalog generation from all bands collected to date as well as a discussion of photometric redshifts and initial and expected science results and goals. To illustrate the scientific potential of this unique dataset, we also present here IRAC color color diagrams.

We have undertaken a systematic search for candidate supernovae from high-redshift Population III stars in a field that has been observed with repeated imaging on a cadence of 2-3 weeks over a 2.2 year baseline, the Spitzer/IRAC Dark Field. The individual epochs reach a typical 5-sigma depth of 1 uJy in IRAC Channel 1 (3.6 um). Requiring a minimum of four epochs coverage, the total effective area searched is 214 sq arcminutes. The unprecedented depth and multi-epochal nature of these data make it ideal for a first foray to detect transient objects which may be candidate luminous Pair Instability Supernovae from the primordial-metallicity first stars. The search was conducted over a broad range of timescales, allowing for different durations of the putative candidates' light curve plateau phases. All candidates were vetted by inspection of the Spitzer imaging data, as well as deep HST/ACS F814W imaging available over the full field. While many resolved-source objects were found with Spitzer variability, no transient objects could plausibly be identified as high-redshift supernovae candidates. The resulting 95% confidence level upper limit is 23/sq deg/yr, for sources with plateau timescales under 400/(1+z) days and brightnesses above ~1 uJy.

We present infrared luminosities, star formation rates, colors, morphologies, locations, and AGN properties of 24 micron-detected sources in photometrically detected high-redshift clusters in order to understand the impact of environment on star formation and AGN evolution in cluster galaxies. We use three newly-identified z=1 clusters selected from the IRAC dark field; the deepest ever mid-IR survey with accompanying, 14 band multiwavelength data including deep HST imaging and deep wide-area Spitzer MIPS 24 micron imaging. We find 90 cluster members with MIPS detections within two virial radii of the cluster centers, of which 17 appear to have spectral energy distributions dominated by AGN and the rest dominated by star formation. We find that 43 of the star forming are luminous infrared galaxies (LIRGs). The majority of sources (81%) are spirals or irregulars. A large fraction (at least 25%) show obvious signs of interactions. The MIPS -detected member galaxies have varied spatial distributions as compared to the MIPS-undetected members with one of the three clusters showing SF galaxies being preferentially located on the cluster outskirts, while the other 2 clusters show no such trend. Both the AGN fraction and the summed SFR of cluster galaxies increases from z=0 to 1, at a rate that is a few times faster in clusters than over the same redshift range in the field. Cluster environment does have an effect on the evolution of both AGN fraction and SFR from redshift one to the present, but does not effect the infrared luminosities or morphologies of the MIPS sample. Star formation happens in the same way regardless of environment making MIPS sources look the same in the cluster and field, however the cluster environment does encourage a more rapid evolution with time as compared to the field.

Using three newly identified galaxy clusters at z~1 (photometric redshift) we measure the evolution of the galaxies within clusters from high redshift to the present day by studying the growth of the red cluster sequence. The clusters are located in the Spitzer Infrared Array Camera (IRAC) Dark Field, an extremely deep mid-infrared survey near the north ecliptic pole with photometry in 18 total bands from X-ray through far-IR. Two of the candidate clusters are additionally detected as extended emission in matching Chandra data in the survey area allowing us to measure their masses to be M_{500}= 6.2 \pm 1.0 \times 10^{13} and 3.6 \pm 1.1 \times 10^{13} solar masses. For all three clusters we create a composite color magnitude diagram in rest-frame B-K using our deep HST and Spitzer imaging. By comparing the fraction of low luminosity member galaxies on the composite red sequence with the corresponding population in local clusters at z=0.1 taken from the COSMOS survey, we examine the effect of a galaxy's mass on its evolution. We find a deficit of faint galaxies on the red sequence in our z~1 clusters which implies that more massive galaxies have evolved in clusters faster than less massive galaxies, and that the less massive galaxies are still forming stars in clusters such that they have not yet settled onto the red sequence.

We have measured the flux, profile, color, and substructure in the diffuse intracluster light (ICL) in a sample of ten galaxy clusters with a range of mass, morphology, redshift, and density. Deep, wide-field observations for this project were made in two bands at the one meter Swope and 2.5 meter du Pont telescope at Las Campanas Observatory. Careful attention in reduction and analysis was paid to the illumination correction, background subtraction, point spread function determination, and galaxy subtraction. ICL flux is detected in both bands in all ten clusters ranging from 7.6 x 10^{10} to 7.0 x 10^{11} h^{-1} solar luminosities in r and 1.4 x 10^{10} to 1.2 x 10^{11} h^{-1} solar luminosities in the B-band. These fluxes account for 6 to 22% of the total cluster light within one quarter of the virial radius in r and 4 to 21% in the B-band. Average ICL B-r colors range from 1.5 to 2.8 mags when k and evolution corrected to the present epoch. In several clusters we also detect ICL in group environments near the cluster center and up to ~1 h^{-1} Mpc distant from the cluster center. Our sample, having been selected from the Abell sample, is incomplete in that it does not include high redshift clusters with low density, low flux, or low mass, and it does not include low redshift clusters with high flux, mass, or density. This bias makes it difficult to interpret correlations between ICL flux and cluster properties. Despite this selection bias, we do find that the presence of a cD galaxy corresponds to both centrally concentrated galaxy profiles and centrally concentrated ICL profiles. This is consistent with ICL either forming from galaxy interactions at the center, or forming at earlier times in groups and later combining in the center.

We are undertaking a program to measure the characteristics of the intracluster light (total flux, profile, color, and substructure) in a sample of 10 galaxy clusters with a range of cluster mass, morphology, and redshift. We present here the methods and results for the first cluster in that sample, A3888. We have identified an intracluster light (ICL) component in A3888 in V and r that contains 13\pm5% of the total cluster light and extends to 700h_{70}^{-1}kpc (~0.3 r_{200}) from the center of the cluster. The ICL color in our smallest radial bin is V-r = 0.3 \pm 0.1, similar to the central cluster ellipticals. The ICL is redder than the galaxies at 400 < r < 700h_{70}^{-1}kpc although the uncertainty in any one radial bin is high. Based on a comparison of V-r color with simple stellar models, the ICL contains a component which formed more than 7 Gyr ago (at z > 1), coupled with a high metallicity (1.0Z_{\odot} < Z_{ICL} \la 2.5Z_{\odot}), and a more centralized component which contains stars formed within the past 5 Gyr (at z ~ 1). The profile of the ICL can be roughly fit by a shallow exponential in the outer regions and a steeper exponential in the central region. We also find a concentration of diffuse light around a small group of galaxies 1.4h_{70}^{-1}Mpc from the center of the cluster. In addition, we find 3 low surface brightness features near the cluster center which are blue (V-r = 0.0) and contain a total flux of 0.1M*. Based on these observations and X-ray and galaxy morphology, we suggest that this cluster is entering a phase of significant merging of galaxy groups in the core, whereupon we expect the ICL fraction to grow significantly with the formation of a cD galaxy as well as the in-fall of groups.

Observations of extragalactic objects need to be corrected for Galactic absorption and this is often accomplished by using the measured 21 cm HI column. However, within the beam of the radio telescope there are variations in the HI column that can have important effects in interpreting absorption line studies and X-ray spectra at the softest energies. We examine the HI and DIRBE/IRAS data for lines of sight out of the Galaxy, which show evidence for HI variations in of up to a factor of three in 1 degree fields. Column density enhancements would preferentially absorb soft X-rays in spatially extended objects and we find evidence for this effect in the ROSAT PSPC observations of two bright clusters of galaxies, Abell 119 and Abell 2142. For clusters of galaxies, the failure to include column density fluctuations will lead to systematically incorrect fits to the X-ray data in the sense that there will appear to be a very soft X-ray excess. This may be one cause of the soft X-ray excess in clusters, since the magnitude of the effect is comparable to the observed values.