Source author record

Kimberly Ward-Duong

Kimberly Ward-Duong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

ALMA Discovery of a Disk around the Planetary-mass Companion SR 12 c

We report an Atacama Large Millimeter/submillimeter Array 0.88 mm (Band 7) continuum detection of the accretion disk around SR 12 c, an $\sim$11 $M_{\rm Jup}$ planetary-mass companion (PMC) orbiting its host binary at 980 au. This is the first submillimeter detection of a circumplanetary disk around a wide PMC. The disk has a flux density of $127 \pm14~μ$Jy and is not resolved by the $\sim$0.1" beam, so the dust disk radius is likely less than 5 au and can be much smaller if the dust continuum is optically thick. If, however, the dust emission is optically thin, then the SR 12 c disk has a comparable dust mass to the circumplanetary disk around PDS 70 c but is about five times lower than that of the $\sim$12 $M_{\rm Jup}$ free-floating OTS 44. This suggests that disks around bound and unbound planetary-mass objects can span a wide range of masses. The gas mass estimated with an accretion rate of $10^{-11}~M_\odot$ yr$^{-1}$ implies a gas-to-dust ratio higher than 100. If cloud absorption is not significant, a nondetection of ${}^{12}$CO(3-2) implies a compact gas disk around SR 12 c. Future sensitive observations may detect more PMC disks at 0.88 mm flux densities of $\lesssim$100 $μ$Jy.

preprint2022arXiv

HST/WFC3 H$α$ Direct-Imaging Detection of a Point-like Source in the Disk Cavity of AB Aur

Accreting protoplanets enable the direct characterization of planet formation. As part of a high-contrast imaging search for accreting planets with the Hubble Space Telescope (HST) Wide Field Camera 3, we present H$α$ images of AB Aurigae (AB Aur), a Herbig Ae/Be star harboring a transition disk. The data were collected in two epochs of direct-imaging observations using the F656N narrow-band filter. After subtracting the point spread function of the primary star, we identify a point-like source located at a P.A. of $182.5^{\circ}\pm1.4^{\circ}$ and a separation of $600\pm22$~mas relative to the host star. The position is consistent with the recently identified protoplanet candidate AB Aur b. The source is visible in two individual epochs separated by ${\sim}50$ days and the H$α$ intensities in the two epochs agree. The H$α$ flux density is $F_ν=1.5\pm0.4$~mJy, $3.2\pm0.9$ times of the optical continuum determined by published HST/STIS photometry. In comparison to PDS 70 b and c, the H$α$ excess emission is weak. The central star is accreting and the stellar H$α$ emission has a similar line-to-continuum ratio as seen in AB Aur b. We conclude that both planetary accretion and scattered stellar light are possible sources of the H$α$ emission, and the H$α$ detection alone does not validate AB Aur b as an accreting protoplanet. Disentangling the origin of the emission will be crucial for probing planet formation in the AB Aur disk.

preprint2022arXiv

Images of Embedded Jovian Planet Formation At A Wide Separation Around AB Aurigae

Direct images of protoplanets embedded in disks around infant stars provide the key to understanding the formation of gas giant planets like Jupiter. Using the Subaru Telescope and Hubble Space Telescope, we find evidence for a jovian protoplanet around AB Aurigae orbiting at a wide projected separation (93 au), likely responsible for multiple planet-induced features in the disk. Its emission is reproducible as reprocessed radiation from an embedded protoplanet. We also identify two structures located at 430-580 au that are candidate sites of planet formation. These data reveal planet formation in the embedded phase and a protoplanet discovery at wide, > 50 au separations characteristic of most imaged exoplanets. With at least one clump-like protoplanet and multiple spiral arms, the AB Aur system may also provide the evidence for a long-considered alternative to the canonical model for Jupiter's formation: disk (gravitational) instability.

preprint2022arXiv

Performance of near-infrared high-contrast imaging methods with JWST from commissioning

The James Webb Space Telescope (JWST) will revolutionize the field of high-contrast imaging and enable both the direct detection of Saturn-mass planets and the characterization of substellar companions in the mid-infrared. While JWST will feature unprecedented sensitivity, angular resolution will be the key factor when competing with ground-based telescopes. Here, we aim to characterize the performance of several extreme angular resolution imaging techniques available with JWST in the 3-5 micron regime based on data taken during commissioning. Firstly, we introduce custom tools to simulate, reduce, and analyze NIRCam and MIRI coronagraphy data and use these tools to extract companion detection limits from on-sky NIRCam round and bar mask coronagraphy observations. Secondly, we present on-sky NIRISS aperture masking interferometry (AMI) and kernel phase imaging (KPI) observations from which we extract companion detection limits using the publicly available fouriever tool. Scaled to a total integration time of one hour and a target of the brightness of AB Dor, we find that NIRISS AMI and KPI reach contrasts of $\sim$7-8 mag at $\sim$70 mas and $\sim$9 mag at $\sim$200 mas. Beyond $\sim$250 mas, NIRCam coronagraphy reaches deeper contrasts of $\sim$13 mag at $\sim$500 mas and $\sim$15 mag at $\sim$2 arcsec. While the bar mask performs $\sim$1 mag better than the round mask at small angular separations $\lesssim$0.75 arcsec, it is the other way around at large angular separations $\gtrsim$1.5 arcsec. Moreover, the round mask gives access to the full 360 deg field-of-view which is beneficial for the search of new companions. We conclude that already during the instrument commissioning, JWST high-contrast imaging in the L- and M-bands performs close to its predicted limits.

preprint2020arXiv

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

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.

preprint2020arXiv

HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window

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.

preprint2020arXiv

Multiband GPI Imaging of the HR 4796A Debris Disk

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

preprint2020arXiv

Multiband Polarimetric Imaging of HR 4796A with the Gemini Planet Imager

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.

preprint2020arXiv

Revised Astrometric Calibration of the Gemini Planet Imager

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.

preprint2019arXiv

An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

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$σ$ 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}$).

preprint2019arXiv

Detection of a low-mass stellar companion to the accelerating A2IV star HR 1645

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.

preprint2019arXiv

First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

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.

preprint2019arXiv

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry

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.

preprint2016arXiv

Adaptive Optics imaging of VHS 1256-1257: A Low Mass Companion to a Brown Dwarf Binary System

Recently, Gauza et al. (2015) reported the discovery of a companion to the late M-dwarf, VHS J125601.92-125723.9 (VHS 1256-1257). The companion's absolute photometry suggests its mass and atmosphere are similar to the HR 8799 planets. However, as a wide companion to a late-type star, it is more accessible to spectroscopic characterization. We discovered that the primary of this system is an equal-magnitude binary. For an age $\sim300$ Myr the A and B components each have a mass of $64.6^{+0.8}_{-2.0}~M_{\mathrm{Jup}}$, and the b component has a mass of $11.2^{+9.7}_{-1.8}$, making VHS 1256-1257 only the third brown dwarf triple system. There exists some tension between the spectrophotometric distance of $17.2\pm2.6$ pc and the parallax distance of $12.7\pm1.0$ pc. At 12.7 pc VHS1256-1257 A and B would be the faintest known M7.5 objects, and are even faint outliers among M8 types. If the larger spectrophotmetric distance is more accurate than the parallax, then the mass of each component increases. In particular, the mass of the b component increases well above the deuterium burning limit to $\sim35~M_{\mathrm{Jup}}$ and the mass of each binary component increases to $73^{+20}_{-17}~M_{\mathrm{Jup}}$. At 17.1 pc, the UVW kinematics of the system are consistent with membership in the AB~Dor moving group. The architecture of the system resembles a hierarchical stellar multiple suggesting it formed via an extension of the star-formation process to low masses. Continued astrometric monitoring will resolve this distance uncertainty and will provide dynamical masses for a new benchmark system.

preprint2016arXiv

Constraints on the architecture of the HD 95086 planetary system with the Gemini Planet Imager

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.

preprint2016arXiv

Discovery and Validation of a High-Density sub-Neptune from the K2 Mission

We report the discovery of BD+20594b, a high density sub-Neptune exoplanet, made using photometry from Campaign 4 of the two-wheeled Kepler (K2) mission, ground-based radial velocity follow-up from HARPS and high resolution lucky and adaptive optics imaging obtained using AstraLux and MagAO, respectively. The host star is a bright ($V=11.04$, $K_s = 9.37$), slightly metal poor ([Fe/H]$=-0.15\pm 0.05$ dex) solar analogue located at $152.1^{+9.7}_{-7.4}$ pc from Earth, for which we find a radius of $R_*=0.928^{+0.055}_{-0.040}R_\odot$ and a mass of $M_* = 0.961^{+0.032}_{-0.029}M_\odot$. A joint analysis of the K2 photometry and HARPS radial velocities reveal that the planet is in a $\approx 42$ day orbit around its host star, has a radius of $2.23^{+0.14}_{-0.11}R_\oplus$, and a mass of $16.3^{+6.0}_{-6.1}M_\oplus$. Although the data at hand puts the planet in the region of the mass-radius diagram where we could expect planets with a pure rock (i.e. magnesium silicate) composition using two-layer models (i.e., between rock/iron and rock/ice compositions), we discuss more realistic three-layer composition models which can explain the high density of the discovered exoplanet. The fact that the planet lies in the boundary between "possibly rocky" and "non-rocky" exoplanets, makes it an interesting planet for future RV follow-up.

preprint2016arXiv

Gemini Planet Imager Observational Calibration XIII: Wavelength Calibration Improvements, Stability, and Nonlinearity

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

preprint2016arXiv

Imaging an 80 AU Radius Dust Ring Around the F5V Star HD 157587

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.

preprint2016arXiv

Spectroscopic characterization of HD 95086 b with the Gemini Planet Imager

We present new $H$ (1.5-1.8 $μ$m) photometric and $K_1$ (1.9-2.2 $μ$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$ ($λ/δλ\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.

preprint2016arXiv

The Orbit and Transit Prospects for $β$ Pictoris b constrained with One Milliarcsecond Astrometry

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ève 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 $β$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$σ$-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 $β$ Pic b firmly rule out the possibility of a transit of the planet at 10-$σ$ significance. However, we confirm that the Hill sphere of $β$ 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.

preprint2016arXiv

The Peculiar Debris Disk of HD 111520 as Resolved by the Gemini Planet Imager

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.

preprint2015arXiv

$β$ Pictoris' inner disk in polarized light and new orbital parameters for $β$ Pictoris b

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

preprint2015arXiv

Astrometric Confirmation and Preliminary Orbital Parameters of the Young Exoplanet 51 Eridani b with the Gemini Planet Imager

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

preprint2015arXiv

Direct imaging of an asymmetric debris disk in the HD 106906 planetary system

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.

preprint2015arXiv

First Scattered-Light Image of the Debris Disk around HD 131835 with the Gemini Planet Imager

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

preprint2015arXiv

Gemini Planet Imager Observations of the AU Microscopii Debris Disk: Asymmetries within One Arcsecond

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.

preprint2015arXiv

The LEECH Exoplanet Imaging Survey: Characterization of the Coldest Directly Imaged Exoplanet, GJ 504 b, and Evidence for Super-Stellar Metallicity

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

preprint2014arXiv

Does the Debris Disk around HD 32297 Contain Cometary Grains?

We present an adaptive optics imaging detection of the HD 32297 debris disk at L' (3.8 \microns) obtained with the LBTI/LMIRcam infrared instrument at the LBT. The disk is detected at signal-to-noise per resolution element ~ 3-7.5 from ~ 0.3-1.1" (30-120 AU). The disk at L' is bowed, as was seen at shorter wavelengths. This likely indicates the disk is not perfectly edge-on and contains highly forward scattering grains. Interior to ~ 50 AU, the surface brightness at L' rises sharply on both sides of the disk, which was also previously seen at Ks band. This evidence together points to the disk containing a second inner component located at $\lesssim$ 50 AU. Comparing the color of the outer (50 $< r$/AU $< 120$) portion of the disk at L' with archival HST/NICMOS images of the disk at 1-2 \microns allows us to test the recently proposed cometary grains model of Donaldson et al. 2013. We find that the model fails to match the disk's surface brightness and spectrum simultaneously (reduced chi-square = 17.9). When we modify the density distribution of the model disk, we obtain a better overall fit (reduced chi-square = 2.9). The best fit to all of the data is a pure water ice model (reduced chi-square = 1.06), but additional resolved imaging at 3.1 \microns is necessary to constrain how much (if any) water ice exists in the disk, which can then help refine the originally proposed cometary grains model.

preprint2014arXiv

Gemini Planet Imager Observational Calibrations I: Overview of the GPI Data Reduction Pipeline

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

preprint2014arXiv

Gemini Planet Imager Observational Calibrations VIII: Characterization and Role of Satellite Spots

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

preprint2014arXiv

High contrast imaging at the LBT: the LEECH exoplanet imaging survey

In Spring 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its $\sim$130-night campaign from the Large Binocular Telescope (LBT) atop Mt Graham, Arizona. This survey benefits from the many technological achievements of the LBT, including two 8.4-meter mirrors on a single fixed mount, dual adaptive secondary mirrors for high Strehl performance, and a cold beam combiner to dramatically reduce the telescope's overall background emissivity. LEECH neatly complements other high-contrast planet imaging efforts by observing stars at L' (3.8 $μ$m), as opposed to the shorter wavelength near-infrared bands (1-2.4 $μ$m) of other surveys. This portion of the spectrum offers deep mass sensitivity, especially around nearby adolescent ($\sim$0.1-1 Gyr) stars. LEECH's contrast is competitive with other extreme adaptive optics systems, while providing an alternative survey strategy. Additionally, LEECH is characterizing known exoplanetary systems with observations from 3-5$μ$m in preparation for JWST.