We present MCSED, a new spectral energy distribution (SED)-fitting code, which mates flexible stellar evolution calculations with the Markov Chain Monte Carlo algorithms of the software package emcee. MCSED takes broad, intermediate, and narrow-band photometry, emission-line fluxes, and/or absorption line spectral indices, and returns probability distributions and co-variance plots for all model parameters. MCSED includes a variety of dust attenuation curves with parameters for varying the UV slopes and bump strengths, a prescription for continuum and PAH emission from dust, models for continuum and line emission from ionized gas, options for fixed and variable stellar metallicity, and a selection of star formation rate (SFR) histories. The code is well-suited for exploring parameter inter-dependencies in sets of galaxies with known redshifts, for which there is multi-band photometry and/or spectroscopy. We apply MCSED to a sample of $\sim2000$ $1.90<z<2.35$ galaxies in the five CANDELS fields, which were selected via their strong [O III] $\lambda5007$ emission, and explore the systematic behavior of their SEDs. We find the galaxies become redder with stellar mass, due to both increasing internal attenuation and a greater population of older stars. The slope of the UV extinction curve also changes with stellar mass, and at least some galaxies exhibit an extinction excess at 2175 Angstroms. Finally, we demonstrate that below $M\lesssim10^9\,M_{\odot}$), the shape of the star-forming galaxy main sequence is highly dependent on the galaxies' assumed SFR history, as calculations which assume a constant SFR produce stellar masses that are $\sim1$ dex smaller than those found using more realistic SFR histories.

Euclid, WFIRST, and HETDEX will make emission-line selected galaxies the largest observed constituent in the $z > 1$ universe. However, we only have a limited understanding of the physical properties of galaxies selected via their Ly$\alpha$ or rest-frame optical emission lines. To begin addressing this problem, we present the basic properties of $\sim 2,000$ AEGIS, COSMOS, GOODS-N, GOODS-S, and UDS galaxies identified in the redshift range $1.90 < z < 2.35$ via their [O II], H$\beta$, and [O III] emission lines. For these $z \sim 2$ galaxies, [O III] is generally much brighter than [O II] and H$\beta$, with typical rest-frame equivalent widths of several hundred Angstroms. Moreover, these strong emission-line systems span an extremely wide range of stellar mass ($\sim 3$ dex), star-formation rate ($\sim 2$ dex), and [O III] luminosity ($\sim 2$ dex). Comparing the distributions of these properties to those of continuum selected galaxies, we find that emission-line galaxies have systematically lower stellar masses and lower optical/UV dust attenuations. These measurements lay the groundwork for an extensive comparison between these rest-frame optical emission-line galaxies and Ly$\alpha$ emitters identified in the HETDEX survey.

We investigate the multiple stellar populations of the globular clusters M3, M5, M13, and M71 using $g^\prime$ and intermediate-band CN-$\lambda 3883$ photometry obtained with the WIYN 0.9-m telescope on Kitt Peak. We find a strong correlation between red giant stars' CN$-g^\prime$ colors and their spectroscopic sodium abundances, thus demonstrating the efficacy of the two-filter system for stellar population studies. In all four clusters, the observed spread in red giant branch CN$-g^\prime$ colors is wider than that expected from photometric uncertainty, confirming the well-known chemical inhomogeneity of these systems. M3 and M13 show clear evidence for a radial dependence in the CN-band strengths of its red giants, while the evidence for such a radial dependence of CN strengths in M5 is ambiguous. Our data suggest that the dynamically old, relatively metal-rich M71 system is well mixed, as it shows no evidence for chemical segregation. Finally, we measure the radial gradients in the integrated CN$-g^\prime$ color of the clusters and find that such gradients are easily detectable in the integrated light. We suggest that photometric observations of color gradients within globular clusters throughout the Local Group can be used to characterize their multiple populations, and thereby constrain the formation history of globular clusters in different galactic environments.

We compare the rest-frame ultraviolet and rest-frame optical morphologies of 2 < z < 3 star-forming galaxies in the GOODS-S field using Hubble Space Telescope WFC3 and ACS images from the CANDELS, GOODS, and ERS programs. We show that the distribution of sizes and concentrations for 1.90 < z < 2.35 galaxies selected via their rest-frame optical emission-lines are statistically indistinguishable from those of Lyman-alpha emitting systems found at z ~ 2.1 and z ~ 3.1. We also show that the z > 2 star-forming systems of all sizes and masses become smaller and more compact as one shifts the observing window from the UV to the optical. We argue that this offset is due to inside-out galaxy formation over the first ~ 2 Gyr of cosmic time.

Exoplanet WASP-14b is a highly irradiated, transiting hot Jupiter. Joshi et al. calculate an equilibrium temperature Teq of 1866 K for zero albedo and reemission from the entire planet, a mass of 7.3 +/- 0.5 Jupiter masses and a radius of 1.28 +/- 0.08 Jupiter radii. Its mean density of 4.6 g/cm3 is one of the highest known for planets with periods less than 3 days. We obtained three secondary eclipse light curves with the Spitzer Space Telescope. The eclipse depths from the best jointly fit model are $0.224\%$ +/- $0.018\%$ at 4.5 {\mu}m and $0.181\%$ +/- $0.022\%$ at 8.0 {\mu}m. The corresponding brightness temperatures are 2212 +/- 94 K and 1590 +/- 116 K. A slight ambiguity between systematic models suggests a conservative 3.6 {\mu}m eclipse depth of $0.19\%$ +/- $0.01\%$ and brightness temperature of 2242 +/- 55 K. Although extremely irradiated, WASP-14b does not show any distinct evidence of a thermal inversion. In addition, the present data nominally favor models with day night energy redistribution less than $~30\%$. The current data are generally consistent with oxygen-rich as well as carbon-rich compositions, although an oxygen-rich composition provides a marginally better fit. We confirm a significant eccentricity of e = 0.087 +/- 0.002 and refine other orbital parameters.

The dayside of HD 149026b is near the edge of detectability by the Spitzer Space Telescope. We report on eleven secondary-eclipse events at 3.6, 4.5, 3 x 5.8, 4 x 8.0, and 2 x 16 microns plus three primary-transit events at 8.0 microns. The eclipse depths from jointly-fit models at each wavelength are 0.040 +/- 0.003% at 3.6 microns, 0.034 +/- 0.006% at 4.5 microns, 0.044 +/- 0.010% at 5.8 microns, 0.052 +/- 0.006% at 8.0 microns, and 0.085 +/- 0.032% at 16 microns. Multiple observations at the longer wavelengths improved eclipse-depth signal-to-noise ratios by up to a factor of two and improved estimates of the planet-to-star radius ratio (Rp/Rs = 0.0518 +/- 0.0006). We also identify no significant deviations from a circular orbit and, using this model, report an improved period of 2.8758916 +/- 0.0000014 days. Chemical-equilibrium models find no indication of a temperature inversion in the dayside atmosphere of HD 149026b. Our best-fit model favors large amounts of CO and CO2, moderate heat redistribution (f=0.5), and a strongly enhanced metallicity. These analyses use BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping, a new technique to model two position-dependent systematics (intrapixel variability and pixelation) by mapping the pixel surface at high resolution. BLISS mapping outperforms previous methods in both speed and goodness of fit. We also present an orthogonalization technique for linearly-correlated parameters that accelerates the convergence of Markov chains that employ the Metropolis random walk sampler. The electronic supplement contains light-curve files and supplementary figures.

The transiting exoplanet WASP-18b was discovered in 2008 by the Wide Angle Search for Planets (WASP) project. The Spitzer Exoplanet Target of Opportunity Program observed secondary eclipses of WASP-18b using Spitzer's Infrared Array Camera (IRAC) in the 3.6 micron and 5.8 micron bands on 2008 December 20, and in the 4.5 micron and 8.0 micron bands on 2008 December 24. We report eclipse depths of 0.30 +/- 0.02%, 0.39 +/- 0.02%, 0.37 +/- 0.03%, 0.41 +/- 0.02%, and brightness temperatures of 3100 +/- 90, 3310 +/- 130, 3080 +/- 140 and 3120 +/- 110 K in order of increasing wavelength. WASP-18b is one of the hottest planets yet discovered - as hot as an M-class star. The planet's pressure-temperature profile most likely features a thermal inversion. The observations also require WASP-18b to have near-zero albedo and almost no redistribution of energy from the day-side to the night side of the planet.

The nearby extrasolar planet GJ 436b--which has been labelled as a 'hot Neptune'--reveals itself by the dimming of light as it crosses in front of and behind its parent star as seen from Earth. Respectively known as the primary transit and secondary eclipse, the former constrains the planet's radius and mass, and the latter constrains the planet's temperature and, with measurements at multiple wavelengths, its atmospheric composition. Previous work using transmission spectroscopy failed to detect the 1.4-\mu m water vapour band, leaving the planet's atmospheric composition poorly constrained. Here we report the detection of planetary thermal emission from the dayside of GJ 436b at multiple infrared wavelengths during the secondary eclipse. The best-fit compositional models contain a high CO abundance and a substantial methane (CH4) deficiency relative to thermochemical equilibrium models for the predicted hydrogen-dominated atmosphere. Moreover, we report the presence of some H2O and traces of CO2. Because CH4 is expected to be the dominant carbon-bearing species, disequilibrium processes such as vertical mixing and polymerization of methane into substances such as ethylene may be required to explain the hot Neptune's small CH4-to-CO ratio, which is at least 10^5 times smaller than predicted.

We observed two secondary eclipses of the exoplanet WASP-12b using the Infrared Array Camera on the Spitzer Space Telescope. The close proximity of WASP-12b to its G-type star results in extreme tidal forces capable of inducing apsidal precession with a period as short as a few decades. This precession would be measurable if the orbit had a significant eccentricity, leading to an estimate of the tidal Love number and an assessment of the degree of central concentration in the planetary interior. An initial ground-based secondary eclipse phase reported by Lopez-Morales et al. (0.510 +/- 0.002) implied eccentricity at the 4.5 sigma level. The spectroscopic orbit of Hebb et al. has eccentricity 0.049 +/- 0.015, a 3 sigma result, implying an eclipse phase of 0.509 +/- 0.007. However, there is a well documented tendency of spectroscopic data to overestimate small eccentricities. Our eclipse phases are 0.5010 +/- 0.0006 (3.6 and 5.8 microns) and 0.5006 +/- 0.0007 (4.5 and 8.0 microns). An unlikely orbital precession scenario invoking an alignment of the orbit during the Spitzer observations could have explained this apparent discrepancy, but the final eclipse phase of Lopez-Morales et al. (0.510 -0.006 / +0.007) is consistent with a circular orbit at better than 2 sigma. An orbit fit to all the available transit, eclipse, and radial-velocity data indicates precession at <1 sigma; a non-precessing solution fits better. We also comment on analysis and reporting for Spitzer exoplanet data in light of recent re-analyses.

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to a silicate-dominated composition as found on Earth, and the atmospheres can also differ from those in the Solar System. The solar C/O is 0.54. Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b that reveals C/O >= 1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapor and enhanced in methane by over two orders of magnitude each compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T > 2,500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has very efficient day- night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.

We report Spitzer/IRAC photometry of the transiting giant exoplanet HAT-P-1b during its secondary eclipse. This planet lies near the postulated boundary between the pM and pL-class of hot Jupiters, and is important as a test of models for temperature inversions in hot Jupiter atmospheres. We derive eclipse depths for HAT-P-1b, in units of the stellar flux, that are: 0.080% +/- 0.008%,[3.6um], 0.135% +/- 0.022%,[4.5um],0.203% +/- 0.031%,[5.8um], and $0.238% +/- 0.040%,[8.0um]. These values are best fit using an atmosphere with a modest temperature inversion, intermediate between the archetype inverted atmosphere (HD209458b) and a model without an inversion. The observations also suggest that this planet is radiating a large fraction of the available stellar irradiance on its dayside, with little available for redistribution by circulation. This planet has sometimes been speculated to be inflated by tidal dissipation, based on its large radius in discovery observations, and on a non-zero orbital eccentricity allowed by the radial velocity data. The timing of the secondary eclipse is very sensitive to orbital eccentricity, and we find that the central phase of the eclipse is 0.4999 +/- 0.0005. The difference between the expected and observed phase indicates that the orbit is close to circular, with a 3-sigma limit of |e cosw| < 0.002.

We report the results of infrared (8 micron) transit and secondary eclipse photometry of the hot Neptune exoplanet, GJ436b using Spitzer. The nearly photon-limited precision of these data allow us to measure an improved radius for the planet, and to detect the secondary eclipse. The transit (centered at HJD = 2454280.78149 +/-0.00016) shows the flat-bottomed shape typical of infrared transits, and it precisely defines the planet-to-star radius ratio (0.0839 +/-0.0005), independent of the stellar properties. However, we obtain the planetary radius, as well as the stellar mass and radius, by fitting to the transit curve simultaneously with an empirical mass-radius relation for M-dwarfs (M=R). We find Rs=Ms=0.47 +/-0.02 in solar units, and Rp=27,600 +/-1170 km (4.33 +/-0.18 Earth radii). This radius significantly exceeds the radius of a naked ocean planet, and requires a gasesous hydrogen-helium envelope. The secondary eclipse occurs at phase 0.587 +/-0.005, proving a significant orbital eccentricity (e=0.15 +/-0.012). The amplitude of the eclipse (5.7 +/-0.8e-4) indicates a brightness temperature for the planet of T=712 +/-36K. If this is indicative of the planet's physical temperature, it suggests the occurrence of tidal heating in the planet. An uncharacterized second planet likely provides ongoing gravitational perturbations, to maintain GJ436b's orbit eccentricity over long time scales.