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.

Atmospheric retrievals are now a standard tool to analyze observations of exoplanet atmospheres. This data-driven approach quantitatively compares atmospheric models to observations in order to estimate atmospheric properties and their uncertainties. In this paper, we introduce a new retrieval package, the PHOENIX Exoplanet Retrieval Analysis (PETRA). PETRA places the PHOENIX atmosphere model in a retrieval framework, allowing us to combine the strengths of a well-tested and widely-used atmosphere model with the advantages of retrieval algorithms. We validate PETRA by retrieving on simulated data for which the true atmospheric state is known. We also show that PETRA can successfully reproduce results from previously published retrievals of WASP-43b and HD 209458b. For the WASP-43b results, we show the effect that different line lists and line profile treatments have on the retrieved atmospheric properties. Lastly, we describe a novel technique for retrieving the temperature structure and $e^{-}$ density in ultra-hot Jupiters using H$^{-}$ opacity, allowing us to probe atmospheres devoid of most molecular features with JWST.

The low-resolution transmission spectra of ultra-hot Jupiters observed shortward of 0.5 microns indicate strong absorption at short-wavelengths. Previous explanations have included scattering, photochemistry, escaping metals, and disequilibrium chemistry. In this Letter, we show that slopes and features shortward of 0.5 microns can be caused by opacity not commonly considered in atmosphere models of exoplanets but guaranteed to be present if conditions are near chemical equilibrium including Fe I, Fe II, Ti I, Ni I, Ca I, Ca II, and SiO. Even relatively trace species (e.g., Cr I and Mn I) can contribute through strong lines in the UV and blue-optical. Using the PHOENIX atmosphere model, we describe how the short-wavelength transit spectrum varies with equilibrium temperature between 1000 K and 4000 K, as well as the effect that the rainout of condensates has at these wavelengths. We define two spectral indices to quantify the strength of the NUV and blue absorption compared to that in the red-optical, finding that the NUV transit depth will significantly exceed the transit depth from Rayleigh scattering alone for all hot Jupiters down to around 1000 K. In the blue-optical, hot Jupiters warmer than 2000 K will have transit depths larger than that from Rayleigh scattering, but below 2000 K, Rayleigh scattering can dominate, if present. We further show that these spectral indices may be used to trace the effects of rainout. We then compare our simulated transit spectra to existing observations of WASP-12b, WASP-33b, WASP-76b, and WASP-121b.

Quantifying the evolution of stellar extreme ultraviolet (EUV, 100 -- 1000 $\overset{\circ}{A}$) emission is critical for assessing the evolution of planetary atmospheres and the habitability of M dwarf systems. Previous studies from the HAbitable Zones and M dwarf Activity across Time (HAZMAT) program showed the far- and near-UV (FUV, NUV) emission from M stars at various stages of a stellar lifetime through photometric measurements from the Galaxy Evolution Explorer (GALEX). The results revealed increased levels of short-wavelength emission that remain elevated for hundreds of millions of years. The trend for EUV flux as a function of age could not be determined empirically because absorption by the interstellar medium prevents access to the EUV wavelengths for the vast majority of stars. In this paper, we model the evolution of EUV flux from early M stars to address this observational gap. We present synthetic spectra spanning EUV to infrared wavelengths of 0.4 $\pm$ 0.05 M$_{\odot}$ stars at five distinct ages between 10 and 5000 Myr, computed with the PHOENIX atmosphere code and guided by the GALEX photometry. We model a range of EUV fluxes spanning two orders of magnitude, consistent with the observed spread in X-ray, FUV, and NUV flux at each epoch. Our results show that the stellar EUV emission from young M stars is 100 times stronger than field age M stars, and decreases as t$^{-1}$ after remaining constant for a few hundred million years. This decline stems from changes in the chromospheric temperature structure, which steadily shifts outward with time. Our models reconstruct the full spectrally and temporally resolved history of an M star's UV radiation, including the unobservable EUV radiation, which drives planetary atmospheric escape, directly impacting a planet's potential for habitability.

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

Recent years have seen increasing interest in the characterization of sub-Neptune sized planets because of their prevalence in the Galaxy, contrasted with their absence in our solar system. HD 97658 is one of the brightest stars hosting a planet of this kind, and we present the transmission spectrum of this planet by combining four HST transits, twelve Spitzer/IRAC transits, and eight MOST transits of this system. Our transmission spectrum has higher signal to noise ratio than that from previous works, and the result suggests that the slight increase in transit depth from wavelength 1.1 to 1.7 microns reported in previous works on the transmission spectrum of this planet is likely systematic. Nonetheless, our atmospheric modeling results are not conclusive as no model provides an excellent match to our data. Nonetheless we find that atmospheres with high C/O ratios (C/O >~ 0.8) and metallicities of >~ 100x solar metallicity are favored. We combine the mid-transit times from all the new Spitzer and MOST observations and obtain an updated orbital period of P=9.489295 +/- 0.000005 d, with a best-fit transit time center at T_0 = 2456361.80690 +/- 0.00038 (BJD). No transit timing variations are found in this system. We also present new measurements of the stellar rotation period (34 +/- 2 d) and stellar activity cycle (9.6 yr) of the host star HD 97658. Finally, we calculate and rank the Transmission Spectroscopy Metric of all confirmed planets cooler than 1000 K and with sizes between 1 and 4 R_Earth. We find that at least a third of small planets cooler than 1000 K can be well characterized using JWST, and of those, HD 97658b is ranked fifth, meaning it remains a high-priority target for atmospheric characterization.

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.

We search for evidence of the cause of the exoplanet radius gap, i.e. the dearth of planets with radii near $1.8\ R_\oplus$. If the cause was photoevaporation, the radius gap should trend with proxies for the early-life high-energy emission of planet-hosting stars. If, alternatively, the cause was core-powered mass loss, no such trends should exist. Critically, spurious trends between the radius gap and stellar properties arise from an underlying correlation with instellation. After accounting for this underlying correlation, we find no trends remain between the radius gap and stellar mass or present-day stellar activity as measured by near-UV emission. We dismiss the nondetection of a radius gap trend with near-UV emission because present-day near-UV emission is unlikely to trace early-life high-energy emission, but we provide a catalog of GALEX near-UV and far-UV emission measurements for general use. We interpret the nondetection of a radius gap trend with stellar mass by simulating photoevaporation with mass-dependent evolution of stellar high-energy emission. The simulation produces an undetectable trend between the radius gap and stellar mass under realistic sources of error. We conclude that no evidence, from this analysis or others in the literature, currently exists that clearly favors either photoevaporation or core powered mass loss as the primary cause of the exoplanet radius gap. However, repeating this analysis once the body of well-characterized $< 4\ R_\oplus$ planets has roughly doubled could confirm or rule out photoevaporation.

We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk.

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 $\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}$).

Correct estimates of stellar extreme ultraviolet (EUV; 100 - 1170 \AA) flux are important for studying the photochemistry and stability of exoplanet atmospheres, as EUV radiation ionizes hydrogen and contributes to the heating, expansion, and potential escape of a planet's upper atmosphere. Contamination from interstellar hydrogen makes observing EUV emission from M stars particularly difficult, and impossible past 100 pc, and necessitates other means to predict the flux in this wavelength regime. We present EUV -- infrared (100 \AA - 5.5 $\mu$m) synthetic spectra computed with the PHOENIX atmospheric code of three early M dwarf planet hosts: GJ 832 (M1.5 V), GJ 176 (M2.5 V), and GJ 436 (M3.5 V). These one-dimensional semiempirical nonlocal thermodynamic equilibrium models include simple temperature prescriptions for the stellar chromosphere and transition region, from where ultraviolet (UV; 100 - 3008 \AA) fluxes originate. We guide our models with Hubble Space Telescope far- and near-UV spectra and discuss the ability to constrain these models using Galaxy Evolution Explorer UV photometry. Our models closely reproduce the observations and predict the unobservable EUV spectrum at a wavelength resolution of < 0.1 \AA. The temperature profiles that best reproduce the observations for all three stars are described by nearly the same set of parameters, suggesting that early M type stars may have similar thermal structures in their upper atmospheres. With an impending UV observation gap and the scarcity of observed EUV spectra for stars less luminous and more distant than the Sun, upper-atmosphere models such as these are important for providing realistic spectra across short wavelengths and for advancing our understanding of the effects of radiation on planets orbiting M stars.

High-contrast medium resolution spectroscopy has been used to detect molecules such as water and carbon monoxide in the atmospheres of gas giant exoplanets. In this work, we show how it can be used to derive radial velocity (RV) measurements of directly imaged exoplanets. Improving upon the traditional cross-correlation technique, we develop a new likelihood based on joint forward modelling of the planetary signal and the starlight background (i.e., speckles). After marginalizing over the starlight model, we infer the barycentric RV of HR 8799 b and c in 2010 yielding -9.2 +- 0.5 km/s and -11.6 +- 0.5 km/s respectively. These RV measurements help to constrain the 3D orientation of the orbit of the planet by resolving the degeneracy in the longitude of ascending node. Assuming coplanar orbits for HR 8799 b and c, but not including d and e, we estimate \Omega = 89 (+27,-17) deg and i = 20.8 (4.5,-3.7) deg.

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.

The discovery of habitable zone (HZ) planets around low-mass stars has highlighted the need for a comprehensive understanding of the radiation environments in which such planets reside. Of particular importance is knowledge of the far-ultraviolet (FUV) radiation, as low-mass stars are typically much more active than solar-type stars and the proximity of their HZs can be one tenth the distance. The vast majority of the flux emitted by low-mass stars at FUV wavelengths occurs in the Lyman-$\alpha$ line at 1216 Angstroms. However, measuring a low-mass star's Lyman-$\alpha$ emission directly is almost always impossible because of the contaminating effects of interstellar hydrogen and geocoronal airglow. We observed Ross 825 (K3) and Ross 1044 (M0), two stars with exceptional radial velocities, with the STIS spectrograph aboard the Hubble Space Telescope (HST). Their radial velocities resulted in significant line shifts, allowing for a more complete view of their Lyman-$\alpha$ line profiles. We provide an updated relation between effective temperature and Lyman-$\alpha$ flux using Gaia DR2 astrometry as well as updated, model-independent relationships between Lyman-$\alpha$ flux and UV flux measurements from the Galaxy Evolution Explorer (GALEX) for low-mass stars. These new relations, in combination with GALEX's considerable spatial coverage, provide substantial predictive power for the Lyman-$\alpha$ environments for thousands of nearby, low-mass stars.

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed.

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 2.9-4.1 micron integral field spectroscopy of the L4+L4 brown dwarf binary HD 130948BC, obtained with the Arizona Lenslets for Exoplanet Spectroscopy (ALES) mode of the Large Binocular Telescope Interferometer (LBTI). The HD 130948 system is a hierarchical triple system, in which the G2V primary is joined by two co-orbiting brown dwarfs. By combining the age of the system with the dynamical masses and luminosities of the substellar companions, we can test evolutionary models of cool brown dwarfs and extra-solar giant planets. Previous near-infrared studies suggest a disagreement between HD 130948BC luminosities and those derived from evolutionary models. We obtained spatially-resolved, low-resolution (R~20) L-band spectra of HD 130948B and C to extend the wavelength coverage into the thermal infrared. Jointly using JHK photometry and ALES L-band spectra for HD 130948BC, we derive atmospheric parameters that are consistent with parameters derived from evolutionary models. We leverage the consistency of these atmospheric quantities to favor a younger age (0.50 \pm 0.07 Gyr) of the system compared to the older age (0.79 \pm 0.22 Gyr) determined with gyrochronology in order to address the luminosity discrepancy.

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.

Ultra-hot Jupiters are the most highly irradiated gas giant planets, with equilibrium temperatures from 2000 to over 4000 K. Ultra-hot Jupiters are amenable to characterization due to their high temperatures, inflated radii, and short periods, but their atmospheres are atypical for planets in that the photosphere possesses large concentrations of atoms and ions relative to molecules. Here we evaluate how the atmospheres of these planets respond to irradiation by stars of different spectral type. We find that ultra-hot Jupiters exhibit temperature inversions that are sensitive to the spectral type of the host star. The slope and temperature range across the inversion both increase as the host star effective temperature increases due to enhanced absorption at short wavelengths and low pressures. The steep temperature inversions in ultra-hot Jupiters around hot stars result in increased thermal dissociation and ionization compared to similar planets around cooler stars. The resulting increase in H$^{-}$ opacity leads to a transit spectrum that has muted absorption features. The emission spectrum, however, exhibits a large contrast in brightness temperature, a signature that will be detectable with both secondary eclipse observations and high-dispersion spectroscopy. We also find that the departures from local thermodynamic equilibrium in the stellar atmosphere can affect the degree of heating caused by atomic metals in the planet's upper atmosphere. Additionally, we further quantify the significance of heating by different opacity sources in ultra-hot Jupiter atmospheres.

Knowing the high-energy radiation environment of a star over a planet's formation and evolutionary period is critical in determining if that planet is potentially habitable and if any biosignatures could be detected, as UV radiation can severely change or destroy a planet's atmosphere. Current efforts for finding a potentially habitable planet are focused on M stars, yet K stars may offer more habitable conditions due to decreased stellar activity and more distant and wider habitable zones (HZ). While M star activity evolution has been observed photometrically and spectroscopically, there has been no dedicated investigation of K-star UV evolution. We present the first comprehensive study of the near-UV, far-UV, and X-ray evolution of K stars. We used members of young moving groups and clusters ranging in age from 10 - 625 Myr combined with field stars and their archived GALEX UV and ROSAT X-ray data to determine how the UV and X-ray radiation evolve. We find that the UV and X-ray flux incident on a HZ planet is 5 - 50 times lower than that of HZ planets around early-M stars and 50 - 1000 times lower than those around late-M stars, due to both an intrinsic decrease in K dwarf stellar activity occurring earlier than for M dwarfs and the more distant location of the K dwarf HZ.

The high energy radiation environment around M dwarf stars strongly impacts the characteristics of close-in exoplanet atmospheres, but these wavelengths are difficult to observe due to geocoronal and interstellar contamination. On account of these observational restrictions, a stellar atmosphere model may be used to compute the stellar extreme ultraviolet (EUV; 100 - 912 \AA) spectrum. We present a case study of the ultra-cool M8 dwarf star, TRAPPIST-1, which hosts seven transiting short-period terrestrial sized planets whose atmospheres will be probed by the James Webb Space Telescope. We construct semi-empirical non-LTE model spectra of TRAPPIST-1 that span EUV to infrared wavelengths (100 \AA\ - 2.5 $\mu$m) using the atmosphere code PHOENIX. These upper-atmosphere models contain prescriptions for the chromosphere and transition region and include newly added partial frequency redistribution capabilities. In the absence of broadband UV spectral observations, we constrain our models using HST Ly$\alpha$ observations from TRAPPIST-1 and GALEX FUV and NUV photometric detections from a set of old M8 stars ($>$1 Gyr). We find that calibrating the models using both data sets separately yield similar FUV and NUV fluxes, and EUV fluxes that range from (1.32 - 17.4) $\times$ 10$^{-14}$ ergs s$^{-1}$ cm$^{-2}$. The results from these models demonstrate that the EUV emission is very sensitive to the temperature structure in the transition region. Our lower activity models predict EUV fluxes similar to previously published estimates derived from semi-empirical scaling relationships, while the highest activity model predicts EUV fluxes a factor of ten higher. Results from this study support the idea that the TRAPPIST-1 habitable zone planets likely do not have much liquid water on their surfaces due to the elevated levels of high energy radiation emitted by the host star.

We present detections of methane in R of $\sim$1300, L band spectra of VHS 1256 b and PSO 318.5, two low gravity, red, late L dwarfs that share the same colors as the HR 8799 planets. These spectra reveal shallow methane features, which indicate VHS 1256 b and PSO 318.5 have photospheres that are out of chemical equilibrium. Directly imaged exoplanets usually have redder near infrared colors than the field-age population of brown dwarfs on a color magnitude diagram. These objects along the L to T transition show reduced methane absorption and evidence of photospheric clouds. Compared to the H and K bands, L band (3 micron - 4 micron) spectroscopy provides stronger constraints on the methane abundances of brown dwarfs and directly imaged exoplanets that have similar effective temperatures as L to T transition objects. When combined with near infrared spectra, the L band extends our conventional wavelength coverage, increasing our understanding of atmospheric cloud structure. Our model comparisons show relatively strong vertical mixing and photospheric clouds can explain the molecular absorption features and continua of VHS 1256 b and PSO 318.5. We also discuss the implications of this work for future exoplanet focused instruments and observations with the James Webb Space Telescope.

M stars are powerful emitters of far-ultraviolet light. Over long timescales, a significant, possibly dominant, fraction of this emission is produced by stellar flares. Characterizing this emission is critical to understanding the atmospheres of the stars producing it and the atmospheric evolution of the orbiting planets subjected to it. Ultraviolet emission is known to be elevated for several hundred million years after M stars form. Whether the same is true of ultraviolet flare activity is a key concern for the evolution of exoplanet atmospheres. Hubble Space Telescope (HST) observations by the HAZMAT program (HAbitable Zones and M dwarf Activity across Time) detected 18 flares on young (40 Myr) early M stars in the Tucana-Horologium association over 10 h of observations, ten having energy $>10^{30}$ erg. These imply flares on young M stars are 100--1000$\times$ more energetic than those occurring at the same rate on "inactive," field-age M dwarfs. However, when energies are are normalized by quiescent emission there is no statistical difference between the young and field-age samples. The most energetic flare observed, dubbed the "Hazflare," emitted an energy of $10^{32.1}$ erg in the FUV, 30$\times$ more energetic than any stellar flare previously observed in the FUV with HST's COS or STIS spectrographs. It was accompanied by $15,500\pm400$ K blackbody emission bright enough to designate it a superflare ($E>10^{33}$ erg), with an estimated bolometric energy of $10^{33.6_{-0.2}^{+0.1}}$ erg. This blackbody emitted 18$^{+2}_{-1}$% of its flux in the FUV (912--1700 \AA) where molecules are generally most sensitive to photolysis. Such hot superflares in young, early M stars could play an important role in the evolution of nascent planetary atmospheres.

We describe the detection of water vapor in the atmosphere of the transiting hot Jupiter KELT-2Ab by treating the star-planet system as a spectroscopic binary with high-resolution, ground-based spectroscopy. We resolve the signal of the planet's motion with deep combined flux observations of the star and the planet. In total, six epochs of Keck NIRSPEC $L$-band observations were obtained, and the full data set was subjected to a cross correlation analysis with a grid of self-consistent atmospheric models. We measure a radial projection of the Keplerian velocity, $K_P$, of 148 $\pm$ 7 km s$^{-1}$, consistent with transit measurements, and detect water vapor at 3.8$\sigma$. We combine NIRSPEC $L$-band data with $Spitzer$ IRAC secondary eclipse data to further probe the metallicity and carbon-to-oxygen ratio of KELT-2Ab's atmosphere. While the NIRSPEC analysis provides few extra constraints on the $Spitzer$ data, it does provide roughly the same constraints on metallicity and carbon-to-oxygen ratio. This bodes well for future investigations of the atmospheres of non-transiting hot Jupiters.

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.

The Star-Planet Activity Research CubeSat (SPARCS) is a NASA-funded astrophysics mission, devoted to the study of the ultraviolet (UV) time-domain behavior in low-mass stars. Given their abundance and size, low-mass stars are important targets in the search for habitable-zone, exoplanets. However, not enough is known about the stars flare and quiescent emission, which powers photochemical reactions on the atmospheres of possible planets. Over its initial 1-year mission, SPARCS will stare at ~10 stars in order to measure short- (minutes) and long- (months) term variability simultaneously in the near-UV (NUV - lam = 280 nm) and far-UV (FUV - lam = 162 nm). The SPARCS payload consists of a 9-cm reflector telescope paired with two high-sensitivity 2D-doped CCDs. The detectors are kept passively cooled at 238K, in order to reduce dark-current contribution. The filters have been selected to provide strong rejection of longer wavelengths, where most of the starlight is emitted. The payload will be integrated within a 6U CubeSat to be placed on a Sun-synchronous terminator orbit, allowing for long observing stares for all targets. Launch is expected to occur not earlier than October 2021.

We present the data reduction pipeline, MEAD, for Arizona Lenslets for Exoplanet Spectroscopy (ALES), the first thermal infrared integral field spectrograph designed for high-contrast imaging. ALES is an upgrade of LMIRCam, the $1-5\,\mu$m imaging camera for the Large Binocular Telescope, capable of observing astronomical objects in the thermal infrared ($3-5\,\mu$m) to produce simultaneous spatial and spectral data cubes. The pipeline is currently designed to perform $L$-band ($2.8-4.2\,\mu$m) data cube reconstruction, relying on methods used extensively by current near-infrared integral field spectrographs. ALES data cube reconstruction on each spectra uses an optimal extraction method. The calibration unit comprises a thermal infrared source, a monochromator and an optical diffuser designed to inject specific wavelengths of light into LBTI to evenly illuminate the pupil plane and ALES lenslet array with monochromatic light. Not only does the calibration unit facilitate wavelength calibration for ALES and LBTI, but it also provides images of monochromatic point spread functions (PSFs). A linear combination of these monochromatic PSFs can be optimized to fit each spectrum in the least-square sense via $\chi^2$ fitting.

The integral field spectrograph configuration of the LMIRCam science camera within the Large Binocular Telescope Interferometer (LBTI) facilitates 2 to 5 $\mu$m spectroscopy of directly imaged gas-giant exoplanets. The mode, dubbed ALES, comprises magnification optics, a lenslet array, and direct-vision prisms, all of which are included within filter wheels in LMIRCam. Our observing approach includes manual adjustments to filter wheel positions to optimize alignment, on/off nodding to track sky-background variations, and wavelength calibration using narrow band filters in series with ALES optics. For planets with separations outside our 1"x1" field of view, we use a three-point nod pattern to visit the primary, secondary and sky. To minimize overheads we select the longest exposure times and nod periods given observing conditions, especially sky brightness and variability. Using this strategy we collected several datasets of low-mass companions to nearby stars.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

Extremely irradiated hot Jupiters, exoplanets reaching dayside temperatures ${>}$2000 K, stretch our understanding of planetary atmospheres and the models we use to interpret observations. While these objects are planets in every other sense, their atmospheres reach temperatures at low pressures comparable only to stellar atmospheres. In order to understand our \textit{a priori} theoretical expectations for the nature of these objects, we self-consistently model a number of extreme hot Jupiter scenarios with the PHOENIX model atmosphere code. PHOENIX is well-tested on objects from cool brown dwarfs to expanding supernovae shells and its expansive opacity database from the UV to far-IR make PHOENIX well-suited for understanding extremely irradiated hot Jupiters. We find several fundamental differences between hot Jupiters at temperatures ${>}$2500 K and their cooler counterparts. First, absorption by atomic metals like Fe and Mg, molecules including SiO and metal hydrides, and continuous opacity sources like H$^-$ all combined with the short-wavelength output of early-type host stars result in strong thermal inversions, without the need for TiO or VO. Second, many molecular species, including H$_2$O, TiO, and VO are thermally dissociated at pressures probed by eclipse observations, biasing retrieval algorithms that assume uniform vertical abundances. We discuss other interesting properties of these objects, as well as future prospects and predictions for observing and characterizing this unique class of astrophysical object, including the first self-consistent model of the hottest known jovian planet, KELT-9b.

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.

Combining visual and spectroscopic orbits of binary stars leads to a determination of the full 3D orbit, individual masses, and distance to the system. We present a full analysis of the evolved binary system $\delta$ Delphini using astrometric data from the MIRC and PAVO instruments on the CHARA long-baseline interferometer, 97 new spectra from the Fairborn Observatory, and 87 unpublished spectra from Lick Observatory. We determine the full set of orbital elements for $\delta$ Del, along with masses of $1.78 \pm 0.07$ $M_{\odot}$ and $1.62 \pm 0.07$ $M_{\odot}$ for each component, and a distance of $63.61 \pm 0.89$ pc. These results are important in two contexts: for testing stellar evolution models and defining the detection capabilities for future planet searches. We find that the evolutionary state of this system is puzzling, as our measured flux ratios, radii, and masses imply a $\sim$ 200 Myr age difference between the components using standard stellar evolution models. Possible explanations for this age discrepancy include mass transfer scenarios with a now ejected tertiary companion. For individual measurements taken over a span of 2 years we achieve $<10$ $\mu$-arcsecond precision on differential position with 10-minute observations. The high precision of our astrometric orbit suggests that exoplanet detection capabilities are within reach of MIRC at CHARA. We compute exoplanet detection limits around $\delta$ Del, and conclude that if this precision is extended to wider systems we should be able to detect most exoplanets $>2$ M$_{J}$ on orbits $>0.75$ AU around individual components of hot binary stars via differential astrometry.

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.

GJ 436b is a prime target for understanding warm Neptune exoplanet atmospheres and a target for multiple JWST GTO programs. Here, we report the first space-based optical transmission spectrum of the planet using two HST/STIS transit observations from 0.53-1.03 microns. We find no evidence for alkali absorption features, nor evidence of a scattering slope longward of 0.53 microns. The spectrum is indicative of moderate to high metallicity (~100-1000x solar) while moderate metallicity scenarios (~100x solar) require aerosol opacity. The optical spectrum also rules out some highly scattering haze models. We find an increase in transit depth around 0.8 microns in the transmission spectra of 3 different sub-Jovian exoplanets (GJ 436b, HAT-P-26b, and GJ 1214b). While most of the data come from STIS, data from three other instruments may indicate this is not an instrumental effect. Only the transit spectrum of GJ 1214b is well fit by a model with stellar plages on the photosphere of the host star. Our photometric monitoring of the host star reveals a stellar rotation rate of 44.1 days and an activity cycle of 7.4 years. Intriguingly, GJ 436 does not become redder as it gets dimmer, which is expected if star spots were dominating the variability. These insights into the nature of the GJ 436 system help refine our expectations for future observations in the era of JWST, whose higher precision and broader wavelength coverage will shed light on the composition and structure of GJ 436b's atmosphere.

We present an optical eclipse observation of the hot Jupiter WASP-12b using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. These spectra allow us to place an upper limit of $A_g < 0.064$ (97.5% confidence level) on the planet's white light geometric albedo across 290--570 nm. Using six wavelength bins across the same wavelength range also produces stringent limits on the geometric albedo for all bins. However, our uncertainties in eclipse depth are $\sim$40% greater than the Poisson limit and may be limited by the intrinsic variability of the Sun-like host star --- the solar luminosity is known to vary at the $10^{-4}$ level on a timescale of minutes. We use our eclipse depth limits to test two previously suggested atmospheric models for this planet: Mie scattering from an aluminum-oxide haze or cloud-free Rayleigh scattering. Our stringent nondetection rules out both models and is consistent with thermal emission plus weak Rayleigh scattering from atomic hydrogen and helium. Our results are in stark contrast with those for the much cooler HD 189733b, the only other hot Jupiter with spectrally resolved reflected light observations; those data showed an increase in albedo with decreasing wavelength. The fact that the first two exoplanets with optical albedo spectra exhibit significant differences demonstrates the importance of spectrally resolved reflected light observations and highlights the great diversity among hot Jupiters.

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.

The upsilon Andromedae system was the first multi-planet system discovered orbiting a main sequence star. We describe the detection of water vapor in the atmosphere of the innermost non-transiting gas giant ups~And~b by treating the star-planet system as a spectroscopic binary with high-resolution, ground-based spectroscopy. We resolve the signal of the planet's motion and break the mass-inclination degeneracy for this non-transiting planet via deep combined flux observations of the star and the planet. In total, seven epochs of Keck NIRSPEC $L$ band observations, three epochs of Keck NIRSPEC short wavelength $K$ band observations, and three epochs of Keck NIRSPEC long wavelength $K$ band observations of the ups~And~system were obtained. We perform a multi-epoch cross correlation of the full data set with an atmospheric model. We measure the radial projection of the Keplerian velocity ($K_P$ = 55 $\pm$ 9 km/s), true mass ($M_b$ = 1.7 $^{+0.33}_{-0.24}$ $M_J$), and orbital inclination \big($i_b$ = 24 $\pm$ 4$^{\circ}$\big), and determine that the planet's opacity structure is dominated by water vapor at the probed wavelengths. Dynamical simulations of the planets in the ups~And~system with these orbital elements for ups~And~b show that stable, long-term (100 Myr) orbital configurations exist. These measurements will inform future studies of the stability and evolution of the ups~And~system, as well as the atmospheric structure and composition of the hot Jupiter.

HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission.

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 optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at H{\alpha}. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

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.

We target the thermal emission spectrum of the non-transiting gas giant HD 88133 b with high-resolution near-infrared spectroscopy, by treating the planet and its host star as a spectroscopic binary. For sufficiently deep summed flux observations of the star and planet across multiple epochs, it is possible to resolve the signal of the hot gas giant's atmosphere compared to the brighter stellar spectrum, at a level consistent with the aggregate shot noise of the full data set. To do this, we first perform a principal component analysis to remove the contribution of the Earth's atmosphere to the observed spectra. Then, we use a cross-correlation analysis to tease out the spectra of the host star and HD 88133 b to determine its orbit and identify key sources of atmospheric opacity. In total, six epochs of Keck NIRSPEC L band observations and three epochs of Keck NIRSPEC K band observations of the HD 88133 system were obtained. Based on an analysis of the maximum likelihood curves calculated from the multi-epoch cross correlation of the full data set with two atmospheric models, we report the direct detection of the emission spectrum of the non-transiting exoplanet HD 88133 b and measure a radial projection of the Keplerian orbital velocity of 40 $\pm$ 15 km/s, a true mass of 1.02$^{+0.61}_{-0.28}M_J$, a nearly face-on orbital inclination of 15${^{+6}_{-5}}^{\circ}$, and an atmosphere opacity structure at high dispersion dominated by water vapor. This, combined with eleven years of radial velocity measurements of the system, provides the most up-to-date ephemeris for HD 88133.

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 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 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.

Radial velocity and transit methods are effective for the study of short orbital period exoplanets but they hardly probe objects at large separations for which direct imaging can be used. We carried out the international deep planet survey of 292 young nearby stars to search for giant exoplanets and determine their frequency. We developed a pipeline for a uniform processing of all the data that we have recorded with NIRC2/Keck II, NIRI/Gemini North, NICI/Gemini South, and NACO/VLT for 14 years. The pipeline first applies cosmetic corrections and then reduces the speckle intensity to enhance the contrast in the images. The main result of the international deep planet survey is the discovery of the HR 8799 exoplanets. We also detected 59 visual multiple systems including 16 new binary stars and 2 new triple stellar systems, as well as 2,279 point-like sources. We used Monte Carlo simulations and the Bayesian theorem to determine that 1.05[+2.80-0.70]% of stars harbor at least one giant planet between 0.5 and 14M_J and between 20 and 300 AU. This result is obtained assuming uniform distributions of planet masses and semi-major axes. If we consider power law distributions as measured for close-in planets instead, the derived frequency is 2.30[+5.95-1.55]%, recalling the strong impact of assumptions on Monte Carlo output distributions. We also find no evidence that the derived frequency depends on the mass of the hosting star, whereas it does for close-in planets. The international deep planet survey provides a database of confirmed background sources that may be useful for other exoplanet direct imaging surveys. It also puts new constraints on the number of stars with at least one giant planet reducing by a factor of two the frequencies derived by almost all previous works.

As followup to our recent detection of a pre-transit signal around HD 189733 b, we obtained full pre-transit phase coverage of a single planetary transit. The pre-transit signal is again detected in the Balmer lines but with variable strength and timing, suggesting that the bow shock geometry reported in our previous work does not describe the signal from the latest transit. We also demonstrate the use of the Ca II H and K residual core flux as a proxy for the stellar activity level throughout the transit. A moderate trend is found between the pre-transit absorption signal in the 2013 data and the Ca II H flux. This suggests that some of the 2013 pre-transit hydrogen absorption can be attributed to varying stellar activity levels. A very weak correlation is found between the Ca II H core flux and the Balmer line absorption in the 2015 transit, hinting at a smaller contribution from stellar activity compared to the 2013 transit. We simulate how varying stellar activity levels can produce changes in the Balmer line transmission spectra. These simulations show that the strength of the 2013 and 2015 pre-transit signals can be reproduced by stellar variability. If the pre-transit signature is attributed to circumplanetary material, its evolution in time can be described by accretion clumps spiraling towards the star, although this interpretation has serious limitations. Further high-cadence monitoring at H-alpha is necessary to distinguish between true absorption by transiting material and short-term variations in the stellar activity level.

We present new astrometric measurements from our ongoing monitoring campaign of the HR 8799 directly imaged planetary system. These new data points were obtained with NIRC2 on the W.M. Keck II 10 meter telescope between 2009 and 2014. In addition, we present updated astrometry from previously published observations in 2007 and 2008. All data were reduced using the SOSIE algorithm, which accounts for systematic biases present in previously published observations. This allows us to construct a self-consistent data set derived entirely from NIRC2 data alone. From this dataset, we detect acceleration for two of the planets (HR 8799b and e) at $>$3$\sigma$. We also assess possible orbital parameters for each of the four planets independently. We find no statistically significant difference in the allowed inclinations of the planets. Fitting the astrometry while forcing coplanarity also returns $\chi^2$ consistent to within 1$\sigma$ of the best fit values, suggesting that if inclination offsets of $\lesssim$20$^{o}$ are present, they are not detectable with current data. Our orbital fits also favor low eccentricities, consistent with predictions from dynamical modeling. We also find period distributions consistent to within 1$\sigma$ with a 1:2:4:8 resonance between all planets. This analysis demonstrates the importance of minimizing astrometric systematics when fitting for solutions to highly undersampled orbits.

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)

GJ 3470b is a warm Neptune-size planet transiting a M dwarf star. Like the handful of other small exoplanets for which transmission spectroscopy has been obtained, GJ 3470b exhibits a flat spectrum in the near- and mid-infrared. Recently, a tentative detection of Rayleigh scattering in its atmosphere has been reported. This signal manifests itself as an observed increase of the planetary radius as a function of decreasing wavelength in the visible. We set out to verify this detection and observed several transits of this planet with the LCOGT network and the Kuiper telescope in four different bands (Sloan g', Sloan i', Harris B and Harris V). Our analysis reveals a strong Rayleigh scattering slope, thus confirming previous results. This makes GJ 3470b the smallest known exoplanet with a detection of Rayleigh scattering. We find that the most plausible scenario is a hydrogen/helium-dominated atmosphere covered by clouds which obscure absorption features in the infrared and hazes which give rise to scattering in the visible. Our results demonstrate the feasibility of exoplanet atmospheric characterization from the ground, even with meter-class telescopes.

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.

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.

We present results from a Hubble Space Telescope (HST) program characterizing the atmospheres of the outer two planets, in the HR8799 system. The images were taken over 15 orbits in three near-infrared medium-band filters - F098M, F127M and F139M - using the Wide Field Camera 3. One of the three filters is sensitive to water absorption band inaccessible from ground-based observations, providing a unique probe of the thermal emission from the atmospheres of these young giant planets. The observations were taken at 30 different spacecraft rolls to enable angular differential imaging, and the full data set was analyzed with the Karhunen-Loeve Image Projection (KLIP) routine, an advanced image processing algorithm adapted to work with HST data. To achieve the required high contrast at sub arcsecond resolution, we utilized the pointing accuracy of HST in combination with an improved pipeline designed to combine the dithered, angular differential imaging data with an algorithm designed to both improve the image resolution and accurately measure the photometry. The results include F127M ($J$) detections of the outer planets, HR8799 b and c and the first detection of HR8799 b in the water-band (F139M) filter. The F127M photometry for HR8799 c agrees well with fitted atmospheric models resolving a long standing difficulty to model the near-IR flux for the planet consistently

Bow shocks are ubiquitous astrophysical phenomena resulting from the supersonic passage of an object through a gas. Recently, pre-transit absorption in UV metal transitions of the hot Jupiter exoplanets HD 189733b and WASP12-b have been interpreted as being caused by material compressed in a planetary bow shock. Here we present a robust detection of a time-resolved pre-transit, as well as in-transit, absorption signature around the hot Jupiter exoplanet HD 189733b using high spectral resolution observations of several hydrogen Balmer lines. The line shape of the pre-transit feature and the shape of the time series absorption provide the strongest constraints on the morphology and physical characteristics of extended structures around an exoplanet. The in-transit measurements confirm the previous exospheric H-alpha detection although the absorption depth measured here is ~50% lower. The pre-transit absorption feature occurs 125 minutes before the predicted optical transit, a projected linear distance from the planet to the stellar disk of 7.2 planetary radii. The absorption strength observed in the Balmer lines indicates an optically thick, but physically small, geometry. We model this signal as the early ingress of a planetary bow shock. If the bow shock is mediated by a planetary magnetosphere, the large standoff distance derived from the model suggests a large equatorial planetary magnetic field strength of 28 G. Better knowledge of exoplanet magnetic field strengths is crucial to understanding the role these fields play in planetary evolution and the potential development of life on planets in the habitable zone.

We used the Magellan adaptive optics system to image the 11 Myr substellar companion 1RXS 1609 B at the bluest wavelengths to date (z' and Ys). Comparison with synthetic spectra yields a higher temperature than previous studies of $T_\mathrm{eff}=2000\pm100\mathrm{K}$ and significant dust extinction of $A_V=4.5^{+0.5}_{-0.7}$ mag. Mass estimates based on the DUSTY tracks gives 0.012-0.015 Msun, making the companion likely a low-mass brown dwarf surrounded by a dusty disk. Our study suggests that 1RXS 1609 B is one of the 25% of Upper Scorpius low-mass members harboring disks, and it may have formed like a star and not a planet out at 320 AU.

Understanding the surface and atmospheric conditions of Earth-size, rocky planets in the habitable zones (HZs) of low-mass stars is currently one of the greatest astronomical endeavors. Knowledge of the planetary effective surface temperature alone is insufficient to accurately interpret biosignature gases when they are observed in the coming decades. The UV stellar spectrum drives and regulates the upper atmospheric heating and chemistry on Earth-like planets, is critical to the definition and interpretation of biosignature gases, and may even produce false-positives in our search for biologic activity. This white paper briefly describes the scientific motivation for panchromatic observations of exoplanetary systems as a whole (star and planet), argues that a future NASA UV/Vis/near-IR space observatory is well-suited to carry out this work, and describes technology development goals that can be achieved in the next decade to support the development of a UV/Vis/near-IR flagship mission in the 2020s.

Absorption lines from water, methane and carbon monoxide are detected in the atmosphere of exoplanet HR8799b. A medium-resolution spectrum presented here shows well-resolved and easily identified spectral features from all three molecules across the K band. The majority of the lines are produced by CO and H2O, but several lines clearly belong to CH4. Comparisons between these data and atmosphere models covering a range of temperatures and gravities yield log mole fractions of H2O between -3.09 and -3.91, CO between -3.30 and -3.72 and CH4 between -5.06 and -5.85. More precise mole fractions are obtained for each temperature and gravity studied. A reanalysis of H-band data, previously obtained at similar spectral resolution, results in a nearly identical water abundance as determined from the K-band spectrum. The methane abundance is shown to be sensitive to vertical mixing and indicates an eddy diffusion coefficient in the range of 10^6 to 10^8 cm^2 s^-1, comparable to mixing in the deep troposphere of Jupiter. The model comparisons also indicate a C/O between ~ 0.58 and 0.7, encompassing previous estimates for a second planet in the same system, HR8799c. Super-stellar C/O could indicate planet formation by core-accretion, however, the range of possible C/O for these planets (and the star) is currently too large to comment strongly on planet formation. More precise values of the bulk properties (e.g., effective temperature and surface gravity) are needed for improved abundance estimates.

Small, cool planets represent the typical end-products of planetary formation. Studying the archi- tectures of these systems, measuring planet masses and radii, and observing these planets' atmospheres during transit directly informs theories of planet assembly, migration, and evolution. Here we report the discovery of three small planets orbiting a bright (Ks = 8.6 mag) M0 dwarf using data collected as part of K2, the new transit survey using the re-purposed Kepler spacecraft. Stellar spectroscopy and K2 photometry indicate that the system hosts three transiting planets with radii 1.5-2.1 R_Earth, straddling the transition region between rocky and increasingly volatile-dominated compositions. With orbital periods of 10-45 days the planets receive just 1.5-10x the flux incident on Earth, making these some of the coolest small planets known orbiting a nearby star; planet d is located near the inner edge of the system's habitable zone. The bright, low-mass star makes this system an excellent laboratory to determine the planets' masses via Doppler spectroscopy and to constrain their atmospheric compositions via transit spectroscopy. This discovery demonstrates the ability of K2 and future space-based transit searches to find many fascinating objects of interest.

We used the Magellan adaptive optics (MagAO) system and its VisAO CCD camera to image the young low mass brown dwarf companion CT Chamaeleontis B for the first time at visible wavelengths. We detect it at r', i', z', and Ys. With our new photometry and Teff~2500 K derived from the shape its K-band spectrum, we find that CT Cha B has Av = 3.4+/-1.1 mag, and a mass of 14-24 Mj according to the DUSTY evolutionary tracks and its 1-5 Myr age. The overluminosity of our r' detection indicates that the companion has significant Halpha emission and a mass accretion rate ~6*10^-10 Msun/yr, similar to some substellar companions. Proper motion analysis shows that another point source within 2" of CT Cha A is not physical. This paper demonstrates how visible wavelength AO photometry (r', i', z', Ys) allows for a better estimate of extinction, luminosity, and mass accretion rate of young substellar companions.

A star's UV emission can greatly affect the atmospheric chemistry and physical properties of closely orbiting planets with the potential for severe mass loss. In particular, the Lyman-alpha emission line at 1216 Angstroms, which dominates the far-ultraviolet spectrum, is a major source of photodissociation of important atmospheric molecules such as water and methane. The intrinsic flux of Lyman-alpha, however, cannot be directly measured due to the absorption of neutral hydrogen in the interstellar medium and contamination by geocoronal emission. To date, reconstruction of the intrinsic Lyman-alpha line based on Hubble Space Telescope spectra has been accomplished for 46 FGKM nearby stars, 28 of which have also been observed by the Galaxy Evolution Explorer (GALEX). Our investigation provides a correlation between published intrinsic Lyman-alpha and GALEX far- and near-ultraviolet chromospheric fluxes for K and M stars. The negative correlations between the ratio of the Lyman-alpha to the GALEX fluxes reveal how the relative strength of Lyman-alpha compared to the broadband fluxes weakens as the FUV and NUV excess flux increase. We also correlate GALEX fluxes with the strong near-ultraviolet Mg II h+k spectral emission lines formed at lower chromospheric temperatures than Lyman-alpha. The reported correlations provide estimates of intrinsic Lyman-alpha and Mg II fluxes for the thousands of K and M stars in the archived GALEX all-sky surveys. These will constrain new stellar upper-atmosphere models for cool stars and provide realistic inputs to models describing exoplanetary photochemistry and atmospheric evolution in the absence of ultraviolet spectroscopy.

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.

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.

It is important to explore the diversity of characteristics of low-mass, low-density planets to understand the nature and evolution of this class of planets. We present a homogeneous analysis of 12 new and 9 previously published broadband photometric observations of the Uranus-sized extrasolar planet GJ 3470b, which belongs to the growing sample of sub-Jovian bodies orbiting M dwarfs. The consistency of our analysis explains some of the discrepancies between previously published results and provides updated constraints on the planetary parameters. Our data are also consistent with previous transit observations of this system. We also provide new spectroscopic measurements of GJ 3470 from 0.33 to 2.42 $\mu$$m$ to aid our analysis. We find $R_{\star}$ = 0.48$\pm$0.04 $R_{\odot}$, $M_{\star}$ = 0.51$\pm$0.06 $M_{\odot}$, and $T_{\rm eff}$ = 3652$\pm$50 K for GJ 3470, along with a rotation period of $20.70\pm{0.15}$ d and an R-band amplitude of 0.01 mag, which is small enough that current transit measurements should not be strongly affected by stellar variability. We also present the most precise orbital ephemeris for this system: T$_{o}$ = 2455983.70472$\pm$0.00021 BJD$_{TDB}$, P = 3.3366487$^{+0.0000043}_{-0.0000033}$ d, and we see no evidence for transit timing variations greater than 1 minute. Our reported planet to star radius ratio is 0.07642$\pm$0.00037. The physical parameters of this planet are $R_{p}$ = 3.88$\pm$0.32 $R_{\oplus}$, and $M_{p}$ = 13.73$\pm$1.61 $M_{\oplus}$. Because of our revised stellar parameters, the planetary radius we present is smaller than previously reported values. We also perform a second analysis of the transmission spectrum of the entire ensemble of transit observations to date, supporting the existence of a H$_{2}$ dominated atmosphere exhibiting a strong Rayleigh scattering slope.

The spectral energy distribution, variability and evolution of the high-energy radiation from an M dwarf planet host is crucial in understanding the planet's atmospheric evolution and habitability and in interpreting the planet's spectrum. The star's extreme-UV (EUV), far-UV (FUV) and near-UV (NUV) emission can chemically modify, ionize, and erode the atmosphere over time. This makes determining the lifetime exposure of such planets to stellar UV radiation critical for both the evolution of a planet's atmosphere and our potential to characterize it. Using the early M star members of nearby young moving groups (YMGs), which sample critical ages in planet formation and evolution, we measure the GALEX NUV and FUV flux as a function of age. The median UV flux remains at a "saturated" level for a few hundred million years, analogous to that observed for X-ray emission. By the age of the Hyades Cluster (650 Myr), we measure a drop in UV flux by a factor of 2-3 followed by a steep drop from old (several Gyrs) field stars. This decline in activity beyond 300 Myr follows roughly t^{-1}. Despite this clear evolution, there remains a wide range of 1-2 orders of magnitude in observed emission levels at every age. These UV data supply the much-needed constraints to M dwarf upper-atmosphere models, which will provide empirically-motivated EUV predictions and more accurate age-dependent UV spectra as inputs to planetary photochemical models.

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.

The study of exoplanetary atmospheres is one of the most exciting and dynamic frontiers in astronomy. Over the past two decades ongoing surveys have revealed an astonishing diversity in the planetary masses, radii, temperatures, orbital parameters, and host stellar properties of exoplanetary systems. We are now moving into an era where we can begin to address fundamental questions concerning the diversity of exoplanetary compositions, atmospheric and interior processes, and formation histories, just as have been pursued for solar system planets over the past century. Exoplanetary atmospheres provide a direct means to address these questions via their observable spectral signatures. In the last decade, and particularly in the last five years, tremendous progress has been made in detecting atmospheric signatures of exoplanets through photometric and spectroscopic methods using a variety of space-borne and/or ground-based observational facilities. These observations are beginning to provide important constraints on a wide gamut of atmospheric properties, including pressure-temperature profiles, chemical compositions, energy circulation, presence of clouds, and non-equilibrium processes. The latest studies are also beginning to connect the inferred chemical compositions to exoplanetary formation conditions. In the present chapter, we review the most recent developments in the area of exoplanetary atmospheres. Our review covers advances in both observations and theory of exoplanetary atmospheres, and spans a broad range of exoplanet types (gas giants, ice giants, and super-Earths) and detection methods (transiting planets, direct imaging, and radial velocity). We close with a discussion of the bright prospects for future studies of exoplanetary atmospheres.

We use high dynamic range, high-resolution L-band spectroscopy to measure the radial velocity variations of the hot Jupiter in the tau Bootis planetary system. The detection of an exoplanet by the shift in the stellar spectrum alone provides a measure of the planet's minimum mass, with the true mass degenerate with the unknown orbital inclination. Treating the tau Boo system as a high flux ratio double-lined spectroscopic binary permits the direct measurement of the planet's true mass as well as its atmospheric properties. After removing telluric absorption and cross-correlating with a model planetary spectrum dominated by water opacity, we measure a 6-sigma detection of the planet at K_p = 111 +- 5 km/s, with a 1-sigma upper limit on the spectroscopic flux ratio of 10^-4. This radial velocity leads to a planetary orbital inclination of i = 45+3-4degrees and a mass of M_P = 5.90+0.35-0.20 M_ Jup. We report the first detection of water vapor in the atmosphere of a non-transiting hot Jupiter, tau Boo b.

We report our spectroscopic investigation of the transiting ice giant GJ 3470b's atmospheric transmission, and the first results of extrasolar planet observations from the new Keck/MOSFIRE spectrograph. We measure a planet/star radius ratio of Rp/Rs = 0.0789 +/- 0.0020 in a bandpass from 2.09-2.36 micron and in six narrower bands across this wavelength range. When combined with existing broadband photometry, these measurements rule out cloud-free atmospheres in chemical equilibrium assuming either solar abundances (5.4 sigma confidence) or a moderate level of metal enrichment (50x solar abundances, 3.8 sigma), confirming previous results that such models are not representative for cool, low-mass, externally irradiated extrasolar planets. Current measurements are consistent with a flat transmission spectrum, which suggests that the atmosphere is explained by high-altitude clouds and haze, disequilibrium chemistry, unexpected abundance patterns, or the atmosphere is extremely metal-rich (>200x solar). Because GJ 3470b's low bulk density sets an upper limit on the planet's atmospheric enrichment of <300x solar, the atmospheric mean molecular weight must be <9. Thus, if the atmosphere is cloud-free its spectral features should be detectable with future observations. Transit observations at shorter wavelengths will provide the best opportunity to discriminate between plausible scenarios. We obtained optical spectroscopy with the GMOS spectrograph, but these observations exhibit large systematic uncertainties owing to thin, persistent cirrus conditions. Finally, we also provide the first detailed look at the steps necessary for well-calibrated MOSFIRE observations, and provide advice for future observations with this instrument.

We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used for a wide range of applications of spectral analysis and stellar parameter synthesis. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. The synthetic spectra cover the wavelength range from 500AA to 50.000AA with resolutions of R=500.000 in the optical and near IR, R=100.000 in the IR and a step size of 0.1AA in the UV. The parameter space covers 2.300K<=Teff<=12.000K, 0.0<=log(g)<=+6.0, -4.0<=[Fe/H]<=+1.0, and -0.2<=[alpha/Fe]<=+1.2. The library is a work in progress and we expect to extend it up to Teff=25.000 K.

Determining the atmospheric structure and chemical composition of an exoplanet remains a formidable goal. Fortunately, advancements in the study of exoplanets and their atmospheres have come in the form of direct imaging - spatially resolving the planet from its parent star - which enables high-resolution spectroscopy of self-luminous planets in Jovian-like orbits. Here, we present a spectrum with numerous, well-resolved, molecular lines from both water and carbon monoxide from a massive planet orbiting less than 40 AU from the star HR 8799. These data reveal the planet's chemical composition, atmospheric structure, and surface gravity, confirming that it is indeed a young planet. The spectral lines suggest an atmospheric carbon-to-oxygn ratio greater than the host star's, providing hints about the planet's formation.

We have monitored the position of the cool Y dwarf WISEPA J182831.08+265037.8 using a combination of ground- and space-based telescopes and have determined its distance to be 11.2$_{-1.0}^{+1.3}$ pc. Its absolute H magnitude, M$_H=22.21^{+0.25}_{-0.22}$ mag, suggests a mass in the range 0.5-20 M$_{Jup}$ for ages of 0.1-10 Gyr with an effective temperature in the range 250-400 K. The broad range in mass is due primarily to the unknown age of the object. Since the high tangential velocity of the object, 51$\pm5$ km s$^{-1}$, is characteristic of an old disk population, a plausible age range of 2-4 Gyr leads to a mass range of 3-6 M$_{Jup}$ based on fits to the (highly uncertain) COND evolutionary models. The range in temperature is due to the fact that no single model adequately represents the 1-5 $\mu$m spectral energy distribution (SED) of the source, failing by factors of up to 5 at either the short or long wavelength portions of the spectral energy distribution. The appearance of this very cold object may be affected by non-equilibrium chemistry or low temperature condensates forming clouds, two atmospheric processes that are known to be important in brown dwarf atmospheres but have proven difficult to model. Finally, we argue that there would have to be a very steep upturn in the number density of late type Y-dwarfs to account for the putative population of objects suggested by by recent microlensing observations. Whether WISE 1828+2650 sits at the low mass end of the brown dwarf population or is the first example of a large number of "free-floating" planets is not yet known.

The hot Jupiter WASP-12b is one of the largest, hottest, and best-studied extrasolar planets. We revisit our recent analysis of WASP-12b's emission spectrum in light of near-infrared spectroscopic measurements which have been claimed to support either a hydride-dominated or carbon-rich atmospheric composition. We show that this new spectrum is still consistent with a featureless blackbody, indicating a nearly isothermal photosphere on the planet's day side. Thus the ensemble of occultation measurements for WASP-12b is still insufficient to constrain the planet's atmospheric composition.

We revisit the atmospheric properties of the extremely hot Jupiter WASP-12b in light of several new developments. First, new narrowband (2.315 micron) secondary eclipse photometry that we present here, which exhibits a planet/star flux ratio of 0.45% +/- 0.06 %, corresponding to a brightness temperature of 3640 K +/- 230 K; second, recent Spitzer/IRAC and Hubble/WFC3 observations; and third, a recently observed star only 1" from WASP-12, which has diluted previous observations and which we further characterize here. We correct past WASP-12b eclipse measurements for the presence of this object, and we revisit the interpretation of WASP-12b's dilution-corrected emission spectrum. The resulting planetary emission spectrum is well-approximated by a blackbody, and consequently our primary conclusion is that the planet's infrared photosphere is nearly isothermal. Thus secondary eclipse spectroscopy is relatively ill-suited to constrain WASP-12b's atmospheric abundances, and transmission spectroscopy may be necessary to achieve this goal.

Breakthrough direct detections of planetary companions orbiting A-type stars confirm the existence of massive planets at relatively large separations, but dedicated surveys are required to estimate the frequency of similar planetary systems. To measure the first estimation of the giant exoplanetary systems frequency at large orbital separation around A-stars, we have conducted a deep-imaging survey of young (8-400 Myr), nearby (19-84 pc) A- and F-stars to search for substellar companions in the 10-300 AU range. The sample of 42 stars combines all A-stars observed in previous AO planet search surveys reported in the literature with new AO observations from VLT/NaCo and Gemini/NIRI. It represents an initial subset of the International Deep Planet Survey (IDPS) sample of stars covering M- to B-stars. The data were obtained with diffraction-limited observations in H- and Ks-band combined with angular differential imaging to suppress the speckle noise of the central stars, resulting in typical 5-sigma detection limits in magnitude difference of 12 mag at 1", 14 mag at 2" and 16 mag at 5" which is sufficient to detect massive planets. A detailed statistical analysis of the survey results is performed using Monte Carlo simulations. Considering the planet detections, we estimate the fraction of A-stars having at least one massive planet (3-14 MJup) in the range 5-320 AU to be inside 5.9-18.8% at 68% confidence, assuming a flat distribution for the mass of the planets. By comparison, the brown dwarf (15-75 MJup) frequency for the sample is 2.0-8.9% at 68% confidence in the range 5-320 AU. Assuming power law distributions for the mass and semimajor axis of the planet population, the AO data are consistent with a declining number of massive planets with increasing orbital radius which is distinct from the rising slope inferred from radial velocity (RV) surveys around evolved A-stars.

We report on a search for H-alpha absorption in four exoplanets. Strong features at H-alpha are observed in the transmission spectra of both HD 189733b and HD 209458b. We attempt to characterize and remove the effects of stellar variability in HD 189733b, and along with an empirical Monte Carlo test the results imply a statistically significant transit-dependent feature of (-8.72+/-1.48)x10^-4 integrated over a 16 Angstrom band relative to the adjacent continuum. We interpret this as the first detection of this line in an exoplanetary atmosphere. A previous detection of Ly-alpha in HD 189733b's atmosphere allows us to calculate an excitation temperature for hydrogen, T_exc=2.6x10^4 K. This calculation depends significantly on certain simplifying assumptions. We explore these assumptions and argue that T_exc is very likely much greater than the radiative equilibrium temperature (the temperature the planet is assumed to be at based on stellar radiation and the planetary distance) of HD 189733b. A large T_exc implies a very low density that is not in thermodynamic equilibrium the planet's lower atmosphere. We argue that the n=2 hydrogen required to cause H-alpha absorption in the atmosphere is created as a result of the greater UV flux at HD 189733b, which has the smallest orbit and most chromospherically active central star in our sample. Though the overall integration of HD 209458b's transmission spectrum over a wide band is consistent with zero, it contains a dramatic, statistically significant feature in the transmission spectrum with reflectional symmetry. We discuss possible physical processes that could cause this feature. Our remaining two targets (HD 147506b and HD 149026b) do not show any clear features, so we place upper limits on their H-alpha absorption levels.

We present rotational velocities for individual components of eleven very low mass (VLM) binaries with spectral types between M7 and L7.5. These results are based on observations taken with the near-infrared spectrograph, NIRSPEC, and the Keck II laser guide star adaptive optics (LGS AO) system. We find that the observed sources tend to be rapid rotators (vsini > 10 km/s), consistent with previous seeing-limited measurements of VLM objects. The two sources with the largest vsini, LP 349-25B and HD 130948C, are rotating at ~30% of their break up speed, and are among the most rapidly rotating VLM objects known. Furthermore, five binary systems, all with orbital semi-major axes <3.5 AU, have component vsini values that differ by greater than 3sigma. To bring the binary components with discrepant rotational velocities into agreement would require the rotational axes to be inclined with respect to each other, and that at least one component is inclined with respect to the orbital plane. Alternatively, each component could be rotating at a different rate, even though they have similar spectral types. Both differing rotational velocities and inclinations have implications for binary star formation and evolution. We also investigate possible dynamical evolution in the triple system HD 130948A-BC. The close binary brown dwarfs B and C have significantly different vsini values. We demonstrate that components B and C could have been torqued into misalignment by the primary star, A, via orbital precession. Such a scenario can also be applied to another triple system in our sample, GJ 569A-Bab. Interactions such as these may play an important role in the dynamical evolution of very low mass binaries. Finally, we note that two of the binaries with large differences in component vsini, LP 349-25AB and 2MASS 0746+20AB, are also known radio sources.

We report the detection of thermal emission from the hot Jupiter WASP-3b in the KS band, using a newly developed guiding scheme for the WIRC instrument at the Palomar Hale 200in telescope. Our new guiding scheme has improved the telescope guiding precision by a factor of ~5-7, significantly reducing the correlated systematics in the measured light curves. This results in the detection of a secondary eclipse with depth of 0.181%\pm0.020% (9-{\sigma}) - a significant improvement in WIRC's photometric precision and a demonstration of the capability of Palomar/WIRC to produce high quality measurements of exoplanetary atmospheres. Our measured eclipse depth cannot be explained by model atmospheres with heat redistribution but favor a pure radiative equilibrium case with no redistribution across the surface of the planet. Our measurement also gives an eclipse phase center of 0.5045\pm0.0020, corresponding to an ecos{\omega} of 0.0070\pm0.0032. This result is consistent with a circular orbit, although it also suggests the planet's orbit might be slightly eccentric. The possible non-zero eccentricity provides insight into the tidal circularization process of the star-planet system, but also might have been caused by a second low-mass planet in the system, as suggested by a previous transit timing variation study. More secondary eclipse observations, especially at multiple wavelengths, are necessary to determine the temperature-pressure profile of the planetary atmosphere and shed light on its orbital eccentricity.

We present a high-precision infrared radial velocity study of late-type stars using spectra obtained with NIRSPEC at the W. M. Keck Observatory. Radial velocity precisions of 50 m/s are achieved for old field mid-M dwarfs using telluric features for precise wavelength calibration. Using this technique, 20 young stars in the {\beta} Pic (age ~12 Myr) and TW Hya (age ~8 Myr) Associations were monitored over several years to search for low mass companions; we also included the chromospherically active field star GJ 873 (EV Lac) in this survey. Based on comparisons with previous optical observations of these young active stars, radial velocity measurements at infrared wavelengths mitigate the radial velocity noise caused by star spots by a factor of ~3. Nevertheless, star spot noise is still the dominant source of measurement error for young stars at 2.3 {\mu}m, and limits the precision to ~77 m/s for the slowest rotating stars (v sin i < 6 km/s), increasing to ~168 m/s for rapidly rotating stars (v sin i > 12 km/s). The observations reveal both GJ 3305 and TWA 23 to be single-lined spectroscopic binaries; in the case of GJ 3305, the motion is likely caused by its 0.09" companion, identified after this survey began. The large amplitude, short-timescale variations of TWA 13A are indicative of a hot Jupiter-like companion, but the available data are insufficient to confirm this. We label it as a candidate radial velocity variable. For the remainder of the sample, these observations exclude the presence of any 'hot' (P < 3 days) companions more massive than 8 MJup, and any 'warm' (P < 30 days) companions more massive than 17 MJup, on average. Assuming an edge-on orbit for the edge-on disk system AU Mic, these observations exclude the presence of any hot Jupiters more massive than 1.8 MJup or warm Jupiters more massive than 3.9 MJup.

We report the tentative detection of the near-infrared emission of the Hot Jupiter WASP-12b with the low-resolution prism on IRTF/SpeX. We find a K-H contrast color of 0.137% +/- 0.054%, corresponding to a blackbody of temperature 2400 (+1500/-500) K and consistent with previous, photometric observations. We also revisit WASP-12b's energy budget on the basis of secondary eclipse observations: the dayside luminosity is a relatively poorly constrained (2.0-4.3) x 10^30 erg/s, but this still allows us to predict a day/night effective temperature contrast of 200-1,000 K (assuming A_B=0). Thus we conclude that WASP-12b probably does not have both a low albedo and low recirculation efficiency. Our results show the promise and pitfalls of using single-slit spectrographs for characterization of extrasolar planet atmospheres, and we suggest future observing techniques and instruments which could lead to further progress. Limiting systematic effects include the use of a too-narrow slit on one night -- which observers could avoid in the future -- and chromatic slit losses (resulting from the variable size of the seeing disk) and variations in telluric transparency -- which observers cannot control. Single-slit observations of the type we present remain the best option for obtaining lambda > 1.7 micron spectra of transiting exoplanets in the brightest systems. Further and more precise spectroscopy is needed to better understand the atmospheric chemistry, structure, and energetics of this, and other, intensely irradiated planet.

We report a new detection of the H-band thermal emission of CoRoT-1b and two confirmation detections of the Ks-band thermal emission of WASP-12b at secondary eclipses. The H-band measurement of CoRoT-1b shows an eclipse depth of 0.145%\pm0.049% with a 3-{\sigma} percentile between 0.033% - 0.235%. This depth is consistent with the previous conclusions that the planet has an isother- mal region with inefficient heat transport from dayside to nightside, and has a dayside thermal inversion layer at high altitude. The two Ks band detections of WASP-12b show a joint eclipse depth of 0.299%\pm0.065%. This result agrees with the measurement of Croll & collaborators, providing independent confirmation of their measurement. The repeatability of the WASP-12b measurements also validates our data analysis method. Our measurements, in addition to a number of previous results made with other telescopes, demonstrate that ground-based observations are becoming widely available for characterization of atmospheres of hot Jupiters.

We obtained over 90 hours of spectroscopic observations of four exoplanetary systems with the Hobby-Eberly Telescope (HET). Observations were taken in transit and out of transit, and we analyzed the differenced spectra---i.e., the transmission spectra---to inspect it for absorption at the wavelengths of the neutral sodium (\ion{Na}{1}) doublet at $\lambda\lambda5889, 5895$ and neutral potassium (\ion{K}{1}) at $\lambda7698$. We used the transmission spectrum at \ion{Ca}{1} $\lambda6122$---which shows strong stellar absorption but is not an alkali metal resonance line that we expect to show significant absorption in these atmospheres---as a control line to examine our measurements for systematic errors. We use an empirical Monte Carlo method to quantity these systematic errors. In a reanalysis of the same dataset using a reduction and analysis pipeline that was derived independently, we confirm the previously seen \ion{Na}{1} absorption in HD 189733b at a level of $(-5.26\pm1.69)\times10^{-4}$ (the average value over a 12 \AA{} integration band to be consistent with previous authors). Additionally, we tentatively confirm the \ion{Na}{1} absorption seen in HD 209458b (independently by multiple authors) at a level of $(-2.63\pm0.81)\times10^{-4}$, though the interpretation is less clear. Furthermore, we find \ion{Na}{1} absorption of $(-3.16\pm2.06)\times10^{-4}$ at $<3\sigma$ in HD 149026b; features apparent in the transmission spectrum are consistent with real absorption and indicate this may be a good target for future observations to confirm. No other results (\ion{Na}{1} in HD 147506b and \ion{Ca}{1} and \ion{K}{1} in all four targets) are significant to $\geq 3\sigma$, although we observe some features that we argue are primarily artifacts.

Multi-wavelength observations/spectroscopy of exoplanetary atmospheres are the basis of the emerging exciting field of comparative exoplanetology. The HR 8799 planetary system is an ideal laboratory to study our current knowledge gap between massive field brown dwarfs and the cold 5-Gyr old Solar system planets. The HR 8799 planets have so far been imaged at J- to L-band, with only upper limits available at M-band. We present here deep high-contrast Keck II adaptive optics M-band observations that show the imaging detection of 3 of the 4 currently known HR 8799 planets. Such detections were made possible due to the development of an innovative LOCI-based background subtraction scheme that is 3 times more efficient than a classical median background subtraction for Keck II AO data, representing a gain in telescope time of up to a factor of 9. These M-band detections extend the broad band photometric coverage out to about 5 microns and provide access to the strong CO fundamental absorption band at 4.5microns. The new M-band photometry shows that the HR 8799 planets are located near the L/T-type dwarf transition, similar to what was found by other studies. We also confirm that the best atmospheric fits are consistent with low surface gravity, dusty and non-equilibrium CO/CH4 chemistry models.

The properties of 2M1207b, a young (~ 8 Myr) planet-mass companion, have lacked a satisfactory explanation for some time. The combination of low luminosity, red near-IR colors, and L-type near-IR spectrum (previously consistent with Teff ~ 1600K) implies an abnormally small radius. Early explanations for the apparent underluminosity of 2M1207b invoked an edge-on disk or the remnant of a recent protoplanetary collision. The discovery of a second planet-mass object (HR8799b) with similar luminosity and colors as 2M1207b indicates that a third explanation, one of a purely atmospheric nature, is more likely. By including clouds, non-equilibrium chemistry, and low-gravity, an atmosphere with effective temperature consistent with evolution cooling-track predictions is revealed. Consequently, 2M1207b, and others like it, requires no new physics to explain nor do they belong to a new class of objects. Instead they most likely represent the natural extension of cloudy substellar atmospheres down to low Teff and log(g). If this atmosphere only explanation for 2M1207b is correct, then very young planet-mass objects with near-IR spectra similar to field T dwarfs may be rare.

We report the detection of thermal emission at 4.5 and 8 micron from the planet WASP-17b. We used Spitzer to measure the system brightness at each wavelength during two occultations of the planet by its host star. By combining the resulting light curves with existing transit light curves and radial velocity measurements in a simultaneous analysis, we find the radius of WASP-17b to be 2.0 Rjup, which is 0.2 Rjup larger than any other known planet and 0.7 Rjup larger than predicted by the standard cooling theory of irradiated gas giant planets. We find the retrograde orbit of WASP-17b to be slightly eccentric, with 0.0012 < e < 0.070 (3 sigma). Such a low eccentricity suggests that, under current models, tidal heating alone could not have bloated the planet to its current size, so the radius of WASP-17b is currently unexplained. From the measured planet-star flux-density ratios we infer 4.5 and 8 micron brightness temperatures of 1881 +/- 50 K and 1580 +/- 150 K, respectively, consistent with a low-albedo planet that efficiently redistributes heat from its day side to its night side.

The nearby star GJ 1214 hosts a planet intermediate in radius and mass between Earth and Neptune, resulting in some uncertainty as to its nature. We have observed this planet, GJ 1214b, during transit with the high-resolution, near-infrared NIRSPEC spectrograph on the Keck II telescope, in order to characterize the planet's atmosphere. By cross-correlating the spectral changes through transit with a suite of theoretical atmosphere models, we search for variations associated with absorption in the planet atmosphere. Our observations are sufficient to rule out tested model atmospheres with wavelength-dependent transit depth variations >5e-4 over the wavelength range 2.1 - 2.4 micron. Our sensitivity is limited by variable slit loss and telluric transmission effects. We find no positive signatures but successfully rule out a number of plausible atmospheric models, including the default assumption of a gaseous, H-dominated atmosphere in chemical equilibrium. Such an atmosphere can be made consistent if the absorption due to methane is reduced. Clouds can also render such an atmosphere consistent with our observations, but only if they lie higher in the atmosphere than indicated by recent optical and infrared measurements. When taken in concert with constraints from other groups, our results support a consensus model in which the atmosphere of GJ 1214b contains significant H and He, but where methane is depleted. If this depletion is the result of photochemical processes, it may also produce a haze that suppresses spectral features in the optical.

Using the integral field spectrograph OSIRIS, on the Keck II telescope, broad near-infrared H and K-band spectra of the young exoplanet HR8799b have been obtained. In addition, six new narrow-band photometric measurements have been taken across the H and K bands. These data are combined with previously published photometry for an analysis of the planet's atmospheric properties. Thick photospheric dust cloud opacity is invoked to explain the planet's red near-IR colors and relatively smooth near-IR spectrum. Strong water absorption is detected, indicating a Hydrogen-rich atmosphere. Only weak CH4 absorption is detected at K band, indicating efficient vertical mixing and a disequilibrium CO/CH4 ratio at photospheric depths. The H-band spectrum has a distinct triangular shape consistent with low surface gravity. New giant planet atmosphere models are compared to these data with best fitting bulk parameters, Teff = 1100K +- 100 and log(g) = 3.5 +-0.5 (for solar composition). Given the observed luminosity, log(Lobs/Lsun) ~ -5.1, these values correspond to a radius of 0.75 Rjup (+0.17, -0.12) and mass ~ 0.72 Mjup (+2.6,-0.6) -- strikingly inconsistent with interior/evolution models. Enhanced metallicity (up to ~ 10 times that of the Sun) along with thick clouds and non-equilibrium chemistry are likely required to reproduce the complete ensemble of spectroscopic and photometric data and the low effective temperatures (< 1000K) required by the evolution models.

Juvenile ultracool dwarfs are late spectral type objects (later than ~M6) with ages between 10 Myr and several 100 Myr. Their age-related properties lie intermediate between very low mass objects in nearby star-forming regions (ages 1-5 Myr) and field stars and brown dwarfs that are members of the disk population (ages 1-5 Gyr). Kinematic associations of nearby young stars with ages from ~10-100 Myr provide sources for juvenile ultracool dwarfs. The lowest mass confirmed members of these groups are late-M dwarfs. Several apparently young L dwarfs and a few T dwarfs are known, but they have not been kinematically associated with any groups. Normalizing the field IMF to the high mass population of these groups suggests that more low mass (mainly late-M and possibly L dwarf) members have yet to be found. The lowest mass members of these groups, along with low mass companions to known young stars, provide benchmark objects with which spectroscopic age indicators for juvenile ultracool dwarfs can be calibrated and evaluated. In this proceeding, we summarize currently used methods for identifying juvenile ultracool dwarfs and discuss the appropriateness and reliability of the most commonly used age indicators.

With an emphasis in detecting Earth-like planets set forth by the 2010 Decadal Survey and in searching for planets around M dwarfs set forth by the 2008 Exoplanet Task Force, radial velocity surveys with infrared echelle spectrometers will have a significant impact on future exoplanet studies. Here, we present the results of an infrared radial velocity survey of a sample of 14 late-M dwarfs with the NIRSPEC echelle spectrometer on the Keck II telescope. Using telluric lines for wavelength calibration, we are able to achieve measurement precisions of 150-300 m/s over a year-long timeframe. While we require more RV epochs to determine whether most of our stars have planetary-mass companions, we have placed upper limits of 5-10 MJ on the masses of planets around a sub-set of our sample. We have also determined the rotational velocities for all the stars in the sample and offer our multi-order, high-resolution spectra over 2.0 to 2.4 micron to the modeling community to better understand the atmospheres of late-M dwarfs.

High-contrast near-infrared imaging of the nearby star HR 8799 has shown three giant planets. Such images were possible due to the wide orbits (> 25 AU) and youth (< 100 Myr) of the imaged planets, which are still hot and bright as they radiate away gravitational energy acquired during their formation. A major area of contention in the extrasolar planet community is whether outer planets (> 10 AU) more massive than Jupiter form via one-step gravitational instabilities or, rather, via a two-step process involving accretion of a core followed by accumulation of a massive outer envelope composed primarily of hydrogen and helium. Here we report the presence of a fourth planet, interior to and about the same mass as the other three. The system, with this additional planet, represents a challenge for current planet formation models as none of them can explain the in situ formation of all four planets. With its four young giant planets and known cold/warm debris belts, the HR 8799 planetary system is a unique laboratory to study the formation and evolution of giant planets at wide > 10 AU separations.

We have conducted a 4030-square-deg near-infrared proper motion survey using multi-epoch data from the Two Micron All-Sky Survey (2MASS). We find 2778 proper motion candidates, 647 of which are not listed in SIMBAD. After comparison to DSS images, we find that 107 of our proper motion candidates lack counterparts at B-, R-, and I-bands and are thus 2MASS-only detections. We present results of spectroscopic follow-up of 188 targets that include the infrared-only sources along with selected optical-counterpart sources with faint reduced proper motions or interesting colors. We also establish a set of near-infrared spectroscopic standards with which to anchor near-infrared classifications for our objects. Among the discoveries are six young field brown dwarfs, five "red L" dwarfs, three L-type subdwarfs, twelve M-type subdwarfs, eight "blue L" dwarfs, and several T dwarfs. We further refine the definitions of these exotic classes to aid future identification of similar objects. We examine their kinematics and find that both the "blue L" and "red L" dwarfs appear to be drawn from a relatively old population. This survey provides a glimpse of the kinds of research that will be possible through time-domain infrared projects such as the UKIDSS Large Area Survey, various VISTA surveys, and WISE, and also through z- or y-band enabled, multi-epoch surveys such as Pan-STARRS and LSST.

We present the first ground-based detection of thermal emission from an exoplanet in the H-band. Using HAWK-I on the VLT, we observed an occultation of WASP-19b by its G8V-type host star. WASP-19b is a Jupiter-mass planet with an orbital period of only 19 h, and thus, being highly irradiated, is expected to be hot. We measure an H-band occultation depth of (0.259 +0.046 -0.044) %, which corresponds to an H-band brightness temperature of T_H = 2580 +/- 125 K. A cloud-free model of the planet's atmosphere, with no redistribution of energy from day-side to night-side, under predicts the planet/star flux density ratio by a factor of two. As the stellar parameters, and thus the level of planetary irradiation, are well-constrained by measurement, it is likely that our model of the planet's atmosphere is too simple.

[ABRIDGED] We present the results of a 3 year monitoring program of a sample of very low mass (VLM) field binaries using both astrometric and spectroscopic data obtained in conjunction with the laser guide star adaptive optics system on the W.M. Keck II 10 m telescope. Fifteen systems have undergone sufficient orbital motion, allowing us to derive their relative orbital parameters and hence their total system mass. These measurements triple the number of masses for VLM objects. Among the 11 systems with both astrometric and spectroscopic measurements, six have sufficient radial velocity variations to allow us to obtain individual component masses. This is the first derivation of the component masses for five of these systems. Altogether, the orbital solutions of these low mass systems show a correlation between eccentricity and orbital period, consistent with their higher mass counterparts. In our primary analysis, we find that there are systematic discrepancies between our dynamical mass measurements and the predictions of theoretical evolutionary models (TUCSON and LYON) with both models either underpredicting or overpredicting the most precisely determined dynamical masses. These discrepancies are a function of spectral type, with late M through mid L systems tending to have their masses underpredicted, while one T type system has its mass overpredicted. These discrepancies imply that either the temperatures predicted by evolutionary and atmosphere models are inconsistent for an object of a given mass, or the mass-radius relationship or cooling timescales predicted by the evolutionary models are incorrect. If these spectral type trends hold into the planetary mass regime, the implication is that the masses of directly imaged extrasolar planets are overpredicted by the evolutionary models.

Tremendous progress in the science of extrasolar planets has been achieved since the discovery of a Jupiter orbiting the nearby Sun-like star 51 Pegasi in 1995. Theoretical models have now reached enough maturity to predict the characteristic properties of these new worlds, mass, radius, atmospheric signatures, and can be confronted with available observations. We review our current knowledge of the physical properties of exoplanets, internal structure and composition, atmospheric signatures, including expected biosignatures for exo-Earth planets, evolution, and the impact of tidal interaction and stellar irradiation on these properties for the short-period planets. We discuss the most recent theoretical achievements in the field and the still pending questions. We critically analyse the different solutions suggested to explain abnormally large radii of a significant fraction of transiting exoplanets. Special attention is devoted to the recently discovered transiting objects in the overlapping mass range between massive planets and low-mass brown dwarfs, stressing the ambiguous nature of these bodies, and we discuss the possible observable diagnostics to identify these two distinct populations. We also review our present understanding of planet formation and critically examine the different suggested formation mechanisms. We expect the present review to provide the basic theoretical background to capture the essential of the physics of exoplanet formation, structure and evolution, and the related observable signatures.

High contrast imaging can find and characterize gas giant planets around nearby young stars and the closest M stars, complementing radial velocity and astrometric searches by exploring orbital separations inaccessible to indirect methods. Ground-based coronagraphs are already probing within 25 AU of nearby young stars to find objects as small as ~ 3 Jupiter masses. This paper compares near-term and future ground-based capabilities with high contrast imaging modes of the James Webb Space Telescope (JWST). Monte Carlo modeling reveals that JWST can detect planets with masses as small as 0.2 MJup across a broad range of orbital separations. We present new calculations for planet brightness as a function of mass and age for specific JWST filters and extending to 0.1 MJup.

We report measurements of the thermal emission of the young and massive planet CoRoT-2b at 4.5 and 8 microns with the Spitzer Infrared Array Camera (IRAC). Our measured occultation depths are 0.510 +- 0.042 % and 0.41 +- 0.11 % at 4.5 and 8 microns, respectively. In addition to the CoRoT optical measurements, these planet/star flux ratios indicate a poor heat distribution to the night side of the planet and are in better agreement with an atmosphere free of temperature inversion layer. Still, the presence of such an inversion is not definitely ruled out by the observations and a larger wavelength coverage is required to remove the current ambiguity. Our global analysis of CoRoT, Spitzer and ground-based data confirms the large mass and size of the planet with slightly revised values (Mp = 3.47 +- 0.22 Mjup, Rp = 1.466 +- 0.044 Rjup). We find a small but significant offset in the timing of the occultation when compared to a purely circular orbital solution, leading to e cos(omega) = -0.00291 +- 0.00063 where e is the orbital eccentricity and omega is the argument of periastron. Constraining the age of the system to be at most of a few hundreds of Myr and assuming that the non-zero orbital eccentricity is not due to a third undetected body, we model the coupled orbital-tidal evolution of the system with various tidal Q values, core sizes and initial orbital parameters. For log(Q_s') = 5 - 6, our modelling is able to explain the large radius of CoRoT-2b if log(Q_p') <= 5.5 through a transient tidal circularization and corresponding planet tidal heating event. Under this model, the planet will reach its Roche limit within 20 Myr at most.

By comparing near-infrared spectra with atmosphere models, we infer the effective temperature, surface gravity, projected rotational velocity, and radial velocity for 21 very-low-mass stars and brown dwarfs. The unique sample consists of two sequences in spectral type from M6-M9, one of 5-10 Myr objects and one of >1 Gyr field objects. A third sequence is comprised of only ~M6 objects with ages ranging from <1 Myr to >1 Gyr. Spectra were obtained in the J band at medium (R~2,000) and high (R~20,000) resolutions with NIRSPEC on the Keck II telescope. Synthetic spectra were generated from atmospheric structures calculated with the PHOENIX model atmosphere code. Using multi-dimensional least-squares fitting and Monte Carlo routines we determine the best-fit model parameters for each observed spectrum and note which spectral regions provide consistent results. We identify successes in the reproduction of observed features by atmospheric models, including pressure-broadened KI lines, and investigate deficiencies in the models, particularly missing FeH opacity, that will need to be addressed in order to extend our analysis to cooler objects. The precision that can be obtained for each parameter using medium- and high- resolution near-infrared spectra is estimated and the implications for future studies of very low mass stars and brown dwarfs are discussed.

We use signal enhancement techniques and a matched filter analysis to search for the K band spectroscopic absorption signature of the close orbiting extrasolar giant planet, HD 189733b. With timeseries observations taken with NIRSPEC at the Keck II telescope, we investigate the relative abundances of H2O and carbon bearing molecules, which have now been identified in the dayside spectrum of HD 189733b. We detect a candidate planet signature with a low level of significance, close to the ~153 km/s velocity amplitude of HD 189733b. However, some systematic variations, mainly due to imperfect telluric line removal, remain in the residual spectral timeseries in which we search for the planetary signal. The robustness of our candidate signature is assessed, enabling us to conclude that it is not possible to confirm the presence of any planetary signal which appears at Fp/F* contrasts deeper than the 95.4 per cent confidence level. Our search does not enable us to detect the planet at a contrast ratio of Fp/F* = 1/1920 with 99.9 per cent confidence. We also investigate the effect of model uncertainties on our ability to reliably recover a planetary signal. The use of incorrect temperature, model opacity wavelengths and model temperature-pressure profiles have important consequences for the least squares deconvolution procedure that we use to boost the S/N ratio in our spectral timeseries observations. We find that mismatches between the empirical and model planetary spectrum may weaken the significance of a detection by ~30-60 per cent, thereby potentially impairing our ability to recover a planetary signal with high confidence.

We present new JHK spectroscopy (R ~ 5000) of GQ Lup b, acquired with the near-infrared integral field spectrograph NIFS and the adaptive optics system ALTAIR at the Gemini North telescope. Angular differential imaging was used in the J and H bands to suppress the speckle noise from GQ Lup A; we show that this approach can provide improvements in signal-to-noise ratio (S/N) by a factor of 2 - 6 for companions located at subarcsecond separations. Based on high quality observations and GAIA synthetic spectra, we estimate the companion effective temperature to Teff = 2400 +/- 100 K, its gravity to log g = 4.0 +/- 0.5, and its luminosity to log(L/L_s) = -2.47 +/- 0.28. Comparisons with the predictions of the DUSTY evolutionary tracks allow us to constrain the mass of GQ Lup b to 8 - 60 MJup, most likely in the brown dwarf regime. Compared with the spectra published by Seifahrt and collaborators, our spectra of GQ Lup b are significantly redder (by 15 - 50%) and do not show important Pa\beta emission. Our spectra are in excellent agreement with the lower S/N spectra previously published by McElwain and collaborators.

We present detailed structure and evolution calculations for the first transiting extrasolar planets discovered by the space-based CoRoT mission. Comparisons between theoretical and observed radii provide information on the internal composition of the CoRoT objects. We distinguish three different categories of planets emerging from these discoveries and from previous ground-based surveys: (i) planets explained by standard planetary models including irradiation, (ii) abnormally bloated planets and (iii) massive objects belonging to the overlapping mass regime between planets and brown dwarfs. For the second category, we show that tidal heating can explain the relevant CoRoT objects, providing non-zero eccentricities. We stress that the usual assumption of a quick circularization of the orbit by tides, as usually done in transit light curve analysis, is not justified a priori, as suggested recently by Levrard et al. (2009), and that eccentricity analysis should be carefully redone for some observations. Finally, special attention is devoted to CoRoT-3b and to the identification of its very nature: giant planet or brown dwarf ? The radius determination of this object confirms the theoretical mass-radius predictions for gaseous bodies in the substellar regime but, given the present observational uncertainties, does not allow an unambiguous identification of its very nature. This opens the avenue, however, to an observational identification of these two distinct astrophysical populations, brown dwarfs and giant planets, in their overlapping mass range, as done for the case of the 8 Jupiter-mass object Hat-P-2b. (abridged)

We present VLT eclipse photometry for the giant planet CoRoT-1b. We observed a transit in the R-band filter and an occultation in a narrow filter centered on 2.09 microns. Our analysis of this new photometry and published radial velocities, in combination with stellar-evolutionary modeling, leads to a planetary mass and radius of 1.07 (+0.13,-0.18) M_Jup and 1.45 (+0.07,-0.13) R_Jup, confirming the very low density previously deduced from CoRoT photometry. The large occultation depth that we measure at 2.09 microns (0.278 (+0.043,-0.066) %) is consistent with thermal emission and is better reproduced by an atmospheric model with no redistribution of the absorbed stellar flux to the night side of the planet.

We present a new code for the calculation of the 1D structure and synthetic spectra of accretion disks. The code is an extension of the general purpose stellar atmosphere code PHOENIX and is therefore capable of including extensive lists of atomic and molecular lines as well as dust in the calculations. We assume that the average viscosity can be represented by a critical Reynolds number in a geometrically thin disk and solve the structure and radiative transfer equations for a number of disk rings in the vertical direction. The combination of these rings provides the total disk structure and spectrum. Since the warm inner regions of protoplanetary disks show a rich molecular spectrum, they are well suited for a spectral analysis with our models. In this paper we test our code by comparing our models with high-resolution VLT CRIRES spectra of the T Tauri star GQ Lup.

Three planets have been directly imaged around the young star HR 8799. The planets are 5-13 Mjup and orbit the star at projected separations of 24-68 AU. While the initial detection occurred in 2007, two of the planets were recovered in a re-analysis of data obtained in 2004. Here we present a detection of the furthest planet of that system, HR 8799 b, in archival HST/NICMOS data from 1998. The detection was made using the locally-optimized combination of images algorithm to construct, from a large set of HST/NICMOS images of different stars taken from the archive, an optimized reference point-spread function image used to subtract the light of the primary star from the images of HR 8799. This new approach improves the sensitivity to planets at small separations by a factor of ~10 compared to traditional roll deconvolution. The new detection provides an astrometry point 10 years before the most recent observations, and is consistent with a Keplerian circular orbit with a~70 AU and low orbital inclination. The new photometry point, in the F160W filter, is in good agreement with an atmosphere model with intermediate clouds and vertical stratification, and thus suggests the presence of significant water absorption in the planet's atmosphere. The success of the new approach used here highlights a path for the search and characterization of exoplanets with future space telescopes, such as the James Webb Space Telescope or a Terrestrial Planet Finder.

Direct imaging of exoplanetary systems is a powerful technique that can reveal Jupiter-like planets in wide orbits, can enable detailed characterization of planetary atmospheres, and is a key step towards imaging Earth-like planets. Imaging detections are challenging due to the combined effect of small angular separation and large luminosity contrast between a planet and its host star. High-contrast observations with the Keck and Gemini telescopes have revealed three planets orbiting the star HR 8799, with projected separations of 24, 38, and 68 astronomical units. Multi-epoch data show counter-clockwise orbital motion for all three imaged planets. The low luminosity of the companions and the estimated age of the system imply planetary masses between 5 and 13 times that of Jupiter. This system resembles a scaled-up version of the outer portion of our Solar System.

In this review, we summarize our present knowledge of the behaviour of the mass-radius relationship from solar-type stars down to terrestrial planets, across the regime of substellar objects, brown dwarfs and giant planets. Particular attention is paid to the identification of the main physical properties or mechanisms responsible for this behaviour. Indeed, understanding the mechanical structure of an object provides valuable information about its internal structure, composition and heat content as well as its formation history. Although the general description of these properties is reasonably well mastered, disagreement between theory and observation in certain cases points to some missing physics in our present modelling of at least some of these objects. The mass-radius relationship in the overlaping domain between giant planets and low-mass brown dwarfs is shown to represent a powerful diagnostic to distinguish between these two different populations and shows once again that the present IAU distinction between these two populations at a given mass has no valid foundation.

Aims. Emission lines in polars show complex profiles with multiple components that are typically ascribed to the accretion stream, threading region, accretion spot, and the irradiated secondary-star. In low-state polars the fractional contribution by the accretion stream, and the accretion spot is greatly reduced offering an opportunity to study the effect of the secondary-star irradiation or stellar activity. We observed VV Pup during an exceptional low-state to study and constrain the properties of the line-forming regions and to search for evidence of chromospheric activity and/or irradiation. Methods. We obtained phase-resolved optical spectra at the ESO VLT+FORS1 with the aim of analyzing the emission line profile and radial velocity as a function of the orbital period. We also tailored irradiated secondary-star models to compare the predicted and the observed emission lines and to establish the nature of the line-forming regions. Results. Our observations and data analysis, when combined with models of the irradiated secondary-star, show that, while the weak low ionization metal lines (FeI and MgI) may be consistent with irradiation processes, the dominant Balmer H emission lines, as well as NaI and HeI, cannot be reproduced by the irradiated secondary-star models. We favor the secondary-star chromospheric activity as the main forming region and cause of the observed H, NaI, and He emission lines, though a threading region very close to the L1 point cannot be excluded.

Using a large sample of optical spectra of late-type dwarfs, we identify a subset of late-M through L field dwarfs that, because of the presence of low-gravity features in their spectra, are believed to be unusually young. From a combined sample of 303 field L dwarfs, we find observationally that 7.6+/-1.6% are younger than 100 Myr. This percentage is in agreement with theoretical predictions once observing biases are taken into account. We find that these young L dwarfs tend to fall in the southern hemisphere (Dec < 0 deg) and may be previously unrecognized, low-mass members of nearby, young associations like Tucana-Horologium, TW Hydrae, beta Pictoris, and AB Doradus. We use a homogeneously observed sample of roughly one hundred and fifty 6300-10000 Angstrom spectra of L and T dwarfs taken with the Low-Resolution Imaging Spectrometer at the W. M. Keck Observatory to examine the strength of the 6708-A Li I line as a function of spectral type and further corroborate the trends noted by Kirkpatrick et al. (2000). We use our low-gravity spectra to investigate the strength of the Li I line as a function of age. The data weakly suggest that for early- to mid-L dwarfs the line strength reaches a maximum for a few 100 Myr, whereas for much older (few Gyr) and much younger (<100 Myr) L dwarfs the line is weaker or undetectable. We show that a weakening of lithium at lower gravities is predicted by model atmosphere calculations, an effect partially corroborated by existing observational data. Larger samples containing L dwarfs of well determined ages are needed to further test this empirically. If verified, this result would reinforce the caveat first cited in Kirkpatrick et al. (2006) that the lithium test should be used with caution when attempting to confirm the substellar nature of the youngest brown dwarfs.

TW Hydrae shows significant radial-velocity variations in the optical regime. They have been attributed to a 10 Jupiter Mass planet orbiting the star at 0.04 AU. In this work, we have tested whether the observed RV variations can be caused by stellar spots. We have also analyzed new optical and infrared data to confirm the signal of the planet companion. We fitted the RV variations of TW Hya using a cool spot model. Our model shows that a cold spot covering 7% of the stellar surface and located at a latitude of 54 deg can reproduce the reported RV variations. The model also predicts a bisector semi-amplitude variation <10 m/s, which is less than the errors of the RV measurements discussed in an earlier publication. The analysis of new optical RV data, with typical errors of 10 m/s, shows a larger RV amplitude that varies depending on the correlation mask used. A slight correlation between the RV variation and the bisector is also observed, although not at a very significant level. The infrared H-band RV curve is almost flat, showing a small variation (<35 m/s) that is not consistent with the optical orbit. All these results support the spot scenario rather than the presence of a hot Jupiter around TW Hya.

We present Spitzer Space Telescope time series photometry of the exoplanet system HD 189733 spanning two times of secondary eclipse, when the planet passes out of view behind the parent star. We estimate the relative eclipse depth in 5 distinct bands and find the planet-to-star flux ratio to be 0.256 +/- 0.014% (3.6 microns), 0.214 +/- 0.020% (4.5 microns), 0.310 +/- 0.034% (5.8 microns), 0.391 +/- 0.022% (8.0 microns), and 0.598 +/- 0.038% (24 microns). For consistency, we re-analyze a previously published time series to deduce a contrast ratio in an additional band, 0.519 +/- 0.020% (16 microns). Our data are strongly inconsistent with a Planck spectrum, and we clearly detect emission near 4 microns as predicted by published theoretical models in which this feature arises from a corresponding opacity window. Unlike recent results for the exoplanet HD 209458b, we find that the emergent spectrum from HD 189733b is best matched by models that do not include an atmospheric temperature inversion. Taken together, these two studies provide initial observational support for the idea that hot Jupiter atmospheres diverge into two classes, in which a thermal inversion layer is present for the more strongly irradiated objects.

We present the discovery of two nearby L dwarfs from our 2MASS proper motion search, which uses multi-epoch 2MASS observations covering ~4700 square degrees of sky. 2MASS J18212815+1414010 and 2MASS J21481628+4003593 were overlooked by earlier surveys due to their faint optical magnitudes and their proximity to the Galactic Plane (10 degrees < |b| < 15 degrees). Assuming that both dwarfs are single, we derive spectrophotometric distances of ~10 pc, thus increasing the number of known L dwarfs within 10 pc to 10. In the near-infrared, 2MASS J21481628+4003593 shows a triangular-shaped H-band spectrum, strong CO absorption, and a markedly red J-Ks color (2.38+/-0.06) for its L6 optical spectral type. 2MASS J18212815+1414010 also shows a triangular-shaped H-band spectrum and a slightly red J-Ks color (1.78+/-0.05) for its L4.5 optical spectral type. Both objects show strong silicate absorption at 9-11 microns. Cumulatively, these features imply an unusually dusty photosphere for both of these objects. We examine several scenarios to explain the underlying cause for their enhanced dust content and find that a metal-rich atmosphere or a low-surface gravity are consistent with these results. 2MASS J18212815+1414010 may be young (and therefore have a low-surface gravity) based on its low tangential velocity of 10 km/s. On the other hand, 2MASS J21481628+4003593 has a high tangential velocity of 62 km/s and is therefore likely old. Hence, high metallicity and low-surface gravity may lead to similar effects.

We have carried out a search for the 2.14 micron spectroscopic signature of the close orbiting extrasolar giant planet, HD 179949b. High cadence time series spectra were obtained with the CRIRES spectrograph at VLT1 on two closely separated nights. Deconvolution yielded spectroscopic profiles with mean S/N ratios of several thousand, enabling the near infrared contrast ratios predicted for the HD 179949 system to be achieved. Recent models have predicted that the hottest planets may exhibit spectral signatures in emission due to the presence of TiO and VO which may be responsible for a temperature inversion high in the atmosphere. We have used our phase dependent orbital model and tomographic techniques to search for the planetary signature under the assumption of an absorption line dominated atmospheric spectrum, where T and V are depleted from the atmospheric model, and an emission line dominated spectrum, where TiO and VO are present. We do not detect a planet in either case, but the 2.120 - 2.174 micron wavelength region covered by our observations enables the deepest near infrared limits yet to be placed on the planet/star contrast ratio of any close orbiting extrasolar giant planet system. We are able to rule out the presence of an atmosphere dominated by absorption opacities in the case of HD 179949b at a contrast ratio of F_p/F_* ~ 1/3350, with 99 per cent confidence.

We examine the uncertainties in current planetary models and we quantify their impact on the planet cooling histories and mass-radius relationships. These uncertainties include (i) the differences between the various equations of state used to characterize the heavy material thermodynamical properties, (ii) the distribution of heavy elements within planetary interiors, (iii) their chemical composition and (iv) their thermal contribution to the planet evolution. Our models, which include a gaseous H/He envelope, are compared with models of solid, gasless Earth-like planets in order to examine the impact of a gaseous envelope on the cooling and the resulting radius. We find that for a fraction of heavy material larger than 20% of the planet mass, the distribution of the heavy elements in the planet's interior affects substantially the evolution and thus the radius at a given age. For planets with large core mass fractions ($\simgr$ 50%), such as the Neptune-mass transiting planet GJ436b, the contribution of the gravitational and thermal energy from the core to the planet cooling history is not negligible, yielding a $\sim$ 10% effect on the radius after 1 Gyr. We show that the present mass and radius determinations of the massive planet Hat-P-2b require at least 200 $\mearth$ of heavy material in the interior, at the edge of what is currently predicted by the core-accretion model for planet formation. We show that if planets as massive as $\sim$ 25 $\mjup$ can form, as predicted by improved core-accretion models, deuterium is able to burn in the H/He layers above the core, even for core masses as large as $\sim$ 100 $\mearth$. We provide extensive grids of planetary evolution models from 10 $\mearth$ to 10 M$_{\rm Jup}$, with various fractions of heavy elements.

We have carried out a multiwavelength observational campaign demonstrating some of the remarkable properties of the infrared-bright variable star BP Psc. Surrounded by a compact dusty, gaseous disk, this little-studied late-G (or early-K) type star emits about 75% of its detected energy flux at infrared wavelengths. Evidence for accretion of gas in conjunction with narrow bi-polar jets and Herbig-Haro objects is apparently consistent with classification of BP Psc as a pre-main sequence star, as postulated in most previous studies. If young, then BP Psc would be one of the nearest and oldest known classical T Tauri stars. However, such an evolutionary classification encounters various problems that are absent or much less severe if BP Psc is instead a luminosity class III post-main sequence star. In this case, it would be the first known example of a first ascent giant surrounded by a massive molecular disk with accompanying rapid gas accretion and prominent jets and HH objects. In this model, the genesis of the massive dusty gaseous disk could be a consequence of the envelopment of a low mass companion star. Properties in the disk may be conducive to the current formation of planets, a gigayear or more after the formation of BP Psc itself.

Using the Keck Interferometer, we spatially resolved the orbit of the pre-main sequence binary, Haro 1-14c, for the first time. We present these interferometric observations along with additional spectroscopic radial velocity measurements of the components. We performed a simultaneous orbit fit to the interferometric visibilities and the radial velocities of Haro 1-14c. Based on a statistical analysis of the possible orbital solutions that fit the data, we determined component masses of M_1 = 0.96 (+0.27/-0.08) Msun and M_2 = 0.33 (+0.09/-0.02) Msun for the primary and secondary, respectively, and a distance to the system of 111 (+19/-18) pc. The distance measurement is consistent with the close distance estimates of the Ophiuchus molecular cloud. Comparing our results with evolutionary tracks suggests an age of 3-4 Myr for Haro 1-14c. With additional interferometric measurements to improve the uncertainties in the masses and distance, we expect the low-mass secondary to provide important empirical data for calibrating the theoretical evolutionary tracks for pre-main sequence stars.

Detailed models are compared to recent infrared observations of the nearby extrasolar planet, HD 189733b. It is demonstrated that atmospheric water is present and that the planet's day side has a non-isothermal structure down to gas pressures of ~ 0.1 bars. Furthermore, model spectra with different amounts of CO are compared to the observations and an atmosphere absent of CO is excluded at roughly 2-sigma. Constraining the CO concentration beyond that is unfortunately not possible with the current Spitzer photometry. However, radically enhanced (or depleted) metal abundances are unlikely and the basic composition of this planet is probably similar to that of its host star. When combined with Spitzer observations, a recent ground-based upper limit for the K-band day side flux allows one to estimate the day-to-night energy redistribution efficiency to be ~ 43%.

We present Spitzer Space Telescope infrared photometry of a secondary eclipse of the hot Neptune GJ436b. The observations were obtained using the 8-micron band of the InfraRed Array Camera (IRAC). The data spanning the predicted time of secondary eclipse show a clear flux decrement with the expected shape and duration. The observed eclipse depth of 0.58 mmag allows us to estimate a blackbody brightness temperature of T_p = 717 +- 35 K at 8 microns. We compare this infrared flux measurement to a model of the planetary thermal emission, and show that this model reproduces properly the observed flux decrement. The timing of the secondary eclipse confirms the non-zero orbital eccentricity of the planet, while also increasing its precision (e = 0.14 +- 0.01). Additional new spectroscopic and photometric observations allow us to estimate the rotational period of the star and to assess the potential presence of another planet.

We obtained 238 spectra of the close orbiting extrasolar giant planet HD 189733b with resolution R ~ 15,000 during one night of observations with the near infrared spectrograph, NIRSPEC, at the Keck II Telescope. We have searched for planetary absorption signatures in the 2.0 - 2.4 micron region where H_2O and CO are expected to be the dominant atmospheric opacities. We employ a phase dependent orbital model and tomographic techniques to search for the planetary absorption signatures in the combined stellar and planetary spectra. Because potential absorption signatures are hidden in the noise of each single exposure, we use a model list of lines to apply a spectral deconvolution. The resulting mean profile possesses a S/N ratio that is 20 times greater than that found in individual lines. Our spectral timeseries thus yields spectral signatures with a mean S/N = 2720. We are unable to detect a planetary signature at a contrast ratio of log_10(F_p/F_*) = -3.40, with 63.8 per cent confidence. Our findings are not consistent with model predictions which nevertheless give a good fit to mid-infrared observations of HD 189733. The 1-sigma result is a factor of 1.7 times less than the predicted 2.185 micron planet/star flux ratio of log_10(F_p/F_*) ~ -3.16.

We present Spitzer Space Telescope infrared photometry of a primary transit of the hot Neptune GJ 436b. The observations were obtained using the 8 microns band of the InfraRed Array Camera (IRAC). The high accuracy of the transit data and the weak limb-darkening in the 8 microns IRAC band allow us to derive (assuming M = 0.44 +- 0.04 Msun for the primary) a precise value for the planetary radius (4.19 +0.21-0.16 Rearth), the stellar radius (0.463 +0.022-0.017 Rsun), the orbital inclination (85.90 +0.19-0.18 degrees) and transit timing (2454280.78186 +0.00015-0.00008 HJD). Assuming current planet models, an internal structure similar to that of Neptune with a small H/He envelope is necessary to account for the measured radius of GJ 436b.

We identify two classes of transiting planet, based on their equilibrium temperatures and Safronov numbers. We examine various possible explanations for the dichotomy. It may reflect the influence of planet or planetesimal scattering in determining when planetary migration stops. Another possibility is that some planets lose more mass to evaporation than others. If this evaporation process preferentially removes Helium from the planet, the consequent reduction in the mean molecular weight may explain why some planets have anomalously large radii.

Tables of mean opacities are often used to compute the transfer of radiation in a variety of astrophysical simulations from stellar evolution models to proto-planetary disks. Often tables, such as Ferguson et al. (2005), are computed with a predetermined set of physical assumptions that may or may not be valid for a specific application. This paper explores the effects of several assumptions of grain physics on the Rosseland mean opacity in an oxygen rich environment. We find that changing the distribution of grain sizes, either the power-law exponent or the shape of the distribution, has a marginal effect on the total mean opacity. We also explore the difference in the mean opacity between solid homogenous grains and grains that are porous or conglomorations of several species. Changing the amount of grain opacity included in the mean by assuming a grain-to-gas ratio significantly affects the mean opacity, but in a predictable way.

We present a search for the near infrared spectroscopic signature of the close orbiting extrasolar giant planet HD 75289b. We obtained ~230 spectra in the wavelength range 2.18 - 2.19 microns using the Phoenix spectrograph at Gemini South. By considering the direct spectrum, derived from irradiated model atmospheres, we search for the absorption profile signature present in the combined star and planet light. Since the planetary spectrum is separated from the stellar spectrum at most phases, we apply a phase dependent orbital model and tomographic techniques to search for absorption signatures. Because the absorption signature lies buried in the noise of a single exposure we apply a multiline deconvolution to the spectral lines available in order to boost the effective S/N ratio of the data. The wavelength coverage of 80 angstroms is expected to contain ~100 planetary lines, enabling a mean line with S/N ratio of ~800 to be achieved after deconvolution. We are nevertheless unable to detect the presence of the planet in the data and carry out further simulations to show that broader wavelength coverage should enable a planet like HD 75289b to be detected with 99.9 per cent (4 sigma) confidence. We investigate the sensitivity of our method and estimate detection tolerances for mismatches between observed and model planetary atmospheres.

Water absorption is identified in the atmosphere of HD209458b by comparing models for the planet's transmitted spectrum to recent, multi-wavelength, eclipse-depth measurements (from 0.3 to 1 microns) published by Knutson et al. (2007). A cloud-free model which includes solar abundances, rainout of condensates, and photoionization of sodium and potassium is in good agreement with the entire set of eclipse-depth measurements from the ultraviolet to near-infrared. Constraints are placed on condensate removal by gravitational settling, the bulk metallicity, and the redistribution of absorbed stellar flux. Comparisons are also made to the Charbonneau et al. (2002) sodium measurements.

We imaged five objects near the star forming clouds of Ophiuchus with the Keck Laser Guide Star AO system. We resolved Allers et al. (2006)'s #11 (Oph 16222-2405) and #16 (Oph 16233-2402) into binary systems. The #11 object is resolved into a 243 AU binary, the widest known for a very low mass (VLM) binary. The binary nature of #11 was discovered first by Allers (2005) and independently here during which we obtained the first spatially resolved R~2000 near-infrared (J & K) spectra, mid-IR photometry, and orbital motion estimates. We estimate for 11A and 11B gravities (log(g)>3.75), ages (5+/-2 Myr), luminosities (log(L/Lsun)=-2.77+/-0.10 and -2.96+/-0.10), and temperatures (Teff=2375+/-175 and 2175+/-175 K). We find self-consistent DUSTY evolutionary model (Chabrier et al. 2000) masses of 17+4-5 MJup and 14+6-5 MJup, for 11A and 11B respectively. Our masses are higher than those previously reported (13-15 MJup and 7-8 MJup) by Jayawardhana & Ivanov (2006b). Hence, we find the system is unlikely a ``planetary mass binary'', (in agreement with Luhman et al. 2007) but it has the second lowest mass and lowest binding energy of any known binary. Oph #11 and Oph #16 belong to a newly recognized population of wide (>100 AU), young (<10 Myr), roughly equal mass, VLM stellar and brown dwarf binaries. We deduce that ~6+/-3% of young (<10 Myr) VLM objects are in such wide systems. However, only 0.3+/-0.1% of old field VLM objects are found in such wide systems. Thus, young, wide, VLM binary populations may be evaporating, due to stellar encounters in their natal clusters, leading to a field population depleted in wide VLM systems.

We present the discovery of 2MASS J21321145+1341584AB as a closely separated (0.066") very low-mass field dwarf binary resolved in the near-infrared by the Keck II Telescope using laser guide star adaptive optics. Physical association is deduced from the angular proximity of the components and constraints on their common proper motion. We have obtained a near-infrared spectrum of the binary and find that it is best described by an L5+/-0.5 primary and an L7.5+/-0.5 secondary. Model-dependent masses predict that the two components straddle the hydrogen burning limit threshold with the primary likely stellar and the secondary likely substellar. The properties of this sytem - close projected separation (1.8+/-0.3 AU) and near unity mass ratio - are consistent with previous results for very low-mass field binaries. The relatively short estimated orbital period of this system (~7-12 yr) makes it a good target for dynamical mass measurements. Interestingly, the system's angular separation is the tightest yet for any very low-mass binary published from a ground-based telescope and is the tightest binary discovered with laser guide star adaptive optics to date.

We have re-analyzed archival HST and Subaru data of the recently discovered planetary mass companion (PMC) GQ Lup B. With these we produce the first R and I band photometry of the companion and fit a radius and effective temperature using detailed model atmospheres. We find an effective temperature of 2335 +/- 100K, and a radius of 0.38 +/- 0.05 Rsol and luminosity of log(L/Lsol) = -2.42 +/- 0.07 (at 140pc). Since we fit wavelengths that span most of the emitted radiation from GQ Lup this luminosity estimate is robust, with uncertainty dominated by the distance uncertainty (+/- 50 pc). The radius obtained for 140pc (0.38Rsol) is significantly larger than the one originally derived and larger than model predictions. The mass of the object is much more model-dependent than the radiative properties, but for the GAIA dusty models we find a mass between 10-20 MJup, in the range of the brown dwarf and PMC deuterium burning boundary. Assuming a distance of 140pc, observations fit to 1sigma the Baraffe evolution model for a ~15 MJup brown dwarf. Additionally, the F606W photometric band is significantly overluminous compared to model predictions and other brown dwarfs. Such overluminosity could be explained by a bright Halpha emission from chromospheric activity, interaction with another undetected companion, or accretion. Assuming that GQ Lup B has a bright Halpha emission line, its Halpha emission strength is 10^(-1.71 +/- 0.10) Lbol, significantly larger than field late-type dwarfs. GQ Lup B might be strongly accreting and still be in its formation phase.

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 review recent theoretical progress aimed at understanding the formation and the early stages of evolution of giant planets, low-mass stars and brown dwarfs. Calculations coupling giant planet formation, within a modern version of the core accretion model, and subsequent evolution yield consistent determinations of the planet structure and evolution. Because of the uncertainties in the initial conditions, however, it is not possible to say whether young planets are faint or bright compared with low-mass young brown dwarfs. We review the effects of irradiation and evaporation on the evolution of short period planets and argue that substantial mass loss may have occurred for these objects. Concerning star formation, geometrical effects in protostar core collapse are examined by comparing 1D and 3D calculations. Spherical collapse is shown to overestimate the core inner density and temperature and thus to yield incorrect initial conditions for PMS or young brown dwarf evolution. Accretion is also shown to occur over a very limited fraction of the protostar surface. Accretion affects the evolution of young brown dwarfs and yields more compact structures for a given mass and age, thus fainter luminosities. This can lead to severe misinterpretations of the mass and/or age of young accreting objects from their location in the HR diagram. We argue that newborn stars and brown dwarfs should appear rapidly over an extended area in the HR diagram, depending on their accretion history, rather than on a well defined birth line. Finally, we suggest that the distinction between planets and brown dwarfs be based on an observational diagnostic, reflecting the different formation mechanisms between these two distinct populations, rather than on an arbitrary, confusing definition.

While following up L dwarf candidates selected photometrically from the Two Micron All Sky Survey, we uncovered an unusual object designated 2MASS J01415823-4633574. Its optical spectrum exhibits very strong bands of vanadium oxide but abnormally weak absorptions by titanium oxide, potassium, and sodium. Morphologically such spectroscopic characteristics fall intermediate between old, field early-L dwarfs (log(g)~5) and very late M giants (log(g)~0), leading us to favor low gravity as the explanation for the unique spectral signatures of this L dwarf. Such a low gravity can be explained only if this L dwarf is much lower in mass than a typical old field L dwarf of similar temperature and is still contracting to its final radius. These conditions imply a very young age. Further evidence of youth is found in the near-infrared spectrum, including a triangular-shaped H-band continuum reminiscent of young brown dwarf candidates discovered in the Orion Nebula Cluster. Using the above information along with comparisons to brown dwarf atmospheric and interior models, our current best estimate is that this L dwarf has an age of 1-50 Myr and a mass of 6-25 M_Jupiter. The location of 2MASS 0141-4633 on the sky coupled with a distance estimate of ~35 pc and the above age estimate suggests that this object may be a brown dwarf member of either the 30-Myr-old Tucana/Horologium Association or the ~12-Myr-old beta Pic Moving Group.

Recently the Spitzer Space Telescope observed the transiting extrasolar planets, TrES-1 and HD209458b. These observations have provided the first estimates of the day side thermal flux from two extrasolar planets orbiting Sun-like stars. In this paper, synthetic spectra from atmospheric models are compared to these observations. The day-night temperature difference is explored and phase-dependent flux densities are predicted for both planets. For HD209458b and TrES-1, models with significant day-to-night energy redistribution are required to reproduce the observations. However, the observational error bars are large and a range of models remains viable.

We compare evolutionary models for close-in exoplanets coupling irradiation and evaporation due respectively to the thermal and high energy flux of the parent star with observations of recently discovered new transiting planets. The models provide an overall good agreement with observations, although at the very limit of the quoted error bars of OGLE-TR-10, depending on its age. Using the same general theory, we show that the three recently detected hot-Neptune planets (GJ436, $\rho$ Cancri, $\mu$ Ara) may originate from more massive gas giants which have undergone significant evaporation. We thus suggest that hot-Neptunes and hot-Jupiters may share the same origin and evolution history. Our scenario provides testable predictions in terms of the mass-radius relationships of these hot-Neptunes.

AB Dor C, a low-mass companion to the young star AB Dor A has recently been imaged. Direct detection and astrometric observations have been used to solve the AB Dor C orbit and obtain a dynamical mass estimate. Near-infrared magnitudes and effective temperatures from models have then been compared with observed values. Models and observations disagree at more than the 2 sigma level. This result prompts Close et al. 2005 to claim that models are overestimating young low-mass object luminosities by roughly a factor of two. This claim is based on the hypothesis that the detected source, AB Dor C, is a single object. Another possible interpretation of the data that was not considered is that AB Dor C is an unresolved binary brown dwarf. Considering that 21 brown dwarfs companions have been found to date around stars and that three of them are confirmed or potential brown dwarf binaries (Gl569Ba & Bb, GL564B & C and Indi Ba & Bb), this suggests that brown dwarf binaries may constitute a significant fraction (14% +/- 8%) of the brown dwarf population around stars. For these reasons, the possibility that AB Dor C is an unresolved binary brown dwarf deserves attention.

Previous computations of low temperature Rosseland and Planck mean opacities from Alexander & Ferguson (1994) are updated and expanded. The new computations include a more complete equation of state with more grain species and updated optical constants. Grains are now explicitly included in thermal equilibrium in the equation of state calculation, which allows for a much wider range of grain compositions to be accurately included than was previously the case. The inclusion of high temperature condensates such as Al$_2$O$_3$ and CaTiO$_3$ significantly affects the total opacity over a narrow range of temperatures before the appearance of the first silicate grains. The new opacity tables are tabulated for temperatures ranging from 30000 K to 500 K with gas densities from 10$^{-4}$ g cm$^{-3}$ to 10$^{-19}$ g cm$^{-3}$. Comparisons with previous Rosseland mean opacity calculations are discussed. At high temperatures, the agreement with OPAL and Opacity Project is quite good. Comparisons at lower temperatures are more divergent as a result of differences in molecular and grain physics included in different calculations. The computation of Planck mean opacities performed with the opacity sampling method are shown to require a very large number of opacity sampling wavelength points; previously published results obtained with fewer wavelength points are shown to be significantly in error. Methods for requesting or obtaining the new tables are provided.

We have obtained time series optical spectra of the cataclysmic variable PQ And in quiescence. The spectra show a white dwarf continuum with narrow Balmer emission superimposed over strong Balmer absorption. The emission lines have blue and red components whose strength changes with time. An analysis of the H$\alpha$ emission line implies a short orbital period below the period gap. Given its lack of accretion disk features, its large and infrequent outbursts, and an orbital period below the period gap, PQ And is probably a low accretion rate object similar to WZ Sge. In addition, white dwarf model fits imply that PQ And is an excellent ZZ Cet candidate.

Model atmospheres have been computed for M dwarfs that are strongly irradiated by nearby hot companions. A variety of primary and secondary spectral types are explored in addition to models specific to four known systems: GD 245, NN Ser, AA Dor, and UU Sge. This work demonstrates that a dramatic temperature inversion is possible on at least one hemisphere of an irradiated M dwarf and the emergent spectrum will be significantly different from an isolated M dwarf or a black body flux distribution. For the first time, synthetic spectra suitable for direct comparison to high-resolution observations of irradiated M dwarfs in non-mass transferring post-common envelope binaries are presented. The effects of departures from local thermodynamic equilibrium on the Balmer line profiles are also discussed.

Brown dwarfs bridge the gap between the stellar and planetary mass regimes. Evolving from conditions very similar to the lowest-mass stars, the atmospheres of older brown dwarfs closely resemble those expected in close-in extrasolar giant planets, and with cooler BDs still being discovered, more and more approach the properties of gas giants at wider separation. Interpreting the spectra of BDs is therefore a crucial step towards understanding and predicting the spectral and thermal properties of EGPs. Essential properties of substellar atmospheres are massive molecular line-blanketing and the condensation of species with decreasing Teff, changing the chemical equilibrium composition and causing absorption from dust grains. More complex details involve the distribution of dust clouds over the surface giving rise to temporal variability, and possible deviations from chemical equilibrium conditions. In the case of close-in EGPs and some BDs in binary systems, the effect of irradiation from the primary significantly affects the spectral properties and thermal evolution.

We include the effect of evaporation in our evolutionary calculations of close-in giant planets, based on a recent model for thermal evaporation taking into account the XUV flux of the parent star (Lammer et al. 2003). Our analysis leads to the existence of a critical mass for a given orbital distance $m_{\rm crit}(a)$ below which the evaporation timescale becomes shorter than the thermal timescale of the planet. For planets with initial masses below $m_{\rm crit}$, evaporation leads to a rapid expansion of the outer layers and of the total planetary radius, speeding up the evaporation process. Consequently, the planet does not survive as long as estimated by a simple application of mass loss rates without following consistently its evolution. We find out that the transit planet HD 209458b might be in such a dramatic phase, although with an extremely small probability. As a consequence, we predict that, after a certain time, only planets above a value $m_{\rm crit}(a)$ should be present at an orbital distance $a$ of a star. For planets with initial masses above $m_{\rm crit}$, evaporation does not affect the evolution of the radius with time.

We review our present understanding of the physical properties of substellar objects, brown dwarfs and irradiated or non-irradiated gaseous exoplanets. This includes a description of their internal properties, mechanical structure and heat content, their atmospheric properties, thermal profile and emergent spectrum, and their evolution, in particular as irradiated companions of a close parent star. The general theory can be used to make predictions in term of detectability for the future observational projects. Special attention is devoted to the evolution of the two presently detected transit planets, HD209458B and OGLE-TR-56B. For this latter, we present a consistent evolution for its recently revised mass and show that we reproduce the observed radius within its error bars. We briefly discuss differences between brown dwarfs and gaseous planets, both in terms of mass function and formation process. We outline several arguments to show that the minimum mass for deuterium burning, recently adopted officially as the limit to distinguish the two types of objects, is unlikely to play any specific role in star formation, so that such a limit is of purely semantic nature and is not supported by a physical justification.

Extending the theory we derived recently for HD209458b to different cases of strongly irradiated gaseous exoplanets, we have calculated the consistent evolution of the new transiting planet, OGLE-TR-56b, for its recently revised mass determination. The theory is shown to successfully reproduce the observed radius, for the proper age of the system. We also examine the dissipation of kinetic energy at the planet's internal adiabat due to atmospheric winds, and place constraints on the efficiency of this process. We show that a fraction $\sim 0.1$ -- 0.5% of the incident flux transformed into thermal energy deposited at the adiabatic level can accommodate the observed radii of both OGLE-TR-56b and HD209458b. The present theory yields quantitative predictions on the evolution of the emergent spectrum and fundamental properties of hot-jupiters. The predictions for radius, luminosity, temperature as a function of the planet's mass and orbital distance can be used as benchmarks for future detections of transit planets.

We present evolutionary models for cool brown dwarfs and extra-solar giant planets. The models reproduce the main trends of observed methane dwarfs in near-IR color-magnitude diagrams. We also present evolutionary models for irradiated planets, coupling for the first time irradiated atmosphere profiles and inner structures. We focus on HD 209458-like systems and show that irradiation effects can substantially affect the radius of sub-jovian mass giant planets. Irradiation effects, however, cannot alone explain the large observed radius of HD 209458b. Adopting assumptions which optimise irradiation effects and taking into account the extension of the outer atmospheric layers, we still find $\sim$ 20% discrepancy between observed and theoretical radii. An extra source of energy seems to be required to explain the observed value of the first transit planet.

We briefly describe the current version of the PHOENIX code. We then present some illustrative results from the modeling of Type Ia and Type II supernovae, hot stars, and irradiated giant planets. Good fits to observations can be obtained, when account is taken for spherically symmetric, line-blanketed, static or expanding atmospheres.

The recent announcement that sodium absorption has been observed in the atmosphere of HD209458b, the only EGP observed to transit its parent star, is the first direct detection of an EGP atmosphere. We explore the possibility that neutral sodium is {\em not} in local thermodynamic equilibrium (LTE) in the outer atmosphere of irradiated EGPs and that the sodium concentration may be underestimated by models that make the LTE assumption. Our results indicate that it may not be necessary to invoke excessive photoionization, low metallicity, or even high altitude clouds to explain the observations.

We have modeled irradiated planets located near a dM5 and a G2 primary star. The impinging radiation field was explicitly included in the solution of the radiative transfer equation and in the computation of the atmospheric structure. We find that large errors in both the thermal and reflected flux will result from models which do not include the impinging radiation in a self-consistent manner. A cool (teff = 500K) and a hot (teff = 1000K) planet were modeled at various orbital separations from both the dM5 and the G2 primary. In all scenarios, we compared the effects of the irradiation in two extreme cases: one where dust clouds form and remain suspended in the atmosphere, and another where dust clouds form but completely settle out of the atmosphere. The atmospheric structure and emergent spectrum strongly depend on the presence or absence of dust clouds. We find that, in the absence of dust opacity, the impinging radiation significantly alters the innermost layers of an EGP atmosphere and that they are actually brighter in the optical than dusty planets. Our models also indicate that the planet-to-star brightness ratio in the optical will be less that 1 x 10^-5 for objects like tau bootis which is consistent with recently reported upper limit values.