We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

High contrast direct imaging of exoplanets can provide many important observables, including measurements of the orbit, spectra that probe the lower layers of the atmosphere, and phase variations of the planet, but cannot directly measure planet radius or mass. Our future understanding of directly imaged exoplanets will therefore rely on extrapolated models of planetary atmospheres and bulk composition, which need robust calibration. We estimate the population of extrasolar planets that could serve as calibrators for these models. Critically, this population of "standard planets" must be accessible to both direct imaging and the transit method, allowing for radius measurement. We show that the search volume of a direct imaging mission eventually overcomes the transit probability falloff with semi-major axis, so that as long as cold planets are not exceedingly rare, the population of transiting planets and directly imageable planets overlaps. Using current extrapolations of Kepler occurrence rates, we estimate that ~8 standard planets could be characterized shortward of 800 nm with an ambitious future direct imaging mission like LUVOIR-A and several dozen could be detected at V band. We show the design space that would expand the sample size and discuss the extent to which ground- and space-based surveys could detect this small but crucial population of planets.

Very little experimental work has been done to explore the properties of photochemical hazes formed in atmospheres with very different compositions or temperatures than that of the outer solar system or of early Earth. With extrasolar planet discoveries now numbering thousands, this untapped phase space merits exploration. This study presents measured chemical properties of haze particles produced in laboratory analogues of exoplanet atmospheres. We used very high resolution mass spectrometry to measure the chemical components of solid particles produced in atmospheric chamber experiments. Many complex molecular species with general chemical formulas C$_w$H$_x$N$_y$O$_z$ were detected. We detect molecular formulas of prebiotic interest in the data, including those for the monosaccharide glyceraldehyde, a variety of amino acids and nucleotide bases, and several sugar derivatives. Additionally, the experimental exoplanetary haze analogues exhibit diverse solubility characteristics, which provide insight into the possibility of further chemical or physical alteration of photochemical hazes in super-Earth and mini-Neptune atmospheres. These exoplanet analogue particles can help us better understand chemical atmospheric processes and suggest a possible source of in situ atmospheric prebiotic chemistry on distant worlds.

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's $M$-band spectrum by including both CO and water clouds in the atmospheric model.

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry.

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

Sulfur gases significantly affect the photochemistry of planetary atmospheres in our Solar System, and are expected to be important components in exoplanet atmospheres. However, sulfur photochemistry in the context of exoplanets is poorly understood due to a lack of chemical-kinetics information for sulfur species under relevant conditions. Here, we study the photochemical role of hydrogen sulfide (H2S) in warm CO2-rich exoplanet atmospheres (800 K) by carrying out laboratory simulations. We find that H2S plays a significant role in photochemistry, even when present in the atmosphere at relatively low concentrations (1.6%). It participates in both gas and solid phase chemistry, leading to the formation of other sulfur gas products (CH3SH/SO, C2H4S/OCS, SO2/S2, and CS2) and to an increase in solid haze particle production and compositional complexity. Our study shows that we may expect thicker haze with small particle sizes (20 to 140 nm) for warm CO2-rich exoplanet atmospheres that possess H2S.

We present the first discovery from the COol Companions ON Ultrawide orbiTS (COCONUTS) program, a large-scale survey for wide-orbit planetary and substellar companions. We have discovered a co-moving system COCONUTS-1, composed of a hydrogen-dominated white dwarf (PSO J058.9855+45.4184; $d=31.5$ pc) and a T4 companion (PSO J058.9869+45.4296) at a $40.6''$ (1280 au) projected separation. We derive physical properties for COCONUTS-1B from (1) its near-infrared spectrum using cloudless Sonora atmospheric models, and (2) its luminosity and the white dwarf's age ($7.3_{-1.6}^{+2.8}$ Gyr) using Sonora evolutionary models. The two methods give consistent temperatures and radii, but atmospheric models infer a lower surface gravity and therefore an unphysically young age. Assuming evolutionary model parameters ($T_{\rm eff}=1255^{+6}_{-8}$ K, $\log{g}=5.44^{+0.02}_{-0.03}$ dex, $R=0.789^{+0.011}_{-0.005}$ R$_{\rm Jup}$), we find cloudless model atmospheres have brighter Y- and J-band fluxes than the data, suggesting condensate clouds have not fully dispersed around 1300 K. The W2 flux (4.6 $\mu$m) of COCONUTS-1B is fainter than models, suggesting non-equilibrium mixing of CO. To investigate the gravity dependence of the L/T transition, we compile all 60 known L6-T6 benchmarks and derive a homogeneous set of temperatures, surface gravities, and masses. As is well-known, young, low-gravity late-L dwarfs have significantly fainter, redder near-infrared photometry and $\approx200-300$ K cooler temperatures than old, high-gravity objects. Our sample now reveals such gravity dependence becomes weaker for T dwarfs, with young objects having comparable near-infrared photometry and $\approx100$ K cooler temperatures compared to old objects. Finally, we find that young objects have a larger amplitude J-band brightening than old objects, and also brighten at H band as they cross the L/T transition.

We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

Among the greatest challenges in understanding ultra-cool brown dwarf and exoplanet atmospheres is the evolution of cloud structure as a function of temperature and gravity. In this study, we present the rotational modulations of GU Psc b -- a rare mid-T spectral type planetary-mass companion at the end of the L/T spectral type transition. Based on the HST/WFC3 1.1-1.67$\rm\, \mu m$ time-series spectra, we observe a quasi-sinusoidal light curve with a peak-to-trough flux variation of 2.7 % and a minimum period of eight hours. The rotation-modulated spectral variations are weakly wavelength-dependent, or largely gray between 1.1-1.67$\rm\,\mu$m. The gray modulations indicate that heterogeneous clouds are present in the photosphere of this low-gravity mid-T dwarf. We place the color and brightness variations of GU Psc b in the context of rotational modulations reported for mid-L to late-T dwarfs. Based on these observations, we report a tentative trend: mid-to-late T dwarfs become slightly redder in $J-H$ color with increasing $J$-band brightness, while L dwarfs become slightly bluer with increasing brightness. If this trend is verified with more T-dwarf samples, it suggests that in addition to the mostly gray modulations, there is a second-order spectral-type dependence on the nature of rotational modulations.

We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4.

We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems.

The characterization of rocky, Earth-like planets is an important goal for future large ground- and space-based telescopes. In support of developing an efficient observational strategy, we have applied Bayesian statistical inference to interpret the albedo spectrum of cloudy true-Earth analogs that include a diverse spread in their atmospheric water vapor mixing ratios. We focus on detecting water-bearing worlds by characterizing their atmospheric water vapor content via the strong 0.94$\,\mu$m H$_2$O absorption feature, with several observational configurations. Water vapor is an essential signpost when assessing planetary habitability, and determining its presence is important in vetting whether planets are suitable for hosting life. We find that R=140 spectroscopy of the absorption feature combined with a same-phase green optical photometric point at $0.525-0.575\,\mu$m is capable of distinguishing worlds with less than $0.1\times$ Earth-like water vapor levels from worlds with $1\times$ Earth-like levels or greater at a signal-to-noise ratio of 5 or better with $2\sigma$ confidence. This configuration can differentiate between $0.01\times$ and $0.1\times$ Earth-like levels when the signal-to-noise ratio is 10 or better at the same confidence. However, strong constraints on the water vapor mixing ratio remained elusive with this configuration even at signal-to-noise of 15. We find that adding the same-phase optical photometric point does not significantly help characterize the H$_2$O mixing ratio, but does enable an upper limit on atmospheric ozone levels. Finally, we find that a 0.94$\,\mu$m photometric point, instead of spectroscopy, combined with the green-optical point, fails to produce meaningful information about atmospheric water content.

Future space telescopes now in the concept and design stage aim to observe reflected light spectra of extrasolar planets. Assessing whether given notional mission and instrument design parameters will provide data suitable for constraining quantities of interest typically requires time consuming retrieval studies in which tens to hundreds of thousands of models are compared to data with a given assumed signal to noise ratio, thereby limiting the rapidity of design iterations. Here we present a novel machine learning approach employing five Multilayer Perceptron's (MLP) trained on model albedo spectra of extrasolar giant planets to estimate a planet's atmospheric metallicity, gravity, effective temperature, and cloud properties given simulated observed spectra. The stand-alone C++ code we have developed can train new MLP's on new training sets, within minutes to hours, depending upon the dimensions of input spectra and desired output. After the MLP is trained, it can classify new input spectra within a second, potentially helping speed observation and mission design planning. Four of the MLP's were tuned to work with various levels of spectral noise. The fifth MLP was developed for cases where the user is uncertain about the noise level. The MLP's were trained using a grid of model spectra that varied in metallicity, gravity, temperature, and cloud properties. We tested the MLP on noisy models and observed spectra of both Jupiter and Saturn. The root mean squared error when applied to noisy models were on the order of the model grid intervals. The results show that a trained MLP is a robust means for reliable in-situ estimations, providing an elegant artificial intelligence that is simple to customize and quick to use.

As part of the PanCET program, we have conducted a spectroscopic study of WASP-79b, an inflated hot Jupiter orbiting an F-type star in Eridanus with a period of 3.66 days. Building on the original WASP and TRAPPIST photometry of Smalley et al (2012), we examine HST/WFC3 (1.125 - 1.650 $\mu$m), Magellan/LDSS-3C (0.6 - 1 $\mu$m) data, and Spitzer data (3.6 and 4.5 $\mu$m). Using data from all three instruments, we constrain the water abundance to be --2.20 $\leq$ log(H$_2$O) $\leq$ --1.55. We present these results along with the results of an atmospheric retrieval analysis, which favor inclusion of FeH and H$^-$ in the atmospheric model. We also provide an updated ephemeris based on the Smalley, HST/WFC3, LDSS-3C, Spitzer, and TESS transit times. With the detectable water feature and its occupation of the clear/cloudy transition region of the temperature/gravity phase space, WASP-79b is a target of interest for the approved JWST Director's Discretionary Early Release Science (DD ERS) program, with ERS observations planned to be the first to execute in Cycle 1. Transiting exoplanets have been approved for 78.1 hours of data collection, and with the delay in the JWST launch, WASP-79b is now a target for the Panchromatic Transmission program. This program will observe WASP-79b for 42 hours in 4 different instrument modes, providing substantially more data by which to investigate this hot Jupiter.

The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.

We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$\sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).

The relative abundance of deuterium and hydrogen is a potent tracer of planet formation and evolution. Jupiter and Saturn have protosolar atmospheric D/H ratios, a relic of substantial gas accretion from the nebula, while the atmospheres of Neptune and Uranus are enhanced in D by accretion of ices into their envelopes. For terrestrial planets, D/H ratios are used to determine the mechanisms of volatile delivery and subsequent atmosphere loss over the lifetime of the planet. Planets and brown dwarfs more massive than 13 M_J quickly fuse their initial D reservoir. Here, we simulate spectra for giant exoplanets and brown dwarfs (2 M_Neptune to ~10 M_Jupiter) from Teff=200--1800 K including both CH3D and HDO to determine the observability of these dominant deuterium isotopologues in mid-infrared thermal emission spectra. Colder objects have stronger molecular features in their spectra, due to the temperature-dependence of molecular cross sections. CH3D is easier to observe than HDO at all temperatures considered, due to the strength of its absorption bands and locations of features at wavelengths with few other strong absorption features. We predict that for nearby cool brown dwarfs, the CH3D abundance can be measured at high signal-to-noise with the James Webb Space Telescope; for objects from 200--800 K closer than 10 pc, a protosolar D/H ratio would be readily observable in 2.5 hours. Moderately young Jupiter-mass planets (100--300 Myr) and young Neptunes (10 Myr) may be discovered with JWST and provide the best targets for detecting deuterium on an exoplanet in the coming decade. Future telescope designs should consider the importance of isotopes for understanding the formation and evolution of planetary atmospheres.

We present constraints on the presence of substellar companions to the nearby ($d\sim7.6$ pc), young ($440\pm40$ Myr) K4Ve star TW Piscis Austrini, the wide ($\sim$0.3 pc) companion to the A4V star Fomalhaut. We combined absolute astrometry from Hipparcos and Gaia with literature radial velocity measurements and dedicated high-contrast imaging observations obtained with Keck/NIRC2 to achieve sensitivity to brown dwarf and planetary-mass companions ($\gtrsim2$ M$_{\rm Jup}$) over many decades of orbital period ($\lesssim 10^3$ yr). The significant astrometric acceleration measured between the Hipparcos and Gaia catalogues, reported previously in the literature, cannot be explained by the orbital motion of TW PsA around the barycenter of the Fomalhaut triple system. Instead, we find that it is consistent with the reflex motion induced by an orbiting substellar companion. The combination of astrometry, radial velocities, and a deep $L^{\prime}$ imaging dataset leads to a constraint on the companion mass of $1.2_{-0.6}^{+0.7}$ M$_{\rm Jup}$. However, the period of the companion is poorly constrained, with a highly multi-modal period posterior distribution due to aliasing with the 24.25 year baseline between Hipparcos and Gaia. If confirmed through continued astrometric or spectroscopic monitoring, or via direct detection, the companion to TW PsA would represent a choice target for detailed atmospheric characterization with high-contrast instruments on the upcoming James Webb Space Telescope and Wide Field Infrared Survey Telescope.

An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay, and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low towards the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations.

Probabilistic relationships between the mass and radius of planets are useful for a variety of purposes, including estimating the yields of planet discovery efforts and the radii of discovered planets given their masses. Previous work on giant planets often do not make the important distinction between cool and (inflated) hot giants. This can lead to needlessly large uncertainties and overestimated radii for the cooler objects. Using hierarchical bayesian modeling and the AIC model selection criterion, we identify an empirical mass-radius relationship for cool giant planets above 15 Earth masses. It takes the form of a quadratic function in log mass, with a power-law dispersion. Parametric uncertainties from MCMC are provided, and we compare our results with the outputs of interior structure models and the radii of the solar system giant planets, which fit comfortably within our relation.

We use the Wide Field Camera 3 on the {\sl Hubble Space Telescope} to spectrophotometrically monitor the young L7.5 companion HD~203030B. Our time series reveal photometric variability at 1.27\,$\mu$m and 1.39\,$\mu$m on time scales compatible with rotation. We find a rotation period of $7.5^{+0.6}_{-0.5}$ h: comparable to those observed in other brown dwarfs and planetary-mass companions younger than 300 Myr. We measure variability amplitudes of $1.1\pm0.3\%$ (1.27\,$\mu$m) and $1.7\pm0.4\%$ (1.39\,$\mu$m), and a phase lag of 56$^\circ\pm$28$^\circ$ between the two light curves. We attribute the difference in photometric amplitudes and phases to a patchy cloud layer that is sinking below the level where water vapor becomes opaque. HD 203030B and the few other known variable young late-L dwarfs are unlike warmer (earlier-type and/or older) L dwarfs, for which variability is much less wavelength-dependent across the 1.1--1.7$\mu$m region. We further suggest that a sinking of the top-most cloud deck below the level where water or carbon monoxide gas become opaque may also explain the often enhanced variability amplitudes of even earlier-type low-gravity L dwarfs. Because these condensate and gas opacity levels are already well-differentiated in T dwarfs, we do not expect the same variability amplitude enhancement in young vs.\ old T dwarfs.

Half of the energy emitted by late-T- and Y-type brown dwarfs emerges at 3.5 < lambda um < 5.5. We present new L' (3.43 < lambda um < 4.11) photometry obtained at the Gemini North telescope for nine late-T and Y dwarfs, and synthesize L' from spectra for an additional two dwarfs. The targets include two binary systems which were imaged at a resolution of 0.25". One of these, WISEP J045853.90+643452.6AB, shows significant motion, and we present an astrometric analysis of the binary using Hubble Space Telescope, Keck Adaptive Optics, and Gemini images. We compare lambda ~4um observations to models, and find that the model fluxes are too low for brown dwarfs cooler than ~700K. The discrepancy increases with decreasing temperature, and is a factor of ~2 at T_eff=500K and ~4 at T_eff=400K. Warming the upper layers of a model atmosphere generates a spectrum closer to what is observed. The thermal structure of cool brown dwarf atmospheres above the radiative-convective boundary may not be adequately modelled using pure radiative equilibrium; instead heat may be introduced by thermochemical instabilities (previously suggested for the L- to T-type transition) or by breaking gravity waves (previously suggested for the solar system giant planets). One-dimensional models may not capture these atmospheres, which likely have both horizontal and vertical pressure/temperature variations.

We observed two full orbital phase curves of the transiting brown dwarf KELT-1b, at 3.6um and 4.5um, using the Spitzer Space Telescope. Combined with previous eclipse data from Beatty et al. (2014), we strongly detect KELT-1b's phase variation as a single sinusoid in both bands, with amplitudes of $964\pm36$ ppm at 3.6um and $979\pm54$ ppm at 4.5um, and confirm the secondary eclipse depths measured by Beatty et al. (2014). We also measure noticeable Eastward hotspot offsets of $28.4\pm3.5$ degrees at 3.6um and $18.6\pm5.2$ degrees at 4.5um. Both the day-night temperature contrasts and the hotspot offsets we measure are in line with the trends seen in hot Jupiters (e.g., Crossfield 2015), though we disagree with the recent suggestion of an offset trend by Zhang et al. (2018). Using an ensemble analysis of Spitzer phase curves, we argue that nightside clouds are playing a noticeable role in modulating the thermal emission from these objects, based on: 1) the lack of a clear trend in phase offsets with equilibrium temperature, 2) the sharp day-night transitions required to have non-negative intensity maps, which also resolves the inversion issues raised by Keating & Cowan (2017), 3) the fact that all the nightsides of these objects appear to be at roughly the same temperature of 1000K, while the dayside temperatures increase linearly with equilibrium temperature, and 4) the trajectories of these objects on a Spitzer color-magnitude diagram, which suggest colors only explainable via nightside clouds.

We present HST near-ultraviolet (NUV) transits of the hot Jupiter WASP-121b, acquired as part of the PanCET program. Time series spectra during two transit events were used to measure the transmission spectra between 2280 and 3070 Angstroms at a resolution of 30,000. Using HST data from 61 STIS visits, we show that data from HST's Pointing Control System can be used to decorrelate the instrument systematic errors (Jitter Decorrelation), which we used to fit the WASP-121b light curves. The NUV spectrum show very strong absorption features, with the NUV white light curve found to be larger than the average optical and near-infrared value at 6-$\sigma$ confidence. We identify and spectrally resolve absorption from the Mg ii doublet in the planetary exosphere at a 5.9-$\sigma$ confidence level. The Mg ii doublet is observed to reach altitudes of $R_{pl}/R_{star}=0.284\pm0.037$ for the 2796 Angstrom line and $0.242\pm0.0431$ in the 2804 Angstrom line, which exceeds the Roche lobe size as viewed in transit geometry ($R_{\rm eqRL}/R_{star}$ = 0.158). We also detect and resolve strong features of the Fe ii UV1 and UV2 multiplets, and observe the lines reaching altitudes of $R_{pl}/R_{star}\approx0.3$. At these high altitudes, the atmospheric Mg ii and Fe ii gas is not gravitationally bound to the planet, and these ionized species may be hydrodynamically escaping or could be magnetically confined. Refractory Mg and Fe atoms at high altitudes also indicates that these species are not trapped into condensate clouds at depth, which places constraints on the deep interior temperature.

WASP-121b is a transiting gas giant exoplanet orbiting close to its Roche limit, with an inflated radius nearly double that of Jupiter and a dayside temperature comparable to a late M dwarf photosphere. Secondary eclipse observations covering the 1.1-1.6 micron wavelength range have revealed an atmospheric thermal inversion on the dayside hemisphere, likely caused by high altitude absorption at optical wavelengths. Here we present secondary eclipse observations made with the Hubble Space Telescope Wide Field Camera 3 spectrograph that extend the wavelength coverage from 1.1 micron down to 0.8 micron. To determine the atmospheric properties from the measured eclipse spectrum, we performed a retrieval analysis assuming chemical equilibrium, with the effects of thermal dissociation and ionization included. Our best-fit model provides a good fit to the data with reduced chi^2=1.04. The data diverge from a blackbody spectrum and instead exhibit emission due to H- shortward of 1.1 micron. The best-fit model does not reproduce a previously reported bump in the spectrum at 1.25 micron, possibly indicating this feature is a statistical fluctuation in the data rather than a VO emission band as had been tentatively suggested. We estimate an atmospheric metallicity of [M/H]=1.09(-0.69,+0.57), and fit for the carbon and oxygen abundances separately, obtaining [C/H]=-0.29(-0.48,+0.61) and [O/H]=0.18(-0.60,+0.64). The corresponding carbon-to-oxygen ratio is C/O=0.49(-0.37,+0.65), which encompasses the solar value of 0.54, but has a large uncertainty.

The reflected light signal from a planet throughout its orbit is a powerful probe of a planet's atmospheric properties. There are a number of planets that are amenable to reflected light phase curve studies with present and future space-based instrumentation and here we assess our ability to characterize these worlds. Using simulated TESS populations we identify the Nine, a set of archetypal exoplanets with the potential to be bright in reflected light, because of their radii and proximity to their star, while still being cool enough to have minimal thermal contamination at optical wavelengths. For each planet we compute albedo spectra for several cloud and atmosphere assumptions (e.g. thermochemical equilibrium, solar composition). We find that in the TESS bandpass the estimated contrast at optical wavelengths is typically <10 ppm except for the brightest, largest, or closest in planets with the highest lofted clouds where contrast can reach a few tens of ppm. Meanwhile, in a bluer bandpass (0.3--0.5 microns) the estimated contrast can be as high as 150 ppm but typically 10-50 ppm. In the temperature range of interest, planets with the highest, most extensive cloud decks are generally darker at bluer wavelengths than cloudless planets because of the low single scattering albedos of their primary condensate constituents. Our models suggest that Neptune-sized planets with relatively low insolation and small semi-major axes are the most conducive to reflected light phase curve studies in TESS.

Here we present the first open-source radiative transfer model for computing the reflected light of exoplanets at any phase geometry, called PICASO: Planetary Intensity Code for Atmospheric Scattering Observations. This code, written in Python, has heritage from a decades old, well-known Fortran model used for several studies of planetary objects within the Solar System and beyond. We have adopted it to include several methodologies for computing both direct and diffuse scattering phase functions, and have added several updates including the ability to compute Raman scattering spectral features. Here we benchmark PICASO against two independent codes and discuss the degree to which the model is sensitive to a user's specification for various phase functions. Then, we conduct a full information content study of the model across a wide parameter space in temperature, cloud profile, SNR and resolving power.

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$\sigma$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\sim\lambda/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $\sigma\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Spectroscopic transiting observations of the atmospheres of hot Jupiters around other stars, first with Hubble Space Telescope and then Spitzer, opened the door to compositional studies of exoplanets. The James Webb Space Telescope will provide such a profound improvement in signal-to-noise ratio that it will enable detailed analysis of molecular abundances, including but not limited to determining abundances of all the major carbon- and oxygen-bearing species in hot Jupiter atmospheres. This will allow determination of the carbon-to-oxygen ratio, an essential number for planet formation models and a motivating goal of the Juno mission currently around Jupiter

Measurements of photometric variability at different wavelengths provide insights into the vertical cloud structure of brown dwarfs and planetary-mass objects. In seven Hubble Space Telescope consecutive orbits, spanning $\sim$10 h of observing time}, we obtained time-resolved spectroscopy of the planetary-mass T8-dwarf Ross 458C using the near-infrared Wide Field Camera 3. We found spectrophotometric variability with a peak-to-peak signal of 2.62$\pm$0.02 % (in the 1.10-1.60~$\mu$m white light curve). Using three different methods, we estimated a rotational period of 6.75$\pm$1.58~h for the white light curve, and similar periods for narrow $J$- and $H$- band light curves. Sine wave fits to the narrow $J$- and $H$-band light curves suggest a tentative phase shift between the light curves with wavelength when we allow different periods between both light curves. If confirmed, this phase shift may be similar to the phase shift detected earlier for the T6.5 spectral type 2MASS J22282889-310262. We find that, in contrast with 2M2228, the variability of Ross~458C shows evidence for a {color trend} within the narrow $J$-band, but gray variations in the narrow $H$-band. The spectral time-resolved variability of Ross 458C might be potentially due to heterogeneous sulfide clouds in the atmosphere of the object. Our discovery extends the study of spectral modulations of condensate clouds to the coolest T dwarfs, planetary-mass companions.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

There are many competing theories and models describing the formation, migration and evolution of exoplanet systems. As both the precision with which we can characterize exoplanets and their host stars, and the number of systems for which we can make such a characterization increase, we begin to see pathways forward for validating these theories. In this white paper we identify predicted, observable correlations that are accessible in the near future, particularly trends in exoplanet populations, radii, orbits and atmospheres with host star age. By compiling a statistically significant sample of well-characterized exoplanets with precisely measured ages, we should be able to begin identifying the dominant processes governing the time-evolution of exoplanet systems.

Gas giants are the most readily detectable exoplanets but fundamental questions about their system architectures, formation, migration, and atmospheres have been unanswerable with the current generation of ground- and space-based facilities. The dominant techniques to detect and characterize giant planets $-$ radial velocities, transits, direct imaging, microlensing, and astrometry $-$ are each isolated to a limited range of planet masses, separations, ages, and temperatures. These windows into the arrangement and physical properties of giant planets have spawned new questions about the timescale and location of their assembly; the distributions of planet mass and orbital separation at young and old ages; the composition and structure of their atmospheres; and their orbital and rotational angular momentum architectures. The ELTs will address these questions by building bridges between these islands of mass, orbital distance, and age. The angular resolution, collecting area, all-sky coverage, and novel instrumentation suite of these facilities are needed to provide a complete map of the orbits and atmospheric evolution of gas giant planets (0.3$-$10 $M_\mathrm{Jup}$) across space (0.1$-$100 AU) and time (1 Myr to 10 Gyr). This white paper highlights the scientific potential of the GMT and TMT to address these outstanding questions, with a particular focus on the role of direct imaging and spectroscopy of large samples of giant planets that will soon be made available with $Gaia$.

This White Paper describes the opportunities for discovery of Jupiter-mass objects with 300K atmospheres. The discovery and characterization of such cold objects is vital for understanding the low-mass terminus of the initial mass function and for optimizing the study of exoplanets by the next generation of large telescopes, space probes and space missions.

Ground-based telescopes coupled with adaptive optics (AO) have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Over the next decade, extreme AO systems on 8-10m telescopes will 1) mitigate risk for WFIRST-CGI by identifying numerous planets the mission can spectrally characterize, 2) validate performance requirements and motivate improvements to atmosphere models needed to unambiguously characterize solar system-analogues from space, and 3) mature novel technological innovations useful for space. Extremely Large Telescopes can deliver the first thermal infrared (10 $\mu m$) images of rocky planets around Sun-like stars and identify biomarkers. These data provide a future NASA direct imaging flagship mission (i.e. HabEx, LUVOIR) with numerous exo-Earth candidates and critical ancillary information to help clarify whether these planets are habitable.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

Directly-imaged planetary-mass companions offer unique opportunities in atmospheric studies of exoplanets. They share characteristics of both brown dwarfs and transiting exoplanets, therefore, are critical for connecting atmospheric characterizations for these objects. Rotational phase mapping is a powerful technique to constrain the condensate cloud properties in ultra-cool atmospheres. Applying this technique to directly-imaged planetary-mass companions will be extremely valuable for constraining cloud models in low mass and surface gravity atmospheres and for determining the rotation rate and angular momentum of substellar companions. Here, we present Hubble Space Telescope Wide Field Camera 3 near-infrared time-resolved photometry for three planetary-mass companions, AB Pic B, 2M0122B, and 2M1207b. Using two-roll differential imaging and hybrid point spread function modeling, we achieve sub-percent photometric precision for all three observations. We find tentative modulations ($<\!\!2\sigma$) for AB Pic B and 2M0122B but cannot reach conclusive results on 2M1207b due to strong systematics. The relatively low significance of the modulation measurements cannot rule out the hypothesis that these planetary-mass companions have the same vertical cloud structures as brown dwarfs. Our rotation rate measurements, combined with archival period measurements of planetary-mass companions and brown dwarfs do not support a universal mass-rotation relation. The high precision of our observations and the high occurrence rates of variable low-surface gravity objects encourage high-contrast time-resolved observations with the James Webb Space Telescope.

Photochemistry induced by stellar UV flux should produce haze particles in exoplanet atmospheres. Recent observations indicate that haze and/or cloud layers exist in the atmospheres of exoplanets. However, photochemical processes in exoplanetary atmospheres remain largely unknown. We performed laboratory experiments with the PHAZER chamber to simulate haze formation in a range of exoplanet atmospheres (hydrogen-rich, water-rich, and carbon dioxide-rich at 300, 400, and 600 K), and observed the gas phase compositional change (the destruction of the initial gas and the formation of new gas species) during these experiments with mass spectrometer. The mass spectra reveal that distinct chemical processes happen in the experiments as a function of different initial gas mixture and different energy sources (plasma or UV photons). We find that organic gas products and O2 are photochemically generated in the experiments, demonstrating that photochemical production is one of the abiotic sources for these potential biosignatures. Multiple simulated atmospheres produce organics and O2 simultaneously, which suggests that even the co-presence of organics and O2 could be false positive biosignature. From the gas phase composition changes, we identify potential precursors (C2H2, HCN, CH2NH, HCHO, etc.) for haze formation, among which complex reactions can take place and produce larger molecules. Our laboratory results indicate that complex atmospheric photochemistry can happen in diverse exoplanet atmospheres and lead to the formation of new gas products and haze particles, including compounds (O2 and organics) that could be falsely identified as biosignatures.

Seismology of the gas giants holds the potential to resolve long-standing questions about their internal structure and rotation state. We construct a family of Saturn interior models constrained by the gravity field and compute their adiabatic mode eigenfrequencies and corresponding Lindblad and vertical resonances in Saturn's C ring, where more than twenty waves with pattern speeds faster than the ring mean motion have been detected and characterized using high-resolution Cassini Visual and Infrared Mapping Spectrometer (VIMS) stellar occultation data. We present identifications of the fundamental modes of Saturn that appear to be the origin of these observed ring waves, and use their observed pattern speeds and azimuthal wavenumbers to estimate the bulk rotation period of Saturn's interior to be $10{\rm h}\, 33{\rm m}\, 38{\rm s}^{+1{\rm m}\, 52{\rm s}}_{-1{\rm m}\, 19{\rm s}}$ (median and 5%/95% quantiles), significantly faster than Voyager and Cassini measurements of periods in Saturn's kilometric radiation, the traditional proxy for Saturn's bulk rotation period. The global fit does not exhibit any clear systematics indicating strong differential rotation in Saturn's outer envelope.

GJ 3470b is a warm Neptune transiting an M-dwarf star at the edge of the evaporation desert. It offers the possibility of investigating how low-mass, close-in exoplanets evolve under the irradiation from their host stars. We observed three transits of GJ 3470b in the Lyman-alpha line with the Hubble Space Telescope (HST) as part of the Panchromatic Comparative Exoplanet Treasury (PanCET) program. Absorption signatures are detected with similar properties in all three independent epochs, with absorption depths of 35+-7% in the blue wing of the line, and 23+-5% in the red wing. The repeatability of these signatures, their phasing with the planet transit, and the radial velocity of the absorbing gas allow us to conclude that there is an extended upper atmosphere of neutral hydrogen around GJ 3470 b. We determine from our observations the stellar radiation pressure and XUV irradiation from GJ 3470 and use them to perform numerical simulations of the upper atmosphere of GJ 3470b with the EVaporating Exoplanets (EVE) code. The unusual redshifted signature can be explained by the damping wings of dense layers of neutral hydrogen that extend beyond the Roche lobe and are elongated in the direction of the planet motion. This structure could correspond to a shocked layer of planetary material formed by the collision of the expanding thermosphere with the wind of the star. The blueshifted signature is well explained by neutral hydrogen atoms escaping at rates of about 1e10 g s-1 that are blown away from the star by its strong radiation pressure and are quickly photoionized, resulting in a smaller exosphere than that of the warm Neptune GJ 436b. The stronger escape from GJ 3470b, however, may have led to the loss of about 4-35% of its current mass over its 2 Gyr lifetime.

Bayesian atmospheric retrieval tools can place constraints on the properties of brown dwarfs and hot Jupiters atmospheres. To fully exploit these methods, high signal-to-noise spectral libraries with well-understood uncertainties are essential. We present a high signal-to-noise spectral library (1.10-1.69 microns) of the thermal emission of 76 brown dwarfs and hot Jupiters. All our spectra have been acquired with the Hubble Space Telescope's Wide Field Camera 3 instrument and its G141 grism. The near-infrared spectral types of these objects range from L4 to Y1. Eight of our targets have estimated masses below the deuterium-burning limit. We analyze the database to identify peculiar objects and/or multiple systems, concluding that this sample includes two very-low-surface-gravity objects and five intermediate-surface-gravity objects. In addition, spectral indices designed to search for composite atmosphere brown dwarfs, indicate that eight objects in our sample are strong candidates to have such atmospheres. None of these objects are overluminous, thus their composite atmospheres are unlikely a companion-induced artifact. Five of the eight confirmed candidates have been reported as photometrically variable, suggesting that composite atmospheric indices are useful in identifying brown dwarfs with strongly heterogeneous cloud covers. We compare hot Jupiters and brown dwarfs in a near-infrared color-magnitude diagram. We confirm that the coldest hot Jupiters in our sample have spectra similar to mid-L dwarfs, and the hottest hot Jupiters have spectra similar to those of M-dwarfs. Our sample provides a uniform dataset of a broad range of ultracool atmospheres, allowing large-scale, comparative studies, and providing a HST legacy spectral library.

We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.

We have observed the eclipsing, post-common envelope white dwarf-brown dwarf binary, SDSS141126.20+200911.1, in the near-IR with the HAWK-I imager, and present here the first direct detection of the dark side of an irradiated brown dwarf in the $H$ band, and a tentative detection in the $K_s$ band. Our analysis of the lightcurves and indicates that the brown dwarf is likely to have an effective temperature of 1300 K, which is not consistent with the effective temperature of 800 K suggested by its mass and radius. As the brown dwarf is already absorbing almost all the white dwarf emission in the $K_s$ band we suggest that this inconsistency may be due to the UV-irradiation from the white dwarf inducing an artificial brightening in the $K_s$ band, similar to that seen for the similar system WD0137-349B, suggesting this brightening may be characteristic of these UV-irradiated binaries.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7\mu m$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using the SPARC/MITgcm, we investigate how thermal dissociation, ionization, H$^-$ opacity and clouds shape the thermal structures and spectral properties of ultra hot Jupiters with a special focus on WASP-121b. We expand our findings to the whole population of ultra hot Jupiters through analytical quantification of the thermal dissociation and its influence on the strength of spectral features. We predict that most molecules are thermally dissociated and alkalies are ionized in the dayside photospheres of ultra hot Jupiters. This includes H$_{\rm 2}$O, TiO, VO, and H$_{\rm 2}$ but not CO, which has a stronger molecular bond. The vertical molecular gradient created by the dissociation significantly weakens the spectral features from water while the $4.5\mu m$ CO feature remain unchanged. The water band in the HST/WFC3 bandpass is further weakened by H$^-$ continuum opacity. Molecules are expected to recombine before reaching the limb, leading to order of magnitude variations of the chemical composition and cloud coverage between the limb and the dayside. Overall, molecular dissociation provides a qualitative understanding of the lack of strong spectral feature of water in the $1-2\mu m$ bandpass observed in most ultra hot Jupiters. Quantitatively, however, our model does not provide a satisfactory match to the WASP-121b emission spectrum. Together with WASP-33b and Kepler-33Ab, they seem the outliers among the population of ultra hot Jupiters in need of a more thorough understanding.

Atmospheric characterization of directly imaged planets has thus far been limited to ground-based observations of young, self-luminous, Jovian planets. Near-term space- and ground- based facilities like \emph{WFIRST} and ELTs will be able to directly image mature Jovian planets in reflected light, a critical step in support of future facilities that aim to directly image terrestrial planets in reflected light (e.g. HabEx, LUVOIR). These future facilities are considering the use of photometry to classify planets. Here, we investigate the intricacies of using colors to classify gas-giant planets by analyzing a grid of 9,120 theoretical reflected light spectra spread across different metallicities, pressure-temperature profiles, cloud properties, and phase angles. We determine how correlated these planet parameters are with the colors in the \emph{WFIRST} photometric bins and other photometric bins proposed in the literature. Then we outline under what conditions giant planet populations can be classified using several supervised multivariate classification algorithms. We find that giant planets imaged in reflected light can be classified by metallicity with an accuracy of $>$90\% if they are \emph{a prior} known to not have significant cloud coverage in the visible part of the atmosphere, and at least 3 filter observations are available. If the presence of clouds is not known \emph{a priori}, directly imaged planets can be more accurately classified by their cloud properties, as oppposed to metallicity or temperature. Furthermore, we are able to distinguish between cloudy and cloud-free populations with $>$90\% accuracy with 3 filter observations. Our statistical pipeline is available on GitHub and can be extended to optimize science yield of future mission concepts.

UV radiation can induce photochemical processes in exoplanet atmospheres and produce haze particles. Recent observations suggest that haze and/or cloud layers could be present in the upper atmospheres of exoplanets. Haze particles play an important role in planetary atmospheres and may provide a source of organic material to the surface which may impact the origin or evolution of life. However, very little information is known about photochemical processes in cool, high-metallicity exoplanetary atmospheres. Previously, we investigated haze formation and particle size distribution in laboratory atmosphere simulation experiments using AC plasma as the energy source. Here, we use UV photons to initiate the chemistry rather than the AC plasma, since photochemistry driven by UV radiation is important for understanding exoplanet atmospheres. We present photochemical haze formation in current UV experiments, we investigated a range of atmospheric metallicities (100x, 1000x, and 10000x solar metallicity) at three temperatures (300 K, 400 K, and 600 K). We find that photochemical hazes are generated in all simulated atmospheres with temperature-dependent production rates: the particles produced in each metallicity group decrease as the temperature increases. The images taken with atomic force microscopy show the particle size (15-190 nm) varies with temperature and metallicity. Our laboratory experimental results provide new insight into the formation and properties of photochemical haze, which could guide exoplanet atmosphere modeling and help to analyze and interpret current and future observations of exoplanets.

We explore KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR 8799, using PyKLIP and show algorithm stability with varying KLIP parameters. We report new and re-reduced spectrophotometry of HR 8799 c, d, and e in H & K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting simulated sources and recovering them over a range of parameters. The K1/K2 spectra for HR 8799 c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR 8799 e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We show that HR 8799 c and d show significant differences in their H & K spectra, but do not find any conclusive differences between d and e or c and e, likely due to large error bars in the recovered spectrum of e. Compared to M, L, and T-type field brown dwarfs, all three planets are most consistent with mid and late L spectral types. All objects are consistent with low gravity but a lack of standard spectra for low gravity limit the ability to fit the best spectral type. We discuss how dedicated modeling efforts can better fit HR 8799 planets' near-IR flux and discuss how differences between the properties of these planets can be further explored.

The coldest brown dwarf, WISE 0855, is the closest known planetary-mass, free-floating object and has a temperature nearly as cold as the solar system gas giants. Like Jupiter, it is predicted to have an atmosphere rich in methane, water, and ammonia, with clouds of volatile ices. WISE 0855 is faint at near-infrared wavelengths and emits almost all its energy in the mid-infrared. Skemer et al. 2016 presented a spectrum of WISE 0855 from 4.5-5.1 micron (M band), revealing water vapor features. Here, we present a spectrum of WISE 0855 in L band, from 3.4-4.14 micron. We present a set of atmosphere models that include a range of compositions (metallicities and C/O ratios) and water ice clouds. Methane absorption is clearly present in the spectrum. The mid-infrared color can be better matched with a methane abundance that is depleted relative to solar abundance. We find that there is evidence for water ice clouds in the M band spectrum, and we find a lack of phosphine spectral features in both the L and M band spectra. We suggest that a deep continuum opacity source may be obscuring the near-infrared flux, possibly a deep phosphorous-bearing cloud, ammonium dihyrogen phosphate. Observations of WISE 0855 provide critical constraints for cold planetary atmospheres, bridging the temperature range between the long-studied solar system planets and accessible exoplanets. JWST will soon revolutionize our understanding of cold brown dwarfs with high-precision spectroscopy across the infrared, allowing us to study their compositions and cloud properties, and to infer their atmospheric dynamics and formation processes.

The H2O abundance of a planetary atmosphere is a powerful indicator of formation conditions. Inferring H2O in the solar system giant planets is challenging, due to condensation depleting the upper atmosphere of water vapour. Substantially warmer hot Jupiter exoplanets readily allow detections of H2O via transmission spectroscopy, but such signatures are often diminished by the presence of clouds made of other species. In contrast, highly scattering H2O clouds can brighten planets in reflected light, enhancing molecular signatures. Here, we present an extensive parameter space survey of the prominence of H2O absorption features in reflection spectra of cool (Teff<400K) giant exoplanetary atmospheres. The impact of effective temperature, gravity, metallicity, and sedimentation efficiency is explored. We find prominent H2O features around 0.94um, 0.83um, and across a wide spectral region from 0.4-0.73um. The 0.94um feature is only detectable where high-altitude water clouds brighten the planet: T~150K, g>~20m/s^2, fsed>~3, m<~10xsolar. In contrast, planets with g<~20m/s^2 and T>~180K display prominent H2O features embedded in the Rayleigh slope from 0.4-0.73um. High fsed enhances H2O features around 0.94um, and enables them to be detected at lower temperatures. High m results in dampened H2O absorption features, due to H2O vapour condensing to form bright optically thick clouds that dominate the continuum. We verify these trends via self-consistent modelling of the low gravity exoplanet HD 192310c, revealing that its reflection spectrum is expected to be dominated by H2O absorption from 0.4-0.73um for m<~10xsolar. Our results demonstrate that H2O is manifestly detectable in reflected light spectra of cool giant planets only marginally warmer than Jupiter, providing an avenue to directly constrain the C/O and O/H ratios of a hitherto unexplored population of exoplanetary atmospheres.

Transmission spectroscopy during planetary transits is expected to be a major source of information on the atmospheres of small (approximately Earth-sized) exoplanets in the next two decades. This technique, however, is intrinsically affected by stellar spectral contamination caused by the fact that stellar photo- and chromospheres are not perfectly homogeneous. Such stellar contamination will often reach or exceed the signal introduced by the planetary spectral features. Finding effective methods to correct stellar contamination -- or at least to quantify its possible range -- for the most important exoplanets is a necessary step for our understanding of exoplanet atmospheres. This will require significantly deepening our understanding of stellar heterogeneity, which is currently limited by the available data.

Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $\mu$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $\mu$m and $R = 50$ spectroscopy from 0.6-0.97 $\mu$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass.

Ground-based telescopes have been playing a leading role in exoplanet direct imaging science and technological development for the past two decades and will continue to have an indispensable role for the next decade and beyond. Extreme adaptive optics (AO) systems will advance focal-plane wavefront control and coronagraphy, augmenting the performance of and mitigating risk for WFIRST-CGI, while validating performance requirements and motivating improvements to atmosphere models needed to unambiguously characterize solar system-analogues with HabEx/LUVOIR. Specialized instruments for Extremely Large Telescopes may deliver the first thermal infrared images of rocky planets around Sun-like stars, providing HabEx/LUVOIR with numerous exo-Earth candidates and key ancillary information that can help clarify whether the planets are habitable.

Exoplanets mass measurements will be a critical next step to assess the habitability of Earth-like planets: a key aspect of the 2020 vision in the previous decadal survey and also central to NASA's strategic priorities. Precision astrometry delivers measurement of exoplanet masses, allowing discrimination of rocky planets from water worlds and enabling much better modeling of their atmosphere improving species retrieval from spectroscopy. The scientific potential of astrometry will be enormous. The intrinsic astrophysical noise floor set by star spots and stellar surface activity is about a factor of ten more benign for astrometry than for the more established technique of Radial Velocity, widening the discovery region and pushing detection thresholds to lower masses than previously possible. On the instrumental side, precision astrometry is limited by optical field distortion and detector calibration issues. Both technical challenges are now being addressed successfully in the laboratory. However, we have identified the need to continue these technology development efforts to achieve sub-microarcsecond astrometry precision necessary for detection and characterization of Earth-like planets around nearby FGK stars. The international community has realized the importance of astrometry, and various astrometry missions have been proposed and under development, with a few high profile missions now operational. We believe that it is vital for the U.S. scientific community to participate in the development of these new technologies and scientific discoveries. We recommend exploring alternatives to incorporate astrometric capabilities into future exoplanet flagship missions such as HABEX and LUVOIR, substantially increasing the scientific return associated with the expected yield of earth-like planets to be recovered.

Condensation clouds in substellar atmospheres have been widely inferred from spectra and photometric variability. Up until now, their horizontally averaged vertical distribution and mean particle size have been largely characterized using models, one of which is the eddy diffusion-sedimentation model from Ackerman & Marley (2001) that relies on a sedimentation efficiency parameter, $f_{\rm sed}$, to determine the vertical extent of clouds in the atmosphere. However, the physical processes controlling the vertical structure of clouds in substellar atmospheres are not well understood. In this work, we derive trends in $f_{\rm sed}$ across a large range of eddy diffusivities ($K_{zz}$), gravities, material properties, and cloud formation pathways by fitting cloud distributions calculated by a more detailed cloud microphysics model. We find that $f_{\rm sed}$ is dependent on $K_{zz}$, but not gravity, when $K_{zz}$ is held constant. $f_{\rm sed}$ is most sensitive to the nucleation rate of cloud particles, as determined by material properties like surface energy and molecular weight. High surface energy materials form fewer, larger cloud particles, leading to large $f_{\rm sed}$ ($>$1), and vice versa for materials with low surface energy. For cloud formation via heterogeneous nucleation, $f_{\rm sed}$ is sensitive to the condensation nuclei flux and radius, connecting cloud formation in substellar atmospheres to the objects' formation environments and other atmospheric aerosols. These insights could lead to improved cloud models that help us better understand substellar atmospheres. For example, we demonstrate that $f_{\rm sed}$ could increase with increasing cloud base depth in an atmosphere, shedding light on the nature of the brown dwarf L/T transition.

Time-resolved observations of brown dwarfs' rotational modulations provide powerful insights into the properties of condensate clouds in ultra-cool atmospheres. Multi-wavelength light curves reveal cloud vertical structures, condensate particle sizes, and cloud morphology, which directly constrain condensate cloud and atmospheric circulation models. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared G141 taken in six consecutive orbits observations of HN Peg B, an L/T transition brown dwarf companion to a G0V type star. The best-fit sine wave to the $1.1-1.7\mu$m broadband light curve has the amplitude of $1.206\pm0.025\%$ and period of $15.4\pm0.5$ hr. The modulation amplitude has no detectable wavelength dependence except in the 1.4 $\mu$m water absorption band, indicating that the characteristic condensate particle sizes are large ($>1\mu$m). We detect significantly ($4.4\sigma$) lower modulation amplitude in the 1.4$\mu$m water absorption band, and find that HN Peg B's spectral modulation resembles those of early T type brown dwarfs. We also describe a new empirical interpolation method to remove spectral contamination from the bright host star. This method may be applied in other high-contrast time-resolved observations with WFC3.

Numerous solar system atmospheres possess aerosols including the characteristic organic hazes of Titan and Pluto. Haze particles substantially impact atmospheric temperatures structures and may provide organic material to the surface of a world, thereby affecting its habitability. Observations of exoplanet atmospheres suggest the presence of aerosols, especially in cooler (<800 K), smaller (<0.3 times Jupiter's mass) exoplanets. It remains unclear if the aerosols muting the spectroscopic features of exoplanet atmospheres are condensate clouds or photochemical hazes, which is difficult to predict from theory alone. We present here the first laboratory haze simulation experiments for atmospheric compositions expected for super-Earths and mini-Neptunes. We explored temperatures from 300 to 600 K and a range of atmospheric metallicities (100x, 1000x, 10000x solar); all simulated atmospheres produced particles, and the cooler (300 and 400 K) 1000x solar metallicity ("H2O-dominated", CH4-rich) experiments exhibited haze production rates higher than our standard Titan simulation (~10 mg/hr versus 7.4 mg/hr for Titan). However the particle production rates varied greatly, with measured rates as low as 0.04 mg/hr (100x solar metallicity, 600 K). Here we show that we should expect some, but not all, super-Earth and mini-Neptune atmospheres to possess a thick photochemically generated haze.

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.

Observations of rotational modulations of brown dwarfs and giant exoplanets allow the characterization of condensate cloud properties. As of now rotational spectral modulations have only been seen in three L-type brown dwarfs. We report here the discovery of rotational spectral modulations in LP261-75B, an L6-type intermediate surface gravity companion to an M4.5 star. As a part of the Cloud Atlas Treasury program we acquired time-resolved Wide Field Camera 3 grism spectroscopy (1.1--1.69~$\mu$m) of LP261-75B. We find gray spectral variations with the relative amplitude displaying only a weak wavelength dependence and no evidence for lower-amplitude modulations in the 1.4~$\mu$m water band than in the adjacent continuum. The likely rotational modulation period is 4.78$\pm$0.95 h, although the rotational phase is not well sampled. The minimum relative amplitude in the white light curve measured over the whole wavelength range is 2.41$\pm$0.14\%. We report an unusual light curve with seemingly three peaks approximately evenly distributed in rotational phase. The spectral modulations suggests that the upper atmosphere cloud properties in {LP261-75B} are similar to two other mid-L dwarfs of typical infrared colors, but differ from that of the extremely red L-dwarf WISE0047.

Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

Brown dwarf spectra are rich in information revealing of the chemical and physical processes operating in their atmospheres. We apply a recently developed atmospheric retrieval tool to an ensemble of late T-dwarf (600-800K) near infrared spectra. With these spectra we are able to place direct constraints the molecular abundances of H$_2$O, CH$_4$, CO, CO$_2$, NH$_3$, H$_2$S, and Na+K, gravity, thermal structure (and effective temperature), photometric radius, and cloud optical depths. We find that ammonia, water, methane, and the alkali metals are present and well constrained in all 11 objects. From the abundance constraints we find no significant trend in the water, methane, or ammonia abundances with temperature, but find a very strong ($>$25$\sigma$) increasing trend in the alkali metal abundances with effective temperature, indicative of alkali rainout. We also find little evidence for optically thick clouds. With the methane and water abundances, we derive the intrinsic atmospheric metallicity and carbon-to-oxygen ratios. We find in our sample, that metallicities are typically sub solar and carbon-to-oxygen ratios are somewhat super solar, different than expectations from the local stellar population. We also find that the retrieved vertical thermal profiles are consistent with radiative equilibrium over the photospheric regions. Finally, we find that our retrieved effective temperatures are lower than previous inferences for some objects and that our radii are larger than expectations from evolutionary models, possibly indicative of un-resolved binaries. This investigation and methodology represents a paradigm in linking spectra to the determination of the fundamental chemical and physical processes governing cool brown dwarf atmospheres.

Remote sensing of the atmospheres of distant worlds motivates a firm understanding of radiative transfer. In this review, we provide a pedagogical cookbook that describes the principal ingredients needed to perform a radiative transfer calculation and predict the spectrum of an exoplanet atmosphere, including solving the radiative transfer equation, calculating opacities (and chemistry), iterating for radiative equilibrium (or not), and adapting the output of the calculations to the astronomical observations. A review of the state of the art is performed, focusing on selected milestone papers. Outstanding issues, including the need to understand aerosols or clouds and elucidating the assumptions and caveats behind inversion methods, are discussed. A checklist is provided to assist referees/reviewers in their scrutiny of works involving radiative transfer. A table summarizing the methodology employed by past studies is provided.

Observations of the highly-eccentric (e~0.9) hot-Jupiter HD 80606b with Spitzer have provided some of best probes of the physics at work in exoplanet atmospheres. By observing HD 80606b during its periapse passage, atmospheric radiative, advective, and chemical timescales can be directly measured and used to constrain fundamental planetary properties such as rotation period, tidal dissipation rate, and atmospheric composition (including aerosols). Here we present three-dimensional general circulation models for HD 80606b that aim to further explore the atmospheric physics shaping HD 80606b's observed Spitzer phase curves. We find that our models that assume a planetary rotation period twice that of the pseudo-synchronous rotation period best reproduce the phase variations observed for HD~80606b near periapse passage with Spitzer. Additionally, we find that the rapid formation/dissipation and vertical transport of clouds in HD 80606b's atmosphere near periapse passage likely shapes its observed phase variations. We predict that observations near periapse passage at visible wavelengths could constrain the composition and formation/advection timescales of the dominant cloud species in HD 80606b's atmosphere. The time-variable forcing experienced by exoplanets on eccentric orbits provides a unique and important window on radiative, dynamical, and chemical processes in planetary atmospheres and an important link between exoplanet observations and theory.

We present the first results from applying the spectral inversion technique in the cloudy L dwarf regime. Our new framework provides a flexible approach to modelling cloud opacity which can be built incrementally as the data requires, and improves upon previous retrieval experiments in the brown dwarf regime by allowing for scattering in two stream radiative transfer. Our first application of the tool to two mid-L dwarfs is able to reproduce their near-infrared spectra far more closely than grid models. Our retrieved thermal, chemical, and cloud profiles allow us to estimate $T_{\rm eff} = 1796^{+23}_{-25}$ K and $\log g = 5.21^{+0.05}_{-0.08}$ for 2MASS J05002100+0330501 and for 2MASSW J2224438-015852 we find $T_{\rm eff} = 1723^{+18}_{-19}$ K and $\log g = 5.31^{+0.04}_{-0.08}$, in close agreement with previous empirical estimates. Our best model for both objects includes an optically thick cloud deck which passes $\tau_{cloud} \geq 1$ (looking down) at a pressure of around 5 bar. The temperature at this pressure is too high for silicate species to condense, and we argue that corundum and/or iron clouds are responsible for this cloud opacity. Our retrieved profiles are cooler at depth, and warmer at altitude than the forward grid models that we compare, and we argue that some form of heating mechanism may be at work in the upper atmospheres of these L dwarfs. We also identify anomalously high CO abundance in both targets, which does not correlate with the warmth of our upper atmospheres or our choice of cloud model, and find similarly anomalous alkali abundance for one of our targets. These anomalies may reflect unrecognised shortcomings in our retrieval model, or inaccuracies in our gas phase opacities.

We present a new matched filter algorithm for direct detection of point sources in the immediate vicinity of bright stars. The stellar Point Spread Function (PSF) is first subtracted using a Karhunen-Lo\'eve Image Processing (KLIP) algorithm with Angular and Spectral Differential Imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the Signal-to-Noise Ratio (SNR) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal SNR loss. We also developed a complete pipeline for the automated detection of point source candidates, the calculation of Receiver Operating Characteristics (ROC), false positives based contrast curves, and completeness contours. We process in a uniform manner more than 330 datasets from the Gemini Planet Imager Exoplanet Survey (GPIES) and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false positive rate. We show that the new forward model matched filter allows the detection of $50\%$ fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false positive rate.

We present spectro-photometry spanning 1-5 $\mu$m of 51 Eridani b, a 2-10 M$_\text{Jup}$ planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new $K1$ (1.90-2.19 $\mu$m) and $K2$ (2.10-2.40 $\mu$m) spectra taken with the Gemini Planet Imager as well as an updated $L_P$ (3.76 $\mu$m) and new $M_S$ (4.67 $\mu$m) photometry from the NIRC2 Narrow camera. The new data were combined with $J$ (1.13-1.35 $\mu$m) and $H$ (1.50-1.80 $\mu$m) spectra from the discovery epoch with the goal of better characterizing the planet properties. 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4-T8) and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud free and patchy/intermediate opacity clouds. Model fits suggest that 51 Eri b has an effective temperature ranging between 605-737 K, a solar metallicity, a surface gravity of $\log$(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the SED. From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 ($\log L/L_{\odot}$), leaving 51 Eri b in the unique position as being one of the only directly imaged planet consistent with having formed via cold-start scenario. Comparisons of the planet SED against warm-start models indicates that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M$_{\oplus}$.

We present H band spectroscopic and Halpha photometric observations of HD 100546 obtained with GPI and MagAO. We detect H band emission at the location of the protoplanet HD 100546b, but show that choice of data processing parameters strongly affects the morphology of this source. It appears point-like in some aggressive reductions, but rejoins an extended disk structure in the majority of the others. Furthermore, we demonstrate that this emission appears stationary on a timescale of 4.6 yrs, inconsistent at the 2sigma level with a Keplerian clockwise orbit at 59 au in the disk plane. The H band spectrum of the emission is inconsistent with any type of low effective temperature object or accreting protoplanetary disk. It strongly suggests a scattered light origin, as it is consistent with the spectrum of the star and the spectra extracted at other locations in the disk. A non detection at the 5sigma level of HD 100546b in differential Halpha imaging places an upper limit, assuming the protoplanet lies in a gap free of extinction, on the accretion luminosity and accretion rate of 1.7E-4 Lsun and MMdot<6.4E-7Mjup^2/yr for 1Rjup. These limits are comparable to the accretion luminosity and rate of TTauri-stars or LkCa 15b. Taken together, these lines of evidence suggest that the H band source at the location of HD 100546b is not emitted by a planetary photosphere or an accreting circumplanetary disk but is a disk feature enhanced by the PSF subtraction process. This non-detection is consistent with the non-detection in the K band reported in an earlier study but does not exclude the possibility that HD 100546b is deeply embedded.

We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.

Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.

Recent work has shown that sulfur hazes may arise in the atmospheres of some giant exoplanets due to the photolysis of H$_{2}$S. We investigate the impact such a haze would have on an exoplanet's geometric albedo spectrum and how it may affect the direct imaging results of WFIRST, a planned NASA space telescope. For temperate (250 K $<$ T$_{\rm eq}$ $<$ 700 K) Jupiter--mass planets, photochemical destruction of H$_{2}$S results in the production of $\sim$1 ppmv of \seight between 100 and 0.1 mbar, which, if cool enough, will condense to form a haze. Nominal haze masses are found to drastically alter a planet's geometric albedo spectrum: whereas a clear atmosphere is dark at wavelengths between 0.5 and 1 $\mu$m due to molecular absorption, the addition of a sulfur haze boosts the albedo there to $\sim$0.7 due to scattering. Strong absorption by the haze shortward of 0.4 $\mu$m results in albedos $<$0.1, in contrast to the high albedos produced by Rayleigh scattering in a clear atmosphere. As a result, the color of the planet shifts from blue to orange. The existence of a sulfur haze masks the molecular signatures of methane and water, thereby complicating the characterization of atmospheric composition. Detection of such a haze by WFIRST is possible, though discriminating between a sulfur haze and any other highly reflective, high altitude scatterer will require observations shortward of 0.4 $\mu$m, which is currently beyond WFIRST's design.

We report here on our search for excess power in photometry of Neptune collected by the K2 mission that may be due to intrinsic global oscillations of the planet Neptune. To conduct this search, we developed new methods to correct for instrumental effects such as intrapixel variability and gain variations. We then extracted and analyzed the time-series photometry of Neptune from 49 days of nearly continuous broadband photometry of the planet. We find no evidence of global oscillations and place an upper limit of $\sim$5 ppm at 1000 \uhz\ for the detection of a coherent signal. With an observed cadence of 1-minute and point-to-point scatter less than 0.01\%, the photometric signal is dominated by reflected light from the Sun, which is in turn modulated by atmospheric variability of Neptune at the 2\% level. A change in flux is also observed due to the increasing distance between Neptune and the K2 spacecraft, and solar variability with convection-driven solar p modes present.

We present results from the first observations of the Hubble Space Telescope (HST) Panchromatic Comparative Exoplanet Treasury (PanCET) program for WASP-101b, a highly inflated hot Jupiter and one of the community targets proposed for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. From a single HST Wide Field Camera 3 (WFC3) observation, we find that the near-infrared transmission spectrum of WASP-101b contains no significant H$_2$O absorption features and we rule out a clear atmosphere at 13{\sigma}. Therefore, WASP-101b is not an optimum target for a JWST ERS program aimed at observing strong molecular transmission features. We compare WASP-101b to the well studied and nearly identical hot Jupiter WASP-31b. These twin planets show similar temperature-pressure profiles and atmospheric features in the near-infrared. We suggest exoplanets in the same parameter space as WASP-101b and WASP-31b will also exhibit cloudy transmission spectral features. For future HST exoplanet studies, our analysis also suggests that a lower count limit needs to be exceeded per pixel on the detector in order to avoid unwanted instrumental systematics.

Starting in December 2014, Kepler K2 observed Neptune continuously for 49 days at a 1-minute cadence. The goals consisted of studying its atmospheric dynamics (Simon et al. 2016), detecting its global acoustic oscillations (Rowe et al., submitted), and those of the Sun, which we report on here. We present the first indirect detection of solar oscillations in intensity measurements. Beyond the remarkable technical performance, it indicates how Kepler would see a star like the Sun. The result from the global asteroseismic approach, which consists of measuring the oscillation frequency at maximum amplitude "nu_max" and the mean frequency separation between mode overtones "Delta nu", is surprising as the nu_max measured from Neptune photometry is larger than the accepted value. Compared to the usual reference nu_max_sun = 3100 muHz, the asteroseismic scaling relations therefore make the solar mass and radius appear larger by 13.8 +/- 5.8 % and 4.3 +/- 1.9 % respectively. The higher nu_max is caused by a combination of the value of nu_max_sun, being larger at the time of observations than the usual reference from SOHO/VIRGO/SPM data (3160 +/- 10 muHz), and the noise level of the K2 time series, being ten times larger than VIRGO's. The peak-bagging method provides more consistent results: despite a low signal-to-noise ratio (SNR), we model ten overtones for degrees l=0,1,2. We compare the K2 data with simultaneous SOHO/VIRGO/SPM photometry and BiSON velocity measurements. The individual frequencies, widths, and amplitudes mostly match those from VIRGO and BiSON within 1 sigma, except for the few peaks with lowest SNR.

Future space-based telescopes, such as the Wide-Field Infrared Survey Telescope (WFIRST), will observe the reflected-light spectra of directly imaged extrasolar planets. Interpretation of such data presents a number of novel challenges, including accounting for unknown planet radius and uncertain stellar illumination phase angle. Here we report on our continued development of Markov Chain Monte Carlo retrieval methods for addressing these issues in the interpretation of such data. Specifically we explore how the unknown planet radius and potentially poorly known observer-planet-star phase angle impacts retrievals of parameters of interest such as atmospheric methane abundance, cloud properties and surface gravity. As expected, the uncertainty in retrieved values is a strong function of signal-to-noise ratio (SNR) of the observed spectra, particularly for low metallicity atmospheres, which lack deep absorption signatures. Meaningful results may only be possible above certain SNR thresholds; for cases across a metallicity range of 1-50 times solar, we find that only an SNR of 20 systematically reproduces close to the correct methane abundance at all phase angles. However, even in cases where the phase angle is poorly known we find that the planet radius can be constrained to within a factor of two. We find that uncertainty in planet radius decreases at phase angles past quadrature, as the highly forward scattering nature of the atmosphere at these geometries limits the possible volume of phase space that relevant parameters can occupy. Finally, we present an estimation of possible improvement that can result from combining retrievals against observations at multiple phase angles.

The Neptune-mass GJ 436b is one of the most-studied transiting exoplanets with repeated measurements of both its thermal emission and transmission spectra. We build on previous studies to answer outstanding questions about this planet, including its potentially high metallicity and tidal heating of its interior. We present new observations of GJ 436b's thermal emission at 3.6 and 4.5 micron, which reduce uncertainties in estimates of GJ 436b's flux at those wavelengths and demonstrate consistency between Spitzer observations spanning more than 7 years. We analyze the Spitzer thermal emission photometry and Hubble WFC3 transmission spectrum in tandem. We use a powerful dual-pronged modeling approach, comparing these data to both self-consistent and retrieval models. We vary the metallicity, intrinsic luminosity from tidal heating, disequilibrium chemistry, and heat redistribution. We also study the effect of clouds and photochemical hazes on the spectra, but do not find strong evidence for either. The self-consistent and retrieval modeling combine to suggest that GJ 436b has a high atmospheric metallicity, with best fits at or above several hundred times solar metallicity, tidal heating warming its interior with best-fit intrinsic effective effective temperatures around 300--350 K, and disequilibrium chemistry. High metal-enrichments (>600x solar) can only occur from the accretion of rocky, rather than icy, material. Assuming Tint~300--350 K, we find that Q'~2x10^5--10^6, larger than Neptune's Q', and implying a long tidal circularization timescale for the planet's orbit. We suggest that Neptune-mass planets may be a more diverse class than previously imagined, with metal-enhancements potentially spanning several orders of magnitude, to perhaps over 1000x solar metallicity. High fidelity observations with instruments like JWST will be critical for characterizing this diversity.

We present phase curves of Jupiter from 0-140 degrees as measured in multiple optical bandpasses by Cassini/ISS during the Millennium flyby of Jupiter in late 2000 to early 2001. Phase curves are of interest for studying the energy balance of Jupiter and understanding the scattering behavior of Jupiter as an exoplanet analog. We find that Jupiter is significantly darker at partial phases than an idealized Lambertian planet by roughly 25% and is not well fit by Jupiter-like exoplanet atmospheric models across all wavelengths. We provide analytic fits to Jupiter's phase function in several Cassini/ISS imaging filter bandpasses. In addition, these observations show that Jupiter's color is more variable with phase angle than predicted by models. Therefore, the color of even a near Jupiter-twin planet observed at a partial phase cannot be assumed to be comparable to that of Jupiter at full phase. We discuss how WFIRST and other future direct-imaging missions can enhance the study of cool giants.

Deciphering the role of clouds is central to our understanding of exoplanet atmospheres, as they have a direct impact on the temperature and pressure structure, and observational properties of the planet. Super-hot Jupiters occupy a temperature regime similar to low mass M-dwarfs, where minimal cloud condensation is expected. However, observations of exoplanets such as WASP-12b (Teq ~ 2500 K) result in a transmission spectrum indicative of a cloudy atmosphere. We re-examine the temperature and pressure space occupied by these super-hot Jupiter atmospheres, to explore the role of the initial Al- and Ti-bearing condensates as the main source of cloud material. Due to the high temperatures a majority of the more common refractory material is not depleted into deeper layers and would remain in the vapor phase. The lack of depletion into deeper layers means that these materials with relatively low cloud masses can become significant absorbers in the upper atmosphere. We provide condensation curves for the initial Al- and Ti-bearing condensates that may be used to provide quantitative estimates of the effect of metallicity on cloud masses, as planets with metal-rich hosts potentially form more opaque clouds because more mass is available for condensation. Increased metallicity also pushes the point of condensation to hotter, deeper layers in the planetary atmosphere further increasing the density of the cloud. We suggest that planets around metal-rich hosts are more likely to have thick refractory clouds, and discuss the implication on the observed spectra of WASP-12b.

We present new spatially resolved astrometry and photometry from the Gemini Planet Imager of the inner binary of the young multiple star system V343 Normae, which is a member of the beta Pictoris moving group. V343 Normae comprises a K0 and mid-M star in a ~4.5 year orbit (AaAb) and a wide 10" M5 companion (B). By combining these data with archival astrometry and radial velocities we fit the orbit and measure individual masses for both components of M_Aa = 1.10 +/- 0.10 M_sun and M_Ab = 0.290 +/- 0.018 M_sun. Comparing to theoretical isochrones, we find good agreement for the measured masses and JHK band magnitudes of the two components consistent with the age of the beta Pic moving group. We derive a model-dependent age for the beta Pic moving group of 26 +/- 3 Myr by combining our results for V343 Normae with literature measurements for GJ 3305, which is another group member with resolved binary components and dynamical masses.

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

Condensate clouds fundamentally impact the atmospheric structure and spectra of exoplanets and brown dwarfs but the connections between surface gravity, cloud structure, dust in the upper atmosphere, and the red colors of some brown dwarfs remain poorly understood. Rotational modulations enable the study of different clouds in the same atmosphere, thereby providing a method to isolate the effects of clouds. Here we present the discovery of high peak-to-peak amplitude (8%) rotational modulations in a low-gravity, extremely red (J-Ks=2.55) L6 dwarf WISEP J004701.06+680352.1 (W0047). Using the Hubble Space Telescope (HST) time-resolved grism spectroscopy we find a best-fit rotational period (13.20$\pm$0.14 hours) with a larger amplitude at 1.1 micron than at 1.7 micron. This is the third largest near-infrared variability amplitude measured in a brown dwarf, demonstrating that large-amplitude variations are not limited to the L/T transition but are present in some extremely red L-type dwarfs. We report a tentative trend between the wavelength dependence of relative amplitude, possibly proxy for small dust grains lofted in the upper atmosphere, and the likelihood of large-amplitude variability. By assuming forsterite as haze particle, we successfully explain the wavelength dependent amplitude with submicron-sized haze particles sizes of around 0.4 {\mu}m. W0047 links the earlier spectral and later spectral type brown dwarfs in which rotational modulations have been observed, the large amplitude variations in this object make this a benchmark brown dwarf for the study of cloud properties close to the L/T transition.

Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler (emcee) and a multimodal nested sampling algorithm (MultiNest) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model, and highlights possible discrepancies in the likelihood maps. As a proof-of-concept, our current atmospheric model contains 1 or 2 cloud layers, methane as a major absorber, and a H$_2$-He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492 c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model-dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets.

We present $H$-band near-infrared polarimetric imaging observations of the F5V star HD~157587 obtained with the Gemini Planet Imager (GPI) that reveal the debris disk as a bright ring structure at a separation of $\sim$80$-$100~AU. The new GPI data complement recent HST/STIS observations that show the disk extending out to over 500~AU. The GPI image displays a strong asymmetry along the projected minor axis as well as a fainter asymmetry along the projected major axis. We associate the minor and major axis asymmetries with polarized forward scattering and a possible stellocentric offset, respectively. To constrain the disk geometry we fit two separate disk models to the polarized image, each using a different scattering phase function. Both models favor a disk inclination of $\sim 70\degr$ and a $1.5\pm0.6$ AU stellar offset in the plane of the sky along the projected major axis of the disk. We find that the stellar offset in the disk plane, perpendicular to the projected major axis is degenerate with the form of the scattering phase function and remains poorly constrained. The disk is not recovered in total intensity due in part to strong adaptive optics residuals, but we recover three point sources. Considering the system's proximity to the galactic plane and the point sources' positions relative to the disk, we consider it likely that they are background objects and unrelated to the disk's offset from the star.

The past decade has seen significant progress on the direct detection and characterization of young, self-luminous giant planets at wide orbital separations from their host stars. Some of these planets show evidence for disequilibrium processes like transport-induced quenching in their atmospheres; photochemistry may also be important, despite the large orbital distances. These disequilibrium chemical processes can alter the expected composition, spectral behavior, thermal structure, and cooling history of the planets, and can potentially confuse determinations of bulk elemental ratios, which provide important insights into planet-formation mechanisms. Using a thermo/photochemical kinetics and transport model, we investigate the extent to which disequilibrium chemistry affects the composition and spectra of directly imaged giant exoplanets. Results for specific "young Jupiters" such as HR 8799 b and 51 Eri b are presented, as are general trends as a function of planetary effective temperature, surface gravity, incident ultraviolet flux, and strength of deep atmospheric convection. We find that quenching is very important on young Jupiters, leading to CO/CH4 and N2/NH3 ratios much greater than, and H2O mixing ratios a factor of a few less than, chemical-equilibrium predictions. Photochemistry can also be important on such planets, with CO2 and HCN being key photochemical products. Carbon dioxide becomes a major constituent when stratospheric temperatures are low and recycling of water via the H2 + OH reaction becomes kinetically stifled. Young Jupiters with effective temperatures <~ 700 K are in a particularly interesting photochemical regime that differs from both transiting hot Jupiters and our own solar-system giant planets.

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 present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo\`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $\beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$\sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $\beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$\sigma$ significance. However, we confirm that the Hill sphere of $\beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

We have monitored photometrically the Y0 brown dwarf WISEP J173835.52+273258.9 (W1738) at both near- and mid-infrared wavelengths. This ~1 Gyr-old 400K dwarf is at a distance of 8pc and has a mass around 5 M_Jupiter. We observed W1738 using two near-infrared filters at lambda~1um, Y and J, on Gemini observatory, and two mid-infrared filters at lambda~4um, [3.6] and [4.5], on the Spitzer observatory. Twenty-four hours were spent on the source by Spitzer on each of June 30 and October 30 2013 UT. Between these observations, around 5 hours were spent on the source by Gemini on each of July 17 and August 23 2013 UT. The mid-infrared light curves show significant evolution between the two observations separated by four months. We find that a double sinusoid can be fit to the [4.5] data, where one sinusoid has a period of 6.0 +/- 0.1 hours and the other a period of 3.0 +/- 0.1 hours. The near-infrared observations suggest variability with a ~3.0 hour period, although only at a <~2 sigma confidence level. We interpret our results as showing that the Y dwarf has a 6.0 +/- 0.1 hour rotation period, with one or more large-scale surface features being the source of variability. The peak-to-peak amplitude of the light curve at [4.5] is 3%. The amplitude of the near-infrared variability, if real, may be as high as 5 to 30%. Intriguingly, this size of variability and the wavelength dependence can be reproduced by atmospheric models that include patchy KCl and Na_2S clouds and associated small changes in surface temperature. The small number of large features, and the timescale for evolution of the features, is very similar to what is seen in the atmospheres of the solar system gas giants.

Over a large range of equilibrium temperatures, clouds shape the transmission spectrum of hot Jupiter atmospheres, yet their composition remains unknown. Recent observations show that the Kepler lightcurves of some hot Jupiters are asymmetric: for the hottest planets, the lightcurve peaks before secondary eclipse, whereas for planets cooler than $\sim1900\rm\,K$, it peaks after secondary eclipse. We use the thermal structure from 3D global circulation models to determine the expected cloud distribution and Kepler lightcurves of hot Jupiters. We demonstrate that the change from an optical lightcurve dominated by thermal emission to one dominated by scattering (reflection) naturally explains the observed trend from negative to positive offset. For the cool planets the presence of an asymmetry in the Kepler lightcurve is a telltale sign of the cloud composition, because each cloud species can produce an offset only over a narrow range of effective temperatures. By comparing our models and the observations, we show that the cloud composition of hot Jupiters likely varies with equilibrium temperature. We suggest that a transition occurs between silicate and manganese sulfide clouds at a temperature near $1600\rm\,K$, analogous to the L/T transition on brown dwarfs. The cold trapping of cloud species below the photosphere naturally produces such a transition and predicts similar transitions for other condensates, including TiO. We predict that most hot Jupiters should have cloudy nightsides, that partial cloudiness should be common at the limb and that the dayside hot spot should often be cloud-free.

We examine the hypothesis that the red near-infrared colors of some L dwarfs could be explained by a "dust haze" of small particles in their upper atmospheres. This dust haze would exist in conjunction with the clouds found in dwarfs with more typical colors. We developed a model which uses Mie theory and the Hansen particle size distributions to reproduce the extinction due to the proposed dust haze. We apply our method to 23 young L dwarfs and 23 red field L dwarfs. We constrain the properties of the dust haze including particle size distribution and column density using Markov-Chain Monte Carlo methods. We find that sub-micron range silicate grains reproduce the observed reddening. Current brown dwarf atmosphere models include large grain (1--100~$\mu m$) dust clouds but not submicron dust grains. Our results provide a strong proof of concept and motivate a combination of large and small dust grains in brown dwarf atmosphere models.

We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection.

The recently discovered brown dwarf WISE 0855 presents our first opportunity to directly study an object outside the Solar System that is nearly as cold as our own gas giant planets. However the traditional methodology for characterizing brown dwarfs---near infrared spectroscopy---is not currently feasible as WISE 0855 is too cold and faint. To characterize this frozen extrasolar world we obtained a 4.5-5.2 $\mu$m spectrum, the same bandpass long used to study Jupiter's deep thermal emission. Our spectrum reveals the presence of atmospheric water vapor and clouds, with an absorption profile that is strikingly similar to Jupiter. The spectrum is high enough quality to allow the investigation of dynamical and chemical processes that have long been studied in Jupiter's atmosphere, but now on an extrasolar world.

We present $Spitzer$/IRAC Ch1 and Ch2 monitoring of six brown dwarfs during 8 different epochs over the course of 20 months. For four brown dwarfs, we also obtained simulataneous $HST$/WFC3 G141 Grism spectra during two epochs and derived light curves in five narrow-band filters. Probing different pressure levels in the atmospheres, the multi-wavelength light curves of our six targets all exhibit variations, and the shape of the light curves evolves over the timescale of a rotation period, ranging from 1.4 h to 13 h. We compare the shapes of the light curves and estimate the phase shifts between the light curves observed at different wavelengths by comparing the phase of the primary Fourier components. We use state-of-the-art atmosphere models to determine the flux contribution of different pressure layers to the observed flux in each filter. We find that the light curves that probe higher pressures are similar and in phase, but are offset and often different from the light curves that probe lower pressures. The phase differences between the two groups of light curves suggest that the modulations seen at lower and higher pressures may be introduced by different cloud layers.

WISE J104915.57-531906.1 is the nearest brown dwarf binary to our Solar system, consisting of two brown dwarfs in the L/T transition: Luhman 16A & B. In this paper we present the first map of Luhman 16A, and maps of Luhman 16B for two epochs. Our maps were created by applying Aeolus, a Markov-Chain Monte Carlo code that maps the top-of-the-atmosphere structure of brown dwarf and other ultracool atmospheres, to light curves of Luhman 16A & B using the Hubble Space Telescope's G141 and G102 grisms. Aeolus retrieved three or four spots in the top-of-the-atmosphere of Luhman 16A & B, with a surface coverage of 19%-32% (depending on an assumed rotational period of 5 hr or 8 hr) or 21%-38.5% (depending on the observational epoch) respectively. The brightness temperature of the spots of the best-fit models was ~200 K hotter than the background top-of-the-atmosphere. We compared our Luhman 16B map with the only previously published map. Interestingly, our map contained a large, cooler (DT~51 K) than the background top-of-the-atmosphere spot that lay at low latitudes, in agreement with the previous Luhman 16B map. Finally, we report the detection of a feature reappearing in Luhman 16B light curves that are separated by tens of hundreds of rotations from each other. We speculate this feature is related to top-of-the-atmosphere structures of Luhman 16B.

We use a 1D model to address photochemistry and possible haze formation in the irradiated warm Jupiter, 51 Eridani b. The intended focus was to be carbon, but sulfur photochemistry turns out to be important. The case for organic photochemical hazes is intriguing but falls short of being compelling. If organic hazes form, they are likeliest to do so if vertical mixing in 51 Eri b is weaker than in Jupiter, and they would be found below the altitudes where methane and water are photolyzed. The more novel result is that photochemistry turns H$_2$S into elemental sulfur, here treated as S$_8$. In the cooler models, S$_8$ is predicted to condense in optically thick clouds of solid sulfur particles, whilst in the warmer models S$_8$ remains a vapor along with several other sulfur allotropes that are both visually striking and potentially observable. For 51 Eri b, the division between models with and without condensed sulfur is at an effective temperature of 700 K, which is within error its actual effective temperature; the local temperature where sulfur condenses is between 280 and 320 K. The sulfur photochemistry we have discussed is quite general and ought to be found in a wide variety of worlds over a broad temperature range, both colder and hotter than the 650-750 K range studied here, and we show that products of sulfur photochemistry will be nearly as abundant on planets where the UV irradiation is orders of magnitude weaker than it is on 51 Eri b.

Many of the directly imaged self-luminous gas giant exoplanets have been found to have cloudy atmospheres. Scattering of the emergent thermal radiation from these planets by the dust grains in their atmospheres should locally give rise to significant linear polarization of the emitted radiation. However, the observable disk averaged polarization should be zero if the planet is spherically symmetric. Rotation-induced oblateness may yield a net non-zero disk averaged polarization if the planets have sufficiently high spin rotation velocity. On the other hand, when a large natural satellite or exomoon transits a planet with cloudy atmosphere along the line of sight, the asymmetry induced during the transit should give rise to a net non-zero, time resolved linear polarization signal. The peak amplitude of such time dependent polarization may be detectable even for slowly rotating exoplanets. Therefore, we suggest that large exomoons around directly imaged self-luminous exoplanets may be detectable through time resolved imaging polarimetry. Adopting detailed atmospheric models for several values of effective temperature and surface gravity which are appropriate for self-luminous exoplanets, we present the polarization profiles of these objects in the infrared during transit phase and estimate the peak amplitude of polarization that occurs during the inner contacts of the transit ingress/egress phase. The peak polarization is predicted to range between 0.1 and 0.3 % in the infrared.

We present new near-infrared spectra, obtained at Gemini Observatory, for two Y dwarfs: WISE J035000.32-565830.2 (W0350) and WISEP J173835.52+273258.9 (W1738). A FLAMINGOS-2 R=540 spectrum was obtained for W0350, covering 1.0 < lambda um < 1.7, and a cross-dispersed GNIRS R=2800 spectrum was obtained for W1738, covering 0.993-1.087 um, 1.191-1.305 um, 1.589-1.631 um, and 1.985-2.175 um, in four orders. We also present revised YJH photometry for W1738, using new NIRI Y and J imaging, and a re-analysis of the previously published NIRI H band images. We compare these data, together with previously published data for late-T and Y dwarfs, to cloud-free models of solar metallicity, calculated both in chemical equilibrium and with disequilibrium driven by vertical transport. We find that for the Y dwarfs the non-equilibrium models reproduce the near-infrared data better than the equilibrium models. The remaining discrepancies suggest that fine-tuning the CH_4/CO and NH_3/N_2 balance is needed. Improved trigonometric parallaxes would improve the analysis. Despite the uncertainties and discrepancies, the models reproduce the observed near-infrared spectra well. We find that for the Y0, W1738, T_eff = 425 +/- 25 K and log g = 4.0 +/- 0.25, and for the Y1, W0350, T_eff = 350 +/- 25 K and log g = 4.0 +/- 0.25. W1738 may be metal-rich. Based on evolutionary models, these temperatures and gravities correspond to a mass range for both Y dwarfs of 3-9 Jupiter masses, with W0350 being a cooler, slightly older, version of W1738; the age of W0350 is 0.3-3 Gyr, and the age of W1738 is 0.15-1 Gyr.

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

We present results from an atmospheric circulation study of nine hot Jupiters that comprise a large transmission spectral survey using the Hubble and Spitzer Space Telescopes. These observations exhibit a range of spectral behavior over optical and infrared wavelengths which suggest diverse cloud and haze properties in their atmospheres. By utilizing the specific system parameters for each planet, we naturally probe a wide phase space in planet radius, gravity, orbital period, and equilibrium temperature. First, we show that our model "grid" recovers trends shown in traditional parametric studies of hot Jupiters, particularly equatorial superrotation and increased day-night temperature contrast with increasing equilibrium temperature. We show how spatial temperature variations, particularly between the dayside and nightside and west and east terminators, can vary by hundreds of K, which could imply large variations in Na, K, CO and CH4 abundances in those regions. These chemical variations can be large enough to be observed in transmission with high-resolution spectrographs, such as ESPRESSO on VLT, METIS on the E-ELT, or with MIRI and NIRSpec aboard JWST. We also compare theoretical emission spectra generated from our models to available Spitzer eclipse depths for each planet, and find that the outputs from our solar-metallicity, cloud-free models generally provide a good match to many of the datasets, even without additional model tuning. Although these models are cloud-free, we can use their results to understand the chemistry and dynamics that drive cloud formation in their atmospheres.

We present the first detection of photometric variability of a spectroscopically-confirmed Y dwarf. The Infrared Array Camera on board the Spitzer Space Telescope was used to obtain times series photometry at 3.6 and 4.5 microns over a twenty four hour period at two different epochs separated by 149 days. Variability is evident at 4.5 um in the first epoch and at 3.6 and 4.5 um in the second epoch which suggests that the underlying cause or causes of this variability change on the timescales of months. The second-epoch [3.6] and [4.5] light curves are nearly sinusoidal in form, in phase, have periods of roughly 8.5 hours, and have semi-amplitudes of 3.5%. We find that a simple geometric spot model with a single bright spot reproduces these observations well. We also compare our measured semi-amplitudes of the second epoch light curves to predictions of the static, one-dimensional, partly cloudy and hot spot models of Morley and collaborators and find that neither set of models can reproduce the observed [3.6] and[4.5] semi-amplitudes simultaneously. More advanced two- or three-dimensional models that include time-dependent phenomena like vertical mixing, cloud formation, and thermal relaxation are therefore sorely needed in order to properly interpret our observations.

Detecting polarized light from self-luminous exoplanets has the potential to provide key information about rotation, surface gravity, cloud grain size, and cloud coverage. While field brown dwarfs with detected polarized emission are common, no exoplanet or substellar companion has yet been detected in polarized light. With the advent of high contrast imaging spectro-polarimeters such as GPI and SPHERE, such a detection may now be possible with careful treatment of instrumental polarization. In this paper, we present 28 minutes of $H$-band GPI polarimetric observations of the benchmark T5.5 companion HD 19467 B. We detect no polarization signal from the target, and place an upper limit on the degree of linear polarization of $p_{\text{CL}99.73\%} \leq 2.4\%$. We discuss our results in the context of T dwarf cloud models and photometric variability.

Observations of Neptune with the Kepler Space Telescope yield a 49-day light curve with 98% coverage at a 1-minute cadence. A significant signature in the light curve comes from discrete cloud features. We compare results extracted from the light curve data with contemporaneous disk-resolved imaging of Neptune from the Keck 10-meter telescope at 1.65 microns and Hubble Space Telescope visible imaging acquired 9 months later. This direct comparison validates the feature latitudes assigned to the K2 light curve periods based on Neptune's zonal wind profile, and confirms observed cloud feature variability. Although Neptune's clouds vary in location and intensity on short and long time scales, a single large discrete storm seen in Keck imaging dominates the K2 and Hubble light curves; smaller or fainter clouds likely contribute to short-term brightness variability. The K2 Neptune light curve, in conjunction with our imaging data, provides context for the interpretation of current and future brown dwarf and extrasolar planet variability measurements. In particular we suggest that the balance between large, relatively stable, atmospheric features and smaller, more transient, clouds controls the character of substellar atmospheric variability. Atmospheres dominated by a few large spots may show inherently greater light curve stability than those which exhibit a greater number of smaller features.

Rotational modulations of brown dwarfs have recently provided powerful constraints on the properties of ultra-cool atmospheres, including longitudinal and vertical cloud structures and cloud evolution. Furthermore, periodic light curves directly probe the rotational periods of ultra-cool objects. We present here, for the first time, time-resolved high-precision photometric measurements of a planetary-mass companion, 2M1207b. We observed the binary system with HST/WFC3 in two bands and with two spacecraft roll angles. Using point spread function-based photometry, we reach a nearly photon-noise limited accuracy for both the primary and the secondary. While the primary is consistent with a flat light curve, the secondary shows modulations that are clearly detected in the combined light curve as well as in different subsets of the data. The amplitudes are 1.36% in the F125W and 0.78% in the F160W filters, respectively. By fitting sine waves to the light curves, we find a consistent period of $10.7^{+1.2}_{-0.6}$ hours and similar phases in both bands. The J- and H-band amplitude ratio of 2M1207b is very similar to a field brown dwarf that has identical spectral type but different J-H color. Importantly, our study also measures, for the first time, the rotation period for a directly imaged extra-solar planetary-mass companion.

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

Several concepts now exist for small, space-based missions to directly characterize exoplanets in reflected light. Here, we develop an instrument noise model suitable for studying the spectral characterization potential of a coronagraph-equipped, space-based telescope. We adopt a baseline set of telescope and instrument parameters appropriate for near-future planned missions like WFIRST-AFTA, including a 2 m diameter primary aperture, an operational wavelength range of 0.4-1.0 um, and an instrument spectral resolution of 70, and apply our baseline model to a variety of spectral models of different planet types, including Earth twins, Jupiter twins, and warm and cool Jupiters and Neptunes. With our exoplanet spectral models, we explore wavelength-dependent planet-star flux ratios for main sequence stars of various effective temperatures, and discuss how coronagraph inner and outer working angle constraints will influence the potential to study different types of planets. For planets most favorable to spectroscopic characterization, we study the integration times required to achieve moderate signal-to-noise ratio spectra. We also explore the sensitivity of the integration times required to either detect the bottom or presence of key absorption bands to coronagraph raw contrast performance, exozodiacal light levels, and the distance to the planetary system. Decreasing detector quantum efficiency at longer visible wavelengths makes the detection of water vapor in the atmospheres of Earth-like planets extremely challenging, and also hinders detections of the 0.89 um methane band. Additionally, most modeled observations have noise dominated by dark current, indicating that improving CCD performance could substantially drive down requisite integration times. Finally, we briefly discuss the extension of our models to a more distant future Large UV-Optical-InfraRed (LUVOIR) mission.

Planets larger than Earth and smaller than Neptune are some of the most numerous in the galaxy, but observational efforts to understand this population have proved challenging because optically thick clouds or hazes at high altitudes obscure molecular features (Kreidberg et al. 2014b). We present models of super Earths that include thick clouds and hazes and predict their transmission, thermal emission, and reflected light spectra. Very thick, lofted clouds of salts or sulfides in high metallicity (1000x solar) atmospheres create featureless transmission spectra in the near-infrared. Photochemical hazes with a range of particle sizes also create featureless transmission spectra at lower metallicities. Cloudy thermal emission spectra have muted features more like blackbodies, and hazy thermal emission spectra have emission features caused by an inversion layer at altitudes where the haze forms. Close analysis of reflected light from warm (~400-800 K) planets can distinguish cloudy spectra, which have moderate albedos (0.05-0.20), from hazy models, which are very dark (0.0-0.03). Reflected light spectra of cold planets (~200 K) accessible to a space-based visible light coronagraph will have high albedos and large molecular features that will allow them to be more easily characterized than the warmer transiting planets. We suggest a number of complementary observations to characterize this population of planets, including transmission spectra of hot (>1000 K) targets, thermal emission spectra of warm targets using the James Webb Space Telescope (JWST), high spectral resolution (R~10^5) observations of cloudy targets, and reflected light spectral observations of directly-imaged cold targets. Despite the dearth of features observed in super Earth transmission spectra to date, different observations will provide rich diagnostics of their atmospheres.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

The re-emergence of the 0.99 $\mu$m FeH feature in brown dwarfs of early- to mid-T spectral type has been suggested as evidence for cloud disruption where flux from deep, hot regions below the Fe cloud deck can emerge. The same mechanism could account for color changes at the L/T transition and photometric variability. We present the first observations of spectroscopic variability of brown dwarfs covering the 0.99 $\mu$m FeH feature. We observed the spatially resolved very nearby brown dwarf binary WISE J104915.57-531906.1 (Luhman 16AB), a late-L and early-T dwarf, with HST/WFC3 in the G102 grism at 0.8-1.15 $\mu$m. We find significant variability at all wavelengths for both brown dwarfs, with peak-to-valley amplitudes of 9.3% for Luhman 16B and 4.5% for Luhman 16A. This represents the first unambiguous detection of variability in Luhman 16A. We estimate a rotational period between 4.5 and 5.5 h, very similar to Luhman 16B. Variability in both components complicates the interpretation of spatially unresolved observations. The probability for finding large amplitude variability in any two brown dwarfs is less than 10%. Our finding may suggest that a common but yet unknown feature of the binary is important for the occurrence of variability. For both objects, the amplitude is nearly constant at all wavelengths except in the deep K I feature below 0.84 $\mu$m. No variations are seen across the 0.99 $\mu$m FeH feature. The observations lend strong further support to cloud height variations rather than holes in the silicate clouds, but cannot fully rule out holes in the iron clouds. We re-evaluate the diagnostic potential of the FeH feature as a tracer of cloud patchiness.

We present $H$-band observations of $\beta$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $\beta$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($\phi=85.27{\deg}^{+0.26}_{-0.19}$) with a position angle $\theta_{PA}=30.35{\deg}^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $\theta_{PA}\approx29{\deg}$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $\beta$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $\Omega=31.75{\deg}\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $\beta$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $\beta$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $\beta$ Pic b; 2) the mutual inclination of the inner disk and $\beta$ Pic $b$ is $4{\deg}$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $\beta$ Pic $b$ or self-stirring by the disk's parent bodies.

Interpreting the spectra of brown dwarfs is key to determining the fundamental physical and chemical processes occurring in their atmospheres. Powerful Bayesian atmospheric retrieval tools have recently been applied to both exoplanet and brown dwarf spectra to tease out the thermal structures and molecular abundances to understand those processes. In this manuscript we develop a significantly upgraded retrieval method and apply it to the SpeX spectral library data of two benchmark late T-dwarfs, Gl570D and HD3651B, to establish the validity of our upgraded forward model parameterization and Bayesian estimator. Our retrieved metallicities, gravities, and effective temperature are consistent with the metallicity and presumed ages of the systems. We add the carbon-to-oxygen ratio as a new dimension to benchmark systems and find good agreement between carbon-to-oxygens ratio derived in the brown dwarfs and the host stars. Furthermore, we have for the first time unambiguously determined the presence of ammonia in the low-resolution spectra of these two late T-dwarfs. We also show that the retrieved results are not significantly impacted by the possible presence of clouds, though some quantities are significantly impacted by uncertainties in photometry. This investigation represents a watershed study in establishing the utility of atmospheric retrieval approaches on brown dwarf spectra.

We use a planetary albedo model to investigate variations in visible wavelength phase curves of exoplanets. The presence of clouds on these exoplanets significantly alters their planetary albedo spectra. We confirm that non-uniform cloud coverage on the dayside of tidally locked exoplanets will manifest as changes to the magnitude and shift of the phase curve. In this work, we first investigate a test case of our model using a Jupiter-like planet, at temperatures consistent to 2.0 AU insolation from a solar type star, to consider the effect of H2O clouds. We then extend our application of the model to the exoplanet Kepler-7b and consider the effect of varying cloud species, sedimentation efficiency, particle size, and cloud altitude. We show that, depending on the observational filter, the largest possible shift of the phase curve maximum will be 2-10 deg for a Jupiter-like planet, and up to 30 deg (0.08 in fractional orbital phase) for hot-Jupiter exoplanets at visible wavelengths as a function of dayside cloud distribution with a uniformly averaged thermal profile. Finally, we tailor our model for comparison with, and confirmation of, the recent optical phase-curve observations of Kepler-7b with the Kepler space telescope. The average planetary albedo can vary between 0.1-0.6 for the 1300 cloud scenarios that were compared to the observations. We observe that smaller particle size and increasing cloud altitude have a strong effect on increasing albedo. In particular, we show that a set of models where Kepler-7b has roughly half of its dayside covered in small-particle clouds high in the atmosphere, made of bright minerals like MgSiO3 and Mg2SiO4, provide the best fits to the observed offset and magnitude of the phase-curve, whereas Fe clouds are found to have too dark to fit the observations.

We explored two aspects of the problem of characterizing cool extrasolar giant planets in scattered optical light with a space based coronagraph. First, for a number of the known radial velocity (RV) giants we computed traditional forward models of their atmospheric structure and clouds, given various input assumptions, and computed model albedo spectra. Such models have been computed before, but mostly for generic planets. Our new models demonstrate that there is likely interesting spectral diversity among those planets that are most favorable for direct detection. Second, we applied a powerful Markov Chain Monte Carlo (MCMC) retrieval technique to synthetic noisy data of cool giants to better understand how well various atmospheric parameters--particularly molecular abundances and cloud properties--could be constrained. This is the first time such techniques have been applied to this problem. The process is time consuming, so only a dozen or so cases could be completed in the limited time available. Nevertheless the results clearly show that even at S/N ~ 5, scientifically interesting and valuable conclusions can be drawn about the properties of giant planet atmospheres from noisy spectra. We find that atmospheric abundances are best constrained when the planet gravity is bounded. Thus direct imaging observations of known radial velocity planets are extremely valuable as limits on a target planet's gravity can be obtained from astrometric imaging and reflectivity and mass-radius relationship arguments (for mass and radius, respectively). Further retrieval studies are clearly warranted and would be valuable to help guide instrument design decisions.

WD0137-349 is a white dwarf-brown dwarf binary system in a 116 minute orbit. We present radial velocity observations and multiwaveband photometry from V, R and I in the optical, to J, H and Ks in the near-IR and [3.6], [4.5], [5.8] and [8.0] microns in the mid-IR. The photometry and lightcurves show variability in all wavebands, with the amplitude peaking at [4.5] microns, where the system is also brightest. Fluxes and brightness temperatures were computed for the heated and unheated atmosphere of the brown dwarf (WD0137-349B) using synthetic spectra of the white dwarf using model atmosphere simulations. We show that the flux from the brown dwarf dayside is brighter than expected in the Ks and [4.5] micron bands when compared to models of irradiated brown dwarfs with full energy circulation and suggest this over-luminosity may be attributed to H2 fluorescence or H3+ being generated in the atmosphere by the UV irradiation.

We present new near-infrared photometry for seven late-type T dwarfs and nine Y-type dwarfs, and lower limit magnitudes for a tenth Y dwarf, obtained at Gemini Observatory. We also present a reanalysis of H-band imaging data from the Keck Observatory Archive, for an eleventh Y dwarf. These data are combined with earlier MKO-system photometry, Spitzer and WISE mid-infrared photometry, and available trigonometric parallaxes, to create a sample of late-type brown dwarfs which includes ten T9-T9.5 dwarfs or dwarf systems, and sixteen Y dwarfs. We compare the data to our models which include updated H_2 and NH_3 opacity, as well as low-temperature condensate clouds. The models qualitatively reproduce the trends seen in the observed colors, however there are discrepancies of around a factor of two in flux for the Y0-Y1 dwarfs, with T_eff~350-400K. At T_eff~400K, the problems could be addressed by significantly reducing the NH_3 absorption, for example by halving the abundance of NH_3 possibly by vertical mixing. At T_eff~350K, the discrepancy may be resolved by incorporating thick water clouds. The onset of these clouds might occur over a narrow range in T_eff, as indicated by the observed small change in 5um flux over a large change in J-W2 color. Of the known Y dwarfs, the reddest in J-W2 are WISEP J182831.08+265037.8 and WISE J085510.83-071442.5. We interpret the former as a pair of identical 300-350K dwarfs, and the latter as a 250K dwarf. If these objects are ~3 Gyrs old, their masses are ~10 and ~5 Jupiter-masses respectively.

We present results from the "Weather on Other Worlds" Spitzer Exploration Science program to investigate photometric variability in L and T dwarfs, usually attributed to patchy clouds. We surveyed 44 L3-T8 dwarfs, spanning a range of $J-K_s$ colors and surface gravities. We find that 14/23 (61%; 95% confidence interval: 41%-78%) of our single L3-L9.5 dwarfs are variable with peak-to-peak amplitudes between 0.2% and 1.5%, and 5/16 (31%; 95% confidence interval: 14%-56%) of our single T0-T8 dwarfs are variable with amplitudes between 0.8% and 4.6%. After correcting for sensitivity, we find that 80% (95% confidence interval: 53%-100%) of L dwarfs vary by >0.2%, and 36% (95% confidence interval: 19%-52%) of T dwarfs vary by >0.4%. Given viewing geometry considerations, we conclude that photospheric heterogeneities causing >0.2% 3-5-micron flux variations are present on virtually all L dwarfs, and probably on most T dwarfs. A third of L dwarf variables show irregular light curves, indicating that L dwarfs may have multiple spots that evolve over a single rotation. Also, approximately a third of the periodicities are on time scales >10 h, suggesting that slowly-rotating brown dwarfs may be common. We observe an increase in the maximum amplitudes over the entire spectral type range, revealing a potential for greater temperature contrasts in T dwarfs than in L dwarfs. We find a tentative association (92% confidence) between low surface gravity and high-amplitude variability among L3-L5.5 dwarfs. Although we can not confirm whether lower gravity is also correlated with a higher incidence of variables, the result is promising for the characterization of directly imaged young extrasolar planets through variability.

We present time-resolved near-infrared spectroscopy of two L5 dwarfs, 2MASS J18212815+1414010 and 2MASS J15074759-1627386, observed with the Wide Field Camera 3 instrument on the Hubble Space Telescope (HST). We study the wavelength dependence of rotation-modulated flux variations between 1.1 $\mu$m and 1.7 $\mu$m. We find that the water absorption bands of the two L5 dwarfs at 1.15 $\mu$m and 1.4 $\mu$m vary at similar amplitudes as the adjacent continuum. This differs from the results of previous HST observations of L/T transition dwarfs, in which the water absorption at 1.4 $\mu$m displays variations of about half of the amplitude at other wavelengths. We find that the relative amplitude of flux variability out of the water band with respect to that in the water band shows a increasing trend from the L5 dwarfs toward the early T dwarfs. We utilize the models of Saumon & Marley (2008) and find that the observed variability of the L5 dwarfs can be explained by the presence of spatially varying high-altitude haze layers above the condensate clouds. Therefore, our observations show that the heterogeneity of haze layers - the driver of the variability - must be located at very low pressures, where even the water opacity is negligible. In the near future, the rotational spectral mapping technique could be utilized for other atomic and molecular species to probe different pressure levels in the atmospheres of brown dwarfs and exoplanets and uncover both horizontal and vertical cloud structures.

The recent discovery and characterization of the diversity of the atmospheres of exoplanets and brown dwarfs calls for the development of fast and accurate analytical models. We quantify the accuracy of the analytical solution derived in paper I for an irradiated, non-grey atmosphere by comparing it to a state-of-the-art radiative transfer model. Then, using a grid of numerical models, we calibrate the different coefficients of our analytical model for irradiated solar-composition atmospheres of giant exoplanets and brown dwarfs. We show that the so-called Eddington approximation used to solve the angular dependency of the radiation field leads to relative errors of up to 5% on the temperature profile. We show that for realistic non-grey planetary atmospheres, the presence of a convective zone that extends to optical depths smaller than unity can lead to changes in the radiative temperature profile on the order of 20% or more. When the convective zone is located at deeper levels (such as for strongly irradiated hot Jupiters), its effect on the radiative atmosphere is smaller. We show that the temperature inversion induced by a strong absorber in the optical, such as TiO or VO is mainly due to non-grey thermal effects reducing the ability of the upper atmosphere to cool down rather than an enhanced absorption of the stellar light as previously thought. Finally, we provide a functional form for the coefficients of our analytical model for solar-composition giant exoplanets and brown dwarfs. This leads to fully analytical pressure-temperature profiles for irradiated atmospheres with a relative accuracy better than 10% for gravities between 2.5m/s^2 and 250 m/s^2 and effective temperatures between 100 K and 3000 K. This is a great improvement over the commonly used Eddington boundary condition.

The binary brown dwarf WISE J104915.57$-$531906.1 (also Luhman 16AB), composed of a late L and early T dwarf, is a prototypical L/T transition flux reversal binary located at only 2 pc distance. Luhman 16B is a known variable whose light curves evolve rapidly. We present spatially resolved spectroscopic time-series of Luhman 16A and B covering 6.5 h using HST/WFC3 at 1.1 to 1.66 $\mu$m. The small, count-dependent variability of Luhman 16A at the beginning of the observations likely stems from instrumental systematics; Luhman 16A appears non-variable above $\approx$0.4%. Its spectrum is well fit by a single cloud layer with intermediate cloud thickness (f_sed=2, Teff=1200 K). Luhman 16B varies at all wavelengths with peak-to-valley amplitudes of 7-11%. The amplitude and light curve shape changes over only one rotation period. The lowest relative amplitude is found in the deep water absorption band at 1.4 $\mu$m, otherwise it mostly decreases gradually from the blue to the red edge of the spectrum. This is very similar to the other two known highly variable early T dwarfs. A two-component cloud model accounts for most of the variability, although small deviations are seen in the water absorption band. We fit the mean spectrum and relative amplitudes with a linear combination of two models of a warm, thinner cloud (Teff=1300 K, fsed=3) and a cooler, thicker cloud (Teff=1000-1100 K, f_sed=1), assuming out-of-equilibrium atmospheric chemistry. A cloud as for Luhman 16A but with holes cannot reproduce the variability of Luhman 16B, indicating more complex cloud evolution through the L/T transition. The projected separation of the binary has decreased by $\approx$0.3'' in 8 months.

The atmosphere of a brown dwarf or extrasolar giant planet controls the spectrum of radiation emitted by the object and regulates its cooling over time. While the study of these atmospheres has been informed by decades of experience modeling stellar and planetary atmospheres, the distinctive characteristics of these objects present unique challenges to forward modeling. In particular, complex chemistry arising from molecule-rich atmospheres, molecular opacity line lists (sometimes running to 10 billion absorption lines or more) multiple cloud-forming condensates, and disequilibrium chemical processes all combine to create a challenging task for any modeling effort. This review describes the process of incorporating these complexities into one-dimensional radiative-convective equilibrium models of sub-stellar objects. We discuss the underlying mathematics as well as the techniques used to model the physics, chemistry, radiative transfer, and other processes relevant to understanding these atmospheres. The review focuses on the process of the creation of atmosphere models and briefly presents some comparisons of model predictions to data. Current challenges in the field and some comments on the future conclude the review.

We address disequilibrum abundances of some simple molecules in the atmospheres of solar composition brown dwarfs and self-luminous extrasolar giant planets using a kinetics-based 1D atmospheric chemistry model. Our approach is to use the full kinetics model to survey the parameter space with effective temperatures between 500 K and 1100 K. In all of these worlds equilibrium chemistry favors CH4 over CO in the parts of the atmosphere that can be seen from Earth, but in most disequilibrium favors CO. The small surface gravity of a planet strongly discriminates against CH4 when compared to an otherwise comparable brown dwarf. If vertical mixing is like Jupiter's, the transition from methane to CO occurs at 500 K in a planet. Sluggish vertical mixing can raise this to 600 K; but clouds or more vigorous vertical mixing could lower this to 400 K. The comparable thresholds in brown dwarfs are $1100\pm100$ K. Ammonia is also sensitive to gravity, but unlike CH4/CO, the NH3/N2 ratio is insensitive to mixing, which makes NH3 a potential proxy for gravity. HCN may become interesting in high gravity brown dwarfs with very strong vertical mixing. Detailed analysis of the CO-CH4 reaction network reveals that the bottleneck to CO hydrogenation goes through methanol, in partial agreement with previous work. Simple, easy to use quenching relations are derived by fitting to the complete chemistry of the full ensemble of models. These relations are valid for determining CO, CH4, NH3, HCN, and CO2 abundances in the range of self-luminous worlds we have studied but may not apply if atmospheres are strongly heated at high altitudes by processes not considered here (e.g., wave breaking).

High altitude clouds and hazes are integral to understanding exoplanet observations, and are proposed to explain observed featureless transit spectra. However, it is difficult to make inferences from these data because of the need to disentangle effects of gas absorption from haze extinction. Here, we turn to the quintessential hazy world -- Titan -- to clarify how high altitude hazes influence transit spectra. We use solar occultation observations of Titan's atmosphere from the Visual and Infrared Mapping Spectrometer (VIMS) aboard NASA's Cassini spacecraft to generate transit spectra. Data span 0.88-5 microns at a resolution of 12-18 nm, with uncertainties typically smaller than 1%. Our approach exploits symmetry between occultations and transits, producing transit radius spectra that inherently include the effects of haze multiple scattering, refraction, and gas absorption. We use a simple model of haze extinction to explore how Titan's haze affects its transit spectrum. Our spectra show strong methane absorption features, and weaker features due to other gases. Most importantly, the data demonstrate that high altitude hazes can severely limit the atmospheric depths probed by transit spectra, bounding observations to pressures smaller than 0.1-10 mbar, depending on wavelength. Unlike the usual assumption made when modeling and interpreting transit observations of potentially hazy worlds, the slope set by haze in our spectra is not flat, and creates a variation in transit height whose magnitude is comparable to those from the strongest gaseous absorption features. These findings have important consequences for interpreting future exoplanet observations, including those from NASA's James Webb Space Telescope.

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

The hot-Jupiter HAT-P-2b has become a prime target for Spitzer Space Telescope observations aimed at understanding the atmospheric response of exoplanets on highly eccentric orbits. Here we present a suite of three-dimensional atmospheric circulation models for HAT-P-2b that investigate the effects of assumed atmospheric composition and rotation rate on global scale winds and thermal patterns. We compare and contrast atmospheric models for HAT-P-2b, which assume one and five times solar metallicity, both with and without TiO/VO as atmospheric constituents. Additionally we compare models that assume a rotation period of half, one, and two times the nominal pseudo-synchronous rotation period. We find that changes in assumed atmospheric metallicity and rotation rate do not significantly affect model predictions of the planetary flux as a function of orbital phase. However, models in which TiO/VO are present in the atmosphere develop a transient temperature inversion between the transit and secondary eclipse events that results in significant variations in the timing and magnitude of the peak of the planetary flux compared with models in which TiO/VO are omitted from the opacity tables. We find that no one single atmospheric model can reproduce the recently observed full orbit phase curves at 3.6, 4.5 and 8.0 microns, which is likely due to a chemical process not captured by our current atmospheric models for HAT-P-2b. Further modeling and observational efforts focused on understanding the chemistry of HAT-P-2b's atmosphere are needed and could provide key insights into the interplay between radiative, dynamical, and chemical processes in a wide range of exoplanet atmospheres.

We present new calculations of Rosseland and Planck gaseous mean opacities relevant to the atmospheres of giant planets and ultracool dwarfs. Such calculations are used in modeling the atmospheres, interiors, formation, and evolution of these objects. Our calculations are an expansion of those presented in Freedman et al. (2008) to include lower pressures, finer temperature resolution, and also the higher metallicities most relevant for giant planet atmospheres. Calculations span 1 microbar to 300 bar, and 75 K to 4000 K, in a nearly square grid. Opacities at metallicities from solar to 50 times solar abundances are calculated. We also provide an analytic fit to the Rosseland mean opacities over the grid in pressure, temperature, and metallicity. In addition to computing mean opacities at these local temperatures, we also calculate them with weighting functions up to 7000 K, to simulate the mean opacities for incident stellar intensities, rather than locally thermally emitted intensities. The chemical equilibrium calculations account for the settling of condensates in a gravitational field and are applicable to cloud-free giant planet and ultracool dwarf atmospheres, but not circumstellar disks. We provide our extensive opacity tables for public use.

We present a new Y dwarf, WISE J030449.03-270508.3, confirmed from a candidate sample designed to pick out low temperature objects from the WISE database. The new object is typed Y0pec following a visual comparison with spectral standards, and lies at a likely distance of 10-17 pc. Its tangential velocity suggests thin disk membership, but it shows some spectral characteristics that suggest it may be metal-poor and/or older than previously identified Y0 dwarfs. Based on trends seen for warmer late type T dwarfs, the Y-band flux peak morphology is indicative of sub-solar metallicity, and the enhanced red wing of the J-band flux peak offers evidence for high gravity and/or low metallicity (with associated model trends suggesting an age closer to ~10 Gyr and mass in the range 0.02-0.03 Mo). This object may thus be extending the population parameter-space of the known Y0 dwarfs.

Brown dwarf spectra contain a wealth of information about their molecular abundances, temperature structure, and gravity. We present a new data driven retrieval approach, previously used in planetary atmosphere studies, to extract the molecular abundances and temperature structure from brown dwarf spectra. The feasibility of the approach is first demonstrated on a synthetic brown dwarf spectrum. Given typical spectral resolutions, wavelength coverage, and noise properties precisions of tens of percent can be obtained for the molecular abundances and 10s-100s K on the temperature profile. The technique is then applied to the well studied brown dwarf, Gl 570D. From this spectral retrieval the spectroscopic radius is constrained to be 0.75 - 0.83 $R_{\mathrm{J}}$, $\log(g)$ to be 5.13 - 5.46 and $T_{\mathrm{eff}}$ to be between 804 and 849 K. Estimates for the range of abundances and allowed temperature profiles are also derived. The results from our retrieval approach are in agreement with the self-consistent grid modeling results of Saumon et al (2006). This new approach will allow us to address issues of compositional differences between brown dwarfs and possibly their formation environments, disequilibrium chemistry, missing physics in current grid modeling approaches as well as a many other issues.

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

We present an analytical model for the transmission spectrum of a transiting exoplanet, showing that a cloud base can produce an observable inflection point in the spectrum. The wavelength and magnitude of the inflection can be used to break the degeneracy between the atmospheric pressure and the abundance of the main cloud material, however, the abundance still depends on cloud particle size. An observed inflection also provides a specific point on the atmospheric P-T profile, giving us a "thermometer" to directly validate or rule out postulated cloud species. We apply the model to the transit spectrum of HD 189733b.

Gas-giant planets emit a large fraction of their light in the mid-infrared ($\gtrsim$3$\mu$m), where photometry and spectroscopy are critical to our understanding of the bulk properties of extrasolar planets. Of particular importance are the L and M-band atmospheric windows (3-5$\mu$m), which are the longest wavelengths currently accessible to ground-based, high-contrast imagers. We present binocular LBT AO images of the HR 8799 planetary system in six narrow-band filters from 3-4$\mu$m, and a Magellan AO image of the 2M1207 planetary system in a broader 3.3$\mu$m band. These systems encompass the five known exoplanets with luminosities consistent with L$\rightarrow$T transition brown dwarfs. Our results show that the exoplanets are brighter and have shallower spectral slopes than equivalent temperature brown dwarfs in a wavelength range that contains the methane fundamental absorption feature (spanned by the narrowband filters and encompassed by the broader 3.3$\mu$m filter). For 2M1207 b, we find that thick clouds and non-equilibrium chemistry caused by vertical mixing can explain the object's appearance. For the HR 8799 planets, we present new models that suggest the atmospheres must have patchy clouds, along with non-equilibrium chemistry. Together, the presence of a heterogeneous surface and vertical mixing presents a picture of dynamic planetary atmospheres in which both horizontal and vertical motions influence the chemical and condensate profiles.

Brown dwarfs of a variety of spectral types have been observed to be photometrically variable. Previous studies have focused on objects at the L/T transition, where the iron and silicate clouds in L dwarfs break up or dissipate. However, objects outside of this transitional effective temperature regime also exhibit variability. Here, we present models for mid-late T dwarfs and Y dwarfs. We present models that include patchy salt and sulfide clouds as well as water clouds for the Y dwarfs. We find that for objects over 375 K, patchy cloud opacity would generate the largest amplitude variability within near-infrared spectral windows. For objects under 375 K, water clouds also become important and generate larger amplitude variability in the mid-infrared. We also present models in which we perturb the temperature structure at different pressure levels of the atmosphere to simulate hot spots. These models show the most variability in the absorption features between spectral windows. The variability is strongest at wavelengths that probe pressure levels at which the heating is the strongest. The most illustrative types of observations for understanding the physical processes underlying brown dwarf variability are simultaneous, multi-wavelength observations that probe both inside and outside of molecular absorption features.

Recent discoveries of variable brown dwarfs have provided us with a new window into their three-dimensional cloud structure. The highest variables are found at the L/T transition, where the cloud cover is thought to break up, but variability has been found to occur also for both cloudy L dwarfs and (mostly) cloud-free mid T dwarfs. We summarize results from recent HST programs measuring the spectral variability of brown dwarfs in the near-infrared and compare to results from ground-based programs. We discuss the patchy cloud structure of L/T transition objects, for which it is becoming increasingly certain that the variability does not arise from cloud holes into the deep hot regions but from varying cloud thickness. We present a new patchy cloud model to explain the spectral variability of 2MASSJ21392676+0220226. We also discuss the curious multi-wavelength variability behavior of the recently discovered very nearby early T dwarf WISE J104915.57-531906.1B (Luhman 16B) and the mid T dwarf 2MASS J22282889-431026.

A number of brown dwarfs are now known to be variable with observed amplitudes as large as 10-30% at some wavelengths. While spatial inhomogeneities in cloud coverage and thickness are likely responsible for much of the observed variability, it is possible that some of the variations arise from atmospheric temperature fluctuations instead of, or in addition to, clouds. To better understand the role that thermal variability might play we present a case study of brown dwarf variability using a newly-developed one-dimensional, time-stepping model of atmospheric thermal structure. We focus on the effects of thermal perturbations, intentionally simplifying the problem through omission of clouds and atmospheric circulation. Model results demonstrate that thermal perturbations occurring deep in the atmosphere (at pressures greater than 10 bar) of a model T-dwarf can be communicated to the upper atmosphere through radiative heating via the windows in near-infrared water opacity. The response time depends on where in the atmosphere a thermal perturbation is introduced. We show that, for certain periodic perturbations, the emission spectrum can have complex, time- and wavelength-dependent behaviors, including phase shifts in times of maximum flux observed at different wavelengths. Since different wavelengths probe different levels in the atmosphere, these variations track a wavelength-dependent set of radiative exchanges happening between different atmospheric levels as a perturbation evolves in time. We conclude that thermal--as well as cloud--fluctuations must be considered as possible contributors to the observed brown dwarf variability.

The formation of clouds affects brown dwarf and planetary atmospheres of nearly all effective temperatures. Iron and silicate condense in L dwarf atmospheres and dissipate at the L/T transition. Minor species such as sulfides and salts condense in mid-late T dwarfs. For brown dwarfs below Teff=450 K, water condenses in the upper atmosphere to form ice clouds. Currently over a dozen objects in this temperature range have been discovered, and few previous theoretical studies have addressed the effect of water clouds on brown dwarf or exoplanetary spectra. Here we present a new grid of models that include the effect of water cloud opacity. We find that they become optically thick in objects below Teff=350-375 K. Unlike refractory cloud materials, water ice particles are significantly non-gray absorbers; they predominantly scatter at optical wavelengths through J band and absorb in the infrared with prominent features, the strongest of which is at 2.8 microns. H2O, NH3, CH4, and H2 CIA are dominant opacity sources; less abundant species such as may also be detectable, including the alkalis, H2S, and PH3. PH3, which has been detected in Jupiter, is expected to have a strong signature in the mid-infrared at 4.3 microns in Y dwarfs around Teff=450 K; if disequilibrium chemistry increases the abundance of PH3, it may be detectable over a wider effective temperature range than models predict. We show results incorporating disequilibrium nitrogen and carbon chemistry and predict signatures of low gravity in planetary- mass objects. Lastly, we make predictions for the observability of Y dwarfs and planets with existing and future instruments including the James Webb Space Telescope and Gemini Planet Imager.

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean-molecular weight (i.e., H2-dominated) and a high mean-molecular weight (i.e. water- and CO2-dominated). We find that atmospheres with a low mean-molecular weight have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day-night temperature variations and hence stronger superrotation. In comparison, atmospheres with a high mean-molecular weight have larger day-night and equator-to-pole temperature variations than low mean-molecular weight atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO2-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and lightcurves for 50x solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50x solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations.

We have determined astrometric positions for 15 WISE-discovered late-type brown dwarfs (6 T8-9 and 9 Y dwarfs) using the Keck II telescope, the Spitzer Space Telescope, and the Hubble Space Telescope. Combining data from 8 to 20 epochs we derive parallactic and proper motions for these objects which put the majority within 15 parsecs. For ages greater than a few Gyr, as suggested from kinematic considerations, we find masses of 10-30 M$_{Jup}$ based on standard models for the evolution of low mass objects with a range of mass estimates for individual objects depending on the model in question. Three of the coolest objects have effective temperatures $\sim$ 350 K and inferred masses of 10-15 M$_{Jup}$. Our parallactic distances confirm earlier photometric estimates (Kirkpatrick et al. 2012) and direct measurements (Marsh et al. 2013, Beichman et al. 2013, Dupuy and Krauss 2013) and suggest that the number of objects with masses below about 15 M$_{Jup}$ must be flat or declining relative to higher mass objects. The masses of the coldest Y dwarfs may be similar to those inferred for recently imaged planet-mass companions to nearby young stars. Objects in this mass range, which appear to be rare in both the interstellar and proto-planetary environments, may both have formed via gravitational fragmentation: the brown dwarfs in interstellar clouds and companion objects in a protoplanetary disk. In both cases, however, the fact that objects in this mass range are relatively infrequent, suggests that this mechanism must be inefficient in both environments.

Some wave features found in the C-ring of Saturn appear to be excited by resonances with normal mode oscillations of the planet. The waves are found at locations in the rings where the ratio of orbital to oscillation frequencies is given by m : m+1 where m is a small integer. I suggest here that it is plausible that ring waves may also be launched at second order resonances, where the frequency ratio would be m : m+2. Indeed, otherwise unassociated wave features are found at such locations in the C-ring. If confirmed the association of planetary modes with additional C-ring wave features would measure additional oscillation frequencies of Saturn and improve the utility of the waves for constraining the internal structure of the planet. Second-order resonances in general do not lie near first order ring resonance locations and thus are not the explanation for the apparent frequency splitting of modes.

It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys.

We present 0.9 - 2.5 um resolved spectra for the ultracool binary WISEPC J121756.91+162640.2AB. The system consists of a pair of brown dwarfs that straddles the currently defined T/Y spectral type boundary. We use synthetic spectra generated by model atmospheres that include chloride and sulfide clouds (Morley et al.), the distance to the system (Dupuy & Kraus), and the radius of each component based on evolutionary models (Saumon & Marley) to determine a probable range of physical properties for the binary. The effective temperature of the T8.5 primary is 550 - 600 K, and that of the Y0 - Y0.5 secondary is 450 K. The atmospheres of both components are either free of clouds or have extremely thin cloud layers. We find that the masses of the primary and secondary are 30 and 22 M_Jup, respectively, and that the age of the system is 4 - 8 Gyr. This age is consistent with astrometric measurements (Dupuy & Kraus) that show that the system has kinematics intermediate between those of the thin and thick disks of the Galaxy. An older age is also consistent with an indication by the H - K colors that the system is slightly metal-poor.

Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planet's atmosphere is vertically compact, perhaps due to a water-rich composition that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted absorption features. Previous models that incorporate clouds have included their effect without a strong physical motivation for their existence. Here, we present model atmospheres of GJ 1214b that include physically-motivated clouds of two types. We model the clouds that form as a result of condensation in chemical equilibrium, as they likely do on brown dwarfs, which include KCl and ZnS for this planet. We also include clouds that form as a result of photochemistry, forming a hydrocarbon haze layer. We use a photochemical kinetics model to understand the vertical distribution and available mass of haze-forming molecules. We model both solar and enhanced-metallicity cloudy models and determine the cloud properties necessary to match observations. In enhanced-metallicity atmospheres, we find that the equilibrium clouds can match the observations of GJ 1214b if they are lofted high into the atmosphere and have a low sedimentation efficiency (fsed=0.1). We find that models with a variety of hydrocarbon haze properties can match the observations. Particle sizes from 0.01 to 0.25 micron can match the transmission spectrum with haze-forming efficiencies as low as 1-5%.

We report the discovery of 76 new T dwarfs from the UKIDSS Large Area Survey (LAS). Near-infrared broad and narrow-band photometry and spectroscopy are presented for the new objects, along with WISE and warm-Spitzer photometry. Proper motions for 128 UKIDSS T dwarfs are presented from a new two epoch LAS proper motion catalogue. We use these motions to identify two new benchmark systems: LHS 6176AB, a T8+M4 pair and HD118865AB, a T5.5+F8 pair. Using age constraints from the primaries and evolutionary models to constrain the radii we have estimated their physical properties from their bolometric luminosity. We compare the colours and properties of known benchmark T dwarfs to the latest model atmospheres and draw two principal conclusions. Firstly, it appears that the H - [4.5] and J - W2 colours are more sensitive to metallicity than has previously been recognised, such that differences in metallicity may dominate over differences in Teff when considering relative properties of cool objects using these colours. Secondly, the previously noted apparent dominance of young objects in the late-T dwarf sample is no longer apparent when using the new model grids and the expanded sample of late-T dwarfs and benchmarks. This is supported by the apparently similar distribution of late-T dwarfs and earlier-type T dwarfs on reduced proper motion diagrams that we present. Finally, we present updated space densities for the late-T dwarfs, and compare our values to simulation predictions and those from WISE.

Photometric monitoring from warm Spitzer reveals that the L3 dwarf DENIS-P J1058.7-1548 varies sinusoidally in brightness with a period of 4.25 +0.26/-0.16 hours and an amplitude of 0.388 +/- 0.043% (peak-to-valley) in the 3.6 micron band, confirming the reality of a 4.31 +/- 0.31 hour periodicity detected in J-band photometry from the SOAR telescope. The J-band variations are a factor of 2.17 +/- 0.35 larger in amplitude than those at 3.6 microns, while 4.5 micron Spitzer observations yield a 4.5 micron/3.6 micron amplitude ratio of only 0.23 +/- 0.15, consistent with zero 4.5 micron variability. This wide range in amplitudes indicates rotationally modulated variability due to magnetic phenomena and/or inhomogeneous cloud cover. Weak H-alpha emission indicates some magnetic activity, but it is difficult to explain the observed amplitudes by magnetic phenomena unless they are combined with cloud inhomogeneities (which might have a magnetic cause). However, inhomogenous cloudcover alone can explain all our observations, and our data align with theory in requiring that the regions with the thickest clouds also have the lowest effective temperature. Combined with published v sin(i) results, our rotation period yields a 95% confidence lower limit of 0.111 solar radii on the radius of DENIS-P J1058.7-1548, suggesting upper limits of 320 Myr and 0.055 solar masses on the age and mass. These limits should be regarded cautiously because of 3-sigma inconsistencies with other data; however, a lower limit of 45 degrees on the inclination is more secure. DENIS-P J1058.7-1548 is only the first of nearly two dozen low-amplitude variables discovered and analyzed by the Weather on Other Worlds project.

Of the over 800 exoplanets detected to date, over half are on non-circular orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable stellar heating, which has implications for the planet's atmospheric dynamical regime. However, little is known about this dynamical regime, and how it may influence observations. Therefore, we present a systematic study of hot Jupiters on highly eccentric orbits using the SPARC/MITgcm, a model which couples a three-dimensional general circulation model with a plane-parallel, two-stream, non-grey radiative transfer model. In our study, we vary the eccentricity and orbit-average stellar flux over a wide range. We demonstrate that the eccentric hot Jupiter regime is qualitatively similar to that of planets on circular orbits; the planets possess a superrotating equatorial jet and exhibit large day-night temperature variations. We show that these day-night heating variations induce momentum fluxes equatorward to maintain the superrotating jet throughout its orbit. As the eccentricity and/or stellar flux is increased, the superrotating jet strengthens and narrows, due to a smaller Rossby deformation radius. For a select number of model integrations, we generate full-orbit lightcurves and find that the timing of transit and secondary eclipse viewed from Earth with respect to periapse and apoapse can greatly affect what we see in infrared (IR) lightcurves; the peak in IR flux can lead or lag secondary eclipse depending on the geometry. For those planets that have large day-night temperature variations and rapid rotation rates, we find that the lightcurves exhibit "ringing" as the planet's hottest region rotates in and out of view from Earth. These results can be used to explain future observations of eccentric transiting exoplanets.

We report the discovery of 87 new T dwarfs uncovered with the Wide-field Infrared Survey Explorer (WISE) and three brown dwarfs with extremely red near-infrared colors that exhibit characteristics of both L and T dwarfs. Two of the new T dwarfs are likely binaries with L7+/-1 primaries and mid-type T secondaries. In addition, our follow-up program has confirmed 10 previously identified T dwarfs and four photometrically-selected L and T dwarf candidates in the literature. This sample, along with the previous WISE discoveries, triples the number of known brown dwarfs with spectral types later than T5. Using the WISE All-Sky Source Catalog we present updated color-color and color-type diagrams for all the WISE-discovered T and Y dwarfs. Near-infrared spectra of the new discoveries are presented, along with spectral classifications. To accommodate later T dwarfs we have modified the integrated flux method of determining spectral indices to instead use the median flux. Furthermore, a newly defined J-narrow index differentiates the early-type Y dwarfs from late-type T dwarfs based on the J-band continuum slope. The K/J indices for this expanded sample show that 32% of late-type T dwarfs have suppressed K-band flux and are blue relative to the spectral standards, while only 11% are redder than the standards. Comparison of the Y/J and K/J index to models suggests diverse atmospheric conditions and supports the possible re-emergence of clouds after the L/T transition. We also discuss peculiar brown dwarfs and candidates that were found not to be substellar, including two Young Stellar Objects and two Active Galactic Nuclei. The coolest WISE-discovered brown dwarfs are the closest of their type and will remain the only sample of their kind for many years to come.

Clouds and hazes are commonplace in the atmospheres of solar system planets and are likely ubiquitous in the atmospheres of extrasolar planets as well. Clouds affect every aspect of a planetary atmosphere, from the transport of radiation, to atmospheric chemistry, to dynamics and they influence - if not control - aspects such as surface temperature and habitability. In this review we aim to provide an introduction to the role and properties of clouds in exoplanetary atmospheres. We consider the role clouds play in influencing the spectra of planets as well as their habitability and detectability. We briefly summarize how clouds are treated in terrestrial climate models and consider the far simpler approaches that have been taken so far to model exoplanet clouds, the evidence for which we also review. Since clouds play a major role in the atmospheres of certain classes of brown dwarfs we briefly discuss brown dwarf cloud modeling as well. We also review how the scattering and extinction efficiencies of cloud particles may be approximated in certain limiting cases of small and large particles in order to facilitate physical understanding. Since clouds play such important roles in planetary atmospheres, cloud modeling may well prove to be the limiting factor in our ability to interpret future observations of extrasolar planets.

We present YJHK photometry, or a subset, for the six Y dwarfs discovered in WISE data by Cushing et al.. The data were obtained using NIRI on the Gemini North telescope. We also present a far-red spectrum obtained using GMOS-North for WISEPC J205628.90+145953.3. We compare the data to Morley et al. (2012) models, which include cloud decks of sulfide and chloride condensates. We find that the models with these previously neglected clouds can reproduce the energy distributions of T9 to Y0 dwarfs quite well, other than near 5um where the models are too bright. This is thought to be because the models do not include departures from chemical equilibrium caused by vertical mixing, which would enhance the abundance of CO, decreasing the flux at 5um. Vertical mixing also decreases the abundance of NH_3, which would otherwise have strong absorption features at 1.03um and 1.52um that are not seen in the Y0 WISEPC J205628.90+145953.3. We find that the five Y0 to Y0.5 dwarfs have 300 < T_eff K < 450, 4.0 < log g < 4.5 and f_sed ~ 3. These temperatures and gravities imply a mass range of 5 - 15 M_Jupiter and ages around 5 Gyr. We suggest that WISEP J182831.08+265037.8 is a binary system, as this better explains its luminosity and color. We find that the data can be made consistent with observed trends, and generally consistent with the models, if the system is composed of a T_eff = 325 K and log g <~ 4.5 primary, and a T_eff = 300 K and log g >~ 4.0 secondary, corresponding to masses of 10 and 7 M_Jupiter and an age around 2 Gyr. If our deconvolution is correct, then the T_eff = 300 K cloud-free model fluxes at K and W2 are too faint by 0.5 - 1.0 magnitudes. We will address this discrepancy in our next generation of models, which will incorporate water clouds and mixing.

Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1-1.7 microns and broadband 4.5 microns. Spectroscopic observations were taken with Wide Field Camera 3 on board the Hubble Space Telescope and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability with a period of 1.4 hr at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ from wavelength to wavelength with a maximum difference of more than half of a rotational period. We compare the spectra with atmospheric models of different cloud prescriptions, from which we determine the pressure levels probed at different wavelengths. We find that the phase lag increases with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that could cause this trend of light curve phase with probed pressure level. These observations are the first to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus are an ideal test case for realistic three dimensional simulations of the atmospheric structure with clouds in brown dwarfs and extrasolar planets.

As brown dwarfs cool, a variety of species condense in their atmospheres, forming clouds. Iron and silicate clouds shape the emergent spectra of L dwarfs, but these clouds dissipate at the L/T transition. A variety of other condensates are expected to form in cooler T dwarf atmospheres. These include Cr, MnS, Na2S, ZnS, and KCl, but the opacity of these optically thinner clouds has not been included in previous atmosphere models. Here, we examine their effect on model T and Y dwarf atmospheres. The cloud structures and opacities are calculated using the Ackerman & Marley (2001) cloud model, which is coupled to an atmosphere model to produce atmospheric pressure-temperature profiles in radiative-convective equilibrium. We generate a suite of models between Teff = 400 and 1300 K, log g=4.0 and 5.5, and condensate sedimentation efficiencies from fsed=2 to 5. Model spectra are compared to two red T dwarfs, Ross 458C and UGPS 0722-05; models that include clouds are found to match observed spectra significantly better than cloudless models. The emergence of sulfide clouds in cool atmospheres, particularly Na2S, may be a more natural explanation for the "cloudy" spectra of these objects, rather than the re-emergence of silicate clouds that wane at the L-to-T transition. We find that sulfide clouds provide a mechanism to match the near- and mid-infrared colors of observed T dwarfs. Our results indicate that including the opacity of condensates in T dwarf atmospheres is necessary to accurately determine the physical characteristics of many of the observed objects.

Jupiter is expected to pulsate in a spectrum of acoustic modes and recent re-analysis of a spectroscopic time series has identified a regular pattern in the spacing of the frequencies \citep{gaulme2011}. This exciting result can provide constraints on gross Jovian properties and warrants a more in-depth theoretical study of the seismic structure of Jupiter. With current instrumentation, such as the SYMPA instrument \citep{schmider2007} used for the \citet{gaulme2011} analysis, we assume that, at minimum, a set of global frequencies extending up to angular degree $\ell=25$ could be observed. In order to identify which modes would best constrain models of Jupiter's interior and thus help motivate the next generation of observations, we explore the sensitivity of derived parameters to this mode set. Three different models of the Jovian interior are computed and the theoretical pulsation spectrum from these models for $\ell\leq 25$ is obtained. We compute sensitivity kernels and perform linear inversions to infer details of the expected discontinuities in the profiles in the Jovian interior. We find that the amplitude of the sound-speed jump of a few percent in the inner/outer envelope boundary seen in two of the applied models should be reasonably inferred with these particular modes. Near the core boundary where models predict large density discontinuities, the location of such features can be accurately measured, while their amplitudes have more uncertainty. These results suggest that this mode set would be sufficient to infer the radial location and strength of expected discontinuities in Jupiter's interior, and place strong constraints on the core size and mass. We encourage new observations to detect these Jovian oscillations.

The near-infrared colors of the planets directly imaged around the A star HR 8799 are much redder than most field brown dwarfs of the same effective temperature. Previous theoretical studies of these objects have concluded that the atmospheres of planets b, c, and d are unusually cloudy or have unusual cloud properties. Some studies have also found that the inferred radii of some or all of the planets disagree with expectations of standard giant planet evolution models. Here we compare the available data to the predictions of our own set of atmospheric and evolution models that have been extensively tested against observations of field L and T dwarfs, including the reddest L dwarfs. Unlike some previous studies we require mutually consistent choices for effective temperature, gravity, cloud properties, and planetary radius. This procedure thus yields plausible values for the masses, effective temperatures, and cloud properties of all three planets. We find that the cloud properties of the HR 8799 planets are not unusual but rather follow previously recognized trends, including a gravity dependence on the temperature of the L to T spectral transition--some reasons for which we discuss. We find the inferred mass of planet b is highly sensitive to whether or not we include the H and K band spectrum in our analysis. Solutions for planets c and d are consistent with the generally accepted constraints on the age of the primary star and orbital dynamics. We also confirm that, like in L and T dwarfs and solar system giant planets, non-equilibrium chemistry driven by atmospheric mixing is also important for these objects. Given the preponderance of data suggesting that the L to T spectral type transition is gravity dependent, we present an exploratory evolution calculation that accounts for this effect. Finally we recompute the the bolometric luminosity of all three planets.

We present multiple-epoch photometric monitoring in the $J$, $H$, and $K_s$ bands of the T1.5 dwarf 2MASS J21392676+0220226 (2M2139), revealing persistent, periodic ($P=7.721\pm$0.005 hr) variability with a peak-to-peak amplitude as high as 26% in the $J$-band. The light curve shape varies on a timescale of days, suggesting that evolving atmospheric cloud features are responsible. Using interpolations between model atmospheres with differing cloud thicknesses to represent a heterogeneous surface, we find that the multi-wavelength variations and the near-infrared spectrum of 2M2139 can be reproduced by either (1)cool, thick cloud features sitting above a thinner cloud layer, or (2)warm regions of low condensate opacity in an otherwise cloudy atmosphere, possibly indicating the presence of holes or breaks in the cloud layer. We find that temperature contrasts between thick and thin cloud patches must be greater than 175 K and as high as 425 K. We also consider whether the observed variability could arise from an interacting binary system, but this scenario is ruled out. 2M2139 joins the T2.5 dwarf SIMP0136 discovered by Artigau and coworkers as the second L/T transition brown dwarf to display large-amplitude variability on rotational timescales, suggesting that the fragmentation of dust clouds at the L/T transition may contribute to the abrupt decline in condensate opacity and $J$-band brightening observed to occur over this regime.

We present new cloudy and cloudless model atmospheres for brown dwarfs using recent ab initio calculations of the line list of ammonia (NH3) and of the collision-induced absorption of molecular hydrogen (H2). We compare the new synthetic spectra with models based on an earlier description of the H2 and NH3 opacities. We find a significant improvement in fitting the nearly complete spectral energy distribution of the T7p dwarf Gliese 570D and in near infrared color-magnitude diagrams of field brown dwarfs. We apply these new models to the identification of NH3 absorption in the H band peak of very late T dwarfs and the new Y dwarfs and discuss the observed trend in the NH3-H spectral index. The new NH3 line list also allows a detailed study of the medium resolution spectrum of the T9/T10 dwarf UGPS J072227.51-054031.2 where we identify several specific features caused by NH3.

We present i and z photometry for 25 T dwarfs and one L dwarf. Combined with published photometry, the data show that the i - z, z - Y and z - J colors of T dwarfs are very red, and continue to increase through to the late-type T dwarfs, with a hint of a saturation for the latest types with T_eff ~ 600 K. We present new 0.7-1.0 um and 2.8-4.2 um spectra for the very late-type T dwarf UGPS J072227.51-054031.2, as well as improved astrometry for this dwarf. Examination of the spectral energy distribution using the new and published data, with Saumon & Marley models, shows that the dwarf has T_eff = 505 +/- 10 K, a mass of 3-11 M_Jupiter and an age between 60 Myr and 1 Gyr. This young age is consistent with the thin disk kinematics of the dwarf. The mass range overlaps with that usually considered to be planetary, despite this being an unbound object discovered in the field near the Sun. This apparently young rapid rotator is also undergoing vigorous atmospheric mixing, as determined by the IRAC and WISE-2 4.5 um photometry and the Saumon & Marley models. The optical spectrum for this 500 K object shows clearly detected lines of the neutral alkalis Cs and Rb, which are emitted from deep atmospheric layers with temperatures of 900-1200 K.

Using two epochs of 4.5um images from the Infrared Array Camera (IRAC) on board the Spitzer Space Telescope, we recently identified a common proper motion companion to the white dwarf WD 0806-661 that is a candidate for the coldest known brown dwarf. To verify its cool nature, we have obtained images of this object at 3.6um with IRAC, at J with HAWK-I on the Very Large Telescope, and in a filter covering the red half of J with FourStar on Magellan. WD 0806-661 B is detected by IRAC but not HAWK-I or FourStar. From these data we measure colors of [3.6]-[4.5]=2.77+/-0.16 and J-[4.5]>7.0 (SNR<3). Based on these colors and its absolute magnitudes, WD 0806-661 B is the coldest companion directly imaged outside of the solar system and is a contender for the coldest known brown dwarf with the Y dwarf WISEP J1828+2650. It is unclear which of these two objects is colder given the available data. A comparison of its absolute magnitude at 4.5um to the predictions of theoretical spectra and evolutionary models suggests that WD 0806-661 B has T=300-345 K.

We present the discovery of seven ultracool brown dwarfs identified with the Wide-field Infrared Survey Explorer (WISE). Near-infrared spectroscopy reveals deep absorption bands of H_2O and CH_4 that indicate all seven of the brown dwarfs have spectral types later than UGPS J072227.51-054031.2, the latest type T dwarf currently known. The spectrum of WISEP J182831.08+265037.8 is distinct in that the heights of the J- and H-band peaks are approximately equal in units of f_lambda, so we identify it as the archetypal member of the Y spectral class. The spectra of at least two of the other brown dwarfs exhibit absorption on the blue wing of the H-band peak that we tentatively ascribe to NH_3. These spectral morphological changes provide a clear transition between the T dwarfs and the Y dwarfs. In order to produce a smooth near-infrared spectral sequence across the T/Y dwarf transition, we have reclassified UGPS J0722-0540 as the T9 spectral standard and tentatively assign WISEP J173835.52+273258.9 as the Y0 spectral standard. In total, six of the seven new brown dwarfs are classified as Y dwarfs: four are classified as Y0, one is classified as Y0 (pec?), and WISEP J1828+2650 is classified as >Y0. We have also compared the spectra to the model atmospheres of Marley and Saumon and infer that the brown dwarfs have effective temperatures ranging from 300 K to 500 K, making them the coldest spectroscopically confirmed brown dwarfs known to date.

It is becoming clear that the atmospheres of the young, self-luminous extrasolar giant planets imaged to date are dusty. Planets with dusty atmospheres may exhibit detectable amounts of linear polarization in the near-infrared, as has been observed from some field L dwarfs. The asymmetry required in the thermal radiation field to produce polarization may arise either from the rotation-induced oblateness or from surface inhomogeneities, such as partial cloudiness. While it is not possible at present to predict the extent to which atmospheric dynamics on a given planet may produce surface inhomogeneities substantial enough to produce net non-zero disk integrated polarization, the contribution of rotation-induced oblateness can be estimated. Using a self-consistent, spatially homogeneous atmospheric model and a multiple scattering polarization formalism for this class of exoplanets, we show that polarization on the order of 1% may arise due to the rotation-induced oblateness of the planets. The degree of polarization for cloudy planets should peak at the same wavelengths at which the planets are brightest in the near-infrared. The observed polarization may be even higher if surface inhomogeneities exist and play a significant role. Polarized radiation from self-luminous gas giant exoplanets, if detected, provides an additional tool to characterize these young planets and a new method to constrain their surface gravity and masses.

We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b can help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically removed at higher altitudes. Atomic species, unsaturated hydrocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and NH2 are enhanced overequilibrium predictions because of quenching and photochemistry. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced before its eventual high-altitude destruction. The enhanced abundances of HCN, CH4, and NH3 in particular are expected to affect the spectral signatures and thermal profiles HD 189733b and other, relatively cool, close-in transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficientss and discuss further observational consequences of these models.

We present near-infrared photometry and spectroscopy, and warm-Spitzer IRAC photometry of the young very cool T dwarf Ross 458C, which we have typed as T8.5p. By applying the fiducial age constraints (< 1Gyr) imposed by the properties of the active M dwarf Ross 458A, we have used these data to determine that Ross 458C has Teff = 695 \pm 60K, log g = 4.0 - 4.7 and an inferred mass of 5-20 MJup. We have compared fits of the near-infrared spectrum and IRAC photometry to the BT Settl and Saumon & Marley model grids, and have found that both sets provide best fits that are consistent with our derived properties, whilst the former provide a marginally closer match to the data for all scenarios explored here. The main difference between the model grids arises in the 4.5{\mu}m region, where the BT Settl models are able to better predict the flux through the IRAC filter, suggesting that non-equilibrium effects on the CO-CO2 ratio are important for shaping the mid-infrared spectra of very cool T dwarfs. We have also revisited the issue of dust opacity in the spectra of Ross 458C that was raised by Burgasser et al (2010). We have found that the BT Settl models which also incorporate a condensate cloud model, provide a better match to the near-infrared spectrum of this target than the Saumon & Marley model with fsed = 2, and we briefly discuss the influence of condensate clouds on T dwarf spectra.

We compute grids of radiative-convective model atmospheres for Jupiter, Saturn, Uranus, and Neptune over a range of intrinsic fluxes and surface gravities. The atmosphere grids serve as an upper boundary condition for models of the thermal evolution of the planets. Unlike previous work, we customize these grids for the specific properties of each planet, including the appropriate chemical abundances and incident fluxes as a function of solar system age. Using these grids, we compute new models of the thermal evolution of the major planets in an attempt to match their measured luminosities at their known ages. Compared to previous work, we find longer cooling times, predominantly due to higher atmospheric opacity at young ages. For all planets, we employ simple "standard" cooling models that feature adiabatic temperature gradients in the interior H/He and water layers, and an initially hot starting point for the calculation of subsequent cooling. For Jupiter we find a model cooling age 10% longer than previous work, a modest quantitative difference. This may indicate that the hydrogen equation of state used here overestimates the temperatures in the deep interior of the planet. For Saturn we find a model cooling age 20% longer than previous work. However, an additional energy source, such as that due to helium phase separation, is still clearly needed. For Neptune, unlike in work from the 1980s and 1990s, we match the measured Teff of the planet with a model that also matches the planet's current gravity field. This is predominantly due to advances in the equation of state of water. This may indicate that the planet possesses no barriers to efficient convection in its deep interior. However, for Uranus, our models exacerbate the well-known problem that Uranus is far cooler than calculations predict, which could imply strong barriers to interior convective cooling.

Mid-infrared data, including Spitzer warm-IRAC [3.6] and [4.5] photometry, is critical for understanding the cold population of brown dwarfs now being found, objects which have more in common with planets than stars. As effective temperature (T_eff) drops from 800 K to 400 K, the fraction of flux emitted beyond 3 microns increases rapidly, from about 40% to >75%. This rapid increase makes a color like H-[4.5] a very sensitive temperature indicator, and it can be combined with a gravity- and metallicity-sensitive color like H-K to constrain all three of these fundamental properties, which in turn gives us mass and age for these slowly cooling objects. Determination of mid-infrared color trends also allows better exploitation of the WISE mission by the community. We use new Spitzer Cycle 6 IRAC photometry, together with published data, to present trends of color with type for L0 to T10 dwarfs. We also use the atmospheric and evolutionary models of Saumon & Marley to investigate the masses and ages of 13 very late-type T dwarfs, which have H-[4.5] > 3.2 and T_eff ~ 500 K to 750 K.

Favored theories of giant planet formation center around two main paradigms, namely the core accretion model and the gravitational instability model. These two formation scenarios support the hypothesis that the giant planet metallicities should be higher or equal to that of the parent star. Meanwhile, spectra of the transiting hot Jupiter HD189733b suggest that carbon and oxygen abundances range from depleted to enriched with respect to the star. Here, using a model describing the formation sequence and composition of planetesimals in the protoplanetary disk, we determine the range of volatile abundances in the envelope of HD189733b that is consistent with the 20--80 Earth-masses of heavy elements estimated to be present in the planet's envelope. We then compare the inferred carbon and oxygen abundances to those retrieved from spectroscopy and we find a range of supersolar values that directly fit both spectra and internal structure models. In some cases, we find that the apparent contradiction between the subsolar elemental abundances and the {mass of heavy elements predicted in HD189733b by internal structure models} can be explained by the presence of large amounts of carbon molecules in the form of polycyclic aromatic hydrocarbons and soots in the upper layers of the envelope, as suggested by recent photochemical models. A diagnostic test that would confirm the presence of these compounds in the envelope is the detection of acetylene.

We report the discovery of the first new ultra-cool brown dwarf found with the Wide-field Infrared Survey Explorer (WISE). The object's preliminary designation is WISEPC J045853.90+643451.9. Follow-up spectroscopy with the LUCIFER instrument on the Large Binocular Telescope indicates that it is a very late-type T dwarf with a spectral type approximately equal to T9. Fits to an IRTF/SpeX 0.8-2.5 micron spectrum to the model atmospheres of Marley and Saumon indicate an effective temperature of approximately 600 K as well as the presence of vertical mixing in its atmosphere. The new brown dwarf is easily detected by WISE, with a signal-to-noise ratio of ~36 at 4.6 microns. Current estimates place it at a distance of 6 to 10 pc. This object represents the first in what will likely be hundreds of nearby brown dwarfs found by WISE that will be suitable for follow up observations, including those with the James Webb Space Telescope. One of the two primary scientific goals of the WISE mission is to find the coolest, closest stars to our Sun; the discovery of this new brown dwarf proves that WISE is capable of fulfilling this objective.

GJ436b is a unique member of the transiting extrasolar planet population being one of the smallest and least irradiated and possessing an eccentric orbit. Because of its size, mass and density, GJ436b could plausibly have an atmospheric metallicity similar to Neptune (20-60 times solar abundances), which makes it an ideal target to study the effects of atmospheric metallicity on dynamics and radiative transfer in an extrasolar planetary atmosphere. We present three-dimensional atmospheric circulation models that include realistic non-gray radiative transfer for 1, 3, 10, 30, and 50 times solar atmospheric metallicity cases of GJ436b. Low metallicity models (1 and 3 times solar) show little day/night temperature variation and strong high-latitude jets. In contrast, higher metallicity models (30 and 50 times solar) exhibit day/night temperature variations and a strong equatorial jet. Spectra and light curves produced from these simulations show strong orbital phase dependencies in the 50 times solar case and negligible variations with orbital phase in the 1 times solar case. Comparisons between the predicted planet/star flux ratio from these models and current secondary eclipse measurements support a high metallicity atmosphere (30-50 times solar abundances) with disequilibrium carbon chemistry at play for GJ436b. Regardless of the actual atmospheric composition of GJ436b, our models serve to illuminate how metallicity influences the atmospheric circulation for a broad range of warm extrasolar planets.

One of the mechanisms suggested for the L to T dwarf spectral type transition is the appearance of relatively cloud-free regions across the disk of brown dwarfs as they cool. The existence of partly cloudy regions has been supported by evidence for variability in dwarfs in the late L to early T spectral range, but no self-consistent atmosphere models of such partly cloudy objects have yet been constructed. Here we present a new approach for consistently modeling partly cloudy brown dwarfs and giant planets. We find that even a small fraction of cloud holes dramatically alter the atmospheric thermal profile, spectra, and photometric colors of a given object. With decreasing cloudiness objects briskly become bluer in J - K and brighten in J band, as is observed at the L/T transition. Model spectra of partly cloudy objects are similar to our models with globally homogenous, but thinner, clouds. Hence spectra alone may not be sufficient to distinguish partial cloudiness although variability and polarization measurements are potential observational signatures. Finally we note that partial cloud cover may be an alternative explanation for the blue L dwarfs.

Condensate clouds are a salient feature of L dwarf atmospheres, but have been assumed to play little role in shaping the spectra of the coldest T-type brown dwarfs. Here we report evidence of condensate opacity in the near-infrared spectrum of the brown dwarf candidate Ross 458C, obtained with the Folded-Port Infrared Echellette (FIRE) spectrograph at the Magellan Telescopes. These data verify the low-temperature nature of this source, indicating a T8 spectral classification, log Lbol/Lsun = -5.62+/-0.03, Teff = 650+/-25 K, and a mass at or below the deuterium burning limit. The data also reveal enhanced emission at K-band associated with youth (low surface gravity) and supersolar metallicity, reflecting the properties of the Ross 458 system (age = 150-800 Myr, [Fe/H] = +0.2 to +0.3). We present fits of FIRE data for Ross 458C, the T9 dwarf ULAS J133553.45+113005.2, and the blue T7.5 dwarf SDSS J141624.08+134826.7B, to cloudless and cloudy spectral models from Saumon & Marley. For Ross 458C we confirm a low surface gravity and supersolar metallicity, while the temperature differs depending on the presence (635 [+25,-35] K) or absence (760 [+70,-45] K) of cloud extinction. ULAS J1335+1130 and SDSS J1416+1348B have similar temperatures (595 [+25,-45] K), but distinct surface gravities (log g = 4.0-4.5 cgs versus 5.0-5.5 cgs) and metallicities ([M/H] ~ +0.2 versus -0.2). In all three cases, cloudy models provide better fits to the spectral data, significantly so for Ross 458C. These results indicate that clouds are an important opacity source in the spectra of young cold T dwarfs, and should be considered when characterizing the spectra of planetary-mass objects in young clusters and directly-imaged exoplanets. The characteristics of Ross 458C suggest it could itself be regarded as a planet, albeit one whose cosmogony does not conform with current planet formation theories.

First generation optical coronagraphic telescopes will obtain images of cool gas and ice giant exoplanets around nearby stars. The albedo spectra of exoplanets at planet-star separations larger than about 1 AU are dominated by reflected light to beyond 1 {\mu}m and are punctuated by molecular absorption features. We consider how exoplanet albedo spectra and colors vary as a function of planet-star separation, metallicity, mass, and observed phase for Jupiter and Neptune analogs from 0.35 to 1 {\mu}m. We model Jupiter analogs with 1x and 3x the solar abundance of heavy elements, and Neptune analogs with 10x and 30x. Our model planets orbit a solar analog parent star at separations of 0.8 AU, 2 AU, 5 AU, and 10 AU. We use a radiative-convective model to compute temperature-pressure profiles. The giant exoplanets are cloud-free at 0.8 AU, have H2O clouds at 2 AU, and have both NH3 and H2O clouds at 5 AU and 10 AU. For each model planet we compute moderate resolution spectra as a function of phase. The presence and structure of clouds strongly influence the spectra. Since the planet images will be unresolved, their phase may not be obvious, and multiple observations will be needed to discriminate between the effects of planet-star separation, metallicity, and phase. We consider the range of these combined effects on spectra and colors. For example, we find that the spectral influence of clouds depends more on planet-star separation and hence temperature than metallicity, and it is easier to discriminate between cloudy 1x and 3x Jupiters than between 10x and 30x Neptunes. In addition to alkalis and methane, our Jupiter models show H2O absorption features near 0.94 {\mu}m. We also predict that giant exoplanets receiving greater insolation than Jupiter will exhibit higher equator to pole temperature gradients than are found on Jupiter and thus may have differing atmospheric dynamics.

We present an analysis of the recently discovered blue L dwarf SDSS J141624.08+134826.7. We extend the spectral coverage of its published spectrum to ~4 microns by obtaining a low-resolution L band spectrum with SpeX on the NASA IRTF. The spectrum exhibits a tentative weak CH4 absorption feature at 3.3 microns but is otherwise featureless. We derive the atmospheric parameters of SDSS J141624.08+134826.7 by comparing its 0.7-4.0 micron spectrum to the atmospheric models of Marley and Saumon which include the effects of both condensate cloud formation and non-equilibrium chemistry due to vertical mixing and find the best fitting model has Teff=1700 K, log g=5.5 [cm s-2], fsed=4, and Kzz=10^4 cm2 s-1. The derived effective temperature is significantly cooler than previously estimated but we confirm the suggestion by Bowler et al. that the peculiar spectrum of SDSS J141624.08+134826.7 is primarily a result of thin condensate clouds. In addition, we find strong evidence of vertical mixing in the atmosphere of SDSS J141624.08+134826.7 based on the absence of the deep 3.3 micron CH4 absorption band predicted by models computed in chemical equilibrium. This result suggests that observations of blue L dwarfs are an appealing way to quantitatively estimate the vigor of mixing in the atmospheres of L dwarfs because of the dramatic impact such mixing has on the strength of the 3.3 micron CH4 band in the emergent spectra of L dwarfs with thin condensate clouds.

Light scattering by atmospheric dust particles is responsible for the polarization observed in some L dwarfs. Whether this polarization arises from an inhomogeneous distribution of dust across the disk or an oblate shape induced by rotation remains unclear. Here we argue that the latter case is plausible and, for many L dwarfs, the more likely one. Furthermore evolutionary models of mature field L dwarfs predict surface gravities ranging from about 200 to 2500 m/s^2 (corresponding to masses of about 15 to 70 times of Jupiter mass). Yet comparison of observed spectra to available synthetic spectra often does not permit more precise determination of the surface gravity of individual field L dwarfs, leading to important uncertainties in their properties. Since rotationally-induced non-sphericity, which gives rise to non-zero disk-integrated polarization, is more pronounced at lower gravities, polarization is a promising low gravity indicator. Here we combine a rigorous multiple scattering analysis with a self-consistent cloudy atmospheric model and observationally inferred rotational velocities and find that the observed optical polarization can be explained if the surface gravity of the polarized objects is about 300 m/s^2 or less, potentially providing a new method for constraining L dwarf masses.

We use a 1D thermochemical and photochemical kinetics model to predict the disequilibrium stratospheric chemistries of warm and hot Jupiters (800 < T < 1200 K). Thermal chemistry and vertical mixing are generally more important than photochemistry. At 1200 K, methane is oxidized to CO and CO2 by OH radicals from thermal decomposition of water. At T < 1000 K, methane is reactive but stable enough to reach the stratosphere, while water is stable enough that OH levels are suppressed by reaction with H2. These trends raise the effective C/O ratio in the reacting gases above unity. Reduced products such as ethylene, acetylene, and hydrogen cyanide become abundant; further polymerization should lead to formation of PAHs (Poly-Aromatic Hydrocarbons) and soots. Parallel shifts are seen in the sulfur chemistry, with CS and CS2 displacing S2 and HS as the interesting disequilibrium products. Although lower temperature is a leading factor favoring hydrocarbons, higher rates of vertical mixing, lower metallicities, and lower incident UV radiation also favor organic synthesis. Acetylene (the first step toward PAHs) formation is especially favored by high eddy diffusion coefficients Kzz > 10^10 cm2/s. In most cases planetary compositions inferred from transit observations will differ markedly from those inferred from reflected or emitted light from the same planet. The peculiar properties of HD 189733b compared to other hot Jupiters - a broadband blue haze, little sign of Na or K, and hints of low metallicity - can be explained by an organic haze if the planet is cool enough. Whether this interpretation applies to HD 189733b itself, many organic-rich warm Jupiters are sure to be discovered in the near future.

We present 7.5-14.2um low-resolution spectroscopy, obtained with the Spitzer Infrared Spectrograph, of the T8.5 dwarf Wolf 940 B, which is a companion to an M4 dwarf with a projected separation of 400 AU. We combine these data with previously published near-infrared spectroscopy and mid-infrared photometry, to produce the spectral energy distribution for the very low-temperature T dwarf. We use atmospheric models to derive the bolometric correction and obtain a luminosity of log L/Lsun = -6.01 +/- 0.05. Evolutionary models are used with the luminosity to constrain the values of effective temperature (T_eff) and surface gravity, and hence mass and age for the T dwarf. We further restrict the allowed range of T_eff and gravity using age constraints implied by the M dwarf primary, and refine the physical properties of the T dwarf by comparison of the observed and modelled spectroscopy and photometry. This comparison indicates that Wolf 940 B has a metallicity within 0.2 dex of solar, as more extreme values give poor fits to the data - lower metallicity produces a poor fit at lambda > 2um while higher metallicity produces a poor fit at lambda < 2um. This is consistent with the independently derived value of [m/H] = +0.24 +/- 0.09 for the primary star, using the Johnson & Apps (2008) M_K:V-K relationship. We find that the T dwarf atmosphere is undergoing vigorous mixing, with an eddy diffusion coefficient K_zz of 10^4 to 10^6 cm^2 s^-1. We derive an effective temperature of 585 K to 625 K, and surface gravity log g = 4.83 to 5.22 (cm s^-2), for an age range of 3 Gyr to 10 Gyr, as implied by the kinematic and H alpha properties of the M dwarf primary. The lower gravity corresponds to the lower temperature and younger age for the system, and the higher value to the higher temperature and older age. The mass of the T dwarf is 24 M_Jupiter to 45 M_Jupiter for the younger to older age limit.

This is a scientific strategy for the detection and characterization of extrasolar planets; that is, planets orbiting other stars. As such, it maps out over a 15-year horizon the techniques and capabilities required to detect and measure the properties of planets as small as Earth around stars as large as our own Sun. It shows how the technology pieces and their development fit together to achieve the primary goal of the strategy: if planets like Earth exist around stars within some tens of light years of our own Solar System, those planets will be found and their basic properties characterized. Essential to this strategy is not only the search for and examination of individual planets, but also a knowledge of the arrangement, or architecture, of planetary systems around as large a number of stars as possible; this is the second goal of the strategy. The final goal of the strategy is the study of disks around stars, important both to understand the implications of the variety of exoplanet systems for planet formation, and to determine how many nearby stars have environments around them clean enough of debris that planets may be sought and, if found, characterized.

We present Spitzer IRAC photometry of twelve very late-type T dwarfs; nine have [3.6], [4.5], [5.8] and [8.0] photometry and three have [3.6] and [4.5] photometry only. We investigate trends with type and color for the planning and interpretation of surveys for the coldest T or Y dwarfs. Brown dwarfs with effective temperature (T_eff) below 700K emit more than half their flux at wavelengths beyond 3um, and the ratio of the mid-infrared flux to the near-infrared flux becomes very sensitive to T_eff. The color H-[4.5] is a good indicator of T_eff with a relatively weak dependence on metallicity and gravity. Conversely, the colors H-K and [4.5]-[5.8] are sensitive to metallicity and gravity. Thus near- and mid-infrared photometry provide useful indicators of the fundamental properties of brown dwarfs, and if temperature and gravity are known, then mass and age can be determined from evolutionary models. There are twelve dwarfs known with H-[4.5]>3.0, and 500< T_eff K < 800, which we examine in detail. The ages of these dwarfs range from very young (0.1-1.0 Gyr) to relatively old (3-12 Gyr). The mass range is possibly as low as 5 Jupiter masses to up to 70 Jupiter masses. The metallicities also span a large range, from [m/H] -0.3 to [m/H] +0.3. The small number of T8-T9 dwarfs found in the UKIRT Infrared Deep Sky Survey to date appear to be mostly young low-mass dwarfs. Extensions to the mid-infrared space missions warm-Spitzer and WISE are desirable in order to obtain the vital mid-infrared data for cold brown dwarfs, and to discover more of these rare objects. (Abridged.)

We compute models of the transmission spectra of planets HD 209458b, HD 189733b, and generic hot Jupiters. We examine the effects of temperature, surface gravity, and metallicity for the generic planets as a guide to understanding transmission spectra in general. We find that carbon dioxide absorption at 4.4 and 15 microns is prominent at high metallicity, and is a clear metallicity indicator. For HD 209458b and HD 189733b, we compute spectra for both one-dimensional and three-dimensional model atmospheres and examine the differences between them. The differences are usually small, but can be large if atmospheric temperatures are near important chemical abundance boundaries. The calculations for the 3D atmospheres, and their comparison with data, serve as constraints on these dynamical models that complement the secondary eclipse and light curve data sets. For HD 209458b, even if TiO and VO gases are abundant on the day side, their abundances can be considerably reduced on the cooler planetary limb. However, given the predicted limb temperatures and TiO abundances, the model's optical opacity is too high. For HD 189733b we find a good match with some infrared data sets and constrain the altitude of a postulated haze layer. For this planet, substantial differences can exist between the transmission spectra of the leading and trailing hemispheres, which is an excellent probe of carbon chemistry. In thermochemical equilibrium, the cooler leading hemisphere is methane-dominated, and the hotter trailing hemisphere is CO-dominated, but these differences may be eliminated by non-equilibrium chemistry due to vertical mixing. It may be possible to constrain the carbon chemistry of this planet, and its spatial variation, with JWST.

The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.

The Terrestrial Planet Finder Coronagraph (TPF-C) mission presented here is an existence proof for a flagship-class internal coronagraph space mission capable of detecting and characterizing Earth-like planets and planetary systems at visible wavelengths around nearby stars, using an existing launch vehicle. TPF-C will use spectroscopy to measure key properties of exoplanets including the presence of atmospheric water or oxygen, powerful signatures in the search for habitable worlds.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

While there have been multiple observational programs aimed at detecting linear polarization of optical radiation emitted by ultracool dwarfs, there has been comparatively less rigorous theoretical analysis of the problem. The general expectation has been that the atmospheres of those substellar-mass objects with condensate clouds would give rise to linear polarization due to scattering. Because of rotation-induced non-sphericity, there is expected to be incomplete cancellation of disk-integrated net polarization and thus a finite polarization. For cloudless objects, however, only molecular Rayleigh scattering will contribute to any net polarization and this limit has not been well studied. Hence in this paper we present a detailed multiple scattering analysis of the polarization expected from those T-dwarfs whose spectra show absence of condensates. For this, we develop and solve the full radiative transfer equations for linearly polarized radiation. Only atomic and molecular Rayleigh scattering are considered to be the source of polarization. We compute the local polarization at different angular directions in a plane-parallel atmospheres calculated for the range of effective temperatures of T dwarfs and then average over the whole surface of the object. The effects of gravity and limb darkening as well as rotation induced non-sphericity are included. It is found that the amount of polarization decreases with the increase in effective temperature. It is also found that significant polarization at any local point in the atmosphere arises only in the optical (B-band). However, the disk integrated polarization--even in the B-band--is negligible. Hence we conclude that, unlike the case for cloudy L dwarfs, polarization of cloudless T-dwarfs by atomic and molecular scattering may not be detectable.

We develop a new 1D photochemical kinetics code to address stratospheric chemistry and stratospheric heating in hot Jupiters. Here we address optically active S-containing species and CO2 at 1200 < T < 2000 K. HS (mercapto) and S2 are highly reactive species that are generated photochemically and thermochemically from H2S with peak abundances between 1-10 mbar. S2 absorbs UV between 240 and 340 nm and is optically thick for metallicities [SH] > 0 at T > 1200 K. HS is probably more important than S2, as it is generally more abundant than S2 under hot Jupiter conditions and it absorbs at somewhat redder wavelengths. We use molecular theory to compute an HS absorption spectrum from sparse available data and find that HS should absorb strongly between 300 and 460 nm, with absorption at the longer wavelengths being temperature sensitive. When the two absorbers are combined, radiative heating (per kg of gas) peaks at 100 microbars, with a total stratospheric heating of about 8 x 10^4 W/m^2 for a jovian planet orbiting a solar-twin at 0.032 AU. Total heating is insensitive to metallicity. The CO2 mixing ratio is a well-behaved quadratic function of metallicity, ranging from 1.6 x 10^-8 to 1.6 x 10^-4 for -0.3 < [M/H] < 1.7. CO2 is insensitive to insolation, vertical mixing, temperature (1200 < T <2000 K), and gravity. The photochemical calculations confirm that CO2 should prove a useful probe of planetary metallicity.

We present global, three-dimensional numerical simulations of HD 189733b and HD 209458b that couple the atmospheric dynamics to a realistic representation of non-gray cloud-free radiative transfer. The model, which we call the Substellar and Planetary Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, and collaborators for the radiation. Like earlier work with simplified forcing, our simulations develop a broad eastward equatorial jet, mean westward flow at higher latitudes, and substantial flow over the poles at low pressure. For HD 189733b, our simulations without TiO and VO opacity can explain the broad features of the observed 8 and 24-micron light curves, including the modest day-night flux variation and the fact that the planet/star flux ratio peaks before the secondary eclipse. Our simulations also provide reasonable matches to the Spitzer secondary-eclipse depths at 4.5, 5.8, 8, 16, and 24 microns and the groundbased upper limit at 2.2 microns. However, we substantially underpredict the 3.6-micron secondary-eclipse depth, suggesting that our simulations are too cold in the 0.1-1 bar region. Predicted temporal variability in secondary-eclipse depths is ~1% at Spitzer bandpasses, consistent with recent observational upper limits at 8 microns. We also show that nonsynchronous rotation can significantly alter the jet structure. For HD 209458b, we include TiO and VO opacity; these simulations develop a hot (>2000 K) dayside stratosphere. Despite this stratosphere, we do not reproduce current Spitzer photometry of this planet. Light curves in Spitzer bandpasses show modest phase variation and satisfy the observational upper limit on day-night phase variation at 8 microns. (abridged)

We present new 5.2-14.5 micron low-resolution spectra of 14 mid-L to mid-T dwarfs, as well as 3.0-4.1 micron spectra for five of these dwarfs. These data are supplemented by existing red and near-infrared spectra (0.6-2.5 micron), as well as red through mid-infrared spectroscopy of seven other L and T dwarfs presented by Cushing et al. (2008). We compare these spectra to those generated from the model atmospheres of Saumon & Marley (2008). The models reproduce the observed spectra well, except in the case of one very red L3.5 dwarf, 2MASS J22244381-0158521. The broad wavelength coverage allows us to constrain almost independently the four parameters used to describe these photospheres in our models: effective temperature (Teff), surface gravity, grain sedimentation efficiency (fsed) and vertical gas transport efficiency (Kzz). The sample of L3.5 to T5.5 dwarfs spans the range 1800 < Teff < 1000 K, with an L-T transition (spectral types L7-T4) that lies between 1400 K and 1100K for dwarfs with typical near-infrared colors; bluer and redder dwarfs can be 100 K warmer or cooler, respectively, when using infrared spectral types. In this model analysis, fsed increases rapidly between types T0 and T4, indicating that increased sedimentation can explain the rapid disappearance of clouds at this stage of brown dwarf evolution. There is a suggestion that the transition to dust-free atmospheres happens at lower temperatures for lower gravity dwarfs.

Significant advances in the discovery and characterization of the planetary systems of nearby stars can be accomplished with a moderate aperture high performance coronagraphic space mission that could be started in the next decade. Its observations would make significant progress in studying terrestrial planets in their habitable zones to giant planets and circumstellar debris disks, also informing the design of a more capable future mission. It is quite exciting that such fundamental exoplanet science can be done with relatively modest capabilities.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee highlights cross-disciplinary science opportunities over the next decade with cold brown dwarfs, sources defined here as having photospheric temperatures less than ~1000 K.

Despite the revolution in our knowledge resulting from the detection of planets around mature stars, we know almost nothing about planets orbiting young stars because rapid rotation and active photospheres preclude detection by radial velocities or transits and because direct imaging has barely penetrated the requisite range of high contrast and angular resolution. Of the techniques presently under consideration for the coming decade, only space-based astrometry offers the prospect of discovering gas giants (100 to >> 300 Mearth), lower mass systems such as icy giants (10 to 100 Mearth), and even a few rocky, super-Earths 300 Mearth) orbiting stars ranging in age from 1 to 100 Myr. Astrometry will complement high contrast imaging which should be able to detect gas giants (1~10 MJup) in orbits from a few to a few hundred AU. An astrometric survey in combination with imaging data for a subsample of objects will allow a detailed physical understanding of the formation and evolution of young gas giant planets impossible to achieve by any one technique.

Newly-formed planetary systems with ages of <10 Myr offer many unique insights into the formation, evolution, and fundamental properties of extrasolar planets. These planets have fallen beyond the limits of past surveys, but as we enter the next decade, we stand on the threshold of several crucial advances in instrumentation and observing techniques that will finally unveil this critical population. In this white paper, we consider several classes of planets (inner gas giants, outer gas giants, and ultrawide planetary-mass companions) and summarize the motivation for their study, the the observational tests that will distinguish between competing theoretical models, and the infrastructure investments and policy choices that will best enable future discovery. We propose that there are two fundamental questions that must be addressed: 1) Do planets form via core accretion, gravitational instability, or a combination of the two methods? 2) What do the atmospheres an interiors of young planets look like, and does the mass-luminosity relation of young planets more closely resemble the "hot start" or "cold start" models? To address these questions, we recommend investment in high-resolution NIR spectrographs (existing and new), support for innovative new techniques and pathfinder surveys for directly-imaged young exoplanets, and continued investment in visible-light adaptive optics to allow full characterization of wide "planetary-mass" companions for calibrating planet evolutionary models. In summary, testing newly proposed planet formation and evolutionary predictions will require the identification of a large population of young (<10 Myr) planets whose orbital, atmospheric, and structural properties can be studied.

We present Spitzer 7.6-14.5um spectra of ULAS J003402.77-005206.7 and ULAS J133553.45+113005.2, two T9 dwarfs with the latest spectral types currently known. We fit synthetic spectra and photometry to the near- through mid-infrared energy distributions of these dwarfs and that of the T8 dwarf 2MASS J09393548-2448279. We also analyse near-infrared data for another T9, CFBD J005910.82-011401.3. We find that the ratio of the mid- to near-infrared fluxes is very sensitive to effective temperature at these low temperatures, and that the 2.2 and 4.5um fluxes are sensitive to metallicity and gravity; there is a degeneracy between these parameters. The 4.5 and 10um fluxes are also sensitive to vertical transport of gas through the atmosphere, which we find to be significant for these dwarfs. The full near- through mid-infrared spectral energy distribution allows us to constrain the effective temperature (K)/gravity (m/s2)/metallicity ([m/H] dex) of ULAS J0034-00 and ULAS J1335+11 to 550-600/ 100-300/ 0.0-0.3 and 500-550/ 100-300/ 0.0-0.3, respectively. These fits imply low masses and young ages for the dwarfs of 5-20 M(Jup) and 0.1-2 Gyr. The fits to 2MASS J0939-24 are in good agreement with the measured distance, the observational data, and the earlier T8 near-infrared spectral type if it is a slightly metal-poor 4-10 Gyr old system consisting of a 500 and 700K, ~25 and ~40 M(Jup), pair, although it is also possible that it is an identical pair of 600K, 30 M(Jup), dwarfs. As no mid-infrared data are available for CFBD J0059-01 its properties are less well constrained; nevertheless it appears to be a 550-600K dwarf with g= 300-2000 m/s2 and [m/H]= 0-0.3 dex. These properties correspond to mass and age ranges of 10-50 M(Jup) and 0.5-10 Gyr for this dwarf.

M band spectra of two late-type T dwarfs, 2MASS J09373487+2931409, and Gliese 570D, confirm evidence from photometry that photospheric CO is present at abundance levels far in excess of those predicted from chemical equilibrium. These new and unambiguous detections of CO, together with an earlier spectroscopic detection of CO in Gliese 229B and existing M band photometry of a large selection of T dwarfs, suggest that vertical mixing in the photosphere drives the CO abundance out of chemical equilibrium and is a common, and likely universal feature of mid-to-late type T dwarfs. The M band spectra allow determinations of the time scale of vertical mixing in the atmosphere of each object, the first such measurements of this important parameter in late T dwarfs. A detailed analysis of the spectral energy distribution of 2MASS J09373487+2931409 results in the following values for metallicity, temperature, surface gravity, and luminosity: [M/H]~-0.3, T_eff=925-975K, log g=5.20-5.47, log L/L_sun=-5.308 +/- 0.027. The age is 3-10 Gyr and the mass is in the range 45-69 M_Jup.

The Super-Earth Explorer is an Off-Axis Space Telescope (SEE-COAST) designed for high contrast imaging. Its scientific objective is to make the physico-chemical characterization of exoplanets possibly down to 2 Earth radii >. For that purpose it will analyze the spectral and polarimetric properties of the parent starlight reflected by the planets, in the wavelength range 400-1250 nm

The characteristics of irradiated solar system planetary atmospheres have been studied for decades, consequently modern planetary science benefits from an exhaustive body of ground- and space-based data. The study of extrasolar planetary atmospheres, by contrast, is still in its infancy and currently rests on a few score of datapoints, mostly of the transiting planets. This short survey aims not to review this dynamic field but rather stresses the importance of a few theoretical concepts and processes for our understanding of exoplanet atmospheres. Topics covered include atmospheric structure and dynamics, cloud processes and photochemistry of planetary atmospheres. Influences on the albedos, spectra, and colors of extrasolar planets are reviewed and caution is urged in the interpretation of exoplanet colors.

The atmospheres of substellar objects contain clouds of oxides, iron, silicates, and other refractory condensates. Water clouds are expected in the coolest objects. The opacity of these `dust' clouds strongly affects both the atmospheric temperature-pressure profile and the emergent flux. Thus any attempt to model the spectra of these atmospheres must incorporate a cloud model. However the diversity of cloud models in atmospheric simulations is large and it is not always clear how the underlying physics of the various models compare. Likewise the observational consequences of different modeling approaches can be masked by other model differences, making objective comparisons challenging. In order to clarify the current state of the modeling approaches, this paper compares five different cloud models in two sets of tests. Test case 1 tests the dust cloud models for a prescribed L, L--T, and T-dwarf atmospheric (temperature T, pressure p, convective velocity vconv)-structures. Test case 2 compares complete model atmosphere results for given (effective temperature Teff, surface gravity log g). All models agree on the global cloud structure but differ in opacity-relevant details like grain size, amount of dust, dust and gas-phase composition. Comparisons of synthetic photometric fluxes translate into an modelling uncertainty in apparent magnitudes for our L-dwarf (T-dwarf) test case of 0.25 < \Delta m < 0.875 (0.1 < \Delta m M 1.375) taking into account the 2MASS, the UKIRT WFCAM, the Spitzer IRAC, and VLT VISIR filters with UKIRT WFCAM being the most challenging for the models. (abr.)

We present new evolution sequences for very low mass stars, brown dwarfs and giant planets and use them to explore a variety of influences on the evolution of these objects. We compare our results with previous work and discuss the causes of the differences and argue for the importance of the surface boundary condition provided by atmosphere models including clouds. The L- to T-type ultracool dwarf transition can be accommodated within the Ackerman & Marley (2001) cloud model by varying the cloud sedimentation parameter. We develop a simple model for the evolution across the L/T transition. By combining the evolution calculation and our atmosphere models, we generate colors and magnitudes of synthetic populations of ultracool dwarfs in the field and in galactic clusters. We focus on near infrared color- magnitude diagrams (CMDs) and on the nature of the ``second parameter'' that is responsible for the scatter of colors along the Teff sequence. Variations in metallicity and cloud parameters, unresolved binaries and possibly a relatively young population all play a role in defining the spread of brown dwarfs along the cooling sequence. We find that the transition from cloudy L dwarfs to cloudless T dwarfs slows down the evolution and causes a pile up of substellar objects in the transition region, in contradiction with previous studies. We apply the same model to the Pleiades brown dwarf sequence. Taken at face value, the Pleiades data suggest that the L/T transition occurs at lower Teff for lower gravity objects. The simulated populations of brown dwarfs also reveal that the phase of deuterium burning produces a distinctive feature in CMDs that should be detectable in ~50-100 Myr old clusters.

We examine the spectra and infrared colors of the cool methane-dominated atmospheres at Teff < 1400 K expected for young gas giant planets. We couple these spectral calculations to an updated version of the Marley et al. (2007) giant planet thermal evolution models that include formation by core accretion-gas capture. These relatively cool "young Jupiters" can be 1-6 magnitudes fainter than predicted by standard cooling tracks that include a traditional initial condition, which may provide a diagnostic of formation. If correct, this would make true Jupiter-like planets much more difficult to detect at young ages than previously thought. Since Jupiter and Saturn are of distinctly super-solar composition, we examine emitted spectra for model planets at both solar metallicity and a metallicity of 5 times solar. These metal-enhanced young Jupiters have lower pressure photospheres than field brown dwarfs of the same effective temperatures arising from both lower surface gravities and enhanced atmospheric opacity. We highlight several diagnostics for enhanced metallicity. A stronger CO absorption band at 4.5 $\mu$m for the warmest objects is predicted. At all temperatures, enhanced flux in $K$ band is expected due to reduced collisional induced absorption by H$_2$. This leads to correspondingly redder near infrared colors, which are redder than solar metallicity models with the same surface gravity by up to 0.7 in $J-K$ and 1.5 in $H-K$. Molecular absorption band depths increase as well, most significantly for the coolest objects. We also qualitatively assess the changes to emitted spectra due to nonequilibrium chemistry.

We report the discovery of a single-lined spectroscopic binary with an Ultra Cool Dwarf (UCD) primary with a spectral type between M8 and L0.5. This system was discovered during the course of an ongoing survey to monitor L dwarfs for radial velocity variations and is the first known small separation (a<1 AU) spectroscopic binary among dwarfs at the M/L transition. Based on radial-velocity measurements with a typical precision of 300 m/s we estimate the orbital parameters of this system to be P=246.73+/-0.49 d, a1 sin(i)=0.159+/-0.003 AU, M2 sin(i)=0.2062 (M1+M2)^(2/3)+/-0.0034 M_{\sun}. Assuming a primary mass of M1=0.08M_{\sun} (based on spectral type), we estimate the secondary minimum mass to be M2 sin(i)=0.054 M_{\sun}. With future photometric, spectroscopic, and interferometric observations it may be possible to determine the dynamical masses of both components directly, making this system one of the best characterized UCD binaries known.

Luhman and collaborators recently discovered an early-T dwarf companion to the G0 dwarf star HN Peg, using Spitzer Infrared Array Camera (IRAC) images. Companionship was established on the basis of the common proper motion inferred from 1998 Two Micron All Sky Survey images and the 2004 IRAC images. In this paper we present new near-infrared imaging data which confirms the common proper motion of the system. We also present new 3 - 4 um spectroscopy of HN Peg B, which provides tighter constraints on both the bolometric luminosity determination and the comparison to synthetic spectra. New adaptive optics imaging data are also presented, which shows the T dwarf to be unresolved, providing limits on the multiplicity of the object. We use the age, distance and luminosity of the solar-metallicity T dwarf to determine its effective temperature and gravity, and compare synthetic spectra with these values, and a range of grain properties and vertical mixing, to the observed 0.8 - 4.0 um spectra and mid-infrared photometry. We find that models with temperature and gravity appropriate for the older end of the age range of the system (0.5 Gyr) can do a reasonable job of fitting the data, but only if the photospheric condensate cloud deck is thin, and if there is significant vertical mixing in the atmosphere. Dwarfs such as HN Peg B, with well-determined metallicity, radius, gravity and temperature will allow development of dynamical atmosphere models, leading to the solution of the puzzle of the L to T dwarf transition.

Ultracool dwarfs exhibit a remarkably varied set of characteristics which hint at the complex physical processes acting in their atmospheres and interiors. Spectra of these objects not only depend upon their mass and effective temperature, but also their atmospheric chemistry, weather, and dynamics. As a consequence divining their mass, metallicity and age solely from their spectra has been a challenge. JWST, by illuminating spectral blind spots and observing objects with constrained masses and ages should finally unearth a sufficient number of ultracool dwarf Rosetta Stones to allow us to decipher the processes underlying the complex brown dwarf cooling sequence. In addition the spectra of objects invisible from the ground, including very low mass objects in clusters and nearby cold dwarfs from the disk population, will be seen for the first time. In combination with other ground- and space-based assets and programs, JWST will usher in a new golden era of brown dwarf science and discovery.

We present global, three-dimensional numerical simulations of the atmospheric circulation on HD 209458b and HD 189733b and calculate the infrared spectra and light curves predicted by these simulations, which we compare with available observations. Radiative heating/cooling is parameterized with a simplified Newtonian relaxation scheme. Our simulations develop day-night temperature contrasts that vary strongly with pressure. At low pressure (<10 mbar), air flows from the substellar point toward the antistellar point, both along the equator and over the poles. At deeper levels, the flow develops an eastward equatorial jet with speeds of 3-4 km/sec, with weaker westward flows at high latitudes. This basic flow pattern is robust to variations in model resolution, gravity, radiative time constant, and initial temperature structure. Nightside spectra show deep absorption bands of H2O, CO, and/or CH4, whereas on the dayside these absorption bands flatten out or even flip into emission. This results from the strong effect of dynamics on the vertical temperature-pressure structure; the temperature decreases strongly with altitude on the nightside but becomes almost isothermal on the dayside. In Spitzer bandpasses, our predicted planet-to-star flux ratios vary by a factor of 2-10 with orbital phase, depending on the wavelength and chemistry. For HD 189733b, where a detailed 8-micron light curve has been obtained, we correctly produce the observed phase offset of the flux maximum, but we do not explain the flux minimum and we overpredict the total flux variation. This discrepancy likely results from the simplifications inherent in the Newtonian relaxation scheme and provides motivation for incorporating realistic radiative transfer in future studies.

Context: L-type ultra-cool dwarfs and brown dwarfs have cloudy atmospheres that could host weather-like phenomena. The detection of photometric or spectral variability would provide insight into unresolved atmospheric heterogeneities, such as holes in a global cloud deck. Aims: It has been proposed that growth of heterogeneities in the global cloud deck may account for the L- to T-type transition as brown dwarf photospheres evolve from cloudy to clear conditions. Such a mechanism is compatible with variability. We searched for variability in the spectra of five L6 to T6 brown dwarfs in order to test this hypothesis. Methods: We obtained spectroscopic time series using VLT/ISAAC, over 0.99-1.13um, and IRTF/SpeX for two of our targets, in J, H and K bands. We search for statistically variable lines and correlation between those. Results: High spectral-frequency variations are seen in some objects, but these detections are marginal and need to be confirmed. We find no evidence for large amplitude variations in spectral morphology and we place firm upper limits of 2 to 3% on broad-band variability, on the time scale of a few hours. The T2 transition brown dwarf SDSS J1254-0122 shows numerous variable features, but a secure variability diagnosis would require further observations. Conclusions: Assuming that any variability arises from the rotation of patterns of large-scale clear and cloudy regions across the surface, we find that the typical physical scale of cloud cover disruption should be smaller than 5-8% of the disk area for four of our targets. The possible variations seen in SDSS J1254-0122 are not strong enough to allow us to confirm the cloud breaking hypothesis.

We highlight the importance of gaseous TiO and VO opacity on the highly irradiated close-in giant planets. The atmospheres of these planets naturally fall into two classes that are somewhat analogous to the M- and L-type dwarfs. Those that are warm enough to have appreciable opacity due to TiO and VO gases we term the ``pM Class'' planets, and those that are cooler we term ``pL Class'' planets. We calculate model atmospheres for these planets, including pressure-temperature profiles, spectra, and characteristic radiative time constants. We show that pM Class planets have hot stratospheres $\sim$2000 K and appear ``anomalously'' bright in the mid infrared secondary eclipse, as was recently found for planets HD 149026b and HD 209458b. This class of planets absorbs incident flux and emits thermal flux from high in their atmospheres. Consequently, they will have large day/night temperature contrasts and negligible phase shifts between orbital phase and thermal emission light curves, because radiative timescales are much shorter than possible dynamical timescales. The pL Class planets absorb incident flux deeper in the atmosphere where atmospheric dynamics will more readily redistribute absorbed energy. This will lead to cooler day sides, warmer night sides, and larger phase shifts in thermal emission light curves. Around a Sun-like primary this boundary occurs at $\sim$0.04-0.05 AU. The eccentric transiting planets HD 147506b and HD 17156b alternate between the classes. Thermal emission in the optical from pM Class planets is significant red-ward of 400 nm, making these planets attractive targets for optical detection. The difference in the observed day/night contrast between ups Andromeda b (pM Class) and HD 189733b (pL Class) is naturally explained in this scenario. (Abridged.)

A test case comparison is presented for different dust cloud model approaches applied in brown dwarfs and giant gas planets. We aim to achieve more transparency in evaluating the uncertainty inherent to theoretical modelling. We show in how far model results for characteristic dust quantities vary due to different assumptions. We also demonstrate differences in the spectral energy distributions resulting from our individual cloud modelling in 1D substellar atmosphere simulations

We present an analysis of the 0.95-14.5 micron spectral energy distributions of nine field ultracool dwarfs with spectral types ranging from L1 to T4.5. Effective temperatures, gravities, and condensate cloud sedimentation efficiencies are derived by comparing the data to synthetic spectra computed from atmospheric models that self-consistently include the formation of condensate clouds. Derived effective temperatures decrease steadily through the L1 to T4.5 spectral types and we confirm that the effective temperatures of ultracool dwarfs at the L/T transition are nearly constant, decreasing by only ~200 K from spectral types L7.5 to T4.5. The two objects in our sample with very red J-Ks colors are best fitted with synthetic spectra that have thick clouds which hints at a possible correlation between the near-infrared colors of L dwarfs and the condensate cloud properties. The fits to the two T dwarfs in our sample (T2 and T4.5) also suggest that the clouds become thinner in this spectral class, in agreement with previous studies. Restricting the fits to narrower wavelength ranges (i.e., individual photometric bands) almost always yields excellent agreement between the data and models. Limitations in our knowledge of the opacities of key absorbers such as FeH, VO, and CH4 at certain wavelengths remain obvious, however. The effective temperatures obtained by fitting the narrower wavelength ranges can show a large scatter compared to the values derived by fitting the full spectral energy distributions; deviations are typically ~200 K and in the worst cases, up to 700 K.

Opacities and chemical abundance data are crucial ingredients of ultracool dwarf and extrasolar giant planet atmosphere models. We report here on the detailed sources of molecular opacity data employed by our group for this application. We also present tables of Rosseland and Planck mean opacities which are of use in some studies of the atmospheres, interiors, and evolution of planets and brown dwarfs. For the tables presented here we have included the opacities of important atomic and molecular species, including the alkali elements, pressure induced absorption by hydrogen, and other significant opacity sources but neglect opacity from condensates. We report for each species how we have assembled molecular line data from a combination of public databases, laboratory data that is not yet in the public databases, and our own numerical calculations. We combine these opacities with abundances computed from a chemical equilibrium model using recently revised solar abundances to compute mean opacities. The chemical equilibrium calculation accounts for the settling of condensates in a gravitational field, and is applicable to ultracool dwarf and extrasolar planetary atmospheres, but not circumstellar disks. We find that the inclusion of alkali atomic opacity substantially increases the mean opacities over those currently in the literature at densities relevant to the atmospheres and interiors of giant planets and brown dwarfs. We provide our opacity tables for public use and discuss their limitations.

Published mid infrared spectra of transiting planets HD 209458b and HD 189733b, obtained during secondary eclipse by the InfraRed Spectrograph (IRS) aboard the Spitzer Space Telescope, are predominantly featureless. In particular these flux ratio spectra do not exhibit an expected feature arising from water vapor absorption short-ward of 10 um. Here we suggest that, in the absence of flux variability, the spectral data for HD 189733b are inconsistent with 8 um-photometry obtained with Spitzer's InfraRed Array Camera (IRAC), perhaps an indication of problems with the challenging reduction of the IRS spectra. The IRAC point, along with previously published secondary eclipse photometry for HD 189733b, are in good agreement with a one-dimensional model of HD 189733b that clearly shows absorption due to water vapor in the emergent spectrum. We are not able to draw firm conclusions regarding the IRS data for HD 209458b, but spectra predicted by 1D and 3D atmosphere models fit the data adequately, without adjustment of the water abundance or reliance on cloud opacity. We argue that the generally good agreement between model spectra and IRS spectra of brown dwarfs with atmospheric temperatures similar to these highly irradiated planets lends confidence in the modeling procedure.

To aid in the physical interpretation of planetary radii constrained through observations of transiting planets, or eventually direct detections, we compute model radii of pure hydrogen-helium, water, rock, and iron planets, along with various mixtures. Masses ranging from 0.01 Earth masses to 10 Jupiter masses at orbital distances of 0.02 to 10 AU are considered. For hydrogen-helium rich planets, our models are the first to couple planetary evolution to stellar irradiation over a wide range of orbital separations (0.02 to 10 AU) through a non-gray radiative-convective equilibrium atmosphere model. Stellar irradiation retards the contraction of giant planets, but its effect is not a simple function of the irradiation level: a planet at 1 AU contracts as slowly as a planet at 0.1 AU. For hydrogen-helium planets, we consider cores up to 90% of the total planet mass, comparable to those of Uranus and Neptune. If "hot Neptunes" have maintained their original masses and are not remnants of more massive planets, radii of 0.30-0.45 times Jupiter's radius are expected. Water planets are ~40-50% larger than rocky planets, independent of mass. Finally, we provide tables of planetary radii at various ages and compositions, and for ice-rock-iron planets we fit our results to analytic functions, which will allow for quick composition estimates, given masses and radii, or mass estimates, given only planetary radii. These results will assist in the interpretation of observations for both the current transiting planet surveys as well as upcoming space missions, including CoRoT and Kepler.

We report on the development of a technique for precise radial-velocity measurements of cool stars and brown dwarfs in the near infrared. Our technique is analogous to the Iodine (I2) absorption cell method that has proven so successful in the optical regime. We rely on telluric CH4 absorption features to serve as a wavelength reference, relative to which we measure Doppler shifts of the CO and H2O features in the spectra of our targets. We apply this technique to high-resolution (R~50,000) spectra near 2.3 micron of nine L dwarfs taken with the Phoenix instrument on Gemini-South and demonstrate a typical precision of 300 m/s. We conduct simulations to estimate our expected precision and show our performance is currently limited by the signal-to-noise of our data. We present estimates of the rotational velocities and systemic velocities of our targets. With our current data, we are sensitive to companions with M sin i > 2MJ in orbits with periods less than three days. We identify no companions in our current data set. Future observations with improved signal-to-noise should result in radial-velocity precision of 100 m/s for L dwarfs.

We use new and published near-IR spectra, with synthetic spectra, to derive physical properties of three of the latest-type T dwarfs. A new R~1700 spectrum of the T7.5 dwarf HD 3651B, with existing data, allows a detailed comparison to the well-studied and very similar dwarf, Gl 570D. We find that HD 3651B has both higher gravity and metallicity than Gl 570D, with Teff=820-830K, log g= 5.4-5.5, [m/H]= +0.2 and Kzz=10^4cm^2/s. Its age is 8-12 Gyr and its implied mass is 60-70 M_Jup. We perform a similar analyis of the T8 and T7.5 dwarfs 2MASS J09393548-2448279 and 2MASS J11145133-2618235 using published data, comparing them to the well-studied T8, 2MASS J04151954-0935066. We find that the two dwarfs have the same Teff as the reference dwarf, and similar gravities, but lower metallicities. The parameters are Teff=725-775K and [m/H]= -0.3; log g=5.3-5.45 for 2MASS J09393548-2448279 and log g=5.0-5.3 for 2MASS J11145133- 261823. The age and mass are ~10Gyr and 60M_Jup for 2MASS J09393548-2448279, and ~5 Gyr and 40M_Jup for 2MASS J11145133-261823. A serious limitation is the incompleteness of the line lists of CH4 and NH3 at lambda <1.7um. Spectra of Saturn and Jupiter, and of laboratory CH4 and NH3 gas, suggest that NH3 features in the Y- and J-bands may be useful as indicators of the next cooler spectral type, and not features in the H- and K-bands as previously thought. However large uncertainties remain, as the abundance of NH3 is likely to be significantly below the chemical equilibrium value, and inclusion of laboratory NH3 opacities predicts band shapes that are discrepant with existing data. It is possible that the T spectral class will have to be extended to low temperatures around 400K, when water clouds condense in the atmosphere [abridged].

We present 10 - 19 um moderate resolution spectra of ten M dwarfs, one L dwarf, and two T dwarf systems obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope. The IRS allows us to examine molecular spectroscopic features/lines at moderate spectral resolution in a heretofore untapped wavelength regime. These R~600 spectra allow for a more detailed examination of clouds, non-equilibrium chemistry, as well as the molecular features of H2O, NH3, and other trace molecular species that are the hallmarks of these objects. A cloud-free model best fits our mid-infrared spectrum of the T1 dwarf epsilon Indi Ba, and we find that the NH3 feature in epsilon Indi Bb is best explained by a non-equilibrium abundance due to vertical transport in its atmosphere. We examined a set of objects (mostly M dwarfs) in multiple systems to look for evidence of emission features, which might indicate an atmospheric temperature inversion, as well as trace molecular species; however, we found no evidence of either.

We present new infrared spectra of the T8 brown dwarf 2MASS J04151954-0935066: 2.9-4.1 micron spectra obtained with the Infrared Camera and Spectrograph on the Subaru Telescope, and 5.2-14.5 micron spectra obtained with the Infrared Spectrograph on the Spitzer Space Telescope. We use these data and models to determine an accurate bolometric luminosity of log L_bol/L_sun=-5.67 and to constrain the effective temperature, gravity, mass and age to 725-775K, log g=5.00-5.37, M=33-58 M_Jupiter and age=3-10Gyr. We perform the same analysis using published 0.6-15 micron spectra for the T7.5 dwarf 2MASS J12171110-0311131, for which we find a metal-rich composition ([Fe/H]~0.3) and log L_bol/L_sun=-5.31, T_eff=850-950K, log g=4.80-5.42, M=25-66M_Jupiter and age=1-10Gyr. These luminosities and effective temperatures straddle those determined with the same method and models for Gl 570D by Saumon et al. (2006) and make 2MASS J04151954-0935066 the coolest and least luminous T dwarf with well-determined properties. We find that synthetic spectra generated by the models reproduce the observed red through mid-infrared spectra of 2MASS J04151954-0935066 and 2MASS J12171110-0311131 very well, except for known discrepancies which are most likely due to the incomplete CH4 opacities. Both objects show evidence of departures from strict chemical equilibrium and we discuss this result in the context of other late T dwarfs where disequilibrium phenomena have been observed.

We present new L' (3.75 um) photometry of six L and T dwarfs, and M' (4.70 um) photometry of ten L and T dwarfs, observed at Gemini Observatory, and new 3.55, 4.49, 5.73 and 7.87 um photometry of nine L and T dwarfs, obtained with the Spitzer Space Telescope. The sample includes unusually blue and red dwarfs from our near-infrared studies. The data are combined with published L', M' and Spitzer photometry of L and T dwarfs, and trends of colors with spectral type and other colors are examined. Model atmospheres by Marley and Saumon are used to generate synthetic colors for ranges of effective temperature, gravity, grain sedimentation efficiency, metallicity and vertical mixing efficiency. We explore how these parameters affect the mid-infrared colors of L and T dwarfs and find that the data are modelled satisfactorily only if substantial vertical mixing occurs in both L- and T-dwarf atmospheres. The location and range of the L and T dwarf sequences in IRAC color-color and color - magnitude diagrams is also only reproduced if this mixing occurs, with a range of efficiency described by K_zz ~ 10^2-10^6 cm^2/s. The colors of the unusually red dwarfs are best reproduced by non-equilibrium models with low sedimentation efficiency, i.e. thick cloud decks, and those of the unusually blue dwarfs by non-equilibrium models with high sedimentation efficiency, i.e. thin cloud decks. The K-L' and Spitzer [3.55]-[4.49] colors can be used as indicators of effective temperature for L and T dwarfs, but care must be taken to include gravity and metallicity effects for late-T dwarfs and vertical mixing for both late-L and T dwarfs.

Traditional thermal evolution models of giant planets employ arbitrary initial conditions selected more for computational expediency than physical accuracy. Since the initial conditions are eventually forgotten by the evolving planet, this approach is valid for mature planets, if not young ones. To explore the evolution at young ages of jovian mass planets we have employed model planets created by one implementation of the core accretion mechanism as initial conditions for evolutionary calculations. We find that young jovian planets are smaller, cooler, and several to 100 times less luminous than predicted by earlier models. Furthermore the time interval during which the young jupiters are fainter than expected depends on the mass of planet. Jupiter mass planets (1 M_J) align with the conventional model luminosity in as little at 20 million years, but 10 M_J planets can take up to 1 billion years to match commonly cited luminosities, given our implementation of the core accretion mechanism. If our assumptions, especially including our treatment of the accretion shock, are correct, then young jovian planets are substantially fainter at young ages than currently believed. These results have important consequences both for detection strategies and for assigning masses to young jovian planets based on observed luminosities.

We explore the infrared spectrum of a three-dimensional dynamical model of planet HD209458b as a function of orbital phase. The dynamical model predicts day-side atmospheric pressure-temperature profiles that are much more isothermal at pressures less than 1 bar than one-dimensional radiative-convective models have found. The resulting day-side thermal spectra are very similar to a blackbody, and only weak water absorption features are seen at short wavelengths. The dayside emission is consequently in significantly better agreement with ground-based and space-based secondary eclipse data than any previous models, which predict strong flux peaks and deep absorption features. At other orbital phases, absorption due to carbon monoxide and methane is also predicted. We compute the spectra under two treatments of atmospheric chemistry: one uses the predictions of equilibrium chemistry, and the other uses non-equilibrium chemistry, which ties the timescales of methane and carbon monoxide chemistry to dynamical timescales. As a function of orbital phase, we predict planet-to-star flux ratios for standard infrared bands and all Spitzer Space Telescope bands. In Spitzer bands, we predict 2-fold to 15-fold variation in planetary flux as a function of orbital phase with equilibrium chemistry, and 2-fold to 4-fold variation with non-equilibrium chemistry. Variation is generally more pronounced in bands from 3-10 $\mu$m than at longer wavelengths. The orbital phase of maximum thermal emission in infrared bands is 15--45 orbital degrees before the time of secondary eclipse. Changes in flux as a function of orbital phase for HD209458b should be observable with Spitzer, given the previously acheived observational error bars.

We have conducted a photometric monitoring program of 3 field late-L brown dwarfs looking for evidence of non-axisymmetric structure or temporal variability in their photospheres. The observations were performed using Spitzer/IRAC 4.5 and 8 micron bandpasses and were designed to cover at least one rotational period of each object. One-sigma RMS (root mean squared) uncertainties of less than 3 mmag at 4.5 micron and around 9 mmag at 8 micron were achieved. Two out of the three objects studied exhibit some modulation in their light curves at 4.5 micron - but not 8 micron - with periods of 7.4 hr and 4.6 hr and peak-to-peak amplitudes of 10 mmag and 8 mmag. Although the lack of detectable 8 micron variation suggests an instrumental origin for the detected variations, the data may nevertheless still be consistent with intrinsic variability since the shorter wavelength IRAC bandpasses probe more deeply into late L dwarf atmospheres than the longer wavelengths. A cloud feature occupying a small percentage (1-2 %) of the visible hemisphere could account for the observed amplitude of variation. If, instead, the variability is indeed instrumental in origin, then our non-variable L dwarfs could be either completely covered with clouds or objects whose clouds are smaller and uniformly distributed. Such scenarios would lead to very small photometric variations. Followup IRAC photometry at 3.6 and 5.8 micron bandpasses should distinguish between the two cases. In any event, the present observations provide the most sensitive search to date for structure in the photospheres of late-L dwarfs at mid-IR wavelengths, and our photometry provides stringent upper limits to the extent to which the photospheres of these transition L dwarfs are structured.

We present a low-resolution (R = 90), 5.5-38 micron spectral sequence of a sample of M, L, and T dwarfs obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope. The spectra exhibit prominent absorption bands of H_2O at 6.27 microns, CH_4 at 7.65 microns, and NH_3 at 10.5 microns and are relatively featureless at lambda > 15 microns. Three spectral indices that measure the strengths of these bands are presented; H_2O absorption features are present throughout the MLT sequence while the CH_4 and NH_3 bands first appear at roughly the L/T transition. Although the spectra are, in general, qualitatively well matched by synthetic spectra that include the formation of spatially homogeneous silicate and iron condensate clouds, the spectra of the mid-type L dwarfs show an unexpected flattening from roughly 9 to 11 microns. We hypothesize that this may be a result of a population of small silicate grains that are not predicted in the cloud models. The spectrum of the peculiar T6 dwarf 2MASS J0937+2931 is suppressed from 5.5-7.5 microns relative to typical T6 dwarfs and may be a consequence of its mildly metal-poor/high surface gravity atmosphere. Finally, we compute bolometric luminosities of a subsample of the M, L, and T dwarfs by combining the IRS spectra with previously published 0.6-4.1 micron spectra and find good agreement with the values of Golimowski et al. who use L'- and M'-band photometry and to account for the flux emitted at lambda > 2.5 microns.

We present the first analysis of an optical to mid-infrared spectrum of the T7.5 dwarf Gliese 570D with model atmospheres, synthetic spectra, and brown dwarf evolution sequences. We obtain precise values for the basic parameters of Gl 570D: Teff=800 - 820K, log g (cm/s^2)=5.09 - 5.23, and log L/Lsun= -5.525 to -5.551. The Spitzer IRS spectrum shows prominent features of ammonia (NH3) that can only be fitted by reducing the abundance of NH3 by about one order of magnitude from the value obtained with chemical equilibrium models. We model departures from chemical equilibrium in the atmosphere of Gl 570D by considering the kinetics of nitrogen and carbon chemistry in the presence of vertical mixing. The resulting model spectrum reproduces the data very well.

The key to understanding an extrasolar giant planet's spectrum--and hence its detectability and evolution--lies with its atmosphere. Now that direct observations of thermal emission from extrasolar giant planets are in hand, atmosphere models can be used to constrain atmospheric composition, thermal structure, and ultimately the formation and evolution of detected planets. We review the important physical processes that influence the atmospheric structure and evolution of extrasolar giant planets and consider what has already been learned from the first generation of observations and modeling. We pay particular attention to the roles of cloud structure, metallicity, and atmospheric chemistry in affecting detectable properties through Spitzer Space Telescope observations of the transiting giant planets. Our review stresses the uncertainties that ultimately limit our ability to interpret EGP observations. Finally we will conclude with a look to the future as characterization of multiple individual planets in a single stellar system leads to the study of comparative planetary architectures.

We investigate the atmosphere and interior of the new transiting planet HD 149026b, which appears to be very rich in heavy elements. We first compute model atmospheres at metallicities ranging from solar to ten times solar, and show how for cases with high metallicity or inefficient redistribution of energy from the day side, the planet may develop a hot stratosphere due to absorption of stellar flux by TiO and VO. The spectra predicted by these models are very different than cooler atmosphere models without stratospheres. The spectral effects are potentially detectable with the Spitzer Space Telescope. In addition the models with hot stratospheres lead to a large limb brightening, rather than darkening. We compare the atmosphere of HD 149026b to other well-known transiting planets, including the recently discovered HD 189733b, which we show have planet-to-star flux ratios twice that of HD 209458 and TrES-1. The methane abundance in the atmosphere of HD 189733b is a sensitive indicator of atmospheric temperature and metallicity and can be constrained with Spitzer IRAC observations. We then turn to interior studies of HD 149026b and use a grid of self-consistent model atmospheres and high-pressure equations of state for all components to compute thermal evolution models of the planet. We estimate that the mass of heavy elements within the planet is in the range of 60 to 93 M_earth. Finally, we discuss trends in the radii of transiting planets with metallicity in light of this new member of the class.

Here we show preliminary calculations of the cooling and contraction of a 2 MJ planet. These calculations, which are being extended to 1-10 MJ, differ from other published "cooling tracks" in that they include a core accretion-gas capture formation scenario, the leading theory for the formation of gas giant planets. We find that the initial post-accretionary intrinsic luminosity of the planet is ~3 times less than previously published models which use arbitrary initial conditions. These differences last a few tens of millions of years. Young giant planets are intrinsically fainter than has been previously appreciated. We also discuss how uncertainties in atmospheric chemistry and the duration of the formation time of giant planets lead to challenges in deriving planetary physical properties from comparison with tabulated model values.

We present new self-consistent atmosphere models for transiting planets TrES-1 and HD209458b. The planets were recently observed with the Spitzer Space Telescope in bands centered on 4.5 and 8.0 $\mu$m, for TrES-1, and 24 $\mu$m, for HD209458b. We find that standard solar metallicity models fit the observations for HD209458b. For TrES-1, which has an T_eff ~300 K cooler, we find that models with a metallicity 3-5 times enhanced over solar abundances can match the 1$\sigma$ error bar at 4.5 $\mu$m and 2$\sigma$ at 8.0$\mu$m. Models with solar abundances that included energy deposition into the stratosphere give fluxes that fall within the 2$\sigma$ error bars in both bands. The best-fit models for both planets assume that reradiation of absorbed stellar flux occurs over the entire planet. For all models of both planets we predict planet/star flux ratios in other Spitzer bandpasses.

We present new JHKL'M' photometry on the MKO system for a large sample of L and T dwarfs identified from SDSS and 2MASS and classified according to the scheme of Geballe et al. (2002). We have compiled a sample of 105 L and T dwarfs that are uniformly classified and measured on a single photometric system. The scattered JHK spectral indices and colors of L dwarfs are likely caused by variations in the altitudes, distributions, and thicknesses of condensate clouds. Scatter in the H-K colors of late T dwarfs probably reflects the sensitivity of the K-band flux to pressure induced H_2 opacity, which itself is sensitive to surface gravity. The M' luminosities of late-T dwarfs are 1.5-2.5 times fainter than predicted under conditions of chemical equilibrium. We have computed L_bol and T_eff for 42 L and T dwarfs whose trigonometric parallaxes have been measured. We find that T_eff ~ 1450 K for types L7-T4, which supports recent models that attribute the changing JHK-band luminosities and spectral features across the L-T transition to rapid changes in the condensate clouds over a narrow range of T_eff. We compute T_eff = 600-750 K for 2MASS 0415-0935 (T9), which supplants Gl 570D as the coolest known brown dwarf.

Significant populations of field L and T dwarfs are now known, and we anticipate the discovery of even cooler dwarfs by Spitzer and ground-based infrared surveys. However, as the number of known L and T dwarfs increases so does the range in their observational properties, and difficulties have arisen in interpreting the observations. Although modellers have made significant advances, the complexity of the very low temperature, high pressure, photospheres means that problems remain such as the treatment of grain condensation as well as incomplete and non-equilibrium molecular chemistry. Also, there are several parameters which control the observed spectral energy distribution - effective temperature, grain sedimentation efficiency, metallicity and gravity - and their effects are not well understood. In this paper, based on a splinter session, we discuss classification schemes for L and T dwarfs, their dependency on wavelength, and the effects of the parameters T_eff, f_sed, [m/H] and log g on optical and infrared spectra. We will discuss the various hypotheses that have been presented for the transition from the dusty L types to the clear atmosphere T types. We conclude with a brief discussion of the spectral class beyond T. Authors of each Section are identified by their initials.

While the precise mechanism responsible for the L to T dwarf transition remains unclear, it is clearly caused by changing cloud characteristics. Here we briefly review data relevant to understanding the nature of the transition and argue that changing atmospheric dynamics produce the transition by opening holes through the global iron and silicate cloud decks. Other possibilities, such as a sudden vertical collapse in these cloud decks are also considered. Any acceptable model of the L to T transition must ultimately connect changing cloud properties to the underlying atmospheric dynamics.

We present the first mid-infrared spectra of brown dwarfs, together with observations of a low-mass star. Our targets are the M3.5 dwarf GJ 1001A, the L8 dwarf DENIS-P J0255-4700, and the T1/T6 binary system epsilon Indi Ba/Bb. As expected, the mid-infrared spectral morphology of these objects changes rapidly with spectral class due to the changes in atmospheric chemistry resulting from their differing effective temperatures and atmospheric structures. By taking advantage of the unprecedented sensitivity of the Infrared Spectrograph on the Spitzer Space Telescope we have detected the 7.8 micron methane and 10 micron ammonia bands for the first time in brown dwarf spectra.

We have compiled L' and M' photometry of 63 single and binary M, L,and T dwarfs obtained at UKIRT using the MKO filter set. This compilation includes new L' data for 8 L dwarfs and 13 T dwarfs and new M' data for 7 L dwarfs, 5 T dwarfs, and Gl 229A. We compute L_bol, BC_K, and T_eff for 42 dwarfs whose flux-calibrated JHK spectra, L' photometry, and trigonometric parallaxes are available, and we estimate these quantities for 9 other dwarfs whose parallaxes and flux-calibrated spectra have been obtained. BC_K is a well-behaved function of near-IR spectral type with a dispersion of ~0.1 mag for types M6-T5; it is significantly more scattered for types T5-T9. T_eff declines steeply and monotonically for types M6-L7 and T4-T9, but is nearly constant at ~1450 K for types L7-T4 with assumed ages of ~3 Gyr. This constant T_eff is evidenced by nearly unchanging values of L'-M' between types L6 and T3. We compare the observed K, L', and M' luminosities of L and T dwarfs in our sample with those predicted by model atmospheres with varying surface gravities and sedimentation efficiencies. The models indicate that the L3-T4.5 dwarfs generally have higher gravities (log g = 5.0-5.5) than the T6-T9 dwarfs (log g = 4.5-5.0). The predicted M' luminosities of T dwarfs are 1.5-2.5 times larger than those derived empirically for the T dwarfs in our sample. This discrepancy is attributed to absorption at 4.5-4.9 um by CO, which is not expected under the condition of chemical equilibrium assumed in the models. We determine that the L3 dwarf Kelu-1 and the T0 dwarf SDSS J0423-0414 are probable binary systems. We compute log(L_bol/L_sun) = -5.73 +/- 0.05 and T_eff = 600-750 K for the T9 dwarf 2MASS J0415-0935, which supplants Gl 570D as the least luminous and coolest brown dwarf presently known.

We present JHK photometry and spectroscopy for a sample of newly identified L and T dwarfs from the Sloan Digital Sky Survey (SDSS) and known dwarfs from the SDSS and the Two Micron All Sky Survey (2MASS). Spectral classification has been carried out using four infrared indices; we identify 9 new L8-9.5 dwarfs and 14 new T dwarfs from SDSS, including the latest SDSS dwarf yet found, the T7 SDSS J1758+4633. We classify 2MASS J0415-0935 as T9, the latest and coolest dwarf found to date. We combine these results with previously published data to produce a sample of 59 L dwarfs and 42 T dwarfs. We compare the near-IR colors and absolute magnitudes of dwarfs near the L-T transition with predictions made by models of the distribution and evolution of photospheric condensates. There is scatter in the spectral indices for L dwarfs, suggesting that these indices are probing different levels of the atmosphere. The near-IR colors of the L dwarfs also scatter, likely due to variations in the altitudes, spatial distributions and thicknesses of the condensate clouds. We identify a group of late L dwarfs that are blue for their spectral type and that have enhanced FeH, H2O and KI absorption, possibly due to an unusually small amount of condensates. The scatter seen in the H-K color for late T dwarfs can be reproduced by models with a range in surface gravity. The variation is probably due to the effect on the K-band flux of pressure-induced H2 opacity. The correlation of H-K color with gravity is supported by the strengths of the J-band KI doublet. Gravity is closely related to mass for field T dwarfs and the gravities implied by the H-K colors indicate that the T dwarfs in our sample have masses in the range 15-75 M_Jupiter, except for SDSS J1110+0116 which possibly has a mass of 10-15 M_Jupiter.

We present an analysis of brown dwarf model spectra in the mid-infrared spectral region (5 - 20 microns), in anticipation of data obtained with the Space Infrared Telescope Facility. The mid-infrared spectra of brown dwarfs are in several ways simpler than those in the near-infrared and yet provide powerful diagnostics of brown dwarf atmospheric physics and chemistry, especially when combined with ground-based data. We discuss the possibility of detection of new molecular species and of the silicate cloud, predict strong observational diagnostics for non-equilibrium chemistry between CO and CH4, and N2 and NH3, and speculate on the possibility of discovering brown dwarf stratospheres.

ESPI has been proposed for direct imaging and spectral analysis of giant planets orbiting solar-type stars. ESPI extends the concept suggested by Nisenson and Papaliolios (2001) for a square aperture apodized telescope that has sufficient dynamic range to directly detect exo-planets. With a 1.5 M square mirror, ESPI can deliver high dynamic range imagery as close as 0.3 arcseconds to bright sources, permitting a sensitive search for exoplanets around nearby stars and a study of their characteristics in reflected light.

A sophisticated approach to condensate opacity is required to properly model the atmospheres of L and T dwarfs. Here we review different models for the treatment of condensates in brown dwarf atmospheres. We conclude that models which include both particle sedimentation and upwards transport of condensate (both gas and particles) provide the best fit for the L dwarf colors. While a globally uniform cloud model fits the L dwarf data, it turns to the blue in J-K too slowly to fit the T dwarfs. Models which include local clearings in the global cloud deck, similar to Jupiter's prominent five-micron hot spots, better reproduce the available photometric data and also account for the observed resurgence of FeH absorption in early type T dwarfs.

Carbon monoxide and ammonia have been detected in the spectrum of Gl 229B at abundances that differ substantially from those obtained from chemical equilibrium. Vertical mixing in the atmosphere is a mechanism that can drive slowly reacting species out of chemical equilibrium. We explore the effects of vertical mixing as a function of mixing efficiency and effective temperature on the chemical abundances in the atmospheres of brown dwarfs and on their spectra. The models compare favorably with the observational evidence and indicate that vertical mixing plays an important role in brown dwarf atmospheres.

We report on the results of an I-band photometric variability survey of eighteen L dwarfs. We find that seven exhibit statistically significant variations above the 95.4% confidence level with root-mean-square scatter (including photometric errors) between 0.010 and 0.083 mag. Another five targets have variability probabilities ~80%, suggesting that these are likely variable objects. Three of the variable objects display significant peaks in a CLEAN periodogram that are several times higher than the noise. The period for 2MASS 0345+25 is clearly not intrinsic to the object and can be dismissed. The periods found for 2MASS 0746+20AB and 2MASS 1300+19 are unique but longer than those periods likely from rotation velocity measurements and they do not represent periodic behavior in the light curve that persists through the entire data set. These observations suggest that we are not observing the rotation modulation of a long-lived albedo feature. Instead, rapid evolution of atmospheric features is likely causing the non-periodic variability. The remaining variable objects show no prominent features in their light curves, suggesting even more rapid evolution of atmospheric features. We argue against the existence of magnetic spots in these atmospheres and favor the idea that non-uniform condensate coverage is responsible for these variations. The magnetic Reynolds number in the atmosphere of L dwarfs is too small to support the formation of magnetic spots. In contrast, silicate and iron clouds are expected to form in the photospheres of L dwarfs. Inhomogeneities in such cloud decks and the evolution of the inhomogeneities can plausibly produce the observed photometric variations.

Clouds of metal-bearing condensates play a critical role in shaping the emergent spectral energy distributions of the coolest classes of low-mass stars and brown dwarfs, L and T dwarfs. Because condensate clouds in planetary atmospheres show distinct horizontal structure, we have explored a model for partly cloudy atmospheres in brown dwarfs. Our model successfully reproduces the colors and magnitudes of both L and T dwarfs for the first time, including the unexpected brightning of the early- and mid-type T dwarfs at J-band, provided that clouds are rapidly removed from the photosphere at T_eff ~ 1200 K. The clearing of cloud layers also explains the surprising persistence and strengthening of gaseous FeH bands in early- and mid-type T dwarfs. The breakup of cloud layers is likely driven by convection in the troposphere, analogous to phenomena observed on Jupiter. Our results demonstrate that planetary-like atmospheric dynamics must be considered when examining the evolution of free-floating brown dwarfs.

The optical and infrared colors of L and T dwarfs are sensitive to cloud sedimentation and chemical equilibrium processes in their atmospheres. The i'-z' vs. J-K color-color diagram provides a window into diverse atmospheric processes mainly because different chemical processes govern each color, and cloud opacity largely affects J-K but not i'-z'. Using theoretical atmosphere models that include for the first time a self-consistent treatment of cloud formation, we present an interpretation of the i'-z' vs. J-K color trends of known L and T dwarfs. We find that the i'-z' color is extremely sensitive to chemical equilibrium assumptions: chemical equilibrium models accounting for cloud sedimentation predict redder i'-z' colors--by up to 2 magnitudes--than models that neglect sedimentation. We explore the previously known J-K color trends where objects first become redder, then bluer with decreasing effective temperature. Only models that include sedimentation of condensates are able to reproduce these trends. We find that the exact track of a cooling brown in J-K (and i'-z') is very sensitive to the details of clouds, in particular to the efficiency of sedimentation of condensates in its atmosphere. We also find that clouds still affect the strength of the J, H, and K band fluxes of even the coolest T dwarfs. In addition, we predict the locus in the i'-z' vs. J-K color-color diagram of brown dwarfs cooler than yet discovered.

We present K- and L-band photometry obtained with the Keck I telescope for a representative sample of L and T dwarfs. These observations were motivated in part by the dominant role water and methane play in shaping the flux near 2 and 3 microns and by the potential use of these bands as indicators of spectral class in the infrared. In addition, these observations aid the determination of the bolometric luminosity of L and T dwarfs. Here we report the K, L' and Ls magnitudes of our objects and the trends observed in the (K-L') and (K-Ls) colors as a function of L- and T-dwarf spectral class. We compare these colors with theoretical models, derive a relationship between effective temperature and L-spectral class, and compare our temperature estimates with others.

The recent discovery of photometric variations in L dwarfs has opened a discussion on the cause of the variations. We argue against the existence of magnetic spots in these atmospheres and favor the idea that non-uniform condensate coverage (i.e. clouds) is responsible for the variations. The magnetic Reynolds number (Rm) in the atmosphere of L dwarfs, which describes how well the gas couples with the magnetic field, is too small (<<1) to support the formation of magnetic spots. In constrast silicate and iron clouds form in the photospheres of L dwarfs. Inhomogeneities in such cloud decks can plausibly produce the observed photometric variations. Further evidence in support of clouds is the tendency for variable L dwarfs to be bluer in J-Ks than the average L dwarf of a given spectral type. This color effect is expected if clear holes appear in an otherwise uniform cloud layer.

We present a method to calculate vertical profiles of particle size distributions in condensation clouds of giant planets and brown dwarfs. The method assumes a balance between turbulent diffusion and sedimentation in horizontally uniform cloud decks. Calculations for the Jovian ammonia cloud are compared with results from previous methods. An adjustable parameter describing the efficiency of sedimentation allows the new model to span the range of predictions made by previous models. Calculations for the Jovian ammonia cloud are consistent with observations. Example calculations are provided for water, silicate, and iron clouds on brown dwarfs and on a cool extrasolar giant planet. We find that precipitating cloud decks naturally account for the characteristic trends seen in the spectra and colors of L- and T-type ultracool dwarfs.

Clouds and hazes are important throughout our solar system and in the atmospheres of brown dwarfs and extrasolar giant planets. Among the brown dwarfs, clouds control the colors and spectra of the L-dwarfs; the disappearance of clouds helps herald the arrival of the T-dwarfs. The structure and composition of clouds will be among the first remote-sensing results from the direct detection of extrasolar giant planets.

We have obtained a good quality R$\sim$400 0.8-2.5 $\mu$m spectrum and accurate photometry of Gl 570D, one of the coolest and least luminous brown dwarfs currently known. The spectrum shows that Gl 570D has deeper absorptions in the strong water and methane bands at 1.12-1.17 $\mu$m, 1.33-1.45 $\mu$m, 1.62-1.88 $\mu$m, and 2.20-2.45 $\mu$m and is both bluer at J-K and redder at $K-L^{\prime}$ than previously observed T dwarfs. Data analysis using model spectra coupled with knowledge of the well-understood primary implies that for the same surface gravity, Gl 570D is about 160 K cooler than Gl 229B. For an age range of 2-5 Gyr Gl 570D has an effective temperature in the range 784-824 K, a log gravity in the range 5.00-5.27 cm s$^{-2}$, and a luminosity in the range 2.88-2.98 $\times$ 10$^{-6}$ L$_\odot$.

One of the main scientific drivers of the Space InfraRed Telescope Facility (SIRTF) is the search for brown dwarfs and extrasolar superplanets. We discuss observational strategies for identification of these objects, and conclude that an optimal strategy is a wide IRAC survey (18 deg$^2$) with a 5 $\sigma$ sensitivity of 3.9 $\mu$Jy in channel 2 (M$\sim19.1^m$). For this sensitivity, we provide estimates of the number of low mass brown dwarfs and isolated planets detected per square degree for power-law mass functions with $\alpha$=1.5, and 1.0. Shallower surveys covering a larger area are inefficient because of large overheads and detector noise. Deeper surveys covering a smaller area become more and more affected by crowding with galaxies. A survey like the one that we propose would determine the field mass function down to a few Jupiter masses through the identification of a large sample of brown dwarfs and isolated planets. The proposed SIRTF survey would also allow the first detection of ultracool substellar objects with temperatures between 700 K and 200 K. The cooling curves of substellar objects with masses less than 20 Jupiters imply that they should spend most of their lifetimes at temperatures below 700 K. Preliminary models indicate that their atmospheres could be dominated by water clouds, which would diminish their optical and near-infrared fluxes. The properties of those objects are still completely unexplored.

We have detected weak absorption features produced by the strong nu(3) methane band at 3.3 microns in two L dwarfs, 2MASSW J1507476-162738 and 2MASSI J0825196+211552, classified by Kirkpatrick et al. (2000) as spectral types L5 and L7.5 respectively. These absorptions occur in objects warmer than any in which methane previously has been detected, and mark the first appearance of methane in the ultracool star-to-brown dwarf spectral sequence.

We present new, high resolution, infrared spectra of the T dwarf Gliese 229B in the J, H, and K bandpasses. We analyze each of these as well as previously published spectra to determine its metallicity and the abundances of NH3 and CO in terms of the surface gravity of Gl 229B, which remains poorly constrained. The metallicity increases with increasing gravity and is below the solar value unless Gl 229B is a high-gravity brown dwarf with log g(cgs) ~ 5.5. The NH3 abundance is determined from both the H and the K band spectra which probe two different levels in the atmosphere. We find that the abundance from the K band data is well below that expected from chemical equilibrium, which we interpret as strong evidence for dynamical transport of NH3 in the atmosphere. This is consistent with the previous detection of CO and provides additional constraints on the dynamics of the atmosphere of this T dwarf.

The recent discovery that the close-in extrasolar giant planet, HD209458b, transits its star has provided a first-of-its-kind measurement of the planet's radius and mass. In addition, there is a provocative detection of the light reflected off of the giant planet, $\tau$ Boo b. Including the effects of stellar irradiation, we estimate the general behavior of radius/age trajectories for such planets and interpret the large measured radii of HD209458b and $\tau$ Boo b in that context. We find that HD209458b must be a hydrogen-rich gas giant. Furthermore, the large radius of close-in gas giant is not due to the thermal expansion of its atmosphere, but to the high residual entropy that remains throughout its bulk by dint of its early proximity to a luminous primary. The large stellar flux does not inflate the planet, but retards its otherwise inexorable contraction from a more extended configuration at birth. This implies either that such a planet was formed near its current orbital distance or that it migrated in from larger distances ($\geq$0.5 A.U.), no later than a few times $10^7$ years of birth.

We identify the pressure--broadened red wings of the saturated potassium resonance lines at 7700 \AA as the source of anomalous absorption seen in the near-infrared spectra of Gliese 229B and, by extension, of methane dwarfs in general. This conclusion is supported by the recent work of Tsuji {\it et al.} 1999, though unlike them we find that dust need not be invoked to explain the spectra of methane dwarfs shortward of 1 micron. We find that a combination of enhanced alkali abundances due to rainout and a more realistic non-Lorentzian theory of resonant line shapes may be all that is needed to properly account for these spectra from 0.5 to 1.0 microns. The WFPC2 $I$ measurement of Gliese 229B is also consistent with this theory. Furthermore, a combination of the blue wings of this K I resonance doublet, the red wings of the Na D lines at 5890 \AA, and, perhaps, the Li I line at 6708 \AA can explain in a natural way the observed WFPC2 $R$ band flux of Gliese 229B. Hence, we conclude that the neutral alkali metals play a central role in the near-infrared and optical spectra of methane dwarfs and that their lines have the potential to provide crucial diagnostics of brown dwarfs. We speculate on the systematics of the near-infrared and optical spectra of methane dwarfs, for a given mass and composition, that stems from the progressive burial with decreasing \teff of the alkali metal atoms to larger pressures and depths.

The atmospheres of extrasolar giant planets are modeled with various effective temperatures and gravities, with and without clouds. Bond albedos are computed by calculating the ratio of the flux reflected by a planet (integrated over wavelength) to the total stellar flux incident on the planet. This quantity is useful for estimating the effective temperature and evolution of a planet. We find it is sensitive to the stellar type of the primary. For a 5 M_Jup planet the Bond albedo varies from 0.4 to 0.3 to 0.06 as the primary star varies from A5V to G2V to M2V in spectral type. It is relatively insensitive to the effective temperature and gravity for cloud--free planets. Water clouds increase the reflectivity of the planet in the red, which increases the Bond albedo. The Bond albedo increases by an order of magnitude for a 13 M_Jup planet with an M2V primary when water clouds are present. Silicate clouds, on the other hand, can either increase or decrease the Bond albedo, depending on whether there are many small grains (the former) or few large grains (the latter).

We present evolutionary sequences for Jupiter and Saturn, based on new nongray model atmospheres, which take into account the evolution of the solar luminosity and partitioning of dense components to deeper layers. The results are used to set limits on the extent to which possible interior phase separation of hydrogen and helium may have progressed in the two planets. When combined with static models constrained by the gravity field, our evolutionary calculations constrain the helium mass fraction in Jupiter to be between 0.20 and 0.27, relative to total hydrogen and helium. This is in agreement with the Galileo determination. The helium mass fraction in Saturn's atmosphere lies between 0.11 and 0.25, higher than the Voyager determination. Based on the discrepancy between the Galileo and Voyager results for Jupiter, and our models, we predict that Cassini measurements will yield a higher atmospheric helium mass fraction for Saturn relative to the Voyager value.

The reflected spectra of extrasolar giant planets are primarily influenced by Rayleigh scattering, molecular absorption, and atmospheric condensates. We present model geometric albedo and phase integral spectra and Bond albedos for planets and brown dwarfs with masses between 0.8 and 70 Jupiter masses. Rayleigh scattering predominates in the blue while molecular absorption removes most red and infrared photons. Thus cloud-free atmospheres, found on giant planets with effective temperatures exceeding about 400 K, are quite dark in reflected light beyond 0.6 microns. In cooler atmospheres first water clouds and then other condensates provide a bright reflecting layer. Only planets with cloudy atmospheres will be detectable in reflected light beyond 1 micron. Thermal emission dominates the near-infrared for warm objects with clear atmospheres. However the presence of other condensates, not considered here, may brighten some planets in reflected near-infrared light and darken them in the blue and UV. Bond albedos, the ratio of the total reflected to incident power, are sensitive to the spectral type of the primary. Most incident photons from early type stars will be Rayleigh scattered, while most incident photons from late type stars will be absorbed. The Bond albedo of a given planet thus may range from 0.4 to 0.05, depending on the primary type. Condensation of a water cloud increases the Bond albedo of a given planet by up to a factor of two. The spectra of cloudy planets are strongly influenced by poorly constrained cloud microphysical properties, particularly particle size and supersaturation. Both Bond and geometric albedos are comparatively less sensitive to variations in planet mass and effective temperature.

We summarize and extend recent work on the theory of extrasolar giant planets (EGPs) and brown dwarfs, paying particular attention to Gliese 229 B, the albedos of EGPs, the compositions of substellar atmospheres, the connections with the giant planets in the solar system, cloud physics, and non-gray spectral synthesis. The role of condensates in altering the optical spectrum of Gliese 229 B is explored, as are the systematics of the reflection spectra from extrasolar giant planets near their primaries. In addition, we discuss the role of convection and disequilibrium chemistry in explaining the anomalous detection of CO in Gliese 229 B. Throughout, we highlight the distinctive chemistry that defines this new class of objects and set goals for future study.

Since October of 1995, employing Doppler spectroscopy astronomers have discovered as many as 20 giant planets and brown dwarfs, including companions to $\tau$ Boo, 51 Peg, $\upsilon$ And, 55 Cnc, $\rho$ CrB, 70 Vir, 16 Cyg B, and 47 UMa. These discoveries have excited the planetary science community, astronomers, and the public at large. Within hours of the announcement of the planet 51 Peg b, the first direct detection of an unimpeachable brown dwarf, Gl229 B, was also announced. Gl229 B is a watershed since it has methane spectral features and a surface temperature below 1000 Kelvin, characteristics unique to objects with substellar masses. During the last two years, building upon our previous experience in the modeling of brown dwarfs and M stars, we published theoretical studies of the evolution and spectra of extrasolar giant planets. We have recently upgraded our capabilities and now generate non-gray spectral and color models of both giant planets and brown dwarfs. This theory will soon encompass objects whose effective temperatures range from 100 K to 4000 K and whose masses span three orders of magnitude. The evolutionary, spectral, and color calculations upon which we have embarked are in direct support of the searches now being planned in earnest with the HST (WFPC2, NICMOS), the IRTF, the MMT 6.5-meter upgrade, the LBT, Keck's I and II, ISO, UKIRT, NGST, the VLT, COROT, DENIS, 2MASS, and SIRTF.

Stratospheric heating, condensation, convective transport of non-equilibrium species, and deep radiative energy transport are important processes in the atmospheres of the solar jovian planets. They likely affect the atmospheres of extrasolar giant planets and brown dwarfs as well. Stratospheric temperatures control thermal fluxes in strong molecular bands, and may dramatically affect observed spectra. Condensation processes affect the appearance of an object, alter the abundances of atmospheric species, and influence the fluxes of both reflected and emitted radiation. Convection can dredge up non-thermochemical equilibrium species from the deep atmosphere to the observable atmosphere. Finally deep radiative zones can limit the effectiveness of deep convection and alter the boundary condition at which the planet radiates energy to space. Here I discuss the role these processes may play in the atmospheres of brown dwarfs and extrasolar giant planets.

We have developed a new non-gray theory of the evolution, spectra, and colors of extrasolar giant planets (EGPs) and brown dwarfs that reveals their exotic nature and uniqueness. We have discovered that the fluxes of such objects for T$_{eff}$s from 1300 K to 100 K can be spectacularly higher in the near infrared bands than black body values and that their infrared colors are anomalously blue. As a consequence, EGPs and brown dwarfs reside in hitherto unoccupied realms of the H--R diagram and may be more easily found with current and planned telescopes than previously imagined.

We present the results of a new series of non-gray calculations of the atmospheres, spectra, colors, and evolution of extrasolar giant planets (EGPs) and brown dwarfs for effective temperatures below 1300 K. This theory encompasses most of the mass/age parameter space occupied by substellar objects and is the first spectral study down to 100 K. These calculations are in aid of the multitude of searches being conducted or planned around the world for giant planets and brown dwarfs and reveal the exotic nature of the class. Generically, absorption by H_2 at longer wavelengths and H_2O opacity windows at shorter wavelengths conspire to redistribute flux blueward. Below 1200 K, methane is the dominant carbon bearing molecule and is a universal diagnostic feature of EGP and brown dwarf spectra. We find that the primary bands in which to search are $Z$ (\sim 1.05 \mic), $J$ (\sim 1.2 \mic), $H$ (\sim 1.6 \mic), $K$ (\sim 2.2 \mic), $M$ (\sim 5 \mic), and $N$ (\sim 10 \mic), that enhancements of the emergent flux over blackbody values, in particular in the near infrared, can be by many orders of magnitude, and that the infrared colors of EGPs and brown dwarfs are much bluer than previously believed. In particular, relative to $J$ and $H$, the $K$ band flux is reduced by CH_4 and H_2 absorption. Furthermore, we derive that for T_{eff}s below 1200 K most or all true metals are sequestered below the photosphere, that an interior radiative zone is a generic feature of substellar objects, and that clouds of H_2O and NH_3 are formed for T_{eff}s below \sim 400 K and \sim 200 K, respectively. This study is done for solar-metallicity objects in isolation and does not include the effects of stellar insolation. Nevertheless, it is an attempt to bridge the gap between the planetary and stellar realms and to develop a non-gray theory of objects from 0.3 \mj (``saturn'') to 70 \mj ($\sim$0.07 \mo). We find that the detection ranges for brown dwarf/EGP discovery of both ground-- and space-based telescopes are larger than previously estimated.

Electron-degenerate, pressure-ionized hydrogen (usually referred to as metallic hydrogen) is the principal constituent of brown dwarfs, the long-sought objects which lie in the mass range between the lowest-mass stars (about eighty times the mass of Jupiter) and the giant planets. The thermodynamics and transport properties of metallic hydrogen are important for understanding the properties of these objects, which, unlike stars, continually and slowly cool from initial nondegenerate (gaseous) states. Within the last year, a brown dwarf (Gliese 229 B) has been detected and its spectrum observed and analyzed, and several examples of extrasolar giant planets have been discovered. The brown dwarf appears to have a mass of about forty to fifty Jupiter masses and is now too cool to be fusing hydrogen or deuterium, although we predict that it will have consumed all of its primordial deuterium. This paper reviews the current understanding of the interrelationship between its interior properties and its observed spectrum, and also discusses the current status of research on the structure of giant planets, both in our solar system and elsewhere.

The recent discoveries of giant planets around nearby stars have galvanized the planetary science community, astronomers, and the public at large. Since {\it direct} detection is now feasible, and is suggested by the recent acquisition of Gl229 B, it is crucial for the future of extrasolar planet searches that the fluxes, evolution, and physical structure of objects from Saturn's mass to 15 Juipter masses be theoretically investigated. We discuss our first attempts to explore the characteristics of extrasolar giant planets (EGPs), in aid of both NASA's and ESA's recent plans to search for such planets around nearby stars.

Theoretical spectra and evolutionary models that span the giant planet--brown dwarf continuum have been computed based on the recent discovery of the brown dwarf, Gliese 229 B. A flux enhancement in the 4--5 micron window is a universal feature from Jovian planets to brown dwarfs. We confirm the existence of methane and water in Gl 229 B's spectrum and find its mass to be 30 to 55 Jovian masses. Although these calculations focus on Gliese 229 B, they are also meant to guide future searches for extra-solar giant planets and brown dwarfs.