Draft version May 11, 2016 A Preprint typeset using LTEX style emulateapj v. 5/2/11 IMAGING EXTRASOLAR GIANT PLANETS Brendan P. Bowler1,2 Draft version May 11, 2016 ABSTRACT High-contrast adaptive optics imaging is a powerful technique to probe the architectures of planetary systems from the outside-in and survey the atmospheres of self-luminous giant planets. Direct imaging has rapidly matured over the past decade and especially the last few years with the advent of high-order adaptive optics systems, dedicated planet-finding instruments with specialized coronagraphs, and innovative observing and post-processing strategies to suppress speckle noise. This review summarizes recent progress in high-contrast imaging with particular emphasis on observational results, discoveries near and below the deuterium-burning limit, and a practical overview of large-scale surveys and dedicated instruments. I conclude with a statistical meta-analysis of deep imaging surveys in the literature. Based on observations of 384 unique and single young (≈5–300 Myr) stars spanning stellar masses between 0.1–3.0 M⊙, the overall occurrence rate of 5–13 MJup companions at orbital distances +0.7 of 30–300 AU is 0.6−0.5% assuming hot-start evolutionary models. The most massive giant planets regularly accessible to direct imaging are about as rare as hot Jupiters are around Sun-like stars. Dividing this sample into individual stellar mass bins does not reveal any statistically-significant +3.7 trend in planet frequency with host mass: giant planets are found around 2.8−2.3% of BA stars, <4.1% of FGK stars, and <3.9% of M dwarfs. Looking forward, extreme adaptive optics systems and the next generation of ground- and space-based telescopes with smaller inner working angles and deeper detection limits will increase the pace of discovery to ultimately map the demographics, composition, evolution, and origin of planets spanning a broad range of masses and ages. Subject headings: planets and satellites: detection — planets and satellites: gaseous planets 1. INTRODUCTION Webb Space Telescope, WFIRST, and the thirty meter- Over the past two decades the orbital architecture of class telescopes. This “classical” period of high-contrast giant planets has expanded from a single order of mag- imaging spanning approximately 2000 to 2015 has set nitude in the Solar System (5–30 AU) to over five or- the stage and baseline expectations for the next gen- ders of magnitude among extrasolar planetary systems eration of instruments and telescopes that will deliver (0.01–5000 AU; Figure 1). High-contrast adaptive op- ultra-high contrasts and reach unprecedented sensitivi- tics (AO) imaging has played a critical role in this ad- ties. In addition to the first images of bona fide extraso- vancement by probing separations beyond ∼10 AU and lar planets, this early phase experienced a number of sur- prising discoveries including planetary-mass companions masses &1 MJup. Uncovering planetary-mass objects at hundreds and thousands of AU has fueled novel theo- orbiting brown dwarfs; planets on ultra-wide orbits be- ries of planet formation and migration, inspiring a more yond 100 AU; enigmatic (and still poorly understood) ob- complex framework for the origin of giant planets in jects like the optically-bright companion to Fomalhaut; which multiple mechanisms (core accretion, dynamical and unexpectedly red, methane-free, and dust-rich atmo- scattering, disk instability, and cloud fragmentation) op- spheres at low surface gravities. Among the most impor- erate on different timescales and orbital separations. In tant results has been the gradual realization that mas- addition to probing unexplored orbital distances, imag- sive planets are exceedingly rare on wide orbits; only a arXiv:1605.02731v1 [astro-ph.EP] 9 May 2016 ing entails directly capturing photons that originated in handful of discoveries have been made despite thousands planetary atmospheres, providing unparalleled informa- of hours spent on hundreds of targets spanning over a tion about the initial conditions, chemical composition, dozen surveys. Although dismaying, these null detec- internal structure, atmospheric dynamics, photospheric tions provide valuable information about the efficiency condensates, and physical properties of extrasolar plan- of planet formation and the resulting demographics at ets. These three science goals — the architecture, for- wide separations. Making use of mostly general, non- mation, and atmospheres of gas giants — represent the optimized facility instruments and early adaptive optics main motivations to directly image and spectroscopically systems has also led to creative observing strategies and characterize extrasolar giant planets. post-processing solutions for PSF subtraction. This pedagogical review summarizes the field of di- Distinguishing giant planets from low-mass brown rect imaging in the era leading up to and transition- dwarfs is a well-trodden intellectual exercise (e.g., ing towards extreme adaptive optics systems, the James Oppenheimeretal. 2000; Basri&Brown 2006; Chabrier et al. 2007; Chabrier et al. 2014). Except [email protected] in the few rare cases where the architectures or abun- 1 McDonald Observatory and the University of Texas at dance patterns of individual systems offer clues about Austin, Department of Astronomy, 2515 Speedway, Stop C1400, Austin, TX 78712 a specific formation route, untangling the origin of 2 McDonald Prize Fellow. imaged planetary-mass companions must necessarily 2 1996 2006 2016 75 ) Jup 10 Mass (M 1 0.01 0.1 1 10 100 0.01 0.1 1 10 100 0.01 0.1 1 10 100 1000 Separation (AU) Fig. 1.— Substellar companions discovered via radial velocities (gray circles) and direct imaging (red circles) as of 1996, 2006, and 2016. Over this twenty year period the number of directly imaged companions below 10 MJup has steadily increased from one (2M1207–3932 b in 2004) to over a dozen. The surprising discovery of planetary companions at extremely wide separations of hundreds to thousands of AU has expanded the architecture of planetary systems to over five orders of magnitude. Note that the radial velocity planets are minimum masses (msini) and the directly imaged companion masses are inferred from evolutionary models. RV-detected planets are from exoplanets.org (Wright et al. 2011; Han et al. 2014) and are supplemented with a compilation of RV-detected brown dwarfs from the literature. Imaged companions are from Deacon et al. (2014) together with other discoveries from the literature. be addressed as a population and in a statistical man- Marois et al. 2010a; Traub & Oppenheimer 2010; ner. This review is limited in scope to self-luminous Mawet et al. 2012b; Davies & Kasper 2012). companions detected in thermal emission at near- and 2. OPTIMAL TARGETS FOR HIGH-CONTRAST IMAGING mid-infrared wavelengths (1–5 µm) with masses between Planets radiatively cool over time by endlessly releas- ≈1–13 MJup with the understanding that multiple formation routes can probably produce objects in this ing the latent heat generated during their formation “planetary” mass regime (see Section 5). Indeed, the and gravitational contraction. Fundamental scaling rela- separations regularly probed in high-contrast imaging tions for the evolution of brown dwarfs and giant plan- surveys— typically tens to hundreds of AU— lie beyond ets can be derived analytically with basic assumptions the regions in protoplanetary disks containing the of a polytropic equation of state and degenerate electron highest surface densities of solids where core accretion gas (Stevenson 1991; Burrows & Liebert 1993). Neglect- operates most efficiently (e.g., Andrews 2015). Direct ing the influence of lithium burning, deuterium burning, imaging has therefore predominantly surveyed the wide and atmospheres, which act as partly opaque wavelength- orbital distances where alternative formation and migra- dependent boundary conditions, substellar objects with tion channels like disk instability, cloud fragmentation, different masses cool in a similar monotonic fashion over and planet-planet scattering are most likely to apply. time: In the future, the most efficient strategy to detect even −5/4 5/2 smaller super-Earths and terrestrial worlds close to their Lbol ∝ t M . (1) host stars will be in reflected light from a dedicated Here Lbol is the bolometric luminosity, t is the object’s space-based optical telescope. age, and M is its mass. By focusing on the optimal targets, early dis- This steep mass-luminosity relationship means that coveries, largest surveys, and statistical results, luminosity tracks are compressed in the brown dwarf this observationally-oriented overview aims to com- regime (≈13–75 MJup) and fan out in the planetary plement recent reviews on giant planet formation regime with significant consequences for high-contrast (Chabrier et al. 2007; Helled et al. 2013; Chabrier et al. imaging. A small gain in contrast in the brown dwarf 2014; Helling et al. 2014), atmospheric models regime results in a large gain in the limiting detectable (Marley et al. 2007a; Helling et al. 2008, Allard et al. mass, whereas the same contrast gain in the planetary 2012; Marley & Robinson 2015), evolutionary mod- regime has a much smaller influence on limiting mass els, (Burrows et al. 2001; Fortney & Nettelmann (Figure 2). It is much more difficult, for example, to 2009), observational results (Absil & Mawet 2009; improve sensitivity from 10 MJup to 1 MJup than from Lagrange 2014; Bailey 2014; Helling & Casewell 2014; 80 MJup to 10 MJup. Moreover, sensitivity to low masses Madhusudhan et al. 2014; Quanz 2015; Crossfield 2015), and close separations is highly dependent on a star’s and high contrast imaging instruments and speckle youth and proximity. In terms of limiting detectable suppression techniques (Guyon et al. 2006; Beuzit et al. planet mass, observing younger and closer stars is equiv- 2007; Oppenheimer & Hinkley 2009; Biller et al. 2008; alent to improving speckle suppression or integrating for 3 100 −10 mag 5 Gyr ) Jup 1 Gyr 500 Myr 10 100 pc 30 pc −5 mag 10 pc 100 Myr −2 mag NICI 30 Myr +2 mag Planet Mass (M 1 +5 mag 10 Myr 5 Myr 10 100 10 100 10 100 Semimajor Axis (AU) Fig.