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.

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

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

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.

HISPEC (High-resolution Infrared Spectrograph for Exoplanet Characterization) is a proposed diffraction-limited spectrograph for the W.M. Keck Observatory, and a pathfinder for the MODHIS facility project (Multi-Object Diffraction-limited High-resolution Infrared Spectrograph) on the Thirty Meter Telescope. HISPEC/MODHIS builds on diffraction-limited spectrograph designs which rely on adaptively corrected single-mode fiber feeds. Seeing-limited high-resolution spectrographs, by virtue of the conservation of beam etendue, grow in volume following a D^3 power law (D is the telescope diameter), and are subject to daunting challenges associated with their large size. Diffraction-limited spectrographs fed by single mode fibers are decoupled from the telescope input, and are orders of magnitude more compact and have intrinsically stable line spread functions. Their efficiency is directly proportional to the performance of the adaptive optics (AO) system. AO technologies have matured rapidly over the past two decades and are baselined for future extremely large telescopes. HISPEC/MODHIS will take R>100,000 spectra of a few objects in a 10" field-of-view sampled at the diffraction limit (~10-50 mas), simultaneously from 0.95 to 2.4 microns (y-K). The scientific scope ranges from exoplanet infrared precision radial velocities, spectroscopy of transiting, close-in, and directly imaged exoplanets (atmospheric composition and dynamics, RM effect, spin measurements, Doppler imaging), brown dwarf characterization, stellar physics/chemistry, proto-planetary disk kinematics/composition, Solar system, extragalactic science, and cosmology. HISPEC/MODHIS features a compact, cost-effective design optimized to fully exploit the existing Keck-AO and future TMT-NFIRAOS infrastructures and boost the scientific reach of Keck Observatory and TMT soon after first light.

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.

We propose a new approach to Bayesian prior probability distributions (priors) that can improve orbital solutions for low-phase-coverage orbits, where data cover less than approximately 40% of an orbit. In instances of low phase coverage such as with stellar orbits in the Galactic center or with directly-imaged exoplanets, data have low constraining power and thus priors can bias parameter estimates and produce under-estimated confidence intervals. Uniform priors, which are commonly assumed in orbit fitting, are notorious for this. We propose a new observable-based prior paradigm that is based on uniformity in observables. We compare performance of this observable-based prior and of commonly assumed uniform priors using Galactic center and directly-imaged exoplanet (HR 8799) data. The observable-based prior can reduce biases in model parameters by a factor of two and helps avoid under-estimation of confidence intervals for simulations with less than about 40% phase coverage. Above this threshold, orbital solutions for objects with sufficient phase coverage such as S0-2, a short-period star at the Galactic center with full phase coverage, are consistent with previously published results. Below this threshold, the observable-based prior limits prior influence in regions of prior dominance and increases data influence. Using the observable-based prior, HR 8799 orbital analyses favor lower eccentricity orbits and provide stronger evidence that the four planets have a consistent inclination around 30 degrees to within 1-sigma. This analysis also allows for the possibility of coplanarity. We present metrics to quantify improvements in orbital estimates with different priors so that observable-based prior frameworks can be tested and implemented for other low-phase-coverage orbits.

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

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

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

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

As we begin to discover rocky planets in the habitable zone of nearby stars with missions like TESS and CHEOPS, we will need quick advancements on instrumentation and observational techniques that will enable us to answer key science questions, such as What are the atmospheric characteristics of habitable zone rocky planets? How common are Earth-like biosignatures in rocky planets?} How similar or dissimilar are those planets to Earth? Over the next decade we expect to have discovered several Earth-analog candidates, but we will not have the tools to study the atmospheres of all of them in detail. Ground-based ELTs can identify biosignatures in the spectra of these candidate exo-Earths and understand how the planets' atmospheres compare to the Earth at different epochs. Transit spectroscopy, high-resolution spectroscopy, and reflected-light direct imaging on ELTs can identify multiple biosignatures for habitable zone, rocky planets around M stars at optical to near-infrared wavelengths. Thermal infrared direct imaging can detect habitable zone, rocky planets around AFGK stars, identifying ozone and motivating reflected-light follow-up observations with NASA missions like HabEx/LUVOIR. Therefore, we recommend that the Astro2020 Decadal Survey Committee support: (1) the search for Earth-like biosignatures on rocky planets around nearby stars as a key science case; (2) the construction over the next decade of ground-based Extremely Large Telecopes (ELTs), which will provide the large aperture and spatial resolution necessary to start revealing the atmospheres of Earth-analogues around nearby stars; (3) the development of instrumentation that optimizes the detection of biosignatures; and (4) the generation of accurate line lists for potential biosignature gases, which are needed as model templates to detect those molecules.

Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe.

Substellar multiplicity is a key outcome of the formation process. The biggest challenge for the next decade will be to distinguish between the formation history, environmental conditions, and dynamical evolution leading to the least massive brown dwarfs and the most massive planets at the tail ends of their mass functions. In this white paper, we advocate for a comprehensive characterization of both the statistical distributions of the population of ultracool dwarf multiple systems and the fundamental properties of their individual components as a function of age. A space-based precision astrometry mission in near-infrared wavelengths would provide the necessary measurements to identify and characterize age-calibrated populations of multiple systems.

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

We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties and compositions, and how they interact with protoplanetary disks. While we have begun to detect protoplanet candidates during the last several years, we are presently only sensitive to the widest separation, highest mass / accretion rate cases. Current observing facilities lack the angular resolution and inner working angle to probe the few-AU orbital separations where giant planet formation is thought to be most efficient. They also lack the contrast to detect accretion rates that would form lower mass gas giants and ice giants. Instruments and telescopes coming online over the next decade will provide high contrast in the inner giant-planet-forming regions around young stars, allowing us to build a protoplanet census and to characterize planet formation in detail for the first time.

Having discovered that Earth-sized planets are common, we are now embarking on a journey to determine if Earth-like planets are also common. Finding Earth-like planets is one of the most compelling endeavors of the 21st century - leading us toward finally answering the question: Are we alone? To achieve this forward-looking goal, we must determine the masses of the planets; the sizes of the planets, by themselves, are not sufficient for the determination of the bulk and atmospheric compositions. Masses, coupled with the radii, are crucial constraints on the bulk composition and interior structure of the planets and the composition of their atmospheres, including the search for biosignatures. Precision radial velocity is the most viable technique for providing essential mass and orbit information for spectroscopy of other Earths. The development of high quality precision radial velocity instruments coupled to the building of the large telescope facilities like TMT and GMT or space-based platforms like EarthFinder can enable very high spectral resolution observations with extremely precise radial velocities on minute timescales to allow for the modeling and removal of radial velocity jitter. Over the next decade, the legacy of exoplanet astrophysics can be cemented firmly as part of humankind's quest in finding the next Earth - but only if we can measure the masses and orbits of Earth-sized planets in habitable zone orbits around Sun-like stars.

High resolution spectroscopy of the lowest-mass stars and brown dwarfs reveals their origins, multiplicity, compositions and physical properties, with implications for the star formation and chemical evolution history of the Milky Way. We motivate the need for high-resolution, infrared spectroscopic surveys to reach these faint sources.

The Orion Nebula Cluster (ONC) is the nearest site of ongoing massive star formation, which allows us to study the kinematics and dynamics of the region in detail and constrain star formation theories. Using HST ACS/WFPC2/WFC3IR and Keck II NIRC2 data, we have measured the proper motions of 701 stars within an $\sim6'\times6'$ field of view around the center of the ONC. We have found more than 10 escaping star candidates, concentrated predominantly at the core of the cluster. The proper motions of the bound stars are consistent with a normal distribution, albeit elongated North-South along the Orion filament, with proper motion dispersions of $(\sigma_{\mu,\alpha^*}, \sigma_{\mu,\delta}) = (0.83\pm0.02,\,1.12\pm0.03)$ mas yr$^{-1}$ or intrinsic velocity dispersions of $(\sigma_{v,\alpha^*}, \sigma_{v,\delta}) = (1.57\pm0.04,\,2.12\pm0.06)$ km s$^{-1}$ assuming a distance of 400 pc to the ONC. The cluster shows no evidence for tangential-to-radial anisotropy. Our velocity dispersion profile agrees with the prediction from the observed stellar + gas density profile from Da Rio et al. (2014), indicating that the ONC is in virial equilibrium. This finding suggests that the cluster was formed with a low star formation efficiency per dynamical timescale based on comparisons with current star formation theories. Our survey also recovered high-velocity IR sources BN, x, and n in the BN/KL region. The estimated location of the first two sources $\sim500$ years ago agrees with that of the radio source I, consistent with their proposed common origin from a multi-stellar disintegration. However, source n appears to have a small proper motion and is unlikely to have been involved in the event.

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.

Direct detection, also known as direct imaging, is a method for discovering and characterizing the atmospheres of planets at intermediate and wide separations. It is the only means of obtaining spectra of non-transiting exoplanets. Characterizing the atmospheres of planets in the <5 AU regime, where RV surveys have revealed an abundance of other worlds, requires a 30-m-class aperture in combination with an advanced adaptive optics system, coronagraph, and suite of spectrometers and imagers - this concept underlies planned instruments for both TMT (the Planetary Systems Imager, or PSI) and the GMT (GMagAO-X). These instruments could provide astrometry, photometry, and spectroscopy of an unprecedented sample of rocky planets, ice giants, and gas giants. For the first time habitable zone exoplanets will become accessible to direct imaging, and these instruments have the potential to detect and characterize the innermost regions of nearby M-dwarf planetary systems in reflected light. High-resolution spectroscopy will not only illuminate the physics and chemistry of exo-atmospheres, but may also probe rocky, temperate worlds for signs of life in the form of atmospheric biomarkers (combinations of water, oxygen and other molecular species). By completing the census of non-transiting worlds at a range of separations from their host stars, these instruments will provide the final pieces to the puzzle of planetary demographics. This whitepaper explores the science goals of direct imaging on 30-m telescopes and the technology development needed to achieve them.

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.

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.

We present adaptive-optics assisted near-infrared high-spectral resolution observations of late-type giants in the nuclear star cluster of the Milky Way. The metallicity and elemental abundance measurements of these stars offer us an opportunity to understand the formation and evolution of the nuclear star cluster. In addition, their proximity to the supermassive black hole ($\sim 0.5$ pc) offers a unique probe of the star formation and chemical enrichment in this extreme environment. We observed two stars identified by medium spectral-resolution observations as potentially having very high metallicities. We use spectral-template fitting with the PHOENIX grid and Bayesian inference to simultaneously constrain the overall metallicity, [M/H], alpha-element abundance [$\alpha$/Fe], effective temperature, and surface gravity of these stars. We find that one of the stars has very high metallicity ([M/H] $> 0.6$) and the other is slightly above solar metallicity. Both Galactic center stars have lines from scandium (Sc), vanadium (V), and yttrium (Y) that are much stronger than allowed by the PHOENIX grid. We find, using the spectral synthesis code Spectroscopy Made Easy, that [Sc/Fe] may be an order of magnitude above solar. For comparison, we also observed an empirical calibrator in NGC6791, the highest metallicity cluster known ([M/H] $\sim 0.4$). Most lines are well matched between the calibrator and the Galactic center stars, except for Sc, V, and Y, which confirms that their abundances must be anomalously high in these stars. These unusual abundances, which may be a unique signature of nuclear star clusters, offer an opportunity to test models of chemical enrichment in this region.

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.

We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf.

At a distance of 2~pc, our nearest brown dwarf neighbor, Luhman 16 AB, has been extensively studied since its discovery 3 years ago, yet its most fundamental parameter -- the masses of the individual dwarfs -- has not been constrained with precision. In this work we present the full astrometric orbit and barycentric motion of Luhman 16 AB and the first precision measurements of the individual component masses. We draw upon archival observations spanning 31 years from the European Southern Observatory (ESO) Schmidt Telescope, the Deep Near-Infrared Survey of the Southern Sky (DENIS), public FORS2 data on the Very Large Telescope (VLT), and new astrometry from the Gemini South Multiconjugate Adaptive Optics System (GeMS). Finally, we include three radial velocity measurements of the two components from VLT/CRIRES, spanning one year. With this new data sampling a full period of the orbit, we use a Markov Chain Monte Carlo algorithm to fit a 16-parameter model incorporating mutual orbit and barycentric motion parameters and constrain the individual masses to be~$27.9^{+1.1}_{-1.0}$~$M_{J}$ for the T dwarf and~$34.2^{+1.3}_{-1.1}$~$M_{J}$ for the L dwarf. Our measurements of Luhman 16 AB's mass ratio and barycentric motion parameters are consistent with previous estimates in the literature utilizing recent astrometry only. The GeMS-derived measurements of the Luhman 16 AB separation in 2014-2015 agree closely with Hubble Space Telescope (HST) measurements made during the same epoch Bedin et al. 2017, and the derived mutual orbit agrees with those measurements to within the HST uncertainties of $0.3 - 0.4$ milliarcseconds.

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 describe a Bayesian rejection sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of MCMC in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to ten long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.

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

The Gemini Planet Imager (GPI) is a next-generation high-contrast imager built for the Gemini Observatory. The GPI exoplanet survey (GPIES) consortium is made up of 102 researchers from 28 institutions in North and South America and Europe. In November 2014, we launched a search for young Jovian planets and debris disks. In this paper, we discuss how we have coordinated the work done by this large team to improve the technical and scientific productivity of the campaign, and describe lessons we have learned that could be useful for future instrumentation-based astronomical surveys. The success of GPIES lies mostly on its decentralized structure, clear definition of policies that are signed by each member, and the heavy use of modern tools for communicating, exchanging information, and processing data.

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.

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.

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

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.

We present astrometric monitoring of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager between 2013 and 2016. A small but significant position angle change is detected at constant separation; the orbital motion is confirmed with literature measurements. Efficient Monte Carlo techniques place preliminary constraints on the orbital parameters of HD 95086 b. With 68% confidence, a semimajor axis of 61.7^{+20.7}_{-8.4} au and an inclination of 153.0^{+9.7}_{-13.5} deg are favored, with eccentricity less than 0.21. Under the assumption of a co-planar planet-disk system, the periastron of HD 95086 b is beyond 51 au with 68% confidence. Therefore HD 95086 b cannot carve the entire gap inferred from the measured infrared excess in the SED of HD 95086. We use our sensitivity to additional planets to discuss specific scenarios presented in the literature to explain the geometry of the debris belts. We suggest that either two planets on moderately eccentric orbits or three to four planets with inhomogeneous masses and orbital properties are possible. The sensitivity to additional planetary companions within the observations presented in this study can be used to help further constrain future dynamical simulations of the planet-disk system.

We present new $H$ (1.5-1.8 $\mu$m) photometric and $K_1$ (1.9-2.2 $\mu$m) spectroscopic observations of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager. The $H$-band magnitude has been significantly improved relative to previous measurements, whereas the low resolution $K_1$ ($\lambda/\delta\lambda \approx 66$) spectrum is featureless within the measurement uncertainties, and presents a monotonically increasing pseudo-continuum consistent with a cloudy atmosphere. By combining these new measurements with literature $L^{\prime}$ photometry, we compare the spectral energy distribution of the planet to other young planetary-mass companions, field brown dwarfs, and to the predictions of grids of model atmospheres. HD 95086 b is over a magnitude redder in $K_1-L^{\prime}$ color than 2MASS J12073346-3932539 b and HR 8799 c and d, despite having a similar $L^{\prime}$ magnitude. Considering only the near-infrared measurements, HD 95086 b is most analogous to the brown dwarfs 2MASS J2244316+204343 and 2MASS J21481633+4003594, both of which are thought to have dusty atmospheres. Morphologically, the spectral energy distribution of HD 95086 b is best fit by low temperature ($T_{\rm eff} =$ 800-1300 K), low surface gravity spectra from models which simulate high photospheric dust content. This range of effective temperatures is consistent with field L/T transition objects, but the spectral type of HD 95086 b is poorly constrained between early L and late T due to its unusual position the color-magnitude diagram, demonstrating the difficulty in spectral typing young, low surface gravity substellar objects. As one of the reddest such objects, HD 95086 b represents an important empirical benchmark against which our current understanding of the atmospheric properties of young extrasolar planets can be tested.

We present the first scattered-light image of the debris disk around HD 131835 in $H$ band using the Gemini Planet Imager. HD 131835 is a $\sim$15 Myr old A2IV star at a distance of $\sim$120 pc in the Sco-Cen OB association. We detect the disk only in polarized light and place an upper limit on the peak total intensity. No point sources resembling exoplanets were identified. Compared to its mid-infrared thermal emission, the disk in scattered light shows similar orientation but different morphology. The scattered-light disk extends from $\sim$75 to $\sim$210 AU in the disk plane with roughly flat surface density. Our Monte Carlo radiative transfer model can well describe the observations with a model disk composed of a mixture of silicates and amorphous carbon. In addition to the obvious brightness asymmetry due to stronger forward scattering, we discover a weak brightness asymmetry along the major axis with the northeast side being 1.3 times brighter than the southwest side at a 3-{\sigma} level.

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

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

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

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.

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.

The Gemini Planet Imager (GPI) has as its science instrument an infrared integral field spectrograph/polarimeter (IFS). Integral field spectrographs are scientificially powerful but require sophisticated data reduction systems. For GPI to achieve its scientific goals of exoplanet and disk characterization, IFS data must be reconstructed into high quality astrometrically and photometrically accurate datacubes in both spectral and polarization modes, via flexible software that is usable by the broad Gemini community. The data reduction pipeline developed by the GPI instrument team to meet these needs is now publicly available following GPI's commissioning. This paper, the first of a series, provides a broad overview of GPI data reduction, summarizes key steps, and presents the overall software framework and implementation. Subsequent papers describe in more detail the algorithms necessary for calibrating GPI data. The GPI data reduction pipeline is open source, available from planetimager.org, and will continue to be enhanced throughout the life of the instrument. It implements an extensive suite of task primitives that can be assembled into reduction recipes to produce calibrated datasets ready for scientific analysis. Angular, spectral, and polarimetric differential imaging are supported. Graphical tools automate the production and editing of recipes, an integrated calibration database manages reference files, and an interactive data viewer customized for high contrast imaging allows for exploration and manipulation of data.

We present the results of both laboratory and on sky astrometric characterization of the Gemini Planet Imager (GPI). This characterization includes measurement of the pixel scale of the integral field spectrograph (IFS), the position of the detector with respect to north, and optical distortion. Two of these three quantities (pixel scale and distortion) were measured in the laboratory using two transparent grids of spots, one with a square pattern and the other with a random pattern. The pixel scale in the laboratory was also estimate using small movements of the artificial star unit (ASU) in the GPI adaptive optics system. On sky, the pixel scale and the north angle are determined using a number of known binary or multiple systems and Solar System objects, a subsample of which had concurrent measurements at Keck Observatory. Our current estimate of the GPI pixel scale is 14.14 $\pm$ 0.01 millarcseconds/pixel, and the north angle is -1.00 $\pm$ 0.03$\deg$. Distortion is shown to be small, with an average positional residual of 0.26 pixels over the field of view, and is corrected using a 5th order polynomial. We also present results from Monte Carlo simulations of the GPI Exoplanet Survey (GPIES) assuming GPI achieves ~1 milliarcsecond relative astrometric precision. We find that with this precision, we will be able to constrain the eccentricities of all detected planets, and possibly determine the underlying eccentricity distribution of widely separated Jovians.

The Gemini Planet Imager (GPI) is a new facility instrument for the Gemini Observatory designed to provide direct detection and characterization of planets and debris disks around stars in the solar neighborhood. In addition to its extreme adaptive optics and corona graphic systems which give access to high angular resolution and high-contrast imaging capabilities, GPI contains an integral field spectrograph providing low resolution spectroscopy across five bands between 0.95 and 2.5 $\mu$m. This paper describes the sequence of processing steps required for the spectro-photometric calibration of GPI science data, and the necessary calibration files. Based on calibration observations of the white dwarf HD 8049B we estimate that the systematic error in spectra extracted from GPI observations is less than 5%. The flux ratio of the occulted star and fiducial satellite spots within coronagraphic GPI observations, required to estimate the magnitude difference between a target and any resolved companions, was measured in the $H$-band to be $\Delta m = 9.23\pm0.06$ in laboratory measurements and $\Delta m = 9.39\pm 0.11$ using on-sky observations. Laboratory measurements for the $Y$, $J$, $K1$ and $K2$ filters are also presented. The total throughput of GPI, Gemini South and the atmosphere of the Earth was also measured in each photometric passband, with a typical throughput in $H$-band of 18% in the non-coronagraphic mode, with some variation observed over the six-month period for which observations were available. We also report ongoing development and improvement of the data cube extraction algorithm.

The Gemini Planet Imager (GPI) is a complex optical system designed to directly detect the self-emission of young planets within two arcseconds of their host stars. After suppressing the starlight with an advanced AO system and apodized coronagraph, the dominant residual contamination in the focal plane are speckles from the atmosphere and optical surfaces. Since speckles are diffractive in nature their positions in the field are strongly wavelength dependent, while an actual companion planet will remain at fixed separation. By comparing multiple images at different wavelengths taken simultaneously, we can freeze the speckle pattern and extract the planet light adding an order of magnitude of contrast. To achieve a bandpass of 20%, sufficient to perform speckle suppression, and to observe the entire two arcsecond field of view at diffraction limited sampling, we designed and built an integral field spectrograph with extremely low wavefront error and almost no chromatic aberration. The spectrograph is fully cryogenic and operates in the wavelength range 1 to 2.4 microns with five selectable filters. A prism is used to produce a spectral resolution of 45 in the primary detection band and maintain high throughput. Based on the OSIRIS spectrograph at Keck, we selected to use a lenslet-based spectrograph to achieve an rms wavefront error of approximately 25 nm. Over 36,000 spectra are taken simultaneously and reassembled into image cubes that have roughly 192x192 spatial elements and contain between 11 and 20 spectral channels. The primary dispersion prism can be replaced with a Wollaston prism for dual polarization measurements. The spectrograph also has a pupil-viewing mode for alignment and calibration.

We present the results of a search for companions to young brown dwarfs in the Taurus and Chamaeleon I star forming regions (1/2-3 Myr). We have used WFPC2 on board HST to obtain F791W and F850LP images of 47 members of these regions that have spectral types of M6-L0 (0.01-0.1 Msun). An additional late-type member of Taurus, FU Tau (M7.25+M9.25), was also observed with adaptive optics at Keck Observatory. We have applied PSF subtraction to the primaries and have searched the resulting images for objects that have colors and magnitudes that are indicative of young low-mass objects. Through this process, we have identified promising candidate companions to 2MASS J04414489+2301513 (rho=0.105"/15 AU), 2MASS J04221332+1934392 (rho=0.05"/7 AU), and ISO 217 (rho=0.03"/5 AU). We reported the discovery of the first candidate in a previous study, showing that it has a similar proper motion as the primary through a comparison of astrometry measured with WFPC2 and Gemini adaptive optics. We have collected an additional epoch of data with Gemini that further supports that result. By combining our survey with previous high-resolution imaging in Taurus, Chamaeleon, and Upper Sco (10 Myr), we measure binary fractions of 14/93 = 0.15+0.05/-0.03 for M4-M6 (0.1-0.3 Msun) and 4/108 = 0.04+0.03/-0.01 for >M6 (<0.1 Msun) at separations of >10 AU. Given the youth and low density of these three regions, the lower binary fraction at later types is probably primordial rather than due to dynamical interactions among association members. The widest low-mass binaries (>100 AU) also appear to be more common in Taurus and Chamaeleon than in the field, which suggests that the widest low-mass binaries are disrupted by dynamical interactions at >10 Myr, or that field brown dwarfs have been born predominantly in denser clusters where wide systems are disrupted or inhibited from forming.

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.

Our general understanding of multiple star and planet formation is primarily based on observations of young multiple systems in low density regions like Tau-Aur and Oph. Since many, if not most, of the stars are born in clusters, observational constraints from young binaries in those environments are fundamental for understanding both the formation of multiple systems and planets in multiple systems throughout the Galaxy. We build upon the largest survey for young binaries in the Orion Nebula Cluster (ONC) which is based on Hubble Space Telescope observations to derive both stellar and circumstellar properties of newborn binary systems in this cluster environment. We present Adaptive Optics spatially-resolved JHKL'-band photometry and K-band R$\sim$\,5000 spectra for a sample of 8 ONC binary systems from this database. We characterize the stellar properties of binary components and obtain a census of protoplanetary disks through K-L' color excess. For a combined sample of ONC binaries including 7 additional systems with NIR spectroscopy from the literature, we derive mass ratio and relative age distributions. We compare the stellar and circumstellar properties of binaries in ONC with those in Tau-Aur and Oph from samples of binaries with stellar properties derived for each component from spectra and/or visual photometry and with a disk census obtained through K-L color excess. The mass ratio distribution of ONC binaries is found to be indistinguishable from that of Tau-Aur and, to some extent, to that of Oph in the separation range 85-560\,AU and for primary mass in the range 0.15 to 0.8\,M$_{\sun}$.A trend toward a lower mass ratio with larger separation is suggested in ONC binaries which is not seen in Tau-Aur binaries.The components of ONC binaries are found to be significantly more coeval than the overall ONC population and as coeval as components of binaries in Tau-Aur and Oph[...]

We present the first high-angular resolution survey for multiple systems among very low-mass stars and brown dwarfs in the Hyades open cluster. Using the Keck\,II adaptive optics system, we observed a complete sample of 16 objects with estimated masses $\lesssim$0.1 Msun. We have identified three close binaries with projected separation $\lesssim$0.11", or $\lesssim$5 AU. A number of wide, mostly faint candidate companions are also detected in our images, most of which are revealed as unrelated background sources based on astrometric and/or photometric considerations. The derived multiplicity frequency, 19+13/-6 % over the 2-350 AU range, and the rarity of systems wider than 10 AU are both consistent with observations of field very low-mass objects. In the limited 3-50 AU separation range, the companion frequency is essentially constant from brown dwarfs to solar-type stars in the Hyades cluster, which is also in line with our current knowledge for field stars. Combining the binaries discovered in this surveys with those already known in the Pleiades cluster reveals that very low-mass binaries in open clusters, as well as in star-forming regions, are skewed toward lower mass ratios ($0.6 \lesssim q \lesssim 0.8$) than are their field counterparts, a result that cannot be accounted for by selection effects. Although the possibility of severe systematic errors in model-based mass estimates for very low-mass stars cannot be completely excluded, it is unlikely to explain this difference. We speculate that this trend indicates that surveys among very low-mass field stars may have missed a substantial population of intermediate mass ratio systems, implying that these systems are more common and more diverse than previously thought.

Studies of solar-type binaries have found coplanarity between the equatorial and orbital planes of systems with $<$40 AU separation. By comparison, the alignment of the equatorial and orbital axes in the substellar regime, and the associated implications for formation theory, are relatively poorly constrained. Here we present the discovery of the rotation period of 3.32 $\pm$ 0.15 hours from 2MASS J0746+20A - the primary component of a tight (2.7 AU) ultracool dwarf binary system (L0+L1.5). The newly discovered period, together with the established period via radio observations of the other component, and the well constrained orbital parameters and rotational velocity measurements, allow us to infer alignment of the equatorial planes of both components with the orbital plane of the system to within 10 degrees. This result suggests that solar-type binary formation mechanisms may extend down into the brown dwarf mass range, and we consider a number of formation theories that may be applicable in this case. This is the first such observational result in the very low mass binary regime. In addition, the detected period of 3.32 $\pm$ 0.15 hours implies that the reported radio period of 2.07 $\pm$ 0.002 hours is associated with the secondary star, not the primary, as was previously claimed. This in turn refutes the claimed radius of 0.78 $\pm$ 0.1 $R_{J}$ for 2MASS J0746+20A, which we demonstrate to be 0.99 $\pm$ 0.03 $R_{J}$.

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.

To image faint substellar companions obscured by the stellar halo and speckles, scattered light from the bright primary star must be removed in hardware or software. We apply the "locally-optimized combination of images" (LOCI) algorithm to 1-minute Keck Observatory snapshots of GKM dwarfs in the Hyades using source diversity to determine the most likely PSF. We obtain a mean contrast of 10^{-2} at 0.01", 10^{-4} at <1", and 10^{-5} at 5". New brown dwarf and low-mass stellar companions to Hyades primaries are found in a third of the 84 targeted systems. This campaign shows the efficacy of LOCI on snapshot imaging as well as on bright wide binaries with off-axis LOCI, reaching contrasts sufficient for imaging 625-Myr late-L/early-T dwarfs purely in post-processing.

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 on new high-resolution imaging and spectroscopy on the multiple T Tauri star system V773 Tau over the 2003 -- 2009 period. With these data we derive relative astrometry, photometry between the A and B components, and radial velocity (RV) of the A-subsystem components. Combining these new data with previously published astrometry and RVs, we update the relative A-B orbit model. This updated orbit model, the known system distance, and A subsystem parameters yields a dynamical mass for the B component for the first time. Remarkably the derived B dynamical mass is in the range of 1.7 -- 3.0 M$_\sun$. This is much higher than previous estimates, and suggests that like A, B is also a multiple stellar system. Among these data, spatially-resolved spectroscopy provide new insight into the nature of the B component. Similar to A, these near-IR spectra indicate that the dominant source in B is of mid-K spectral type. If B is in fact a multiple star system as suggested by the dynamical mass estimate, the simplest assumption is that B is composed of similar $\sim$ 1.2 M$_\sun$ PMS stars in a close ($<$ 1 AU) binary system. This inference is supported by line-shape changes in near-IR spectroscopy of B, tentatively interpreted as changing RV among components in V773 Tau B. Relative photometry indicate that B is highly variable in the near-IR. The most likely explanation for this variability is circum-B material resulting in variable line-of-sight extinction. The distribution of this material must be significantly affected by both the putative B multiplicity, and the A-B orbit.

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.

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.

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 identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7.

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

Owing to its extremely high luminosity and long duration, SN2006gy radiated more energy in visual light than any other known SN. Two hypotheses to explain its high luminosity -- that it was powered by shock interaction with CSM as implied by its Type IIn spectrum, or that it was fueled by radioactive decay from a large mass of 56Ni synthesized in a pair-instability SN -- predicted different late-time properties. Here we present observations of SN2006gy obtained more than a year after discovery. We were unable to detect it at visual wavelengths, but clear near-IR K and H-band detections show that it is still at least as luminous as the peak of a normal SN II. We also present spectra giving an upper limit to the late-time Ha luminosity of about 1e39 erg/s. Based on the weak late-time Ha, X-ray, and radio emission, combined with the difficulty of explaining the shift to IR wavelengths, we can rule out ongoing CSM interaction as the primary late-time power source. Instead, we propose that the evolution of SN2006gy is consistent with one of two possible scenarios: (1) A pair-instability SN plus modest CSM interaction, where the radioactive decay luminosity shifts to the IR because of dust formation. (2) An IR echo, where radiation emitted during peak luminosity heats a pre-existing dust shell at radii near 1 light year, requiring the progenitor star to have ejected another 10 Msun shell about 1500 yr before the SN.

We present a review of precursor observing programs for the SIM PlanetQuest Key project devoted to detecting Jupiter mass planets around young stars. In order to ensure that the stars in the sample are free of various sources of astrometric noise that might impede the detection of planets, we have initiated programs to collect photometry, high contrast images, interferometric data and radial velocities for stars in both the Northern and Southern hemispheres. We have completed a high contrast imaging survey of target stars in Taurus and the Pleiades and found no definitive common proper motion companions within one arcsecond (140 AU) of the SIM targets. Our radial velocity surveys have shown that many of the target stars in Sco-Cen are fast rotators and a few stars in Taurus and the Pleiades may have sub-stellar companions. Interferometric data of a few stars in Taurus show no signs of stellar or sub-stellar companions with separations of <5 mas. The photometric survey suggests that approximately half of the stars initially selected for this program are variable to a degree (1 sigma>0.1 mag) that would degrade the astrometric accuracy achievable for that star. While the precursor programs are still a work in progress, we provide a comprehensive list of all targets ranked according to their viability as a result of the observations taken to date. By far, the observable that moves the most targets from the SIM-YSO program is photometric variability.

We surveyed thirteen very low mass (VLM; M < 0.2 M_sun) objects in the Taurus star-forming region using near-infrared diffraction-limited imaging techniques on the W.M. Keck I 10 m telescope. Of these thirteen, five were found to be binary, with separations ranging from 0.04" to 0.6" and flux ratios from 1.4 to 3.7. In all cases, the companions are likely to be physically associated with the primaries (probability > 4-sigma). Using the theoretical models of Baraffe et al. (1998), we find that all five new companions, as well as one of the primaries, are likely brown dwarfs. The discovery of these systems therefore increases the total number of known, young VLM binaries by ~50%. These new systems, along with other young VLM binaries from the literature, have properties that differ significantly from older field VLM binaries in that the young systems have wider separations and lower mass ratios, supporting the idea that VLM binaries undergo significant dynamical evolution ~5 - 10 Myr after their formation. The range of separations of these binaries, four of which are over 30 AU, argues against the ejection scenario of brown dwarf formation. While several of the young, VLM binaries discovered in this study have lower binding energies than the previously suggested minimum for VLM binaries, the apparent minimum is still significantly higher than that found among higher mass binaries. We suggest that this discrepancy may be due to the small mass of a VLM binary relative to the average perturbing star, leading to more substantial changes in their binding energy over time.

We present the results of a five year monitoring campaign of the close binary TWA 5Aab in the TW Hydrae association, using speckle and adaptive optics on the W.M. Keck 10 m telescopes. These measurements were taken as part of our ongoing monitoring of pre-main sequence (PMS) binaries in an effort to increase the number of dynamically determined PMS masses and thereby calibrate the theoretical PMS evolutionary tracks. Our observations have allowed us to obtain the first determination of this system's astrometric orbit. We find an orbital period of 5.94 +- 0.09 years and a semi-major axis of 0.066" +- 0.005". Combining these results with a kinematic distance, we calculate a total mass of 0.71 +- 0.14 M_sun (D/44 pc)^3. for this system. This mass measurement, as well as the estimated age of this system, are consistent to within 2$\sigma$ of all theoretical models considered. In this analysis, we properly account for correlated uncertainties, and show that while these correlations are generally ignored, they increase the formal uncertainties by up to a factor of five and therefore are important to incorporate. With only a few more years of observation, this type of measurement will allow the theoretical models to be distinguished.

The goal of this study is to obtain accurate estimates for the individual masses of the components of the tight binary system T Tau S in order to settle the ongoing debate on the nature of T Tau Sa, a so-called infrared companion. We take advantage of the fact that T Tau S belongs to a triple system composed of two hierarchical orbits to simultaneously analyze the motion of T Tau Sb in the rest frames of T Tau Sa and T Tau N. With this method, it is possible to pinpoint the location of the center of mass of T Tau S and, thereby, to determine individual masses for T Tau Sa and T Tau Sb with no prior assumption about the mass/flux ratio of the system. This improvement over previous studies of the system results in much better constraints on orbital parameters. We find individual masses of 2.73+/-0.31 Msun for T Tau Sa and of 0.61+/-0.17 Msun for T Tau Sb (in agreement with its early-M spectral type), including the uncertainty on the distance to the system. These are among the most precise estimates of the mass of any Pre-Main Sequence star, a remarkable result since this is the first system in which individual masses of T Tauri stars can be determined from astrometry only. This model-independent analysis confirms that T Tau Sa is an intermediate-mass star, presumably a very young Herbig Ae star, that may possess an almost edge-on disk.

The explosion that results in a cosmic gamma-ray burst (GRB) is thought to produce emission from two physical processes -- the activity of the central engine gives rise to the high-energy emission of the burst through internal shocking and the subsequent interaction of the flow with the external environment produces long-wavelength afterglow. While afterglow observations continue to refine our understanding of GRB progenitors and relativistic shocks, gamma-ray observations alone have not yielded a clear picture of the origin of the prompt emission nor details of the central engine. Only one concurrent visible-light transient has been found and was associated with emission from an external shock. Here we report the discovery of infrared (IR) emission contemporaneous with a GRB, beginning 7.2 minutes after the onset of GRB 041219a. Our robotic telescope acquired 21 images during the active phase of the burst, yielding the earliest multi-colour observations of any long-wavelength emission associated with a GRB. Analysis of an initial IR pulse suggests an origin consistent with internal shocks. This opens a new possibility to study the central engine of GRBs with ground-based observations at long wavelengths.

HD 199143 and HD 358623 (BD-17 6128) are two sets of binary stars which are physically associated and 48 pc from Earth. We present heliocentric radial velocities and high lithium abundances which establish these stars as members of the 12 Myr-old Beta Pictoris Moving Group. We also present mid-IR photometric measurements which show no firm evidence for warm dust around all four stars.

If enough of their Lyman limit continuum escapes, star-forming galaxies could be significant contributors to the cosmic background of ionizing photons. To investigate this possibility, we obtained the first deep imaging in the far ultraviolet of eleven bright blue galaxies at intermediate redshift (z=1.1--1.4). NO Lyman continuum emission was detected. Sensitive, model-independent, upper limits of typically 2 x 10**-19 erg/sec/cm2/Ang were obtained for the ionizing flux escaping from these normal galaxies. This corresponds to lower limits on the observed ratio of 1500 to 700Ang flux of 150 up to 1000. Based on a wide range of stellar synthesis models, this suggests that less than 6%, down to less than 1%, of the available ionizing flux emitted by hot stars is escaping these galaxies. The magnitude of this spectral break at the Lyman l imit confirms that the basic premise of `Lyman break' searches for galaxies at high redshift can also be applied at intermediate redshifts. This implies that the integrated contribution of galaxies to the UV cosmic background at z around 1.2 is less than 15%, and may be less than 2%.

A significant population of nearby stars have strong far-infrared excesses, now known to be due to circumstellar dust in regions analogous to the Kuiper Belt of our solar system, though orders of magnitude more dense. Recent sub-mm and mm imaging of these systems resolves the circumstellar dust and reveals complex structures, often in the form of rings with azimuthal non-axisymmetric variations. This structure might well be due to the presence of embedded brown dwarfs or planets. We have carried out deep adaptive optics imaging of two nearby stars with such asymmetric dust: Epsilon Eridani and Vega. Ten and seven candidate companions were seen in and near the dust rings of Epsilon Eridani and Vega respectively, but second-epoch proper motion measurements indicate that all are background objects. Around these two stars we can thus exclude planetary companions at spatial scales comparable to the radius of the dust structures to a level of MK=24, corresponding to 5 Jupiter masses, for Epsilon Eridani, and MK=19-21, corresponding to 6-8 Jupiter masses, for Vega.