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

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

Direct imaging of Earth-like planets from space requires dedicated observatories, combining large segmented apertures with instruments and techniques such as coronagraphs, wavefront sensors, and wavefront control in order to reach the high contrast of 10^10 that is required. The complexity of these systems would be increased by the segmentation of the primary mirror, which allows for the larger diameters necessary to image Earth-like planets but also introduces specific patterns in the image due to the pupil shape and segmentation and making high-contrast imaging more challenging. Among these defects, the phasing errors of the primary mirror are a strong limitation to the performance. In this paper, we focus on the wavefront sensing of segment phasing errors for a high-contrast system, using the COronagraphic Focal plane wave-Front Estimation for Exoplanet detection (COFFEE) technique. We implemented and tested COFFEE on the High-contrast imaging for Complex Aperture Telescopes (HiCAT) testbed, in a configuration without any coronagraph and with a classical Lyot coronagraph, to reconstruct errors applied on a 37 segment mirror. We analysed the quality and limitations of the reconstructions. We demonstrate that COFFEE is able to estimate correctly the phasing errors of a segmented telescope for piston, tip, and tilt aberrations of typically 100nm RMS. We also identified the limitations of COFFEE for the reconstruction of low-order wavefront modes, which are highly filtered by the coronagraph. This is illustrated using two focal plane mask sizes on HiCAT. We discuss possible solutions, both in the hardware system and in the COFFEE optimizer, to mitigate these issues.

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.

Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

Space coronagraphs are projected to detect exoplantes that are at least 1e10 times dimmer than their host stars. Yet, the actual detection threshold depends on the instrument's wavefront stability and varies by an order of magnitude with the choice of observation strategy and post-processing method. In this paper the authors consider the performance of the previously introduced observation strategy (dark hole maintenance) and post-processing algorithm (electric field order reduction) in the presence of various realistic effects. In particular, it will be shown that under some common assumptions, the telescope's averaged pointing jitter translates into an additional light source incoherent with the residual light from the star (speckles), and that jitter "modes" can be identified in post-processing and distinguished from a planet signal. We also show that the decrease in contrast due to drift of voltages in deformable mirror actuators can be mitigated by recursive estimation of the electric field in the high-contrast region of the image (dark hole) using Electric Field Conjugation (EFC). Moreover, this can be done even when the measured intensity is broadband, as long as it is well approximated by an incoherent sum of monochromatic intensities. Finally, we assess the performance of closed-loop vs. open-loop observation scenarios through a numerical simulation of the Wide-Field Infra-Red Survey Telescope (WFIRST). In particular, we compare the post-processing factors of Angular Differential Imaging (ADI) with and without Electric Field Order Reduction (EFOR), which we extended to account for possible telescope rolls and the presence of pointing jitter. For all observation parameters considered in this paper, close-loop dark hole maintenance resulted in significantly higher post-processing accuracy.

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.

Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

We present new results concerning the sub-stellar binary population in the Orion Nebula Cluster (ONC). Using the Karhunen-Lo\`{e}ve Image Projection (KLIP) algorithm, we have reprocessed images taken with the IR channel of the Wide Field Camera 3 mounted on the Hubble Space Telescope to unveil faint close companions in the wings of the stellar PSFs. Starting with a sample of $1392$ bona-fide not saturated cluster members, we detect $39$ close-pairs cluster candidates with separation $0.16''-0.77''$. The primary masses span a range M$_p$ $\sim 0.015-1.27$ M$_{\odot}$ whereas for the companions we derive M$_c$ $\sim 0.004-0.54$ M$_{\odot}$. Of these $39$ binary systems, $18$ were already known while the remaining $21$ are new detections. Correcting for completeness and combining our catalog with previously detected ONC binaries, we obtain an overall binary fraction of $11.5\% \pm 0.9\%$. Compared to other star forming regions, our multiplicity function is $\sim 2$ smaller than e.g. Taurus, while compared to the binaries in the field we obtain comparable values. We analyze the mass function of the binaries, finding differences between the mass distribution of binaries and single stars and between primary and companion mass distributions. The mass ratio shows a bottom-heavy distribution with median value of $M_c/M_p \sim 0.25$. Overall our results suggest that ONC binaries may represent a template for the typical population of field binaries, supporting the hypothesis that the ONC may be regarded as a most typical star forming region in the Milky Way.

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

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

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

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

Further advances in exoplanet detection and characterisation require sampling a diverse population of extrasolar planets. One technique to detect these distant worlds is through the direct detection of their thermal emission. The so-called direct imaging technique, is suitable for observing young planets far from their star. These are very low signal-to-noise-ratio (SNR) measurements and limited ground truth hinders the use of supervised learning approaches. In this paper, we combine deep generative and discriminative models to bypass the issues arising when directly training on real data. We use a Generative Adversarial Network to obtain a suitable dataset for training Convolutional Neural Network classifiers to detect and locate planets across a wide range of SNRs. Tested on artificial data, our detectors exhibit good predictive performance and robustness across SNRs. To demonstrate the limits of the detectors, we provide maps of the precision and recall of the model per pixel of the input image. On real data, the models can re-confirm bright source detections.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

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

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

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

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

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

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

Alpha Centauri A is the closest solar-type star to the Sun and offers an excellent opportunity to detect the thermal emission of a mature planet heated by its host star. The MIRI coronagraph on JWST can search the 1-3 AU (1"-2") region around alpha Cen A which is predicted to be stable within the alpha Cen AB system. We demonstrate that with reasonable performance of the telescope and instrument, a 20 hr program combining on-target and reference star observations at 15.5 um could detect thermal emission from planets as small as ~5 RE. Multiple visits every 3-6 months would increase the geometrical completeness, provide astrometric confirmation of detected sources, and push the radius limit down to ~3 RE. An exozodiacal cloud only a few times brighter than our own should also be detectable, although a sufficiently bright cloud might obscure any planet present in the system. While current precision radial velocity (PRV) observations set a limit of 50-100 ME at 1-3 AU for planets orbiting alpha Cen A, there is a broad range of exoplanet radii up to 10 RE consistent with these mass limits. A carefully planned observing sequence along with state-of-the-art post-processing analysis could reject the light from alpha Cen A at the level of ~10^-5 at 1"-2" and minimize the influence of alpha Cen B located 7-8" away in the 2022-2023 timeframe. These space-based observations would complement on-going imaging experiments at shorter wavelengths as well as PRV and astrometric experiments to detect planets dynamically. Planetary demographics suggest that the likelihood of directly imaging a planet whose mass and orbit are consistent with present PRV limits is small, ~5%, and possibly lower if the presence of a binary companion further reduces occurrence rates. However, at a distance of just 1.34 pc, alpha Cen A is our closest sibling star and certainly merits close scrutiny.

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.

Direct imaging of exo-Earths and search for life is one of the most exciting and challenging objectives for future space observatories. Segmented apertures in space will be required to reach the needed large diameters beyond the capabilities of current or planned launch vehicles. These apertures present additional challenges for high-contrast coronagraphy, not only in terms of static phasing but also in terms of their stability. The Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS) was developed to model the effects of segment-level optical aberrations on the final image contrast. In this paper, we extend the original PASTIS propagation model from a purely analytical to a semi-analytical method, in which we substitute the use of analytical images with numerically simulated images. The inversion of this model yields a set of orthonormal modes that can be used to determine segment-level wavefront tolerances. We present results in the case of segment-level piston error applied to the baseline coronagraph design of LUVOIR A, with minimum and maximum wavefront error constraint between 56 pm and 290 pm per segment. The analysis is readily generalizable to other segment-level aberrations modes, and can also be expanded to establish stability tolerances for these missions.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

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

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

The expected yield of potentially Earth-like planets is a useful metric for designing future exoplanet-imaging missions. Recent yield studies of direct-imaging missions have focused primarily on yield methods and trade studies using "toy" models of missions. Here we increase the fidelity of these calculations substantially, adopting more realistic exoplanet demographics as input, an improved target list, and a realistic distribution of exozodi levels. Most importantly, we define standardized inputs for instrument simulations, use these standards to directly compare the performance of realistic instrument designs, include the sensitivity of coronagraph contrast to stellar diameter, and adopt engineering-based throughputs and detector parameters. We apply these new high-fidelity yield models to study several critical design trades: monolithic vs segmented primary mirrors, on-axis vs off-axis secondary mirrors, and coronagraphs vs starshades. We show that as long as the gap size between segments is sufficiently small, there is no difference in yield for coronagraph-based missions with monolithic off-axis telescopes and segmented off-axis telescopes, assuming that the requisite engineering constraints imposed by the coronagraph can be met in both scenarios. We show that there is currently a factor of ~2 yield penalty for coronagraph-based missions with on-axis telescopes compared to off-axis telescopes, and note that there is room for improvement in coronagraph designs for on-axis telescopes. We also reproduce previous results in higher fidelity showing that the yields of coronagraph-based missions continue to increase with aperture size while the yields of starshade-based missions turnover at large apertures if refueling is not possible. Finally, we provide absolute yield numbers with uncertainties that include all major sources of astrophysical noise to guide future mission design.

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.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

Exoplanet imaging and spectroscopy are now routinely achieved by dedicated instruments on large ground-based observatories (e.g. Gemini/GPI, VLT/SPHERE, or Subaru/SCExAO). In addition to extreme adaptive optics (ExAO) and post-processing methods, these facilities make use of the most advanced coronagraphs to suppress light of an observed star and enable the observation of circumstellar environments. The Apodized Pupil Lyot Coronagraph (APLC) is one of the leading coronagraphic baseline in the current generation of instruments. This concept combines a pupil apodization, an opaque focal plane mask (FPM), and a Lyot stop. APLC can be optimized for a range of applications and designs exist for on-axis segmented aperture telescopes at $10^{10}$ contrast in broadband light. In this communication, we propose novel designs to push the limits of this concept further by modifying the nature of the FPM from its standard opaque mask to a smaller size occulting spot surrounded by circular phase shifting zones. We present the formalism of this new concept which solutions find two possible applications: 1) upgrades for the current generation of ExAO coronagraphs since these solutions remain compatible with the existing designs and will provide better inner working angle, contrast and throughput, and 2) coronagraphy at $10^{10}$ contrast for future flagship missions such as LUVOIR, with the goal to increase the throughput of the existing designs for the observation of Earth-like planets around nearby stars.

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

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

Rigorously quantifying the information in high contrast imaging data is important for informing follow-up strategies to confirm the substellar nature of a point source, constraining theoretical models of planet-disk interactions, and deriving planet occurrence rates. However, within the exoplanet direct imaging community, non-detections have almost exclusively been defined using a frequentist detection threshold (i.e. contrast curve) and associated completeness. This can lead to conceptual inconsistencies when included in a Bayesian framework. A Bayesian upper limit is such that the true value of a parameter lies below this limit with a certain probability. The associated probability is the integral of the posterior distribution with the upper limit as the upper bound. In summary, a frequentist upper limit is a statement about the detectability of planets while a Bayesian upper limit is a statement about the probability of a parameter to lie in an interval given the data. The latter is therefore better suited for rejecting hypotheses or theoretical models based on their predictions. In this work we emphasize that Bayesian statistics and upper limits are more easily interpreted and typically more constraining than the frequentist approach. We illustrate the use of Bayesian analysis in two different cases: 1) with a known planet location where we also propose to use model comparison to constrain the astrophysical nature of the point source and 2) gap-carving planets in TW Hya. To finish, we also mention the problem of combining radial velocity and direct imaging observations.

We present the results of a deep high-contrast imaging search for planets around Vega. Vega is an ideal target for high-contrast imaging because it is bright, nearby, and young with a face-on two-belt debris disk which may be shaped by unseen planets. We obtained $J-$ and $H-$band data on Vega with the coronagraphic integral-field spectrograph Project 1640 (P1640) at Palomar Observatory. Two nights of data were obtained in 2016, in poor seeing conditions, and two additional nights in more favorable conditions in 2017. In total, we obtained 5.5 hours of integration time on Vega in moderate to good seeing conditions (<1.5"). We did not detect any low mass companions in this system. Our data present the most sensitive contrast limits around Vega at very small separations (2-15 au) thus far, allowing us to place new constraints on the companions which may be sculpting the Vega system. In addition to new constraints, as the deepest data obtained with P1640, these observations form the final legacy of the now decommissioned instrument.

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 Optimal Optical Coronagraph (OOC) Workshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this first installment of a series of three papers summarizing the outcomes of the OOC workshop, we present an overview of design methods and optical performance metrics developed for coronagraph instruments. The design and optimization of coronagraphs for future telescopes has progressed rapidly over the past several years in the context of space mission studies for Exo-C, WFIRST, HabEx, and LUVOIR as well as ground-based telescopes. Design tools have been developed at several institutions to optimize a variety of coronagraph mask types. We aim to give a broad overview of the approaches used, examples of their utility, and provide the optimization tools to the community. Though it is clear that the basic function of coronagraphs is to suppress starlight while maintaining light from off-axis sources, our community lacks a general set of standard performance metrics that apply to both detecting and characterizing exoplanets. The attendees of the OOC workshop agreed that it would benefit our community to clearly define quantities for comparing the performance of coronagraph designs and systems. Therefore, we also present a set of metrics that may be applied to theoretical designs, testbeds, and deployed instruments. We show how these quantities may be used to easily relate the basic properties of the optical instrument to the detection significance of the given point source in the presence of realistic noise.

The Optimal Optical Coronagraph (OOC) Workshop held at the Lorentz Center in September 2017 in Leiden, the Netherlands, gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this second installment of a series of three papers summarizing the outcomes of the OOC workshop (see also~\citenum{ruane2018,snik2018}), we present an overview of common path wavefront sensing/control and Coherent Differential Imaging techniques, highlight the latest results, and expose their relative strengths and weaknesses. We layout critical milestones for the field with the aim of enhancing future ground/space based high contrast imaging platforms. Techniques like these will help to bridge the daunting contrast gap required to image a terrestrial planet in the zone where it can retain liquid water, in reflected light around a G type star from space.

The Optimal Optical CoronagraphWorkshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. This contribution is the final part of a series of three papers summarizing the outcomes of the workshop, and presents an overview of novel optical technologies and systems that are implemented or considered for high-contrast imaging instruments on both ground-based and space telescopes. The overall objective of high contrast instruments is to provide direct observations and characterizations of exoplanets at contrast levels as extreme as 10^-10. We list shortcomings of current technologies, and identify opportunities and development paths for new technologies that enable quantum leaps in performance. Specifically, we discuss the design and manufacturing of key components like advanced deformable mirrors and coronagraphic optics, and their amalgamation in "adaptive coronagraph" systems. Moreover, we discuss highly integrated system designs that combine contrast-enhancing techniques and characterization techniques (like high-resolution spectroscopy) while minimizing the overall complexity. Finally, we explore extreme implementations using all-photonics solutions for ground-based telescopes and dedicated huge apertures for space telescopes.

The Gemini Planet Imager (GPI) is the dedicated high-contrast imaging facility, located on Gemini South, designed for the direct detection and characterization of young Jupiter mass exoplanets. In 2019, Gemini is considering moving GPI from Gemini South to Gemini North. Analysis of GPI's as-built performance has highlighted several key areas of improvement to its detection capabilities while leveraging its current capabilities as a facility class instrument. We present the proposed upgrades which include a pyramid wavefront sensor, broadband low spectral resolution prisms and new apodized-pupil Lyot coronagraph designs all of which will enhance the current science capabilities while enabling new science programs.

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.

Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and co-rotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10-yr timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 $L'$-band observations in 2017. MWC~758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of ${0.6^\circ}^{+3.3^\circ}_{-0.6^\circ}\, \mathrm{yr}^{-1}$ at $3\sigma$ for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best fit orbit is at $89$ au ($0.59"$), with a $3\sigma$ lower limit of 30 au ($0.20"$). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.

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

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

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

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

Pandeia is the exposure time calculator (ETC) system developed for the James Webb Space Telescope (JWST) that will be used for creating JWST proposals. It includes a simulation-hybrid Python engine that calculates the two-dimensional pixel-by-pixel signal and noise properties of the JWST instruments. This allows for appropriate handling of realistic point spread functions, MULTIACCUM detector readouts, correlated detector readnoise, and multiple photometric and spectral extraction strategies. Pandeia includes support for all the JWST observing modes, including imaging, slitted/slitless spectroscopy, integral field spectroscopy, and coronagraphy. Its highly modular, data-driven design makes it easily adaptable to other observatories. An implementation for use with WFIRST is also available.

BANYAN $\Sigma$ is a new Bayesian algorithm to identify members of young stellar associations within 150 pc of the Sun. It includes 27 young associations with ages in the range ~1-800 Myr, modelled with multivariate Gaussians in 6-dimensional XYZUVW space. It is the first such multi-association classification tool to include the nearest sub-groups of the Sco-Cen OB star-forming region, the IC 2602, IC 2391, Pleiades and Platais 8 clusters, and the $\rho$ Ophiuchi, Corona Australis, and Taurus star-formation regions. A model of field stars is built from a mixture of multivariate Gaussians based on the Besan\c{c}on Galactic model. The algorithm can derive membership probabilities for objects with only sky coordinates and proper motion, but can also include parallax and radial velocity measurements, as well as spectrophotometric distance constraints from sequences in color-magnitude or spectral type-magnitude diagrams. BANYAN $\Sigma$ benefits from an analytical solution to the Bayesian marginalization integrals that makes it more accurate and significantly faster than its predecessor BANYAN II. A contamination versus hit rate analysis is presented and demonstrates that BANYAN $\Sigma$ achieves a better classification performance than other moving group classification tools, especially in terms of cross-contamination between young associations. An updated list of bona fide members in the 27 young associations, augmented by the Gaia-DR1 release, are presented. This new tool will make it possible to analyze large data sets such as the upcoming Gaia-DR2 to identify new young stars. IDL and Python versions of BANYAN $\Sigma$ are made available with this publication. (shortened)

The physical properties of brown dwarf companions found to orbit nearby, solar-type stars can be benchmarked against independent measures of their mass, age, chemical composition, and other parameters, offering insights into the evolution of substellar objects. The TRENDS high-contrast imaging survey has recently discovered a (mass/age/metallicity) benchmark brown dwarf orbiting the nearby (d=18.69+/-0.19 pc), G8V/K0V star HD 4747. We have acquired follow-up spectroscopic measurements of HD 4747 B using the Gemini Planet Imager to study its spectral type, effective temperature, surface gravity, and cloud properties. Observations obtained in the H-band and K1-band recover the companion and reveal that it is near the L/T transition (T1+/-2). Fitting atmospheric models to the companion spectrum, we find strong evidence for the presence of clouds. However, spectral models cannot satisfactorily fit the complete data set: while the shape of the spectrum can be well-matched in individual filters, a joint fit across the full passband results in discrepancies that are a consequence of the inherent color of the brown dwarf. We also find a $2\sigma$ tension in the companion mass, age, and surface gravity when comparing to evolutionary models. These results highlight the importance of using benchmark objects to study "secondary effects" such as metallicity, non-equilibrium chemistry, cloud parameters, electron conduction, non-adiabatic cooling, and other subtleties affecting emergent spectra. As a new L/T transition benchmark, HD 4747 B warrants further investigation into the modeling of cloud physics using higher resolution spectroscopy across a broader range of wavelengths, polarimetric observations, and continued Doppler radial velocity and astrometric monitoring.

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

We apply the vectorized Non-negative Matrix Factorization (NMF) method to post-processing of direct imaging data for exoplanetary systems such as circumstellar disks. NMF is an iterative approach, which first creates a non-orthogonal and non-negative basis of components using given reference images, then models a target with the components. The constructed model is then rescaled with a factor to compensate for the contribution from a disk. We compare NMF with existing methods (classical reference differential imaging method, and the Karhunen-Lo\`eve image projection algorithm) using synthetic circumstellar disks, and demonstrate the superiority of NMF: with no need for prior selection of references, NMF can detect fainter circumstellar disks, better preserve low order disk morphology, and does not require forward modeling. As an application to a well-known disk example, we process the archival Hubble Space Telescope (HST) STIS coronagraphic observations of HD~181327 with different methods and compare them. NMF is able to extract some circumstellar material inside the primary ring for the first time. In the appendix, we mathematically investigate the stability of NMF components during iteration, and the linearity of NMF modeling.

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

The increasing complexity of the aperture geometry of the future space- and ground based-telescopes will limit the performance of the next generation of coronagraphic instruments for high contrast imaging of exoplanets. We propose here a new closed-loop optimization technique using two deformable mirrors to correct for the effects of complex apertures on coronagraph performance, alternative to the ACAD technique previously developed by our group. This technique, ACAD-OSM, allows the use of any coronagraphs designed for continuous apertures, with complex, segmented, apertures, maintaining high performance in contrast and throughput. We show the capabilities of this technique on several pupil geometries (segmented LUVOIR type aperture, WFIRST, ELTs) for which we obtained high contrast levels with several deformable mirror setups (size, number of actuators, separation between them), coronagraphs (apodized pupil Lyot and vortex coronagraphs) and spectral bandwidths, which will help us present recommendations for future coronagraphic instruments. We show that this active technique handles, without any revision to the algorithm, changing or unknown optical aberrations or discontinuities in the pupil, including optical design misalignments, missing segments and phase errors.

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

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

We present two state-of-the-art models of the solar system, one corresponding to the present day and one to the Archean Eon 3.5 billion years ago. Each model contains spatial and spectral information for the star, the planets, and the interplanetary dust, extending to 50 AU from the sun and covering the wavelength range 0.3 to 2.5 micron. In addition, we created a spectral image cube representative of the astronomical backgrounds that will be seen behind deep observations of extrasolar planetary systems, including galaxies and Milky Way stars. These models are intended as inputs to high-fidelity simulations of direct observations of exoplanetary systems using telescopes equipped with high-contrast capability. They will help improve the realism of observation and instrument parameters that are required inputs to statistical observatory yield calculations, as well as guide development of post-processing algorithms for telescopes capable of directly imaging Earth-like planets.

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

The current generation of ground-based coronagraphic instruments uses deformable mirrors to correct for phase errors and to improve contrast levels at small angular separations. Improving these techniques, several space and ground based instruments are currently developed using two deformable mirrors to correct for both phase and amplitude errors. However, as wavefront control techniques improve, more complex telescope pupil geome- tries (support structures, segmentation) will soon be a limiting factor for these next generation coronagraphic instruments. In this paper we discuss fundamental limits associated with wavefront control with deformable mirrors in high contrast coronagraph. We start with an analytic prescription of wavefront errors, along with their wave- length dependence, and propagate them through coronagraph models. We then consider a few wavefront control architectures, number of deformable mirrors and their placement in the optical train of the instrument, and algorithms that can be used to cancel the starlight scattered by these wavefront errors over a finite bandpass. For each configuration we derive the residual contrast as a function of bandwidth and of the properties of the incoming wavefront. This result has consequences when setting the wavefront requirements, along with the wavefront control architecture of future high contrast instrument both from the ground and from space. In particular we show that these limits can severely affect the effective Outer Working Angle that can be achieved by a given coronagraph instrument.

The current generation of ground-based coronagraphic instruments uses deformable mirrors to correct for phase errors and to improve contrast levels at small angular separations. Improving these techniques, several space and ground based instruments are currently developed using two deformable mirrors to correct for both phase and amplitude errors. However, as wavefront control techniques improve, more complex telescope pupil geometries will soon be a limiting factor for these next generation coronagraphic instru- ments. The technique presented in this proceeding, the Active Correction of Aperture Discontinuities method, is taking advantage of the fact that most future coronagraphic instruments will include two deformable mirrors, and is proposing to find the shapes and actuator movements to correct for the effect introduced by these complex pupil geometries. We here present the results of the parametric analysis realized on the WFIRST pupil for which we obtained high contrast levels with several deformable mirror setups (size, separation between them), coronagraphs (Vortex charge 2, vortex charge 4, APLC) and spectral bandwidths. However, because contrast levels and separation are not the only metrics to maximize the scientific return of an instrument, we also included in this study the influence of these deformable mirror shapes on the throughput of the instrument and sensitivity to pointing jitters. Finally, we present results obtained on another potential space based telescope segmented aperture. The main result of this proceeding is that we now obtain comparable performance than the coronagraphs previously designed for WFIRST. First result from the parametric analysis strongly suggest that the 2 deformable mirror set up (size and distance between them) has a important impact on the performance in contrast and throughput of the final instrument.

High-contrast imaging and spectroscopy provide unique constraints for exoplanet formation models as well as for planetary atmosphere models. Instrumentation techniques in this field have greatly improved over the last two decades, with the development of stellar coronagraphy, in parallel with specific methods of wavefront sensing and control. Next generation space- and ground-based telescopes will allow the characterization cold solar-system like planets for the first time and maybe even in situ detection of bio-markers. However, the growth of primary mirror diameters, necessary for these detection, comes with an increase of their complexity (segmentation, secondary mirror features). These discontinuities in the aperture can greatly limit the performance of coronagraphic instruments. In this context, we introduced a new technique, Active Correction of Aperture Discontinuities - Optimized Stroke Minimization (ACAD-OSM), to correct for the diffractive effects of aperture discontinuities in the final image plane of a coronagraph, using deformable mirrors. In this paper, we present several tools that can be used to optimize the performance of this technique for its application to future large missions. In particular, we analyze the influence of the deformable setup (size and separating distance) and found that there is an optimal point for this setup, optimizing the performance of the instrument in contrast and throughput while minimizing the strokes applied to the deformable mirrors. These results will help us design future coronagraphic instruments to obtain the best performance.

Several coronagraph designs have been proposed over the last two decades to directly image exoplanets. Among these designs, the vector vortex coronagraphs provide theoretically perfect starlight cancellation along with small inner working angles when deployed on telescopes with unobstructed pupils. However, current and planned space missions and ground-based extremely large telescopes present complex pupil geometries, including secondary mirror central obscurations, that prevent vortex coronagraphs from rejecting on-axis sources entirely. Recent solutions combining the vortex phase mask with a ring-apodized pupil have been proposed to circumvent this issue, but provide a limited throughput for vortex charges $>2$. We present a family of pupil plane apodizations that compensate for pupil geometries with circularly symmetric central obstructions caused by on-axis secondary mirrors for charge 2, 4, and 6 vector vortex coronagraphs. These apodizations are derived analytically and allow the vortex coronagraph to retain theoretically perfect nulling in the presence of central obscurations. For a charge 4 vortex, we design polynomial apodization functions assuming a greyscale apodizing filter that represent a substantial gain in throughput over the ring-apodized vortex coronagraph design, while for a charge 6 vortex, we design polynomial apodized vortex coronagraphs that have $\gtrsim 70\%$ total energy throughput for the entire range of central obscuration sizes studied. We propose methods for optimizing apodizations produced with either greyscale apodizing filters or shaped mirrors. We conclude by demonstrating how this design may be combined with apodizations numerically optimized for struts and segment gaps in telescope pupils to design terrestrial exoplanet imagers for complex pupils.

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

Future searches for biomarkers on habitable exoplanets will rely on telescope instruments that achieve extremely high contrast at small planet-to-star angular separations. Coronagraphy is a promising starlight suppression technique, providing excellent contrast and throughput for off-axis sources on clear apertures. However, the complexity of space- and ground-based telescope apertures goes on increasing over time, owing to the combination of primary mirror segmentation, secondary mirror, and support structures. These discontinuities in the telescope aperture limit the coronagraph performance. In this paper, we present ACAD-OSM, a novel active method to correct for the diffractive effects of aperture discontinuities in the final image plane of a coronagraph. Active methods use one or several deformable mirrors that are controlled with an interaction matrix to correct for the aberrations in the pupil. However, they are often limited by the amount of aberrations introduced by aperture discontinuities. This algorithm relies on the recalibration of the interaction matrix during the correction process to overcome this limitation. We first describe the ACAD-OSM technique and compare it to the previous active methods for the correction of aperture discontinuities. We then show its performance in terms of contrast and off-axis throughput for static aperture discontinuities (segmentation, struts) and for some aberrations evolving over the life of the instrument (residual phase aberrations, artifacts in the aperture, misalignments in the coronagraph design). This technique can now obtain the earth-like planet detection threshold of 10^(-10) contrast on any given aperture over at least a 10% spectral bandwidth, with several coronagraph designs.

Gravitational deflection of starlight around the Sun during the 1919 total solar eclipse provided measurements that confirmed Einstein's general theory of relativity. We have used the Hubble Space Telescope to measure the analogous process of astrometric microlensing caused by a nearby star, the white dwarf Stein 2051 B. As Stein 2051 B passed closely in front of a background star, the background star's position was deflected. Measurement of this deflection at multiple epochs allowed us to determine the mass of Stein 2051 B -- the sixth nearest white dwarf to the Sun -- as 0.675 +/- 0.051 solar masses. This mass determination provides confirmation of the physics of degenerate matter and lends support to white dwarf evolutionary theory.

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.

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

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 have obtained new images of the protoplanetary disk orbiting TW Hya in visible, total intensity light with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST), using the newly commissioned BAR5 occulter. These HST/STIS observations achieved an inner working angle $\sim$0.2\arcsec, or 11.7~AU, probing the system at angular radii coincident with recent images of the disk obtained by ALMA and in polarized intensity near-infrared light. By comparing our new STIS images to those taken with STIS in 2000 and with NICMOS in 1998, 2004, and 2005, we demonstrate that TW Hya's azimuthal surface brightness asymmetry moves coherently in position angle. Between 50~AU and 141~AU we measure a constant angular velocity in the azimuthal brightness asymmetry of 22.7$^\circ$~yr$^{-1}$ in a counter-clockwise direction, equivalent to a period of 15.9~yr assuming circular motion. Both the (short) inferred period and lack of radial dependence of the moving shadow pattern are inconsistent with Keplerian rotation at these disk radii. We hypothesize that the asymmetry arises from the fact that the disk interior to 1~AU is inclined and precessing due to a planetary companion, thus partially shadowing the outer disk. Further monitoring of this and other shadows on protoplanetary disks potentially opens a new avenue for indirectly observing the sites of planet formation.

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

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

We present the first resolved image of the debris disk around the 16+/-8 Myr old star, HD 114082. The observation was made in the H-band using the SPHERE instrument. The star is at a distance of 92+/-6 pc in the Lower Centaurus Crux association. Using a Markov Chain Monte Carlo analysis, we determined that the debris is likely in the form of a dust ring with an inner edge of 27.7+2.8/-3.5 au, position angle -74+0.5/-1.5 deg, and an inclination with respect to the line of sight of 6.7+3.8/-0.4 deg. The disk imaged in scattered light has a surface density declining with radius like ~r^(-4), steeper than expected for grain blowout by radiation pressure. We find only marginal evidence (2 sigma) of eccentricity, and rule out planets more massive than 1.0 Mjup orbiting within 1 au of the ring's inner edge, since such a planet would have disrupted the disk. The disk has roughly the same fractional disk luminosity (Ldisk/L*=3.3x10^(-3)) as HR4796A and Beta Pictoris, however it was not detected by previous instrument facilities most likely because of its small angular size (radius~0.4"), low albedo (~0.2) and low scattering efficiency far from the star due to high scattering anisotropy. With the arrival of extreme adaptive optics systems like SPHERE and GPI, the morphology of smaller, fainter and more distant debris disks are being revealed, providing clues to planet-disk interactions in young protoplanetary systems.

HR8799 is orbited by at least four giant planets, making it a prime target for the recently commissioned Spectro-Polarimetric High-contrast Exoplanet REsearch (VLT/SPHERE). As such, it was observed on five consecutive nights during the SPHERE science verification in December 2014. We aim to take full advantage of the SPHERE capabilities to derive accurate astrometric measurements based on H-band images acquired with the Infra-Red Dual-band Imaging and Spectroscopy (IRDIS) subsystem, and to explore the ultimate astrometric performance of SPHERE in this observing mode. We also aim to present a detailed analysis of the orbital parameters for the four planets. We report the astrometric positions for epoch 2014.93 with an accuracy down to 2.0 mas, mainly limited by the astrometric calibration of IRDIS. For each planet, we derive the posterior probability density functions for the six Keplerian elements and identify sets of highly probable orbits. For planet d, there is clear evidence for nonzero eccentricity ($e \simeq 0.35$), without completely excluding solutions with smaller eccentricities. The three other planets are consistent with circular orbits, although their probability distributions spread beyond $e = 0.2$, and show a peak at $e \simeq 0.1$ for planet e. The four planets have consistent inclinations of about $30\deg$ with respect to the sky plane, but the confidence intervals for the longitude of ascending node are disjoint for planets b and c, and we find tentative evidence for non-coplanarity between planets b and c at the $2 \sigma$ level.

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

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

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

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

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

Current state-of-the-art high contrast imaging instruments take advantage of a number of elegant coronagraph designs to suppress starlight and image nearby faint objects, such as exoplanets and circumstellar disks. The ideal performance and complexity of the optical systems depends strongly on the shape of the telescope aperture. Unfortunately, large primary mirrors tend to be segmented and have various obstructions, which limit the performance of most conventional coronagraph designs. We present a new family of vortex coronagraphs with numerically-optimized gray-scale apodizers that provide the sensitivity needed to directly image faint exoplanets with large, segmented aperture telescopes, including the Thirty Meter Telescope (TMT) as well as potential next-generation space telescopes.

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.

A new class of high-contrast image analysis algorithms that empirically fit and subtract systematic noise has lead to recent discoveries of faint exoplanet /substellar companions and scattered light images of circumstellar disks. These methods are extremely efficient at enhancing the detectability of faint astrophysical signal, but they do generally create systematic biases in their observed properties. This paper provides a general solution for this outstanding problem. We present the analytical derivation of a linear expansion that captures the impact of astrophysical over-subtraction and/or self-subtraction these image analysis techniques. We examine the general case for which the reference images of the astrophysical scene move azimuthally and/or radially across the field of view as a result of the observation strategy. Our new method is based on perturbing the covariance matrix underlying any least-squares speckles problem, and propagating this perturbation through the data analysis algorithm. We then demonstrate practical applications of this new algorithm. We first consider the case of the spectral extraction of faint point sources in IFS data and illustrate, using public Gemini Planet Imager commissioning data, that our novel perturbation-based Karhunen-Lo\`eve Image Processing Forward Modeling (KLIP-FM) can indeed alleviate algorithmic biases. We then apply KLIP-FM to the problem associated with the detection of point sources. We show how it decreases the rate of false negatives (e.g missed planets) while keeping the rate of false positives unchanged when compared to classical least-squares fitting methods. This can potentially have important consequences on the design of follow-up strategies of ongoing direct imaging surveys.

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

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

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

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

The new frontier in the quest for the highest contrast levels in the focal plane of a coronagraph is now the correction of the large diffractive artifacts effects introduced at the science camera by apertures of increasing complexity. The coronagraph for the WFIRST/AFTA mission will be the first of such instruments in space with a two Deformable Mirrors wavefront control system. Regardless of the control algorithm for these multi Deformable Mirrors, they will have to rely on quick and accurate simulation of the propagation effects introduced by the out-of-pupil surface. In the first part of this paper, we present the analytical description of the different approximations to simulate these propagation effects. In Annex A, we prove analytically that, in the special case of surfaces inducing a converging beam, the Fresnel method yields high fidelity for simulations of these effects. We provide numerical simulations showing this effect. In the second part, we use these tools in the framework of the Active Compensation of Aperture Discontinuities technique (ACAD) applied to pupil geometries similar to WFIRST-AFTA. We present these simulations in the context of the optical layout of the High-contrast imager for Complex Aperture Telescopes, which will test ACAD on a optical bench. The results of this analysis show that using the ACAD method, an apodized pupil lyot coronagraph and the performance of our current deformable mirrors, we are able to obtain, in numerically simulations, a dark hole with an AFTA-like pupil. Our numerical simulation shows that we can obtain contrast better than $2.10^{-9}$ in monochromatic light and better than 3.e-8 with 10% bandwidth between 5 and 14 lambda/D.

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

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

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

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

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

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

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

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

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

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

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

We present results from a Hubble Space Telescope (HST) program characterizing the atmospheres of the outer two planets, in the HR8799 system. The images were taken over 15 orbits in three near-infrared medium-band filters - F098M, F127M and F139M - using the Wide Field Camera 3. One of the three filters is sensitive to water absorption band inaccessible from ground-based observations, providing a unique probe of the thermal emission from the atmospheres of these young giant planets. The observations were taken at 30 different spacecraft rolls to enable angular differential imaging, and the full data set was analyzed with the Karhunen-Loeve Image Projection (KLIP) routine, an advanced image processing algorithm adapted to work with HST data. To achieve the required high contrast at sub arcsecond resolution, we utilized the pointing accuracy of HST in combination with an improved pipeline designed to combine the dithered, angular differential imaging data with an algorithm designed to both improve the image resolution and accurately measure the photometry. The results include F127M ($J$) detections of the outer planets, HR8799 b and c and the first detection of HR8799 b in the water-band (F139M) filter. The F127M photometry for HR8799 c agrees well with fitted atmospheric models resolving a long standing difficulty to model the near-IR flux for the planet consistently

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.

The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r~0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

HR 4796A is surrounded by a well-structured and very bright circumstellar disc shaped like an annulus with many interesting features: very sharp inner and outer edges, brightness asymmetries, centre offset, and suspected distortions in the ring. We aim to constrain the properties of the dust surrounding the star HR 4796A, in particular the grain size and composition. We also want to confirm and refine the morphological parameters derived from previous scattered light observations, and reveal the dust spatial extent in regions unexplored so far due to their proximity to the star. We have obtained new images in polarised light of the binary system HR 4796A and B in the Ks and Lp band with the NaCo instrument at the Very Large Telescope (VLT). In addition, we revisit two archival data sets obtained in the Lp band with that same instrument and at 2.2 microns with the NICMOS instrument on the Hubble Space Telescope. We analyse these observations with simulations using the radiative transfer code MCFOST to investigate the dust properties. We explore a grid of models with various dust compositions and sizes in a Bayesian approach. We detect the disc in polarised light in the Ks band and reveal for the first time the innermost regions down to $0.3"$ along the semi-minor axis. We measure a polarised fraction of $29\%\pm8\%$ in the two disc ansae, with a maximum occurring more than $13^\circ$ westwards from the ansae. A very pronounced brightness asymmetry between the north-west and south-east side is detected. This contradicts the asymmetry previously reported in all images of the disc in unpolarised light at wavelengths smaller than or equal to 2.2 microns and is inconsistent with the predicted scattered light from spherical grains using the Mie theory. Our modelling suggests the north-west side is most likely inclined towards the Earth, contrary to previous conclusions.

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 Apodized Pupil Lyot Coronagraph (APLC) is a diffraction suppression system installed in the recently deployed instruments Palomar/P1640, Gemini/GPI, and VLT/SPHERE to allow direct imaging and spectroscopy of circumstellar environments. Using a prolate apodization, the current implementations offer raw contrasts down to $10^{-7}$ at 0.2 arcsec from a star over a wide bandpass (20\%), in the presence of central obstruction and struts, enabling the study of young or massive gaseous planets. Observations of older or lighter companions at smaller separations would require improvements in terms of inner working angle (IWA) and contrast, but the methods originally used for these designs were not able to fully explore the parameter space. We here propose a novel approach to improve the APLC performance. Our method relies on the linear properties of the coronagraphic electric field with the apodization at any wavelength to develop numerical solutions producing coronagraphic star images with high-contrast region in broadband light. We explore the parameter space by considering different aperture geometries, contrast levels, dark-zone sizes, bandpasses, and focal plane mask sizes. We present an application of these solutions to the case of Gemini/GPI with a design delivering a $10^{-8}$ raw contrast at 0.19 arcsec and offering a significantly reduced sensitivity to low-order aberrations compared to the current implementation. Optimal solutions have also been found to reach $10^{-10}$ contrast in broadband light regardless of the telescope aperture shape (in particular the central obstruction size), with effective IWA in the $2-3.5\lambda/D$ range, therefore making the APLC a suitable option for the future exoplanet direct imagers on the ground or in space.

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 the first observations obtained with the L'-band AGPM vortex coronagraph recently installed on LBTI/LMIRCam. The AGPM (Annular Groove Phase Mask) is a vector vortex coronagraph made from diamond subwavelength gratings. It is designed to improve the sensitivity and dynamic range of high-resolution imaging at very small inner working angles, down to 0.09 arcseconds in the case of LBTI/LMIRCam in the L' band. During the first hours on sky, we observed the young A5V star HR\,8799 with the goal to demonstrate the AGPM performance and assess its relevance for the ongoing LBTI planet survey (LEECH). Preliminary analyses of the data reveal the four known planets clearly at high SNR and provide unprecedented sensitivity limits in the inner planetary system (down to the diffraction limit of 0.09 arcseconds).

The new planet finder for the Very Large Telescope (VLT), the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE), just had its first light in Paranal. A dedicated instrument for the direct detection of planets, SPHERE, is composed of a polametric camera in visible light, the Zurich IMager POLarimeter (ZIMPOL), and two near-infrared sub-systems: the Infra-Red Dual-beam Imager and Spectrograph (IRDIS), a multi-purpose camera for imaging, polarimetry, and long-slit spectroscopy, and the integral field spectrograph (IFS), an integral field spectrograph. We present the results obtained from the analysis of data taken during the laboratory integration and validation phase, after the injection of synthetic planets. Since no continuous field rotation could be performed in the laboratory, this analysis presents results obtained using reduction techniques that do not use the angular differential imaging (ADI) technique. To perform the simulations, we used the instrumental point spread function (PSF) and model spectra of L and T-type objects scaled in contrast with respect to the host star. We evaluated the expected error in astrometry and photometry as a function of the signal to noise of companions, after spectral differential imaging (SDI) reduction for IRDIS and spectral deconvolution (SD) or principal component analysis (PCA) data reductions for IFS. We deduced from our analysis, for example, that $\beta$Picb, a 12~Myr old planet of $\sim$10~\MJ and semi-major axis of 9--10 AU, would be detected with IRDIS with a photometric error of 0.16~mag and with a relative astrometric position error of 1.1~mas. With IFS, we could retrieve a spectrum with error bars of about 0.15~mag on each channel and astrometric relative position error of 0.6~mas. For a fainter object such as HR8799d, a 13~\MJ planet at a distance of 27~AU, IRDIS could obtain a relative astrometric error of 3~mas.

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.

High dynamic-range imagers aim to block out or null light from a very bright primary star to make it possible to detect and measure far fainter companions; in real systems a small fraction of the primary light is scattered, diffracted, and unocculted. We introduce S4, a flexible data-driven model for the unocculted (and highly speckled) light in the P1640 spectroscopic coronograph. The model uses Principal Components Analysis (PCA) to capture the spatial structure and wavelength dependence of the speckles but not the signal produced by any companion. Consequently, the residual typically includes the companion signal. The companion can thus be found by filtering this error signal with a fixed companion model. The approach is sensitive to companions that are of order a percent of the brightness of the speckles, or up to $10^{-7}$ times the brightness of the primary star. This outperforms existing methods by a factor of 2-3 and is close to the shot-noise physical limit.

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.

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

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

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

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

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

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

In this paper, we review the impact of small sample statistics on detection thresholds and corresponding confidence levels (CLs) in high contrast imaging at small angles. When looking close to the star, the number of resolution elements decreases rapidly towards small angles. This reduction of the number of degrees of freedom dramatically affects CLs and false alarm probabilities. Naively using the same ideal hypothesis and methods as for larger separations, which are well understood and commonly assume Gaussian noise, can yield up to one order of magnitude error in contrast estimations at fixed CL. The statistical penalty exponentially increases towards very small inner working angles. Even at 5-10 resolution elements from the star, false alarm probabilities can be significantly higher than expected. Here we present a rigorous statistical analysis which ensures robustness of the CL, but also imposes a substantial limitation on corresponding achievable detection limits (thus contrast) at small angles. This unavoidable fundamental statistical effect has a significant impact on current coronagraphic and future high contrast imagers. Finally, the paper concludes with practical recommendations to account for small number statistics when computing the sensitivity to companions at small angles and when exploiting the results of direct imaging planet surveys.

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

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

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

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

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

With a clear circular aperture, the vortex coronagraph perfectly cancels an on-axis point source and offers a 0.9 or 1.75 lambda/D inner working angle for topological charge 2 or 4, respectively. Current and near-future large telescopes are on-axis, however, and the diffraction effects of the central obscuration, and the secondary supports are strong enough to prevent the detection of companions 1e-3 - 1e-5 as bright as, or fainter than, their host star. Recent advances show that a ring apodizer can restore the performance of this coronagraph by compensating for the diffraction effects of a circular central obscuration in a 1D modeling of the pupil. We extend this work and optimize apodizers for arbitrary apertures in 2D in order to tackle the diffraction effects of the spiders and other noncircular artefacts in the pupil. We use a numerical optimization scheme to compute hybrid coronagraph designs that combine the advantages of the vortex coronagraph (small in IWA) and of shaped pupils coronagraphs (robustness to central obscuration and pupil asymmetric structures). We maximize the apodizer transmission, while constraints are set on the extremum values of the electric field that is computed in chosen regions of the Lyot plane through closed form expressions. Optimal apodizers are computed for topological charges 2 and 4 vortex coronagraphs and for telescope apertures with 10-30% central obscurations and 0-1% thick spiders. We characterize the impacts of the obscuration ratio and the thickness of the spiders on the throughput and the IWA for the two topological charges.

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.

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.

The vortex coronagraph (VC) is a new generation small inner working angle (IWA) coronagraph currently offered on various 8-meter class ground-based telescopes. On these observing platforms, the current level of performance is not limited by the intrinsic properties of actual vortex devices, but by wavefront control residuals and incoherent background (e.g. thermal emission of the sky) or the light diffracted by the imprint of the secondary mirror and support structures on the telescope pupil. In the particular case of unfriendly apertures (mainly large central obscuration) when very high contrast is needed (e.g. direct imaging of older exoplanets with extremely large telescopes or space- based coronagraphs), a simple VC, as most coronagraphs, can not deliver its nominal performance because of the contamination due to the diffraction from the obscured part of the pupil. Here we propose a novel yet simple concept that circumvents this problem. We combine a vortex phase mask in the image plane of a high-contrast instrument with a single pupil-based amplitude ring apodizer, tailor designed to exploit the unique convolution properties of the VC at the Lyot-stop plane. We show that such a ring-apodized vortex coronagraph (RAVC) restores the perfect attenuation property of the VC regardless of the size of the central obscuration, and for any (even) topological charge of the vortex. More importantly the RAVC maintains the IWA and conserves a fairly high throughput, which are signature properties of the VC.

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 report low resolution near infrared spectroscopic observations of the eruptive star FU Orionis using the Integral Field Spectrograph Project 1640 installed at the Palomar Hale telescope. This work focuses on elucidating the nature of the faint source, located 0.5" south of FU Ori, and identified in 2003 as FU Ori S. We first use our observations in conjunction with published data to demonstrate that the two stars are indeed physically associated and form a true binary pair. We then proceed to extract J and H band spectro-photometry using the damped LOCI algorithm, a reduction method tailored for high contrast science with IFS. This is the first communication reporting the high accuracy of this technique, pioneered by the Project 1640 team, on a faint astronomical source. We use our low resolution near infrared spectrum in conjunction with 10.2 micron interferometric data to constrain the infrared excess of FU Ori S. We then focus on estimating the bulk physical properties of FU Ori S. Our models lead to estimates of an object heavily reddened, A_V =8-12, with an effective temperature of ~ 4000-6500 K . Finally we put these results in the context of the FU Ori N-S system and argue that our analysis provides evidence that FU Ori S might be the more massive component of this binary system

We present a new method to achieve high-contrast images using segmented and/or on-axis telescopes. Our approach relies on using two sequential Deformable Mirrors to compensate for the large amplitude excursions in the telescope aperture due to secondary support structures and/or segment gaps. In this configuration the parameter landscape of Deformable Mirror Surfaces that yield high contrast Point Spread Functions is not linear, and non-linear methods are needed to find the true minimum in the optimization topology. We solve the highly non-linear Monge-Ampere equation that is the fundamental equation describing the physics of phase induced amplitude modulation. We determine the optimum configuration for our two sequential Deformable Mirror system and show that high-throughput and high contrast solutions can be achieved using realistic surface deformations that are accessible using existing technologies. We name this process Active Compensation of Aperture Discontinuities (ACAD). We show that for geometries similar to JWST, ACAD can attain at least 10^-7 in contrast and an order of magnitude higher for both the future Extremely Large Telescopes and on-axis architectures reminiscent of HST. We show that the converging non-linear mappings resulting from our Deformable Mirror shapes actually damp near-field diffraction artifacts in the vicinity of the discontinuities. Consequently, ACAD is a true broadband solution to the problem of high-contrast imaging with segmented and/or on-axis apertures. We finally show that once the non-linear solution is found, fine tuning with linear methods used in wavefront control can be applied to further contrast by another order of magnitude. Generally speaking, the ACAD technique can be used to significantly improve a broad class of telescope designs for a variety of problems.

Small-angle coronagraphy is technically and scientifically appealing because it enables the use of smaller telescopes, allows covering wider wavelength ranges, and potentially increases the yield and completeness of circumstellar environment - exoplanets and disks - detection and characterization campaigns. However, opening up this new parameter space is challenging. Here we will review the four posts of high contrast imaging and their intricate interactions at very small angles (within the first 4 resolution elements from the star). The four posts are: choice of coronagraph, optimized wavefront control, observing strategy, and post-processing methods. After detailing each of the four foundations, we will present the lessons learned from the 10+ years of operations of zeroth and first-generation adaptive optics systems. We will then tentatively show how informative the current integration of second-generation adaptive optics system is, and which lessons can already be drawn from this fresh experience. Then, we will review the current state of the art, by presenting world record contrasts obtained in the framework of technological demonstrations for space-based exoplanet imaging and characterization mission concepts. Finally, we will conclude by emphasizing the importance of the cross-breeding between techniques developed for both ground-based and space-based projects, which is relevant for future high contrast imaging instruments and facilities in space or on the ground.

We describe a new method to achieve point spread function (PSF) subtractions for high- contrast imaging using Principal Component Analysis (PCA) that is applicable to both point sources or extended objects (disks). Assuming a library of reference PSFs, a Karhunen-Lo`eve transform of theses references is used to create an orthogonal basis of eigenimages, on which the science target is projected. For detection this approach provides comparable suppression to the Locally Optimized Combination of Images (LOCI) algorithm, albeit with increased robustness to the algorithm parameters and speed enhancement. For characterization of detected sources the method enables forward modeling of astrophysical sources. This alleviates the biases in the astrometry and photometry of discovered faint sources, which are usually associated with LOCI- based PSF subtractions schemes. We illustrate the algorithm performance using archival Hubble Space Telescope (HST) images, but the approach may also be considered for ground-based data acquired with Angular Differential Imaging (ADI) or integral-field spectrographs (IFS).

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.

High-contrast imaging instruments are now being equipped with integral field spectrographs (IFS) to facilitate the detection and characterization of faint substellar companions. Algorithms currently envisioned to handle IFS data, such as the Locally Optimized Combination of Images (LOCI) algorithm, rely upon aggressive point-spread-function (PSF) subtraction, which is ideal for initially identifying companions but results in significantly biased photometry and spectroscopy due to unwanted mixing with residual starlight. This spectro-photometric issue is further complicated by the fact that algorithmic color response is a function of the companion's spectrum, making it difficult to calibrate the effects of the reduction without using iterations involving a series of injected synthetic companions. In this paper, we introduce a new PSF calibration method, which we call "damped LOCI", that seeks to alleviate these concerns. By modifying the cost function that determines the weighting coefficients used to construct PSF reference images, and also forcing those coefficients to be positive, it is possible to extract companion spectra with a precision that is set by calibration of the instrument response and transmission of the atmosphere, and not by post-processing. We demonstrate the utility of this approach using on-sky data obtained with the Project 1640 IFS at Palomar. Damped-LOCI does not require any iterations on the underlying spectral type of the companion, nor does it rely upon priors involving the chromatic and statistical properties of speckles. It is a general technique that can readily be applied to other current and planned instruments that employ IFS's.

The past decade has seen a significant growth in research targeted at space based observatories for imaging exo-solar planets. The challenge is in designing an imaging system for high-contrast. Even with a perfect coronagraph that modifies the point spread function to achieve high-contrast, wavefront sensing and control is needed to correct the errors in the optics and generate a "dark hole". The high-contrast imaging laboratory at Princeton University is equipped with two Boston Micromachines Kilo-DMs. We review here an algorithm designed to achieve high-contrast on both sides of the image plane while minimizing the stroke necessary from each deformable mirror (DM). This algorithm uses the first DM to correct for amplitude aberrations and the second DM to create a flat wavefront in the pupil plane. We then show the first results obtained at Princeton with this correction algorithm, and we demonstrate a symmetric dark hole in monochromatic light.

The purpose of this paper is to present the results of a theoretical study pertaining to the feasibility of PIAA units using Deformable Mirrors. We begin by reviewing the general derivation of the design equations driving PIAA. We then show how to solve these equations for square apertures and show the performance of pure PIAA systems in the ray optics regime. We tie these design equations into the study of edge diffraction effects, and provide a general expression for the field after a full propagation through a PIAA coronagraph. Third, we illustrate how a combination of pre and post apodisers yields to a contrast of 10^10 even in the presence of diffractive effects, for configuration with neither wavefront errors or wavefront control. Finally we present novel PIAA configurations over square apertures which circumvent the constraints on the manufacturing of PIAA optics by inducing the apodisation with two square Deformable Mirrors (DM). Such solutions rely on pupil size smaller than currently envisioned static PIAA solutions and thus require aggressive pre and post-apodizing screens in order to mitigate for diffractive effect between the two mirrors. As a result they are associated to significant loss in performance, throughput in particular.

The specification of polishing requirements for the optics in coronagraphs dedicated to exo-planet detection requires careful and accurate optical modelling. Numerical representations of the propagation of aberrations through the system as well as simulations of the broadband wavefront compensation system using multiple DMs are critical when one devises an error budget for such a class of instruments. In this communication we introduce an analytical tool that serves this purpose for Phase Induced Amplitude Apodisation (PIAA) coronagraphs. We first start by deriving the analytical form of the propagation of a harmonic ripple through a PIAA unit. Using this result we derive the chromaticity of the field at any plane in the optical train of a telescope equipped with such a coronagraph. Finally we study the chromatic response of a sequential DM wavefront actuator correcting such a corrugated field and thus quantify the requirements on the manufacturing of PIAA mirrors

HR 8799 is currently the only multiple-planet system that has been detected with direct imaging, with four giant planets orbiting at large separations from this young late A star. Orbital motion provides insight into the stability, and possible formation mechanisms of this planetary system. Dynamical studies can also provide constraints on the planets' masses, which help calibrate evolutionary models. Yet, measuring the orbital motion is a very difficult task because the long-period orbits (50-500 yr) require long time baselines and high-precision astrometry. This paper studies the three planets HR 8799b, c and d in the archival data set of HR 8799 obtained with the HST NICMOS coronagraph in 1998. The detection of all three planets is made possible by a careful optimization of the LOCI algorithm. This work confirms previous astrometry for planet b, and presents new detections and astrometry for c and d. These HST images provide a ten-year baseline with the discovery images from 2008, and therefore offer a unique opportunity to constrain their orbital motion now. Recent dynamical studies of this system show the existence of a few possible stable solutions involving mean motion resonances, where the interaction between c and d plays a major role. We study the compatibility of a few of these stable scenarios (1d:1c, 1d:2c, or 1d:2c:4d) with the new astrometric data from HST. In the hypothesis of a 1d:2c:4b mean motion resonance our best orbit fit is close to the stable solution previously identified for a three-planet system, and involves low eccentricity for planet d (ed = 0.10) and moderate inclination of the system (i = 28.0 deg), assuming a coplanar system, circular orbits for b and c, and exact resonance with integer period ratios. Under these assumptions, we can place strong constraints on the inclination of the system (27.3 - 31.4 deg) and on the eccentricity for d ed < 0.46.

The vortex coronagraph is one of the most promising coronagraphs for high contrast imaging because of its simplicity, small inner working angle, high throughput, and clear off-axis discovery space. However, as with most coronagraphs, centrally-obscured on-axis telescopes degrade contrast. Based on the remarkable ability of vortex coronagraphs to move light between the interior and exterior of pupils, we propose a method, based on multiple vortices, that, without sacrificing throughput, reduces the residual light leakage to (a/A)^n, with n >=4, and a and A being the radii of the central obscuration and primary mirror, respectively. This method thus enables high contrasts to be reached even with an on-axis telescope.

Project 1640 is a high-contrast imaging instrument recently commissioned at Palomar observatory. A combination of a coronagraph with an integral field spectrograph (IFS), Project 1640 is designed to detect and characterize extrasolar planets, brown dwarfs, and circumstellar material orbiting nearby stars. In this paper, we present our data processing techniques for improving upon instrument raw sensitivity via the removal of quasi-static speckles. Our approach utilizes the chromatic image diversity provided by the IFS in combination with the locally-optimized combination of images (LOCI) algorithm to suppress the intensity of residual contaminating light in close angular proximity to target stars. We describe the Project 1640 speckle suppression pipeline (PSSP) and demonstrate the ability to detect companions with brightness comparable to and below that of initial speckle intensities using on-sky commissioning data. Our preliminary results indicate that suppression factors of at least one order of magnitude are consistently possible, reaching $5\sigma$ contrast levels of $2.1\times10^{-5}$ at $1\arcsec$ in the H-band in 20 minutes of on-source integration time when non-common-path errors are reasonably well-calibrated. These results suggest that near-infrared contrast levels of order $\approx10^{-7}$ at subarcsecond separations will soon be possible for Project 1640 and similarly designed instruments that receive a diffraction-limited beam corrected by adaptive optics (AO) systems employing deformable mirrors with high actuator-density.

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.

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.