Spatially Resolved Spectroscopy of the Exoplanet HR 8799 C∗
Total Page:16
File Type:pdf, Size:1020Kb
Spatially resolved spectroscopy of the exoplanet HR 8799 c∗ M. Janson1,4, C. Bergfors2, M. Goto2, W. Brandner2, D. Lafreni`ere3 [email protected] ABSTRACT HR 8799 is a multi-planet system detected in direct imaging, with three companions known so far. Here, we present spatially resolved VLT/NACO 3.88–4.10 µm spectroscopy of the middle planet, HR 8799 c, which has an estimated mass of ∼10 MJup, temperature of ∼1100 K and projected separation of 38 AU. The spectrum shows some differences in the continuum from existing theoretical models, particularly longwards of 4µm, implying that detailed cloud structure or non-equilibrium conditions may play an important role in the physics of young exoplanetary atmospheres. Subject headings: planetary systems — techniques: spectroscopic 1. Introduction spectral signal from the stellar one. During the past year, several planet candi- Spectroscopy of extrasolar planets will be a key dates from direct imaging have been reported (e.g. tool towards understanding their formation and Lafreniere et al. 2008, Kalas et al. 2008, Lagrange atmospheres, and could eventually lead to identifi- et al. 2009). Of particular interest here is the cation of chemical biomarkers (e.g. Kaltenegger et HR 8799 planetary system (Marois et al. 2008), al. 2007), thus addressing the question of life else- where the three companions have been confirmed where in the universe. However, the large bright- through common proper motion, orbital motion, ness contrast and small angular separation be- and detection of thermal radiation in several pho- tween a star and planet makes such spectroscopy tometric bands. This system is in many ways rem- technologically challenging. In the special geomet- iniscent of a scaled-up version of our own Solar ric case where the orbital plane of the planet is system. HR 8799 itself is an A5 star of 1.5 M . aligned along our line of sight, a spectroscopic sig- Sun HR 8799 b, c, and d have estimated masses of 7 nal can be temporally resolved by comparing the M , 10 M , and 10 M and separations of combined star and planet light in and out of tran- Jup Jup Jup 68 AU, 38 AU and 24 AU, respectively. The mass sit or secondary eclipse. This has been successfully estimates are based on an age estimate of 60 Myr done with the Spitzer and Hubble space telescopes for HR 8799. Although much higher ages are diffi- (e.g. Grillmair et al. 2007; Swain et al. 2009). cult to exclude categorically, several age indicators However, the vast majority of exoplanets do not (such as the massive debris disk) and arguments share this fortunate geometry. For these planets, from dynamical stability mutually indicate such a it is necessary to spatially resolve the planetary young age (e.g. Reidemeister et al. 2009). Mod- eling of infrared excess by Spitzer (e.g. Chen et *Based on observations collected at the European Southern Observatory, Chile (ESO No. 084.C-0072) and al. 2009) has shown that a Kuiper-analogous de- at the Subaru telescope, which is operated by the National bris belt exists outside of HR 8799 b, and a debris Astronomical Observatory of Japan. belt analogous to the Asteroid belt exists inside of 1 University of Toronto, Toronto, Canada HR 8799 d. The planets appear to be largely co- 2 Max-Planck-Institute for Astronomy, Heidelberg, Ger- planar with an orbital plane close to face-on, and many have low eccentricities. 3University of Montreal, Montreal, Canada 4Reinhardt fellow 1 2. Observations and Data Analysis to optimize the slit alignment. A 150 second se- quence without any high-contrast technique ap- Spectroscopic runs were performed with IRCS plied showed the presence the planet, confirming at Subaru and NACO at the VLT during the sec- that HR 8799 c indeed dominates the stellar PSF ond half of 2009 with the purpose of detecting the ′ at its separation (see Fig. 1). The slit was aligned spectroscopic L -band signature of the planets, in along the star-planet axis to ensure an adequate particular HR 8799 b and c. The IRCS run was slit positioning and a good PSF characterization of aimed primarily at HR 8799 b, and proved the the star, as well as provide an optimal correction feasibility of the technique, but the thermal back- of telluric and instrumental transmission effects. ground was much higher than predicted, probably The price to pay for this is some ghost images due to dust on the image rotator or the high am- on the detector from internal reflections. Pairwise bient humidity during the run. Hence, no planet nodding with 5” nod size was applied for back- spectroscopy could be extracted from the IRCS ground subtraction. Observations were also taken data. For the NACO run, one out of four half- with a slit rotation of 180o for PSF reference sub- nights was clouded out and two had thin cirrus traction. Flat field frames were produced with the clouds. During one of the nights with cirrus, HR slit and grism in, pointing at the thermal sky. 8799 b was included in the slit. Hence, some The data reduction was performed with ded- NACO data of HR 8799 b exists, but given the icated procedures written in IDL. The flat field ambient conditions we focused mainly on getting frame reproduced the expected vertical fringes a strong detection on HR 8799 c during the rest that occur in 4 µm spectroscopy. The fringe pat- of the time. The HR 8799 c spectroscopy is the tern was constant but shifted laterally with time, topic of this Letter. so the flat frame had to be digitally shifted to The technique employed was slit spectroscopy ′ each science frame. Since doing so would elimi- in the L -band range. This choice is based on nate pixel-to-pixel variation correction on the de- the combination of facts that the star-to-planet tector, both the image and flat frames first had to contrast is more favorable there (9.5 mag for HR be corrected for such variations. This was done 8799 c, as compared to 11.6 in H-band), and that with an imaging L′-band flat which was also ac- an exceptional adaptive optics performance can quired during the night. The spectral dependency be achieved in this range (upwards of 85 % at 4 of the pixel response was negligible, and so this µm, see Janson et al. 2008). For these reasons, a correction could be done satisfactorily. The shift planet like HR 8799 c dominates the stellar PSF ′ of the spectroscopic flat was determined individu- at its angular separation. A disadvantage of L - ally for each science frame by cross-correlation of band with respect to shorter wavelengths is the the fringe pattern with minimization of the resid- higher thermal background, but previous observa- uals. The flat- and dark-corrected frames were tions had shown that this was not a major limita- pairwise subtracted to remove the thermal back- tion. The widest NACO slit of 172 mas was used, ground. The spectrum of the star was traced to which allows for the companion PSF to be entirely an absolute center using a Gaussian fit on the PSF included in the slit, and along with the high PSF core individually for each wavelength, and rela- stability, allows to ensure that PSF chromaticity tive centering was checked using cross-correlation will have a negligible effect over the relevant wave- on the core with minimization of the residuals. A length range. The spectral resolution, limited by trace was then fitted over all wavelengths with 21 the FWHM, was approximately 900. The pixel pixel median box filtering, and all spectra were scale was 27 mas/pixel. Individual integration shifted to a common center. A collapse of all the times were 5 sec, and 20 coadds per frame were frames revealed the spectrum of the planet, at the taken. After excluding bad frames, 102 frames expected position and with the expected average were used for the target characterization and 94 brightness (see Fig. 2). The short-wavelength tail for the PSF characterization, giving total integra- of the planetary spectrum interfered with a resid- tion times of 10200 sec and 9400 sec, respectively. ual ghost feature, which could not be entirely re- ′ A quick-look L -band image of the HR 8799 moved since the PSF reference frames were not system was taken in the beginning of the night taken at identical dithering positions as the sci- 2 ence frames. Hence, we count this range as being wavelength end, and likewise a significant devia- of insufficient quality and use only the > 3.88 µm tion from the model. range where the planetary spectrum is clean. Hence, there are features in the spectrum that The star and planet spectra were extracted the COND model cannot explain. Part of the using 4-pixel apertures (108 mas, close to the problem is probably dust in the atmosphere – FWHM of the PSF). The background for the COND assumes dust-free conditions. However, a planet was estimated on the corresponding flux in DUSTY model (Chabrier et al. 2000) with the the PSF reference frame, averaged over 7 neigh- same range of parameters can also not reproduce boring pixels in the wavelength direction so as not the spectrum (see Fig. 4). Models with interme- to increase the background noise. The background diate abundances of dust have also been unable residual spectrum (non-averaged over wavelength) to reproduce the full set of observed properties of was also used to estimate the error by taking the the HR 8799 planets (Marois et al. 2008). The standard deviation over 21 neighboring pixels in results therefore imply that a more detailed treat- the wavelength direction.