Source author record

Christian Marois

Christian Marois 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

64works
4topics
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

64 published item(s)

preprint2022arXiv

Deep orbital search for additional planets in the HR 8799 system

The HR 8799 system hosts four massive planets orbiting 15 and 80 AU. Studies of the system's orbital stability and its outer debris disk open the possibility of additional planets, both interior to and exterior to the known system. Reaching a sufficient sensitivity to search for interior planets is very challenging due to the combination of bright quasi static speckle noise close to the stellar diffraction core and relatively fast orbital motion. In this work, we present a deep L-band imaging campaign using NIRC2 at Keck comprising 14 observing sequences. We further re-reduce archival data for a total of 16.75 hours, one of the largest uniform datasets of a single direct imaging target. Using a Bayesian modeling technique for detecting planets in images while compensating for plausible orbital motion, we then present deep limits on the existence of additional planets in the HR 8799 system. The final combination shows a tentative candidate, consistent with 4-7 $M_{jup}$ at 4-5 AU, detected with an equivalent false alarm probability better than $3σ$. This analysis technique is widely applicable to archival data and to new observations from upcoming missions that revisit targets at multiple epochs.

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

Focal Plane Wavefront Sensing with the FAST TGV Coronagraph

The continual push to directly image exoplanets at lower masses and closer separations orbiting around bright stars remains limited by both quasi-static and residual adaptive optics (AO) aberration. In previous papers we have proposed a modification of the self-coherent camera (SCC) design to address both of these limitations, called the Fast Atmospheric SCC Technique (FAST). In this paper we introduce an additional modification to the FAST focal plane mask design, including the existing Tip/tilt and Gaussian components and adding a charge four Vortex (TGV) component. In addition to boosting SCC fringe signal-to-noise ratio (S/N) as in our previous design, we show that the FAST TGV mask is also optimized to reach high contrast at separations closer to the star. In this paper we use numerical simulations to consider the performance improvement on correcting quasi-static aberration using this new mask compared to the previously proposed Tip/tilt+Gaussian mask. Using active deformable mirror control to generate a calibrated half dark hole improves contrast by a factor of about 200 at 2 - 5 $λ/D$ and up to a factor of 10 at 5 - 20 $λ/D$. The new methodology presented in this paper, now simultaneously considering both contrast and fringe S/N, opens the door to a new ideology of coronagraph design, where the coronagraph is now considered in duality as both a diffraction attenuator and a wavefront sensor.

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

preprint2020arXiv

The Gemini Planet Imager view of the HD 32297 debris disk

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

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

LRP2020: Signposts of planet formation in protoplanetary disks

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

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

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 of a Substellar Companion to the Nearby Debris Disk Host HR 2562

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

preprint2016arXiv

High contrast imaging of exoplanets on ELTs using a super-Nyquist wavefront control scheme

One of the key science goals for extremely large telescopes (ELTs) is the detailed characterization of already known directly imaged exoplanets. The typical adaptive optics (AO) Nyquist control region for ELTs is ~0.4 arcseconds, placing many already known directly imaged planets outside the DM control region and not allowing any standard wavefront control scheme to remove speckles that would allow higher SNR images/spectra to be acquired. This can be fixed with super-Nyquist wavefront control (SNWFC), using a sine wave phase plate to allow for wavefront control outside the central DM Nyquist region. We demonstrate that SNWFC is feasible through a simple, deterministic, non-coronagraphic, super-Nyquist speckle nulling technique in the adaptive optics laboratory at the National Research Council of Canada. We also present results in simulation of how SNWFC using the self coherent camera (SCC) can be used for high contrast imaging. This technique could be implemented on future high contrast imaging instruments to improve contrast outside the standard central dark hole for higher SNR characterization of exoplanets.

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

Large collaboration in observational astronomy: the Gemini Planet Imager exoplanet survey case

The Gemini Planet Imager (GPI) is a next-generation high-contrast imager built for the Gemini Observatory. The GPI exoplanet survey (GPIES) consortium is made up of 102 researchers from 28 institutions in North and South America and Europe. In November 2014, we launched a search for young Jovian planets and debris disks. In this paper, we discuss how we have coordinated the work done by this large team to improve the technical and scientific productivity of the campaign, and describe lessons we have learned that could be useful for future instrumentation-based astronomical surveys. The success of GPIES lies mostly on its decentralized structure, clear definition of policies that are signed by each member, and the heavy use of modern tools for communicating, exchanging information, and processing data.

preprint2016arXiv

Planet detection down to a few $λ$/D: an RSDI/TLOCI approach to PSF subtraction

Most current high contrast imaging point spread function (PSF) subtraction algorithms use some form of a least-squares noise minimization to find exoplanets that are, before post-processing, often hidden below the instrumental speckle noise. In the current standard PSF subtraction algorithms, a set of reference images is derived from the target image sequence to subtract each target image, using Angular and/or Simultaneous Spectral Differential Imaging (ADI, SSDI, respectively). However, to avoid excessive exoplanet self-subtraction, ADI and SSDI (in the absence of a strong spectral feature) severely limit the available number of reference images at small separations. This limits the performance of the least-squares algorithm, resulting in lower sensitivity to exoplanets at small angular separations. Possible solutions are to use additional reference images by acquiring longer sequences, use SSDI if the exoplanet is expected to show strong spectral features, or use images acquired on other targets. The latter option, known as Reference Star Differential Imaging (RSDI), which relies on the use of reference images that are highly correlated to the target image, has been ineffective in previous ground-based high contrast imaging surveys. We present the results of work to optimize PSF subtraction with the GPIES reference library using a least-squares algorithm designed to minimize speckle noise and maximize planet throughput, thus maximizing the planet signal to noise ratio (SNR). Using December 2014 51 Eri GPI data in the inner 100 mas to 300 mas annulus, we find no apparent improvement in SNR when using RSDI and/or our optimization scheme. This result, while still being investigated, seems to show that current algorithms on ADI+SSDI data sets are optimized, and that limited gains can be achieved by using a PSF archive.

preprint2016arXiv

Searching for the HR 8799 Debris Disk with HST/STIS

We present a new algorithm for space telescope high contrast imaging of close-to-face-on planetary disks called Optimized Spatially Filtered (OSFi) normalization. This algorithm is used on HR 8799 Hubble Space Telescope (HST) coronagraphic archival data, showing an over-luminosity after reference star point spread function (PSF) subtraction that may be from the inner disk and/or planetesimal belt components of this system. The PSF-subtracted radial profiles in two separate epochs from 2011 and 2012 are consistent with one another, and self-subtraction shows no residual in both epochs. We explore a number of possible false-positive scenarios that could explain this residual flux, including telescope breathing, spectral differences between HR 8799 and the reference star, imaging of the known warm inner disk component, OSFi algorithm throughput and consistency with the standard spider normalization HST PSF subtraction technique, and coronagraph misalignment from pointing accuracy. In comparison to another similar STIS dataset, we find that the over-luminosity is likely a result of telescope breathing and spectral difference between HR 8799 and the reference star. Thus, assuming a non-detection, we derive upper limits on the HR 8799 dust belt mass in small grains. In this scenario, we find that the flux of these micron-sized dust grains leaving the system due to radiation pressure is small enough to be consistent with measurements of other debris disk halos.

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 InfraRed Imaging Spectrograph (IRIS) for TMT: latest science cases and simulations

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

preprint2016arXiv

The International Deep Planet Survey II: The frequency of directly imaged giant exoplanets with stellar mass

Radial velocity and transit methods are effective for the study of short orbital period exoplanets but they hardly probe objects at large separations for which direct imaging can be used. We carried out the international deep planet survey of 292 young nearby stars to search for giant exoplanets and determine their frequency. We developed a pipeline for a uniform processing of all the data that we have recorded with NIRC2/Keck II, NIRI/Gemini North, NICI/Gemini South, and NACO/VLT for 14 years. The pipeline first applies cosmetic corrections and then reduces the speckle intensity to enhance the contrast in the images. The main result of the international deep planet survey is the discovery of the HR 8799 exoplanets. We also detected 59 visual multiple systems including 16 new binary stars and 2 new triple stellar systems, as well as 2,279 point-like sources. We used Monte Carlo simulations and the Bayesian theorem to determine that 1.05[+2.80-0.70]% of stars harbor at least one giant planet between 0.5 and 14M_J and between 20 and 300 AU. This result is obtained assuming uniform distributions of planet masses and semi-major axes. If we consider power law distributions as measured for close-in planets instead, the derived frequency is 2.30[+5.95-1.55]%, recalling the strong impact of assumptions on Monte Carlo output distributions. We also find no evidence that the derived frequency depends on the mass of the hosting star, whereas it does for close-in planets. The international deep planet survey provides a database of confirmed background sources that may be useful for other exoplanet direct imaging surveys. It also puts new constraints on the number of stars with at least one giant planet reducing by a factor of two the frequencies derived by almost all previous works.

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 PDS 66 Circumstellar Disk as seen in Polarized Light with the Gemini Planet Imager

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

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

Characterizing the Atmospheres of the HR8799 Planets with HST/WFC3

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

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

Simultaneous Detection of Water, Methane and Carbon Monoxide in the Atmosphere of Exoplanet HR8799b

Absorption lines from water, methane and carbon monoxide are detected in the atmosphere of exoplanet HR8799b. A medium-resolution spectrum presented here shows well-resolved and easily identified spectral features from all three molecules across the K band. The majority of the lines are produced by CO and H2O, but several lines clearly belong to CH4. Comparisons between these data and atmosphere models covering a range of temperatures and gravities yield log mole fractions of H2O between -3.09 and -3.91, CO between -3.30 and -3.72 and CH4 between -5.06 and -5.85. More precise mole fractions are obtained for each temperature and gravity studied. A reanalysis of H-band data, previously obtained at similar spectral resolution, results in a nearly identical water abundance as determined from the K-band spectrum. The methane abundance is shown to be sensitive to vertical mixing and indicates an eddy diffusion coefficient in the range of 10^6 to 10^8 cm^2 s^-1, comparable to mixing in the deep troposphere of Jupiter. The model comparisons also indicate a C/O between ~ 0.58 and 0.7, encompassing previous estimates for a second planet in the same system, HR8799c. Super-stellar C/O could indicate planet formation by core-accretion, however, the range of possible C/O for these planets (and the star) is currently too large to comment strongly on planet formation. More precise values of the bulk properties (e.g., effective temperature and surface gravity) are needed for improved abundance estimates.

preprint2014arXiv

Archival Legacy Investigations of Circumstellar Environments: Overview and First Results

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

preprint2014arXiv

Fundamental limitations of high contrast imaging set by small sample statistics

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

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 IX: Least-Squares Inversion Flux Extraction

The Gemini Planet Imager (GPI) is an instrument designed to directly image planets and circumstellar disks from 0.9 to 2.5 microns (the $YJHK$ infrared bands) using high contrast adaptive optics with a lenslet-based integral field spectrograph. We develop an extraction algorithm based on a least-squares method to disentangle the spectra and systematic noise contributions simultaneously. We utilize two approaches to adjust for the effect of flexure of the GPI optics which move the position of light incident on the detector. The first method is to iterate the extraction to achieve minimum residual and the second is to cross-correlate the detector image with a model image in iterative extraction steps to determine an offset. Thus far, this process has made clear qualitative improvements to the cube extraction by reducing the Moiré pattern. There are also improvements to the automated routines for finding flexure offsets which are reliable to with $\sim0.5$ pixel accuracy compared to pixel accuracy prior. Further testing and optimization will follow before implementation into the GPI pipeline.

preprint2014arXiv

Gemini Planet Imager Observational Calibrations V: Astrometry and Distortion

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

preprint2014arXiv

Gemini Planet Imager Observational Calibrations VI: Photometric and Spectroscopic Calibration for the Integral Field Spectrograph

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

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

Gemini Planet Imager Spectroscopy of the HR 8799 planets c and d

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

preprint2014arXiv

GPI PSF subtraction with TLOCI: the next evolution in exoplanet/disk high-contrast imaging

To directly image exoplanets and faint circumstellar disks, the noisy stellar halo must be suppressed to a high level. To achieve this feat, the angular differential imaging observing technique and the least-squares Locally Optimized Combination of Images (LOCI) algorithm have now become the standard in single band direct imaging observations and data reduction. With the development and commissioning of new high-order high-contrast adaptive optics equipped with integral field units, the image subtraction algorithm needs to be modified to allow the optimal use of polychromatic images, field-rotated images and archival data. A new algorithm, TLOCI (for Template LOCI), is designed to achieve this task by maximizing a companion signal-to-noise ratio instead of simply minimizing the noise as in the original LOCI algorithm. The TLOCI technique uses an input spectrum and template Point Spread Functions (PSFs, generated from unocculted and unsaturated stellar images) to optimize the reference image least-squares coefficients to minimize the planet self-subtraction, thus maximizing its throughput per wavelength, while simultaneously providing a maximum suppression of the speckle noise. The new algorithm has been developed using on-sky GPI data and has achieved impressive contrast. This paper presents the TLOCI algorithm, on-sky performance, and will discuss the challenges in recovering the planet spectrum with high fidelity.

preprint2014arXiv

Near-Infrared Detection and Characterization of the Exoplanet HD 95086 b with the Gemini Planet Imager

HD 95086 is an intermediate-mass debris-disk-bearing star. VLT/NaCo $3.8 μm$ observations revealed it hosts a $5\pm2 \mathrm{M}_{Jup}$ companion (HD 95086 b) at $\simeq 56$ AU. Follow-up observations at 1.66 and 2.18 $μm$ yielded a null detection, suggesting extremely red colors for the planet and the need for deeper direct-imaging data. In this Letter, we report H- ($1.7 μm$) and $\mathrm{K}_1$- ($2.05 μm$) band detections of HD 95086 b from Gemini Planet Imager (GPI) commissioning observations taken by the GPI team. The planet position in both spectral channels is consistent with the NaCo measurements and we confirm it to be comoving. Our photometry yields colors of H-L'= $3.6\pm 1.0$ mag and K$_1$-L'=$2.4\pm 0.7$ mag, consistent with previously reported 5-$σ$ upper limits in H and Ks. The photometry of HD 95086 b best matches that of 2M 1207 b and HR 8799 cde. Comparing its spectral energy distribution with the BT-SETTL and LESIA planet atmospheric models yields T$_{\mathrm{eff}}\sim$600-1500 K and log g$\sim$2.1-4.5. Hot-start evolutionary models yield M=$5\pm2$ M$_{Jup}$. Warm-start models reproduce the combined absolute fluxes of the object for M=4-14 M$_{Jup}$ for a wide range of plausible initial conditions (S$_{init}$=8-13 k$_{B}$/baryon). The color-magnitude diagram location of HD 95086 b and its estimated T$_{\mathrm{eff}}$ and log g suggest that the planet is a peculiar L-T transition object with an enhanced amount of photospheric dust.

preprint2014arXiv

Polarimetry with the Gemini Planet Imager: Methods, Performance at First Light, and the Circumstellar Ring around HR 4796A

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

preprint2014arXiv

The First H-band Spectrum of the Massive Gas Giant Planet beta Pictoris b with the Gemini Planet Imager

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

preprint2014arXiv

The Gemini Planet Imager: First Light

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

preprint2014arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: Overview of innovative science programs

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

preprint2013arXiv

Detection of Carbon Monoxide and Water Absorption Lines in an Exoplanet Atmosphere

Determining the atmospheric structure and chemical composition of an exoplanet remains a formidable goal. Fortunately, advancements in the study of exoplanets and their atmospheres have come in the form of direct imaging - spatially resolving the planet from its parent star - which enables high-resolution spectroscopy of self-luminous planets in Jovian-like orbits. Here, we present a spectrum with numerous, well-resolved, molecular lines from both water and carbon monoxide from a massive planet orbiting less than 40 AU from the star HR 8799. These data reveal the planet's chemical composition, atmospheric structure, and surface gravity, confirming that it is indeed a young planet. The spectral lines suggest an atmospheric carbon-to-oxygn ratio greater than the host star's, providing hints about the planet's formation.

preprint2013arXiv

Fomalhaut b: Independent Analysis of the Hubble Space Telescope Public Archive Data

The nature and even the existence of a putative planet-mass companion ("Fomalhaut b") to Fomalhaut has been debated since 2008. In the present paper we reanalyze the multi-epoch ACS/STIS/WFC3 Hubble Space Telescope (HST) optical/near infrared images on which the discovery and some other claims were based. We confirm that the HST images do reveal an object in orbit around Fomalhaut but the detailed results from our analysis differ in some ways from previous discussions. In particular, we do not confirm flux variability over a two-year interval at 0.6 microns wavelength and we detect Fomalhaut b for the first time at the short wavelength of 0.43 microns. We find that the HST image of Fomalhaut b at $0.8\,μ$m may be extended beyond the PSF. We cannot determine from our astrometry if Fomalhaut b will cross or not the dust ring. The optical through mid-infrared spectral energy distribution (SED) of Fomalhaut b cannot be explained as due to direct or scattered radiation from a massive planet. We consider two models to explain the SED: (1) a large circumplanetary disk around an unseen planet and (2) the aftermath of a collision during the past 50-150 years of two Kuiper Belt-like objects of radii 50 km.

preprint2013arXiv

Resolved Imaging of the HR 8799 Debris Disk with Herschel

We present Herschel far-infrared and submillimeter maps of the debris disk associated with the HR 8799 planetary system. We resolve the outer disk emission at 70, 100, 160 and 250 um and detect the disk at 350 and 500 um. A smooth model explains the observed disk emission well. We observe no obvious clumps or asymmetries associated with the trapping of planetesimals that is a potential consequence of planetary migration in the system. We estimate that the disk eccentricity must be <0.1. As in previous work by Su et al. (2009), we find a disk with three components: a warm inner component and two outer components, a planetesimal belt extending from 100 - 310 AU, with some flexibility (+/- 10 AU) on the inner edge, and the external halo which extends to ~2000 AU. We measure the disk inclination to be 26 +/- 3 deg from face-on at a position angle of 64 deg E of N, establishing that the disk is coplanar with the star and planets. The SED of the disk is well fit by blackbody grains whose semi-major axes lie within the planetesimal belt, suggesting an absence of small grains. The wavelength at which the spectrum steepens from blackbody, 47 +/- 30 um, however, is short compared to other A star debris disks, suggesting that there are atypically small grains likely populating the halo. The PACS longer wavelength data yield a lower disk color temperature than do MIPS data (24 and 70 um), implying two distinct halo dust grain populations.

preprint2012arXiv

High-contrast imaging in the Hyades with snapshot LOCI

To image faint substellar companions obscured by the stellar halo and speckles, scattered light from the bright primary star must be removed in hardware or software. We apply the "locally-optimized combination of images" (LOCI) algorithm to 1-minute Keck Observatory snapshots of GKM dwarfs in the Hyades using source diversity to determine the most likely PSF. We obtain a mean contrast of 10^{-2} at 0.01", 10^{-4} at <1", and 10^{-5} at 5". New brown dwarf and low-mass stellar companions to Hyades primaries are found in a third of the 84 targeted systems. This campaign shows the efficacy of LOCI on snapshot imaging as well as on bright wide binaries with off-axis LOCI, reaching contrasts sufficient for imaging 625-Myr late-L/early-T dwarfs purely in post-processing.

preprint2012arXiv

Review of small-angle coronagraphic techniques in the wake of ground-based second-generation adaptive optics systems

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

preprint2012arXiv

The VAST Survey -- II. Orbital motion monitoring of A-type star multiples

As a part of our ongoing Volume-limited A-Star (VAST) adaptive optics survey, we have obtained observations of 26 binary systems with projected separations <100 AU, 13 of which have sufficient historical measurements to allow for refinement of their orbital elements. For each system with an estimated orbit, the dynamical system mass obtained was compared with the system mass estimated from mass-magnitude relations. Discrepancies between the dynamical and theoretical system mass can be explained by the presence of a previously unresolved spectroscopic component, or by a non-solar metallicity of the system. Using this approach to infer the presence of additional companions, a lower limit to the fraction of binaries, triples, and quadruples can be estimated as 39, 46, and 15 per cent, for systems with at least one companion within 100 AU. The fraction of multiple systems with three or more components shows a relative increase compared to the fraction for Solar-type primaries resolved in previous volume-limited surveys. The observations have also revealed a pair of potentially young ($<$100 Myr) M-dwarf companions, which would make an ideal benchmark for the theoretical models during the pre-Main Sequence contraction phase for M-dwarfs. In addition to those systems with orbit fits, we report 13 systems for which further orbital monitoring observations are required, 11 of which are newly resolved as a part of the VAST survey.

preprint2011arXiv

A Substellar Companion to the Dusty Pleiades Star HD 23514

With adaptive optics imaging at Keck observatory, we have discovered a substellar companion to the F6 Pleiades star HD 23514, one of the dustiest main-sequence stars known to date (L_{IR}/L_{*}~2%). This is one of the first brown dwarfs discovered as a companion to a star in the Pleiades. The 0.06 M$_\odot$ late-M secondary has a projected separation of ~360 AU. The scarcity of substellar companions to stellar primaries in the Pleiades combined with the extremely dusty environment make this a unique system to study.

preprint2011arXiv

Clouds and Chemistry in the Atmosphere of Extrasolar Planet HR8799b

Using the integral field spectrograph OSIRIS, on the Keck II telescope, broad near-infrared H and K-band spectra of the young exoplanet HR8799b have been obtained. In addition, six new narrow-band photometric measurements have been taken across the H and K bands. These data are combined with previously published photometry for an analysis of the planet's atmospheric properties. Thick photospheric dust cloud opacity is invoked to explain the planet's red near-IR colors and relatively smooth near-IR spectrum. Strong water absorption is detected, indicating a Hydrogen-rich atmosphere. Only weak CH4 absorption is detected at K band, indicating efficient vertical mixing and a disequilibrium CO/CH4 ratio at photospheric depths. The H-band spectrum has a distinct triangular shape consistent with low surface gravity. New giant planet atmosphere models are compared to these data with best fitting bulk parameters, Teff = 1100K +- 100 and log(g) = 3.5 +-0.5 (for solar composition). Given the observed luminosity, log(Lobs/Lsun) ~ -5.1, these values correspond to a radius of 0.75 Rjup (+0.17, -0.12) and mass ~ 0.72 Mjup (+2.6,-0.6) -- strikingly inconsistent with interior/evolution models. Enhanced metallicity (up to ~ 10 times that of the Sun) along with thick clouds and non-equilibrium chemistry are likely required to reproduce the complete ensemble of spectroscopic and photometric data and the low effective temperatures (< 1000K) required by the evolution models.

preprint2011arXiv

Determination of the inclination of the multi-planet hosting star HR8799 using asteroseismology

Direct imaging of the HR8799 system was a major achievement in the study of exoplanets. HR8799 is a $γ$\,Doradus variable and asteroseismology can provide an independent constraint on the inclination. Using 650 high signal-to-noise, high resolution, full visual wavelength spectroscopic observations obtained over two weeks at Observatoire de Haute Provence (OHP) with the SOPHIE spectrograph we find that the main frequency in the radial velocity data is 1.9875 d$^{-1}$. This frequency corresponds to the main frequency as found in previous photometric observations. Using the FAMIAS software to identify the pulsation modes, we find this frequency is a prograde $\ell$=1 sectoral mode and obtain the constraint that inclination $i\gtrsim$40$^{\circ}$.

preprint2011arXiv

Experimental Design for the Gemini Planet Imager

The Gemini Planet Imager (GPI) is a high performance adaptive optics system being designed and built for the Gemini Observatory. GPI is optimized for high contrast imaging, combining precise and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with integral field and polarimetry capabilities. The primary science goal for GPI is the direct detection and characterization of young, Jovian-mass exoplanets. For plausible assumptions about the distribution of gas giant properties at large semi-major axes, GPI will be capable of detecting more than 10% of gas giants more massive than 0.5 M_J around stars younger than 100 Myr and nearer than 75 parsecs. For systems younger than 1 Gyr, gas giants more massive than 8 M_J and with semi-major axes greater than 15 AU are detected with completeness greater than 50%. A survey targeting young stars in the solar neighborhood will help determine the formation mechanism of gas giant planets by studying them at ages where planet brightness depends upon formation mechanism. Such a survey will also be sensitive to planets at semi-major axes comparable to the gas giants in our own solar system. In the simple, and idealized, situation in which planets formed by either the "hot-start" model of Burrows et al. (2003) or the core accretion model of Marley et al. (2007), a few tens of detected planets are sufficient to distinguish how planets form.

preprint2011arXiv

M-band Imaging of the HR 8799 Planetary System Using an Innovative LOCI-based Background Subtraction Technique

Multi-wavelength observations/spectroscopy of exoplanetary atmospheres are the basis of the emerging exciting field of comparative exoplanetology. The HR 8799 planetary system is an ideal laboratory to study our current knowledge gap between massive field brown dwarfs and the cold 5-Gyr old Solar system planets. The HR 8799 planets have so far been imaged at J- to L-band, with only upper limits available at M-band. We present here deep high-contrast Keck II adaptive optics M-band observations that show the imaging detection of 3 of the 4 currently known HR 8799 planets. Such detections were made possible due to the development of an innovative LOCI-based background subtraction scheme that is 3 times more efficient than a classical median background subtraction for Keck II AO data, representing a gain in telescope time of up to a factor of 9. These M-band detections extend the broad band photometric coverage out to about 5 microns and provide access to the strong CO fundamental absorption band at 4.5microns. The new M-band photometry shows that the HR 8799 planets are located near the L/T-type dwarf transition, similar to what was found by other studies. We also confirm that the best atmospheric fits are consistent with low surface gravity, dusty and non-equilibrium CO/CH4 chemistry models.

preprint2011arXiv

The VAST Survey - I. Companions and the unexpected X-ray detection of B6-A7 stars

With an adaptive optics imaging survey of 148 B6-A7 stars, we have tested the hypothesis that unresolved lower-mass companions are the source of the unexpected X-ray detections of stars in this spectral type range. The sample is composed of 63 stars detected in X-rays within the ROSAT All-Sky Survey and 85 stars that form a control sample; both subsets have the same restricted distribution of spectral type, age, X-ray sensitivity and separation coverage. A total of 68 companion candidates are resolved with separations ranging from 0.3" to 26.2", with 23 new detections. The multiple star frequency of the X-ray sample based on companions resolved within the ROSAT error ellipse is found to be 43 (+6,-6)%. The corresponding control sample multiple star frequency is three times lower at 12 (+4,-3)% -- a difference of 31\pm7%. These results are presented in the first of a series of papers based on our Volume-limited A-Star (VAST) survey -- a comprehensive study of the multiplicity of A-type stars.

preprint2011arXiv

The Young Planet-mass Object 2M1207b: a cool, cloudy, and methane-poor atmosphere

The properties of 2M1207b, a young (~ 8 Myr) planet-mass companion, have lacked a satisfactory explanation for some time. The combination of low luminosity, red near-IR colors, and L-type near-IR spectrum (previously consistent with Teff ~ 1600K) implies an abnormally small radius. Early explanations for the apparent underluminosity of 2M1207b invoked an edge-on disk or the remnant of a recent protoplanetary collision. The discovery of a second planet-mass object (HR8799b) with similar luminosity and colors as 2M1207b indicates that a third explanation, one of a purely atmospheric nature, is more likely. By including clouds, non-equilibrium chemistry, and low-gravity, an atmosphere with effective temperature consistent with evolution cooling-track predictions is revealed. Consequently, 2M1207b, and others like it, requires no new physics to explain nor do they belong to a new class of objects. Instead they most likely represent the natural extension of cloudy substellar atmospheres down to low Teff and log(g). If this atmosphere only explanation for 2M1207b is correct, then very young planet-mass objects with near-IR spectra similar to field T dwarfs may be rare.

preprint2010arXiv

The highest resolution near-IR spectrum of the imaged planetary mass companion 2M1207 b

Direct-imaging searches for planets reveal wide orbit planets amenable to spectroscopy, and their atmospheres represent an important comparison to the irradiated atmospheres of Hot Jupiters. Using AO integral field spectroscopy of 2M1207 b, the shape of the continuum emission over the J, H, and K bands from the atmosphere of this young, planetary mass companion is measured in order to compare with atmospheric and evolutionary models, and objects of similar temperature in young clusters and the field. The 2M1207 b spectrum has the highest spectral resolution (R~300-1500) and largest wavelength coverage, including the first J-band spectrum, for this benchmark object. The high signal-to-noise of the data allow a clear identification of signatures of low surface gravity, and comparison with a grid of AMES-Dusty models reveals a best-fit effective temperature of Teff=1600 K with a preferred surface gravity of log g=4.5. The J-band flux is depressed relative to nearly all L-type objects, and the detailed shape of the absorption features across the H-band exhibit differences from the model predictions. The possible origins of 2M1207 b and its low luminosity are examined with the new data and analysis which suggest that extinction from a disk with large grains is a viable scenario and is preferred over scatttering off an optically thick disk. The 2M1207 b spectrum presents an important comparison for the types of features which may be present in upcoming spectra of the atmospheres of planets imaged in orbit around stellar primaries.

preprint2009arXiv

Near Infrared Observations of GQ Lup b Using the Gemini Integral Field Spectrograph NIFS

We present new JHK spectroscopy (R ~ 5000) of GQ Lup b, acquired with the near-infrared integral field spectrograph NIFS and the adaptive optics system ALTAIR at the Gemini North telescope. Angular differential imaging was used in the J and H bands to suppress the speckle noise from GQ Lup A; we show that this approach can provide improvements in signal-to-noise ratio (S/N) by a factor of 2 - 6 for companions located at subarcsecond separations. Based on high quality observations and GAIA synthetic spectra, we estimate the companion effective temperature to Teff = 2400 +/- 100 K, its gravity to log g = 4.0 +/- 0.5, and its luminosity to log(L/L_s) = -2.47 +/- 0.28. Comparisons with the predictions of the DUSTY evolutionary tracks allow us to constrain the mass of GQ Lup b to 8 - 60 MJup, most likely in the brown dwarf regime. Compared with the spectra published by Seifahrt and collaborators, our spectra of GQ Lup b are significantly redder (by 15 - 50%) and do not show important Paβemission. Our spectra are in excellent agreement with the lower S/N spectra previously published by McElwain and collaborators.

preprint2009arXiv

Pre-Discovery 2007 Image of the HR 8799 Planetary System

We present a pre-discovery H-band image of the HR 8799 planetary system that reveals all three planets in August 2007. The data were obtained with the Keck adaptive optics system, using angular differential imaging and a coronagraph. We confirm the physical association of all three planets, including HR 8799d, which had only been detected in 2008 images taken two months apart, and whose association with HR 8799 was least secure until now. We confirm that the planets are 2-3 mag fainter than field brown dwarfs of comparable near-infrared colors. We note that similar under-luminosity is characteristic of young substellar objects at the L/T spectral type transition, and is likely due to enhanced dust content and non-equilibrium CO/CH_4 chemistry in their atmospheres. Finally, we place an upper limit of 18 mag per square arc second on the >120 AU H-band dust-scattered light from the HR 8799 debris disk. The upper limit on the integrated scattered light flux is 1e-4 times the photospheric level, 24 times fainter than the debris ring around HR 4796A.

preprint2007arXiv

A new algorithm for point spread function subtraction in high-contrast imaging: a demonstration with angular differential imaging

Direct imaging of exoplanets is limited by bright quasi-static speckles in the point spread function (PSF) of the central star. This limitation can be reduced by subtraction of reference PSF images. We have developed an algorithm to construct an optimized reference PSF image from a set of reference images. This image is built as a linear combination of the reference images available and the coefficients of the combination are optimized inside multiple subsections of the image independently to minimize the residual noise within each subsection. The algorithm developed can be used with many high-contrast imaging observing strategies relying on PSF subtraction, such as angular differential imaging (ADI), roll subtraction, spectral differential imaging, reference star observations, etc. The performance of the algorithm is demonstrated for ADI data. It is shown that for this type of data the new algorithm provides a gain in sensitivity by up to a factor 3 at small separation over the algorithm used in Marois et al. (2006).

preprint2007arXiv

Improving the speckle noise attenuation of simultaneous spectral differential imaging with a focal plane holographic diffuser

Direct exoplanet detection is limited by speckle noise in the point spread function (PSF) of the central star. This noise can be reduced by subtracting PSF images obtained simultaneously in adjacent narrow spectral bands using a multi-channel camera (MCC), but only to a limit imposed by differential optical aberrations in the MCC. To alleviate this problem, we suggest the introduction of a holographic diffuser at the focal plane of the MCC to convert the PSF image into an incoherent illumination scene that is then re-imaged with the MCC. The re-imaging is equivalent to a convolution of the scene with the PSF of each spectral channel of the camera. Optical aberrations in the MCC affect only the convolution kernel of each channel and not the PSF globally, resulting in better correlated images. We report laboratory measurements with a dual channel prototype (1.575 micron and 1.625 micron) to validate this approach. A speckle noise suppression factor of 12-14 was achieved, an improvement by a factor ~5 over that obtained without the holographic diffuser. Simulations of realistic exoplanet populations for three representative target samples show that the increase in speckle noise attenuation achieved in the laboratory would roughly double the number of planets that could be detected with current adaptive optics systems on 8-m telescopes.