We demonstrate a path to hitherto unachievable differential photometric precisions from the ground, both in the optical and near-infrared (NIR), using custom-fabricated beam-shaping diffusers produced using specialized nanofabrication techniques. Such diffusers mold the focal plane image of a star into a broad and stable top-hat shape, minimizing photometric errors due to non-uniform pixel response, atmospheric seeing effects, imperfect guiding, and telescope-induced variable aberrations seen in defocusing. This PSF reshaping significantly increases the achievable dynamic range of our observations, increasing our observing efficiency and thus better averages over scintillation. Diffusers work in both collimated and converging beams. We present diffuser-assisted optical observations demonstrating $62^{+26}_{-16}$ppm precision in 30 minute bins on a nearby bright star 16-Cygni A (V=5.95) using the ARC 3.5m telescope---within a factor of $\sim$2 of Kepler's photometric precision on the same star. We also show a transit of WASP-85-Ab (V=11.2) and TRES-3b (V=12.4), where the residuals bin down to $180^{+66}_{-41}$ppm in 30 minute bins for WASP-85-Ab---a factor of $\sim$4 of the precision achieved by the K2 mission on this target---and to 101ppm for TRES-3b. In the NIR, where diffusers may provide even more significant improvements over the current state of the art, our preliminary tests have demonstrated $137^{+64}_{-36}$ppm precision for a $K_S =10.8$ star on the 200" Hale Telescope. These photometric precisions match or surpass the expected photometric precisions of TESS for the same magnitude range. This technology is inexpensive, scalable, easily adaptable, and can have an important and immediate impact on the observations of transits and secondary eclipses of exoplanets.

Saturn's moon Enceladus offers a unique opportunity in the search for life and habitable environments beyond Earth, a key theme of the National Research Council's 2013-2022 Decadal Survey. A plume of water vapor and ice spews from Enceladus's south polar region. Cassini data suggest that this plume, sourced by a liquid reservoir beneath the moon's icy crust, contain organics, salts, and water-rock interaction derivatives. Thus, the ingredients for life as we know it-- liquid water, chemistry, and energy sources-- are available in Enceladus's subsurface ocean. We have only to sample the plumes to investigate this hidden ocean environment. We present a New Frontiers class, solar-powered Enceladus orbiter that would take advantage of this opportunity, Testing the Habitability of Enceladus's Ocean (THEO). Developed by the 2015 Jet Propulsion Laboratory Planetary Science Summer School student participants under the guidance of TeamX, this mission concept includes remote sensing and in situ analyses with a mass spectrometer, a sub-mm radiometer-spectrometer, a camera, and two magnetometers. These instruments were selected to address four key questions for ascertaining the habitability of Enceladus's ocean within the context of the moon's geological activity: (1) How are the plumes and ocean connected? (2) Are the abiotic conditions of the ocean suitable for habitability? (3) How stable is the ocean environment? (4) Is there evidence of biological processes? By taking advantage of the opportunity Enceladus's plumes offer, THEO represents a viable, solar-powered option for exploring a potentially habitable ocean world of the outer solar system.

We report secondary eclipse photometry of the hot Jupiter HAT-P-32Ab, taken with Hale/WIRC in H and Ks bands and with Spitzer/IRAC at 3.6 and 4.5 micron. We carried out adaptive optics imaging of the planet host star HAT-P-32A and its companion HAT-P-32B in the near-IR and the visible. We clearly resolve the two stars from each other and find a separation of 2.923" +/- 0. 004" and a position angle 110.64 deg +/- 0.12 deg. We measure the flux ratios of the binary in g' r' i' z' and H & Ks bands, and determine Teff = 3565 +/- 82 K for the companion star, corresponding to an M1.5 dwarf. We use PHOENIX stellar atmosphere models to correct the dilution of the secondary eclipse depths of the hot Jupiter due to the presence of the M1.5 companion. We also improve the secondary eclipse photometry by accounting for the non-classical, flux-dependent nonlinearity of the WIRC IR detector in the H band. We measure planet-to-star flux ratios of 0.090 +/- 0.033%, 0.178 +/- 0.057%, 0.364 +/- 0.016%, and 0.438 +/- 0.020% in the H, Ks, 3.6 and 4.5 micron bands, respectively. We compare these with planetary atmospheric models, and find they prefer an atmosphere with a temperature inversion and inefficient heat redistribution. However, we also find that the data are equally well-described by a blackbody model for the planet with Tp = 2042 +/- 50 K. Finally, we measure a secondary eclipse timing offset of 0.3 +/- 1.3 min from the predicted mid-eclipse time, which constrains e = 0.0072 +0.0700/-0.0064 when combined with RV data and is more consistent with a circular orbit.

Kepler-13Ab (= KOI-13.01) is a unique transiting hot Jupiter. It is one of very few known short-period planets orbiting a hot A-type star, making it one of the hottest planets currently known. The availability of Kepler data allows us to measure the planet's occultation (secondary eclipse) and phase curve in the optical, which we combine with occultations observed by warm Spitzer at 4.5 mic and 3.6 mic and a ground-based occultation observation in the Ks band (2.1 mic). We derive a day-side hemisphere temperature of 2,750 +- 160 K as the effective temperature of a black body showing the same occultation depths. Comparing the occultation depths with one-dimensional planetary atmosphere models suggests the presence of an atmospheric temperature inversion. Our analysis shows evidence for a relatively high geometric albedo, Ag= 0.33 +0.04 -0.06. While measured with a simplistic method, a high Ag is supported also by the fact that the one-dimensional atmosphere models underestimate the occultation depth in the optical. We use stellar spectra to determine the dilution, in the four wide bands where occultation was measured, due to the visual stellar binary companion 1.15 +- 0.05" away. The revised stellar parameters measured using these spectra are combined with other measurements leading to revised planetary mass and radius estimates of Mp = 4.94 - 8.09 Mjup and Rp = 1.406 +- 0.038 Rjup. Finally, we measure a Kepler mid-occultation time that is 34.0 +- 6.9 s earlier than expected based on the mid-transit time and the delay due to light travel time, and discuss possible scenarios.

We report secondary eclipse photometry of two hot Jupiters, WASP-48b and HAT-P-23b, at 3.6 and 4.5 um taken with the InfraRed Array Camera aboard the Spitzer Space Telescope during the warm Spitzer mission and in the H and Ks bands with the Wide Field IR Camera at the Palomar 200-inch Hale Telescope. WASP-48b and HAT-P-23b are Jupiter-mass and twice Jupiter-mass objects orbiting an old, slightly evolved F star and an early G dwarf star, respectively. In the H, Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.047% +/- 0.016%, 0.109% +/- 0.027%, 0.176% +/- 0.013%, and 0.214% +/- 0.020% for WASP-48b. In the Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.234% +/- 0.046%, 0.248% +/- 0.019%, and 0.309% +/- 0.026% for HAT-P-23b. For WASP-48b and HAT-P-23b, respectively, we measure delays of 2.6 +/- 3.9 minutes and 4.0 +/- 2.4 minutes relative to the predicted times of secondary eclipse for circular orbits, placing 2-sigma upper limits on |e cos(omega)| of 0.0053 and 0.0080, both of which are consistent with circular orbits. The dayside emission spectra of these planets are well-described by blackbodies with effective temperatures of 2158 +/- 100 K (WASP-48b) and 2154 +/- 90 K (HAT-P-23b), corresponding to moderate recirculation in the zero albedo case. Our measured eclipse depths are also consistent with one-dimensional radiative transfer models featuring varying degrees of recirculation and weak thermal inversions or no inversions at all. We discuss how the absence of strong temperature inversions on these planets may be related to the activity levels and metallicities of their host stars.

Titan's moment of inertia, calculated assuming hydrostatic equilibrium from gravity field data obtained during the Cassini-Huygens mission, implies an internal mass distribution that may be incompatible with complete differentiation. This suggests that Titan may have a mixed ice/rock core, possibly consistent with slow accretion in a gas-starved disk, which may initially spare Titan from widespread ice melting and subsequent differentiation. A partially differentiated Titan, however, must still efficiently remove radiogenic heat over geologic time. We argue that compositional heterogeneity in the major Saturnian satellites indicates that Titan formed from planetesimals with disparate densities. The resulting compositional anomalies would quickly redistribute to form a vertical density gradient that would oppose thermal convection. We use elements of the theory of double-diffusive convection to create a parameterized model for the thermal evolution of ice/rock mixtures with a stabilizing compositional gradient. Simulations are performed for a wide range of initial conditions to account for large uncertainties in material properties and accretionary processes. Ultimately, for realistic density gradients, double-diffusive convection in the ice/rock interior can delay, but not prevent, ice melting and differentiation, even if a substantial fraction of potassium is leached from the rock component. Consequently, Titan is not partially differentiated.

The ongoing discovery of terrestrial exoplanets accentuates the importance of studying planetary evolution for a wide range of initial conditions. We perform thermal evolution simulations for generic terrestrial planets with masses ranging from that of Mars to 10 Earth-masses in the stagnant-lid regime, the most natural mode of convection with strongly temperature- dependent viscosity. Given considerable uncertainty surrounding the dependency of mantle rheology on pressure, we choose to focus on the end-member case of pressure-independent potential viscosity, where viscosity does not change with depth along an adiabatic temperature gradient. We employ principal component analysis and linear regression to capture the first-order systematics of possible evolutionary scenarios from a large number of simulation runs. With increased planetary mass, crustal thickness and the degree of mantle processing are both predicted to decrease, and such size effects can also be derived with simple scaling analyses. The likelihood of plate tectonics is quantified using a mantle rheology that takes into account both ductile and brittle deformation mechanisms. Confirming earlier scaling analyses, the effects of lithosphere hydration dominate the effects of planetary mass. The possibility of basalt-eclogite phase transition in the planetary crust is found to increase with planetary mass, and we suggest that massive terrestrial planets may escape the stagnant-lid regime through the formation of a self-destabilizing dense eclogite layer.