Search

We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 K to 1800 K for effective temperature and from 3.0 to 5.0 for log($g$), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of sub-micron dust particles, similar to dereddening approaches applied to the peculiar red field L-dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L-dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from four years of astrometric monitoring is consistent with a $\sim$30-year period, orbital inclination of 147$^{\circ}$, and semimajor axis of 10 au, well within the estimated disk inner radius of $\sim$50 au. As one of very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion-disk interaction and provide a benchmark for substellar and planetary atmospheric study., Accepted for publication in AJ (October 15, 2020). 35 pages, 20 figures, with tables and appendices presented in their entirety

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

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

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

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

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers., Comment: 48 pages, 23 figures. This paper is part of a coordinated effort between groups. Accepted for publication in ApJ

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., 16 pages, 8 figures, accepted at 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

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

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies., Comment: 29 pages, 12 figures, accepted for publication in ApJ

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

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.

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

We report 885$\mu$m ALMA continuum flux densities for 24 Taurus members spanning the stellar/substellar boundary, with spectral types from M4 to M7.75. Of the 24 systems, 22 are detected at levels ranging from 1.0-55.6 mJy. The two non-detections are transition disks, though other transition disks in the sample are detected. Converting ALMA continuum measurements to masses using standard scaling laws and radiative transfer modeling yields dust mass estimates ranging from $\sim$0.3-20M$_{\oplus}$. The dust mass shows a declining trend with central object mass when combined with results from submillimeter surveys of more massive Taurus members. The substellar disks appear as part of a continuous sequence and not a distinct population. Compared to older Upper Sco members with similar masses across the substellar limit, the Taurus disks are brighter and more massive. Both Taurus and Upper Sco populations are consistent with an approximately linear relationship in $M_{dust}$ to $M_{star}$, although derived power-law slopes depend strongly upon choices of stellar evolutionary model and dust temperature relation. The median disk around early M-stars in Taurus contains a comparable amount of mass in small solids as the average amount of heavy elements in Kepler planetary systems on short-period orbits around M-dwarf stars, with an order of magnitude spread in disk dust mass about the median value. Assuming a gas:dust ratio of 100:1, only a small number of low-mass stars and brown dwarfs have a total disk mass amenable to giant planet formation, consistent with the low frequency of giant planets orbiting M-dwarfs., Comment: 41 pages and 32 figures, with all tables and appendices presented here in their entirety. Accepted for publication in AJ (November 26, 2017)

With the discovery of Y dwarfs by the WISE mission, the population of field brown dwarfs now extends to objects with temperatures comparable to those of Solar System planets. To investigate the atmospheres of these newly identified brown dwarfs, we have conducted a pilot study monitoring an initial sample of three late T-dwarfs (T6.5, T8 and T8.5) and one Y-dwarf (Y0) for infrared photometric variability at multiple epochs. With J-band imaging, each target was observed for a period of 1.0h to 4.5h per epoch, which covers a significant fraction of the expected rotational period. These measurements represent the first photometric monitoring for these targets. For three of the four targets (2M1047, Ross 458C and WISE0458), multi-epoch monitoring was performed, with the time span between epochs ranging from a few hours to ~2 years. During the first epoch, the T8.5 target WISE0458 exhibited variations with a remarkable min-to-max amplitude of 13%, while the second epoch light curve taken ~2 years later did not note any variability to a 3% upper limit. With an effective temperature of ~600 K, WISE0458 is the coldest variable brown dwarf published to-date, and combined with its high and variable amplitude makes it a fascinating target for detailed follow-up. The three remaining targets showed no significant variations, with a photometric precision between 0.8% and 20.0%, depending on the target brightness. Combining the new results with previous multi-epoch observations of brown dwarfs with spectral types of T5 or later, the currently identified variables have locations on the colour-colour diagram better matched by theoretical models incorporating cloud opacities rather than cloud-free atmospheres. This preliminary result requires further study to determine if there is a definitive link between variability among late-T dwarfs and their location on the colour-colour diagram., Comment: 9 pages, 6 figures, 3 tables, accepted for publication in MNRAS

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

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

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

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

We present the results of ALMA band 7 observations of dust and CO gas in the disks around 7 objects with spectral types ranging between M5.5 and M7.5 in Upper Scorpius OB1, and one M3 star in Ophiuchus. We detect unresolved continuum emission in all but one source, and the $^{12}$CO J=3-2 line in two sources. We constrain the dust and gas content of these systems using a grid of models calculated with the radiative transfer code MCFOST, and find disk dust masses between 0.1 and 1 M$_\oplus$, suggesting that the stellar mass / disk mass correlation can be extrapolated for brown dwarfs with masses as low as 0.05 M$_\odot$. The one disk in Upper Sco in which we detect CO emission, 2MASS J15555600, is also the disk with warmest inner disk as traced by its H - [4.5] photometric color. Using our radiative transfer grid, we extend the correlation between stellar luminosity and mass-averaged disk dust temperature originally derived for stellar mass objects to the brown dwarf regime to $\langle T_{dust} \rangle \approx 22 (L_{*} /L_{\odot})^{0.16} K$, applicable to spectral types of M5 and later. This is slightly shallower than the relation for earlier spectral type objects and yields warmer low-mass disks. The two prescriptions cross at 0.27 L$_\odot$, corresponding to masses between 0.1 and 0.2 M$_\odot$ depending on age., 9 pages,6 figures, accepted to ApJ on 26/01/2016

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$��$ 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., Accepted for publication in the Astrophysical Journal; 15 pages, 7 figures

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., Comment: 29 pages, 3 figures, 2 tables, and Supplementary Materials. published in Science Express on Aug 13 2015. List of authors and the magnitudes of the star were correted

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., 14 pages, 4 figures, accepted by ApJ Letters

With Herschel/PACS 134 low mass members of the Taurus star-forming region spanning the M4-L0 spectral type range and covering the transition from low mass stars to brown dwarfs were observed. Combining the new Herschel results with other programs, a total of 150 of the 154 M4-L0 Taurus members members have observations with Herschel. Among the 150 targets, 70um flux densities were measured for 7 of the 7 ClassI objects, 48 of the 67 ClassII members, and 3 of the 76 ClassIII targets. For the detected ClassII objects, the median 70um flux density level declines with spectral type, however, the distribution of excess relative to central object flux density does not change across the stellar/substellar boundary in the M4-L0 range. Connecting the 70um TBOSS values with the results from K0-M3 ClassII members results in the first comprehensive census of far-IR emission across the full mass spectrum of the stellar and substellar population of a star-forming region, and the median flux density declines with spectral type in a trend analogous to the flux density decline expected for the central objects. SEDs were constructed for all TBOSS targets covering the optical to far-IR range and extending to the submm/mm for a subset of sources. Based on an initial exploration of the impact of different physical parameters; inclination, scale height and flaring have the largest influence on the PACS flux densities. From the 24um to 70um spectral index of the SEDs, 5 new candidate transition disks were identified. The steep 24um to 70um slope for a subset of 8 TBOSS targets may be an indication of truncated disks in these systems.Two examples of mixed pair systems that include secondaries with disks were measured. Finally, comparing the TBOSS results with a Herschel study of Ophiuchus brown dwarfs reveals a lower fraction of disks around the Taurus substellar population., 64 pages, 33 figures, 12 tables, accepted for publication in A&A

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 ∼7% and scales linearly with central star brightness and that the shape of the satellite spot spectrum varies on the 2% level., Ground-Based and Airborne Instrumentation for Astronomy V, June 22-26, 2014, Series: Proceedings of SPIE; no. 9147

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., Comment: 9 pages, 8 Figures, 1 table, Accepted to ApJ

We present the preliminary findings of an investigation of the multiplicity of debris disk stars identified within our Volume-limited A-star (VAST) multiplicity survey. Previous studies have produced conflicting results regarding the multiplicity fraction of debris disk-hosting stars compared with non-excess stars. By combining our large-scale volume-limited AO survey of A-type stars with the all-sky WISE catalogue, we have investigated the frequency of binary companions to a large sample of A-type stars with and without measured 22μm excess. The results of this study will allow for a greater understanding of the interaction between a companion star and a circumstellar debris disk, informing future study into the formation and stability of planetary-mass companions within binary systems., 299th Symposium of the International Astronomical Union, June 2-7, 2013, Victoria, Canada, Series: Proceedings of the International Astronomical Union

With the combination of Herschel PACS (far-IR), submillimeter (submm) and millimeter (mm) ground-based observations, we are leading an investigation of all members with spectral types M4 and onwards in the Taurus star-forming region. This complete census spans the stellar to substellar boundary (M6.25), and the wavelength range covers the transition from optically thick to optically thin emission. From our Hershel PACS observations of ~135 sources, we obtain a detection rate of ~40% at 70 μm and a detection rate of ~20% at 160 μm, and provide the first far-IR measurements for the majority of these members. With our complementary submm and mm observations, best fit SED models from the radiative transfer code MCFOST will be used to infer disk properties such as scale height, mass, outer radius and maximum dust grain size. These comprehensive population statistics of disks are critical for testing star/brown dwarf and planet formation models around these later type members of Taurus.