We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

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.

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

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

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.

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

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

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

Ultra-cool brown dwarfs offer a unique window into understanding substellar atmospheric physics and chemistry. Their strong molecular absorption bands at infrared wavelengths, Jupiter-like radii, cool temperatures, and lack of complicating stellar irradiation, make them ideal test-beds for understanding Jovian-like atmospheres. Here we report the findings of a uniform atmospheric retrieval analysis on a set of 14 Y and T-dwarfs observed with the Hubble Space Telescope Wide Field Camera 3 instrument. From our retrieval analysis, we find the temperature-structures to be largely consistent with radiative-convective equilibrium in most objects. We also determine the abundances of water, methane, and ammonia and upper limits on the alkali metals sodium and potassium. The constraints on water and methane are consistent with predictions from chemical equilibrium models, while those of ammonia may be affected by vertical disequilibrium mixing, consistent with previous works. Our key result stems from the constraints on the alkali metal abundances where we find their continued depletion with decreasing effective temperature, consistent with the trend identified in a previous retrieval analysis on a sample of slightly warmer late T-dwarfs in Line et al. (2017). These constraints show that the previously observed Y-J color trend across the T/Y transition is most likely due to the depletion of these metals in accordance with predictions from equilibrium condensate rainout chemistry. Finally, we simulate future James Webb Space Telescope observations of ultra-cool dwarfs and find that the NIRSpec PRISM offers the best chance at developing high-precision constraints on fundamental atmospheric characteristics.

We present the results of the largest $L^{\prime}$ ($3.8~\mu$m) direct imaging survey for exoplanets to date, the Large Binocular Telescope Interferometer (LBTI) Exozodi Exoplanet Common Hunt (LEECH). We observed 98 stars with spectral types from B to M. Cool planets emit a larger share of their flux in $L^{\prime}$ compared to shorter wavelengths, affording LEECH an advantage in detecting low-mass, old, and cold-start giant planets. We emphasize proximity over youth in our target selection, probing physical separations smaller than other direct imaging surveys. For FGK stars, LEECH outperforms many previous studies, placing tighter constraints on the hot-start planet occurrence frequency interior to $\sim20$ au. For less luminous, cold-start planets, LEECH provides the best constraints on giant-planet frequency interior to $\sim20$ au around FGK stars. Direct imaging survey results depend sensitively on both the choice of evolutionary model (e.g., hot- or cold-start) and assumptions (explicit or implicit) about the shape of the underlying planet distribution, in particular its radial extent. Artificially low limits on the planet occurrence frequency can be derived when the shape of the planet distribution is assumed to extend to very large separations, well beyond typical protoplanetary dust-disk radii ($\lesssim50$ au), and when hot-start models are used exclusively. We place a conservative upper limit on the planet occurrence frequency using cold-start models and planetary population distributions that do not extend beyond typical protoplanetary dust-disk radii. We find that $\lesssim90\%$ of FGK systems can host a 7 to 10 $M_{\mathrm{Jup}}$ planet from 5 to 50 au. This limit leaves open the possibility that planets in this range are common.

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.

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.

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 report 885$\mu$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from $\sim$0.3-20M$_{\oplus}$. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in $M_{dust}$ to $M_{star}$, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early M-stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas:dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M-dwarfs.

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.

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 conducted a 12-month monitoring campaign of 33 T Tauri stars (TTS) in Taurus. Our goal was to monitor objects that possess a disk but have a weak Halpha line, a common accretion tracer for young stars, to determine whether they host a passive circumstellar disk. We used medium-resolution optical spectroscopy to assess the objects' accretion status and to measure the Halpha line. We found no convincing example of passive disks; only transition disk and debris disk systems in our sample are non-accreting. Among accretors, we find no example of flickering accretion, leading to an upper limit of 2.2% on the duty cycle of accretion gaps assuming that all accreting TTS experience such events. Combining literature results with our observations, we find that the reliability of traditional Halpha-based criteria to test for accretion is high but imperfect, particularly for low-mass TTS. We find a significant correlation between stellar mass and the full width at 10 per cent of the peak (W10%) of the Halpha line that does not seem to be related to variations in free-fall velocity. Finally, our data reveal a positive correlation between the Halpha equivalent width and its W10%, indicative of a systematic modulation in the line profile whereby the high-velocity wings of the line are proportionally more enhanced than its core when the line luminosity increases. We argue that this supports the hypothesis that the mass accretion rate on the central star is correlated with the Halpha W10% through a common physical mechanism.

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.

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.

Using the Gemini Planet Imager (GPI), we have resolved the circumstellar debris disk around HD 111520 at a projected range of ~30-100 AU in both total and polarized $H$-band intensity. The disk is seen edge-on at a position angle of ~165$^{\circ}$ along the spine of emission. A slight inclination or asymmetric warping are covariant and alters the interpretation of the observed disk emission. We employ 3 point spread function (PSF) subtraction methods to reduce the stellar glare and instrumental artifacts to confirm that there is a roughly 2:1 brightness asymmetry between the NW and SE extension. This specific feature makes HD 111520 the most extreme examples of asymmetric debris disks observed in scattered light among similar highly inclined systems, such as HD 15115 and HD 106906. We further identify a tentative localized brightness enhancement and scale height enhancement associated with the disk at ~40 AU away from the star on the SE extension. We also find that the fractional polarization rises from 10 to 40% from 0.5" to 0.8" from the star. The combination of large brightness asymmetry and symmetric polarization fraction leads us to believe that an azimuthal dust density variation is causing the observed asymmetry.

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 results of ALMA band 7 observations of dust and CO gas in the disks around 7 objects with spectral types ranging between M5.5 and M7.5 in Upper Scorpius OB1, and one M3 star in Ophiuchus. We detect unresolved continuum emission in all but one source, and the $^{12}$CO J=3-2 line in two sources. We constrain the dust and gas content of these systems using a grid of models calculated with the radiative transfer code MCFOST, and find disk dust masses between 0.1 and 1 M$_\oplus$, suggesting that the stellar mass / disk mass correlation can be extrapolated for brown dwarfs with masses as low as 0.05 M$_\odot$. The one disk in Upper Sco in which we detect CO emission, 2MASS J15555600, is also the disk with warmest inner disk as traced by its H - [4.5] photometric color. Using our radiative transfer grid, we extend the correlation between stellar luminosity and mass-averaged disk dust temperature originally derived for stellar mass objects to the brown dwarf regime to $\langle T_{dust} \rangle \approx 22 (L_{*} /L_{\odot})^{0.16} K$, applicable to spectral types of M5 and later. This is slightly shallower than the relation for earlier spectral type objects and yields warmer low-mass disks. The two prescriptions cross at 0.27 L$_\odot$, corresponding to masses between 0.1 and 0.2 M$_\odot$ depending on age.

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

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 the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$\deg$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet.

We have observed and spatially resolved a set of seven A-type stars in the nearby Ursa Major moving group with the Classic, CLIMB, and PAVO beam combiners on the CHARA Array. At least four of these stars have large rotational velocities ($v \sin i$ $\gtrsim$ 170 $\mathrm{km~s^{-1}}$) and are expected to be oblate. These interferometric measurements, the stars' observed photometric energy distributions, and $v \sin i$ values are used to computationally construct model oblate stars from which stellar properties (inclination, rotational velocity, and the radius and effective temperature as a function of latitude, etc.) are determined. The results are compared with MESA stellar evolution models (Paxton et al. 2011, 2013) to determine masses and ages. The value of this new technique is that it enables the estimation of the fundamental properties of rapidly rotating stars without the need to fully image the star. It can thus be applied to stars with sizes comparable to the interferometric resolution limit as opposed to those that are several times larger than the limit. Under the assumption of coevality, the spread in ages can be used as a test of both the prescription presented here and the MESA evolutionary code for rapidly rotating stars. With our validated technique, we combine these age estimates and determine the age of the moving group to be 414 $\pm$ 23 Myr, which is consistent with, but much more precise than previous estimates.

We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband imaging polarimetry were obtained during the commissioning of GPI. In our broadband imaging polarimetry observations, we detect the disk only in total intensity and find asymmetries in the morphology of the disk between the southeast and northwest sides. The southeast side of the disk exhibits a bump at 1$''$ (10 AU projected separation) that is three times more vertically extended and three times fainter in peak surface brightness than the northwest side at similar separations. This part of the disk is also vertically offset by 69$\pm$30 mas to the northeast at 1$''$ when compared to the established disk mid-plane and consistent with prior ALMA and Hubble Space Telescope/STIS observations. We see hints that the southeast bump might be a result of detecting a horizontal sliver feature above the main disk that could be the disk backside. Alternatively when including the morphology of the northwest side, where the disk mid-plane is offset in the opposite direction $\sim$50 mas between 0$.''$4 and 1$.''$2, the asymmetries suggest a warp-like feature. Using our integral field spectroscopy data to search for planets, we are 50% complete for $\sim$4 $M_\mathrm{Jup}$ planets at 4 AU. We detect a source, resolved only along the disk plane, that could either be a candidate planetary mass companion or a compact clump in the disk.

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

Context. Astrometric monitoring of directly-imaged exoplanets allows the study of their orbital parameters and system architectures. Because most directly-imaged planets have long orbital periods (>20 AU), accurate astrometry is challenging when based on data acquired on timescales of a few years and usually with different instruments. The LMIRCam camera on the LBT is being used for the LEECH survey to search for and characterize young and adolescent exoplanets in L' band, including their system architectures. Aims. We first aim to provide a good astrometric calibration of LMIRCam. Then, we derive new astrometry, test the predictions of the orbital model of 8:4:2:1 mean motion resonance proposed by Go\'zdziewski & Migaszewski, and perform new orbital fitting of the HR 8799 bcde planets. We also present deep limits on a putative fifth planet interior to the known planets. Methods. We use observations of HR 8799 and the Theta1 Ori C field obtained during the same run in October 2013. Results. We first characterize the distortion of LMIRCam. We determine a platescale and a true north orientation for the images of 10.707 +/- 0.012 mas/pix and -0.430 +/- 0.076 deg, respectively. The errors on the platescale and true north orientation translate into astrometric accuracies at a separation of 1 of 1.1 mas and 1.3 mas, respectively. The measurements for all planets are usually in agreement within 3 sigma with the ephemeris predicted by Go\'zdziewski & Migaszewski. The orbital fitting based on the new astrometric measurements favors an architecture for the planetary system based on 8:4:2:1 mean motion resonance. The detection limits allow us to exclude a fifth planet slightly brighter/more massive than HR 8799 b at the location of the 2:1 resonance with HR 8799 e (~9.5 AU) and about twice as bright as HR 8799 cde at the location of the 3:1 resonance with HR 8799 e (~7.5 AU).

We present a large-scale, volume-limited companion survey of 245 late-K to mid-M (K7-M6) dwarfs within 15 pc. Infrared adaptive optics (AO) data were analysed from the Very Large Telescope, Subaru Telescope, Canada-France-Hawaii Telescope, and MMT Observatory to detect close companions to the sample from $\sim$1 au to 100 au, while digitised wide-field archival plates were searched for wide companions from $\sim$100 au to 10,000 au. With sensitivity to the bottom of the main sequence over a separation range of 3 au to 10,000 au, multiple AO and wide-field epochs allow us to confirm candidates with common proper motions, minimize background contamination, and enable a measurement of comprehensive binary statistics. We detected 65 co-moving stellar companions and find a companion star fraction of $23.5 \pm 3.2$ per cent over the 3 au to 10,000 au separation range. The companion separation distribution is observed to rise to a higher frequency at smaller separations, peaking at closer separations than measured for more massive primaries. The mass ratio distribution across the $q = 0.2 - 1.0$ range is flat, similar to that of multiple systems with solar-type primaries. The characterisation of binary and multiple star frequency for low-mass field stars can provide crucial comparisons with star forming environments and hold implications for the frequency and evolutionary histories of their associated disks and planets.

With the discovery of Y dwarfs by the WISE mission, the population of field brown dwarfs now extends to objects with temperatures comparable to those of Solar System planets. To investigate the atmospheres of these newly identified brown dwarfs, we have conducted a pilot study monitoring an initial sample of three late T-dwarfs (T6.5, T8 and T8.5) and one Y-dwarf (Y0) for infrared photometric variability at multiple epochs. With J-band imaging, each target was observed for a period of 1.0h to 4.5h per epoch, which covers a significant fraction of the expected rotational period. These measurements represent the first photometric monitoring for these targets. For three of the four targets (2M1047, Ross 458C and WISE0458), multi-epoch monitoring was performed, with the time span between epochs ranging from a few hours to ~2 years. During the first epoch, the T8.5 target WISE0458 exhibited variations with a remarkable min-to-max amplitude of 13%, while the second epoch light curve taken ~2 years later did not note any variability to a 3% upper limit. With an effective temperature of ~600 K, WISE0458 is the coldest variable brown dwarf published to-date, and combined with its high and variable amplitude makes it a fascinating target for detailed follow-up. The three remaining targets showed no significant variations, with a photometric precision between 0.8% and 20.0%, depending on the target brightness. Combining the new results with previous multi-epoch observations of brown dwarfs with spectral types of T5 or later, the currently identified variables have locations on the colour-colour diagram better matched by theoretical models incorporating cloud opacities rather than cloud-free atmospheres. This preliminary result requires further study to determine if there is a definitive link between variability among late-T dwarfs and their location on the colour-colour diagram.

We report the discovery of a wide co-moving substellar companion to the nearby ($D=67.5\pm1.1$ pc) A3V star $\zeta$ Delphini based on imaging and follow-up spectroscopic observations obtained during the course of our Volume-limited A-Star (VAST) multiplicity survey. $\zeta$ Del was observed over a five-year baseline with adaptive optics, revealing the presence of a previously-unresolved companion with a proper motion consistent with that of the A-type primary. The age of the $\zeta$ Del system was estimated as $525\pm125$ Myr based on the position of the primary on the colour-magnitude and temperature-luminosity diagrams. Using intermediate-resolution near-infrared spectroscopy, the spectrum of $\zeta$ Del B is shown to be consistent with a mid-L dwarf (L$5\pm2$), at a temperature of $1650\pm200$ K. Combining the measured near-infrared magnitude of $\zeta$ Del B with the estimated temperature leads to a model-dependent mass estimate of $50\pm15$ M$_{\rm Jup}$, corresponding to a mass ratio of $q=0.019\pm0.006$. At a projected separation of $910\pm14$ au, $\zeta$ Del B is among the most widely-separated and extreme-mass ratio substellar companions to a main-sequence star resolved to-date, providing a rare empirical constraint of the formation of low-mass ratio companions at extremely wide separations.

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.

With Herschel/PACS 134 low mass members of the Taurus star-forming region spanning the M4-L0 spectral type range and covering the transition from low mass stars to brown dwarfs were observed. Combining the new Herschel results with other programs, a total of 150 of the 154 M4-L0 Taurus members members have observations with Herschel. Among the 150 targets, 70um flux densities were measured for 7 of the 7 ClassI objects, 48 of the 67 ClassII members, and 3 of the 76 ClassIII targets. For the detected ClassII objects, the median 70um flux density level declines with spectral type, however, the distribution of excess relative to central object flux density does not change across the stellar/substellar boundary in the M4-L0 range. Connecting the 70um TBOSS values with the results from K0-M3 ClassII members results in the first comprehensive census of far-IR emission across the full mass spectrum of the stellar and substellar population of a star-forming region, and the median flux density declines with spectral type in a trend analogous to the flux density decline expected for the central objects. SEDs were constructed for all TBOSS targets covering the optical to far-IR range and extending to the submm/mm for a subset of sources. Based on an initial exploration of the impact of different physical parameters; inclination, scale height and flaring have the largest influence on the PACS flux densities. From the 24um to 70um spectral index of the SEDs, 5 new candidate transition disks were identified. The steep 24um to 70um slope for a subset of 8 TBOSS targets may be an indication of truncated disks in these systems.Two examples of mixed pair systems that include secondaries with disks were measured. Finally, comparing the TBOSS results with a Herschel study of Ophiuchus brown dwarfs reveals a lower fraction of disks around the Taurus substellar population.

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.

In Spring 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its $\sim$130-night campaign from the Large Binocular Telescope (LBT) atop Mt Graham, Arizona. This survey benefits from the many technological achievements of the LBT, including two 8.4-meter mirrors on a single fixed mount, dual adaptive secondary mirrors for high Strehl performance, and a cold beam combiner to dramatically reduce the telescope's overall background emissivity. LEECH neatly complements other high-contrast planet imaging efforts by observing stars at L' (3.8 $\mu$m), as opposed to the shorter wavelength near-infrared bands (1-2.4 $\mu$m) of other surveys. This portion of the spectrum offers deep mass sensitivity, especially around nearby adolescent ($\sim$0.1-1 Gyr) stars. LEECH's contrast is competitive with other extreme adaptive optics systems, while providing an alternative survey strategy. Additionally, LEECH is characterizing known exoplanetary systems with observations from 3-5$\mu$m in preparation for JWST.

We present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses a new integral field polarimetry architecture to provide high contrast linear polarimetry with minimal systematic biases between the orthogonal polarizations. We describe the design, data reduction methods, and performance of polarimetry with GPI. Point spread function subtraction via differential polarimetry suppresses unpolarized starlight by a factor of over 100, and provides sensitivity to circumstellar dust reaching the photon noise limit for these observations. In the case of the circumstellar disk around HR 4796A, GPI's advanced adaptive optics system reveals the disk clearly even prior to PSF subtraction. In polarized light, the disk is seen all the way in to its semi-minor axis for the first time. The disk exhibits surprisingly strong asymmetry in polarized intensity, with the west side >9 times brighter than the east side despite the fact that the east side is slightly brighter in total intensity. Based on a synthesis of the total and polarized intensities, we now believe that the west side is closer to us, contrary to most prior interpretations. Forward scattering by relatively large silicate dust particles leads to the strong polarized intensity on the west side, and the ring must be slightly optically thick in order to explain the lower brightness in total intensity there. These findings suggest that the ring is geometrically narrow and dynamically cold, perhaps shepherded by larger bodies in the same manner as Saturn's F ring.

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.

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) combines extreme adaptive optics, an integral field spectrograph, and a high performance coronagraph to directly image extrasolar planets in the near-infrared. Because the coronagraph blocks most of the light from the star, it prevents the properties of the host star from being measured directly. Instead, satellite spots, which are created by diffraction from a square grid in the pupil plane, can be used to locate the star and extract its spectrum. We describe the techniques implemented into the GPI Data Reduction Pipeline to measure the properties of the satellite spots and discuss the precision of the reconstructed astrometry and spectrophotometry of the occulted star. We find the astrometric precision of the satellite spots in an $H$-band datacube to be $0.05$ pixels and is best when individual satellite spots have a signal to noise ratio (SNR) of $> 20$. In regards to satellite spot spectrophotometry, we find that the total flux from the satellite spots is stable to $\sim 7\%$ and scales linearly with central star brightness and that the shape of the satellite spot spectrum varies on the $2\%$ level.

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.

Using SofI on the 3.5m New Technology Telescope, we have conducted an extensive near-infrared monitoring survey of an unbiased sample of 69 brown dwarfs spanning the L0 and T8 spectral range, with at least one example of each spectral type. Each target was observed for a 2-4 hour period in the Js-band, and the median photometric precision of the data is ~0.7%. A total of 14 brown dwarfs were identified as variables with min-to-max amplitudes ranging from 1.7% to 10.8% over the observed duration. All variables satisfy a statistical significance threshold with a p-value <5% based on comparison with the median reference star light curve. Approximately half of the variables show sinusoidal amplitude variations similar to 2M2139, and the remainder shows short timescale evolving light curves similar to SIMP0136. The L/T transition has been suggested to be a region of a higher degree of variability if patchy clouds are present and this survey was designed to test the patchy cloud model with photometric monitoring of both the L/T transition and non-transition brown dwarfs. Considering the targets identified as variable with the same statistical threshold, the measured variability frequency of 13^{+10}_{-4}% for the L7 -- T4 transition region is indistinguishable from that of the earlier spectral types (32^{+11}_{-8}%), the later spectral types (13^{+10}_{-4}%), or the combination of all non-transition region brown dwarfs (21^{+7}_{-5}%). The variables are not concentrated at the transition, in a specific colour, or with binaries. We note that of the systems previously monitored for variability only ~60% maintained the state of variability (variable or constant), with the remaining switching states. The 14 variables include nine newly identified variables that will provide important systems for follow-up multi-wavelength monitoring to further investigate brown dwarf atmosphere physics.

We present results from our ongoing program to map the visual orbits of pre-main sequence binaries in the Taurus star forming region using adaptive optics imaging at the Keck Observatory. We combine our results with measurements reported in the literature to analyze the orbital motion for each binary. We present preliminary orbits for DF Tau, T Tau S, ZZ Tau, and the Pleiades binary HBC 351. Seven additional binaries show curvature in their relative motion. Currently, we can place lower limits on the orbital periods for these systems; full solutions will be possible with more orbital coverage. Five other binaries show motion that is indistinguishable from linear motion. We suspect that these systems are bound and might show curvature with additional measurements in the future. The observations reported herein lay critical groundwork toward the goal of measuring precise masses for low-mass pre-main sequence stars.

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.

In February 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its 100-night campaign from the Large Binocular Telescope atop Mount Graham in Arizona. LEECH nearly complements other high-contrast planet imaging efforts by observing stars in L' band (3.8 microns) as opposed to the shorter wavelength near-infrared bands (1-2.3 microns). This part of the spectrum offers deeper mass sensitivity for intermediate age (several hundred Myr-old) systems, since their Jovian-mass planets radiate predominantly in the mid-infrared. In this proceedings, we present the science goals for LEECH and a preliminary contrast curve from some early data.

With a combination of adaptive optics imaging and a multi-epoch common proper motion search, we have conducted a large volume-limited (D $\le$ 75 pc) multiplicity survey of A-type stars, sensitive to companions beyond 30 au. The sample for the Volume-limited A-STar (VAST) survey consists of 435 A-type stars: 363 stars were observed with adaptive optics, 228 stars were searched for wide common proper motion companions and 156 stars were measured with both techniques. The projected separation coverage of the VAST survey extends from 30 to 45,000 au. A total of 137 stellar companions were resolved, including 64 new detections from the VAST survey, and the companion star fraction, projected separation distribution and mass ratio distribution were measured. The separation distribution forms a log-normal distribution similar to the solar-type binary distribution, but with a peak shifted to a significantly wider value of 387 (+132,-98) au. Integrating the fit to the distribution over the 30 to 10,000 au observed range, the companion star fraction for A-type stars is estimated as 33.8%+-2.6%. The mass ratio distribution of closer (<125 au) binaries is distinct from that of wider systems, with a flat distribution for close systems and a distribution that tends towards smaller mass ratios for wider binaries. Combining this result with previous spectroscopic surveys of A-type stars gives an estimate of the total companion star fraction of 68.9%+-7.0%. The most complete assessment of higher order multiples was estimated from the 156-star subset of the VAST sample with both adaptive optics and common proper motion measurements, combined with a literature search for companions, yielding a lower limit on the frequency of single, binary, triple, quadruple and quintuple A-type star systems of 56.4 (-4.0,+3.8), 32.1 (-3.5,+3.9), 9.0 (-1.8,+2.8), 1.9 (-0.6,+1.8) and 0.6 (-0.2,+1.4) per cent, respectively.

A set of six debris disk candidates identified with IRAS or WISE excesses were observed at either 350 um or 450 um with the CSO. Five of the targets - HIP 51658, HIP 68160, HIP 73512, HIP 76375, and HIP 112460 - have among the largest measured excess emission from cold dust from IRAS in the 25-100 um bands. Single temperature blackbody fits to the excess dust emission of these sources predict 350-450 um fluxes above 240 mJy. The final target - HIP 73165 - exhibits weak excess emission above the stellar photosphere from WISE measurements at 22 um, indicative of a population of warm circumstellar dust. None of the six targets were detected, with 3 sigma upper limits ranging from 51-239 mJy. These limits are significantly below the expected fluxes from SED fitting. Two potential causes of the null detections were explored - companion stars and contamination. To investigate the possible influence of companion stars, imaging data were analyzed from new AO data from the MMT and archival HST, NIRI, and POSS/2MASS data. The images are sensitive to all stellar companions beyond a radius of 1-94 AU. One target is identified as a binary system, but with a separation too large to impact the disk. While the gravitational effects of a companion do not appear to provide an explanation for the submm upper limits, the majority of the IRAS excess targets show evidence for contaminating sources, based on investigation of higher resolution WISE and archival Spitzer and Herschel images. Finally, the exploratory submm measurements of the WISE excess source suggest that the hot dust present around these targets is not matched by a comparable population of colder, outer dust. More extensive and more sensitive Herschel observations of WISE excess sources will build upon this initial example to further define the characteristics of warm debris disks sources.

We have measured the multiplicity fractions and separation distributions of seven young star-forming regions using a uniform sample of young binaries. Both the multiplicity fractions and separation distributions are similar in the different regions. A tentative decline in the multiplicity fraction with increasing stellar density is apparent, even for binary systems with separations too close (19-100au) to have been dynamically processed. The separation distributions in the different regions are statistically indistinguishable over most separation ranges, and the regions with higher densities do not exhibit a lower proportion of wide (300-620au) relative to close (62-300au) binaries as might be expected from the preferential destruction of wider pairs. Only the closest (19-100au) separation range, which would be unaffected by dynamical processing, shows a possible difference in separation distributions between different regions. The combined set of young binaries, however, shows a distinct difference when compared to field binaries, with a significant excess of close (19-100au) systems among the younger binaries. Based on both the similarities and differences between individual regions, and between all seven young regions and the field, especially over separation ranges too close to be modified by dynamical processing, we conclude that multiple star formation is not universal and, by extension, the star formation process is not universal.

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.

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

SPICES (Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems) is a five-year M-class mission proposed to ESA Cosmic Vision. Its purpose is to image and characterize long-period extrasolar planets and circumstellar disks in the visible (450 - 900 nm) at a spectral resolution of about 40 using both spectroscopy and polarimetry. By 2020/22, present and near-term instruments will have found several tens of planets that SPICES will be able to observe and study in detail. Equipped with a 1.5 m telescope, SPICES can preferentially access exoplanets located at several AUs (0.5-10 AU) from nearby stars ($<$25 pc) with masses ranging from a few Jupiter masses to Super Earths ($\sim$2 Earth radii, $\sim$10 M$_{\oplus}$) as well as circumstellar disks as faint as a few times the zodiacal light in the Solar System.

As a part of our ongoing Volume-limited A-Star (VAST) adaptive optics survey, we have obtained observations of 26 binary systems with projected separations <100 AU, 13 of which have sufficient historical measurements to allow for refinement of their orbital elements. For each system with an estimated orbit, the dynamical system mass obtained was compared with the system mass estimated from mass-magnitude relations. Discrepancies between the dynamical and theoretical system mass can be explained by the presence of a previously unresolved spectroscopic component, or by a non-solar metallicity of the system. Using this approach to infer the presence of additional companions, a lower limit to the fraction of binaries, triples, and quadruples can be estimated as 39, 46, and 15 per cent, for systems with at least one companion within 100 AU. The fraction of multiple systems with three or more components shows a relative increase compared to the fraction for Solar-type primaries resolved in previous volume-limited surveys. The observations have also revealed a pair of potentially young ($<$100 Myr) M-dwarf companions, which would make an ideal benchmark for the theoretical models during the pre-Main Sequence contraction phase for M-dwarfs. In addition to those systems with orbit fits, we report 13 systems for which further orbital monitoring observations are required, 11 of which are newly resolved as a part of the VAST survey.

With a uniform VLT SINFONI data set of nine targets, we have developed an empirical grid of J,H,K spectra of the atmospheres of objects estimated to have very low substellar masses of \sim5-20 MJup and young ages of \sim1-50 Myr. Most of the targets are companions, objects which are especially valuable for comparison with atmosphere and evolutionary models, as they present rare cases in which the age is accurately known from the primary. Based on the sample youth, all objects are expected to have low surface gravity, and this study investigates the critical early phases of the evolution of substellar objects. The spectra are compared with grids of five different theoretical atmosphere models. This analysis represents the first systematic model comparison with infrared spectra of young brown dwarfs. The fits to the full JHK spectra of each object result in a range of best fit effective temperatures of +/-150-300K whether or not the full model grid or a subset restricted to lower log(g) values is used. This effective temperature range is significantly larger than the uncertainty typically assigned when using a single model grid. Fits to a single wavelength band can vary by up to 1000K using the different models. Since the overall shape of these spectra is governed more by the temperature than surface gravity, unconstrained model fits did not find matches with low surface gravity or a trend in log(g) with age. This suggests that empirical comparison with spectra of unambiguously young objects targets (such as these SINFONI data) may be the most reliable method to search for indications of low surface gravity and youth. For two targets, the SINFONI data are a second epoch and the data show no variations in morphology over time. The analysis of two other targets, AB Pic B and CT Cha B, suggests that these objects may have lower temperatures, and consequently lower masses, than previously estimated.

We have collated multiplicity data for five clusters (Taurus, Chamaeleon I, Ophiuchus, IC348, and the Orion Nebula Cluster). We have applied the same mass ratio (flux ratios of delta K <= 2.5) and primary mass cuts (~0.1-3.0 Msun) to each cluster and therefore have directly comparable binary statistics for all five clusters in the separation range 62-620 au, and for Taurus, Chamaeleon I, and Ophiuchus in the range 18-830 au. We find that the trend of decreasing binary fraction with cluster density is solely due to the high binary fraction of Taurus, the other clusters show no obvious trend over a factor of nearly 20 in density. With N-body simulations we attempt to find a set of initial conditions that are able to reproduce the density, morphology and binary fractions of all five clusters. Only an initially clumpy (fractal) distribution with an initial total binary fraction of 73 per cent (17 per cent in the range 62-620 au) is able to reproduce all of the observations (albeit not very satisfactorily). Therefore, if star formation is universal the initial conditions must be clumpy and with a high (but not 100 per cent) binary fraction. This could suggest that most stars, including M-dwarfs, form in binaries.

Dynamical interactions between planets and debris disks may sculpt the disk structure and impact planetary orbits, but only a few systems with both imaged planets and spatially resolved debris disks are known. With the Caltech Submm Observatory (CSO), we have observed the HR 8799 debris disk at 350{\mu}m. The 350{\mu}m map is the first spatially resolved measurement of the debris disk encircling the HR 8799 planetary system at this wavelength. Both the flux and size of the emission are consistent with a Kuiper belt of dust extending from ~100-300 AU. Although the resolution of the current map is limited, the map shows an indication of offset asymmetric emission, and several scenarios for this possibility are explored with radiative transfer calculations of a star-disk system and N-body numerical simulations of planet-disk interactions with parameters representative of the HR 8799 system.

With an adaptive optics imaging survey of 148 B6-A7 stars, we have tested the hypothesis that unresolved lower-mass companions are the source of the unexpected X-ray detections of stars in this spectral type range. The sample is composed of 63 stars detected in X-rays within the ROSAT All-Sky Survey and 85 stars that form a control sample; both subsets have the same restricted distribution of spectral type, age, X-ray sensitivity and separation coverage. A total of 68 companion candidates are resolved with separations ranging from 0.3" to 26.2", with 23 new detections. The multiple star frequency of the X-ray sample based on companions resolved within the ROSAT error ellipse is found to be 43 (+6,-6)%. The corresponding control sample multiple star frequency is three times lower at 12 (+4,-3)% -- a difference of 31\pm7%. These results are presented in the first of a series of papers based on our Volume-limited A-Star (VAST) survey -- a comprehensive study of the multiplicity of A-type stars.

We report the CO(J=1-0) observations of M51 using both the Combined Array for Research in Millimeter Astronomy (CARMA) and the Nobeyama 45m telescope (NRO45). We describe a procedure for the combination of interferometer and single-dish data. In particular, we discuss (1) the joint imaging and deconvolution of heterogeneous data, (2) the weighting scheme based on the root-mean-square (RMS) noise in the maps, (3) the sensitivity and uv-coverage requirements, and (4) the flux recovery of a combined map. We generate visibilities from the single-dish map and calculate the noise of each visibility based on the RMS noise. Our weighting scheme, though it is applied to discrete visibilities in this paper, is applicable to grids in uv-space, and this scheme may advance in future software development. For a realistic amount of observing time, the sensitivities of the NRO45 and CARMA visibility data sets are best matched by using the single dish baselines only up to 4-6 kilo-lambda (about 1/4-1/3 of the dish diameter). The synthesized beam size is determined to conserve the flux between synthesized beam and convolution beam. The superior uv-coverage provided by the combination of CARMA long baseline data with 15 antennas and NRO45 short spacing data results in the high image fidelity, which is evidenced by the excellent overlap between even the faint CO emission and dust lanes in an optical HST image and PAH emission in an Spitzer 8 micron image.

Direct-imaging searches for planets reveal wide orbit planets amenable to spectroscopy, and their atmospheres represent an important comparison to the irradiated atmospheres of Hot Jupiters. Using AO integral field spectroscopy of 2M1207 b, the shape of the continuum emission over the J, H, and K bands from the atmosphere of this young, planetary mass companion is measured in order to compare with atmospheric and evolutionary models, and objects of similar temperature in young clusters and the field. The 2M1207 b spectrum has the highest spectral resolution (R~300-1500) and largest wavelength coverage, including the first J-band spectrum, for this benchmark object. The high signal-to-noise of the data allow a clear identification of signatures of low surface gravity, and comparison with a grid of AMES-Dusty models reveals a best-fit effective temperature of Teff=1600 K with a preferred surface gravity of log g=4.5. The J-band flux is depressed relative to nearly all L-type objects, and the detailed shape of the absorption features across the H-band exhibit differences from the model predictions. The possible origins of 2M1207 b and its low luminosity are examined with the new data and analysis which suggest that extinction from a disk with large grains is a viable scenario and is preferred over scatttering off an optically thick disk. The 2M1207 b spectrum presents an important comparison for the types of features which may be present in upcoming spectra of the atmospheres of planets imaged in orbit around stellar primaries.

We report the highest-fidelity observations of the spiral galaxy M51 in CO emission, revealing the evolution of giant molecular clouds (GMCs) vis-a-vis the large-scale galactic structure and dynamics. The most massive GMCs (so-called GMAs) are first assembled and then broken up as the gas flow through the spiral arms. The GMAs and their H2 molecules are not fully dissociated into atomic gas as predicted in stellar feedback scenarios, but are fragmented into smaller GMCs upon leaving the spiral arms. The remnants of GMAs are detected as the chains of GMCs that emerge from the spiral arms into interarm regions. The kinematic shear within the spiral arms is sufficient to unbind the GMAs against self-gravity. We conclude that the evolution of GMCs is driven by large-scale galactic dynamics --their coagulation into GMAs is due to spiral arm streaming motions upon entering the arms, followed by fragmentation due to shear as they leave the arms on the downstream side. In M51, the majority of the gas remains molecular from arm entry through the inter-arm region and into the next spiral arm passage.

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

Four Ophiuchus binaries, two Class I systems and two Class II systems, with separations of ~450-1100 AU, were observed with the Owens Valley Radio Observatory (OVRO) millimeter interferometer. In each system, the 3 mm continuum maps show dust emission at the location of the primary star, but no emission at the position of the secondary. This result is different from observations of less evolved Class 0 binaries, in which dust emission is detected from both sources. The nondetection of secondary disks is, however, similar to the dust distribution seen in wide Class II Taurus binaries. The combined OVRO results from the Ophiuchus and Taurus binaries suggest that secondary disk masses are significantly lower than primary disk masses by the Class II stage, with initial evidence that massive secondary disks are reduced by the Class I stage. Although some of the secondaries retain hot inner disk material, the early dissipation of massive outer disks may negatively impact planet formation around secondary stars. Masses for the circumprimary disks are within the range of masses measured for disks around single T Tauri stars and, in some cases, larger than the minimum mass solar nebula. More massive primary disks are predicted by several formation models and are broadly consistent with the observations. Combining the 3 mm data with previous 1.3 mm observations, the dust opacity power-law index for each primary disk is estimated. The opacity index values are all less than the scaling for interstellar dust, possibly indicating grain growth within the circumprimary disks.

With the Navy Prototype Optical Interferometer (NPOI), the binary system Theta 1 Orionis C, the most massive member of the Trapezium, was spatially resolved over a time period extending from February 2006 to March 2007. The data show significant orbital motion over the 14 months, and, after combining the NPOI data with previous measurements of the system from the literature, the observations span 10 years of the orbit. Our results indicate that the secondary did not experience an unusually close periastron passage this year, in contradiction to the prediction of a recently published, highly eccentric ~11 year orbit. Future observations of this source will be required to improve the orbital solution. Possible implications of the results in terms of system distance are discussed, although a main conclusion of this work is that a definitive orbit solution will require more time to obtain sufficient phase coverage, and that the interaction effects expected at periastron did not occur in 2007.

Close pre-main-sequence binary stars are expected to clear central holes in their protoplanetary disks, but the extent to which material can flow from the circumbinary disk across the gap onto the individual circumstellar disks has been unclear. In binaries with eccentric orbits, periodic perturbation of the outer disk is predicted to induce mass flow across the gap, resulting in accretion that varies with the binary period. This accretion may manifest itself observationally as periodic changes in luminosity. Here we present a search for such periodic accretion in the pre-main-sequence spectroscopic binary UZ Tau E. We present BVRI photometry spanning three years; we find that the brightness of UZ Tau E is clearly periodic, with a best-fit period of 19.16 +/- 0.04 days. This is consistent with the spectroscopic binary period of 19.13 days, refined here from analysis of new and existing radial velocity data. The brightness of UZ Tau E shows significant random variability, but the overall periodic pattern is a broad peak in enhanced brightness, spanning more than half the binary orbital period. The variability of the H-alpha line is not as clearly periodic, but given the sparseness of the data, some periodic component is not ruled out. The photometric variations are in good agreement with predictions from simulations of binaries with orbital parameters similar to those of UZ Tau E, suggesting that periodic accretion does occur from circumbinary disks, replenishing the inner disks and possibly extending the timescale over which they might form planets.

It is now well established that a large fraction of the low-mass stars are binaries or higher order multiples. Similarly a sizable fraction have giant planets. In contrast to these, the situation for brown dwarf companions is complicated: While close systems seem to be extremely rare, wide systems are possibly more common. In this paper, we present new results on a survey for low-mass companions in the Hyades. After measuring precisely the radial velocity of 98 Hyades dwarf stars for 5 years, we have selected all stars that show low-amplitude long-period trends. With AO-observations of these 14 stars we found companion candidates around nine of them, where one star has two companions. The two companions of HIP 16548 have masses between 0.07 to 0.08 Mo, and are thus either brown dwarfs or very low mass stars. In the case of HAN 172 we found a companion with a mass between 0.08 to 0.10 Mo, which is again between a star and a brown dwarf. The other seven stars all have stellar companions. In two additional cases, the RV-variations are presumably caused by stellar activity, and in another case the companion could be a short-period binary. The images of the remaining two stars are slightly elongated, which might imply that even these are binaries. Because at least 12 of the 14 stars showing low-amplitude RV trends turn out to have companions with a mass greater than 70 MJupiter, or are just active, we finally estimate the number of companions with masses between 10 MJupiter and 70 MJupiter within 8 AU of the host stars in the Hyades as less equal 2%.

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.

A combination of high-resolution and wide-field imaging reveals two binary stars and one triple star system among the sample of the first 11 stars with planets detected by radial velocity variations. High resolution speckle or adaptive optics (AO) data probe subarcsecond scales down to the diffraction limit of the Keck 10-m or Lick 3-m, and direct images or AO images are sensitive to a wider field, extending to 10" or 38", depending upon the camera. One of the binary systems -- HD 114762 -- was not previously known to be a spatially resolved multiple system; additional data taken with the combination of Keck adaptive optics and NIRSPEC are used to characterize the new companion. The second binary system -- Tau Boo -- was a known multiple with two conflicting orbital solutions; the current data will help constrain the discrepant estimates of periastron time and separation. Another target -- 16 Cyg B -- was a known common proper motion binary, but the current data resolve a new third component, close to the wide companion 16 Cyg A. Both the HD 114762 and 16 Cyg B systems harbor planets in eccentric orbits, while the Tau Boo binary contains an extremely close planet in a tidally-circularized orbit. Although the sample is currently small, the proportion of binary systems is comparable to that measured in the field over a similar separation range. Incorporating the null result from another companion search project lowers the overall fraction of planets in binary systems, but the detections in our survey reveal that planets can form in binaries separated by less than 50 AU.

Two hundred and forty-two members of the Praesepe and alpha Persei clusters ranging from B (~5Msun) to early-M (~0.5 Msun) have been surveyed with high angular resolution imaging. The 39 binary and 1 quadruple systems detected encompass separations from 0".053 to 7".28; 28 of the systems are new detections and there are 9 candidate substellar companions. The main results from the survey are: * Over the projected separation range of 26-581 AU and magnitude differences of DeltaK < 4.0 mag, the companion star fraction (CSF) for alpha Persei is 0.09 +/- 0.03 and for Praesepe is 0.10 +/- 0.03. This fraction is consistent with the field G-dwarf value, implying that there is not a systematic decline in multiplicity with age on timescales of at 4 +1/-1.5 AU, significantly smaller value than both the field G-dwarf and the nearby T Tauri distributions. A simple population synthesis model suggests that the G-dwarf binary population is a combination of ~30% dark cloud and ~70% giant molecular cloud systems. * An exploration of the binary star properties reveals a cluster CSF that increases with decreasing target mass and a cluster mass ratio distribution that rises more sharply for higher mass stars, but is independent of binary separation. These observational trends are consistent with several models of capture in small clusters and simulations of accretion following fragmentation in a cluster environment. * Among the cluster A stars, there is a higher fraction of binaries in the subset with X-ray detections, consistent with the hypothesis that lower mass companions are the true source of X-ray emission. * In alpha Persei, the rotational velocities for solar-type binaries with separations less that 60 AU the rotational evolution of young stars.

We have undertaken an adaptive optics (AO) imaging survey of extrasolar planetary systems and stars showing interesting radial velocity trends from high precision radial velocity searches. AO increases the resolution and dynamic range of an image, substantially improving the detectability of faint close companions. This survey is sensitive to objects less luminous than the bottom of the main sequence at separations as close as 1 arcsec. We have detected stellar companions to the planet bearing stars HD 114762 and Tau Boo. We have also detected a companion to the non-planet bearing star 16 Cyg A.

We report the discovery of five T Tauri star systems, two of which are resolved binaries, in the vicinity of the nearest known region of recent star formation, the TW Hydrae Association. The newly discovered systems display the same signatures of youth (namely high X-ray flux, large Li abundance and strong chromospheric activity) and the same proper motion as the original five members. These similarities firmly establish the group as a bona fide T Tauri association, unique in its proximity to Earth and its complete isolation from any known molecular clouds. At an age of ~10 Myr and a distance of ~50 pc, the association members are excellent candidates for future studies of circumstellar disk dissipation and the formation of brown dwarfs and planets. Indeed, as an example, our speckle imaging revealed a faint, very likely companion 2" north of CoD-33 7795 (TWA 5). Its color and brightness suggest a spectral type ~M8.5 which, at an age of ~10^7 years, implies a mass ~20 M(Jupiter).

A 2.2 micron speckle imaging survey of 167 bright (K $< 8.5$ mag) Hyades members reveals a total of 33 binaries with separations spanning 0."044 to 1."34 and magnitude differences as large as 5.5 mag. Of these binaries, 13 are new detections. Combining the results of this survey with previous radial velocity, optical speckle, and direct imaging Hyades surveys, the detected multiplicity of the sample is: 98 singles, 59 binaries, and 10 triples. A statistical analysis of this sample investigates a variety of multiple star formation and evolution theories. Accounting for the inability to detect high flux ratio binaries results in an implied companion star fraction (csf) of 0.30 +/- 0.06 in this separation range 0."1 to 1."07 (5 to 50 AU). The Hyades csf is intermediate between the values derived from observations of T Tauri stars (T Tauri csf=0.40 +/- 0.08) and solar neighborhood G-dwarfs (G-dwarf csf=0.14 +/- 0.03). This result allows for an evolution of the csf from an initially high value for the pre-main sequence to that found for main sequence stars. Within the Hyades, the csf is independent of the radial distance from the cluster center and the primary star mass. The distribution of mass ratios is best fit by a power law q^(-1.3 +/- 0.3) and shows no dependence on the primary mass, binary separation, or the radial distance from the cluster center. Overall, the Hyades data are consistent with scale-free fragmentation, but inconsistent with capture in small clusters and disk-assisted capture in small clusters. Without testable predictions, scale-dependent fragmentation and disk fragmentation cannot be assessed with the Hyades data.