More than a dozen young stars host spiral arms in their surrounding protoplanetary disks. The excitation mechanisms of such arms are under debate. The two leading hypotheses -- companion-disk interaction and gravitational instability (GI) -- predict distinct motion for spirals. By imaging the MWC 758 spiral arm system at two epochs spanning ${\sim}5$ yr using the SPHERE instrument on the Very Large Telescope (VLT), we test the two hypotheses for the first time. We find that the pattern speeds of the spirals are not consistent with the GI origin. Our measurements further evince the existence of a faint "missing planet" driving the disk arms. The average spiral pattern speed is $0.\!^\circ22\pm0.\!^\circ03$ yr$^{-1}$, pointing to a driver at $172_{-14}^{+18}$ au around a $1.9$ $M_\odot$ central star if it is on a circular orbit. In addition, we witness time varying shadowing effects on a global scale that are likely originated from an inner disk.

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.

To characterize circumstellar systems in high contrast imaging, the fundamental step is to construct a best point spread function (PSF) template for the non-circumstellar signals (i.e., star light and speckles) and separate it from the observation. With existing PSF construction methods, the circumstellar signals (e.g., planets, circumstellar disks) are unavoidably altered by over-fitting and/or self-subtraction, making forward modeling a necessity to recover these signals. We present a forward modeling--free solution to these problems with data imputation using sequential non-negative matrix factorization (DI-sNMF). DI-sNMF first converts this signal separation problem to a "missing data" problem in statistics by flagging the regions which host circumstellar signals as missing data, then attributes PSF signals to these regions. We mathematically prove it to have negligible alteration to circumstellar signals when the imputation region is relatively small, which thus enables precise measurement for these circumstellar objects. We apply it to simulated point source and circumstellar disk observations to demonstrate its proper recovery of them. We apply it to Gemini Planet Imager (GPI) K1-band observations of the debris disk surrounding HR 4796A, finding a tentative trend that the dust is more forward scattering as the wavelength increases. We expect DI-sNMF to be applicable to other general scenarios where the separation of signals is needed.

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.

HiCAT is a high-contrast imaging testbed designed to provide complete solutions in wavefront sensing, control and starlight suppression with complex aperture telescopes. The pupil geometry of such observatories includes primary mirror segmentation, central obstruction, and spider vanes, which make the direct imaging of habitable worlds very challenging. The testbed alignment was completed in the summer of 2014, exceeding specifications with a total wavefront error of 12nm rms over a 18mm pupil. The installation of two deformable mirrors for wavefront control is to be completed in the winter of 2015. In this communication, we report on the first testbed results using a classical Lyot coronagraph. We also present the coronagraph design for HiCAT geometry, based on our recent development of Apodized Pupil Lyot Coronagraph (APLC) with shaped-pupil type optimizations. These new APLC-type solutions using two-dimensional shaped-pupil apodizer render the system quasi-insensitive to jitter and low-order aberrations, while improving the performance in terms of inner working angle, bandpass and contrast over a classical APLC.

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.

Instrumentation techniques in the field of direct imaging of exoplanets have greatly advanced over the last two decades. Two of the four NASA-commissioned large concept studies involve a high-contrast instrument for the imaging and spectral characterization of exo-Earths from space: LUVOIR and HabEx. This whitepaper describes the status of 8 optical testbeds in the US and France currently in operation to experimentally validate the necessary technologies to image exo-Earths from space. They explore two complementary axes of research: (i) coronagraph designs and manufacturing and (ii) active wavefront correction methods and technologies. Several instrument architectures are currently being analyzed in parallel to provide more degrees of freedom for designing the future coronagraphic instruments. The necessary level of performance has already been demonstrated in-laboratory for clear off-axis telescopes (HabEx-like) and important efforts are currently in development to reproduce this accomplishment on segmented and/or on-axis telescopes (LUVOIR-like) over the next two years.

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.

Segmented telescopes are a possibility to enable large-aperture space telescopes for the direct imaging and spectroscopy of habitable worlds. However, the complexity of their aperture geometry, due to the central obstruction, support structures and segment gaps, makes high-contrast imaging challenging. The High-contrast Imager for Complex Aperture Telescopes (HiCAT) testbed was designed to study and develop solutions for such telescope pupils using wavefront control and coronagraphic starlight suppression. The testbed design has the flexibility to enable studies with increasing complexity for telescope aperture geometries: off-axis telescopes, on-axis telescopes with central obstruction and support structures - e.g. the Wide Field Infrared Survey Telescope (WFIRST) - to on-axis segmented telescopes, including various concepts for a Large UV, Optical, IR telescope (LUVOIR). In the past year, HiCAT has made significant hardware and software updates to accelerate the development of the project. In addition to completely overhauling the software that runs the testbed, we have completed several hardware upgrades, including the second and third deformable mirror, and the first custom Apodized Pupil Lyot Coronagraph (APLC) optimized for the HiCAT aperture, which is similar to one of the possible geometries considered for LUVOIR. The testbed also includes several external metrology features for rapid replacement of parts, and in particular the ability to test multiple apodizers readily, an active tip-tilt control system to compensate for local vibration and air turbulence in the enclosure. On the software and operations side, the software infrastructure enables 24/7 automated experiments that include routine calibration tasks and high-contrast experiments. We present an overview and status update of the project, on the hardware and software side, and describe results obtained with APLC WFC.

High contrast imaging using coronagraphy is one of the main avenues to enable the search for life on extrasolar Earth analogs. The HiCAT testbed aims to demonstrate coronagraphy and wavefront control for segmented on-axis space telescopes as envisioned for a future large UV optical IR mission (LUVOIR). Our software infrastructure enables 24/7 automated operation of high-contrast imaging experiments while monitoring for safe operating parameters, along with graceful shutdown processes for unsafe conditions or unexpected errors. The infrastructure also includes a calibration suite that can run nightly to catch regressions and track optical performance changes over time, and a testbed simulator to support software development and testing, as well as optical modeling necessary for high-contrast algorithms. This paper presents a design and implementation of testbed control software to leverage continuous integration whether the testbed is available or not.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

The James Webb Space Telescope (JWST) Optical Simulation Testbed (JOST) is a hardware simulator for wavefront sensing and control designed to produce JWST-like images. A model of the JWST three mirror anas- tigmat is realized with three lenses in the form of a Cooke triplet, which provides JWST-like optical quality over a field equivalent to a NIRCam module. An Iris AO hexagonally segmented mirror stands in for the JWST primary. This setup successfully produces images extremely similar to expected JWST in-flight point spread functions (PSFs), and NIRCam images from cryotesting, in terms of the PSF morphology and sampling relative to the diffraction limit. The segmentation of the primary mirror into subapertures introduces complexity into wavefront sensing and control (WFS&C) of large space based telescopes like JWST. JOST provides a platform for independent analysis of WFS&C scenarios for both commissioning and maintenance activities on such ob- servatories. We present an update of the current status of the testbed including both single field and wide-field alignment results. We assess the optical quality of JOST over a wide field of view to inform the future imple- mentation of different wavefront sensing algorithms including the currently implemented Linearized Algorithm for Phase Diversity (LAPD). JOST complements other work at the Makidon Laboratory at the Space Telescope Science Institute, including the High-contrast imager for Complex Aperture Telescopes (HiCAT) testbed, that investigates coronagraphy for segmented aperture telescopes. Beyond JWST we intend to use JOST for WFS&C studies for future large segmented space telescopes such as LUVOIR.

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

The Gemini Planet Imager Exoplanet Survey (GPIES) is a direct imaging campaign designed to search for young, self-luminous, giant exoplanets. To date, GPIES has observed nearly 500 targets, and generated over 30,000 individual exposures using its integral field spectrograph (IFS) instrument. The GPIES team has developed a campaign data system with a database incorporating all of the metadata for all individual raw data products, including environmental conditions and instrument performance metrics. The same database also indexes metadata associated with multiple levels of reduced data products, including contrast measures for individual images and combined image sequences, which serve as the primary metric of performance for the final science products. The database is also used to track telemetry products from the adaptive optics subsystem, and associate these with corresponding IFS data. Here, we discuss several data exploration and visualization projects enabled by the GPIES database. Of particular interest are any correlations between instrument performance and environmental or operating conditions. We show single and multiple-parameter fits of single-image and observing sequence contrast as functions of various seeing measures, and discuss automated outlier rejection and other fitting concerns. Supervised learning techniques are employed in order to partition the space of raw (single image) to final (full sequence) contrast in order to better predict the value of the final data set from the first few completed observations. Finally, we discuss the particular features of the database design that aid in performing these analyses, and suggest potential future upgrades and refinements.

Accurately predicting the performance of coronagraphs and tolerancing optical surfaces for high-contrast imaging requires a detailed accounting of diffraction effects. Unlike simple Fraunhofer diffraction modeling, near and far-field diffraction effects, such as the Talbot effect, are captured by plane-to-plane propagation using Fresnel and angular spectrum propagation. This approach requires a sequence of computationally intensive Fourier transforms and quadratic phase functions, which limit the design and aberration sensitivity parameter space which can be explored at high-fidelity in the course of coronagraph design. This study presents the results of optimizing the multi-surface propagation module of the open source Physical Optics Propagation in PYthon (POPPY) package. This optimization was performed by implementing and benchmarking Fourier transforms and array operations on graphics processing units, as well as optimizing multithreaded numerical calculations using the NumExpr python library where appropriate, to speed the end-to-end simulation of observatory and coronagraph optical systems. Using realistic systems, this study demonstrates a greater than five-fold decrease in wall-clock runtime over POPPY's previous implementation and describes opportunities for further improvements in diffraction modeling performance.

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 present Hubble Space Telescope Near-Infrared Camera and Multi-Object Spectrometer coronagraphic imaging polarimetry of the TW Hydrae protoplanetary disk. These observations simultaneously measure the total and polarized intensity, allowing direct measurement of the polarization fraction across the disk. In accord with the self-shadowing hypothesis recently proposed by Debes et al., we find that the total and polarized intensity of the disk exhibits strong azimuthal asymmetries at projected distances consistent with the previously reported bright and dark ring-shaped structures (~45-99 au). The sinusoidal-like variations possess a maximum brightness at position angles near ~268-300 degrees and are up to ~28% stronger in total intensity. Furthermore, significant radial and azimuthal variations are also detected in the polarization fraction of the disk. In particular, we find that regions of lower polarization fraction are associated with annuli of increased surface brightness, suggesting that the relative proportion of multiple-to-single scattering is greater along the ring and gap structures. Moreover, we find strong (~20%) azimuthal variation in the polarization fraction along the shadowed region of the disk. Further investigation reveals that the azimuthal variation is not the result of disk flaring effects, but instead from a decrease in the relative contribution of multiple-to-single scattering within the shadowed region. Employing a two-layer scattering surface, we hypothesize that the diminished contribution in multiple scattering may result from shadowing by an inclined inner disk, which prevents direct stellar light from reaching the optically thick underlying surface component.

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

We present the first scattered-light images of two debris disks around the F8 star HD 104860 and the F0V star HD 192758, respectively $\sim45$ and $\sim67$ pc away. We detected these systems in the F110W and F160W filters through our re-analysis of archival Hubble Space Telescope NICMOS data with modern starlight subtraction techniques. Our image of HD 104860 confirms the morphology previously observed by Herschel in thermal emission with a well-defined ring at radius $\sim114$ au inclined $\sim58$ degrees. Although the outer edge profile is consistent with dynamical evolution models, the sharp inner edge suggests sculpting by unseen perturbers. Our images of HD 192758 reveal a disk at radius $\sim95$ au inclined by $\sim59$ degrees, never resolved so far. These disks have low scattering albedos of 10% and 13% respectively, inconsistent with water ice grain compositions. They are reminiscent of several other disks with similar inclination and scattering albedos: Fomalhaut, HD 92945, HD 202628, and HD 207129. They are also very distinct from brighter disks in the same inclination bin, which point to different compositions between these two populations. Varying scattering albedo values can be explained by different grain porosities, chemical compositions, or grain size distributions, which may indicate distinct formation mechanisms or dynamical processes at work in these systems. Finally, these faint disks with large infrared excesses may be representative of an underlying population of systems with low albedo values. Searches with more sensitive instruments on HST or on the James Webb Space Telescope and using state-of-the art starlight-subtraction methods may help discover more of such faint systems.

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.

The optically and IR bright, and starlight-scattering, HR 4796A ring-like debris disk is one of the most (and best) studied exoplanetary debris systems. The presence of a yet-undetected planet has been inferred (or suggested) from the narrow width and inner/outer truncation radii of its r = 1.05" (77 au) debris ring. We present new, highly sensitive, Hubble Space Telescope (HST) visible-light images of the HR 4796A circumstellar debris system and its environment over a very wide range of stellocentric angles from 0.32" (23 au) to ~ 15" (1100 au). These very high contrast images were obtained with the Space Telescope Imaging Spectrograph (STIS) using 6-roll PSF-template subtracted coronagraphy suppressing the primary light of HR 4796A and using three image plane occulters and simultaneously subtracting the background light from its close angular proximity M2.5V companion. The resulting images unambiguously reveal the debris ring embedded within a much larger, morphologically complex, and bi-axially asymmetric exoring scattering structure. These images at visible wavelengths are sensitive to, and map, the spatial distribution, brightness, and radial surface density of micron size particles over 5 dex in surface brightness. These particles in the exo-ring environment may be unbound from the system and interacting with the local ISM. Herein we present a new morphological and photometric view of the larger than prior seen HR 4796A exoplanetary debris system with sensitivity to small particles at stellocentric distances an order of magnitude greater than has previously been observed.

We present new HST ACS observations and detailed models for a recently discovered edge-on protoplanetary disk around ESO H$\alpha$ 569 (a low-mass T Tauri star in the Cha I star forming region). Using radiative transfer models we probe the distribution of the grains and overall shape of the disk (inclination, scale height, dust mass, flaring exponent and surface/volume density exponent) by model fitting to multiwavelength (F606W and F814W) HST observations together with a literature compiled spectral energy distribution. A new tool set was developed for finding optimal fits of MCFOST radiative transfer models using the MCMC code emcee to efficiently explore the high dimensional parameter space. It is able to self-consistently and simultaneously fit a wide variety of observables in order to place constraints on the physical properties of a given disk, while also rigorously assessing the uncertainties in those derived properties. We confirm that ESO H$\alpha$ 569 is an optically thick nearly edge-on protoplanetary disk. The shape of the disk is well described by a flared disk model with an exponentially tapered outer edge, consistent with models previously advocated on theoretical grounds and supported by millimeter interferometry. The scattered light images and spectral energy distribution are best fit by an unusually high total disk mass (gas+dust assuming a ratio of 100:1) with a disk-to-star mass ratio of 0.16.

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.

As the performance of coronagraphs improves, the achievable contrast is more and more dependent of the shape of the pupil. The future generation of space and ground based coronagraphic instruments will have to achieve high contrast levels on on-axis and/or segmented telescopes. To correct for the high amplitude aberrations introduced by secondary mirror structures and segmentation of the primary mirror, we explore a two deformable mirror (DM) method. The major difficulty of several DM methods is the non-linear relation linking actuator strokes to the point spread function in the coronagraph focal plane. The Active Compensation of Aperture Discontinuities (ACAD) method is achieving this minimization by solving a non linear differential Monge Ampere equation. Once this open loop method have reached the minimum, a close-loop stroke minimization method can be applied to correct for phase and amplitude aberrations to achieve the ultimate contrast. In this paper, I describe the results of the parametric analysis that that I have undertaken on this method. After recalling the principle of the method, I will described the explored parameter space (deformable mirror set-up, shape of the pupil, bandwidth, coronagraph designs). I will precisely described the way I simulated the Vortex coronagraph for this numerical simulation. Finally I will present the preliminary results of this parametric analysis for space telescope pupils only.

In the past 20 years, the Hubble Space Telescope (HST) STIS coronagraphic instrument has observed more than 100 stars, obtaining more than 4,000 readouts since its installment on HST in 1997 and the numbers are still increasing. We reduce the whole STIS coronagraphic archive at the most commonly observed positions (Wedge A0.6 and A1.0) with new post-processing methods, and present our results here. We are able to recover all of the 32 previously reported circumstellar disks, and obtain better contrast close to the star. For some of the disks, our results are limited by the over subtraction of the methods, and therefore the major regions of the disks can be recovered except the faintest regions. We also explain our efforts in the calibration of its new BAR5 occulting position, enabling STIS to explore inner regions as close as 0.2".

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.

In order to look for signs of on-going planet formation in young disks, we carried out the first J-band polarized emission imaging of the Herbig Ae/Be stars HD 150193, HD 163296, and HD 169142 using the Gemini Planet Imager (GPI), along with new H band observations of HD 144432. We confirm the complex "double ring" structure for the nearly face-on system HD 169142 first seen in H-band, finding the outer ring to be substantially redder than the inner one in polarized intensity. Using radiative transfer modeling, we developed a physical model that explains the full spectral energy distribution (SED) and J- and H-band surface brightness profiles, suggesting that the differential color of the two rings could come from reddened starlight traversing the inner wall and may not require differences in grain properties. In addition, we clearly detect an elongated, off-center ring in HD 163296 (MWC 275), locating the scattering surface to be 18 AU above the midplane at a radial distance of 77 AU, co-spatial with a ring seen at 1.3mm by ALMA linked to the CO snow line. Lastly, we report a weak tentative detection of scattered light for HD 150193 (MWC 863) and a non-detection for HD 144432; the stellar companion known for each of these targets has likely disrupted the material in the outer disk of the primary star. For HD 163296 and HD 169142, the prominent outer rings we detect could be evidence for giant planet formation in the outer disk or a manifestation of large-scale dust growth processes possibly related to snow-line chemistry.

We present the first scattered-light images of the debris disk around 49 ceti, a ~40 Myr A1 main sequence star at 59 pc, famous for hosting two massive dust belts as well as large quantities of atomic and molecular gas. The outer disk is revealed in reprocessed archival Hubble Space Telescope NICMOS F110W images, as well as new coronagraphic H band images from the Very Large Telescope SPHERE instrument. The disk extends from 1.1" (65 AU) to 4.6" (250 AU), and is seen at an inclination of 73degr, which refines previous measurements at lower angular resolution. We also report no companion detection larger than 3 M_Jup at projected separations beyond 20 AU from the star (0.34"). Comparison between the F110W and H-band images is consistent with a grey color of 49 ceti's dust, indicating grains larger than >2microns. Our photometric measurements indicate a scattering efficiency / infrared excess ratio of 0.2-0.4, relatively low compared to other characterized debris disks. We find that 49 ceti presents morphological and scattering properties very similar to the gas-rich HD 131835 system. From our constraint on the disk inclination we find that the atomic gas previously detected in absorption must extend to the inner disk, and that the latter must be depleted of CO gas. Building on previous studies, we propose a schematic view of the system describing the dust and gas structure around 49 ceti and hypothetic scenarios for the gas nature and origin.

Uncovering the ingredients and the architecture of planetary systems is a very active field of research that has fuelled many new theories on giant planet formation, migration, composition, and interaction with the circumstellar environment. We aim at discovering and studying new such systems, to further expand our knowledge of how low-mass companions form and evolve. We obtained high-contrast H-band images of the circumstellar environment of the F5V star HD206893, known to host a debris disc never detected in scattered light. These observations are part of the SPHERE High Angular Resolution Debris Disc Survey (SHARDDS) using the InfraRed Dual-band Imager and Spectrograph (IRDIS) installed on VLT/SPHERE. We report the detection of a source with a contrast of 3.6x10^{-5} in the H-band, orbiting at a projected separation of 270 milliarcsecond or 10 au, corresponding to a mass in the range 24 to 73 Mjup for an age of the system in the range 0.2 to 2 Gyr. The detection was confirmed ten months later with VLT/NaCo, ruling out a background object with no proper motion. A faint extended emission compatible with the disc scattered light signal is also observed. The detection of a low-mass companion inside a massive debris disc makes this system an analog of other young planetary systems such as beta Pictoris, HR8799 or HD95086 and requires now further characterisation of both components to understand their interactions.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

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 improvements to the wavelength calibration for the lenslet-based Integral Field Spectrograph (IFS), that serves as the science instrument for the Gemini Planet Imager (GPI). The GPI IFS features a 2.7" $\times$ 2.7" field of view and a 190 $\times$ 190 lenslet array (14.1 mas/lenslet) with spectral resolving power ranging from R $\sim$ 35 to 78. A unique wavelength solution is determined for each lenslet characterized by a two-dimensional position, an n-dimensional polynomial describing the spectral dispersion, and the rotation of the spectrum with respect to the detector axis. We investigate the non-linearity of the spectral dispersion across all Y, J, H, and K bands through both on-sky arc lamp images and simulated IFS images using a model of the optical path. Additionally, the 10-hole non-redundant masking mode on GPI provides an alternative measure of wavelength dispersion within a datacube by cross-correlating reference PSFs with science images. This approach can be used to confirm deviations from linear dispersion in the reduced datacubes. We find that the inclusion of a quadratic term provides a factor of 10 improvement in wavelength solution accuracy over the linear solution and is necessary to achieve uncertainties of a few hundredths of a pixel in J band to a few thousands of a pixel in the K bands. This corresponds to a wavelength uncertainty of $\sim$0.2 nm across all filters.

The Gemini Planet Imager (GPI) has been designed for the direct detection and characterization of exoplanets and circumstellar disks. GPI is equipped with a dual channel polarimetry mode designed to take advantage of the inherently polarized light scattered off circumstellar material to further suppress the residual seeing halo left uncorrected by the adaptive optics. We explore how recent advances in data reduction techniques reduce systematics and improve the achievable contrast in polarimetry mode. In particular, we consider different flux extraction techniques when constructing datacubes from raw data, division by a polarized flat-field and a method for subtracting instrumental polarization. Using observations of unpolarized standard stars we find that GPI's instrumental polarization is consistent with being wavelength independent within our errors. In addition, we provide polarimetry contrast curves that demonstrate typical performance throughout the GPIES campaign.

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

The Gemini Planet Imager is a high-contrast near-infrared instrument specifically designed to image exoplanets and circumstellar disks over a narrow field of view. We use science data and AO telemetry taken during the first 1.5 yr of the GPI Exoplanet Survey to quantify the performance of the AO system. In a typical 60 sec H-band exposure, GPI achieves a 5$\sigma$ raw contrast of 10$^{-4}$ at 0.4"; typical final 5$\sigma$ contrasts for full 1 hr sequences are more than 10 times better than raw contrasts. We find that contrast is limited by bandwidth wavefront error over much of the PSF. Preliminary exploratory factor analysis can explain 60-70% of the variance in raw contrasts with combinations of seeing and wavefront error metrics. We also examine the effect of higher loop gains on contrast by comparing wavefront error maps reconstructed from AO telemetry to concurrent IFS images. These results point to several ways that GPI performance could be improved in software or hardware.

The HD 61005 debris disk ("The Moth") stands out from the growing collection of spatially resolved circumstellar disks by virtue of its unusual swept-back morphology, brightness asymmetries, and dust ring offset. Despite several suggestions for the physical mechanisms creating these features, no definitive answer has been found. In this work, we demonstrate the plausibility of a scenario in which the disk material is shaped dynamically by an eccentric, inclined planet. We present new Keck NIRC2 scattered-light angular differential imaging of the disk at 1.2-2.3 microns that further constrains its outer morphology (projected separations of 27-135 AU). We also present complementary Gemini Planet Imager 1.6 micron total intensity and polarized light detections that probe down to projected separations less than 10 AU. To test our planet-sculpting hypothesis, we employed secular perturbation theory to construct parent body and dust distributions that informed scattered-light models. We found that this method produced models with morphological and photometric features similar to those seen in the data, supporting the premise of a planet-perturbed disk. Briefly, our results indicate a disk parent body population with a semimajor axis of 40-52 AU and an interior planet with an eccentricity of at least 0.2. Many permutations of planet mass and semimajor axis are allowed, ranging from an Earth mass at 35 AU to a Jupiter mass at 5 AU.

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 have imaged GM Aur with HST, detected its disk in scattered light at 1400A and 1650A, and compared these with observations at 3300A, 5550A, 1.1 microns, and 1.6 microns. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300A shows no evidence of the 24AU radius cavity that has been previously observed in sub-mm observations. Comparison with dust grain opacity models indicates the surface of the entire disk is populated with sub-micron grains. We have compiled an SED from 0.1 microns to 1 mm, and used it to constrain a model of the star+disk system that includes the sub-mm cavity using the Monte Carlo Radiative Transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of sub-microns grains interior to the sub-mm cavity wall. We suggest one explanation for this which could be due to a planet of mass <9 Jupiter masses interior to 24 AU. A unique cylindrical structure is detected in the FUV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semi-minor axis, but does not resemble an accretion-driven jet. The structure is limb-brightened and extends 190 +/- 35 AU above the disk midplane. The inner radius of the limb-brightening is 40 +/- 10 AU, just beyond the sub-millimeter cavity wall.

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

We present new Hubble Space Telescope observations of three a priori known starlight-scattering circumstellar debris systems (CDSs) viewed at intermediate inclinations around nearby close-solar analog stars: HD 207129, HD 202628, and HD 202917. Each of these CDSs possesses ring-like components that are more-massive analogs of our solar system's Edgeworth- Kuiper belt. These systems were chosen for follow-up observations to provide higher-fidelity and better sensitivity imaging for the sparse sample of solar-analog CDSs that range over two decades in systemic ages with HD 202628 and HD 202917 (both ~ 2.3 Gyr) currently the oldest CDSs imaged in visible or near-IR light. These deep (10 - 14 ksec) observations, with six-roll point-spread-function template subtracted visible-light coronagraphy using the Space Telescope Imaging Spectrograph, were designed to better reveal their angularly large, diffuse/low surface brightness, debris rings, and for all targets probe their exo-ring environments for starlight-scattering materials that present observational challenges for current ground-based facilities and instruments. Contemporaneously also observing with a narrower occulter position, these observations additionally probe the CDS endo-ring environments seen to be relatively devoid of scatterers. We discuss the morphological, geometrical, and photometric properties of these CDSs also in the context of other FGK-star hosted CDSs we have previously imaged as a homogeneously observed ensemble. From this combined sample we report a general decay in quiescent disk F_disk/F_star optical brightness ~ t^-0.8, similar to what is seen in at thermal IR wavelengths, and CDSs with a significant diversity in scattering phase asymmetries, and spatial distributions of their starlight-scattering grains.

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

We report the discovery of a scattering component around the HD 141569 A circumstellar debris system, interior to the previously known inner ring. The discovered inner disk component, obtained in broadband optical light with HST/STIS coronagraphy, was imaged with an inner working angle of 0".25, and can be traced from 0".4 (~46 AU) to 1".0 (~116 AU) after deprojection using i=55deg. The inner disk component is seen to forward scatter in a manner similar to the previously known rings, has a pericenter offset of ~6 AU, and break points where the slope of the surface brightness changes. It also has a spiral arm trailing in the same sense as other spiral arms and arcs seen at larger stellocentric distances. The inner disk spatially overlaps with the previously reported warm gas disk seen in thermal emission. We detect no point sources within 2" (~232 AU), in particular in the gap between the inner disk component and the inner ring. Our upper limit of 9+/-3 M_J is augmented by a new dynamical limit on single planetary mass bodies in the gap between the inner disk component and the inner ring of 1 M_J, which is broadly consistent with previous estimates.

Some circumstellar disks, called transitional or hybrid disks, present characteristics of both protoplanetary disks (significant amount of gas) and debris disks (evolved structures around young main-sequence stars, composed of second generation dust, from collisions between planetesimals). Therefore, they are an ideal astrophysical laboratory to witness the last stages of planet formation. The circumstellar disk around HD~141569A was intensively observed and resolved in the past from space but also from the ground but the recent implementation of high contrast imaging systems opens new opportunities to re-analyze this object. We analyzed Gemini archival data from the Near-Infrared Coronagraphic Imager (NICI) obtained in 2011 in the H band, using several angular differential imaging techniques (classical ADI, LOCI, KLIP). These images reveal the complex structures of this disk with an unprecedented resolution. We also include archival Hubble Space Telescope (HST) images as an independent dataset to confirm these findings. Using an analysis of the inner edge of the disk, we show that the inner disk is almost axisymmetrical. The measurement of an offset towards the East observed by previous authors is likely due to the fact that the Eastern part of this disk is wider and more complex in substructure. Our precise re-analysis of the eastern side shows several structures including a splitting of the disk and a small finger detached from the inner edge on the southeast. Finally, we find that the arc at 250 AU is unlikely to be a spiral, at least not at the inclination derived from the first ring, but instead could be interpreted as a third belt at a different inclination. If the very symmetrical inner disk edge is carved by a companion, the data presented here put additional constraints on its position.

We introduce a new class of solutions for Apodized Pupil Lyot Coronagraphs (APLC) with segmented aperture telescopes to remove broadband diffracted light from a star with a contrast level of $10^{10}$. These new coronagraphs provide a key advance to enabling direct imaging and spectroscopy of Earth twins with future large space missions. Building on shaped pupil (SP) apodization optimizations, our approach enables two-dimensional optimizations of the system to address any aperture features such as central obstruction, support structures or segment gaps. We illustrate the technique with a design that could reach $10^{10}$ contrast level at 34\,mas for a 12\,m segmented telescope over a 10\% bandpass centered at a wavelength $\lambda_0=$500\,nm. These designs can be optimized specifically for the presence of a resolved star, and in our example, for stellar angular size up to 1.1\,mas. This would allow probing the vicinity of Sun-like stars located beyond 4.4\,pc, therefore fully retiring this concern. If the fraction of stars with Earth-like planets is $\eta_{\Earth}=0.1$, with 18\% throughput, assuming a perfect, stable wavefront and considering photon noise only, 12.5 exo-Earth candidates could be detected around nearby stars with this design and a 12\,m space telescope during a five-year mission with two years dedicated to exo-Earth detection (one total year of exposure time and another year of overheads). Our new hybrid APLC/SP solutions represent the first numerical solution of a coronagraph based on existing mask technologies and compatible with segmented apertures, and that can provide contrast compatible with detecting and studying Earth-like planets around nearby stars. They represent an important step forward towards enabling these science goals with future large space missions.

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 the first images of four debris disks observed in scattered light around the young (4--250 Myr old) M dwarfs TWA 7 and TWA 25, the K6 star HD 35650, and the G2 star HD 377. We obtained these images by reprocessing archival Hubble Space Telescope NICMOS coronagraph data with modern post-processing techniques as part of the Archival Legacy Investigation of Circumstellar Environments (ALICE) program. All four disks appear faint and compact compared with other debris disks resolved in scattered light. The disks around TWA 25, HD 35650, and HD 377 appear very inclined, while TWA 7's disk is viewed nearly face-on. The surface brightness of HD 35650's disk is strongly asymmetric. These new detections raise the number of disks resolved in scattered light around M and late-K stars from one (the AU Mic system) to four. This new sample of resolved disks enables comparative studies of heretofore scarce debris disks around low-mass stars relative to solar-type stars.

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 coronagraphic instrument currently proposed for the WFIRST-AFTA mission will be the first example of a space-based coronagraph optimized for extremely high contrasts that are required for the direct imaging of exoplanets reflecting the light of their host star. While the design of this instrument is still in progress, this early stage of development is a particularly beneficial time to consider the operation of such an instrument. In this paper, we review current or planned operations on the Hubble Space Telescope (HST) and the James Webb Space Telescope (JWST) with a focus on which operational aspects will have relevance to the planned WFIRST-AFTA coronagraphic instrument. We identify five key aspects of operations that will require attention: 1) detector health and evolution, 2) wavefront control, 3) observing strategies/post-processing, 4) astrometric precision/target acquisition, and 5) polarimetry. We make suggestions on a path forward for each of these items.

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.

The ALICE program, for Archival Legacy Investigation of Circumstellar Environment, is currently conducting a virtual survey of about 400 stars, by re-analyzing the HST-NICMOS coronagraphic archive with advanced post-processing techniques. We present here the strategy that we adopted to identify detections and potential candidates for follow-up observations, and we give a preliminary overview of our detections. We present a statistical analysis conducted to evaluate the confidence level on these detection and the completeness of our candidate search.

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

The imaging channel on the Mid-Infrared Instrument (MIRI) is equipped with four coronagraphs that provide high contrast imaging capabilities for studying faint point sources and extended emission that would otherwise be overwhelmed by a bright point-source in its vicinity. Such bright sources might include stars that are orbited by exoplanets and circumstellar material, mass-loss envelopes around post-main-sequence stars, the near-nuclear environments in active galaxies, and the host galaxies of distant quasars. This paper describes the coronagraphic observing modes of MIRI, as well as performance estimates based on measurements of the MIRI flight model during cryo-vacuum testing. A brief outline of coronagraphic operations is also provided. Finally, simulated MIRI coronagraphic observations of a few astronomical targets are presented for illustration.

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

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

The newly commissioned Gemini Planet Imager (GPI) combines extreme adaptive optics, an advanced coronagraph, precision wavefront control and a lenslet-based integral field spectrograph (IFS) to measure the spectra of young extrasolar giant planets between 0.9-2.5 um. Each GPI detector image, when in spectral model, consists of ~37,000 microspectra which are under or critically sampled in the spatial direction. This paper demonstrates how to obtain high-resolution microlens PSFs and discusses their use in enhancing the wavelength calibration, flexure compensation and spectral extraction. This method is generally applicable to any lenslet-based integral field spectrograph including proposed future instrument concepts for space missions.

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 is an extreme AO instrument with an integral field spectrograph (IFS) operating in Y, J, H, and K bands. Both the Gemini telescope and the GPI instrument are very complex systems. Our goal is that the combined telescope and instrument system may be run by one observer operating the instrument, and one operator controlling the telescope and the acquisition of light to the instrument. This requires a smooth integration between the two systems and easily operated control interfaces. We discuss the definition of the software and hardware interfaces, their implementation and testing, and the integration of the instrument with the telescope environment.

The Gemini Planet Imager (GPI) Extreme Adaptive Optics Coronograph contains an interferometric mode: a 10-hole non-redundant mask (NRM) in its pupil wheel. GPI operates at $Y, J, H$, and $K$ bands, using an integral field unit spectrograph (IFS) to obtain spectral data at every image pixel. NRM on GPI is capable of imaging with a half resolution element inner working angle at moderate contrast, probing the region behind the coronagraphic spot. The fine features of the NRM PSF can provide a reliable check on the plate scale, while also acting as an attenuator for spectral standard calibrators that would otherwise saturate the full pupil. NRM commissioning data provides details about wavefront error in the optics as well as operations of adaptive optics control without pointing control from the calibration system. We compare lab and on-sky results to evaluate systematic instrument properties and examine the stability data in consecutive exposures. We discuss early on-sky performance, comparing images from integration and tests with the first on-sky images, and demonstrate resolving a known binary. We discuss the status of NRM and implications for future science with this mode.

The Gemini Planet Imager (GPI) is an instrument designed to directly image planets and circumstellar disks from 0.9 to 2.5 microns (the $YJHK$ infrared bands) using high contrast adaptive optics with a lenslet-based integral field spectrograph. We develop an extraction algorithm based on a least-squares method to disentangle the spectra and systematic noise contributions simultaneously. We utilize two approaches to adjust for the effect of flexure of the GPI optics which move the position of light incident on the detector. The first method is to iterate the extraction to achieve minimum residual and the second is to cross-correlate the detector image with a model image in iterative extraction steps to determine an offset. Thus far, this process has made clear qualitative improvements to the cube extraction by reducing the Moir\'{e} pattern. There are also improvements to the automated routines for finding flexure offsets which are reliable to with $\sim0.5$ pixel accuracy compared to pixel accuracy prior. Further testing and optimization will follow before implementation into the GPI pipeline.

The Gemini Planet Imager instrument's adaptive optics (AO) subsystem was designed specifically to facilitate high-contrast imaging. It features several new technologies, including computationally efficient wavefront reconstruction with the Fourier transform, modal gain optimization every 8 seconds, and the spatially filtered wavefront sensor. It also uses a Linear-Quadratic-Gaussian (LQG) controller (aka Kalman filter) for both pointing and focus. We present on-sky performance results from verification and commissioning runs from December 2013 through May 2014. The efficient reconstruction and modal gain optimization are working as designed. The LQG controllers effectively notch out vibrations. The spatial filter can remove aliases, but we typically use it oversized by about 60% due to stability problems.

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.

The Gemini Planet Imager is a newly commissioned facility instrument designed to measure the near-infrared spectra of young extrasolar planets in the solar neighborhood and obtain imaging polarimetry of circumstellar disks. GPI's science instrument is an integral field spectrograph that utilizes a HAWAII-2RG detector with a SIDECAR ASIC readout system. This paper describes the detector characterization and calibrations performed by the GPI Data Reduction Pipeline to compensate for effects including bad/hot/cold pixels, persistence, non-linearity, vibration induced microphonics and correlated read noise.

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.

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

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

We present the wavelength calibration for the lenslet-based Integral Field Spectrograph (IFS) that serves as the science instrument for the Gemini Planet Imager (GPI). The GPI IFS features a 2.7" x 2.7" field of view and a 190 x 190 lenslet array (14.3 mas/lenslet) operating in $Y$, $J$, $H$, and $K$ bands with spectral resolving power ranging from $R$ $\sim$ 35 to 78. Due to variations across the field of view, a unique wavelength solution is determined for each lenslet characterized by a two-dimensional position, the spectral dispersion, and the rotation of the spectrum with respect to the detector axes. The four free parameters are fit using a constrained Levenberg-Marquardt least-squares minimization algorithm, which compares an individual lenslet's arc lamp spectrum to a simulated arc lamp spectrum. This method enables measurement of spectral positions to better than 1/10th of a pixel on the GPI IFS detector using Gemini's facility calibration lamp unit GCAL, improving spectral extraction accuracy compared to earlier approaches. Using such wavelength calibrations we have measured how internal flexure of the spectrograph with changing zenith angle shifts spectra on the detector. We describe the methods used to compensate for these shifts when assembling datacubes from on-sky observations using GPI.

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.

An Atmospheric Dispersion Corrector (ADC) uses a double-prism arrangement to nullify the vertical chromatic dispersion introduced by the atmosphere at non-zero zenith distances. The ADC installed in the Gemini Planet Imager (GPI) was first tested in August 2012 while the instrument was in the laboratory. GPI was installed at the Gemini South telescope in August 2013 and first light occurred later that year on November 11th. In this paper, we give an overview of the characterizations and performance of this ADC unit obtained in the laboratory and on sky, as well as the structure of its control software.

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

Searching for nearby habitable worlds with direct imaging and spectroscopy will require a telescope large enough to provide angular resolution and sensitivity to planets around a significant sample of stars. Segmented telescopes are a compelling option to obtain such large apertures. However, these telescope designs have a complex geometry (central obstruction, support structures, segmentation) that makes high-contrast imaging more challenging. We are developing a new high-contrast imaging testbed at STScI to provide an integrated solution for wavefront control and starlight suppression on complex aperture geometries. We present our approach for the testbed optical design, which defines the surface requirements for each mirror to minimize the amplitude-induced errors from the propagation of out-of-pupil surfaces. Our approach guarantees that the testbed will not be limited by these Fresnel propagation effects, but only by the aperture geometry. This approach involves iterations between classical ray-tracing optical design optimization, and end-to-end Fresnel propagation with wavefront control (e.g. Electric Field Conjugation / Stroke Minimization). The construction of the testbed is planned to start in late Fall 2013.

Spatially resolved scattered-light images of circumstellar (CS) debris in exoplanetary systems constrain the physical properties and orbits of the dust particles in these systems. They also inform on co-orbiting (but unseen) planets, systemic architectures, and forces perturbing starlight-scattering CS material. Using HST/STIS optical coronagraphy, we have completed the observational phase of a program to study the spatial distribution of dust in ten CS debris systems, and one "mature" protoplanetrary disk all with HST pedigree, using PSF-subtracted multi-roll coronagraphy. These observations probe stellocentric distances > 5 AU for the nearest stars, and simultaneously resolve disk substructures well beyond, corresponding to the giant planet and Kuiper belt regions in our Solar System. They also disclose diffuse very low-surface brightness dust at larger stellocentric distances. We present new results inclusive of fainter disks such as HD92945 confirming, and better revealing, the existence of a narrow inner debris ring within a larger diffuse dust disk. Other disks with ring-like sub-structures, significant asymmetries and complex morphologies include: HD181327 with a posited spray of ejecta from a recent massive collision in an exo-Kuiper belt; HD61005 suggested interacting with the local ISM; HD15115 & HD32297, discussed also in the context of environmental interactions. These disks, and HD15745, suggest debris system evolution cannot be treated in isolation. For AU Mic's edge-on disk, out-of-plane surface brightness asymmetries at > 5 AU may implicate one or more planetary perturbers. Time resolved images of the MP Mus proto-planetary disk provide spatially resolved temporal variability in the disk illumination. These and other new images from our program enable direct inter-comparison of the architectures of these exoplanetary debris systems in the context of our own Solar System.

We present a new high-contrast imaging testbed designed to provide complete solutions in wavefront sensing, control and starlight suppression with complex aperture telescopes. The testbed was designed to enable a wide range of studies of the effects of such telescope geometries, with primary mirror segmentation, central obstruction, and spiders. The associated diffraction features in the point spread function make high-contrast imaging more challenging. In particular the testbed will be compatible with both AFTA-like and ATLAST-like aperture shapes, respectively on-axis monolithic, and on-axis segmented telescopes. The testbed optical design was developed using a novel approach to define the layout and surface error requirements to minimize amplitude-induced errors at the target contrast level performance. In this communication we compare the as-built surface errors for each optic to their specifications based on end-to-end Fresnel modeling of the testbed. We also report on the testbed optical and optomechanical alignment performance, coronagraph design and manufacturing, and preliminary first light results.

We are currently conducting a comprehensive and consistent re-processing of archival HST-NICMOS coronagraphic surveys using advanced PSF subtraction methods, entitled the Archival Legacy Investigations of Circumstellar Environments program (ALICE, HST/AR 12652). This virtual campaign of about 400 targets has already produced numerous new detections of previously unidentified point sources and circumstellar structures. We present five newly spatially resolved debris disks revealed in scattered light by our analysis of the archival data. These images provide new views of material around young solar-type stars at ages corresponding to the period of terrestrial planet formation in our solar system. We have also detected several new candidate substellar companions, for which there are ongoing followup campaigns (HST/WFC3 and VLT/SINFONI in ADI mode). Since the methods developed as part of ALICE are directly applicable to future missions (JWST, AFTA coronagraph) we emphasize the importance of devising optimal PSF subtraction methods for upcoming coronagraphic imaging missions. We describe efforts in defining direct imaging high-level science products (HLSP) standards that can be applicable to other coronagraphic campaigns, including ground-based (e.g., Gemini Planet Imager), and future space instruments (e.g., JWST). ALICE will deliver a first release of HLSPs to the community through the MAST archive at STScI in 2014.

The James Webb Space Telescope (JWST) Optical Simulation Testbed (JOST) is a tabletop workbench to study aspects of wavefront sensing and control for a segmented space telescope, including both commissioning and maintenance activities. JOST is complementary to existing optomechanical testbeds for JWST (e.g. the Ball Aerospace Testbed Telescope, TBT) given its compact scale and flexibility, ease of use, and colocation at the JWST Science & Operations Center. We have developed an optical design that reproduces the physics of JWST's three-mirror anastigmat using three aspheric lenses; it provides similar image quality as JWST (80% Strehl ratio) over a field equivalent to a NIRCam module, but at HeNe wavelength. A segmented deformable mirror stands in for the segmented primary mirror and allows control of the 18 segments in piston, tip, and tilt, while the secondary can be controlled in tip, tilt and x, y, z position. This will be sufficient to model many commissioning activities, to investigate field dependence and multiple field point sensing & control, to evaluate alternate sensing algorithms, and develop contingency plans. Testbed data will also be usable for cross-checking of the WFS&C Software Subsystem, and for staff training and development during JWST's five- to ten-year mission.

The James Webb Space Telescope (JWST) Optical Simulation Testbed (JOST) is a tabletop experiment designed to reproduce the main aspects of wavefront sensing and control (WFSC) for JWST. To replicate the key optical physics of JWST's three-mirror anastigmat (TMA) design at optical wavelengths we have developed a three-lens anastigmat optical system. This design uses custom lenses (plano-convex, plano-concave, and bi-convex) with fourth-order aspheric terms on powered surfaces to deliver the equivalent image quality and sampling of JWST NIRCam at the WFSC wavelength (633~nm, versus JWST's 2.12~micron). For active control, in addition to the segmented primary mirror simulator, JOST reproduces the secondary mirror alignment modes with five degrees of freedom. We present the testbed requirements and its optical and optomechanical design. We study the linearity of the main aberration modes (focus, astigmatism, coma) both as a function of field point and level of misalignments of the secondary mirror. We find that the linearity with the transmissive design is similar to what is observed with a traditional TMA design, and will allow us to develop a linear-control alignment strategy based on the multi-field methods planned for JWST.

We have spatially resolved five debris disks (HD 30447, HD 35841, HD 141943, HD 191089, and HD 202917) for the first time in near-infrared scattered light by reanalyzing archival Hubble Space Telescope (HST)/NICMOS coronagraphic images obtained between 1999 and 2006. One of these disks (HD 202917) was previously resolved at visible wavelengths using HST/Advanced Camera for Surveys. To obtain these new disk images, we performed advanced point-spread function subtraction based on the Karhunen-Loeve Image Projection (KLIP) algorithm on recently reprocessed NICMOS data with improved detector artifact removal (Legacy Archive PSF Library And Circumstellar Environments Legacy program). Three of the disks (HD 30447, HD 35841, and HD 141943) appear edge-on, while the other two (HD 191089 and HD 202917) appear inclined. The inclined disks have been sculpted into rings; in particular, the disk around HD 202917 exhibits strong asymmetries. All five host stars are young (8-40 Myr), nearby (40-100 pc) F and G stars, and one (HD 141943) is a close analog to the young sun during the epoch of terrestrial planet formation. Our discoveries increase the number of debris disks resolved in scattered light from 19 to 23 (a 21% increase). Given their youth, proximity, and brightness (V = 7.2 to 8.5), these targets are excellent candidates for follow-up investigations of planet formation at visible wavelengths using the HST/STIS coronagraph, at near-infrared wavelengths with the Gemini Planet Imager (GPI) and Very Large Telescope (VLT)/SPHERE, and at thermal infrared wavelengths with the James Webb Space Telescope NIRCam and MIRI coronagraphs.

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.

PDS 144 is a pair of Herbig Ae stars that are separated by 5.35" on the sky. It has previously been shown to have an A2Ve Herbig Ae star viewed at 83\circ inclination as its northern member and an A5Ve Herbig Ae star as its southern member. Direct imagery revealed a disk occulting PDS 144 N - the first edge-on disk observed around a Herbig Ae star. The lack of an obvious disk in direct imagery suggested PDS 144 S might be viewed face-on or not physically associated with PDS 144 N. Multi-epoch Hubble Space Telescope imagery of PDS 144 with a 5 year baseline demonstrates PDS 144 N & S are comoving and have a common proper motion with TYC 6782-878-1. TYC 6782-878-1 has previously been identified as a member of Upper Sco sub-association A at d = 145 \pm 2 pc with an age of 5-10 Myr. Ground-based imagery reveals jets and a string of Herbig-Haro knots extending 13' (possibly further) which are aligned to within 7\circ \pm 6\circ on the sky. By combining proper motion data and the absence of a dark mid-plane with radial velocity data, we measure the inclination of PDS 144 S to be i = 73\circ \pm 7\circ. The radial velocity of the jets from PDS 144 N & S indicates they, and therefore their disks, are misaligned by 25\circ \pm 9\circ. This degree of misalignment is similar to that seen in T Tauri wide binaries.

We report new observations of the unusually active, high proper motion L5e dwarf 2MASS J13153094-2649513. Optical spectroscopy with Magellan/MagE reveals persistent nonthermal emission, with narrow H I Balmer, Na I and K I lines all observed in emission. Low-resolution near-infrared spectroscopy with IRTF/SpeX indicates the presence of a low-temperature companion, which is resolved through multi-epoch laser guide star adaptive optics imaging at Keck. The comoving companion is separated by 338 \pm 4 mas, and its relative brightness (\Delta K_s = 5.09 \pm 0.10) makes this system the second most extreme flux ratio very low-mass binary identified to date. Resolved near-infrared spectroscopy with Keck/OSIRIS identifies this companion as a T7 dwarf. The absence of Li I absorption in combined-light optical spectroscopy constrains the system age to >~0.8-1.0 Gyr, while the system's kinematics and unusually low mass ratio (M_2/M_1 = 0.3-0.6) suggests that it is even older. A coevality test of the components also indicates an older age, but reveals discrepancies between evolutionary and atmosphere model fits of the secondary which are likely attributable to poor reproduction of its near-infrared spectrum. With a projected separation of 6.6 \pm 0.9 AU, the 2MASS J13153094-2649513 system is too widely separated for mass exchange or magnetospheric interactions to be powering its persistent nonthermal emission. Rather, the emission is probably chromospheric in nature, signaling an inversion in the age-activity relation in which strong magnetic fields are maintained by relatively old and massive ultracool dwarfs.

We are undertaking a multi-frequency Expanded Very Large Array (EVLA) survey of edge-on protoplanetary disks to probe the growth of solids in each disk, sedimentation of such material into the disk midplane, and the connection of these phenomena to the planet formation process. The projection of edge-on disk systems along our line of sight enables a study of the vertical stratification of large grains with fewer model dependencies than would be required for disks that are more face-on. Robust studies of the spatial distribution of grains up to ~1 cm in size are possible with the wavelength range and sensitivity of the EVLA. In this contribution we describe target selection and observational strategies. First results concerning the Class 0 source IRAS04368+2557 (L1527 IRS) are presented, including a study of this source's 8.46 GHz continuum variability over short and long time baselines and an indication that its protoplanetary disk may have a dearth of pebble-sized grains.

We present high-resolution ($0\farcs2$) near-infrared images of the area surrounding the Herbig Be star LkH$\alpha$ 234 taken with the Coronagraphic Imager with Adaptive Optics (CIAO) and the adaptive optics on the Subaru Telescope. The near-infrared ($J$, $H$, $K$, $L'$ and $M'$ bands) images reveal circumstellar structures around LkH$\alpha$ 234 in detail. Eight young stellar object (YSO) candidates (Object B, C, D, E, F, G, NW1, and NW2) were detected at 2$\arcsec$--11$\arcsec$ from LkH$\alpha$ 234. Objects B and C are likely variable stars, which is consistent with a young evolutionary status. Three objects (LkH$\alpha$ 234, NW1, and NW2) were identified in the 11.8 $\micron$ (SiC) and 17.65 $\micron$ images obtained with the Long Wavelength Spectrometer (LWS) on the Keck Telescope. NW1 and NW2 are thought to be embedded young stars. We suggest that NW1 and not LkH$\alpha$ 234 is the source illuminating the reflection nebula west of LkH$\alpha$ 234, although Object G may be another candidate. In our images, these objects are located at the center of the 2 $\micron$ polarization, and NW1 resides at the center of a monopolar cavity. The SED of NW1 suggests that it is a YSO with a spectral type of B6--7. Object F and G were discovered by our observations at $1\farcs9$ and $2\farcs3$, respectively, from LkH$\alpha$ 234, and their proximity to LkH$\alpha$ 234 suggests that they may be its lower-mass companions.

The nearby star Alpha Oph (Ras Alhague) is a rapidly rotating A5IV star spinning at ~89% of its breakup velocity. This system has been imaged extensively by interferometric techniques, giving a precise geometric model of the star's oblateness and the resulting temperature variation on the stellar surface. Fortuitously, Alpha Oph has a previously known stellar companion, and characterization of the orbit provides an independent, dynamically-based check of both the host star and the companion mass. Such measurements are crucial to constrain models of such rapidly rotating stars. In this study, we combine eight years of Adaptive Optics imaging data from the Palomar, AEOS, and CFHT telescopes to derive an improved, astrometric characterization of the companion orbit. We also use photometry from these observations to derive a model-based estimate of the companion mass. A fit was performed on the photocenter motion of this system to extract a component mass ratio. We find masses of 2.40^{0.23}_{0.37} solar masses and 0.85^{0.06}_{0.04} solar masses for Alpha Oph A and Alpha Oph B, respectively. Previous orbital studies of this system found a mass too high for this system, inconsistent with stellar evolutionary calculations. Our measurements of the host star mass are more consistent with these evolutionary calculations, but with slightly higher uncertainties. In addition to the dynamically-derived masses, we use IJHK photometry to derive a model-based mass for Alpha Oph B, of 0.77 +/- 0.05 solar masses marginally consistent with the dynamical masses derived from our orbit. Our model fits predict a periastron passage on 2012 April 19, with the two components having a ~50 milliarcsec separation from March to May 2012. A modest amount of interferometric and radial velocity data during this period could provide a mass determination of this star at the few percent level.

The Z CMa binary is understood to undergo both FU Orionis (FUOR) and EX Orionis (EXOR) type outbursts. While the SE component has been spectro- scopically identified as an FUOR, the NW component, a Herbig Be star, is the source of the EXOR outbursts. The system has been identified as the source of a large outflow, however, previous studies have failed to identify the driver. Here we present adaptive optics (AO) assisted [FeII] spectro-images which reveal for the first time the presence of two jets. Observations made using OSIRIS at the Keck Observatory show the Herbig Be star to be the source of the parsec-scale outflow, which within 2'' of the source shows signs of wiggling and the FUOR to '' be driving a ~ 0.4 jet. The wiggling of the Herbig Be star's jet is evidence for an additional companion which could in fact be generating the EXOR outbursts, the last of which began in 2008 (Grankin & Artemenko 2009). Indeed the dy- namical scale of the wiggling corresponds to a time-scale of 4-8 years which is in agreement with the time-scale of these outbursts. The spectro-images also show a bow-shock shaped feature and possible associated knots. The origin of this structure is as of yet unclear. Finally interesting low velocity structure is also observed. One possibility is that it originates in a wide-angle outflow launched from a circumbinary disk.

BP Psc is a remarkable emission-line field star that is orbited by a dusty disk and drives a parsec-scale system of jets. We report the detection by the Chandra X-ray Observatory of a weak X-ray point source coincident with the centroids of optical/IR and submillimeter continuum emission at BP Psc. As the star's photosphere is obscured throughout the visible and near-infrared, the Chandra X-ray source likely represents the first detection of BP Psc itself. The X-rays most likely originate with magnetic activity at BP Psc and hence can be attributed either to a stellar corona or to star-disk interactions. The log of the ratio of X-ray to bolometric luminosity (log(L_X/L_{bol}) lies in the range -5.8 to -4.2. This is smaller than log(L_X/L_{bol}) ratios typical of low-mass, pre-main sequence stars, but is well within the log(L_X/L_{bol}) range observed for rapidly-rotating (FK Com-type) G giant stars. Hence, the Chandra results favor an exotic model wherein the disk/jet system of BP Psc is the result of its very recently engulfing a companion star or giant planet, as the primary star ascended the giant branch.

The Lyot project used an optimized Lyot coronagraph with Extreme Adaptive Optics at the 3.63m Advanced Electro-Optical System telescope (AEOS) to observe 86 stars from 2004 to 2007. In this paper we give an overview of the survey results and a statistical analysis of the observed nondetections around 58 of our targets to place constraints on the population of substellar companions to nearby stars. The observations did not detect any companion in the substellar regime. Since null results can be as important as detections, we analyzed each observation to determine the characteristics of the companions that can be ruled out. For this purpose we use a Monte Carlo approach to produce artificial companions, and determine their detectability by comparison with the sensitivity curve for each star. All the non-detection results are combined using a Bayesian approach and we provide upper limits on the population of giant exoplanets and brown dwarfs for this sample of stars. Our nondetections confirm the rarity of brown dwarfs around solar-like stars and we constrain the frequency of massive substellar companions (M>40Mjup) at orbital separation between and 10 and 50 AU to be <20%.

Through the combination of high-order Adaptive Optics and coronagraphy, we report the discovery of a faint stellar companion to the A3V star zeta Virginis. This companion is ~7 magnitudes fainter than its host star in the H-band, and infrared imaging spanning 4.75 years over five epochs indicates this companion has common proper motion with its host star. Using evolutionary models, we estimate its mass to be 0.168+/-.016 solar masses, giving a mass ratio for this system q = 0.082. Assuming the two objects are coeval, this mass suggests a M4V-M7V spectral type for the companion, which is confirmed through integral field spectroscopic measurements. We see clear evidence for orbital motion from this companion and are able to constrain the semi-major axis to be greater than 24.9 AU, the period > 124$ yrs, and eccentricity > 0.16. Multiplicity studies of higher mass stars are relatively rare, and binary companions such as this one at the extreme low end of the mass ratio distribution are useful additions to surveys incomplete at such a low mass ratio. Moreover, the frequency of binary companions can help to discriminate between binary formation scenarios that predict an abundance of low-mass companions forming from the early fragmentation of a massive circumstellar disk. A system such as this may provide insight into the anomalous X-ray emission from A stars, hypothesized to be from unseen late-type stellar companions. Indeed, we calculate that the presence of this M-dwarf companion easily accounts for the X-ray emission from this star detected by ROSAT.

With the Keck Interferometer, we have studied at 2 um the innermost regions of several nearby, young, dust depleted "transitional" disks. Our observations target five of the six clearest cases of transitional disks in the Taurus/Auriga star-forming region (DM Tau, GM Aur, LkCa 15, UX Tau A, and RY Tau) to explore the possibility that the depletion of optically thick dust from the inner disks is caused by stellar companions rather than the more typical planet-formation hypothesis. At the 99.7% confidence level, the observed visibilities exclude binaries with flux ratios of at least 0.05 and separations ranging from 2.5 to 30 mas (0.35 - 4 AU) over >= 94% of the area covered by our measurements. All targets but DM Tau show near-infrared excess in their SED higher than our companion flux ratio detection limits. While a companion has previously been detected in the candidate transitional disk system CoKu Tau/4, we can exclude similar mass companions as the typical origin for the clearing of inner dust in transitional disks and of the near-infrared excess emission. Unlike CoKu Tau/4, all our targets show some evidence of accretion. We find that all but one of the targets are clearly spatially resolved, and UX Tau A is marginally resolved. Our data is consistent with hot material on small scales (0.1 AU) inside of and separated from the cooler outer disk, consistent with the recent SED modeling. These observations support the notion that some transitional disks have radial gaps in their optically thick material, which could be an indication for planet formation in the habitable zone (~ a few AU) of a protoplanetary disk.

Using the NICMOS coronagraph, we have obtained high-contrast 2.0 micron imaging polarimetry and 1.1 micron imaging of the circumstellar disk around AB Aurigae on angular scales of 0.3-3 arcsec (40-550 AU). Unlike previous observations, these data resolve the disk in both total and polarized intensity, allowing accurate measurement of the spatial variation of polarization fraction across the disk. Using these observations we investigate the apparent "gap" in the disk reported by Oppenheimer et al. 2008. In polarized intensity, the NICMOS data closely reproduces the morphology seen by Oppenheimer et al., yet in total intensity we find no evidence for a gap in either our 1.1 or 2.0 micron images. We find instead that region has lower polarization fraction, without a significant decrease in total scattered light, consistent with expectations for back-scattered light on the far side of an inclined disk. Radiative transfer models demonstrate this explanation fits the observations. Geometrical scattering effects are entirely sufficient to explain the observed morphology without any need to invoke a gap or protoplanet at that location.

We demonstrate the versatility of a dual imaging polarimeter working in tandem with a Lyot coronagraph and Adaptive Optics to suppress the highly static speckle noise pattern--the greatest hindrance to ground-based direct imaging of planets and disks around nearby stars. Using a double difference technique with the polarimetric data, we quantify the level of speckle suppression, and hence improved sensitivity, by placing an ensemble of artificial faint companions into real data, with given total brightness and polarization. For highly polarized sources within 0.5 arcsec, we show that we achieve 3 to 4 magnitudes greater sensitivity through polarimetric speckle suppression than simply using a coronagraph coupled to a high-order Adaptive Optics system. Using such a polarimeter with a classical Lyot coronagraph at the 3.63m AEOS telescope, we have obtained a 6.5 sigma detection in the H-band of the 76 AU diameter circumstellar debris disk around the star HR 4796A. Our data represent the first definitive, ground-based, near-IR polarimetric image of the HR 4796A debris disk and clearly show the two outer ansae of the disk, evident in Hubble Space Telescope NICMOS/STIS imaging. We derive a lower limit to the fractional linear polarization of 29% caused by dust grains in the disk. In addition, we fit simple morphological models of optically thin disks to our data allowing us to constrain the dust disk scale height to 2.5{+5.0}_{-1.3} AU and scattering asymmetry parameter (g=0.20^{+.07}_{-.10}). These values are consistent with several lines of evidence suggesting that the HR 4796A disk is dominated by a micron-sized dust population, and are indeed typical of disks in transition between those surrounding the Herbig Ae stars to those associated with Vega-like stars.

The mid-infrared properties of pre-planetary disks are sensitive to the temperature and flaring profiles of disks for the regions where planet formation is expected to occur. In order to constrain theories of planet formation, we have carried out a mid-infrared (wavelength 10.7 microns) size survey of young stellar objects using the segmented Keck telescope in a novel configuration. We introduced a customized pattern of tilts to individual mirror segments to allow efficient sparse-aperture interferometry, allowing full aperture synthesis imaging with higher calibration precision than traditional imaging. In contrast to previous surveys on smaller telescopes and with poorer calibration precision, we find most objects in our sample are partially resolved. Here we present the main observational results of our survey of 5 embedded massive protostars, 25 Herbig Ae/Be stars, 3 T Tauri stars, 1 FU Ori system, and 5 emission-line objects of uncertain classification. The observed mid-infrared sizes do not obey the size-luminosity relation found at near-infrared wavelengths and a companion paper will provide further modelling analysis of this sample. In addition, we report imaging results for a few of the most resolved objects, including complex emission around embedded massive protostars, the photoevaporating circumbinary disk around MWC 361A, and the subarcsecond binaries T Tau, FU Ori and MWC 1080.

Non-redundant masking (NRM) is a high contrast, high resolution technique relevant to future space missions concerned with extrasolar planetary system and star formation, as well as general high angular resolution galactic and extragalactic astronomy. NRM enables the highest angular resolution science possible given the telescope's diameter and operating wavelength. It also provides precise information on a telescope's optical state. We must assess NRM contrast limits realistically to understand the science yield of NRM in space, and, simultaneously, develop NRM science for planet and star formation and extragalactic science in the UV-NIR, to help steer high resolution space-based astronomy in the coming decade.

MWC 778 is an unusual and little-studied young stellar object located in the IC 2144 nebula. Recent spectroscopy by Herbig and Vacca (2008) suggested the presence of an edge-on circumstellar disk around it. We present near-infrared adaptive optics imaging polarimetry and mid-infrared imaging which directly confirm the suspected nearly-edge-on disk around MWC 778 (i ~ 70-80 degrees) plus reveal a more extensive envelope pierced by bipolar outflow cavities. In addition, our mid-infrared images and near-infrared polarization maps detect a spiral-shaped structure surrounding MWC 778, with arms that extend beyond 6" on either side of the star. Although MWC 778 has previously been classified as an Herbig Ae/Be star, the properties of its central source (including its spectral type) remain fairly uncertain. Herbig & Vacca (2008) suggested an F or G spectral type based on the presence of metallic absorption lines in the optical spectrum, which implies that MWC 778 may belong to the fairly rare class of Intermediate-Mass T Tauri Stars (IMTTSs) which are the evolutionary precursors to Herbig Ae/Be objects. Yet its integrated bolometric luminosity, > 750 L_sun (for an assumed distance of 1 kpc) is surprisingly high for an F or G spectral type, even for an IMTTS. We speculate on several possible explanations for this discrepancy, including its true distance being much closer than 1 kpc, the presence of a binary companion, and/or a non-stellar origin for the observed absorption lines.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of planetary systems and star formation. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of cosmology. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of stars and stellar evolution. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of Galactic science. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of extragalactic astrophysics. Community-based White Paper to Astro2010 in response to the call for such papers.

We present a panchromatic study, involving a multiple technique approach, of the circumstellar disc surrounding the T Tauri star IM Lupi (Sz 82). We have undertaken a comprehensive observational study of IM Lupi using photometry, spectroscopy, millimetre interferometry and multi-wavelength imaging. For the first time, the disc is resolved from optical and near-infrared wavelengths in scattered light, to the millimetre regime in thermal emission. Our data-set, in conjunction with existing photometric data, provides an extensive coverage of the spectral energy distribution, including a detailed spectrum of the silicate emission bands. We have performed a simultaneous modelling of the various observations, using the radiative transfer code MCFOST, and analysed a grid of models over a large fraction of the parameter space via Bayesian inference. We have constructed a model that can reproduce all of the observations of the disc. Our analysis illustrates the importance of combining a wide range of observations in order to fully constrain the disc model, with each observation providing a strong constraint only on some aspects of the disc structure and dust content. Quantitative evidence of dust evolution in the disc is obtained: grain growth up to millimetre-sized particles, vertical stratification of dust grains with micrometric grains close to the disc surface and larger grains which have settled towards the disc midplane, and possibly the formation of fluffy aggregates and/or ice mantles around grains.

We have upgraded IRCAL, the near-infrared science camera of the Lick Observatory adaptive optics system, to add a dual-channel imaging polarimetry mode. This mode uses an optically contacted YLF (LiYF_4) Wollaston prism to provide simultaneous images in perpendicular linear polarizations, providing high resolution, high dynamic range polarimetry in the near infrared. We describe the design and construction of the polarimeter, discuss in detail the data reduction algorithms adopted, and evaluate the instrument's on-the-sky performance. The IRCAL polarimeter is capable of reducing the stellar PSF halo by about two orders of magnitude, thereby increasing contrast for studies of faint circumstellar dust-scattered light. We discuss the various factors that limit the achieved contrast, and present lessons applicable to future high contrast imaging polarimeters.

The young star AB Aurigae is surrounded by a complex combination of gas-rich and dust dominated structures. The inner disk which has not been studied previously at sufficient resolution and imaging dynamic range seems to contain very little gas inside a radius of least 130 astronomical units (AU) from the star. Using adaptive-optics coronagraphy and polarimetry we have imaged the dust in an annulus between 43 and 302 AU from the star, a region never seen before. An azimuthal gap in an annulus of dust at a radius of 102 AU, along with a clearing at closer radii inside this annulus, suggests the formation of at least one small body at an orbital distance of about 100 AU. This structure seems consistent with crude models of mean motion resonances, or accumulation of material at two of the Lagrange points relative to the putative object and the star. We also report a low significance detection of a point source in this outer annulus of dust. This source may be an overdensity in the disk due to dust accreting onto an unseen companion. An alternate interpretation suggests that the object's mass is between 5 and 37 times the mass of Jupiter. The results have implications for circumstellar disk dynamics and planet formation.

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

Observations of outflows associated with pre-main-sequence stars reveal details about morphology, binarity and evolutionary states of young stellar objects. We present molecular line data from the Berkeley-Illinois-Maryland Association array and Five Colleges Radio Astronomical Observatory toward the regions containing the Herbig Ae/Be stars LkHa 198 and LkHa 225S. Single dish observations of 12CO 1-0, 13CO 1-0, N2H+ 1-0 and CS 2-1 were made over a field of 4.3' x 4.3' for each species. 12CO data from FCRAO were combined with high resolution BIMA array data to achieve a naturally-weighted synthesized beam of 6.75'' x 5.5'' toward LkHa 198 and 5.7'' x 3.95'' toward LkHa 225S, representing resolution improvements of factors of approximately 10 and 5 over existing data. By using uniform weighting, we achieved another factor of two improvement. The outflow around LkHa 198 resolves into at least four outflows, none of which are centered on LkHa 198-IR, but even at our resolution, we cannot exclude the possibility of an outflow associated with this source. In the LkHa 225S region, we find evidence for two outflows associated with LkHa 225S itself and a third outflow is likely driven by this source. Identification of the driving sources is still resolution-limited and is also complicated by the presence of three clouds along the line of sight toward the Cygnus molecular cloud. 13CO is present in the environments of both stars along with cold, dense gas as traced by CS and (in LkHa 225S) N2H+. No 2.6 mm continuum is detected in either region in relatively shallow maps compared to existing continuum observations.

We have used the integral field spectrograph OSIRIS and laser guide star adaptive optics at Keck Observatory to obtain high angular resolution (0.06"), moderate spectral resolution (R ~ 3800) images of the bipolar jet from the Herbig Ae star LkHa 233, seen in near-IR [Fe II] emission at 1.600 & 1.644 microns. This jet is narrow and tightly collimated, with an opening angle of only 9 degrees, and has an average radial velocity of ~ 100 km/s. The jet and counterjet are asymmetric, with the red-shifted jet much clumpier than its counterpart at the angular resolution of our observations. The observed properties are in general similar to jets seen around T Tauri stars, though it has a relatively large mass flux of (1.2e-7 +- 0.3e-7) M_sun/year, near the high end of the observed mass flux range around T Tauri stars. We also spatially resolve an inclined circumstellar disk around LkHa 233, which obscures the star from direct view. By comparison with numerical radiative transfer disk models, we estimate the disk midplane to be inclined i = 65 +- 5 degrees relative to the plane of the sky. Since the star is seen only in scattered light at near-infrared wavelengths, we detect only a small fraction of its intrinsic flux. Because previous estimates of its stellar properties did not account for this, either LkHa 233 must be located closer than the previously believed, or its true luminosity must be greater than previously supposed, consistent with its being a ~4 M_sun star near the stellar birthline.

The major obstacle to the direct detection of companions to nearby stars is the overwhelming brightness of the host star. Current instruments employing the combination of adaptive optics (AO) and coronagraphy can typically detect objects within 2'' of the star that are 10^{4-5} times fainter. Correlated speckle noise is one of the biggest obstacles limiting such high-contrast imaging. We have obtained a series of 284 8 s, AO-corrected, coronagraphically occulted H-band images of the star Vega at the 3.63 m AEOS telescope located on Haleakala, Hawaii. This dataset is unique for studying the temporal behavior of speckle noise, and represents the first time such a study on highly corrected _coronagraphic_ AO images has been carried out in a quantitative way. We find the speckle pattern to be highly stable in both position and time in our data. This is due to the fact that the AO system corrects disturbances to the stellar wave front at the level where the instrumental wave front errors dominate. Because of this, we find that our detection limit is not significantly improved simply with increased exposure time alone. However, we are able to improve our dynamic range by 1.5-2 magnitudes through subtraction of static/quasi-static speckles in two rotating frames: the telescope pupil frame and the deformable mirror frame. Furthermore, from our data, we are able to constrain the mass of any purported companion to Vega to be less than 45 M_J at 8 AU and less than 30 M_J at 16 AU, radii not previously probed at these sensitivities.

We have discovered an optically thick, edge-on circumstellar disk around a Herbig Ae star in the binary system PDS 144, providing the first intermediate-mass analog of HK Tau and similar T Tauris. This system consists of a V ~ 13 mag. primary and a fainter companion, with spectra of both stars showing evidence for circumstellar disks and accretion; both stars were classified as Herbig Aes by the Pico dos Dias survey. In Lick adaptive optics polarimetry, we resolved extended polarized light scattered from dust around the northern star. Followup Keck adaptive optics and mid-infrared observations show that this star is entirely hidden by an optically thick disk at all wavelengths from 1.2 to 11.7 microns. The disk major axis subtends ~ 0.8" on the sky, corresponding to ~ 800 AU at a distance of 1000 pc. Bright "wings" extend 0.3" above and below the disk ansae, due most likely to scattering from the edges of an outflow cavity in a circumstellar envelope. We discuss the morphology of the disk and the spectral energy distributions of the two PDS 144 stars, present preliminary disk models, and identify a number of open questions regarding this fascinating system.

Adaptive optics (AO) systems have significantly improved astronomical imaging capabilities over the last decade, and are revolutionizing the kinds of science possible with 4-5m class ground-based telescopes. A thorough understanding of AO system performance at the telescope can enable new frontiers of science as observations push AO systems to their performance limits. We look at recent advances with wave front reconstruction (WFR) on the Advanced Electro-Optical System (AEOS) 3.6 m telescope to show how progress made in improving WFR can be measured directly in improved science images. We describe how a "waffle mode" wave front error (which is not sensed by a Fried geometry Shack-Hartmann wave front sensor) affects the AO point-spread function (PSF). We model details of AEOS AO to simulate a PSF which matches the actual AO PSF in the I-band, and show that while the older observed AEOS PSF contained several times more waffle error than expected, improved WFR techniques noticeably improve AEOS AO performance. We estimate the impact of these improved WFRs on H-band imaging at AEOS, chosen based on the optimization of the Lyot Project near-infrared coronagraph at this bandpass.

We have used laser guide star adaptive optics and a near-infrared dual-channel imaging polarimeter to observe light scattered in the circumstellar environment of Herbig Ae/Be stars on scales of 100-300 AU. We discover a strongly polarized, biconical nebula 10 arcseconds in diameter (6000 AU) around the star LkHa 198, and also observe a polarized jet-like feature associated with the deeply embedded source LkHa 198-IR. The star LkHa 233 presents a narrow, unpolarized dark lane consistent with an optically thick circumstellar disk blocking our direct view of the star. These data show that the lower-mass T Tauri and intermediate mass Herbig Ae/Be stars share a common evolutionary sequence.

We describe the symmetries present in the point-spread function (PSF) of an optical system either located in space or corrected by an adaptive o to Strehl ratios of about 70% and higher. We present a formalism for expanding the PSF to arbitrary order in terms of powers of the Fourier transform of the residual phase error, over an arbitrarily shaped and apodized entrance aperture. For traditional unapodized apertures at high Strehl ratios, bright speckles pinned to the bright Airy rings are part of an antisymmetric perturbation of the perfect PSF, arising from the term that is first order in the residual phase error. There are two symmetric second degree terms. One is negative at the center, and, like the first order term, is modulated by the perfect image's field strength -- it reduces to the Marechal approximation at the center of the PSF. The other is non-negative everywhere, zero at the image center, and can be responsible for an extended halo -- which limits the dynamic range of faint companion detection in the darkest portions of the image. In regimes where one or the other term dominates the speckles in an image, the symmetry of the dominant term can be exploited to reduce the effect of those speckles, potentially by an order of magnitude or more. We demonstrate the effects of both secondary obscuration and pupil apodization on the structure of residual speckles, and discuss how these symmetries can be exploited by appropriate telescope and instrument design, observing strategies, and filter bandwidths to improve the dynamic range of high dynamic range AO and space-based observations. Finally, we show that our analysis is relevant to high dynamic range coronagraphy.

We present results from the flight of two prototype CZT detectors on a scientific balloon payload in September 2000. The first detector, referred to as ``CZT1,'' consisted of a 10 mm x 10 mm x 2 mm CZT crystal with a single gold planar electrode readout. This detector was shielded by a combination of a passive collimator surrounded by plastic scintillator and a thick BGO crystal in the rear. The second detector, ``CZT2,'' comprised two 10 mm x 10 mm x 5 mm CZT crystals, one made of eV Products high pressure Bridgman material and the other of IMARAD horizontal Bridgman material, each fashioned with a 4 x 4 array of gold pixels on a 2.5 mm pitch. The pixellated detectors were flip-chip-mounted side by side and read out by a 32-channel ASIC. This detector was also shielded by a passive/plastic collimator in the front, but used only additional passive/plastic shielding in the rear. Both experiments were flown from Ft. Sumner, NM on September 19, 2000 on a 24 hour balloon flight. CZT1 recorded a non-vetoed background level at 100 keV of ~1e-3 cts/cm2/s/keV. Raising the BGO threshold from 50 keV to ~1 MeV produced only an 18% increase in this level. CZT2 recorded a background at 100 keV of ~4e-3 cts/cm2/s/keV in the eV Products detector and ~6e-3 cts/cm2/s/keV in the IMARAD detector. Both CZT1 and CZT2 spectra were in basic agreement with Monte Carlo simulations, though both recorded systematically higher count rates at high energy than predicted. No lines were observed, indicating that neutron capture reactions, at least those producing decay lines at a few 100 keV, are not significant components of the CZT background. Comparison of the CZT1 and CZT2 spectra indicates that passive/plastic shielding may provide adequately low background levels for many applications.

Low resolution spectra of 23 stars towards the bulge globular clusters NGC 6528 and NGC 6553 are analysed. Radial velocities and atmospheric parameters are derived in order to check their membership in the clusters. Effective temperatures were obtained from photometric data for stars with Teff > 3800 K, whereas for cooler stars, they were derived from equivalent widths of TiO bands. Calibrations of W(TiO) as a function of stellar parameters based on a grid of synthetic spectra are presented. Metallicities were derived from a comparison of the observed spectra to a grid of synthetic spectra. The sample comprises evolutionary stages from the Red Giant Branch to the Horizontal Branch, with parameters in the range 3200 <= Teff <= 5000 K and -0.5 <= log g <= 2.4. The mean metallicities obtained for NGC 6528 and NGC 6553 are [Fe/H] ~ -0.5 and -0.7, in both cases with [Mg/Fe] = +0.3; assuming the same overabundance for the $\alpha$ elements O, Mg, Si, S, Ca and Ti, this gives [Z/Zsolar] = -0.25 and -0.45. Membership verification by means of low resolution spectra is a crucial step in preparing targets for high resolution spectroscopy with 8m class telescopes.

From an initial sample of a few hundred stars selected for an Hipparcos program, a subsample of 33 objects with direct determination of their effective temperature, bolometric magnitude and [Fe/H] has been derived. An accurate observational (logTeff, Mbol) HR diagram has been obtained, and compared to theoretical isochrones. Attention has been concentrated on the cool part of the HR diagram, for which evolutionary effects are negligible. Whereas solar metallicity stars show a reasonable agreement with theoretical expectation, a clear discrepancy occurs for stars more metal-poor than [Fe/H] = -0.5. This discrepancy is reduced to practically zero, if corrections to LTE abundances are taken into account for Fe, and if sedimentation effects are included in stellar evolution computations.