Search

The Near-Infrared Imager and Slitless Spectrograph (NIRISS) is the science module of the Canadian-built Fine Guidance Sensor (FGS) onboard the James Webb Space Telescope (JWST). NIRISS has four observing modes: 1) broadband imaging featuring seven of the eight NIRCam broadband filters, 2) wide-field slitless spectroscopy (WFSS) at a resolving power of $\sim$150 between 0.8 and 2.2 $\mu$m, 3) single-object cross-dispersed slitless spectroscopy (SOSS) enabling simultaneous wavelength coverage between 0.6 and 2.8 $\mu$m at R$\sim$700, a mode optimized for exoplanet spectroscopy of relatively bright ($J<6.3$) stars and 4) aperture masking interferometry (AMI) between 2.8 and 4.8 $\mu$m enabling high-contrast ($\sim10^{-3}-10^{-4}$) imaging at angular separations between 70 and 400 milliarcsec for relatively bright ($M<8$) sources. This paper presents an overview of the NIRISS instrument, its design, its scientific capabilities, and a summary of in-flight performance. NIRISS shows significantly better response shortward of $\sim2.5\,\mu$m resulting in 10-40% sensitivity improvement for broadband and low-resolution spectroscopy compared to pre-flight predictions. Two time-series observations performed during instrument commissioning in the SOSS mode yield very stable spectro-photometry performance within $\sim$10% of the expected noise. The first space-based companion detection of the tight binary star AB Dor AC through AMI was demonstrated.

Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit., Comment: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the orbital fit and, in combination with absolute astrometry from GAIA, place a 3$\sigma$ upper limit of 18.5 $M_J$ on the companion's mass. To interpret the ALMA observations, we used radiative transfer modeling to determine the disk properties. We find that the disk is well resolved and nearly edge on. While the misalignment angle between the disk and the orbit is weakly constrained due to the short orbital arc available, the data strongly support a (near) coplanar geometry for the system. Furthermore, we find that the models that describe the ALMA data best have an inner radius that is close to the companion's semi-major axis. Including a posteriori knowledge of the system's SED further narrows the constraints on the disk's inner radius and place it at a location that is in reasonable agreement with, possibly interior to, predictions from existing dynamical models of disk truncation by an interior substellar companion. HR\,2562 has the potential over the next few years to become a new testbed for dynamical interaction between a debris disk and a substellar companion., Comment: 22 pages, 17 figures

The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$\mu$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$\mu$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles., Comment: 21 pages, 7 figures. Accepted for Publication in PASP

The direct detection of exoplanets with high-contrast instruments can be boosted with high spectral resolution. For integral field spectrographs yielding hyperspectral data, this means that the field of view consists of diffracted starlight spectra and a spatially localized planet. Analysis usually relies on cross-correlation with theoretical spectra. In a purely blind-search context, this supervised strategy can be biased with model mismatch and/or be computationally inefficient. Using an approach that is inspired by the remote-sensing community, we aim to propose an alternative to cross-correlation that is fully data-driven, which decomposes the data into a set of individual spectra and their corresponding spatial distributions. This strategy is called spectral unmixing. We used an orthogonal subspace projection to identify the most distinct spectra in the field of view. Their spatial distribution maps were then obtained by inverting the data. These spectra were then used to break the original hyperspectral images into their corresponding spatial distribution maps via non-negative least squares. The performance of our method was evaluated and compared with a cross-correlation using simulated hyperspectral data with medium resolution from the ELT/HARMONI integral field spectrograph. We show that spectral unmixing effectively leads to a planet detection solely based on spectral dissimilarities at significantly reduced computational cost. The extracted spectrum holds significant signatures of the planet while being not perfectly separated from residual starlight. The sensitivity of the supervised cross-correlation is three to four times higher than with unsupervised spectral unmixing, the gap is biased toward the former because the injected and correlated spectrum match perfectly. The algorithm was furthermore vetted on real data obtained with VLT/SINFONI of the beta Pictoris system., Comment: 13 pages, 8 figures

Medium-resolution integral-field spectrographs (IFS) coupled with adaptive-optics such as Keck/OSIRIS, VLT/MUSE, or SINFONI are appearing as a new avenue for enhancing the detection and characterization capabilities of young, gas giant exoplanets at large heliocentric distances (>5 au). We analyzed K-band VLT/SINFONI medium-resolution (R_lambda~5577) observations of the young giant exoplanet HIP 65426 b. Our dedicated IFS data analysis toolkit (TExTRIS) optimized the cube building, star registration, and allowed for the extraction of the planet spectrum. A Bayesian inference with the nested sampling algorithm coupled with the self-consistent forward atmospheric models BT-SETTL15 and Exo-REM using the ForMoSA tool yields Teff=1560 +/- 100K, log(g)<4.40dex, [M/H]=0.05 +/- 0.22dex, and an upper limit on the C/O ratio (<0.55). The object is also re-detected with the so-called "molecular mapping" technique. The technique yields consistent atmospheric parameters, but the loss of the planet pseudo-continuum in the process degrades or modifies the constraints on these parameters. The solar to sub-solar C/O ratio suggests an enrichment by solids at formation if the planet was formed beyond the water snowline (>20 au) by core-accretion. However, a formation by gravitational instability can not be ruled out. The metallicity is compatible with the bulk enrichment of massive Jovian planets from the Bern planet population models. Finally, we measure a radial velocity of 26 +/- 15km/s compatible with our revised measurement on the star. This is the fourth imaged exoplanet for which a radial velocity can be evaluated, illustrating the potential of such observations for assessing the coevolution of imaged systems belonging to star forming regions, such as HIP 65426., Comment: 21 pages, 16 figures, 3 Appendix

We present a velocimetric and spectropolarimetric analysis of 27 observations of the 22-Myr M1 star AU Microscopii (Au Mic) collected with the high-resolution $YJHK$ (0.98-2.35 $\mu$m) spectropolarimeter SPIRou from 2019 September 18 to November 14. Our radial velocity (RV) time-series exhibits activity-induced fluctuations of 45 m/s RMS, about three times smaller than those measured in the optical domain, that we filter using Gaussian Process Regression. We report a 3.9$\sigma$-detection of the recently-discovered 8.46-d transiting planet AU Mic b, with an estimated mass of $17.1^{+4.7}_{-4.5}$ M$_{\odot}$ and a bulk density of $1.3 \pm 0.4$ g/cm$^{-3}$, inducing a RV signature of semi-amplitude $K=8.5^{+2.3}_{-2.2}$ m/s in the spectrum of its host star. A consistent detection is independently obtained when we simultaneously image stellar surface inhomogeneities and estimate the planet parameters with Zeeman-Doppler Imaging (ZDI). Using ZDI, we invert the time series of unpolarized and circularly-polarized spectra into surface brightness and large-scale magnetic maps. We find a mainly poloidal and axisymmetric field of 475 G, featuring, in particular, a dipole of 450 G tilted at 19{\deg} to the rotation axis. Moreover, we detect a strong differential rotation of d$\Omega = 0.167 \pm 0.009$ rad/d shearing the large-scale field, about twice stronger than that shearing the brightness distribution, suggesting that both observables probe different layers of the convective zone. Even though we caution that more RV measurements are needed to accurately pin down the planet mass, AU Mic b already appears as a prime target for constraining planet formation models, studying the interactions with the surrounding debris disk, and characterizing its atmosphere with upcoming space- and ground-based missions., Comment: 19 pages, 20 figures, MNRAS, in press

Nearly 15 years of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star $\beta$ Pictoris. The $\delta$ Scuti pulsations of the star, overwhelming planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than in the literature to model the activity of the star. We show that a Gaussian Process (GP) with a stochastically driven damped harmonic oscillator kernel can model the $\delta$ Scuti pulsations. It provides similar results as parametric models but with a simpler framework, using only 3 hyperparameters. It also enables to model poorly sampled RV data, that were excluded from previous analysis, hence extending the RV baseline by nearly five years. Altogether, the orbit and the mass of both planets can be constrained from RV only, which was not possible with the parametric modelling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for $\beta$ Pic b favour a low eccentricity of $0.029^{+0.061}_{-0.024}$ and a relatively short period of $21.1^{+2.0}_{-0.8}$ yr. The orbit of $\beta$ Pic c is eccentric with $0.206^{+0.074}_{-0.063}$ with a period of $3.36\pm0.03$ yr. We find model-independent masses of $11.7\pm1.4$ and $8.5\pm0.5$ M$_{Jup}$ for $\beta$ Pic b and c, respectively, assuming coplanarity. The mass of $\beta$ Pic b is consistent with the hottest start evolutionary models, at an age of $25\pm3$ Myr. A direct direction of $\beta$ Pic c would provide a second calibration measurement in a coeval system., Comment: 25 pages, 18 figures, Accepted to AJ, Typos Corrected

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., Comment: Accepted for publication in AJ. 3 tables, 11 figures

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., Comment: arXiv resubmission with typographical corrections. Accepted for publication in AJ. 19 figures, 7 tables

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., Comment: Accepted for publication in the Astronomical Journal

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., Comment: 25 pages, 12 figures

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., Comment: 29 Pages, 23 Figures. Accepted to AJ

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., Comment: 20 pages, 9 figures

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., Comment: 14 pages, 9 figures. Accepted for publication in the Astronomical Journal

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}$)., Comment: 17 pages, 11 figures. Accepted for publication in the Astronomical Journal

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., Comment: 36 pages, 17 figures, accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems. Updated version includes revisions made during the referee process

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., Comment: 28 Pages, 13 Figures, Accepted to JATIS

We present the first statistical analysis of exoplanet direct imaging surveys combining adaptive optics imaging at small separations with deep seeing-limited observations at large separations allowing us to study the entire orbital separation domain from 5 to 5000~au simultaneously. Our sample of 344 stars includes only confirmed members of nearby young associations and is based on all AO direct-imaging detection limits readily available online, with addition of our own previous seeing limited surveys. Assuming that the companion distribution in mass and semi-major axis follows a power law distribution and adding a dependence on the mass of the host star, such as $d^2n\propto fM^{\alpha}a^{\beta} (M_\star/M_{\odot})^{\gamma}$d$ M $d$a$, we constrain the parameters to obtained $\alpha=-0.18^{+0.77}_{-0.65}$, $\beta=-1.43^{+0.23}_{-0.24}$, and $\gamma=0.62^{+0.56}_{-0.50}$,at a 68\% confidence level, and we obtain $f=0.11^{+0.11}_{-0.05}$, for the overall planet occurrence rate for companions with masses between 1 to 20~\mj\ in the range 5--5000~au. Thus, we find that occurrence of companions is negatively correlated with semi-major axis and companion mass (marginally) but is positively correlated with the stellar host mass. Our inferred mass distribution is in good agreement with other distributions found previously from direct imaging surveys for planets and brown dwarfs, but is shallower as a function of mass than the distributions inferred by radial velocity surveys of gas giants in the 1--3\,au range. This may suggest that planets at these wide and very-wide separations represent the low-mass tail of the brown dwarfs and stellar companion distribution rather than an extension of the distribution of the inner planets., Comment: 26 pages, 8 figures Accepted to Astronomical Journal Updated Table 4

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., Comment: 29 pages, 18 figures, ApJ accepted

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., Comment: Accepted to AJ, 22 pages, 14 figures

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., Comment: 52 pages, 18 figures. AJ in press

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., Comment: 21 pages, 10 figures, accepted to AJ

We report results from the Wide-orbit Exoplanet search with InfraRed Direct imaging (WEIRD), a survey designed to search for Jupiter-like companions on very wide orbits (1000 to 5000 AU) around young stars ($<$120 Myr) that are known members of moving groups in the solar neighborhood ($<$70 pc). Sharing the same age, distance, and metallicity as their host while being on large enough orbits to be studied as "isolated" objects make such companions prime targets for spectroscopic observations and valuable benchmark objects for exoplanet atmosphere models. The search strategy is based on deep imaging in multiple bands across the near-infrared domain. For all 177 objects of our sample, $z_{ab}^\prime$, $J$, [3.6] and [4.5] images were obtained with CFHT/MegaCam, GEMINI/GMOS, CFHT/WIRCam, GEMINI/Flamingos-2, and $Spitzer$/IRAC. Using this set of 4 images per target, we searched for sources with red $z_{ab}^\prime$ and $[3.6]-[4.5]$ colors, typically reaching good completeness down to 2Mjup companions, while going down to 1Mjup for some targets, at separations of $1000-5000$ AU. The search yielded 4 candidate companions with the expected colors, but they were all rejected through follow-up proper motion observations. Our results constrain the occurrence of 1-13 Mjup planetary-mass companions on orbits with a semi-major axis between 1000 and 5000 AU at less than 0.03, with a 95\% confidence level., Comment: 55 pages, 16 figures, accepted to AJ

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., Comment: 14 pages, 6 figures, 1 table

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., Comment: Accepted to AJ. 22 pages, 11 figures, 3 tables

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., Comment: Accepted to AJ, 25 pages, 16 Figures

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., Comment: 21 pages, 3 figures, accepted in JATIS

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., Comment: 35 pages, 20 figures, Published in AJ

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., Comment: ApJ accepted

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}$., Comment: 27 pages, 19 figures, Accepted for publication in The Astronomical Journal

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., Comment: Accepted to AJ, 10 pages, 6 figures, 1 table

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., Comment: accepted to AJ

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., Comment: 32 pages, 28 figures, Accepted to AJ

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., Comment: Accepted for publication in The Astronomical Journal

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., Comment: 21 pages, 11 figures, Accepted to AJ

Atmospheric magnetic fields in stars with convective envelopes heat stellar chromospheres. Starting with the historical Mount Wilson monitoring program, CaH&K lines have been widely used to trace stellar magnetic activity, and as a proxy for rotation period and consequently for stellar age. Monitoring stellar activity has also become essential in filtering out false-positives due to magnetic activity in extra-solar planet surveys. The Ca H&Kemission is traditionally quantified through the R'HK-index, which compares the chromospheric flux in the doublet to the overall bolometric flux of the star. Much work has been done to characterize this index for FGK-dwarfs, but M-dwarfs were left out of these analyses and no calibration of their Ca ii H&K emission to an R'HK exists to date. We set out to characterize the magnetic activity of the low and very low-mass stars by providing a calibration of the R'HK-index that extends to the realm of M-dwarfs, and by evaluating the relation between R'HK and the rotation period. We calibrated the bolometric and photospheric factors for M-dwarfs to properly transform the S-index into the R'HK. We monitored magnetic activity through the Ca ii H&K emission lines in the HARPS M-dwarf sample. The R'HK index, like the fractional X-ray luminosity L_X/Lbol, shows a saturated correlation with rotation, with saturation setting in around a ten days rotation period. Above that period, slower rotators show weaker Ca ii activity, as expected. Under that period, the R'HK index saturates to ~10^-4. Stellar mass modulates the Ca ii activity, with R'HK showing a constant basal activity above 0.6Msun and then decreasing with mass between 0.6Msun and the fully-convective limit of 0.35Msun. Short-term variability of the activity correlates with its mean level, stars with higher R'HK index show larger R'HK variability, as previously observed for earlier spectral types., Comment: Accepted for publication in A&A. 21 pages, 10 figures, 6 tables

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., Comment: 12 pages, 7 figures. Accepted to AJ

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., Comment: Updated with proper references

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., Comment: 10 pages, 4 figures. Accepted for publication in ApJ Letters

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., Comment: 18 pages, 7 figures. Accepted to AJ. v2 corrects an error in plotting the full size of the Hill sphere in Figure 6

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., Comment: Accepted for publication in ApJL

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., Comment: 20 pages, 12 figures, accepted for publication in ApJ