We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models.

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

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

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

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

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

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

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

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

An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay, and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low towards the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations.

The James Webb Space Telescope will offer high-angular resolution observing capability in the near-infrared with masking interferometry on NIRISS, and coronagraphic imaging on NIRCam & MIRI. Full aperture kernel-phase based interferometry complements these observing modes, probing for companions at small separations while preserving the telescope throughput. Our goal is to derive both theoretical and operational contrast detection limits for the kernel-phase analysis of JWST NIRISS full-pupil observations by using tools from hypothesis testing theory, applied to observations of faint brown dwarfs with this instrument, but the tools and methods introduced here are applicable in a wide variety of contexts. We construct a statistically independent set of observables from aberration-robust kernel phases. Three detection tests based on these observable quantities are designed and analysed, all guaranteeing a constant false alarm rate for small phase aberrations. One of these tests, the Likelihood Ratio or Neyman-Pearson test, provides a theoretical performance bound for any detection test. The operational detection method considered here is shown to exhibit only marginal power loss with respect to the theoretical bound. In principle, for the test set to a false alarm probability of 1%, companion at contrasts reaching 10^3 at separations of 200 mas around objects of magnitude 14.1 are detectable. With JWST NIRISS, contrasts of up to 10^4 at separations of 200 mas could be ultimately achieved, barring significant wavefront drift. The proposed detection method is close to the ultimate bound and offers guarantees over the probability of making a false detection for binaries, as well as over the error bars for the estimated parameters of the binaries detectable by JWST NIRISS. This method is not only applicable to JWST NIRISS but to any imaging system with adequate sampling.

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.

A coronagraphic starlight suppression system situated on a future flagship space observatory offers a promising avenue to image Earth-like exoplanets and search for biomarkers in their atmospheric spectra. One NASA mission concept that could serve as the platform to realize this scientific breakthrough is the Large UV/Optical/IR Surveyor (LUVOIR). Such a mission would also address a broad range of topics in astrophysics with a multiwavelength suite of instruments. The apodized pupil Lyot coronagraph (APLC) is one of several coronagraph design families that the community is assessing as part of NASAs Exoplanet Exploration Program Segmented aperture coronagraph design and analysis (SCDA) team. The APLC is a Lyot-style coronagraph that suppresses starlight through a series of amplitude operations on the on-axis field. Given a suite of seven plausible segmented telescope apertures, we have developed an object-oriented software toolkit to automate the exploration of thousands of APLC design parameter combinations. This has enabled us to empirically establish relationships between planet throughput and telescope aperture geometry, inner working angle, bandwidth, and contrast level. In parallel with the parameter space exploration, we have investigated several strategies to improve the robustness of APLC designs to fabrication and alignment errors. We also investigate the combination of APLC with wavefront control or complex focal plane masks to improve inner working angle and throughput. Preliminary scientific yield evaluations based on design reference mission simulations indicate the APLC is a very competitive concept for surveying the local exoEarth population with a mission like LUVOIR.

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.

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 imaging and spectroscopy of habitable worlds will require large-aperture space-based telescopes, to increase the collecting area and the angular resolution. These large telescopes will necessarily use segmented primaries to fit in a rocket. However, these massively segmented mirrors make high-contrast performance very difficult to achieve and stabilize, compared to more common monolithic primaries. Despite space telescopes operating in a friendlier environment than ground-based telescopes, remaining vibrations and resonant modes on the segments can still deteriorate the performance. In this context, we present the Pair-based Analytical model for Segmented Telescopes Imaging from Space (PASTIS) that enables the establishment of a comprehensive error budget, both in term of segment alignment and stability. Using this model, one may evaluate the influence of the segment cophasing and surface quality evolution on the final images and contrasts, and set up requirements for any given mission. One can also identify the dominant modes of a given geometry for a given coronagraphic instrument and design the feedback control systems accordingly. In this paper, we first develop and validate this analytical model by comparing its outputs to the images and contrasts predicted by an end-to-end simulation. We show that the contrasts predicted using PASTIS are accurate enough compared to the end-to-end propagation results, at the exo-Earth detection level. Second, we develop a method for a fast and efficient error budget in term of segment manufacturing and alignment that takes into account the disparities of the segment effects on the final performance. This technique is then applied on a specific aperture to provide static and quasi-static requirements on each segment for local aberrations. Finally we discuss potential application of this new technique to future missions.

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

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

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.

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

We present the first near infra-red spectrum of the faint white dwarf companion HD 114174 B, obtained with Project 1640. Our spectrum, covering the Y, J and H bands, combined with previous TRENDS photometry measurements, allows us to place further constraints on this companion. We suggest two possible scenarios; either this object is an old, low mass, cool H atmosphere white dwarf with Teff ~ 3800 K or a high mass white dwarf with Teff > 6000 K, potentially with an associated cool (Teff ~ 700 K) brown dwarf or debris disk resulting in an infra-red excess in the L band. We also provide an additional astrometry point for 2014 June 12 and use the modelled companion mass combined with the RV and direct imaging data to place constraints on the orbital parameters for this companion.

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.

The nearby Sun-like star GJ 758 hosts a cold substellar companion, GJ 758 B, at a projected separation of $\lesssim$30 AU, previously detected in high-contrast multi-band photometric observations. In order to better constrain the companion's physical characteristics, we acquired the first low-resolution ($R \sim 50$) near-infrared spectrum of it using the high-contrast hyperspectral imaging instrument Project 1640 on Palomar Observatory's 5-m Hale telescope. We obtained simultaneous images in 32 wavelength channels covering the $Y$, $J$, and $H$ bands ($\sim$952-1770 nm), and used data processing techniques based on principal component analysis to efficiently subtract chromatic background speckle-noise. GJ 758 B was detected in four epochs during 2013 and 2014. Basic astrometric measurements confirm its apparent northwest trajectory relative to the primary star, with no clear signs of orbital curvature. Spectra of SpeX/IRTF observed T dwarfs were compared to the combined spectrum of GJ 758 B, with ${\chi}^2$ minimization suggesting a best fit for spectral type T7.0$\pm$1.0, but with a shallow minimum over T5-T8. Fitting of synthetic spectra from the BT-Settl13 model atmospheres gives an effective temperature $T_{\text{eff}}=741 \pm 25$ K and surface gravity $\log g = 4.3 \pm 0.5$ dex (cgs). Our derived best-fit spectral type and effective temperature from modeling of the low-resolution spectrum suggest a slightly earlier and hotter companion than previous findings from photometric data, but do not rule out current results, and confirm GJ 758 B as one of the coolest sub-stellar companions to a Sun-like star to date.

We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.

Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.

We present new spatially resolved astrometry and photometry from the Gemini Planet Imager of the inner binary of the young multiple star system V343 Normae, which is a member of the beta Pictoris moving group. V343 Normae comprises a K0 and mid-M star in a ~4.5 year orbit (AaAb) and a wide 10" M5 companion (B). By combining these data with archival astrometry and radial velocities we fit the orbit and measure individual masses for both components of M_Aa = 1.10 +/- 0.10 M_sun and M_Ab = 0.290 +/- 0.018 M_sun. Comparing to theoretical isochrones, we find good agreement for the measured masses and JHK band magnitudes of the two components consistent with the age of the beta Pic moving group. We derive a model-dependent age for the beta Pic moving group of 26 +/- 3 Myr by combining our results for V343 Normae with literature measurements for GJ 3305, which is another group member with resolved binary components and dynamical masses.

We present $H$-band near-infrared polarimetric imaging observations of the F5V star HD~157587 obtained with the Gemini Planet Imager (GPI) that reveal the debris disk as a bright ring structure at a separation of $\sim$80$-$100~AU. The new GPI data complement recent HST/STIS observations that show the disk extending out to over 500~AU. The GPI image displays a strong asymmetry along the projected minor axis as well as a fainter asymmetry along the projected major axis. We associate the minor and major axis asymmetries with polarized forward scattering and a possible stellocentric offset, respectively. To constrain the disk geometry we fit two separate disk models to the polarized image, each using a different scattering phase function. Both models favor a disk inclination of $\sim 70\degr$ and a $1.5\pm0.6$ AU stellar offset in the plane of the sky along the projected major axis of the disk. We find that the stellar offset in the disk plane, perpendicular to the projected major axis is degenerate with the form of the scattering phase function and remains poorly constrained. The disk is not recovered in total intensity due in part to strong adaptive optics residuals, but we recover three point sources. Considering the system's proximity to the galactic plane and the point sources' positions relative to the disk, we consider it likely that they are background objects and unrelated to the disk's offset from the star.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo\`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $\beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$\sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $\beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$\sigma$ significance. However, we confirm that the Hill sphere of $\beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

Wavefront estimation using in-focus image data is critical to many applications. This data is invariant to a sign flip with complex conjugation of the complex amplitude in the pupil, making for a non-unique solution. Information from an in-focus image taken through a non-redundant pupil mask (NRM) can break this ambiguity, enabling the true aberration to be determined. We demonstrate this by priming a full pupil Gerchberg-Saxton phase retrieval with NRM fringe phase information. We apply our method to measure simulated aberrations on the segmented James Webb Space Telescope (JWST) mirror using full pupil and NRM data from its Near Infrared Imager and Slitless Spectrograph (NIRISS).

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.

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.

$\mu$ Her is a nearby quadruple system with a G-subgiant primary and several low mass companions arranged in a 2+2 architecture. While the BC components have been well characterized, the Ab component has been detected astrometrically and with direct imaging but there has been some confusion over its nature, in particular whether the companion is stellar or substellar. Using near-infrared spectroscopy we are able to estimate the spectral type of the companion as a M4$\pm$1V star. In addition, we have measured the astrometry of the system for over a decade. We combined the astrometry with archival radial velocity measurements to compute an orbit of the system. From the combined orbit, we are able to compute the mass sum of the system. Using the estimated mass of the primary, we estimate the mass of the secondary as 0.32 M_sun, which agrees with the estimated spectral type. Our computed orbit is preliminary due to the incomplete orbital phase coverage, but it should be sufficient to predict ephemerides over the next decade.

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 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 James Webb Space Telescope (JWST) will allow observations with a unique combination of spectral, spatial, and temporal resolution for the study of outer planet satellites within our Solar System. We highlight the infrared spectroscopy of icy moons and temporal changes on geologically active satellites as two particularly valuable avenues of scientific inquiry. While some care must be taken to avoid saturation issues, JWST has observation modes that should provide excellent infrared data for such studies.

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.

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

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

We present $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.

HD 177830 is an evolved K0IV star with two known exoplanets. In addition to the planetary companions it has a late-type stellar companion discovered with adaptive optics imagery. We observed the binary star system with the PHARO near-IR camera and the Project 1640 coronagraph. Using the Project 1640 coronagraph and integral field spectrograph we extracted a spectrum of the stellar companion. This allowed us to determine that the spectral type of the stellar companion is a M4$\pm$1V. We used both instruments to measure the astrometry of the binary system. Combining these data with published data, we determined that the binary star has a likely period of approximately 800 years with a semi-major axis of 100-200 AU. This implies that the stellar companion has had little or no impact on the dynamics of the exoplanets. The astrometry of the system should continue to be monitored, but due to the slow nature of the system, observations can be made once every 5-10 years.

HD 19467 B is presently the only directly imaged T dwarf companion known to induce a measurable Doppler acceleration around a solar type star. We present spectroscopy measurements of this important benchmark object taken with the Project 1640 integral field unit at Palomar Observatory. Our high-contrast R~30 observations obtained simultaneously across the $JH$ bands confirm the cold nature of the companion as reported from the discovery article and determine its spectral type for the first time. Fitting the measured spectral energy distribution to SpeX/IRTF T dwarf standards and synthetic spectra from BT-Settl atmospheric models, we find that HD 19467 B is a T5.5+/-1 dwarf with effective temperature Teff=$978^{+20}_{-43}$ K. Our observations reveal significant methane absorption affirming its substellar nature. HD 19467 B shows promise to become the first T dwarf that simultaneously reveals its mass, age, and metallicity independent from the spectrum of light that it emits.

The high angular resolution technique of non-redundant masking (NRM) or aperture masking interferometry (AMI) has yielded images of faint protoplanetary companions of nearby stars from the ground. AMI on James Webb Space Telescope (JWST)'s Near Infrared Imager and Slitless Spectrograph (NIRISS) has a lower thermal background than ground-based facilites and does not suffer from atmospheric instability. NIRISS AMI images are likely to have 90 - 95% Strehl ratio between 2.77 and 4.8 micron. In this paper we quantify factors that limit the raw point source contrast of JWST NRM. We develop an analytic model of the NRM point spread function which includes different optical path delays (pistons) between mask holes and fit the model parameters with image plane data. It enables a straightforward way to exclude bad pixels, is suited to limited fields of view, and can incorporate effects such as intra-pixel sensitivity variations. We simulate various sources of noise to estimate their effect on the standard deviation of closure phase, sigma_CP (a proxy for binary point source contrast). If sigma_CP < 10^-4 radians --- a contrast ratio of 10 magnitudes --- young accreting gas giant planets (e.g. in the nearby Taurus star-forming region) could be imaged with JWST NIRISS. We show the feasibility of using NIRISS' NRM with the sub-Nyquist sampled F277W, which would enable some exoplanet chemistry characterization. In the presence of small piston errors, the dominant sources of closure phase error (depending on pixel sampling, and filter bandwidth) are flat field errors and unmodeled variations in intra-pixel sensitivity. The in-flight stability of NIRISS will determine how well these errors can be calibrated by observing a point source. Our results help develop efficient observing strategies for space-based NRM.

We present an analysis of the orbital motion of the four sub-stellar objects orbiting HR8799. Our study relies on the published astrometric history of this system augmented with an epoch obtained with the Project 1640 coronagraph + Integral Field Spectrograph (IFS) installed at the Palomar Hale telescope. We first focus on the intricacies associated with astrometric estimation using the combination of an Extreme Adaptive Optics system (PALM-3000), a coronagraph and an IFS. We introduce two new algorithms. The first one retrieves the stellar focal plane position when the star is occulted by a coronagraphic stop. The second one yields precise astrometric and spectro-photometric estimates of faint point sources even when they are initially buried in the speckle noise. The second part of our paper is devoted to studying orbital motion in this system. In order to complement the orbital architectures discussed in the literature, we determine an ensemble of likely Keplerian orbits for HR8799bcde, using a Bayesian analysis with maximally vague priors regarding the overall configuration of the system. While the astrometric history is currently too scarce to formally rule out coplanarity, HR8799d appears to be misaligned with respect to the most likely planes of HR8799bce orbits. This misalignment is sufficient to question the strictly coplanar assumption made by various authors when identifying a Laplace resonance as a potential architecture. Finally, we establish a high likelihood that HR8799de have dynamical masses below 13 M_Jup using a loose dynamical survival argument based on geometric close encounters. We illustrate how future dynamical analyses will further constrain dynamical masses in the entire system.

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.

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.

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

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

The Gemini Planet Imager (GPI) 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.

JWST's Near-Infrared Imager and Slitless Spectrograph (NIRISS) includes an Aperture Masking Interferometry (AMI) mode designed to be used between 2.7{\mu}m and 4.8{\mu}m. At these wavelengths, it will have the highest angular resolution of any mode on JWST, and, for faint targets, of any existing or planned infrastructure. NIRISS AMI is uniquely suited to detect thermal emission of young massive planets and will permit the characterization of the mid-IR flux of exoplanets discovered by the GPI and SPHERE adaptive optics surveys. It will also directly detect massive planets found by GAIA through astrometric accelerations, providing the first opportunity ever to get both a mass and a flux measurement for non-transiting giant planets. NIRISS AMI will also enable the study of the nuclear environment of AGNs.

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.

Due to feedback from accretion onto supermassive black holes (SMBHs), Active Galactic Nuclei (AGNs) are believed to play a key role in LambdaCDM cosmology and galaxy formation. However, AGNs' extreme luminosities and the small angular size of their accretion flows create a challenging imaging problem. We show James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) Aperture Masking Interferometry (AMI) mode will enable true imaging (i.e. without any requirement of prior assumptions on source geometry) at ~65 mas angular resolution at the centers of AGNs. This is advantageous for studying complex extended accretion flows around SMBHs, and in other areas of angular-resolution-limited astrophysics. By simulating data sequences incorporating expected sources of noise, we demonstrate that JWST-NIRISS AMI mode can map extended structure at a pixel-to-pixel contrast of ~10^{-2} around an L=7.5 point source, using short exposure times (minutes). Such images will test models of AGN feedback, fuelling and structure (complementary with ALMA observations), and are not currently supported by any ground-based IR interferometer or telescope. Binary point source contrast with NIRISS is ~10^{-4} (for observing binary nuclei in merging galaxies), significantly better than current ground-based optical or IR interferometry. JWST-NIRISS' seven-hole non-redundant mask has a throughput of 15%, and utilizes NIRISS' F277W (2.77\micron), F380M (3.8\micron), F430M (4.3\micron), and F480M (4.8\micron) filters. NIRISS' square pixels are 65mas per side, with a field of view ~2\arcmin x 2\arcmin. We also extrapolate our results to AGN science enabled by non-redundant masking on future 2.4m and 16m space telescopes working at long-UV to near-IR wavelengths.

The Project 1640 instrument on the 200-inch Hale telescope at Palomar Observatory is a coronagraphic instrument with an integral field spectrograph at the back end, designed to find young, self-luminous planets around nearby stars. To reach the necessary contrast for this, the PALM-3000 adaptive optics system corrects for fast atmospheric speckles, while CAL, a phase-shifting interferometer in a Mach-Zehnder configuration, measures the quasistatic components of the complex electric field in the pupil plane following the coronagraphic stop. Two additional sensors measure and control low-order modes. These field measurements may then be combined with a system model and data taken separately using a white-light source internal to the AO system to correct for both phase and amplitude aberrations. Here, we discuss and demonstrate the procedure to maintain a half-plane dark hole in the image plane while the spectrograph is taking data, including initial on-sky performance.

Kappa Andromedae is a B9IVn star at 52 pc for which a faint substellar companion separated by 55 AU was recently announced. In this work, we present the first spectrum of the companion, "kappa And B," using the Project 1640 high-contrast imaging platform. Comparison of our low-resolution YJH-band spectra to empirical brown dwarf spectra suggests an early-L spectral type. Fitting synthetic spectra from PHOENIX model atmospheres to our observed spectrum allows us to constrain the effective temperature to ~2000K, as well as place constraints on the companion surface gravity. Further, we use previously reported log(g) and effective temperature measurements of the host star to argue that the kappa And system has an isochronal age of 220 +/- 100 Myr, older than the 30 Myr age reported previously. This interpretation of an older age is corroborated by the photometric properties of kappa And B, which appear to be marginally inconsistent with other 10-100 Myr low-gravity L-dwarfs for the spectral type range we derive. In addition, we use Keck aperture masking interferometry combined with published radial velocity measurements to rule out the existence of any tight stellar companions to kappa And A that might be responsible for the system's overluminosity. Further, we show that luminosity enhancements due to a nearly "pole-on" viewing angle coupled with extremely rapid rotation is unlikely. Kappa And A is thus consistent with its slightly evolved luminosity class (IV) and we propose here that kappa And, with a revised age of 220 +/- 100 Myr, is an interloper to the 30 Myr Columba association with which it was previously associated. The photometric and spectroscopic evidence for kappa And B combined with our re-assesment of the system age implies a substellar companion mass of 50^{+16}_{-13} Jupiter Masses, consistent with a brown dwarf rather than a planetary mass companion.

We obtained spectra, in the wavelength range \lambda = 995 - 1769 nm, of all four known planets orbiting the star HR 8799. Using the suite of instrumentation known as Project 1640 on the Palomar 5-m Hale Telescope, we acquired data at two epochs. This allowed for multiple imaging detections of the companions and multiple extractions of low-resolution (R ~ 35) spectra. Data reduction employed two different methods of speckle suppression and spectrum extraction, both yielding results that agree. The spectra do not directly correspond to those of any known objects, although similarities with L and T-dwarfs are present, as well as some characteristics similar to planets such as Saturn. We tentatively identify the presence of CH_4 along with NH_3 and/or C_2H_2, and possibly CO_2 or HCN in varying amounts in each component of the system. Other studies suggested red colors for these faint companions, and our data confirm those observations. Cloudy models, based on previous photometric observations, may provide the best explanation for the new data presented here. Notable in our data is that these presumably co-eval objects of similar luminosity have significantly different spectra; the diversity of planets may be greater than previously thought. The techniques and methods employed in this paper represent a new capability to observe and rapidly characterize exoplanetary systems in a routine manner over a broad range of planet masses and separations. These are the first simultaneous spectroscopic observations of multiple planets in a planetary system other than our own.

We present adaptive optics photometry and spectra in the JHKL-bands along with high spectral resolution K-band spectroscopy for each component of the Z Canis Majoris system. Our high angular resolution photometry of this very young (<1 Myr) binary, comprised of an FU Ori object and a Herbig Ae/Be star, were gathered shortly after the 2008 outburst while our high resolution spectroscopy was gathered during a quiescent phase. Our photometry conclusively determine that the outburst was due solely to the embedded Herbig Ae/Be member, supporting results from earlier works, and that the optically visible FU Ori component decreased slightly (~30%) in luminosity during the same period, consistent with previous works on the variability of FU Ori type systems. Further, our high-resolution K-band spectra definitively demonstrate that the 2.294 micron CO absorption feature seen in composite spectra of the system is due solely to the FU Ori component, while a prominent CO emission feature at the same wavelength, long suspected to be associated with the innermost regions of a circumstellar accretion disk, can be assigned to the Herbig Ae/Be member. These findings are in contrast to previous analyses (e.g. Malbet et al 2010, Benisty et al. 2010) of this complex system which assigned the CO emission to the FU Ori component.

We used the Project 1640 near-infrared coronagraph and integral field spectrograph to observe 19 young solar type stars. Five of these stars are known binary stars and we detected the late-type secondaries and were able to measure their JH spectra with a resolution of R\sim30. The reduced, extracted, and calibrated spectra were compared to template spectra from the IRTF spectral library. With this comparison we test the accuracy and consistency of spectral type determination with the low-resolution near-infrared spectra from P1640. Additionally, we determine effective temperature and surface gravity of the companions by fitting synthetic spectra calculated with the PHOENIX model atmosphere code. We also present several new epochs of astrometry of each of the systems. Together these data increase our knowledge and understanding of the stellar make up of these systems. In addition to the astronomical results, the analysis presented helps validate the Project 1640 data reduction and spectral extraction processes and the utility of low-resolution, near-infrared spectra for characterizing late-type companions in multiple systems.

We describe a new instrument that forms the core of a long-term high contrast imaging program at the 200-inch Hale Telescope at Palomar Observatory. The primary scientific thrust is to obtain images and low-resolution spectroscopy of brown dwarfs and young Jovian mass exoplanets in the vicinity of stars within 50 parsecs of the Sun. The instrument is a microlens-based integral field spectrograph integrated with a diffraction limited, apodized-pupil Lyot coronagraph, mounted behind the Palomar adaptive optics system. The spectrograph obtains imaging in 23 channels across the J and H bands (1.06 - 1.78 microns). In addition to obtaining spectra, this wavelength resolution allows suppression of the chromatically dependent speckle noise, which we describe. We have recently installed a novel internal wave front calibration system that will provide continuous updates to the AO system every 0.5 - 1.0 minutes by sensing the wave front within the coronagraph. The Palomar AO system is undergoing an upgrade to a much higher-order AO system ("PALM-3000"): a 3388-actuator tweeter deformable mirror working together with the existing 241-actuator mirror. This system will allow correction with subapertures as small as 8cm at the telescope pupil using natural guide stars. The coronagraph alone has achieved an initial dynamic range in the H-band of 2 X 10^-4 at 1 arcsecond, without speckle noise suppression. We demonstrate that spectral speckle suppression is providing a factor of 10-20 improvement over this bringing our current contrast at an arcsecond to ~2 X 10^-5. This system is the first of a new generation of apodized pupil coronagraphs combined with high-order adaptive optics and integral field spectrographs (e.g. GPI, SPHERE, HiCIAO), and we anticipate this instrument will make a lasting contribution to high contrast imaging in the Northern Hemisphere for years.

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

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

The A5V star Alcor has an M3-M4 dwarf companion, as evidenced by a novel astrometric technique. Imaging spectroscopy combined with adaptive optics coronagraphy allowed for the detection and spectrophotometric characterization of the point source at a contrast of ~6 J- and H-band magnitudes and separation of 1" from the primary star. The use of an astrometric pupil plane grid allowed us to determine the projected separations between the companion and the coronagraphically occulted primary star to <=3 milliarcsecond precision at two observation epochs. Our measurements demonstrate common parallactic and proper motion over the course of 103 days, significantly shorter than the period of time needed for most companion confirmations through proper motion measurements alone. This common parallax method is potentially more rigorous than common proper motion, ensuring that the neighboring bodies lie at the same distance, rather than relying on the statistical improbability that two objects in close proximity to each other on the sky move in the same direction. The discovery of a low-mass (~0.25M_sun) companion around a bright (V = 4.0), nearby (d = 25 pc) star highlights a region of binary star parameter space that to date has not been fully probed.

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.

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.

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 present a new method for numerical propagation through Lyot-style coronagraphs using finite occulting masks. Standard methods for coronagraphic simulations involve Fast Fourier Transforms (FFT) of very large arrays, and computing power is an issue for the design and tolerancing of coronagraphs on segmented Extremely Large Telescopes (ELT) in order to handle both the speed and memory requirements. Our method combines a semi-analytical approach with non-FFT based Fourier transform algorithms. It enables both fast and memory-efficient computations without introducing any additional approximations. Typical speed improvements based on computation costs are of about ten to fifty for propagations from pupil to Lyot plane, with thirty to sixty times less memory needed. Our method makes it possible to perform numerical coronagraphic studies even in the case of ELTs using a contemporary commercial laptop computer, or any standard commercial workstation computer.

The Gemini Planet (GPI) imager is an "extreme" adaptive optics system being designed and built for the Gemini Observatory. GPI combines precise and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with integral field and polarimetry capabilities. GPI's primary science goal is the direct detection and characterization of young, Jovian-mass exoplanets. For systems younger than 2 Gyr exoplanets more massive than 6 MJ and semimajor axes beyond 10 AU are detected with completeness greater than 50%. GPI will also discover faint debris disks, explore icy moons and minor planets in the solar system, reveal high dynamic range main-sequence binaries, and study mass loss from evolved stars. This white paper explains the role of GPI in exoplanet discovery and characterization and summarizes our recommendations to the NSF-NASA-DOE Astronomy and Astrophysics Advisory Committee ExoPlanet Task Force.

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 propose a solution to the problem of astrometric and photometric calibration of coronagraphic images with a simple optical device which, in theory, is easy to use. Our design uses the Fraunhofer approximation of Fourier optics. Placing a periodic grid of wires (we use a square grid) with known width and spacing in a pupil plane in front of the occulting coronagraphic focal plane mask produces fiducial images of the obscured star at known locations relative to the star. We also derive the intensity of these fiducial images in the coronagraphic image. These calibrator images can be used for precise relative astrometry, to establish companionship of other objects in the field of view through measurement of common proper motion or common parallax, to determine orbits, and to observe disk structure around the star quantitatively. The calibrator spots also have known brightness, selectable by the coronagraph designer, permitting accurate relative photometry in the coronagraphic image. This technique, which enables precision exoplanetary science, is relevant to future coronagraphic instruments, and is particularly useful for `extreme' adaptive optics and space-based coronagraphy.

We study the way Lyot coronagraphs with unapodized entrance pupils respond to small, low order phase aberrations. This study is applicable to ground-based adaptive optics coronagraphs operating at 90% and higher Strehl ratios, as well as to some space-based coronagraphs with intrinsically higher Strehl ratio imaging. We utilize a second order expansion of the monochromatic point-spread function (written as a power spectrum of a power series in the phase aberration over clear aperture) to derive analytical expressions for the response of a `band-limited' Lyot coronagraph (BLC) to small, low order, phase aberrations. The BLC possesses a focal plane mask with an occulting spot whose opacity profile is a spatially band-limited function rather than a hard-edged, opaque disk. The BLC is, to first order, insensitive to tilt and astigmatism. Undersizing the stop in the re-imaged pupil plane (the Lyot plane) following the focal plane mask can alleviate second order effects of astigmatism, at the expense of system throughput and angular resolution. The optimal degree of such undersizing depends on individual instrument designs and goals. Our analytical work engenders physical insight, and complements existing numerical work on this subject. Our methods can be extended to treat the passage of higher order aberrations through band-limited Lyot coronagraphs, by using our polynomial decomposition or an analogous Fourier approach.

We study the way Lyot coronagraphs with unapodized entrance pupils respond to small, low order phase aberrations. This study is applicable to ground-based adaptive optics coronagraphs operating at 90% and higher Strehl ratios, as well as to some space-based coronagraphs with intrinsically higher Strehl ratio imaging. We utilize a second order expansion of the monochromatic point-spread function (written as a power spectrum of a power series in the phase aberration over clear aperture) to derive analytical expressions for the response of a `band-limited' Lyot coronagraph (BLC) to small, low order, phase aberrations. The BLC possesses a focal plane mask with an occulting spot whose opacity profile is a spatially band-limited function rather than a hard-edged, opaque disk. The BLC is, to first order, insensitive to tilt and astigmatism. Undersizing the stop in the re-imaged pupil plane (the Lyot plane) following the focal plane mask can alleviate second order effects of astigmatism, at the expense of system throughput and angular resolution. The optimal degree of such undersizing depends on individual instrument designs and goals. Our analytical work engenders physical insight, and complements existing numerical work on this subject. Our methods can be extended to treat the passage of higher order aberrations through band-limited Lyot coronagraphs, by using our polynomial decomposition or an analogous Fourier approach.

In principle, suppression of on-axis stellar light by a coronagraph is easier on an unobscured aperture telescope than on one with an obscured aperture. Recent designs such as the apodized pupil Lyot coronagraph, the `band-limited' Lyot coronagraph, and several variants of phase mask coronagraphs work best on unobscured circular aperture telescopes. These designs are developed to enable the discovery and characterization of nearby jovian or even terrestrial exoplanets. All of today's major space-based and adaptive optics-equipped ground based telescopes are obscured aperture systems, with a secondary mirror held in place by secondary support `spider' vanes. The presence of a secondary obscuration can be dealt with by ingenious coronagraph designs, but the spider vanes themselves cause diffracted light that can hamper the search for jovian exoplanets around nearby stars. We look at the problem of suppressing spider vane diffraction in Lyot coronagraphs, including apodized pupil and band-limited designs. We show how spider vane diffraction can be reduced drastically, and in fact contained, in the final coronagraphic image, within one resolution element of the geometric image of the focal plane mask's occulting spot. This makes adaptive optics coronagraphic searches for exojupiters possible on the next generation of adaptive optics systems being developed for 8-10 m class telescopes such as Gemini and the Very Large Telescopes.

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 present a study of Lyot style (i.e., classical, band-limited, and Gaussian occulter) coronagraphy on extremely large, highly-segmented telescopes. We show that while increased telescope diameter is always an advantage for high dynamic range science (assuming wavefront errors have been corrected sufficiently well), segmentation itself sets a limit on the performance of Lyot coronagraphs. Diffraction from inter-segment gaps sets a floor to the achievable extinction of on-axis starlight with Lyot coronagraphy. We derive an analytical expression for the manner in which coronagraphic suppression of an on-axis source decreases with increasing gap size when the segments are placed in a spatially periodic array over the telescope aperture, regardless of the details of the arrangement. A simple Lyot stop masking out pupil edges produces good extinction of the central peak in the point-spread function (PSF), but leaves satellite images caused by inter-segment gaps essentially unaffected. Masking out the bright segment gaps in the Lyot plane with a reticulated mask reduces the satellite images'intensity to a contrast of 5x10^{-9} on a 30 m telescope with 10 mm gaps, at the expense of an increase in the brightness of the central peak. The morphology of interesting targets will dictate which Lyot stop geometry is preferable: the reticulated Lyot stop produces a conveniently uni-modal PSF, whereas a simple Lyot stop produces an extended array of satellite spots. A cryogenic reticulate Lyot stop will also benefit both direct and coronagraphic mid-IR imaging.

The direct detection of extrasolar planets by imaging means is limited by the large flux of light from the host star being scattered into the region of interest by a variety of processes, including diffraction. Coronagraphs are devices that suppress the undesirable scattering of light caused by diffraction. In a coronagraph the sensitivity limit for high dynamic range is limited by the propagation of errors introduced by the imperfect optical system to the final image. In this paper we develop theory and simulations to understand how such errors propagate in a coronagraph. We describe the response of classical and band-limited Lyot coronagraphs to small and large errors in the placement of the central star, and identify ways of making such coronagraphs more robust to small guiding errors. We also uncover features of the decentered PSF that can lead to spurious detection of companions, especially with aggressive, high dynamic range coronagraphs dedicated to companion searches aimed at finding extrasolar terrestrial or Jovian planets.

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.

The useful dynamic range of an image in the diffraction limited regime is usually limited by speckles caused by residual phase errors in the optical system forming the image. The technique of speckle decorrelation involves introducing many independent realizations of additional phase error into a wavefront during one speckle lifetime, changing the instantaneous speckle pattern. A commonly held assumption is that this results in the speckles being `moved around' at the rate at which the additional phase screens are applied. The intention of this exercise is to smooth the speckles out into a more uniform background distribution during their persistence time, thereby enabling companion detection around bright stars to be photon noise limited rather than speckle-limited. We demonstrate analytically why this does not occur, and confirm this result with numerical simulations. We show that the original speckles must persist, and that the technique of speckle decorrelation merely adds more noise to the original speckle noise, thereby degrading the dynamic range of the image.

We summarize the theory of coronagraphic optics, and identify a dimensionless fine-tuning parameter, F, which we use to describe the Lyot stop size in the natural units of the coronagraphic optical train and the observing wavelength. We then present simulations of coronagraphs matched to adaptive optics (AO) systems on the Calypso 1.2m, Palomar Hale 5m and Gemini 8m telescopes under various atmospheric conditions, and identify useful parameter ranges for AO coronagraphy on these telescopes. Our simulations employ a tapered, high-pass filter in spatial frequency space to mimic the action of adaptive wavefront correction. We test the validity of this representation of AO correction by comparing our simulations with recent K-band data from the 241-channel Palomar Hale AO system and its dedicated PHARO science camera in coronagraphic mode.