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

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

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

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

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

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

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

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.

Starlight that becomes linearly polarized by magnetically aligned dust grains provides a viable diagnostic of the interstellar magnetic field (ISMF). A survey is underway to map the local ISMF using data collected at eight observatories in both hemispheres. Two approaches are used to obtain the magnetic structure: statistically evaluating magnetic field directions traced by multiple polarization position angles, and least-squares fits that provide the dipole component of the magnetic field. We find that the magnetic field in the circumheliospheric interstellar medium (CHM), which drives winds of interstellar gas and dust through the heliosphere, drapes over the heliopause and influences polarization measurements. We discover a polarization band that can be described with a great circle that traverses the heliosphere nose and ecliptic poles. A gap in the band appears in a region coinciding both with the highest heliosheath pressure, found by IBEX, and the center of the Loop I superbubble. The least-squares analysis finds a magnetic dipole component of the polarization band with the axis oriented toward the ecliptic poles. The filament of dust around the heliosphere and the warm helium breeze flowing through the heliosphere trace the same magnetic field directions. Regions along the polarization band near the heliosphere nose have magnetic field orientations within 15 degrees of sightlines. Regions in the IBEX ribbon have field directions within 40 degrees of the plane of the sky. Several spatially coherent magnetic filaments are within 15 pc. Most of the low frequency radio emissions detected by the two Voyager spacecraft follow the polarization band. The geometry of the polarization band is compared to the Local Interstellar Cloud, the Cetus Ripple, the BICEP2 low opacity region, Ice Cube IC59 galactic cosmic ray data, and Cassini results.

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

We present H and K band imaging polarimetry for the PDS 66 circumstellar disk obtained during the commissioning of the Gemini Planet Imager (GPI). Polarization images reveal a clear detection of the disk in to the 0.12'' inner working angle (IWA) in H band, almost 3 times as close to the star as the previous HST observations with NICMOS and STIS (0.35'' effective IWA). The centro-symmetric polarization vectors confirm that the bright inner disk detection is due to circumstellar scattered light. A more diffuse disk extends to a bright outer ring centered at 80 AU. We discuss several physical mechanisms capable of producing the observed ring + gap structure. GPI data confirm enhanced scattering on the East side of the disk which is inferred to be nearer to us. We also detect a lateral asymmetry in the South possibly due to shadowing from material within the inner working angle. This likely corresponds to a temporally variable azimuthal asymmetry observed in HST/STIS coronagraphic imaging.

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

We present a multiple scattering vector radiative transfer model which produces disk integrated, full phase polarized light curves for reflected light from an exoplanetary atmosphere. We validate our model against results from published analytical and computational models and discuss a small number of cases relevant to the existing and possible near-future observations of the exoplanet HD 189733b. HD 189733b is arguably the most well observed exoplanet to date and the only exoplanet to be observed in polarized light, yet it is debated if the planet's atmosphere is cloudy or clear. We model reflected light from clear atmospheres with Rayleigh scattering, and cloudy or hazy atmospheres with Mie and fractal aggregate particles. We show that clear and cloudy atmospheres have large differences in polarized light as compared to simple flux measurements, though existing observations are insufficient to make this distinction. Futhermore, we show that atmospheres that are spatially inhomogeneous, such as being partially covered by clouds or hazes, exhibit larger contrasts in polarized light when compared to clear atmospheres. This effect can potentially be used to identify patchy clouds in exoplanets. Given a set of full phase polarimetric measurements, this model can constrain the geometric albedo, properties of scattering particles in the atmosphere and the longitude of the ascending node of the orbit. The model is used to interpret new polarimetric observations of HD 189733b in a companion paper.

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.

The interstellar magnetic field (ISMF) near the heliosphere is a basic part of the solar neighborhood that can only be studied using polarized starlight. Results of an ongoing survey of polarized starlight are analyzed with the goal of linking the interstellar magnetic field that shapes the heliosphere to the nearby field in interstellar space. New results for the direction of the nearby ISMF, based on a merit function that utilizes polarization position angles, identify several magnetic components. The dominant interstellar field, B_pol, is aligned with the direction L,B= 36.2,49.0 (+/-16.0) degrees and is within 8 degrees of the IBEX Ribbon ISMF direction. Stars tracing B_pol have the same mean distance as stars that do not trace B_pol, but show weaker polarizations consistent with lower column densities of polarizing grains. The variations in the polarization position angle directions indicate a low level of magnetic turbulence. B_pol is found after excluding polarizations that trace a separate magnetic structure that apparently is due to interstellar dust deflected around the heliosphere. Local interstellar cloud velocities relative to the LSR increase with the angles between the LSR velocities and ISMF, indicating that the kinematics of local interstellar material is ordered by the ISMF. Polarization and color excess data are consistent with an extension of Loop I to the solar vicinity. Polarizations are consistent with previous findings of more efficient grain alignment in low column density sightlines. Optical polarization and color excess data indicate the presence of nearby interstellar dust in the BICEP2 field. Color excess E(B-V) indicates an optical extinction of A_V about 0.59 mag in the BICEP2 field, while the polarization data indicate that A_V is larger than 0.09 mag. The IBEX Ribbon ISMF extends to the boundaries of the BICEP2 region.

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.

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

We present 50 nights of polarimetric observations of HD 189733 in $B$ band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of order 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a $3 \sigma$ upper limit to the planet-induced polarization of the system of 58 ppm in $B$ band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope (Wiktorowicz 2009). A physically-motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a $3 \sigma$ upper limit to the geometric albedo of HD 189733b of $A_g < 0.37$. This value is consistent with the value $A_g = 0.226 \pm 0.091$ derived from occultation observations with HST STIS (Evans et al. 2013), but it is inconsistent with the large $A_g = 0.61 \pm 0.12$ albedo reported by (Berdyugina et al. 2011).

A recently discovered filament of polarized starlight that traces a coherent magnetic field is shown to have several properties that are consistent with an origin in the outer heliosheath of the heliosphere: (1) The magnetic field that provides the best fit to the polarization position angles is directed within 6.7+-11 degrees of the observed upwind direction of the flow of interstellar neutral helium gas through the heliosphere. (2) The magnetic field is ordered; the component of the variation of the polarization position angles that can be attributed to magnetic turbulence is small. (3) The axis of the elongated filament can be approximated by a line that defines an angle of 80+/-14 degrees with the plane that is formed by the interstellar magnetic field vector and the vector of the inflowing neutral gas (the "BV" plane). We propose that this polarization feature arises from aligned interstellar dust grains in the outer heliosheath where the interstellar plasma and magnetic field are deflected around the heliosphere. The proposed outer heliosheath location of the polarizing grains requires confirmation by modeling grain-propagation through three-dimensional MHD heliosphere models that simultaneously calculate torques on asymmetric dust grains interacting with the heliosphere.

From a single, 3.8-hour observation of asteroid (4) Vesta at $13.7^\circ$ phase angle with the POLISH2 polarimeter at the Lick Observatory Shane 3-m telescope, we confirm rotational modulation of linear polarization in $B$ and $V$ bands. We measure the peak-to-peak modulation in degree of linear polarization to be $\Delta P = (294 \pm 35) \times 10^{-6}$ (ppm) and time-averaged $\Delta P / P = 0.0577 \pm 0.0069$. After rotating the plane of linear polarization to the scattering plane, asteroidal rotational modulation is detected with $12 \sigma$ confidence and observed solely in Stokes $Q/I$. POLISH2 simultaneously measures Stokes $I$, $Q$, $U$ (linear polarization), and $V$ (circular polarization), but we detect no significant circular polarization with a $1 \sigma$ upper limit of 140 ppm in $B$ band. Circular polarization is expected to arise from multiple scattering of sunlight by rough surfaces, and it has previously been detected in nearly all other classes of Solar System bodies save asteroids. Subsequent observations may be compared with surface albedo maps from the Dawn Mission, which may allow identification of compositional variation across the asteroidal surface. These results demonstrate the high accuracy achieved by POLISH2 at the Lick 3-m telescope, which is designed to directly detect scattered light from spatially unresolvable exoplanets.

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

The local interstellar magnetic field affects both the heliosphere and the surrounding cluster of interstellar clouds (CLIC). Measurements of linearly polarized starlight provide the only test of the magnetic field threading the CLIC. Polarization measurements of the CLIC magnetic field show multiple local magnetic structures, one of which is aligned with the magnetic field traced by the center of the "ribbon" of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX). Comparisons between the bulk motion of the CLIC through the local standard of rest, the magnetic field direction, the geometric center of Loop I, and the polarized dust bridge extending from the heliosphere toward the North Polar Spur direction all suggest that the CLIC is part of the rim region of the Loop I superbubble.

In contrast to photometric transits, whose peak signal occurs at mid-transit due to occultation of the brightest region of the disk, polarimetric transits provide a signal upon ingress and egress due to occultation of the polarized stellar limb. Limb polarization, the bright corollary to limb darkening, arises from the $90^\circ$ scattering angle and low optical depth experienced by photons at the limb. In addition to the ratio $R_{\rm p} / R_*$, the amplitude of a polarimetric transit is expected to be controlled by the strength and width of the stellar limb polarization profile, which depend on the scattering-to-total opacity ratio at the stellar limb. We present a short list of the systems providing the highest expected signal-to-noise ratio for detection of this effect, and we draw particular attention to HD 80606b. This planet is spin/orbit misaligned, has a three-hour ingress, and has a bright parent star, which make it an attractive target. We report on test observations of an HD 80606b ingress with the POLISH2 polarimeter at the Lick Observatory Shane 3-m telescope. We conclude that unmodeled telescope systematic effects prevented polarimetric detection of this event. We outline a roadmap for further refinements of exoplanet polarimetry, whose eventual success will require a further factor of ten reduction in systematic noise.

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

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

We present on-sky polarimetric observations with the Gemini Planet Imager (GPI) obtained at straight Cassegrain focus on the Gemini South 8-m telescope. Observations of polarimetric calibrator stars, ranging from nearly unpolarized to strongly polarized, enable determination of the combined telescope and instrumental polarization. We find the conversion of Stokes $I$ to linear and circular instrumental polarization in the instrument frame to be $I \rightarrow (Q_{\rm IP}, U_{\rm IP}, P_{\rm IP}, V_{\rm IP}) = (-0.037 \pm 0.010\%, +0.4338 \pm 0.0075\%, 0.4354 \pm 0.0075\%, -6.64 \pm 0.56\%)$. Such precise measurement of instrumental polarization enables $\sim 0.1\%$ absolute accuracy in measurements of linear polarization, which together with GPI's high contrast will allow GPI to explore scattered light from circumstellar disk in unprecedented detail, conduct observations of a range of other astronomical bodies, and potentially even study polarized thermal emission from young exoplanets. Observations of unpolarized standard stars also let us quantify how well GPI's differential polarimetry mode can suppress the stellar PSF halo. We show that GPI polarimetry achieves cancellation of unpolarized starlight by factors of 100-200, reaching the photon noise limit for sensitivity to circumstellar scattered light for all but the smallest separations at which the calibration for instrumental polarization currently sets the limit.

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.

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 key indicator of the galactic environment of the Sun is provided by the magnetic field in the interstellar medium (ISM), which influences the shape of the heliosphere. The direction of the nearby interstellar magnetic field (ISMF) is determined from starlight polarized in the ISM. The local ISMF direction is found from the ISMF direction that provides the best fit to the polarization position angles of nearby stars, using weighted fits to the data. New polarization observations are included in the analysis. The best-fitting ISMF is close to the magnetic field direction traced by the center of the Ribbon of energetic neutral atoms, discovered by the Interstellar Boundary Explorer spacecraft. Both the magnetic field and kinematics of the local ISM are consistent with a scenario where the local ISM is a fragment of the Loop I superbubble. An ordered component of the local ISMF is found in a region where PlanetPol data show that polarization increases with distance. It extends to within 8 parsecs of the Sun and implies a weak curvature in the nearby ISMF. Variations from the ordered component indicate turbulence of +/-23 deg. The local ISMF is generally uniform in direction over spatial scales of 8-200 parsecs so that it appears similar to interarm magnetic fields. The best-fitting ISMF direction also agrees with the position of tail-in spatial asymmetries in GeV-TeV galactic cosmic rays. The peculiar geometrical relation between the CMB dipole moment, the heliosphere nose, and local ISMF is supported by these new results. Radiative torques are not likely to play a role in grain alignment for these polarizations.

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.

The Gemini Planet Imager (GPI) is a high performance adaptive optics system being designed and built for the Gemini Observatory. GPI is optimized for high contrast imaging, combining precise and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with integral field and polarimetry capabilities. The primary science goal for GPI is the direct detection and characterization of young, Jovian-mass exoplanets. For plausible assumptions about the distribution of gas giant properties at large semi-major axes, GPI will be capable of detecting more than 10% of gas giants more massive than 0.5 M_J around stars younger than 100 Myr and nearer than 75 parsecs. For systems younger than 1 Gyr, gas giants more massive than 8 M_J and with semi-major axes greater than 15 AU are detected with completeness greater than 50%. A survey targeting young stars in the solar neighborhood will help determine the formation mechanism of gas giant planets by studying them at ages where planet brightness depends upon formation mechanism. Such a survey will also be sensitive to planets at semi-major axes comparable to the gas giants in our own solar system. In the simple, and idealized, situation in which planets formed by either the "hot-start" model of Burrows et al. (2003) or the core accretion model of Marley et al. (2007), a few tens of detected planets are sufficient to distinguish how planets form.

Differential astrometry measurements from the Palomar High-precision Astrometric Search for Exoplanet Systems have been combined with lower precision single-aperture measurements covering a much longer timespan (from eyepiece measurements, speckle interferometry, and adaptive optics) to determine improved visual orbits for 20 binary stars. In some cases, radial velocity observations exist to constrain the full three-dimensional orbit and determine component masses. The visual orbit of one of these binaries---alpha Com (HD 114378)---shows that the system is likely to have eclipses, despite its very long period of 26 years. The next eclipse is predicted to be within a week of 2015 January 24.

The Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES) monitored 51 sub-arcsecond binary systems to determine precision binary orbits, study the geometries of triple and quadruple star systems, and discover previously unknown faint astrometric companions as small as giant planets. PHASES measurements made with the Palomar Testbed Interferometer (PTI) from 2002 until PTI ceased normal operations in late 2008 are presented. Infrared differential photometry of several PHASES targets were measured with Keck Adaptive Optics and are presented.

Variations in the spatial configuration of the interstellar magnetic field (ISMF) near the Sun can be constrained by comparing the ISMF direction at the heliosphere found from the Interstellar Boundary Explorer spacecraft (IBEX) observations of a 'Ribbon' of energetic neutral atoms (ENAs), with the ISMF direction derived from optical polarization data for stars within ~40 pc. Using interstellar polarization observations towards ~30 nearby stars within 90 deg of the heliosphere nose, we find that the best fits to the polarization position angles are obtained for a magnetic pole directed towards ecliptic coordinates of lambda, beta 263 deg, 37 deg (or galactic coordinates of L,B 38 deg, 23deg), with uncertainties of +/- 35 deg, based on the broad minimum of the best fits and the range of data quality. This magnetic pole is 33 deg from the magnetic pole that is defined by the center of the arc of the ENA Ribbon. The IBEX ENA ribbon is seen in sightlines that are perpendicular to the ISMF as it drapes over the heliosphere. The similarity of the polarization and Ribbon directions for the local ISMF suggest that the local field is coherent over scale sizes of tens of parsecs. The ISMF vector direction is nearly perpendicular to the flow of local ISM through the local standard of rest, supporting a possible local ISM origin related to an evolved expanding magnetized shell. The local ISMF direction is found to have a curious geometry with respect to the cosmic microwave background dipole moment.

Using the POLISH instrument, I am unable to reproduce the large-amplitude polarimetric observations of Berdyugina et al. (2008) to the >99.99% confidence level. I observe no significant polarimetric variability in the HD 189733 system, and the upper limit to variability from the exoplanet is Delta_P < 7.9 x 10^(-5) with 99% confidence in the 400 nm to 675 nm wavelength range. Berdyugina et al. (2008) report polarized, scattered light from the atmosphere of the HD 189733b hot Jupiter with an amplitude of two parts in 10^4. Such a large amplitude is over an order of magnitude larger than expected given a geometric albedo similar to other hot Jupiters. However, my non-detection of polarimetric variability phase-locked to the orbital period of the exoplanet, and the lack of any significant variability, shows that the polarimetric modulation reported by Berdyugina et al. (2008) cannot be due to the exoplanet.

We present commissioning data for the POLISH instrument obtained on the Hale 5-m telescope. The goal of this high precision polarimeter is to constrain orbital inclination of high mass X-ray binaries and to therefore obtain independent mass estimates for their black hole companions. We have obtained photon shot noise limited precision on standard stars, and we have measured the polarization of bright stars at the part per million level on a nightly basis. Systematic effects have been reduced to less than 1% of the measured polarization for polarized sources and to the part per million level for weakly polarized sources. The high sensitivity of this instrument to asymmetry suggests that valuable contributions will be made in many other fields, including studies of extrasolar planets, debris disks, and stellar astrophysics.

mu Orionis was identified by spectroscopic studies as a quadruple star system. Seventeen high precision differential astrometry measurements of mu Ori have been collected by the Palomar High-precision Astrometric Search for Exoplanet Systems (PHASES). These show both the motion of the long period binary orbit and short period perturbations superimposed on that caused by each of the components in the long period system being themselves binaries. The new measurements enable the orientations of the long period binary and short period subsystems to be determined. Recent theoretical work predicts the distribution of relative inclinations between inner and outer orbits of hierarchical systems to peak near 40 and 140 degrees. The degree of coplanarity of this complex system is determined, and the angle between the planes of the A-B and Aa-Ab orbits is found to be 136.7 +/- 8.3 degrees, near the predicted distribution peak at 140 degrees; this result is discussed in the context of the handful of systems with established mutual inclinations. The system distance and masses for each component are obtained from a combined fit of the PHASES astrometry and archival radial velocity observations. The component masses have relative precisions of 5% (component Aa), 15% (Ab), and 1.4% (each of Ba and Bb). The median size of the minor axes of the uncertainty ellipses for the new measurements is 20 micro-arcseconds. Updated orbits for delta Equulei, kappa Pegasi, and V819 Herculis are also presented.

Aptly named, ice giants such as Uranus and Neptune contain significant amounts of water. While this water cannot be present near the cloud tops, it must be abundant in the deep interior. We investigate the likelihood of a liquid water ocean existing in the hydrogen-rich region between the cloud tops and deep interior. Starting from an assumed temperature at a given upper tropospheric pressure (the photosphere), we follow a moist adiabat downward. The mixing ratio of water to hydrogen in the gas phase is small in the photosphere and increases with depth. The mixing ratio in the condensed phase is near unity in the photosphere and decreases with depth; this gives two possible outcomes. If at some pressure level the mixing ratio of water in the gas phase is equal to that in the deep interior, then that level is the cloud base. Alternately, if the mixing ratio of water in the condensed phase reaches that in the deep interior, then the surface of a liquid ocean will occur. We find that Neptune is both too warm (photospheric temperature too high) and too dry (mixing ratio of water in the deep interior too low) for liquid oceans to exist at present. To have a liquid ocean, Neptune's deep interior water to gas ratio would have to be higher than current models allow, and the density at 19 kbar would have to be ~ 0.8 g/cm^3. Such a high density is inconsistent with gravitational data obtained during the Voyager flyby. As Neptune cools, the probability of a liquid ocean increases. Extrasolar "hot Neptunes," which presumably migrate inward toward their parent stars, cannot harbor liquid water oceans unless they have lost almost all of the hydrogen and helium from their deep interiors.

The Stellar Planet Survey (STEPS) is an ongoing astrometric search for giant planets and brown dwarfs around a sample of ~30 M-dwarfs. We have discovered several low-mass companions by measuring the motion of our target stars relative to their reference frames. The highest mass discovery thus far is G 78-28B, a companion to the M-dwarf G 78-28A. The orbital period is 4.18 +/- 0.03 y, the system mass is 0.565 +/- 0.055 Msolar, and the semi-major axis is 2.19 +/- 0.10 AU. Imaging observations with the Keck laser guide star adaptive optics (LGSAO) and the Palomar AO instruments resolved the system and also yielded JHK-band delta magnitudes. We use the orbital solution, light ratios, and mass-luminosity relationships to derive component masses of MA = 0.370 +/- 0.034 Msolar and MB = 0.195 +/- 0.021 Msolar. G 78-28B is of type M4 V based upon its colors and mass. We also discovered GJ 231.1C, a companion to GJ 231.1B, with STEPS and imaged the companion with LGSAO and Palomar AO, but the orbital period is longer than our observing baseline; thus the system parameters are less constrained. In GJ 231.1BC the masses are MB = 0.25 +/- 0.06 Msolar and MC =0.12 +/- 0.02 Msolar. The inferred spectral type of GJ 231.1C is M5 V. We demonstrate the results of the current state of mass estimation techniques with our data.

We examine 14 hot white dwarfs for signs of infrared excess using the Spitzer Space Telescope. Twelve of the objects are massive white dwarfs which have been suggested to be the result of binary mergers. The remaining two objects are undermassive white dwarfs which again may be the result of mergers or the inspiral of a substellar companion. In no case do we find any evidence for significant infrared excesses out to wavelengths of 8 microns. This places strong constraints on the presence of orbiting dust and weaker constraints on the presence of close substellar companions.

kappa Pegasi is a well-known, nearby triple star system. It consists of a ``wide'' pair with semi-major axis 235 milli-arcseconds, one component of which is a single-line spectroscopic binary (semi-major axis 2.5 milli-arcseconds). Using high-precision differential astrometry and radial velocity observations, the masses for all three components are determined and the relative inclinations between the wide and narrow pairs' orbits is found to be 43.8 +/- 3.0 degrees, just over the threshold for the three body Kozai resonance. The system distance is determined to 34.60 +/- 0.21 parsec, and is consistent with trigonometric parallax measurements.