# The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics from 10 to 100 au

Eric L. Nielsen, Robert J. de Rosa, Bruce Macintosh, Jason J. Wang, Jean-Baptiste Ruffio, Eugene Chiang, Mark S. Marley, Didier Saumon, Dmitry Savransky, S. Mark Ammons, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Eric L. Nielsen, Robert J. de Rosa, Bruce Macintosh, Jason J. Wang, Jean-Baptiste Ruffio, et al. The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics from 10 to 100 au. The Astronomical Journal, 2019, 158, 10.3847/1538-3881/ab16e9. insu-03703851

HAL Id: insu-03703851

https://insu.hal.science/insu-03703851v1

Submitted on 4 Aug 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics from 10 to 100 au

Eric L. Nielsen<fnref n="1" /> ID, Robert J. De Rosa<fnref n="1" /> ID, Bruce Macintosh<fnref n="1" /> ID, Jason J. Wang<fnref n="2" />,<fnref n="3" />,<fnref n="43" /> ID, Jean-Baptiste Ruffio<fnref n="1" /> ID, Eugene Chiang<fnref n="3" /> ID, Mark S. Marley<fnref n="4" /> ID, Didier Saumon<fnref n="5" /> ID, Dmitry Savransky<fnref n="6" /> ID, S. Mark Ammons<fnref n="7" /> ID, Vanessa P. Bailey<fnref n="8" /> ID, Travis Barman<fnref n="9" /> ID, Célia Blain<fnref n="10" />, Joanna Bulger<fnref n="11" /> ID, Adam Burrows<fnref n="12" /> ID, Jeffrey Chilcote<fnref n="1" />,<fnref n="13" /> ID, Tara Cotten<fnref n="14" /> ID, Ian Czekala<fnref n="1" />,<fnref n="3" />,<fnref n="44" /> ID, Rene Doyon<fnref n="15" />, Gaspard Duchêne<fnref n="3" />,<fnref n="16" /> ID, Thomas M. Esposito<fnref n="3" /> ID, Daniel Fabrycky<fnref n="17" /> ID, Michael P. Fitzgerald<fnref n="18" /> ID, Katherine B. Follette<fnref n="19" /> ID, Jonathan J. Fortney<fnref n="20" /> ID, Benjamin L. Gerard<fnref n="10" />,<fnref n="21" /> ID, Stephen J. Goodsell<fnref n="22" /> ID, James R. Graham<fnref n="3" />, Alexandra Z. Greenbaum<fnref n="23" /> ID, Pascale Hibon<fnref n="24" /> ID, Sasha Hinkley<fnref n="25" /> ID, Lea A. Hirsch<fnref n="1" />, Justin Hom<fnref n="26" />, Li-Wei Hung<fnref n="27" /> ID, Rebekah Ilene Dawson<fnref n="28" /> ID, Patrick Ingraham<fnref n="29" /> ID, Paul Kalas<fnref n="3" />,<fnref n="30" /> ID, Quinn Konopacky<fnref n="31" /> ID, James E. Larkin<fnref n="18" /> ID, Eve J. Lee<fnref n="32" /> ID, Jonathan W. Lin<fnref n="3" />, Jérôme Maire<fnref n="31" />, Franck Marchis<fnref n="30" /> ID, Christian Marois<fnref n="10" />,<fnref n="21" /> ID, Stanimir Metchev<fnref n="33" />,<fnref n="34" /> ID, Maxwell A. Millar-Blanchaer<fnref n="8" />,<fnref n="42" /> ID, Katie M. Morzinski<fnref n="35" /> ID, Rebecca Oppenheimer<fnref n="36" /> ID, David Palmer<fnref n="7" />, Jennifer Patience<fnref n="26" />, Marshall Perrin<fnref n="37" /> ID, Lisa Poyneer<fnref n="7" />, Laurent Pueyo<fnref n="37" />, Roman R. Rafikov<fnref n="38" /> ID, Abhijith Rajan<fnref n="37" /> ID, Julien Rameau<fnref n="15" /> ID, Fredrik T. Rantakyrö<fnref n="39" /> ID, Bin Ren<fnref n="40" /> ID, Adam C. Schneider<fnref n="26" /> ID, Anand Sivaramakrishnan<fnref n="37" /> ID, Inseok Song<fnref n="14" /> ID, Remi Soummer<fnref n="37" /> ID, Melisa Tallis<fnref n="1" /> ID, Sandrine Thomas<fnref n="29" /> ID, Kimberly Ward-Duong<fnref n="26" /> ID, and Schuyler Wolff<fnref n="41" /> ID

## Abstract

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey. This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semimajor axis, and host

<footnote n="1">¹ Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA</footnote>

<footnote n="2">² Department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA</footnote>

<footnote n="3">³ Department of Astronomy, University of California, Berkeley, CA 94720, USA</footnote>

<footnote n="4">⁴ NASA Ames Research Center, Mountain View, CA 94035, USA</footnote>

<footnote n="5">⁵ Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM 87545, USA</footnote>

<footnote n="6">⁶ Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY 14853, USA</footnote>

<footnote n="7">⁷ Lawrence Livermore National Laboratory, Livermore, CA 94551, USA</footnote>

<footnote n="8">⁸ Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA</footnote>

<footnote n="9">⁹ Lunar and Planetary Laboratory, University of Arizona, Tucson AZ 85721, USA</footnote>

<footnote n="10">¹⁰ National Research Council of Canada Herzberg, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada</footnote>

<footnote n="11">¹¹ Subaru Telescope, NAOJ, 650 North A'ohoku Place, Hilo, HI 96720, USA</footnote>

<footnote n="12">¹² Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA</footnote>

<footnote n="13">¹³ Department of Physics, University of Notre Dame, 225 Nieuwland Science Hall, Notre Dame, IN, 46556, USA</footnote>

<footnote n="14">¹⁴ Department of Physics and Astronomy, University of Georgia, Athens, GA 30602, USA</footnote>

<footnote n="15">¹⁵ Institut de Recherche sur les Exoplanètes, Département de Physique, Université de Montréal, Montréal QC, H3C 3J7, Canada</footnote>

<footnote n="16">¹⁶ Univ. Grenoble Alpes/CNRS, IPAG, F-38000 Grenoble, France</footnote>

<footnote n="17">¹⁷ Department of Astronomy & Astrophysics University of Chicago 5640 S. Ellis Avenue, Chicago, IL 60615, USA</footnote>

<footnote n="18">¹⁸ Department of Physics & Astronomy, University of California, Los Angeles, CA 90095, USA</footnote>

<footnote n="19">¹⁹ Physics and Astronomy Department, Amherst College, 21 Merrill Science Drive, Amherst, MA 01002, USA</footnote>

<footnote n="20">²⁰ Department of Astronomy, UC Santa Cruz, 1156 High Street, Santa Cruz, CA 95064, USA</footnote>

<footnote n="21">²¹ University of Victoria, 3800 Finnerty Road, Victoria, BC, V8P 5C2, Canada</footnote>

<footnote n="22">²² Gemini Observatory, 670 N. A'ohoku Place, Hilo, HI 96720, USA</footnote>

<footnote n="23">²³ Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA</footnote>

<footnote n="24">²⁴ European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile</footnote>

<footnote n="25">²⁵ University of Exeter, Astrophysics Group, Physics Building, Stocker Road, Exeter, EX4 4QL, UK</footnote>

<footnote n="26">²⁶ School of Earth and Space Exploration, Arizona State University, P.O. Box 871404, Tempe, AZ 85287, USA</footnote>

<footnote n="27">²⁷ Natural Sounds and Night Skies Division, National Park Service, Fort Collins, CO 80525, USA</footnote>

<footnote n="28">²⁸ Department of Astronomy & Astrophysics, Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA 16802, USA</footnote>

<footnote n="29">²⁹ Large Synoptic Survey Telescope, 950N Cherry Avenue, Tucson, AZ 85719, USA</footnote>

<footnote n="30">³⁰ SETI Institute, Carl Sagan Center, 189 Bernardo Avenue, Mountain View CA 94043, USA</footnote>

<footnote n="31">³¹ Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA 92093, USA</footnote>

<footnote n="32">³² TAPIR, Walter Burke Institute for Theoretical Physics, Mailcode 350-17, Caltech, Pasadena, CA 91125, USA</footnote>

<footnote n="33">³³ Department of Physics and Astronomy, Centre for Planetary Science and Exploration, The University of Western Ontario, London, ON N6A 3K7, Canada</footnote>

<footnote n="34">³⁴ Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA</footnote>

<footnote n="35">³⁵ Steward Observatory, 933 N. Cherry Avenue, University of Arizona, Tucson, AZ 85721, USA</footnote>

<footnote n="36">³⁶ Department of Astrophysics, American Museum of Natural History, New York, NY 10024, USA</footnote>

<footnote n="37">³⁷ Space Telescope Science Institute, Baltimore, MD 21218, USA</footnote>

<footnote n="38">³⁸ Centre for Mathematical Sciences, Department of Applied Mathematics and Theoretical Physics, University of Cambridge, UK</footnote>

<footnote n="39">³⁹ Gemini Observatory, Casilla 603, La Serena, Chile</footnote>

<footnote n="40">⁴⁰ Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA</footnote>

<footnote n="41">⁴¹ Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands</footnote>

<footnote>Received 2019 January 15; revised 2019 March 14; accepted 2019 April 2; published 2019 June 12</footnote>

<footnote n="42">⁴² NASA Hubble Fellow.</footnote>

<footnote n="43">⁴³ 51 Pegasi b Fellow.</footnote>

<footnote n="44">⁴⁴ NASA Hubble Fellowship Program Sagan Fellow.</footnote>
stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars $M_* > 1.5 M_\odot$ more likely to host planets with masses between 2 and $13 M_{\text{Jup}}$ and semimajor axes of 3–100 au at 99.92% confidence. We fit a double power-law model in planet mass ($m$) and semimajor axis ($a$) for planet populations around high-mass stars ($M_* > 1.5 M_\odot$) of the form $d^2N/(dm\,da) \propto m^\alpha a^\beta$, finding $\alpha = -2.4 \pm 0.8$ and $\beta = -2.0 \pm 0.5$, and an integrated occurrence rate of $9_{-4}^{+5}\%$ between 5–$13 M_{\text{Jup}}$ and 10–100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with $0.8_{-0.5}^{+0.8}\%$ of stars hosting a brown dwarf companion between 13–80$M_{\text{Jup}}$ and 10–100 au. Brown dwarfs also appear to be distributed differently in mass and semimajor axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semimajor axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the radial velocity method, our results are consistent with a peak in occurrence of giant planets between ~1 and 10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

**Key words:** instrumentation: adaptive optics – planetary systems – planets and satellites: detection

*Supporting material:* machine-readable table

# 1. Introduction

Since giant extrasolar planets ($\gtrsim 1 M_{\text{Jup}}$) are typically the easiest to detect for most search techniques, there is a long record of studies that investigate the characteristics of the underlying population. A principal scientific goal is to distinguish between different planet formation and evolution theories by measuring planet occurrence rates with respect to the properties of the planets (e.g., mass, semimajor axis, eccentricity) and host stars (e.g., mass and metallicity). A correlation between giant planet frequency and stellar metallicity for FGK stars was among the first empirical findings using radial velocity (RV) data (Gonzalez 1997; Fischer & Valenti 2005). Tabachnik & Tremaine (2002) and Cumming et al. (2008) analyzed data sets of RV-detected planets and completeness to derive occurrence rates for giant planets, fitting power-law distributions to the population in mass and semimajor axis. Johnson et al. (2007) and Johnson et al. (2010a) observed a correlation between RV giant planet occurrence rate and stellar host mass, with more massive stars more likely to host giant planets with $a < 2.5$ au, though there is uncertainty in deriving a main sequence mass for subgiants (Lloyd 2011, 2013; Johnson & Wright 2013). Results from asteroseismology suggest spectroscopically derived stellar masses are too high (Johnson et al. 2014; Campante et al. 2017), but not significantly so (North et al. 2017; Stello et al. 2017). Thus the underlying correlation appears robust: for RV-detected giant planets the occurrence rate is proportional to $M_*^\gamma$, with $\gamma = 1.05_{-0.24}^{+0.28}$ (Ghezzi et al. 2018).

Direct imaging in the infrared targets young, nearby stars to detect self-luminous giant planets at high contrast at separations of $\lesssim 1''$ from their parent stars. This technique is sensitive to giant planets at wider separations ($\gtrsim 5$ au), and has thus far been detecting planets primarily around higher-mass stars ($\sim 1.5–1.8 M_\odot$, e.g., Lagrange et al. 2010; Marois et al. 2010), despite the lower intrinsic luminosity of lower-mass stars making it easier to image planetary-mass companions. While higher-mass stars hosting most directly imaged planets is suggestive that the RV trend of rising giant planet fraction with increasing stellar mass holds at wider separations, it has yet to be robustly demonstrated (Wahhaj et al. 2013; Bowler 2016; Galicher et al. 2016). Lannier et al. (2016) found that while imaged intermediate mass-ratio companions (e.g., brown dwarfs) appeared equally common around higher-mass stars (spectral types A and F) and lower-mass stars (M dwarfs), for low-mass-ratio companions (such as giant planets) there was a 74.5%

probability that the distributions around the higher-mass and lower-mass stars were different. With new observations and a larger statistical sample, we reexamine this question.

Brown dwarfs, intermediate in mass between stars and giant planets, have typically been detected in greater numbers by imaging surveys compared to giant planets, despite relatively low brown dwarf occurrence rates (Metchev & Hillenbrand 2009). Recent advances in high contrast imaging instrumentation (Beuzit et al. 2008; Macintosh et al. 2014), observing techniques (Marois et al. 2006), and data reduction (Lafrenière et al. 2007; Soummer et al. 2012; Wang et al. 2015; Pueyo 2016) have led to detections of planets orbiting HR 8799 (Marois et al. 2008, 2010), $\beta$ Pictoris (Lagrange et al. 2009), HD 95086 (Rameau et al. 2013a, 2013b), 51 Eridani (Macintosh et al. 2015), and HIP 65426 (Chauvin et al. 2018), among others. Brandt et al. (2014) proposed that such wide-separation (10–100 au) giant planets and brown dwarfs are part of the same underlying population. This work tests this hypothesis with new observations.

# 2. Survey Description

The Gemini Planet Imager (GPI) is an instrument optimized to directly detect planets at high contrast within $\sim 1''$ radius from their parent star. Coupling a woofer/tweeter extreme adaptive optics system (Poyneer et al. 2014; Bailey et al. 2016) with an apodized Lyot coronagraph (Soummer et al. 2007; Savransky et al. 2014) and an integral field spectrograph (Chilcote et al. 2012; Larkin et al. 2014), GPI routinely reaches contrasts between $\sim 10^{-5}$ and $10^{-6}$ at $0.4''$ from bright stars ($m_I \lesssim 8.0$ mag) after postprocessing (Macintosh et al. 2014; Ruffio et al. 2017a). The coronagraphic inner working angle (IWA) is $2.8\lambda/D$, or $0.12''$ in H band, though in practice the contrast near the IWA after postprocessing depends on additional factors such as the total sky rotation achieved during the observations of each star.

The Gemini Planet Imager Exoplanet Survey (GPIES) utilizes GPI for a 600-star survey, begun in late 2014 and set to conclude in 2019 (Macintosh et al. 2018). Based at Gemini South at Cerro Pachon, Chile, GPIES is targeting the youngest, closest stars in the sky in a systematic survey for giant planets between 10 and 100 au. In particular, the goal of the planet-finding campaign is to directly measure the planet occurrence rate, determine how planet occurrence correlates with stellar properties, probe the atmospheres of giant planets and brown
dwarfs, and determine the orbital properties of substellar companions. The debris disk portion of GPIES is described in depth in T. M. Esposito et al. (2019, in preparation).

Here we present the results from the first 300 stars observed by GPIES. As the first half of the planned survey, 300 stars represents the natural midcourse point to consider the occurrence rate of giant planets partway through the survey. The 300th star was observed on UT 2016-09-21, and in this work we consider only GPIES campaign observations on or before this date. Candidate companions from direct imaging observations can be reduction artifacts, background stars, or bound companions. Second epoch observations confirm that a source is astrophysical and not spurious, and common proper motion indicates whether the star is a bound companion or a background star (e.g., De Rosa et al. 2015; Nielsen et al. 2017). For these 300 stars we have completed second epoch observations of all but 7 of the 55 candidate companions, and most of these remaining 7 have spectra consistent with background stars; follow-up observations are ongoing for the stars observed after this date. Thus the 300-star sample represents an essentially complete sample where candidate companions have been classified into PSF reduction artifacts, background stars, or bound companions.

## 2.1. Target Selection and Properties

A sample of young, nearby stars was constructed by combining members of kinematic associations (de Zeeuw et al. 1999; Zuckerman et al. 2001, 2011) with X-ray selected FGKM stars within 100 pc that also satisfy the observing limits ($m_I \le 9.0$ mag, $\delta \le 25^\circ$). Candidate B and A stars were selected for youth using evolutionary tracks (Siess et al. 2000). We obtained echelle spectra for approximately 2000 candidate stars to further restrict the sample to the ~800 most promising candidates by measuring lithium abundance and chromospheric activity indicators. Apparent binaries with angular separation $0.02-3.0$ and $\Delta\text{mag} \le 5$ (including those discovered during the course of the campaign) were then removed because the presence of the companion degrades closed-loop performance. The target list was updated in late 2016 to replace newly resolved binaries with stars from the original sample and newly identified moving group members. The final GPIES sample consisted of 602 nearby stars. Almost half of the final sample are members of kinematic associations; 152 are from nearby moving groups and 104 from the more distant Scorpius-Centaurus OB2 association. Of the remaining 346 stars not in known moving groups or associations, 92 were B- and A-type stars selected based on their position on the color-magnitude diagram, and 254 were solar-type stars selected based on the strength of chromospheric activity indicators (e.g., Ca H&K). A few stars in the sample were selected for other indicators of planetary systems such as debris disk structure. For example, even though the age of Fomalhaut is >400 Myr (Mamajek et al. 2013), the stellocentric offset of the outer belt (Kalas et al. 2005) and eccentric orbit of Fomalhaut b could be the dynamical outcome of a yet-to-be-detected Fomalhaut c in the system (Kalas et al. 2013).

The distributions of various target properties for the first 300 stars observed during the GPIES campaign are shown in Figure 1. Parallax measurements were obtained from the Gaia DR2 catalog for 251 stars (Gaia Collaboration et al. 2018). We used Hipparcos parallaxes from van Leeuwen (2007) for the remaining 49 stars because these were either too bright for Gaia or had worse astrometric precision in the DR2. Spectral types were obtained from SIMBAD, compiled from a number

of literature sources. The properties for each target are given in the observing log in Table 4.

The ages of the stars in the 300-star sample were estimated using a number of methods: 145 based on membership of a nearby kinematic moving group or association (Zuckerman et al. 2001, 2011), 109 based on activity indicators in their optical spectra or measured X-ray luminosity, and 48 based on their position on the color-magnitude diagram. For these 109 stars, the process to derive their ages from X-ray and spectroscopic indicators, as well as the new high-resolution spectra obtained for a subset, will be discussed in more detail in an upcoming paper. The ages of the kinematic groups are derived from isochrone fitting (e.g., Naylor & Jeffries 2006), lithium depletion boundary analysis (e.g., Binks & Jeffries 2014), trace-back analysis (e.g., Riedel et al. 2017) and eclipsing or resolved spectroscopic binaries (e.g., Nielsen et al. 2016). Applying these various techniques to many coeval stars simultaneously provides tight constraints on the ages of these groups, and as such these ages are preferred.

For the stars not in a known kinematic group, we instead estimated their ages by the measured decrease in lithium abundance and chromospheric activity indicators for lower-mass stars, and the evolution across the color-magnitude diagram for higher-mass stars. We used a combination of lithium abundance, Ca II H/K and $H\alpha$ emission line strength, and X-ray luminosity to assign an age to the lower-mass stars (F-type and later) based on a comparison to similar measures of stars within open clusters and kinematic associations (Mamajek & Hillenbrand 2008; Soderblom 2010). These methods cannot be applied to the more massive stars within the sample (B through early F-type). For example, lithium depletion only occurs efficiently in stars that have convective envelopes deep enough to transport lithium from the photosphere to regions where the temperature is sufficient for lithium burning to occur. As the depth of the convective envelope is inversely proportional to the mass of the star, this process only occurs efficiently in lower-mass stars (Castro et al. 2016). Early-type stars also do not possess the magnetic field necessary to generate the chromospheric activity indicators observed in later-type stars. Instead, we rely on the rapid evolution of these more massive stars across the color-magnitude diagram to determine their ages. The position of early-type (B and A-type) stars on the color-magnitude diagram is primarily dictated by their age and mass. These stars rapidly evolve onto the zero-age main sequence (ZAMS) where they spend roughly one-third of their main sequence lifetime at their ZAMS location on the HR diagram, before the rising helium abundance in their cores causes them to become cooler and more luminous. Secondary effects include metallicity, rotation, binarity, and extinction, which may need to be considered when attempting to determine the age of a star from its position on the color-magnitude diagram. Table 4 presents our best estimate for the ages and masses of each target star, the median of the distribution if posteriors are available. We are currently exploring the production of posteriors for the 109 stars that rely on X-ray, calcium, and lithium age indicators, and will present the results of this analysis and age and mass posteriors for the full sample in an upcoming paper.

### 2.1.1. Higher-mass Stars: Ages and Masses

We combined the stellar evolutionary model MESA (Paxton et al. 2010), used to predict how the fundamental properties of
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Properties of the first 300 GPIES target stars. While lower-mass stars (<1.5 $M_\odot$) are generally nearby (82% are within 60 pc), higher-mass stars are made up of three groups: nearby moving group stars, nearby field stars, and Sco-Cen stars, with 47% of higher-mass stars beyond 60 pc. For these 300 stars, the distributions have median values of 125 Myr, 45.2 pc, F6, and 1.34 $M_\odot$.</figcaption></figure>

a star evolves with time, with the ATLAS9 model atmospheres (Castelli & Kurucz 2004), used to predict the emergent flux from the stellar photosphere, to estimate both the ages and masses of higher-mass stars and the masses of lower-mass stars based on their position on the color–magnitude diagram. We first generated a grid of evolutionary models spanning a range of masses ($0.6 \le M/M_\odot \le 9.8$), metallicities ($-0.5 \le [M/H] \le 2.0$), and initial rotation rates ($0.0 \le \Omega/\Omega_c \le 0.7$, where $\Omega/\Omega_c$ denotes the rotation rate as a fraction of the critical rotation rate; see Paxton et al. 2013), assuming a solar abundance of $X_\odot = 0.7154$, $Y_\odot = 0.2703$, and $Z_\odot = 0.0142$ (Asplund et al. 2009). These models were generated using the same parameters used to create the MESA Isochrones and Stellar Tracks (MIST) model grid (Dotter 2016; Choi et al. 2016), and our models are consistent with their grid at the rotation rates they computed. Rotation is treated by initializing solid-body rotation at the ZAMS, and we exclude the pre-main sequence (PMS) portion of the evolutionary track given the discontinuity seen at the ZAMS for the rapidly rotating stars. For higher-mass stars, the PMS is very brief (~10 Myr), and it is far more likely that a star consistent with the location of a PMS A-type star on the color–magnitude diagram is actually an older star evolving away from the ZAMS. As with the MIST model grid, stars with $M < 1.2 M_\odot$ were only computed without rotation.

We collected optical photometry ($G_B$, $G$, $G_R$) from the Gaia catalog (Gaia Collaboration et al. 2018), combining the catalog uncertainties with the systematic uncertainties reported in Evans et al. (2018). For the handful of stars too bright for Gaia, we used Tycho2 photometry ($B_T$, $V_T$; Høg et al. 2000) instead. We fit these measurements to synthetic photometry derived from a combination of our evolutionary model and the ATLAS9 model atmospheres. We used the emcee affine-invariant Markov Chain Monte Carlo (MCMC) sampler (Foreman-Mackey et al. 2013) to sample the posterior distributions of star mass $M$, age $t$, metallicity $[M/H]$, initial rotation $\Omega/\Omega_c$, inclination of the rotation axis $i$, and parallax $\pi$. We used the following priors:

1. Uniform for $t$ and $\cos i$

2. A Kroupa (2001) initial mass function for $M$

3. A normal distribution for $[M/H]$ ($-0.05 \pm 0.11$; Nielsen et al. 2013)

4. A mass-dependent Maxwellian distribution for the rotation rate $v$ (Zorec & Royer 2012)

5. A combination of a normal distribution (Gaia or Hipparcos measurement and uncertainty) and an assumed uniform space density with an exponential drop-off in distance (Bailer-Jones 2015) for $\pi$

If the star is part of a known kinematic group, we used a normal distribution described by the age and corresponding uncertainty of that group as the age prior. We also use $v \sin i$ measurements as an additional prior, when available. At each step in the MCMC chain, the rotation rate and radius of the star are combined to give equatorial velocity. This is combined with the inclination angle, also an MCMC parameter, to give $v \sin i$ for the star at that step in the chain. The prior on this derived quantity has a form of a Gaussian centered on the measured value with a sigma of the reported error. For known spectroscopic binaries, each stellar component of the system was included in the fit. Three parameters ($M, \Omega/\Omega_c, i$) were added for each additional star in the system, with a single age, metallicity, and parallax describing the system. The mass ratio measured from the orbit fit to the RVs of both components in the double-lined systems was used as an additional prior. For systems that either eclipsed or contained three stars, the mass measurements for each star were used as priors instead. The masses reported in Table 4 for these spectroscopic systems are the total system mass, rather than the mass of the brightest component.

At each step in the MCMC chain the evolutionary model grid was interpolated to determine the luminosity, equatorial radius, and current rotation rate of the star. A two-dimensional model of the star was constructed using the formalism described in Lovekin et al. (2006) and Espinosa Lara & Rieutord (2011) to account for the latitudinal dependence of the effective temperature and surface gravity, and for the angle between the rotation axis and the observer. We sampled the surface of the star equally in latitude and longitude ($90 \times 180$) to model the star. The emergent flux at each point within each
bandpass was computed from the ATLAS9 model atmosphere, given the local effective temperature and surface gravity, the metallicity of the star, and the viewing angle between the surface normal and the observer. The total flux of the star in each bandpass was calculated by weighting the flux from each point by the projected surface area of each segment as viewed from Earth; segments not visible from Earth have zero projected surface area. The total flux was scaled by the distance to the star and compared to the measured photometry. The MCMC chains were advanced until the median and $1\sigma$ credible region derived from the last half of the chain were no longer evolving.

The ages we derive for these higher-mass stars are consistent with previous works using this method. For example, we derive an age of $750_{-190}^{+170}$ for Fomalhaut, consistent with $520 \pm 100$ Myr from Nielsen et al. (2013) and $520_{-60}^{+100}$ from David & Hillenbrand (2015). Mamajek (2012) finds an age for the system of $440 \pm 40$ Myr, based largely on an analysis of the age of the companion K4V star TW PsA (Fomalhaut B). This is similar to the two previous analyses of Fomalhaut itself, and about ~1.5σ lower than our estimate. This age discrepancy is largely a result of different metallicity values used in these analyses, with the Siess et al. (2000) models used by Nielsen et al. (2013) assuming a solar metallicity of $Z_{\odot} = 0.02$, David & Hillenbrand (2015) using the PARSEC model grid which assumes $Z_{\odot} = 0.0152$, while Mamajek (2012) used a $Z_{\odot} = 0.017$ and the Bertelli et al. (2008) model grid. We use a lower value of $Z_{\odot} = 0.0142$ along with the model grid discussed previously. By increasing the metallicity from 0.0142 to 0.017 to match Mamajek (2012), we recover a more consistent age of $550 \pm 70$ Myr. Similarly, increasing the metallicity further to $Z_{\odot} = 0.02$ gives a still lower age, of $440 \pm 70$ Myr. As we note previously, the Asplund et al. (2009) value of $Z_{\odot} = 0.0142$ more closely matches modern estimates of the protosolar abundance.

For stars in young moving groups, this method tends to overestimate the ages of these higher-mass stars, due to the fact that there is essentially no movement across the HR diagram during the first third of a star's main sequence life. For example, the absolute V magnitude of a 1.8 $M_{\odot}$ star will change less than 0.05 mags from the ZAMS (~15 Myr) until ~240 Myr, given the models of Siess et al. (2000). Thus $\beta$ Pictoris, despite being the namesake of the 26 $\pm 3$ Myr $\beta$ Pic moving group, has a poorly constrained age based on its color-magnitude diagram position alone, with a $1\sigma$ upper limit of <345 Myr, and <510 Myr at $2\sigma$. This is similar to the age derived by Nielsen et al. (2013) of $110_{-30}^{+20}$ Myr, and $520_{-170}^{+400}$ Myr from David & Hillenbrand (2015).

Thus, while there is variation between the different literature age determinations using this method, the precise set of models used and especially assumed protosolar metallicity affects the final shape of the age posterior. Nevertheless, the ages used for the higher-mass stars in our sample not in moving groups are derived in a consistent manner, rather than relying on a heterogeneous compilation of literature estimates.

### 2.1.2. Lower-mass Stars: Masses

The technique described in Section 2.1.1 was used for stars bluer than $B-V = 0.35$ (earlier than F1). We used a similar framework to estimate the masses of the lower-mass stars within the sample. Rotation of these stars can safely be ignored; observed $v \sin i$ values are significantly lower than for early-

type stars (Gallet & Bouvier 2013), and as such there is a negligible dependence of the effective temperature and surface gravity on stellar latitude. For these stars, the *emcee* sampler was used to sample the posterior distributions of star mass *M*, age *t*, metallicity [M/H], and parallax $\pi$. As for the higher-mass stars, an additional mass term was included in the fit for each stellar companion within known spectroscopic binary systems. A prior on age was applied for these stars, either for the group age for members of moving groups or associations, or the age distribution derived by activity indicators. At each step the effective temperature, surface gravity, and radius of the star was calculated by interpolating the evolutionary model. The flux from the star in each bandpass was then computed from an interpolation of the ATLAS9 model atmosphere grid, scaled by the square of the ratio of the radius of the star and the distance. These lower-mass stars not in moving groups or associations are assigned a Gaussian age uncertainty with $\sigma$ of 20%. We are in the process of deriving more realistic age uncertainties for individual FGKM stars based on calcium emission, X-ray emission, and lithium absorption, and will present updated age and mass posteriors in a future paper. For this work, then, we assume a single value of age and mass for each star when deriving completeness. For a future paper detailing the final description of the GPIES sample, however, we will incorporate full posteriors in age and mass.

Extinction is not treated for either the higher-mass or lower-mass stars due to the relative proximity of the stars within the Local Bubble. The overall sample is also relatively free of binaries by construction, although there are undoubtedly spectroscopic binaries that remain undiscovered in the sample. The effect of a binary companion would be to make a star appear more luminous and redder, and therefore typically appear older within our analysis. An overestimate of the age of the star would lead to an underestimate of the sensitivity of the observations of that star to planetary-mass companions, and an overestimate in the overall planet frequency in the final analysis. This bias would be most significant for unresolved binaries with a near equal flux ratio, the type of systems that are most amenable for detection via RV. The effect of this bias on the analysis presented below is expected to be small, given that the B and A-type stars with ages determined using this technique only constitute 16% of the sample, and only a small fraction of them are expected to be undiscovered equal-mass spectroscopic binaries.

There are a total of 12 spectroscopic binaries remaining in the sample, and their properties and estimates of the component masses from the spectral energy distribution (SED) fit described previously are given in Table 1. There are several physical binaries and multiples wider than $3''$ (the limit used in constructing the sample)—for example, GJ 3305 orbits 51 Eridani at $66''$7 (Feigelson et al. 2006). We expect these very wide systems to have less of a dynamical influence on the disk and planet formation. A more detailed analysis of the effect binarity has on planet formation requires a larger, dedicated sample that is complete to both single stars and binaries, with early work on this front being carried out by the SPOTS survey (Asensio-Torres et al. 2018). In this work, we consider very close SBs and wide-orbit binaries the same as single stars, though for SBs we adopt the combined mass of the system as the stellar host mass.
<table><caption>**Table 1**
Known Spectroscopic Binaries</caption><thead><tr><th>Name</th><th>Type</th><th>P<br>(days)</th><th>q</th><th>M<sub>1</sub><br>(M<sub>☉</sub>)</th><th>M<sub>2</sub><br>(M<sub>☉</sub>)</th><th>References</th></tr></thead><tbody><tr><td>CC Eri</td><td>SB2</td><td>1.56</td><td>0.54 ± 0.02</td><td>0.54<sup>+0.02</sup><sub>-0.03</sub></td><td>0.29 ± 0.02</td><td>Amado et al. (2000)</td></tr><tr><td>HR 784</td><td>SB1</td><td>37.09</td><td>...</td><td>1.21<sup>+0.02</sup><sub>-0.04</sub></td><td>0.60<sup>+0.11</sup><sub>-0.13</sub></td><td>Gorynya & Tokovinin (2018)</td></tr><tr><td>ζ Hor</td><td>SB2</td><td>12.93</td><td>0.88 ± 0.03</td><td>1.82<sup>+0.12</sup><sub>-0.13</sub></td><td>1.58<sup>+0.04</sup><sub>-0.30</sub></td><td>Sahade & Hernández (1964)</td></tr><tr><td>HR 3395</td><td>SB1</td><td>14.30</td><td>...</td><td>1.24<sup>+0.05</sup><sub>-0.05</sub></td><td>0.84<sup>+0.07</sup><sub>-0.10</sub></td><td>Abt & Willmarth (2006)</td></tr><tr><td>4 Sex</td><td>SB2</td><td>3.05</td><td>0.945 ± 0.002</td><td>1.42<sup>+0.01</sup><sub>-0.02</sub></td><td>1.35<sup>+0.01</sup><sub>-0.01</sub></td><td>Torres et al. (2003)</td></tr><tr><td>HD 142315</td><td>SB1</td><td>1.264</td><td>...</td><td>2.15<sup>+0.08</sup><sub>-0.06</sub></td><td>2.01<sup>+0.06</sup><sub>-0.08</sub></td><td>Levato et al. (1987)</td></tr><tr><td>ζ TrA</td><td>SB1</td><td>12.98</td><td>...</td><td>1.13 ± 0.01</td><td>0.40<sup>+0.05</sup><sub>-0.03</sub></td><td>Skuljan et al. (2004)</td></tr><tr><td>V824 Ara</td><td>SB2</td><td>1.68</td><td>0.909 ± 0.014</td><td>1.18 ± 0.01</td><td>1.08 ± 0.01</td><td>Strassmeier & Rice (2000)</td></tr><tr><td>V4200 Sgr</td><td>SB2</td><td>46.82</td><td>0.983 ± 0.001</td><td>0.88 ± 0.01</td><td>0.87 ± 0.01</td><td>Fekel et al. (2017)</td></tr><tr><td>Peacock</td><td>SB1</td><td>11.75</td><td>...</td><td>5.31<sup>+0.29</sup><sub>-0.27</sub></td><td>4.45<sup>+0.36</sup><sub>-0.63</sub></td><td>Luyten (1936)</td></tr><tr><td>ι Del</td><td>SB1</td><td>11.04</td><td>...</td><td>1.93<sup>+0.10</sup><sub>-0.24</sub></td><td>0.87<sup>+0.76</sup><sub>-0.21</sub></td><td>Harper (1935)</td></tr><tr><td>42 Cap</td><td>SB2</td><td>13.17</td><td>0.727 ± 0.005</td><td>1.94 ± 0.03</td><td>1.41 ± 0.02</td><td>Fekel (1997)</td></tr></tbody></table>

## 2.2. Data Acquisition

GPIES has been operating in priority visitor mode, where members of the GPIES team operate GPI to take campaign observations, since the start of the campaign in late 2014. An automated target-picker was used to select the best star to observe, given the ages and distances of the remaining unobserved stars, and the current hour angle and environmental conditions (McBride et al. 2011). Targets with known substellar companions were not prioritized. The instrumental setup was similar for each target. A typical sequence consists of 38 1-minute coronagraphic exposures taken with GPI's integral field spectrograph at a low spectral resolution ($\lambda/\Delta\lambda = 50$) in the H band (1.5–1.8 $\mu$m). The Cassegrain rotator was fixed for all GPI observations, causing astrophysical signals to rotate in the IFS images as the target being observed transits overhead. Observations were conducted under program codes GS-2014B-Q-500, GS-2014B-Q-501, GS-2015A-Q-500, GS-2015A-Q-501, GS-2015B-Q-500, and GS-2015B-Q-501. Observations in 2016 were conducted under the 2015B semester program codes. An observing log is given for the first 300 stars observed in the campaign—and any subsequent observations taken to confirm or reject candidate companions—in Table 4.

## 3. Data Reduction and Contrast Curves

GPIES data are processed automatically, with the reduction steps described in detail in Wang et al. (2018b); we briefly describe them here. Each frame of integral field spectroscopy data from GPI is processed with the GPI Data Reduction Pipeline (DRP; Perrin et al. 2014, 2016). The GPI DRP produces dark frames and wavelength solutions for each lenslet using calibration data taken during the daytime. Immediately before each observing sequence, a snapshot argon arc lamp is taken and is processed by the GPI DRP to measure the instrument flexure between the observing sequence and the daytime calibration data. For each science frame, the GPI DRP subtracts the dark frame, corrects for bad pixels, compensates the wavelength solution for flexure, constructs 3D spectral data cubes, corrects for optical distortion, and measures the locations and fluxes of the satellite spots that are used for calibration.

A typical observing sequence consists of 38 one-minute exposures, each of which is reduced into a spectral datacube with 37 wavelength channels, resulting in a total of 1406

individual spectral slices, with some correlation between adjacent wavelength channels. The stellar PSF (i.e., speckles) is subtracted from each of these spectral slices. With GPI, we are able to utilize both angular differential imaging (ADI; Marois et al. 2006) and spectral differential imaging (SDI; Racine et al. 1999; Marois et al. 2000; Sparks & Ford 2002) to disentangle the potential planets from the speckles. The speckle subtraction and the planet detection are performed using an open-source Python package `pyKLIP` (Wang et al. 2015), which includes an implementation of the Karhunen–Loève (K–L) Image Projection algorithm (KLIP; Soummer et al. 2012; Pueyo et al. 2015). We follow the forward model matched filter (FMMF) approach described in Ruffio et al. (2017a) to which a few improvements, which are described in the following, were made.

The first step of the reduction consists in removing the time-averaged uncorrected atmospheric speckles, which result in a smooth stellar halo steeply rising near the edge of the coronagraphic mask. When the wind is strong, the halo becomes lopsided and aligned with the direction of the wind, sometimes denoted wind-butterfly. Consequently, the halo cannot be efficiently subtracted with a classical spatial high-pass filter and a Gaussian kernel, which was the approached taken in Ruffio et al. (2017a). In this work, we first subtract the radial intensity profile in the image, which removes the symmetric component of the stellar halo. Then, the lopsided component is subtracted from its projection onto a set of principal components computed from a library of GPI images containing a strong wind-butterfly.

The second processing step consists in subtracting the speckles individually in small sectors of each image. The results presented in Ruffio et al. (2017a) suffered from artifacts related to the edge of the sectors. In this work, we used tight overlapping sectors centered at each separation in the image therefore mitigating the edge-effects.

Finally, flux, standard deviation, and signal-to-noise ratio (S/N) maps are computed according to Ruffio et al. (2017a) using a matched filter accounting for the heteroskedasticity (i.e., non-uniformity) of the noise and the bias in the planet signal induced by the speckle subtraction (Pueyo 2016). The matched filter assumes a specific spectrum for the planet; however, it was demonstrated that two templates (a cool T-type and a warm L-type spectral template) were sufficient to recover
<table><caption>**Table 2**

Detected Companion Properties</caption><thead><tr><th>Name</th><th>Epoch</th><th>ρ<sup>a</sup><br>(arcsec)</th><th>a<sub>proj</sub><sup>a</sup><br>(au)</th><th>ΔH<br>(mag)</th><th>Mass<sup>b</sup><br>(M<sub>Jup</sub>)</th><th>Astrometry References</th></tr></thead><tbody><tr><td>51 Erid b</td><td>20141218</td><td>0.451 ± 0.002</td><td>13.43 ± 0.09</td><td>14.4 ± 0.2</td><td>2.6 ± 0.3</td><td>De Rosa et al. (2015)</td></tr><tr><td>β Pic b</td><td>20151106</td><td>0.249 ± 0.001</td><td>4.84 ± 0.02</td><td>9.3 ± 0.1</td><td>12.9 ± 0.5</td><td>Wang et al. (2016)</td></tr><tr><td>HD 95086 b</td><td>20160229</td><td>0.621 ± 0.005</td><td>53.65 ± 0.46</td><td>13.7 ± 0.2</td><td>2.6 ± 0.4</td><td>Rameau et al. (2016)</td></tr><tr><td>HR 8799c</td><td>20160919</td><td>0.944 ± 0.001</td><td>38.99 ± 0.15</td><td>12.0 ± 0.1</td><td>8.3 ± 0.6</td><td>Wang et al. (2018a)</td></tr><tr><td>HR 8799 d</td><td>20160919</td><td>0.674 ± 0.001</td><td>27.85 ± 0.11</td><td>12.0 ± 0.1</td><td>8.3 ± 0.6</td><td>Wang et al. (2018a)</td></tr><tr><td>HR 8799 e</td><td>20160919</td><td>0.385 ± 0.002</td><td>15.89 ± 0.09</td><td>11.6 ± 0.1</td><td>9.2 ± 0.6</td><td>Wang et al. (2018a)</td></tr><tr><td>HD 984 B</td><td>20150830</td><td>0.216 ± 0.001</td><td>9.90 ± 0.04</td><td>6.5 ± 0.1</td><td>48 ± 3</td><td>Johnson-Groh et al. (2017)</td></tr><tr><td>HR 2562 B</td><td>20160125</td><td>0.619 ± 0.003</td><td>21.07 ± 0.11</td><td>11.7 ± 0.1</td><td>26<sup>+9</sup><sub>-13</sub></td><td>Konopacky et al. (2016)</td></tr><tr><td>PZ Tel B</td><td>20150730</td><td>0.502 ± 0.001</td><td>23.65 ± 0.08</td><td>5.3 ± 0.1</td><td>64 ± 5</td><td>this work</td></tr></tbody></table>

**Notes.**

<sup>a</sup> Projected separation at the first detection in the GPIES campaign.

<sup>b</sup> Mass derived from absolute *H* magnitude.

most spectral types (See Ruffio et al. 2017b). The entire survey was therefore reduced twice, once for each template.

The speckle subtraction has several free parameters that were optimized in Ruffio et al. (2017b) for GPI. We used 30 K-L modes to model the stellar PSF. For each subtraction zone, the K-L modes were generated from the 150 most correlated reference images selected from the subset of the observing sequence in which a planet in that zone would move by at least 0.7 pixel in *H* band due to a combination of ADI and SDI

Planet detection thresholds are defined from the standard deviation map after calibration (in unit of planet-to-star flux ratio). The standard deviation map calibration is described in Ruffio et al. (2017a). It guarantees that the S/N map has a standard deviation of one and corrects for any algorithm related flux losses (i.e., algorithm throughput) using extensive planet injection and recovery. Then, the 1D standard deviation as a function of separation is computed as the root-mean-square of the standard deviation map in 4 pixel concentric annuli. Detection thresholds are generally derived from a false positive rate chosen a priori. However, the estimation of planet occurrence rates, the goal of this work, requires that all the candidates above the detection threshold be vetted. The vetting process generally requires follow-up observations, which is a limited resource due to telescope time constraints. As a result, we define the detection threshold of σ as the boundary used over the GPIES campaign for candidate follow-up. This 1D detection threshold is the contrast curve. Unfortunately, while the stellar halo subtraction generally improved GPIES' sensitivity at intermediate separations, it introduced a large number of false positives at the very edge of the coronagraphic mask when using the T-type spectral template. We therefore increased the detection threshold to 8σ below <0.′3 in these reductions. Full contrast curves for the entire GPIES campaign, including the subset discussed here, will be presented along with the final statistical analysis of the campaign in a future paper.

## 4. Companion Detections

Nine stellar companions around seven host stars have been detected within the first half of the GPIES survey: three brown dwarfs and six planetary-mass companions. The measured position, contrast, and derived mass for each companion is given in Table 2. The first GPIES campaign image and corresponding sensitivity curves are given in

Figure 2 for the imaged planets, and Figure 3 for the brown dwarfs. 51 Eridani b and HR 2562 were new discoveries from these images, while β Pictoris b, HD 95086 b, HR 8799 cde, HD 984 B, and PZ Tel B were detections in our campaign images of known substellar companions.

### 4.1. 51 Eridani

51 Eridani b was discovered early in the GPIES campaign (Macintosh et al. 2015) around the F0IV β Pictoris moving group member (Bell et al. 2015) 51 Eridani (hereafter 51 Eri). The *H*-band spectrum measured in the discovery epoch exhibited strong methane and water absorption, consistent with both predictions for low-mass exoplanets around young stars, and also the observed spectra of low-temperature T-type brown dwarfs. Follow-up astrometric observations confirmed that 51 Eri b was in a bound orbit around 51 Eri (De Rosa et al. 2015), with a negligible probability of it being an interloping field T-dwarf in chance alignment with 51 Eri.

Subsequent spectroscopic observations with GPI (Rajan et al. 2017) and Very Large Telescope (VLT)/SPHERE (Samland et al. 2017), combined with thermal-IR photometry obtained with Keck/NIRC2 (Macintosh+15, Rajan+17), have been used to characterize the atmospheric properties of the planet. The observed spectral energy distribution is consistent with a photospheric cloud fraction between 75% and 90% (Rajan et al. 2017) and 100% (Samland et al. 2017), depending on the specifics of the cloud model used in the analysis. Both fits find a similar effective temperature (T<sub>eff</sub> = 600–760 K), but a range of values for the surface gravity (log g = 3.5–4.5 (dex)). The estimated mass of the planet is strongly dependent on the assumptions made regarding its formation. Under an assumption of a high-entropy (so-called “hot-start”) formation scenario (Marley et al. 2007), the luminosity of the planet corresponds to a mass of 1–2 M<sub>Jup</sub> (Rajan et al. 2017). Uniquely for young directly imaged companions, the luminosity of 51 Eri b is also consistent with predictions of low-entropy, “cold-start” formation scenarios (Marley et al. 2007; Fortney et al. 2008), with a significantly higher model-dependent mass of 2–12 M<sub>Jup</sub>. Measurements of the reflex motion induced by the orbiting planet, either via spectroscopy or astrometry, are required to differentiate between these two scenarios.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Top row: FMMF S/N maps of the four planetary systems resolved within GPIES after image reduction and PSF subtraction. The 51 Eri S/N map was calculated using a T-type template, the others with an L-type template. Bottom row: corresponding completeness maps for the four stars computed using the procedure described in Section 5.1. Contours are plotted for the 90% (solid), 50% (dashed), and 10% (dotted) completeness levels. The location of each companion is based on the mass derived from the H-band photometry and the projected separation of the companion at the epoch of the campaign observation.</figcaption></figure>

## 4.2. β Pictoris

The exoplanet β Pictoris b (β Pic b) was discovered around the A6V namesake of the β Pictoris moving group using VLT/NaCo in 2003, later confirmed to be a bound companion in follow-up observations taken with the same instrument in 2008 (Lagrange et al. 2009, 2010). The planet is orbiting interior to a large circumstellar disk (Smith & Terrile 1984) and has a measurable effect on the morphology of the inner disk (Lagrange et al. 2012; Millar-Blanchaer et al. 2015). β Pic b was observed several times with GPI—as a first light target during commissioning, as a part of an astrometric monitoring program to constrain the orbit of the system, and finally as part of the GPIES Campaign. Due to its proximity and youth, the primary was a highly ranked GPIES target, and was selected by our automated target picker early in the campaign, resulting in the detection of β Pic b. Milli-arcsecond astrometry obtained with GPI was critical for constraining the orbital parameters of the planet (Wang et al. 2016). Observed variability in the light curve of the host star β Pic was interpreted as being caused by a transit of either the planet or by material within its Hill sphere (Lecavelier Des Etangs & Vidal-Madjar 2009). Wang et al. (2016) used GPI astrometry to tightly constrain the inclination of the orbit of the planet and excluded a transit at the $10\sigma$ level, although a transit of the Hill sphere of the planet was still a possibility. Precise timings on the transit of the planet's Hill sphere were made, allowing for a unique opportunity to probe the circumplanetary environment of a young exoplanet.

The low contrast between β Pic b and its host star compared to other directly imaged exoplanets makes it an ideal target for atmospheric characterization. Chilcote et al. (2017) used a combination of GPI spectroscopy and literature photometry (compiled and calibrated by Morzinski et al. 2015) to measure

the spectral energy distribution of the planet between 0.9 and 5 µm, wavelengths spanning the bulk of the emergent flux of the planet. The data were consistent with models that incorporated the condensation of dust into thick vertically extended clouds within the photosphere, and matched the observed spectral energy distributions of isolated low-surface gravity brown dwarfs (Chilcote et al. 2017). Unlike 51 Eri b, the bolometric luminosity of β Pic b is only consistent with the predictions of the “hot-start” luminosity models, and inconsistent with the Fortney et al. (2008) “cold-start” models. However, this is not a criticism of core accretion; see, for example, Berardo et al. (2017), Owen & Menou (2016) for how core accretion can lead to hot/warm starts, and also our Section 6.3. Chilcote et al. (2017) derived a model-dependent mass of $12.7 \pm 0.3 M_{Jup}$, consistent with the model-independent mass derived from *Hipparchos* and *Gaia* astrometry of the host star of $11 \pm 2 M_{Jup}$ (Snellen & Brown 2018) and $13 \pm 3 M_{Jup}$ (Dupuy et al. 2019).

## 4.3. HD 95086

HD 95086 b was discovered around the A8III member of the 15 Myr (Pecaut & Mamajek 2016a) Lower Centaurus Crux subgroup of the Scorpius-Centaurus OB2 association (Rizzuto et al. 2011) using VLT/NaCo in 2011 (Rameau et al. 2013a, 2013b), one of a growing number of exoplanets that have been resolved around Sco-Cen members (Bailey et al. 2014; Chauvin et al. 2017). The planet was noted for its unusually red infrared colors that were ascribed to a photosphere dominated by thick clouds (Meshkat et al. 2013). HD 95086 was observed as a part of the GPIES campaign to characterize the atmosphere of the planet, monitor its orbital motion, and to search for additional interior companions. Like
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Top row: PyKLIP S/N maps of the three systems with a brown dwarf companion resolved from the first 300 GPIES stars after image reduction and PSF subtraction. Bottom row: corresponding completeness maps for the three stars computed using the procedure described in Section 5.1. Contours are plotted for the 90% (solid), 50% (dashed), and 10% (dotted) completeness levels. The location of each companion is based on the mass derived from the H-band photometry and the projected separation of the companion at the epoch of the campaign observation.</figcaption></figure>

β Pic, the target star was selected for observations for our planet-search program by our automated target picker, resulting in a detection of the planet. The low luminosity and very red infrared colors of the planet combined with its distance of 86 pc make it a challenging target for spectroscopic observations, especially at shorter wavelengths. De Rosa et al. (2016) reported the first spectroscopic measurement of the planet obtained using GPI, a spectrum covering the blue half of the K band. This was combined with literature photometry at *J* and *L'* to reveal a spectral energy distribution consistent with model atmospheres incorporating significant amounts of photospheric dust, a result confirmed in a later analysis incorporating SPHERE observations of the planet (Chauvin et al. 2018).

Monitoring the orbital motion of the planet can provide insight into the architecture of the HD 95086 system. The presence of two large circumstellar dust rings was inferred from the spectral energy distribution of the star (Moór et al. 2013; Su et al. 2017); the outer ring was subsequently resolved with millimeter interferometric observations with ALMA (Su et al. 2017). HD 95086 b lies between these two rings and, despite only sampling a small fraction of the full orbit of the planet with GPI and literature astrometry, Rameau et al. (2016) were able to constrain the eccentricity of the orbit to be near-circular, under the assumption of co-planarity with the outer ring. While the dynamical impact of the planet on the outer ring is still uncertain (Su et al. 2017), Rameau et al. (2016) demonstrated

that it is unlikely that the planet is responsible for clearing the gap in to the predicted outer radius of the inner disk, and posited this as indirect evidence of the existence of additional lower-mass companions that are below present detection limits, a finding also confirmed by Chauvin et al. (2018).

### 4.4. HR 8799

The HR 8799 system is currently the only example of a multiple planet system detected via direct imaging (Marois et al. 2008, 2010; Konopacky & Barman 2018). The host star is a λ Bootis star that exhibits γ Doradus variability (Gray & Kaye 1999; Zerbi et al. 1999). It is a probable member of the 42<sup>+6</sup><sub>-4</sub> Myr (Bell et al. 2015) Columba moving group (Zuckerman et al. 2011; Malo et al. 2013) and has a spectral type of either A5 and F0, depending on which set of spectral lines are used for typing (Gray et al. 2003). The star is also host to two circumstellar dust rings, at ~10 and ~100–300 au, with the planets occupying the gap between them, and a large diffuse dust halo at ~300–1000 au (Su et al. 2009). While the inner ring is only inferred from the spectral energy distribution of the star, the outer ring has been resolved with ALMA and the SMA (Booth et al. 2016; Wilner et al. 2018).

The HR 8799 system has been observed multiple times with GPI: as an early science target during the commissioning of the instrument (Ingraham et al. 2014), and during the campaign to both characterize the atmospheres of the known planets
(Greenbaum et al. 2018) at multiple bands, and as part of the planet search in *H* when chosen by the automated target-picker. The campaign observations resulted in the detection of three of the four planets, HR 8799 cde, with HR 8799 b outside the field of view of GPI in our nominal search mode. Astrometry obtained with GPI was also used in conjunction with measurements obtained with Keck/NIRCC2 to investigate the dynamical properties and long-term stability of the system (Wang et al. 2018a). By utilizing *N*-body simulations, the combination of astrometry and dynamical limits placed tighter constraints on the orbital parameters and the masses of the individual planets, which when combined with evolutionary models give a mass for HR 8799 b of $5.8 \pm 0.5 M_{\text{Jup}}$, and $7.2_{-0.7}^{+0.6} M_{\text{Jup}}$ for planets c, d, and e.

### 4.5. Brown Dwarf Companions

In addition to the four planetary systems discovered during the first half of the GPIES campaign, three brown dwarf companions were resolved: PZ Tel B, HD 984 B, and HR 2562 B. The companions to PZ Tel (Biller et al. 2010; Mugrauer et al. 2010) and HD 984 (Meshkat et al. 2015) were known before the GPIES observations of these stars, and were detected when their host stars were selected by our automated target-picker, while HR 2562 B was first discovered with GPI (Konopacky et al. 2016) as part of our planet search campaign. PZ Tel and HD 984 are probable members of the *β Pictoris* and Columba moving groups, respectively (Malo et al. 2013), while HR 2562 is thought to be <1 Gyr based on photometric (Casagrande et al. 2011) and spectroscopic (Pace 2013) analyses. The three host stars have a later spectral type than the planet host stars discussed previously; G9IV (1.1 $M_{\odot}$) for PZ Tel, F7V (1.3 $M_{\odot}$) for HD 984, and F5V (1.3 $M_{\odot}$) for HR 2562.

A brown dwarf companion was later detected around the campaign target HD 206893 (initially discovered with the VLT/SPHERE by Milli et al. 2016), but these observations were obtained after the 300th star was observed and the companion is therefore not included in any of the analyses presented below for the first half of the GPIES campaign.

## 5. Survey Modeling and Inference

Having observed 300 stars and detected six planets and three brown dwarfs, GPIES is one of the largest, deepest direct imaging surveys, with a high yield for substellar companions. Using this rich data set, we proceed to measure the overall occurrence rate of giant planets and trends in the underlying population.

### 5.1. Survey Completeness to Date

We determine the completeness to planets for each of our observations using the Monte Carlo procedure described in Nielsen et al. (2008, 2013), Nielsen & Close (2010). For each target star, completeness to substellar companions is determined over a grid that is uniform in log space in mass and semimajor axis. At each grid point, $10^4$ simulated substellar companions, each with the same value of mass and semimajor axis, are randomly assigned orbital parameters: inclination angle from a $\sin(i)$ distribution, eccentricity from a linear distribution fit to wider-separation RV planets (Nielsen et al. 2008), $P(e) = 2.1 - 2.2e$ with $0 \le e \le 0.95$, and argument of periastron and epoch of periastron passage from a uniform

distribution. Since contrast curves are one-dimensional, position angle of nodes is not simulated. The projected separation (from the orbital parameters and mass and distance of the host star) and the $\Delta H$ (from the mass of the companion, the age and absolute magnitude of the host star, and luminosity models) for each companion are then computed. A simulated companion that lies above the contrast curve is considered detectable, while one that lies below (or outside the field of view) is undetectable. If more than one contrast curve is available for a single star, the simulated companions generated at the first epoch are advanced forward in their orbits to the next epoch, and compared to the next contrast curve; a companion is detectable if it is above at least one contrast curve.

As noted previously, for the FMMF contrast curves described in Section 2.2, we use an $8\sigma$ threshold inside of $0.1''3$, and a $6\sigma$ threshold outside. FMMF contrasts were produced for both a T-type and L-type template for each star (Ruffio et al. 2017a). The T-type reduction assumes significant methane absorption in the *H* band, so that images in redder channels may be subtracted from images in bluer channels without subtracting the planet from itself. To stay consistent, the same evolutionary models that produce *H* band flux are used to give temperature as a function of mass and age, and planets hotter than 1100 K are assumed undetectable using the T-type curve. The L-type curve, which does not suffer from the possibility of self-subtraction, is used for simulated companions of all temperatures. These FMMF contrast curves are defined out to $1.1''7$, though beyond $1.1''1$ some parts of the field go beyond the edge of the detector, and observations are not complete over all position angles. We account for this with the same method as Nielsen et al. (2013), by recording the fractional completeness as a function of separation, and using a uniform random variable to reject simulated planets that would not fall on the detector.

We utilize the CIFIST2011 BT-Settl atmosphere models (Caffau et al. 2011; Allard 2014; Baraffe et al. 2015)$^{45}$, tied to the COND “hot-start” evolutionary model grid (Baraffe et al. 2003), to give absolute *H* magnitude and temperature as a function of mass and age. From these, we produce completeness maps for each star. These maps give the fractional completeness to companions as a function of mass and semimajor axis (completeness at a given value of mass and semimajor axis is not binary since orbital parameters will set whether a given companion is detectable for a certain contrast curve). These completeness maps are shown for stars with detected companions in Figures 2 and 3. Summing the maps across all target stars produces a map representing our depth of search (e.g., Lunine et al. 2008) for the survey, as given in Figure 4. The contours give the number of stars to which the survey is complete to companions of a given mass and semimajor axis. Since target stars cover a large range of distances (3.2–176.7 pc), the limited field of view of GPIES (~0.1"15–1.1"1, beyond which not all position angles lie on the detector) means that completeness plots for individual stars do not precisely line up. Individual completeness plots never reach 100% for a similar reason, as some simulated planets at each semimajor axis can fall within the IWA or outside the outer working angle, and so go undetectable. As such the maximum value reached on the plot is 246 stars, though all 300 stars have some fractional sensitivity to substellar companions. The

<footnote>$^{45}$ https://phoenix.ens-lyon.fr/Grids/BT-Settl/CIFIST2011/COLORS/colmag.BT-Settl.server.2MASS.Vega</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Depth of search for the first 300 stars observed by GPIES, showing the number of stars to which the survey is complete for planets and brown dwarfs as a function of mass and semimajor axis. Overplotted in red circles are the detected companions, plotted at the projected separation they were first imaged by GPIES.</figcaption></figure>

detected companions are indicated by red dots, plotted at their inferred mass and projected separation at the first GPIES epoch, as given in Table 2. In this figure, and in the subsequent analysis, we assume masses for each companion inferred from the BT-Settl models, the age of the system, and the *H*-band magnitude, even if a more robust mass is available from full SED fitting, to maintain consistency in the completeness calculation. Most of these assumed masses are close to the published mass, with the notable exception of HD 95086, where the *H*-band magnitude and age suggest a mass of $2.6 \pm 0.4 M_{\odot}$, compared to $4.4 \pm 0.8 M_{\odot}$ inferred from the *Ks*-band magnitude and a bolometric correction (De Rosa et al. 2016). HD 95086 is a particularly unusual case, as its dusty atmosphere leads to a redder spectrum than most models predict. We note that the analysis below does not change significantly whether a mass of 2.6 or $4.4 M_{\odot}$ is used.

We also investigated the effect of utilizing full age posteriors for each star, rather than using a single age for each target. We considered the AB Dor moving group target star AN Sex, and computed the completeness plot for a single age of 149 Myr (as used in this work), for a disjoint Gaussian, and for a symmetric Gaussian. The disjoint Gaussian has a $\sigma$ of 19 Myr below 149 Myr, and a $\sigma$ of 51 Myr above, to match the Bell et al. (2015) age for the AB Dor moving group of $149_{-19}^{+51}$ Myr. The symmetric Gaussian, computed for comparison, is given a $\sigma$ of 20 Myr (13% age uncertainty). As expected, assuming younger ages result in more detectable planets compared to older ages; however, the effect on the completeness plot is marginal. The 20% completeness contour, for example, moves at most 7% in planet mass, with a median shift of 1%, for the disjoint Gaussian. For the symmetric Gaussian, the effect is even smaller, maximum of 5% and median 0.2%. In the symmetric case, the contours move toward smaller planets (becoming

more sensitive when using a symmetric Gaussian age posterior compared to a single age), while the completeness plot is less sensitive for the disjoint Gaussian. Similar results are found for the symmetric Gaussian even if the value of $\sigma$ is doubled. We thus expect a minor effect on our results by including age posteriors for each target star, but this effect will be explored more in depth in a future paper on the statistical constraints from the full survey, once full age posteriors are available for all of our stars.

Utilizing mass posteriors is expected to have an even smaller effect on the completeness plots. Since simulated planets are generated as a function of semimajor axis, stellar mass has no effect on the completeness for a star observed at only a single epoch (273 out of our 300 stars), as planets are randomly assigned a mean anomaly, which does not require the period to be known. For stars observed at multiple epochs, however, period is needed to determine the amount of orbital motion between the two observations, but given the slow orbital speeds of these wide-separation planets, and that the observations described here span less than two years, this is a minor effect as well, especially as we typically reach $\lesssim 5\%$ precision on stellar mass on our targets.

## 5.2. A Correlation between Stellar Mass and Planet Occurrence Rate

A shared feature of the planetary companions described in Section 4 is that all orbit the higher-mass stars in our sample: all host stars lie between 1.55 and $1.75 M_{\odot}$ (Figure 5). As shown in Figure 1, there is no clear mass bias in our sample toward higher masses; in fact, 177 out of 300 stars (59%) have masses below $1.5 M_{\odot}$. This trend is also the opposite from what we would expect for observational biases: lower-mass stars are intrinsically fainter, and thus the same achievable contrast would correspond to a lower-mass detectable planet.

We investigate the strength of this trend by comparing occurrence rates of wide-separation giant planets around higher-mass stars and lower-mass stars, here defined to have a mass determination above and below $1.5 M_{\odot}$. We consider planets with semimajor axis between 3 and 100 au, and planet mass between 2 and $13 M_{\odot}$. These limits are chosen to encapsulate a region of planet parameter space to which our observations are most sensitive and encompass the detected planets in the GPIES sample, which are between 2.6 and $12.9 M_{\odot}$ and a projected separation from 4.84 to 53.65 au (as shown in Figure 4). While the depth of search is computed for semimajor axis, the most direct measurement we have of each detected planet is projected separation, so we have chosen a broad enough range of semimajor axes so that planets at these projected separations are unlikely to fall outside this range. Indeed, orbital monitoring of each of these planets (De Rosa et al. 2015; Wang et al. 2016, 2018a; Rameau et al. 2016) place them solidly in this range. We also begin by assuming planets are uniformly distributed over this range in log space in both semimajor axis and planet mass ($\propto a^{-1}m^{-1}$). This is not too dissimilar to the power-law distributions for small-separation ($a < 3.1$ au) giant planets found by Cumming et al. (2008) of $\propto a^{-0.61}m^{-1.31}$, obtained by converting the period distribution of $\propto P^{-0.74}$ to semimajor axis for solar-type stars. For each star, then, we integrate the completeness to planets over this range with uniform weight given to each log bin in semimajor axis and mass.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Depth of search for the sample divided by stellar mass. All three brown dwarfs appear around lower-mass stars, which is not especially surprising given those stars comprise almost 60% of the sample. The fact that all six planets (and so all four planet-hosting stars) are among higher-mass stars, however, is more striking, and suggests a correlation between stellar mass and occurrence rate of intermediate period giant planets.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Left: posterior probability of frequency of giant planet systems (3–100 au, 2–13 $M_{Jup}$) for higher-mass stars (blue) and lower-mass stars (red). Planets with these parameters are more common around higher-mass stars in this sample. Filled circles represent the median of each distribution, and error bars give the 1 and 2σ confidence intervals. Right: posterior for the difference between the two frequencies, with the lower-mass star occurrence rate subtracted from the higher-mass star occurrence rate. The histogram (blue when the higher-mass star occurrence rate is larger, red otherwise) shows the results from Monte Carlo draws from the two posteriors, and the black dashed line gives the posterior computed using Equation (1). At 99.92% confidence, the occurrence rate is larger for higher-mass stars (>1.5 $M_{\odot}$) than lower-mass stars (<1.5 $M_{\odot}$).</figcaption></figure>

Integrating over the entire range from 3 to 100 au and 2 to 13 $M_{Jup}$ produces a completeness to massive giant planets of 17.5 stars among the higher-mass (>1.5 $M_{\odot}$) sample, and 43.3 among the lower-mass sample. These two samples have 4 and 0 planetary systems detected, respectively, as HR 8799 counts as a single planetary system. We infer the frequency of planetary systems for each subsample using a Bayesian approach with a Poisson likelihood ($L = \frac{e^{-\lambda}\lambda^k}{k!}$). The expected number of planetary systems ($\lambda$) is the number of stars to which the subsample is complete multiplied by the frequency, and the measured number ($k$) is the number of detected planetary

systems. The prior is taken to be the Jeffrey's prior on the rate parameter of a Poisson distribution, $\text{Prior}(\lambda) \propto \sqrt{\frac{1}{\lambda}}$. We compute the posterior over a well-sampled regular grid of frequency for both subsamples, plotting the results in Figure 6.

Given no planets were detected around lower-mass stars (<1.5 $M_{\odot}$), our measurement of the frequency for these stars represents an upper limit, <6.9% at 95% confidence. With four detected systems, the fraction of higher-mass stars with giant planet systems is $24^{+13}_{-10}\%$ (68% confidence interval). To evaluate the confidence with which we can conclude that the higher-mass (HM) star frequency is larger than the lower-
mass (LM) star frequency, we compute the posterior of the difference, $\delta = f_{\text{HM}} - f_{\text{LM}}$, given by

$$P(\delta|\text{data}) \propto \int_{0}^{\infty} P_{\text{HM}}(f|\text{data})P_{\text{LM}}(f - \delta|\text{data})df, \quad (1)$$

where $P_{\text{HM}}$ and $P_{\text{LM}}$ are the posteriors on the frequency ($f$) for higher- and lower-mass stars. We confirm this expression with a Monte Carlo method, where a set of two frequencies are randomly generated from the two posteriors and differenced. The two methods generated identical results, as shown in the right panel of Figure 6. For 99.92% of the samples, wide-separation giant planets around higher-mass stars are more common than around lower-mass stars, a >3σ result.

We note that this result is not strongly dependent on our choice of limits in planet mass, semimajor axis, and stellar mass. When increasing the upper semimajor axis limit from 100 to 300 au, the significance of the result that wide-separation giant planets are more common around higher-mass stars drops marginally from 99.92% to 99.85%. Since we are averaging over the depth of search plot in Figure 4, expanding the range to areas of lower sensitivity reduces the effective number of stars to which the survey is sensitive, and for lower-mass stars this number drops from 43.3 to 36.2, and higher-mass stars drops from 17.5 to 17.4, when expanding the semimajor axis range from 100 to 300 au. Increasing this range still further from 3–100 to 1–1000 au reduces the effective number of lower-mass stars to 24.9 and higher-mass stars to 12.0, and the significance in this case remains above 3σ at 99.84%.

Similarly, dropping the planet mass range considered from 2–13 $M_{\text{Jup}}$ to 1–13 $M_{\text{Jup}}$ has a negligible effect, with the significance remaining at 99.92%. Dropping the 12 SBs from the sample is also a minor effect, removing mainly higher-mass stars, and raising the significance slightly to 99.95%. Stellar mass has a more significant effect, since moving the boundary between higher-mass and lower-mass stars results in fewer lower-mass stars in the sample. The significance drops to 99.70% when the boundary is moved to 1.35 $M_{\odot}$ (an effective number of lower-mass stars of 37.9), and 96.99% at 1.10 $M_{\odot}$ (24.3 stars). Thus the larger frequency of wide-separation giant planets around higher-mass stars is relatively robust to specific choices of which range of parameter space we choose.

### 5.3. Wide-separation Giant Planet Occurrence Rate around Higher-mass Stars

Following the example of Cumming et al. (2008), we consider a model of planet distributions defined by power laws in mass and semimajor axis. We adopt a functional form of our model of

$$\frac{d^2N}{dm da} = f C_1 m^\alpha a^\beta, \quad (2)$$

where $m$ and $a$ are planet mass and semimajor axis. We define this equation over a limited range, $[m_1 \le m \le m_2]$ and $[a_1 \le a \le a_2]$. The frequency ($f$) is the occurrence rate of planets in this range, and so the normalization constant $C_1$ is then given by

$$C_1 = \left[ \int_{m_1}^{m_2} \int_{a_1}^{a_2} m^\alpha a^\beta dm da \right]^{-1}. \quad (3)$$

Thus, if $\alpha \ne -1$ and $\beta \ne -1$,

$$C_1 = \frac{\alpha + 1}{m_2^{(\alpha+1)} - m_1^{(\alpha+1)}} \frac{\beta + 1}{a_2^{(\beta+1)} - a_1^{(\beta+1)}}, \quad (4)$$

otherwise, the two terms become $[\ln(a_2) - \ln(a_1)]^{-1}$ and $[\ln(m_2) - \ln(m_1)]^{-1}$.

In order to constrain the three free parameters of this model ($f$, $\alpha$, and $\beta$), we adopt a Bayesian approach utilizing Monte Carlo completeness calculations. Bayes' equation then becomes

$$P(f, \alpha, \beta|\text{data}) \propto P(\text{data}|f, \alpha, \beta) P(f) P(\alpha) P(\beta), \quad (5)$$

and in standard Bayesian terminology, the first term is the posterior, the second the likelihood, and the final three terms are the priors. We follow a similar method to Kraus et al. (2008), dividing the two-dimensional observational space (mass versus projected separation) into a series of bins. Similar Bayesian methods were also used by Biller et al. (2013), Wahhaj et al. (2013), and Brandt et al. (2014) to fit a power-law model to imaged planets. We use 41 bins in mass, logarithmically spaced between 1 and 100 $M_{\text{Jup}}$, and 81 bins in separation, logarithmically spaced between 1 and 1000 au. In each bin, we calculate the expected and actual number of planets in that bin; as a result, we choose a Poisson distribution for the likelihood in each bin, then find the product of probabilities across all the bins,

$$P(\text{data}|f, \alpha, \beta) = \prod_i \prod_j \frac{e^{-E_{i,j}} E_{i,j}^{O_{i,j}}}{O_{i,j}!}, \quad (6)$$

where $i$ and $j$ are indices for bins in mass and projected separation, and $E_{i,j}$ and $O_{i,j}$ are the expected and observed number of planets in each bin, respectively.

Unlike Kraus et al. (2008), our model and likelihood use different parameters: Equation (2) is in terms of mass and semimajor axis, while Equation (6) is in terms of mass and projected separation, since directly imaged planets are detected with a particular separation, but an (often unknown) semimajor axis. Using Monte Carlo simulations, however, we are able to project our modeled planets from semimajor axis into separation. As in Section 5.1, we inject simulated planets at a grid of mass and semimajor axis, and for each grid point, save the fraction of simulated planets that are above the contrast curve (detectable with our observations). Additionally, we also sort the detectable planets into projected separation bins, tracking the projected separation of each detectable planet. The completeness map for each star, then, is a three-dimensional array across mass, projected separation, and semimajor axis. Marginalizing over the separation dimension, and summing over all stars, gives the depth of search shown in Figure 4. If we instead marginalize over semimajor axis, we derive a completeness map in mass and projected separation, the same parameters that define our observed planets.

For a given set of ($f$, $\alpha$, $\beta$), we follow this procedure for evaluating Equation (5):

1. For each star, generate a three-dimensional completeness map using Monte Carlo simulations.

2. Sum the completeness maps over all stars.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Posteriors on our double power-law model of giant planet populations between 5–13 $M_{\text{Jup}}$ and 10–100 au, for stars $>1.5 M_{\odot}$. We find a relatively high occurrence rate of $9^{+5}_{-4}\%$ over these ranges, and negative power laws that point to a rise in planet frequency at smaller masses and smaller semimajor axis (SMA). Red, blue, and green contours give 1, 2, and $3\sigma$ contours on probability.</figcaption></figure>

3. Weight each point in the map using Equation (2) and the values of ($f, \alpha, \beta$), giving the expected number of planets detected in each bin.

4. Marginalize the map over semimajor axis, producing a map of expected planets versus mass and separation.

5. Produce a corresponding observed map of detected planets versus mass and separation.

6. Compute the likelihood from Equation (6).

7. Multiply by priors.

We then utilize a Metropolis-Hastings MCMC procedure (as in Nielsen et al. 2014) to explore this parameter space and generate posteriors on frequency and power-law indices. For priors, we utilize uniform priors for $\alpha$ and $\beta$, and a Jeffrey's prior for $f$ appropriate to a Poisson distribution, $\text{Prior}(f) \propto \sqrt{\frac{1}{f}}$. Our results are not overly dependent on this choice of prior, we find similar posteriors with a prior that is uniform in $\ln(f)$. We first restrict ourselves to a region of parameter space where we have both good completeness and our detections, $[1 M_{\text{Jup}} \le m \le 13 M_{\text{Jup}}]$ and $[1 \text{ au} \le a \le 100 \text{ au}]$, and $M_{*} > 1.5 M_{\odot}$. This range encompasses all six planets detected by the GPIES survey. Given our low sensitivity to low-mass, short-period planets, there is a strong degeneracy between frequency

and power-law indices, as very negative indices (that place large numbers of planets at small separations and small masses) result in larger frequencies. Fits with large occurrence rates place most of these planets in the region of parameter space to which GPIES is least sensitive. In order to report an occurrence rate with high confidence, we define occurrence over a region of high completeness, and so we quote not the frequency over the entire range, but instead a modified frequency over the truncated range $[5 M_{\text{Jup}} \le m \le 13 M_{\text{Jup}}]$ and $[10 \text{ au} \le a \le 100 \text{ au}]$. Thus while 51 Eridani b, HD 95086 b, and $\beta$ Pic are excluded from this modified frequency, they fall within the broader range, and thus are all used to constrain the model. The resulting posteriors are shown in Figure 7 and Table 3.

We find a relatively high giant planet occurrence rate: $9^{+5}_{-4}\%$, for stars with mass $>1.5 M_{\odot}$, and planets in the range $5 M_{\text{Jup}} \le m \le 13 M_{\text{Jup}}$ and $10 \text{ au} \le a \le 100 \text{ au}$. Given the very negative values for $\alpha$ and $\beta$, if we were to extrapolate beyond the region of high completeness to lower planet masses and smaller semimajor axes, we would infer a much larger occurrence rate with significantly larger uncertainties. Vigan et al. (2012) combined the IDPS sample of 42 stars of AF spectral types with three literature stars and, assuming values of $\alpha = \beta = 0$, found a fraction of high-mass stars with at least one $3–14 M_{\text{Jup}}$
<table><caption>**Table 3**
Constraints on Substellar Object Populations from the GPIES 300-star Sample, between 10 and 100 au</caption><thead><tr><th></th><th>f</th><th>α</th><th>β</th><th>γ</th></tr></thead><tbody><tr><td colspan="5">Planets around high-mass stars: 5 M<sub>Jup</sub> &lt; m &lt; 13 M<sub>Jup</sub>, 1.5 M<sub>☉</sub> &lt; M<sub>*</sub> &lt; 5 M<sub>☉</sub></td></tr><tr><td>Median</td><td>8.9%</td><td>-2.37</td><td>-1.99</td><td>...</td></tr><tr><td>68% CI</td><td>5.3%-13.9%</td><td>-3.16 to -1.51</td><td>-2.47 to -1.44</td><td>...</td></tr><tr><td>95% CI</td><td>2.9%-20.6%</td><td>-3.90 to -0.56</td><td>-2.92 to -0.81</td><td>...</td></tr><tr><td colspan="5">Planets around all stars: 5 M<sub>Jup</sub> &lt; m &lt; 13 M<sub>Jup</sub>, 0.2 M<sub>☉</sub> &lt; M<sub>*</sub> &lt; 5 M<sub>☉</sub></td></tr><tr><td>Median</td><td>3.5%</td><td>-2.277</td><td>-1.68</td><td>2.03</td></tr><tr><td>68% CI</td><td>2.1%-5.4%</td><td>-2.98 to -1.48</td><td>-2.18 to -1.10</td><td>1.05-3.07</td></tr><tr><td>95% CI</td><td>1.1%-8.0%</td><td>-3.65 to -0.54</td><td>-2.63 to -0.38</td><td>0.09-4.16</td></tr><tr><td colspan="5">Brown dwarfs around all stars: 13 M<sub>Jup</sub> &lt; m &lt; 80 M<sub>Jup</sub>, 0.2 M<sub>☉</sub> &lt; M<sub>*</sub> &lt; 5 M<sub>☉</sub></td></tr><tr><td>Median</td><td>0.8%</td><td>-0.47</td><td>-0.65</td><td>-0.85</td></tr><tr><td>68% CI</td><td>0.3%-1.6%</td><td>-1.53 to -0.67</td><td>-1.69 to -0.89</td><td>-2.33-0.59</td></tr><tr><td>95% CI</td><td>0.1%-2.8%</td><td>-2.60 to -1.96</td><td>-2.52 to -3.86</td><td>-3.76-2.01</td></tr><tr><td colspan="5">Substellar objects around all stars: 5 M<sub>Jup</sub> &lt; m &lt; 80 M<sub>Jup</sub>, 0.2 M<sub>☉</sub> &lt; M<sub>*</sub> &lt; 5 M<sub>☉</sub></td></tr><tr><td>Median</td><td>4.7%</td><td>-2.09</td><td>-1.44</td><td>1.15</td></tr><tr><td>68% CI</td><td>3.2%-6.7%</td><td>-2.49 to -1.71</td><td>-1.92 to -0.91</td><td>0.35-1.96</td></tr><tr><td>95% CI</td><td>2.0%-9.1%</td><td>-2.89 to -1.33</td><td>-2.36 to -0.28</td><td>-0.48-2.75</td></tr></tbody></table>

**Note.** The power-law model has a number of substellar companions per star of *f* within the given mass and semimajor axis limits, with power-law indexes of *α*, *β*, and *γ* for companion mass, semimajor axis, and stellar host mass, respectively.

planet between 5 and 320 au of $8.7_{-2.8}^{+10.1}\%$. Nielsen et al. (2013) examined the 70 B and A stars observed by the Gemini NICI Planet-Finding Campaign, and with a null result for planets placed $2\sigma$ upper limits on the fraction of high-mass stars with planets. These limits have <10% of high-mass stars having a 10–13 $M_{Jup}$ planet with semimajor axis between 38 and 650 au, and <20% with a 4–13 $M_{Jup}$ planet between 59 and 460 au. Rameau et al. (2013) considered 37 AF stars, including detections of planets around HR 8799 and $\beta$ Pic, finding a planet fraction of $7.4_{-2.4}^{+3.6}\%$ from 5 to 320 au, between 3 and 14 $M_{Jup}$, assuming $\alpha = \beta = 0$. Bowler (2016) considered multiple direct imaging surveys, performing a meta-analysis on 384 unique target stars. This sample was then broken up by spectral type, into subsamples of M, FGK, and BA. For the BA sample (roughly equivalent to our $>1.5 M_{\odot}$ sample), by assuming $\alpha = \beta = -1$, Bowler (2016) find a planet fraction of 5–13 $M_{Jup}$ 30–300 au planets of $2.8_{-2.3}^{+3.7}\%$.

We compare our results to large imaging surveys with detected planets that explicitly measure occurrence rate in the top left panel of Figure 8, where we plot each survey as a line showing the semimajor axis distribution. End points are taken to be the range of semimajor axes considered by each survey. We utilize the value of $\alpha$ from each survey to find the expected occurrence rate over the GPIES mass range (5–13 $M_{Jup}$). In most cases, authors assumed a value of $\beta$ (and $\alpha$) of -1 or 0, as we did in Section 5.2, rather than fitting for the power-law index directly. Vigan et al. (2012) and Rameau et al. (2013) find similar results, as expected given both had similar numbers of stars and detected planets. The Bowler (2016) meta-analysis included these surveys and others, for a total of 110 BA stars, and finds a lower occurrence rate, and assumes a distribution that is flat in log space, rather than one flat in linear space. In all cases, the

GPIES result is similar in fraction of planets at ~100 au, but rises more steeply toward lower semimajor axis. This is as we expect, given our finding of a more negative value of $\beta$ than -1 or 0, placing more planets close to the star.

We also plot possible extrapolations of RV occurrence rates for giant planets. Johnson et al. (2010a) reported a value of 20% planet fraction for retired A stars, for semimajor axis less than 2.5 au, and $K < 20$ m s$^{-1}$ (which they note corresponds to a lower-mass limit of $1.12 M_{Jup}$ for a stellar host of $1.6 M_{\odot}$). Johnson et al. (2010a) do not fit for $\alpha$, and do not give the upper mass for planets in their sample. We assume both are similar to the Cumming et al. (2008) sample for planets around FGK stars, and use the latter values for upper mass limit of $10 M_{Jup}$ and $\alpha = -1.31$, in order to convert the occurrence rate for all planets to 5–13 $M_{Jup}$. Similarly, while Johnson et al. (2010a) give an upper limit of 2.5 au, they do not state the minimum semimajor axis of planets in their sample, though we speculate that this minimum value is 0.08 au, for HD 102956 b (Johnson et al. 2010b). As there is no estimate for the value of $\beta$ for RV planets around higher-mass stars, we investigate three different values, flat in log space, the Cumming et al. (2008) value, and flat in linear space ($\beta$ of -1, -0.61, and 0, respectively), and plot these as gray lines in Figure 8.

For integrated occurrence rate, the Johnson et al. (2010b) value of 20% (<2.5 au) represents a rise from our value of 11.4% (10–100 au), though for different mass ranges. A more detailed study of the distributions of RV-detected giant planets around higher-mass stars is required to definitively determine the nature of the distribution of high-mass planets within 10 au, and whether RV and imaged giant planets around higher-mass stars follow a continuous distribution, but with a break in the planet mass and semimajor axis power-law indexes at ~3 au.
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Comparison of the occurrence rates of 5–13 $M_{Jup}$ planets (top two panels), 13–80 $M_{Jup}$ brown dwarfs (bottom left), and 5–80 $M_{Jup}$ substellar companions (bottom right) as a function of semimajor axis for the GPIES results and other surveys. The top left panel considers only higher-mass stars from GPIES, while all 300 stars are used in the other three panels with a stellar mass power-law term (in these three panels, the occurrence rate is given for $1 M_{\odot}$). The end points of each line represent the range of semimajor axis reported (10–100 au for GPIES), with the slope giving the value of $\beta$, the semimajor axis power-law index, and the normalization showing the planet occurrence rate. We utilize the value of the companion mass power law ($\alpha$) in each survey to find the occurrence rate over the GPIES range in each panel, and in the last three panels we use our value of $\gamma$ (the power-law index for stellar host mass) to convert the GPIES occurrence rate from 1.75 $M_{\odot}$ to 1.0 $M_{\odot}$, to more closely match analyses. Solid lines denote fits to the value of $\beta$ (GPIES; Cumming et al. 2008; Brandt et al. 2014; Fernandes et al. 2019), while dashed lines are used for authors that assumed power-law indexes ($\alpha$ and $\beta$ of -1 or 0) rather than fitting for them directly. For clarity, we omit error bars and only display the median values for power-law indexes and occurrence rates.</figcaption></figure>

## 5.4. Occurrence Rate with a Stellar Mass Dependence

We can instead fit the entire GPIES 300-star sample by including a term that allows the occurrence rate to vary as a function of stellar mass. We thus introduce an additional term into Equation (2)

$$ \frac{d^2N}{dm da} = f C_1 m^{\alpha} a^{\beta} \left( \frac{M_*}{1.75 M_{\odot}} \right)^{\gamma} \quad (7) $$

so that the overall occurrence rate $f$ now varies with stellar mass. We choose $1.75 M_{\odot}$ as the normalization location, given this is where most of our planet-hosting stars lie, and so occurrence rate is most constrained at this mass. Similarly, Equation (5) becomes

$$ P(f, \alpha, \beta, \gamma|\text{data}) \propto P(\text{data}|f, \alpha, \beta, \gamma) \times P(f) P(\alpha) P(\beta) P(\gamma). \quad (8) $$

We solve this equation as before, but with one additional change: instead of summing the completeness map over all stars, we instead introduce an extra dimension in stellar mass, and bin the sample by stellar mass. This produces a four-dimensional completeness map, as a function of companion mass, companion separation, companion semimajor axis, and stellar host mass. Similarly, the map of expected planets is now over three dimensions, companion mass, companion separation, and stellar host mass.

The results from the MCMC fit are plotted in Figure 9. As expected from the prevalence of detected planets around higher-mass stars, the distribution for $\gamma$ is peaked at relatively high values, $\gamma = 2.0 \pm 1.0$. Interestingly, though likely coincidentally, this value of $\gamma$ is close to the negative of the Salpeter IMF power-law index of -2.35 ($\frac{dN_*}{dM_*} \propto M_*^{-2.35}$; Salpeter 1955), which would mean that for a volume-limited sample of young stars, while the number of giant planets per
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Fit to the entire GPIES sample with a triple power-law model of wide-separation giant planet populations between 5–13 $M_{\text{Jup}}$ and 10–100 au, with the occurrence rate ($f$) defined for 1.75 $M_{\oplus}$ stars. Given that all planets were detected around higher-mass stars, a large power-law index for stellar mass dependence is preferred.</figcaption></figure>

star would increase with increasing stellar mass, the number of giant planets per stellar mass bin would only change slightly, decreasing as the cube root of the stellar mass ($\frac{dN_p}{dM_*} \propto M_*^{-0.35}$). When we compare this fit to our earlier fit to only higher-mass stars (Figure 7), the planet mass and semimajor axis power-law indexes, $\alpha$ and $\beta$, are similar for both fits. Occurrence rate, $f$, is significantly lower for the fit utilizing a stellar mass dependence. While a more positive value of $\gamma$ (~4) better explains the additional 177 lower-mass stars without planets in our sample, such a distribution begins to overpredict the number of detectable planets around the highest mass stars in our sample ($\gtrsim 2 M_{\oplus}$). Thus the best fit is found by combining a moderate value of $\gamma$ (2.0 ± 1.0) and a lower overall occurrence rate. While this fit incorporates the entire 300-star sample, we find the fit exclusively to the 123 higher-mass stars presented above to be more robust, as it centers around the detected planets. Multiple models could account for the stellar mass dependence of giant planet occurrence seen in Figure 5. Here we have introduced a scaling of overall occurrence rate with stellar mass, but such a result could also be explained by a more negative power-law fit to either companion mass or semimajor axis for lower-mass stars compared to $>1.5 M_{\oplus}$ stars, or some combination of

varying power-law indices and occurrence rate as a function of stellar mass. With no detections around our lower-mass star sample, we are unable to determine which model best represents these wider-separation giant planets.

Galicher et al. (2016) described statistical results from the IDPS survey of 292 stars and two previous surveys, analyzing a sample of 356 stars with spectral types between B and M. Assuming values of $\alpha = \beta = 0$, a planet fraction was inferred of $1.05_{-0.70}^{+2.80}\%$ between 0.5 and $14 M_{\text{Jup}}$ and 20–300 au, based on detections of planetary systems around HR 8799 and HIP 30034 (AB Pic). An analysis of the VLT/NACO large program by Vigan et al. (2017) found similar results to IDPS with similar assumptions and analysis methods. The survey observed 85 stars with NACO, which was then combined with previously published samples to reach a final sample of 199 FGK stars. This final sample contained two detections of brown dwarfs, GSC 08047-00232 B and PZ Tel B, along with AB Pic b. With $\alpha = \beta = 0$, the occurrence rate of $0.5–14 M_{\text{Jup}}$, 20–300 au planets was found to be between 0.85% and 3.65% at $1\sigma$, with a best value of 1.15%. From the Bowler (2016) meta-analysis, assuming $\alpha = \beta = -1$, an occurrence rate of $0.6_{-0.5}^{+0.7}\%$ was derived for planets between 5–13 $M_{\text{Jup}}$ and 30–300 au around stars between B and M spectral types.
The top right panel of Figure 8 compares direct measurements of the planet occurrence rates to our findings for GPIES, where again our steeper slope places more planets at smaller semimajor axis. We have used our median value of $\gamma$ to convert these occurrence rates from 1.75 to $1 M_\odot$, to better match literature values. As before, we use the values of $\alpha$ assumed from each reference to convert the given mass ranges to the GPIES range of 5–13 $M_{Jup}$. The occurrence rates estimates for different imaging surveys give about the same value at $\sim 100$ au. Cumming et al. (2008) and Fernandes et al. (2019) analyze RV surveys and find values of $\beta = -0.61 \pm 0.15$ and $\beta = -0.025_{-0.22}^{+0.30}$, respectively, at small separations. While Cumming et al. (2008) fit a single power law, Fernandes et al. (2019) find evidence for a broken power law, with a break at $\sim 3$ au. At these wider separations the Fernandes et al. (2019) value of $\beta = -1.975_{-0.30}^{+0.22}$ (from their symmetric epos model) is a good match to our value of $\beta = -1.68_{-0.50}^{+0.57}$, though we note that while these RV surveys contain mainly FGK stars, all our planet detections are orbiting more massive stars, and so our GPIES results for $1 M_\odot$ come from assuming a continuous distribution jointly constrained by non-detections on these lower-mass stars and six detections around higher-mass stars in the sample.

A number of analyses have also considered the giant planet occurrence rate at these separations for specific categories of target stars. Biller et al. (2013) considered 80 stars that are members of moving groups that were observed by the Gemini NICI Planet-Finding Campaign, and used a Bayesian technique to fit a variant of Equation (2), where the upper semimajor axis cut-off was a free parameter. This analysis included $\beta$ Pic b (which was not detected in the first NICI campaign observation of the star), and AB Pic b, but the small number of detections meant that limits could not be placed on the values of $\alpha$ and $\beta$, though a maximum in probability was reached for a planet fraction of 4%, from 10 to 150 au. Wahhaj et al. (2013) performed a similar analysis on a combination of the NICI debris disk subsample and the debris disk hosts in the IDPS high-mass sample (Vigan et al. 2012), with detections of $\beta$ Pic b and HR 8799 bcd. This analysis found $2\sigma$ limits that placed $\beta$ approximately between $-2.5$ and $-1$, and $\alpha$ between $-1$ and $2.5$, with a peak value for the fraction of stars with planets of $\sim 10\%$. Wahhaj et al. (2013) also fit a model with a stellar mass dependence to the average planet multiplicity ($F_P$), which is a separate term from the fraction of stars with planets. These results are suggestive of a rise of planet multiplicity with stellar mass, but with a large degree of uncertainty. For their relation $F_P \propto M_*^\gamma$, the $1\sigma$ results place $\gamma$ between 0.3 and 3.2, and a wider $2\sigma$ range between $-0.5$ and $4.7$, which encompasses the value of 0, corresponding to no stellar mass dependence.

Bryan et al. (2016) examined the occurrence rate of wide-separation giant planets in systems with a previously known inner RV planet. By combining long-period RV trends with (generally non-detections from) AO imaging, the possible mass and semimajor axis of the companions could be partially constrained. Their high occurrence rate ($62.1_{-5.7}^{+3.4}\%$ for substellar companions between 1 and $20 M_{Jup}$ between 5 and 100 au) is noticeably higher than what we find with GPIES, suggesting that the presence of a close-in giant planet substantially boosts the likelihood of a wider substellar companion. While they also fit values of $\alpha$ and $\beta$, Bryan et al. (2016) note that these posteriors are strongly dependent on the choice of mass and semimajor axis limits, as those limits truncate the probability distributions for each companion.

A recent paper by Stone et al. (2018) examined the constraints on occurrence rate from a sample of 98 stars observed in the LBT LEECH survey. While a number of planets were observed in the characterization sample, the statistical sample had no planet detections. For all stars, using the COND models Baraffe et al. (2003), they find a $2\sigma$ upper limit of 25% on the occurrence rate of planets between 4–14 $M_{Jup}$ and 5–100 au, consistent with our results.

### 5.4.1. Brown Dwarf Companions

In addition to detecting six planets, the GPIES campaign has also detected three brown dwarfs, and so we apply our Bayesian framework to inferring the underlying population of brown dwarfs. As before, we fit the four-parameter model of Equation (7) to brown dwarfs in the sample, restricting our fit to semimajor axis values between 1 and 100 au, and companion mass values between 13 and $100 M_{Jup}$.

The depth of search, as shown by Figure 4, is much more uniform over this mass range, and represents a much deeper sensitivity to these companions: the three brown dwarfs all lie within or near the 160-star contour. Figure 10 displays the results from the brown dwarf analysis from the GPIES 300-star sample.

There is a tail to the posterior at large values of semimajor axis power law, which is truncated by the introduction of a top-hat prior, uniform between values of $-5$ and $5$, and $0$ elsewhere. This tail results from the limited number of brown dwarfs, and our model being defined in semimajor axis while our detections are made in separation. A very positive power-law index would place most brown dwarfs at the largest possible semimajor axis allowed, 100 au in our model. Even if all brown dwarfs have a semimajor axis of 100 au, there is a non-zero probability of detecting three companions at projected separations between 9.9 and 23.7 au.

As for planets, we again truncate the range for reporting the occurrence rate, 10–100 au and $13–80 M_{Jup}$. As expected given the deeper sensitivity and fewer detected objects, the brown dwarf occurrence rate is noticeably lower than for giant planets, $0.8_{-0.5}^{+0.8}\%$.

This occurrence rate is consistent with previous literature measurements of the brown dwarf companion frequency. Metchev & Hillenbrand (2009) conducted a deep survey of 100 stars (out of a larger 266-star sample), and determined an occurrence rate of brown dwarfs ($13–75 M_{Jup}$) orbiting between 28 and 1590 au of $3.2_{-1.7}^{+3.1}\%$, assuming $\alpha = 0$ and $\beta = -1$. The occurrence rate of a few percent for wide-separation brown dwarfs has been noted by a number of studies, including Oppenheimer et al. (2001), Brandt et al. (2014), Uyama et al. (2017), and Cheetham et al. (2015).

In the bottom left panel of Figure 8, we compare the semimajor axis distributions of GPIES and Metchev & Hillenbrand (2009). Our result is for a lower occurrence rate than that of Metchev & Hillenbrand (2009), with a more negative power-law slope. This decreasing of slope as we move closer to the star is consistent with results from RV surveys for an even smaller brown dwarf occurrence rate at closer separations, with a brown dwarf frequency $<1\%$ within 5 yr periods (Grether & Lineweaver 2006).

Additionally, we find no significant evidence for a stellar mass dependence of the brown dwarf occurrence rate. We find a value for $\gamma$, the stellar mass power-law, of $-0.9 \pm 1.5$. Thus a value of 0, representing no stellar mass dependence, lies
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Posteriors on the same parameterized model used to fit planets from the GPIES survey, but this time for brown dwarfs. A significantly lower occurrence rate is observed, due to the higher sensitivity to these masses coupled with a detection rate lower by a factor of two.</figcaption></figure>

within the $1\sigma$ confidence interval. A similar conclusion was reached by Lannier et al. (2016), who computed companion frequency for stars from the MASSIVE survey of 58 M stars (which detected two substellar objects) and the survey of 59 mostly dusty AF stars of Rameau et al. (2013). Focusing on companions with an intermediate mass ratio of 0.01–0.05 (corresponding to 10.5 to 52 $M_{Jup}$ for a solar-type primary), Lannier et al. (2016) find no statistical difference between the higher-mass and lower-mass sample. Similarly, the Gemini NICI Planet-Finding Campaign surveyed young 70 B and A stars, finding three brown dwarfs (Nielsen et al. 2013), while Bowler et al. (2015) found four brown dwarfs from a survey of 78 single M stars, again suggesting no strong stellar mass dependence in the wide-separation brown dwarf occurrence rate.

### 5.4.2. A Single Population of Substellar Companions at Wide Separations

Finally, we consider the entire range of substellar objects, and attempt to fit planets and brown dwarfs with the same distribution. GPIES has detected a total of 9 substellar objects, spanning a factor of 25 in companion mass, a factor of 11 in projected separation, and a factor of 1.6 in stellar mass. Thus if a single power-law distribution can describe all these objects,

the GPIES 300-star sample provides excellent constraints on the parameters of this distribution.

Figure 11 displays the results from fitting the entire GPIES substellar object yield with Equation (7). Compared to the other subsamples, the power-law indices are most constrained with this fit, as expected from the larger number of objects to fit. It is also not surprising that the medians of the posteriors of the three power-law indices for the fit to all substellar objects are intermediate to posteriors for the fits to only planets and to only brown dwarfs.

A similar analysis was undertaken by Brandt et al. (2014) based on the SEEDS sample and previously published surveys. Over a similar range to ours, 5–70 $M_{Jup}$ and 10–100 au, Brandt et al. (2014) found values of $\alpha = -0.65 \pm 0.60$ and $\beta = -0.85 \pm 0.39$, compared to the values we derive here of $\alpha = -2.1 \pm 0.4$ and $\beta = -1.4 \pm 0.5$. Thus we are finding significantly more negative power laws from our GPIES sample, placing more planets at smaller semimajor axis, and smaller masses. Brandt et al. (2014) also found an overall occurrence rate of substellar companions in this regime of $1.7_{-0.5}^{+1.1}\%$. While we explicitly fit for $\gamma$, Brandt et al. (2014) do not, but instead adopted a fixed value of $\gamma = 1$. The comparison is shown in the bottom right panel of Figure 8, with the two distributions consistent at ~100 au, but the number
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Fit to the four-parameter model using all nine substellar companions detected in the GPIES 300 star sample.</figcaption></figure>

of planets more steeply rising for GPIES when moving to smaller separations. Vigan et al. (2017) also computed the occurrence rates for substellar companions (0.5–75 $M_{Jup}$, 20–300 au) of $2.1^{+1.95}_{-0.60}\%$, assuming $\alpha = \beta = 0$. With the same assumed distribution, Galicher et al. (2016) derived an occurrence rate of 0.35%, with 95% confidence between 0.10% and 1.95%, for substellar companions between 5–70 $M_{Jup}$ and 10–100 au.

The different values found of $\alpha$, $\beta$, and $f$ by GPIES and Brandt et al. (2014) can be accounted for by differences in the sensitivity and detected companions between the two samples. Brandt et al. (2014) found one planetary-mass companion (GJ 504 b), five brown dwarfs, and one intermediate object crossing the planet/brown dwarf boundary ($\kappa$ And b). By contrast, GPIES detected more lower-mass companions, finding six planets and three brown dwarfs. This leads to a more negative companion mass power law, predicting more lower-mass companions. In addition, the better IWA of GPI has led to three detections of substellar objects closer than 20 au (51 Eri b, $\beta$ Pic b, and HD 984 B), which in turn drives the semimajor axis power law more negative as well. Finally, the larger number of detections, as well as power laws placing more planets at smaller masses (with lower sensitivity), are responsible for the larger occurrence rate we infer.

### 5.5. Cold Start

In addition to the “hot-start” evolutionary models of giant planets we have considered up until this point, there is an alternate set of “cold start” models. In the cold start models (Marley et al. 2007; Fortney et al. 2008), as the planet forms most of the gravitational potential energy of the infalling gas is radiated away very early on ($\lesssim 5$ Myr) in a shock, so that young planets have comparatively less internal energy to radiate away compared to their hot-start counterparts, and so are cooler and less luminous at early ages. At older ages ($\gtrsim 300$ Myr), the cold start and hot-start luminosity tracks converge.

The direct observational consequence of the cold start models is fainter, more difficult to detect planets than would be predicted from the hot-start models. Additionally, these models predict a pile-up in luminosity at the youngest ages, where planets between 2 and 10 $M_{Jup}$ maintain luminosities between about 1 and $3 \times 10^{-6}$ L$_\odot$ for the first 30 Myr of their evolution (Fortney et al. 2008).

We use cold start evolutionary models (Fortney et al. 2008), and couple them with the AMES-COND atmospheric models (Allard et al. 2001; Baraffe et al. 2003) by matching values of luminosity and temperature predicted by the cold start evolutionary grid with gridpoints in the AMES-COND models to extract *H* magnitude. We then produce a depth of search plot
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Depth of search for cold start luminosity models. Of the nine substellar objects detected by GPIES, only 51 Eri b has a luminosity low enough to potentially be a cold start planet. The Sonora 2018 grid (right) predicts generally brighter giant planets, with the exception of a limited range of parameter space: between ~10–50 Myr and $\lesssim 3 M_{\text{Jup}}$, the Fortney et al. (2008) models are slightly brighter. As a result, the GPIES results are overall more sensitive to giant planets with the Sonora 2018 grid, but a handful of very young nearby stars result in the extension of the four-star contour to the lowest masses for Fortney et al. (2008).</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Comparison of the two sets of cold start models we consider here, Fortney et al. (2008) and the Sonora 2018 grid (M. S. Marley et al. 2019, in preparation). The Sonora 2018 grid predicts generally higher temperatures and larger luminosities for giant planets at young ages.</figcaption></figure>

from the GPIES survey, given in the left panel of Figure 12. There are currently no cold start models for brown dwarf masses, and of the six planets detected by GPIES, five have luminosities too large to have formed according to the cold start models. Only 51 Eri b could potentially be a cold start planet. However, due to the pile-up in luminosity of this model grid, the estimated mass can only be constrained to be between 4 and $10 M_{\text{Jup}}$ (Macintosh et al. 2015).

As there is only one possible detection of a cold start planet, rather than attempt to fit a power-law distribution, we instead solve for the occurrence rate of cold start giant planets with an assumed power-law distribution. Following the method

described in Section 5.2, we adopt a uniform distribution in log space, $\alpha = \beta = -1$. The average value over the depth of search between 2 and $13 M_{\text{Jup}}$, between 3 and 100 au shown in Figure 12 is then 10.9 stars. From Bayes’ equation with a Poisson likelihood and a Jeffrey’s prior, and one detection over this range, we find an occurrence rate of $15_{-9}^{+15}\%$. This is intermediate to the value we found for hot-start giant planets around high-mass stars when assuming $\alpha = \beta = -1, 24_{-10}^{+14}\%$ and the $2\sigma$ upper limit for low-mass stars of <6.9% (we caution that 51 Eridani contributed to these occurrence rates for both sets of models). It is not yet certain which set of evolutionary models 51 Eri b follows, and so this cold start
occurrence rate could be overestimated if it is indeed a hot-start planet. Additionally, these hot-start and cold start models represent extremes in evolutionary models; there also exist “warm start” models (Spiegel & Burrows 2012; Mordasini et al. 2017) that give mass as a function of the initial entropy of the planet (e.g., Figure 17 in Rajan et al. 2017). Direct measurements of the mass of 51 Eri b from Gaia astrometry of the host star could break the degeneracy between the different models.

We also consider a cold start variation of the new Sonora 2018 grid (M. S. Marley et al. 2019, in preparation) with cloud-free atmospheres. This new set of cold start models predicts that young planets are generally hotter by ~200 K and ~5× more luminous than the Fortney et al. (2008) models, as shown in Figure 13. Key differences between the two model sets include the assumed helium abundance (Y = 0.24 and 0.28 for M. S. Marley et al. 2019, in preparation) and Fortney et al. 2008, respectively), and the mass of the baryon in the equation of state, with Fortney et al. (2008) using 1 amu and Marley et al. (2019) using the mass of a hydrogen atom.

When using the Sonora 2018 grid, the mass of 51 Eri b is better defined, no longer lying in the overlap region of multiple mass tracks, but at a higher mass of $5.2 M_{\text{Jup}}$. As before, no other detected GPIES planets are low luminosity enough to lie in the range of *H* magnitudes predicted for 1–10 $M_{\text{Jup}}$ planets. The depth of search with the Sonora 2018 grid is given in the right panel of Figure 12. The generally brighter planets results in deeper sensitivity to high-mass planets, though still not as deep as for the hot-start BT-Settl models. For the lower-mass giant planets ($\lesssim 3 M_{\text{Jup}}$), however, the Fortney et al. (2008) models predict slightly brighter planets between ~10 and 50 Myr. Due to a small number of nearby stars with ages $\lesssim 50$ Myr, then, the four-star contour of the depth of search reaches 1.5 $M_{\text{Jup}}$ for the Fortney et al. (2008) models, compared to 3 $M_{\text{Jup}}$ for Sonora 2018. Additionally, the Sonora 2018 grid has GPIES able to reach planets at larger semimajor axis, since a high luminosity planet means it can be detected around more distant target stars, so cold start planets can be detected around more of the GPIES sample. Repeating the occurrence rate calculation above gives a slightly lower value, given the greater sensitivity, with $11_{-6}^{+11\%}$ of stars hosting a Sonora 2018 cold start planet between 2–13 $M_{\text{Jup}}$ and 3–100 au.

Compared to hot-start planets, potential cold start planets are very rare among the directly imaged exoplanets. 51 Eri b and GJ 504 b are the only two imaged planets have low enough luminosity to be consistent with cold start models, and so how such a population depends on planet mass, semimajor axis, and stellar mass is not yet well understood. Unlike previous generations of direct imaging exoplanet surveys, GPIES is reaching higher contrasts within 1", however surveys such as NICI (Biller et al. 2013; Liu et al. 2010; Nielsen et al. 2013; Wahhaj et al. 2013) and IDPS (Vigan et al. 2012; Galicher et al. 2016) reached higher contrasts than GPIES at $\gtrsim 2"$, further from the glare of the central star. For young, nearby moving group stars this corresponds to sensitivity to cold start planets to separations of $\gtrsim 50$ au. The relative dearth of detections could indicate a lower occurrence rate than hot-start planets, or that cold start planets are found predominantly at a much smaller semimajor axis. If 51 Eri b is indeed a cold start planet, and the occurrence rate of ~10% we infer here is accurate, then future surveys, with an upgraded GPI (Chilcote

et al. 2018) or an upgraded SPHERE, or an ELT high-contrast imager, could allow us to better study these planets.

## 6. Discussion

At 300 stars, and with detections of six planets and three brown dwarfs, the early GPIES Campaign data is allowing us to probe the distribution of directly imaged substellar companions at lower masses and closer separations than has been previously possible.

### 6.1. Do Wide-separation Brown Dwarfs and Giant Planets Follow the Same Distribution?

The 13 $M_{\text{Jup}}$ boundary that notionally separates giant planets from brown dwarfs represents the onset of deuterium burning, and a long-standing question has been whether there is a significant difference in formation mechanism between these two classes of objects (e.g., Schlaufman 2018). Brown dwarfs and perhaps also giant planets of the highest masses are thought to represent the low-mass outcome of the process by which stars form (i.e., gravitational instability; Forgan & Rice 2013; Kratter & Lodato 2016) or turbulent fragmentation (e.g., Hopkins 2013). Objects of lower mass (Jupiters and below) may have formed from the “bottom up” via core accretion (i.e., by coagulation of solids; e.g., Goldreich et al. 2004; Ormel 2017) and subsequent accretion of gas (e.g., Harris 1978; Mizuno et al. 1978; Mizuno 1980; Pollack et al. 1996). Evidence for this dichotomy is seen for close-in companions observed by RV and transit surveys, where there exists a brown dwarf desert lacking companions near 30 $M_{\text{Jup}}$ (e.g., Grether & Lineweaver 2006; Triaud et al. 2017). Such a deep minimum in the occurrence rate versus companion mass distribution has not been seen at wider separations (Gizis et al. 2001; Metchev & Hillenbrand 2009), though Kratter et al. (2010) and Currie et al. (2011) suggest a possible gap in substellar companion occurrence rate between a mass ratio of ~0.01–0.02, from 10–100 au. Such a lower mass limit is consistent with simulations which find that the lowest mass objects formed by disk fragmentation should be ~10 $M_{\text{Jup}}$ (Thies et al. 2010). Brandt et al. (2014) found that a single population can fit both brown dwarfs and giant planets at wide separations.

While turbulent fragmentation from a molecular cloud—a process integral to the formation of stars—has been invoked to explain isolated planet-mass companions in systems with primaries not too different in mass (e.g., 2MASSJ 1207334-393254 and 2MASS J11193254-1137466; Chauvin et al. 2004; Best et al. 2017), such a mechanism is not naturally applied to other directly imaged planets, which show evidence for having formed in a disk. The orbit of *β* Pic b is nearly coplanar with the star’s debris disk (Lagrange et al. 2010; Macintosh et al. 2014; Nielsen et al. 2014); HD 95086 b appears to exhibit this same coplanarity with its star’s disk (Rameau et al. 2016); and the HR 8799 planets are likely coplanar with each other and with the surrounding debris disk (Wang et al. 2018a). Coplanarity points to these planets forming in a disk, either by core accretion or gravitational instability, rather than by turbulent fragmentation from a cloud.

In Section 5.3 we fitted for the parameters to Equation (7) for three cases: using all substellar companions, only giant planets, and only brown dwarfs. In Figure 14 we overplot the posteriors on the four parameters for brown dwarfs alone and planets alone. Since the posteriors overlap, we cannot completely rule
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Comparison of population parameter fits to planets and brown dwarfs, based on our GPIES sample. At 1–2σ significance, brown dwarfs and giant planets appear to follow different underlying distributions between 10 and 100 au.</figcaption></figure>

out the possibility that the two populations were drawn from the same underlying distribution. Nevertheless, when we follow the same procedure as in Section 5.2, we find the probabilities that the values of α and β are larger for brown dwarfs at 87.3% and 74.3%, while the probabilities that γ and f are smaller for brown dwarfs are 94.4% and 95.1%, respectively. Thus, we find it likely that planets and brown dwarfs at 10–100 au separation follow different underlying distributions—planets obey a more bottom-heavy mass distribution, and are concentrated toward smaller SMA—at a confidence level between 1 and 2σ.

The greatest degree of overlap in the posteriors is for β, while the other three parameters show a more significant difference for planets and brown dwarfs. In particular, for wide-separation giant planets this distribution is more likely to have lower companion masses and higher stellar host mass, with a higher overall occurrence rate compared to brown dwarfs.

## 6.2. Comparing to Demographics of RV-detected Planets

The RV planet-search method represents one of the more robust ways to measure the underlying distributions of giant planets with periods ≲10 yr. While *Kepler*'s transiting planets

greatly outnumber RV-detected planets, the RV method does not suffer as much bias toward longer periods, especially with precision RV surveys operating for multiple decades (e.g., Bonfils et al. 2013; Fischer et al. 2014; Astudillo-Defru et al. 2017). Cumming et al. (2008, C08) presented an analysis of planets orbiting a uniform sample of Keck data of FGK target stars, and fitted a double power-law model to mass *m* and period *P* (similar in form to our Equation (2)), finding that $d^2N/(dP dm) \propto P^{-0.74}m^{-1.31}$, with an overall planet occurrence rate of 10.5% between 2 and 2000 days and 0.3 and 10 $M_{Jup}$. This period distribution is equivalent to a semimajor axis (a) distribution of $dN/da \propto a^{-0.61}$ over a range of 0.031–3.1 au for solar-mass stars. Past direct imaging surveys have sought to determine the extent to which this distribution holds for wider-separation planets (Nielsen & Close 2010; Nielsen et al. 2013; Brandt et al. 2014). While our occurrence rate *f* gives the number of planets per star, C08 only consider the planet with the largest Doppler amplitude in each system, and therefore calculate the number of planetary systems per star. However, C08 note that this is not a strong bias, as only ~10% of their systems have multiple detected RV planets, and we proceed to use the C08 occurrence rate as an estimate of the number of planets per star.
<figure><img src="image_18.png" /><figcaption>**Figure 15.** Multiple mock GPIES surveys in a universe following the Cumming et al. (2008) population of giant planets, extended out to 100 au, with no occurrence rate dependence on stellar host mass. Cumulative histograms give the GPIES-detected planets (in red) and the mock surveys (in blue), for mass, separation, and stellar mass. Thicker blue lines show the median and $1\sigma$ confidence region for these cumulative distributions. The final panel notes the predicted number of planets from these mock surveys, compared to the six planets that were actually detected. While the Cumming et al. (2008) distribution is consistent in number of planets and mass distribution, it somewhat overpredicts the separations of detected planets, and requires additional model parameters to account for the strong stellar mass dependence seen in the GPIES planet occurrence rate.</figcaption></figure>

Figure 7 is inconsistent with C08’s value of $\alpha$ (planet mass index) of –1.31 at the $1\sigma$ level, and excludes their value of $\beta$ (semimajor axis index) of –0.61 at $\sim 2\sigma$. If we insist on the C08 values of $\alpha, \beta$, and $f$, and extend the distribution out to wider separations, we would predict a planet fraction of 4.7% between 5–13 $M_{Jup}$ and 10–100 au—this is also disfavored by GPIES at the $1\sigma$ level. A more significant difference, though, is in host stellar mass, since the Cumming et al. (2008) sample consisted entirely of planets around FGK stars, while the GPIES planets are exclusively orbiting higher-mass stars.

We further investigate the extent to which the C08 RV distribution is consistent with the GPIES results by setting up a Monte Carlo simulation where planets are drawn from the C08 power laws and occurrence rate, but extrapolated out to an

upper semimajor axis cut-off of 100 au, to account for planets seen out to ~70 au (HD 95086 b, HR 8799 b). For now, we assume that the Cumming et al. (2008) distribution fitted to FGK stars applies to stars of all masses, so that planet frequency is not a function of stellar mass ($\gamma = 0$), and $\alpha$ and $\beta$ are the same for lower-mass and higher-mass stars. We use the same method that generated the GPIES completeness maps to produce an ensemble of simulated planets drawn from a distribution using the C08 values of $\alpha, \beta$, and $f$ around each target star, and determine which planets are detectable given the contrast curve for that star (Nielsen et al. 2008; Nielsen & Close 2010). We then record the mass, projected separation, and stellar host mass for each detectable planet from the ensemble. Finally, we conduct 10,000 mock GPIES surveys,
<figure><img src="image_19.png" /><figcaption>**Figure 16.** Another way to compare the Cumming et al. (2008) distributions with GPIES data, this time distinguishing between stars below $1.5 M_\odot$ (left) and above (right). We keep the mass power-law index ($\alpha$) and normalization ($f$) between 0.031 and 3.1 au fixed at the C08 values, and allow the semimajor axis power-law index ($\beta$) and the semimajor axis upper cut-off to find their most probable values given the GPIES data. The color bar gives the probability of a given combination of $\beta$ and cut-off, with white areas being ruled out, and darker regions allowed. The value of $\beta$ exhibited by RV planets is -0.61 (Cumming et al. 2008) and is shown by a vertical solid line. This value (and the C08 values of $\alpha = -1.31$ and $f = 10.5\%$ that are adopted by construction) could apply to our sample of lower-mass stars if the semimajor axis distribution is cut off at less than 34 au ($2\sigma$) or less than 15 au ($1\sigma$). For higher-mass stars, when adopting the Cumming et al. (2008) values of $\alpha, \beta$, and $f$, a value of semimajor axis cut-off of 264 au is consistent at $1\sigma$, while a cut-off of 85 au is needed to be consistent at $2\sigma$. We emphasize that this analysis is distinct from the Bayesian method used in Figure 7, where the cut-off was fixed, and the values of $\alpha, \beta$, and $f$ were all allowed to float. Additionally, while the MCMC analysis explicitly fitted for the observed separations of the planets, this method considers only occurrence rate.</figcaption></figure>

randomly assigning each star a chance of having a planet in the given mass and semimajor axis range, then using the per-star completeness to determine whether that planet would have been detected.

Figure 15 compares these mock surveys to the results from GPIES. Shown in red are cumulative histograms of the six planets detected by GPIES in planet mass, projected separation, and stellar host mass. Overlaid are these same parameters as measured from the mock surveys. There is generally good agreement in planet mass, as expected from the aforementioned overlap between our measured value of $\alpha$ and that of Cumming et al. (2008). Separation also appears consistent, though the GPIES distribution tends to slightly smaller values of separation than the median of the mock survey, again as we would expect from the GPIES-only fit for $\beta$ finding a more negative value than for the RV planets. The largest departure, however, is clearly for stellar mass, where the median mock survey finds planets around stars over a range of stellar masses (by construction), whereas GPIES observes all planet hosts to be the higher-mass stars in the sample. So while the overall occurrence rate of planets appears similar to predictions of the C08 distribution extended to 100 au, the breakdown by stellar mass is significantly different.

To investigate the dependence of these results on the adopted upper cut-off in semimajor axis, which is currently set to 100 au, and to better account for the difference in stellar masses, we divide our sample into lower-mass stars ($<1.5 M_\odot$) and higher-mass stars ($\ge 1.5 M_\odot$), and conduct a second series of mock surveys for the two subsamples. This time we keep the value of $\alpha$ and the total occurrence rate within 3.1 au constant using the Cumming et al. (2008) values, and vary the value of the semimajor axis power-law index ($\beta$) and the upper cut-off, which sets the maximum semimajor axis planets in the distribution may have. For each combination of these two parameters, we compute the expected number of planets GPIES

should have detected using our Monte Carlo simulations. We then use the Poisson distribution to derive the probability that we achieved our actual number of detections in each subsample, displaying the results in Figure 16.

For lower-mass stars ($<1.5 M_\odot$), with a null result for planets, we can only set upper limits for the two parameters, by taking the probability of detecting 0 planets given the expectation value. For the Cumming et al. (2008) value of $\beta = -0.61$, we find the upper cut-off must be less than 15 and 34 au at 1 and $2\sigma$, respectively. Thus it is plausible that for solar-mass stars the Cumming et al. (2008) distribution extends into the giant planet region of our own solar system, as our GPIES sensitivity does not extend inside 10 au.

For higher-mass stars ($>1.5 M_\odot$), since we have detected six planets, we can better constrain the values of $\beta$ and the upper cut-off within this modified Cumming et al. (2008) distribution. We again use the Poisson distribution to compute the probability of detecting six planets around the GPIES higher-mass stars at each point in the grid, then draw contours encircling 68%, 95%, and 99.7% of the total probability. We set a prior that imposes a lower limit on the semimajor axis cut-off of 70 au, corresponding to HR 8799 b (Wang et al. 2018b). Though we only detected HR 8799 cde given the field of view of GPI, we assert that a distribution fitted to the inner three planets of the system should not exclude the outer planet b. The Cumming et al. (2008) value of $\beta = -0.61$ intersects the $1\sigma$ contour at 264 au, and intersects the $2\sigma$ contour at a semimajor axis cut-off of 85 au. Such a large value for the upper cut-off is difficult to reconcile with the relative lack of detections of planets around higher-mass stars beyond ~100 au by previous generations of imaging surveys. An extrapolation of the C08 distribution would predict an equal number of planets between 10 and 80 au (where all six GPIES planets reside) and between 80 and 264 au. The stated errors of C08 in the period
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Mock GPIES surveys, again based on the Cumming et al. (2008) population of giant planets, extended out to 100 au, but only for the 123 higher-mass stars in the survey, and with planet occurrence rate *f* boosted by a factor of two relative to the value reported by C08. With *f* thus boosted, the number of actually detected planets is consistent with the mock surveys. The planet mass distributions also appear consistent. However, there is a significant discrepancy in the separation distribution, with GPIES observing planets at closer separations than the majority of the mock surveys.</figcaption></figure>

distribution ($\frac{dN}{d\ln P} \propto P^{0.26 \pm 0.1}$), if they are Gaussian, would correspond to errors in $\beta$ for semimajor axis of $-0.61 \pm 0.15$, and so a $1\sigma$ value of $-0.46$ would be more consistent with our results, reaching the $1\sigma$ value in the right panel of Figure 16 at 79 au. Adopting such a value for $\beta$ would make the upper limits for lower-mass stars more stringent, moving from 15 and 34 au at 1 and $2\sigma$ to 12 and 23 au, respectively. Previously we had fitted for all three parameters in Equation (2) simultaneously ($\alpha, \beta, f$) while fixing the value of the upper cut-off to 100 au. Here, as we fix $\alpha$ and $f$ to the Cumming et al. (2008) values and let $\beta$ and the upper cut-off float, we reach a similar conclusion, namely that GPIES excludes C08's value of $\beta$ at $\sim 1\sigma$, and that adjusting the value of the upper cut-off is not a straightforward path to removing this discrepancy.

The upper cut-off analysis above considers only the occurrence rate when evaluating the consistency of the extended Cumming et al. (2008) distribution with our GPIES data. In order to consider not just occurrence rate but also the properties of the detected planets, we repeat the mock surveys on only the 123 stars above $1.5 M_\odot$, assuming an upper semimajor axis cut-off of 100 au, and the C08 values of $\alpha$ and $\beta$. Given this modified distribution and our completeness to planets around our 123 higher-mass stars, the mean number of planets found is only 2.6. To more closely match our actual observed number of six planets, we increase the <3.1 au planet occurrence rate from the Cumming et al. (2008) value by a factor of two (which for a $1.5 M_\odot$ star would represent a value of $\gamma$ near our median value from Figure 9 of 2.03). This moves the mean number of detections up to 5.2, with 42% of the mock
surveys now producing six or more detected planets. As before, we assume there is no stellar mass dependence in planet occurrence within this subsample of higher-mass stars ($\gamma = 0$).

Results from this third set of mock surveys restricted to higher-mass stars ($>1.5 M_\odot$) are plotted in Figure 17. While the planet mass distribution and overall occurrence rate appear consistent with the observed GPIES detections, there remain significant differences for the projected separation and stellar host mass distributions. Again, the Cumming et al. (2008) separation distributions place more planets at wider separation than are observed. We saw the same disagreement with the Cumming et al. (2008) distribution in our full power-law fit (Figure 7), where we found a more negative value of $\beta$ from the GPIES sample. Indeed, 80.2% of mock surveys have a detected planet at a projected separation beyond 53.65 au, the value for HD 95086 b. We could try to reduce this discrepancy by boosting $f$ still further and decreasing the upper cut-off in the semimajor axis distribution. But lowering the cut-off would exclude, without physical justification, observed planets like HD 95086 b in our sample at 55 au, as well as the widest-separation planet in the HR 8799 system, planet b at 70 au. Though HR 8799 b does not appear in our sample, the packed nature of the system (Wang et al. 2018a) leads us to reject a model fit to the inner three planets that does not fit the fourth planet.

Figure 17 also shows that despite our restriction to $>1.5 M_\odot$ stars, a mismatch persists in the stellar mass distributions: the mock surveys predict planets to be detected around more $2 M_\odot$ stars than are actually observed. For such high-mass stars in the GPIES sample, there is no obvious bias where more massive stars are systematically older or more distant. The persistent discrepancy in stellar mass distributions likely arises because the C08 distribution places about as many planets at large orbital distances as at small ones, and those at large distances are detectable in the mock surveys of higher-mass stars. In the actual GPI survey, as the absolute H magnitude decreases by ~1 mag from 1.5 to $2 M_\odot$, planets like 51 Eri b, $\beta$ Pic b, and HR 8799 e that cluster at small orbital separations become undetectable around higher-mass stars.

In sum, we find that the semimajor axis distribution from C08 of RV-discovered giant planets is a poor fit to directly imaged giant planets at wider separations. For lower-mass stars ($<1.5 M_\odot$), while we cannot exclude the possibility that the Cumming et al. (2008) RV population extends to orbital distances of 15–34 au, we do not have any data to support this model, either. For higher-mass stars, the data from GPIES place more planets at smaller semimajor axes than at larger ones; GPIES prefers a semimajor axis power-law index of $\beta \simeq -2$, which is discrepant from the C08 value of $-0.61$ at $2\sigma$.

Fernandes et al. (2019) analyze RV and *Kepler* giant planet occurrence rates, and observe a turnover in occurrence rate as a function of semimajor axis, with a peak at ~3 au. Their symmetric *epos* fit gives a value of the planet mass power law $\alpha = -1.45 \pm 0.05$ (similar to the C08 value of $\alpha = -1.31 \pm 0.2$), and a broken power-law fit to semimajor axis, with $\beta = -0.025_{-0.22}^{+0.30}$ for small periods and $\beta = -1.975_{-0.30}^{+0.22}$ for large periods, again converting a period power-law index into semimajor axis. At small semimajor axis, then, Fernandes et al. (2019) find a more positive value of power-law index compared to the C08 value of $\beta = -0.61 \pm 0.15$, and a value of $\beta$ similar to our fit to planets around higher-mass stars

($\beta = -2.0_{-0.5}^{+0.6}$) and planets around stars with a stellar mass term ($\beta = -1.7_{-0.5}^{+0.6}$). For the mass power-law index, we find a more negative value of $\alpha = -2.3_{-0.6}^{+0.8}$, and while the $1\sigma$ confidence intervals overlap, the offset could be a sign of a different planet mass distribution at larger separations, or between higher-mass and lower-mass stars. Fernandes et al. (2019) find the location of the period break from their fit at $P = 1581_{-392}^{+894}$ days, which corresponds to a semimajor axis for solar-type stars of $2.66_{-0.46}^{+0.93}$ au. The asymmetric *epos* fit gives similar values for $\alpha, \beta$ within the break, and overall occurrence rate, as well as consistent values for $\beta$ beyond the break and the location of the break, though with significantly larger error bars: $\beta = -2.8_{-1.9}^{+1.4}$ and a break at $3.2_{-1.4}^{+1.1}$ au. We find support for the conclusions of Fernandes et al. (2019), that the power-law index $\beta$ measured by C08 cannot continue to large values of semimajor axis, but must turn over within $\lesssim 10$ au (modulo the difference in stellar mass between the planet hosts in the C08 and Fernandes et al. 2019 samples and GPIES).

Finally, we combine RV results for small-separation planets and GPIES results at larger separations to estimate the total number of planets within 100 au, as shown in Figure 18. In the top two panels, we couple GPIES and the C08 occurrence rate to calculate the number of planets as a function of semimajor axis between 1 and 13 $M_{Jup}$ around $1 M_\odot$ stars. Draws are taken from the MCMC chain of our fit to giant planets including a stellar mass term, and combined with Monte Carlo draws of $\alpha, \beta$, and occurrence rate from C08, assuming each is Gaussian distributed and that the three parameters are independent, thus producing two power-law distributions. We attempt to make the combined distribution continuous by finding a point where the two power laws cross between 1 and 20 au; if the two do not intersect in this range, the handover point is taken to be 10 au, with a discontinuity in the final combined distribution. Finally, we integrate the distribution to find the posterior on the total number of giant planets. We find a total occurrence rate (number of planets per star) of $0.09 \pm 0.02$ from 0.031 to 10 au, and $0.04_{-0.03}^{+0.09}$ from 10 to 100 au, producing a total occurrence rate within 100 au of $0.14_{-0.05}^{+0.10}$. Repeating this analysis with the Fernandes et al. (2019) *epos* symmetric distribution produces similar results, with occurrence rates of $0.07_{-0.02}^{+0.03}$, $0.04_{-0.03}^{+0.09}$, and $0.12_{-0.05}^{+0.10}$, respectively.

As GPIES only detected planets around higher-mass stars, the posterior for giant planet occurrence approaches zero between 10 and 100 au, though the GPIES data do set a clear upper limit given the lack of detections of giant planets around these lower-mass stars. By assuming this power-law model, we find a $2\sigma$ upper limit on the occurrence rate of 1–13 $M_{Jup}$ planets around $1 M_\odot$ stars of 39.0% and 38.1% for the C08 and Fernandes et al. (2019) distributions, respectively.

For higher-mass stars, we repeat this analysis, this time using the Johnson et al. (2010a) occurrence rate, stellar mass power law ($\gamma$) and errors, and setting the mass power-law $\alpha$ equal to the C08 value and uncertainties, and allowing $\beta$ to follow a uniform distribution between -1 and 0. We then draw from our MCMC fits to the high-mass star GPIES sample and continue as before. Extrapolating the Johnson et al. (2010a) results gives an occurrence rate within 10 au for $1.75 M_\odot$ stars of $0.42_{-0.12}^{+0.20}$. Given the large occurrence rate determined by GPIES for 5–13 $M_{Jup}$, and the very negative value of $\gamma$, the total occurrence rate for 1–13 $M_{Jup}$ planets between 10 and 100 au orbiting higher-mass stars is even larger, $1.0_{-0.6}^{+1.7}$, with a total occurrence rate within 100 au of $1.5_{-0.7}^{+1.7}$ for these higher-mass stars. Given the negative value of the planet mass power-law
<figure><img src="image_21.png" /><figcaption>**Figure 18.** Estimates of the total occurrence rate of giant planets between 1–13 $M_{Jup}$ and 0.031–100 au, from combining RV occurrence rates and power laws for close-in giant planets with those presented in this work for wider-separation giant planets. Occurrence rates from GPIES are extrapolated to 1 $M_{⊙}$ to compare to the C08 and Fernandes et al. (2019) results for solar-type stars, as well as extrapolated down to 1 $M_{Jup}$, using the power laws from Table 3. From this analysis, we estimate the fraction of giant planets within 100 au of solar-type stars to be $0.14_{-0.05}^{+0.10}$ and $0.12_{-0.05}^{+0.10}$ for the C08 and Fernandes et al. (2019) distributions, while for higher-mass stars ($>1.5 M_{⊙}$), we find a significantly larger occurrence rate of $1.5_{-0.7}^{+1.7}$. As we have limited sensitivity to $\lesssim 3 M_{Jup}$ planets, future observations will be required to determine if the power laws seen by GPIES do indeed extend to these low planet masses, and thus more robustly determine the giant planet occurrence rate.</figcaption></figure>

index ($\alpha = -2.4_{-0.8}^{+0.9}$) found from GPIES, many of these predicted planets have small masses. With GPIES having limited sensitivity to $\lesssim 3 M_{Jup}$ planets orbiting higher-mass stars, we cannot rule out a break in the power-law distribution at lower planet masses, which would result in a smaller occurrence rate.

## 6.3. Connecting to Planet Formation Theory

Our finding that giant planets outnumber brown dwarfs at more than 10-to-1 around $>1.5 M_{⊙}$ stars argues against these

planets having formed by gravitational instability. Fragmentation of marginally Toomre-unstable disks (Toomre 1964; Goldreich & Lynden-Bell 1965) is thought to produce a mass distribution weighted toward larger and not smaller masses. When fragments initially condense, they may already be at least as massive as Jupiter (e.g., Figure 3 of Kratter & Lodato 2016; see also Rafikov 2005). In a population synthesis study predicated on gravitational instability, Forgan & Rice (2013) found the overwhelming majority of clumps that survived tidal disruption to lie above the deuterium burning limit. These initial fragments are expected to accrete still more
mass from their strongly self-gravitating, dynamically active parent disks, which themselves may still be accreting from infalling progenitor clouds (Kratter et al. 2010).

In addition, the requirement that disks rapidly cool in order to fragment (Gammie 2001) predicts that collapse occurs at large stellocentric distances where material is less optically thick (e.g., outside at least ~50 au in a model inspired by the HR 8799 system; Figure 3 of Kratter et al. 2010; see also Matzner & Levin 2005 and Rafikov 2009). The generic predictions of gravitational instability—a top-heavy mass distribution extending to brown dwarfs and stars, and large orbital separations (modulo orbital migration)—are not borne out by our 2–13 $M_{Jup}$ planets located between 10 and 60 au. Instead we infer mass distributions (using the fit to planets around higher-mass stars) that are bottom-heavy ($dN/dm \propto m^\alpha$ with $\alpha = -2.4 \pm 0.8$) and separation distributions that favor small semimajor axes ($dN/da \propto a^\beta$ with $\beta = -2.0 \pm 0.5$).

Formation by gravitational instability remains an option for GPIES brown dwarfs, for which we infer mass and separation indices of $\alpha = -0.5 \pm 1.1$ (more top-heavy) and $\beta = -0.7_{-1.5}^{+1.0}$ (weighted more toward large distances). The three brown dwarfs in our survey are located at orbital distances of 10–30 au, somewhat smaller than the values cited above for fragmentation. However, previous surveys with wider fields of view have found brown dwarfs at larger separations—for example, HR 7329 B at 200 au (Lowrance et al. 2000), $\zeta$ Del B at 910 au (De Rosa et al. 2014), and HIP 79797 Ba/Bb, a pair of brown dwarfs in a 3 au orbit separated by 370 au from the primary star (Huélamo et al. 2010; Nielsen et al. 2013). Moe & Kratter (2018) present mechanisms (e.g., Lidov-Kozai cycles) for stellar/brown dwarf binary systems to decrease their separation over time.

Another observed property of brown dwarf companions that we have highlighted is their lack of preference for either higher-mass or lower-mass stars as host. This, too, is consistent with formation within gravitationally unstable disks. If the disk gas surface density $\Sigma_g \propto M_*$ and the sound speed $c_s \propto T^{1/2} \propto M_*^{1/4}$ (where we have taken the gas temperature $T \propto L_*^{1/4}$ and used the PMS luminosity-mass relation $L_* \propto M_*^2$), then Toomre's $Q \propto 1/M_*^{1/4}$. This weak dependence on stellar mass is made even weaker when we fold in the large scatter in observed gas disk masses for given $M_*$ (e.g., Andrews et al. 2018). Deeper analyses that combine the $Q$ criterion with the requirement of fast cooling reveal that disks gravitationally fragment at orbital locations that are remarkably insensitive to system parameters related to stellar mass (Matzner & Levin 2005; Rafikov 2009).

Our finding that giant planets at 10–100 au are more commonly hosted by A stars than by lower-mass stars extends the trend determined by prior demographic studies that closer-in gas giants are more common around FGK stars than around M dwarfs (Clanton & Gaudi 2014; Mulders et al. 2015; Ghezzi et al. 2018). The correlation with stellar mass presents a major constraint on theory. Having argued above against planet formation by gravitational instability (a “top-down” scenario), we ask here whether the trend with stellar mass is consistent with core accretion, a “bottom-up” scenario whereby giants nucleate from sufficiently massive rocky cores. Note that the high luminosities of directly imaged planets, though consistent with hot-start models and inconsistent with cold start models (Section 5.5), do not necessarily rule out core accretion. Cold and hot-start models differ according to the planet’s assumed initial entropy, which can vary significantly within core

accretion, depending on the shock dynamics of runaway gas accretion (see, e.g., Figure 12 of Berardo et al. 2017).

In a prescient analysis, Laughlin et al. (2004) anticipated fewer giant planets around lower-mass stars, arguing that core formation times would be longer around M dwarfs, and potentially too long compared to gas disk dispersal times, because their disks have longer orbital periods and are expected to have fewer solids. We update this reasoning within the emerging paradigm of “pebble accretion,” in which cores grow by accreting small solid particles aerodynamically entrained in disk gas (see Ormel 2017 and Johansen & Lambrechts 2017 for reviews). As particles drift radially inward through the gas disk (e.g., Weidenschilling 1977), a fraction of them are accreted by the protocore. The final core mass is a monotonic increasing function of $\int_{\text{drift}} \Sigma_s dt$: the disk solid surface density $\Sigma_s$ (local to the protocore) integrated over the time it takes solids to drift past the protocore’s orbit and drain out from the disk (see Equations (36), (A7), and (A9) of Lin et al. 2018). There are dependencies of the final core mass on variables apart from this integral, but they are unlikely to change the qualitative conclusion that pebble accretion breeds more massive cores around more massive stars.<fnref n="46" /> The surface density $\Sigma_s$ scales with the total solid disk mass, which should increase with $M_*$ (a linear dependence is commonly assumed). The drift time $t_{\text{drift}}$ also increases with $M_*$, as inferred empirically from millimeter-wave dust continuum observations showing that disks are systematically larger around higher-mass stars (Pascucci et al. 2016, in particular their Section 6.2; see also Andrews et al. 2018). Thus pebble/core accretion gets at least the sign right with respect to the correlation of giant planet occurrence with host star mass, and it remains to work out more quantitative models that can satisfy the constraints reported in our Table 3.<fnref n="47" />

What about semimajor axis trends in the context of pebble accretion? The preference reported here for giant planets to be located at the lower end of the 10–100 au range makes sense insofar as protocores at smaller orbital distances access larger reservoirs of solids lying exterior to their orbits and drifting past. Pebble accretion probabilities also tend to be higher closer in because disk gas is denser there and more efficiently drags pebbles onto cores: accretion at small distances is more likely to be in the so-called “settling” regime, with particle trajectories more strongly gravitationally focused, as opposed to the less efficient “hyperbolic” regime that obtains farther out (Ormel & Klahr 2010; Lin et al. 2018, in particular their Figure 15 which contrasts fast settling accretion at 5 au with slow hyperbolic accretion at 40 au).

We can also try connecting the observed preference of directly imaged giant planets for smaller semimajor axes

<footnote n="46">⁴⁶ In Equations (A7) and (A9) of Lin et al. (2018) pertaining to the “2D” limit where the pebble accretion radius exceeds the solid disk scale height, the extra variable dependencies only strengthen the trend with stellar mass. In the case of their Equation (36) (“3D”), the extra variables mute the trend but do not reverse it, unless the particle Stokes number $St$ (aerodynamic stopping time normalized to the orbit time) decreases more rapidly than $1/M_*^{3/2}$. Pascucci et al. (2016) support $St \propto 1/M_*^{1/2}$ for fragmentation-limited particle growth (Birnstiel et al. 2012); this relation steepens to $1/M_*$ if particle sizes are fixed and gas disk masses scale linearly with stellar mass.</footnote>

<footnote n="47">⁴⁷ This is not to say that there are not outstanding problems intrinsic to the theory. Lin et al. (2018) point out that pebble drift times in Jupiter-breeding scenarios do not match observations. There may also be insufficient mass in pebbles (Manara et al. 2018). Size distributions of accreting particles are crucial but uncertain; the origin of seed cores is tied to the perennial problem of planetesimal formation. See Johansen & Lambrechts (2017) for some current thinking.</footnote>
(which we reiterate applies to distances of 10–100 au) to the observed preference of giant planets discovered by RVs for larger semimajor axes, which applies to distances of 0.03–3 au (e.g., Cumming et al. 2008; Fernandes et al. 2019). Transit surveys report a similar rise in giant planets out to ~300 days (e.g., Dong & Zhu 2013; Petigura et al. 2018). Assuming these semimajor axis trends can be stitched together—keeping in mind the caveat that they derive from disjoint samples of higher-mass and lower-mass stars—we infer the giant planet occurrence rate peaks between ~3 and 10 au. Remarkably, a similar conclusion has been found by Fernandes et al. (2019) in a combined RV-transit analysis of FGK stars: they find evidence for a turnover in the giant planet occurrence rate at ~2–3 au. A giant planet occurrence rate that peaks at intermediate distances aligns at least qualitatively with pebble/core accretion theory. The aforementioned gains in growth with smaller orbital distance are eventually reversed, at the shortest distances, by other competing physical effects. Cores stop accreting pebbles when they gravitationally perturb the surrounding disk gas so strongly that the inward drift of pebbles stalls (Lambrechts et al. 2014). Such “pebble isolation” sets an upper mass limit that decreases as the orbital distance decreases (Bitsch et al. 2018; Fung & Lee 2018); cores do not grow past sub-Earth masses inside ~0.1 au (see Figure 12 of Lin et al. 2018, who rule out in situ formation of hot Jupiters by pebble accretion on this basis). Another effect stifling growth at the shortest distances is that disk gas becomes hotter and more opaque there, and thus cools and accretes onto cores more slowly (Lee & Chiang 2015, 2016). All the above theoretical reasons, together with the many arguments pointing to inward migration of hot Jupiters by Kozai cycles/tidal friction and perhaps also disk torques (Dawson & Johnson 2018), lead us to posit a “sweet spot” for gas giant formation at stellocentric distances that are neither too large nor too small. From Figures 8 and 18, this sweet spot appears to be at ~2–10 au.

## 7. Conclusions

From the first 300 stars observed out of the planned 600-star survey, reaching contrasts of $10^6$ within 1" radius, GPIES is one of the largest and deepest direct imaging surveys for exoplanets conducted to date. Our analysis of the data shows that there is a clear stellar mass dependence on planet occurrence rate, with stars $>1.5 M_\odot$ more likely to host giant planets (5–13 $M_{Jup}$) at wide separations (semimajor axes 10–100 au) than lower-mass stars. Around higher-mass stars, the total occurrence rate of such planets is $9_{-4}^{+5}\%$. When fitting our data for higher-mass stars ($>1.5 M_\odot$) with a power-law model of planet populations where $d^2N/(dm\,da) \propto m^\alpha a^\beta$, we find negative power-law indices of $\alpha \simeq -2$ and $\beta \simeq -2$, values that favor smaller masses and smaller semimajor axes. All the statistical trends that we have uncovered—a bottom-heavy companion mass distribution, orbital distances concentrated more toward 10 au than 100 au, and a strong preference for higher-mass host stars—appear consistent with a “bottom-up” formation scenario: core/pebble accretion.

By comparison to giant planets, brown dwarfs appear to exhibit different statistics: a more top-heavy mass distribution ($\alpha \simeq -0.47$), orbital distances that are somewhat more concentrated toward 100 than 10 au ($\beta \simeq -0.65$), and no particular preference for host stellar mass. These are consistent with a “top-down” formation scenario: gravitational instability and fragmentation of circumstellar disks. Our finding of a

brown dwarf occurrence rate on the order of 1%, independent of stellar mass, is consistent with previous literature analyses of brown dwarfs from direct imaging surveys.

While planets around GPIES stars are more common around higher-mass stars ($>1.5 M_\odot$), it is interesting to note that the most massive planet host in the sample is $\beta$ Pictoris ($1.73 M_\odot$), yet there are 76 stars more massive than $1.8 M_\odot$ in our 300-star sample. Repeating the analysis of Section 5.2, but dividing the higher-mass sample between 1.5 and $1.8 M_\odot$ and $>1.8 M_\odot$, we find at 99.7% confidence that the intermediate stellar mass bin has more planets than the highest stellar mass bin. But this result is dependent on our choice of bin boundary—when we extend the dividing line from 1.8 to 2.1 to $2.4 M_\odot$, our confidence decreases from 99.7% to 97.1% to 87.1%, respectively. It is thus difficult to know if we are observing a real effect or whether we are fine-tuning our boundary based on a small number of detections. For comparison, the limit we used here to separate higher-mass and lower-mass stars, $1.5 M_\odot$, must be reduced to $1.04 M_\odot$ (with only 67 stars out of 300 below this limit) before the confidence with which we can conclude that higher-mass stars have more planetary systems than lower-mass stars drops below 95%. If the frequency of wide-separation giant planet systems does indeed drop off beyond $\sim 2 M_\odot$, this would suggest that the mechanism that formed the giant planets in the GPIES sample becomes less efficient at stellar masses above this value. Jones et al. (2016) observe reminiscent behavior for close-in giant planets around giant stars, with a planet occurrence rate that increases out to $\sim 2.3 M_\odot$, but then decreases at higher stellar masses. Jones et al. (2016) model this as a Gaussian, though like us they have fewer targets in their sample at the high stellar mass end. Additionally, their RV targets have evolved off the main sequence, while ours are predominantly on the main sequence, and so it is unclear how to disentangle formation effects from post-main sequence evolution effects. Nevertheless, further observations of additional higher-mass stars ($\gtrsim 2 M_\odot$) can test whether this trend is real.

There are a limited number of bright, young, nearby stars available to high-contrast direct imagers like GPI and SPHERE. As a result, at the current achievable contrasts there are unlikely to be a significant number of new detections of imaged giant planets following the completion of the GPIES and SPHERE/SHINE (Feldt et al. 2017) campaigns. An improvement in performance from instrument upgrades (Chilcote et al. 2018), however, could unlock mass/separation phase space around these target stars that is currently inaccessible. For example, 51 Eridani was observed multiple times by planet-hunting surveys (e.g., Heinze et al. 2010; Biller et al. 2013; Rameau et al. 2013), yet the planet 51 Eri b was not discovered until it was imaged with the more advanced GPI (Macintosh et al. 2015). Additional planets discovered with upgraded GPI and SPHERE would allow us to test the robustness of the trends with stellar mass and planetary mass uncovered in this paper. The Gaia mission will identify astrometric signatures of planets with $\gtrsim 10$ yr orbits that may be confirmed with direct imaging, which can also boost the number of companions in this separation range.

Alternatively, high-contrast imaging at longer wavelengths offers the opportunity to expand the list of potential imaging targets to older stars (~300 Myr–1 Gyr). While a $5 M_{Jup}$ planet is 100 times brighter at 30 Myr than 600 Myr at H band, this ratio is closer to a factor of 10 at $M'$ (BT-Settl, Baraffe et al.
2015). For a uniform star formation rate, there should be nine times more stars between 100 Myr and 1 Gyr than there are stars younger than 100 Myr. *James Webb Space Telescope*, GMT, TMT, and ELT will be well-suited to probe wide-separation giant planet populations around these intermediate-age stars at 5–10 µm. Going forward, combining different detection methods to produce a full accounting of giant planet and brown dwarf populations from 0.01 to 1000 au will provide still more constraints on formation mechanisms.

We thank the referee for helpful comments that improved the quality of this manuscript. Based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). This research has made use of the SIMBAD and VizieR databases, operated at CDS, Strasbourg, France. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562.

J.R., R.D., and D.L. acknowledge support from the Fonds de Recherche du Québec. J.R.M.'s work was performed in part under contract with the California Institute of Technology (Caltech)/Jet Propulsion Laboratory (JPL) funded by NASA

through the Sagan Fellowship Program executed by the NASA Exoplanet Science Institute. Support for M.M.B.'s work was provided by NASA through Hubble Fellowship grant no. 51378.01-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. Supported by NSF grants AST-1411868 (E.L.N., K.B. F., B.M., J.P., and J.H.), AST-141378 (G.D.), and AST-1518332 (R.D.R., J.J.W., T.M.E., J.R.G., P.G.K.). Supported by NASA grants NNX14AJ80G (E.L.N., S.C.B., B.M., F.M., and M.P.), NNX15AC89G and NNX15AD95G (B.M., J.E.W., T.M.E., R.J.D.R., G.D., J.R.G., P.G.K.), NN15AB52I (D.S.), and NNX16AD44G (K.M.M). K.W.D. is supported by an NRAO Student Observing Support Award SOSPA3-007. J.J. W. is supported by the Heising-Simons Foundation 51 Pegasi b postdoctoral fellowship. E.C. thanks Courtney Dressing for helpful discussions. Portions of this work were performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. This work benefited from NASA's Nexus for Exoplanet System Science (NExSS) research coordination network sponsored by NASA's Science Mission Directorate. The Center for Exoplanets and Habitable Worlds is supported by the Pennsylvania State University, the Eberly College of Science, and the Pennsylvania Space Grant Consortium.

**Facility:** Gemini: South (GPI).

**Software:** Astropy (Astropy Collaboration et al. 2013), emcee (Foreman-Mackey et al. 2013), GPI DRP (Perrin et al. 2016), IDL Astronomy Library (Landsman 1993), pyKLIP (Wang et al. 2015).

# Appendix
## Properties of Target Stars

Table 4 describes the 300 target stars discussed here, including age, distance, mass, and details of the GPIES observations.
<table><caption>Table 4
300-star Sample Properties and Observing Log</caption><thead><tr><th>Name</th><th>R.A.<br>(hh:mm:ss.ss)</th><th>Decl.<br>(dd:mm:ss.ss)</th><th>SpT</th><th>Ref.</th><th>d<br>(pc)</th><th>Ref.</th><th>I<br>(mag)</th><th>H<br>(mag)</th><th>Moving Group</th><th>Age<br>(Myr)</th><th>Ref.</th><th>Mass<br>(M<sub>☉</sub>)</th><th>UT Date</th><th>t<sub>exp</sub> × n<sub>coadd</sub><br>(s)</th><th>t<sub>int</sub><br>(minutes)</th><th>ΔPA<br>(deg)</th></tr></thead><tbody><tr><td>BU Psc</td><td>00:02:02.59</td><td>-02:45:58.2</td><td>A8/9V</td><td>1</td><td>64.70 ± 0.38</td><td>2</td><td>6.5</td><td>6.3</td><td>...</td><td>554</td><td>3</td><td>1.49</td><td>20160920</td><td>59.6 × 1</td><td>34.8</td><td>21.3</td></tr><tr><td>HD 105</td><td>00:05:52.54</td><td>-41:45:11.0</td><td>G0V</td><td>4</td><td>38.845 ± 0.084</td><td>2</td><td>6.9</td><td>6.2</td><td>TucHor</td><td>45</td><td>5</td><td>1.14</td><td>20141110</td><td>59.6 × 1</td><td>38.8</td><td>31.0</td></tr><tr><td>HR 9</td><td>00:06:50.08</td><td>-23:06:27.1</td><td>F3V</td><td>4</td><td>39.960 ± 0.098</td><td>2</td><td>5.8</td><td>5.3</td><td>β Pic</td><td>26</td><td>6</td><td>1.49</td><td>20141109</td><td>59.6 × 1</td><td>34.8</td><td>44.4</td></tr><tr><td>HD 987</td><td>00:13:52.83</td><td>-74:41:17.5</td><td>G8V</td><td>4</td><td>45.847 ± 0.051</td><td>2</td><td>8.0</td><td>7.1</td><td>TucHor</td><td>45</td><td>5</td><td>0.94</td><td>20151105</td><td>59.6 × 1</td><td>28.8</td><td>16.6</td></tr><tr><td>HD 984</td><td>00:14:10.25</td><td>-07:11:56.8</td><td>F7V</td><td>1</td><td>45.91 ± 0.12</td><td>2</td><td>6.8</td><td>6.2</td><td>Columba</td><td>42</td><td>5</td><td>1.28</td><td>20150830</td><td>59.6 × 1</td><td>21.9</td><td>15.1</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20160919</td><td>59.6 × 1</td><td>13.9</td><td>9.5</td></tr><tr><td>HD 1237 A</td><td>00:16:12.70</td><td>-79:51:04.0</td><td>G8V</td><td>4</td><td>17.5649 ± 0.0073</td><td>2</td><td>5.9</td><td>5.0</td><td>Carina-Near</td><td>200</td><td>7</td><td>1.01</td><td>20150831</td><td>59.6 × 1</td><td>34.8</td><td>13.6</td></tr><tr><td>HD 1466</td><td>00:18:26.10</td><td>-63:28:39.0</td><td>F8V</td><td>4</td><td>42.970 ± 0.055</td><td>2</td><td>6.9</td><td>6.2</td><td>TucHor</td><td>45</td><td>5</td><td>1.19</td><td>20150703</td><td>59.6 × 1</td><td>42.7</td><td>21.0</td></tr><tr><td>9 Cet</td><td>00:22:51.79</td><td>-12:12:34.0</td><td>G2.5V</td><td>8</td><td>21.310 ± 0.023</td><td>2</td><td>5.8</td><td>5.0</td><td>...</td><td>300</td><td>9</td><td>1.14</td><td>20150901</td><td>59.6 × 1</td><td>36.8</td><td>31.2</td></tr><tr><td>κ Phe</td><td>00:26:12.12</td><td>-43:40:47.7</td><td>A5IVn</td><td>10</td><td>23.809 ± 0.085</td><td>11</td><td>3.7</td><td>3.6</td><td>...</td><td>578</td><td>3</td><td>1.83</td><td>20150702</td><td>52.4 × 1</td><td>41.0</td><td>47.5</td></tr><tr><td>λ<sup>1</sup> Phe</td><td>00:31:24.98</td><td>-48:48:12.7</td><td>A1Va</td><td>12</td><td>52.97 ± 0.67</td><td>11</td><td>4.5</td><td>4.8</td><td>...</td><td>473</td><td>3</td><td>2.18</td><td>20150704</td><td>59.6 × 1</td><td>41.8</td><td>35.3</td></tr><tr><td>β<sup>1</sup> Tuc</td><td>00:31:32.67</td><td>-62:57:29.6</td><td>B8/A0</td><td>13</td><td>41.41 ± 0.34</td><td>11</td><td>4.4</td><td>4.7</td><td>TucHor</td><td>45</td><td>5</td><td>2.44</td><td>20151106</td><td>59.6 × 1</td><td>37.8</td><td>18.6</td></tr><tr><td>HD 3221</td><td>00:34:51.20</td><td>-61:54:58.0</td><td>K4Ve</td><td>4</td><td>44.499 ± 0.050</td><td>2</td><td>8.4</td><td>6.7</td><td>TucHor</td><td>45</td><td>5</td><td>1.06</td><td>20160918</td><td>59.6 × 1</td><td>35.8</td><td>19.2</td></tr><tr><td>HD 3375</td><td>00:36:00.80</td><td>-59:43:01.9</td><td>F2V</td><td>13</td><td>57.099 ± 0.095</td><td>2</td><td>6.4</td><td>6.1</td><td>...</td><td>30</td><td>9</td><td>1.51</td><td>20151105</td><td>59.6 × 1</td><td>33.8</td><td>22.5</td></tr><tr><td>HD 3821 A</td><td>00:40:47.48</td><td>-07:13:56.6</td><td>G1V</td><td>4</td><td>26.496 ± 0.033</td><td>2</td><td>6.4</td><td>5.6</td><td>...</td><td>300</td><td>9</td><td>1.05</td><td>20160921</td><td>59.6 × 1</td><td>34.8</td><td>24.5</td></tr><tr><td>HD 3888</td><td>00:40:51.58</td><td>-53:12:35.7</td><td>F7V</td><td>13</td><td>44.177 ± 0.060</td><td>2</td><td>6.8</td><td>6.2</td><td>...</td><td>70</td><td>9</td><td>1.24</td><td>20150801</td><td>59.6 × 1</td><td>12.9</td><td>8.6</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20150831</td><td>59.6 × 1</td><td>36.8</td><td>25.6</td></tr><tr><td>ksi Phe</td><td>00:41:46.30</td><td>-56:30:05.2</td><td>A3VpSrCr Ksn</td><td>14</td><td>68.43 ± 0.42</td><td>2</td><td>5.5</td><td>5.6</td><td>...</td><td>767</td><td>3</td><td>1.91</td><td>20150729</td><td>59.6 × 1</td><td>17.9</td><td>12.0</td></tr><tr><td>HD 5294</td><td>00:54:59.27</td><td>+24:06:01.0</td><td>G5</td><td>15</td><td>29.917 ± 0.045</td><td>2</td><td>6.7</td><td>5.9</td><td>...</td><td>300</td><td>9</td><td>1.03</td><td>20160921</td><td>59.6 × 1</td><td>35.8</td><td>12.3</td></tr><tr><td>HIP 6276</td><td>01:20:32.27</td><td>-11:28:03.7</td><td>G9V</td><td>4</td><td>35.332 ± 0.060</td><td>2</td><td>7.7</td><td>6.7</td><td>AB Dor</td><td>149</td><td>5</td><td>0.94</td><td>20160920</td><td>59.6 × 1</td><td>35.8</td><td>30.4</td></tr><tr><td>HD 8558</td><td>01:23:21.14</td><td>-57:28:50.5</td><td>G7V</td><td>4</td><td>45.314 ± 0.053</td><td>2</td><td>7.8</td><td>6.9</td><td>TucHor</td><td>45</td><td>5</td><td>1.01</td><td>20160919</td><td>59.6 × 1</td><td>37.8</td><td>23.8</td></tr><tr><td>CC Phe</td><td>01:28:08.66</td><td>-52:38:19.1</td><td>K1V</td><td>4</td><td>39.844 ± 0.044</td><td>2</td><td>8.3</td><td>6.9</td><td>TucHor</td><td>45</td><td>5</td><td>0.86</td><td>20160918</td><td>59.6 × 1</td><td>35.8</td><td>24.9</td></tr><tr><td>48 Cet</td><td>01:29:36.10</td><td>-21:37:45.5</td><td>A0V</td><td>16</td><td>71.38 ± 1.48</td><td>11</td><td>4.8</td><td>5.0</td><td>...</td><td>553</td><td>3</td><td>2.27</td><td>20141112</td><td>59.6 × 1</td><td>39.8</td><td>63.3</td></tr><tr><td>49 Cet</td><td>01:34:37.72</td><td>-15:40:34.9</td><td>A1V</td><td>16</td><td>57.07 ± 0.33</td><td>2</td><td>5.4</td><td>5.5</td><td>Argus</td><td>50</td><td>17</td><td>2.02</td><td>20141111</td><td>59.6 × 1</td><td>16.9</td><td>17.4</td></tr><tr><td>EX Cet</td><td>01:37:35.37</td><td>-06:45:36.7</td><td>K0/1V</td><td>1</td><td>24.044 ± 0.025</td><td>2</td><td>6.9</td><td>5.9</td><td>TucHor</td><td>45</td><td>5</td><td>0.88</td><td>20141215</td><td>59.6 × 1</td><td>6.0</td><td>3.6</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20150901</td><td>59.6 × 1</td><td>30.8</td><td>24.5</td></tr><tr><td>HD 10472</td><td>01:40:24.06</td><td>-60:59:56.6</td><td>F2IV/V</td><td>13</td><td>71.15 ± 0.15</td><td>2</td><td>7.2</td><td>6.7</td><td>TucHor</td><td>45</td><td>5</td><td>1.44</td><td>20151106</td><td>59.6 × 1</td><td>37.8</td><td>19.5</td></tr><tr><td>DK Cet</td><td>01:57:48.91</td><td>-21:54:05.0</td><td>G4V</td><td>4</td><td>41.411 ± 0.092</td><td>2</td><td>7.4</td><td>6.6</td><td>TucHor</td><td>45</td><td>5</td><td>1.08</td><td>20141109</td><td>59.6 × 1</td><td>34.8</td><td>10.1</td></tr><tr><td>HD 13183</td><td>02:07:17.98</td><td>-53:11:56.4</td><td>G7V</td><td>4</td><td>49.535 ± 0.081</td><td>2</td><td>7.9</td><td>7.0</td><td>TucHor</td><td>45</td><td>5</td><td>1.05</td><td>20141215</td><td>59.6 × 1</td><td>38.8</td><td>33.8</td></tr><tr><td>HD 13246</td><td>02:07:26.02</td><td>-59:40:45.8</td><td>F7V</td><td>4</td><td>45.610 ± 0.056</td><td>2</td><td>7.0</td><td>6.3</td><td>TucHor</td><td>45</td><td>5</td><td>1.23</td><td>20141218</td><td>59.6 × 1</td><td>38.8</td><td>27.4</td></tr><tr><td>φ Eri</td><td>02:16:30.60</td><td>-51:30:44.0</td><td>B8IV</td><td>18</td><td>47.12 ± 0.27</td><td>11</td><td>3.5</td><td>4.0</td><td>TucHor</td><td>45</td><td>5</td><td>3.58</td><td>20150831</td><td>59.6 × 1</td><td>37.8</td><td>27.0</td></tr><tr><td>HIP 11152</td><td>02:23:26.64</td><td>+22:44:06.7</td><td>M1V</td><td>19</td><td>27.193 ± 0.045</td><td>2</td><td>9.3</td><td>7.6</td><td>β Pic</td><td>26</td><td>6</td><td>0.50</td><td>20160918</td><td>59.6 × 1</td><td>37.8</td><td>12.5</td></tr><tr><td>HD 15115</td><td>02:26:16.20</td><td>+06:17:33.6</td><td>F4IV</td><td>20</td><td>49.00 ± 0.10</td><td>2</td><td>6.3</td><td>5.9</td><td>Columba</td><td>42</td><td>5</td><td>1.45</td><td>20141216</td><td>59.6 × 1</td><td>62.6</td><td>35.2</td></tr><tr><td>CC Eri</td><td>02:34:22.60</td><td>-43:47:47.0</td><td>K7V</td><td>4</td><td>11.5373 ± 0.0049</td><td>2</td><td>7.1</td><td>5.1</td><td>...</td><td>30</td><td>9</td><td>0.84<sup>a</sup></td><td>20141109</td><td>59.6 × 1</td><td>34.8</td><td>32.0</td></tr><tr><td>HD 16699 A</td><td>02:38:44.29</td><td>-52:57:03.1</td><td>F8V</td><td>4</td><td>57.876 ± 0.093</td><td>2</td><td>7.2</td><td>6.7</td><td>...</td><td>27</td><td>9</td><td>1.28</td><td>20160919</td><td>59.6 × 1</td><td>33.8</td><td>22.6</td></tr><tr><td>ε Hyi</td><td>02:39:35.23</td><td>-68:16:00.8</td><td>B9Va</td><td>18</td><td>46.55 ± 0.20</td><td>11</td><td>3.9</td><td>4.4</td><td>TucHor</td><td>45</td><td>5</td><td>2.70</td><td>20150729</td><td>59.6 × 1</td><td>47.7</td><td>28.5</td></tr><tr><td>HR 789</td><td>02:39:47.93</td><td>-42:53:29.8</td><td>A1Vb</td><td>10</td><td>40.47 ± 0.64</td><td>2</td><td>4.5</td><td>4.6</td><td>β Pic</td><td>26</td><td>6</td><td>2.17</td><td>20141110</td><td>59.6 × 1</td><td>39.8</td><td>43.6</td></tr><tr><td>HR 784</td><td>02:40:12.42</td><td>-09:27:10.4</td><td>F8VFe-0.4</td><td>10</td><td>21.93 ± 0.12</td><td>2</td><td>5.3</td><td>4.7</td><td>...</td><td>125</td><td>9</td><td>1.81<sup>a</sup></td><td>20141214</td><td>59.6 × 1</td><td>40.8</td><td>32.7</td></tr><tr><td>ζ Hor</td><td>02:40:39.60</td><td>-54:33:00.0</td><td>F6V</td><td>4</td><td>48.94 ± 0.28</td><td>2</td><td>4.7</td><td>4.3</td><td>...</td><td>30</td><td>9</td><td>3.34<sup>a</sup></td><td>20150831</td><td>59.6 × 1</td><td>21.9</td><td>15.7</td></tr><tr><td>84 Cet B</td><td>02:41:13.78</td><td>-00:41:42.0</td><td>...</td><td>...</td><td>22.59 ± 0.43</td><td>11</td><td>6.6</td><td>4.5</td><td>...</td><td>124</td><td>9</td><td>0.73</td><td>20160920</td><td>59.6 × 1</td><td>37.8</td><td>25.5</td></tr><tr><td>t Hor</td><td>02:42:33.47</td><td>-50:48:01.1</td><td>F8V</td><td>4</td><td>17.327 ± 0.016</td><td>2</td><td>4.8</td><td>4.3</td><td>...</td><td>300</td><td>9</td><td>1.24</td><td>20141111</td><td>59.6 × 1</td><td>39.8</td><td>29.6</td></tr><tr><td>HD 17207</td><td>02:45:01.13</td><td>-22:09:58.7</td><td>F7V</td><td>10</td><td>38.860 ± 0.069</td><td>2</td><td>6.7</td><td>6.0</td><td>...</td><td>125</td><td>9</td><td>1.20</td><td>20141112</td><td>59.6 × 1</td><td>38.8</td><td>58.9</td></tr><tr><td>HR 826</td><td>02:46:45.10</td><td>-21:38:22.2</td><td>F3IV/V</td><td>16</td><td>48.198 ± 0.084</td><td>2</td><td>6.1</td><td>5.6</td><td>...</td><td>300</td><td>9</td><td>1.53</td><td>20141109</td><td>59.6 × 1</td><td>30.8</td><td>69.5</td></tr><tr><td></td><td>02:47:27.21</td><td>+19:22:21.2</td><td>G5V</td><td>21</td><td>33.065 ± 0.056</td><td>2</td><td>7.0</td><td>5.7</td><td>AB Dor</td><td>149</td><td>5</td><td>0.97</td><td>20150830</td><td>59.6 × 1</td><td>35.8</td><td>14.4</td></tr></tbody></table>
<table><caption>Table 4
(Continued)</caption><thead><tr><th>Name</th><th>R.A.<br>(hh:mm:ss.ss)</th><th>Decl.<br>(dd:mm:ss.s)</th><th>SpT</th><th>Ref.</th><th>d<br>(pc)</th><th>Ref.</th><th>I<br>(mag)</th><th>H<br>(mag)</th><th>Moving Group</th><th>Age<br>(Myr)</th><th>Ref.</th><th>Mass<br>(M<sub>☉</sub>)</th><th>UT Date</th><th>t<sub>exp</sub> × n<sub>coadd</sub><br>(s)</th><th>t<sub>int</sub><br>(minutes)</th><th>ΔPA<br>(deg)</th></tr></thead><tbody><tr><td colspan="16"><strong>HD 17332 B</strong></td></tr><tr><td>VZ Hor</td><td>02:51:53.00</td><td>-61:37:05.0</td><td>K1V</td><td>4</td><td>40.126 ± 0.048</td><td>2</td><td>7.9</td><td>6.7</td><td>...</td><td>300</td><td>9</td><td>0.99</td><td>20141216</td><td>59.6 × 1</td><td>38.8</td><td>18.6</td></tr><tr><td>EP Eri</td><td>02:52:32.13</td><td>-12:46:11.0</td><td>K1V</td><td>16</td><td>10.3584 ± 0.0047</td><td>2</td><td>5.1</td><td>4.2</td><td>...</td><td>126</td><td>9</td><td>0.92</td><td>20160920</td><td>59.6 × 1</td><td>37.8</td><td>26.8</td></tr><tr><td>HD 18809</td><td>03:00:19.70</td><td>-37:27:16.0</td><td>G4V</td><td>4</td><td>48.89 ± 0.20</td><td>2</td><td>7.8</td><td>6.9</td><td>...</td><td>125</td><td>9</td><td>1.05</td><td>20141111</td><td>59.6 × 1</td><td>39.8</td><td>58.4</td></tr><tr><td>HD 20232</td><td>03:13:50.19</td><td>-38:48:33.4</td><td>A2/3III/IV</td><td>22</td><td>79.47 ± 0.19</td><td>2</td><td>6.5</td><td>6.4</td><td>...</td><td>429</td><td>3</td><td>1.66</td><td>20141215</td><td>59.6 × 1</td><td>39.8</td><td>53.9</td></tr><tr><td>κ<sup>1</sup> Cet</td><td>03:19:21.70</td><td>+03:22:12.7</td><td>G5V</td><td>8</td><td>9.140 ± 0.023</td><td>11</td><td>4.2</td><td>3.0</td><td>...</td><td>300</td><td>9</td><td>1.05</td><td>20141218</td><td>59.6 × 1</td><td>38.8</td><td>20.5</td></tr><tr><td>ε Eri</td><td>03:32:55.84</td><td>-09:27:29.7</td><td>K2V</td><td>8</td><td>3.2161 ± 0.0017</td><td>11</td><td>2.8</td><td>1.9</td><td>...</td><td>500</td><td>9</td><td>0.86</td><td>20141110</td><td>14.5 × 4</td><td>35.0</td><td>29.2</td></tr><tr><td>HR 1139</td><td>03:43:33.84</td><td>-10:29:08.2</td><td>A5III-IV</td><td>14</td><td>41.81 ± 0.22</td><td>2</td><td>5.3</td><td>5.2</td><td>...</td><td>534</td><td>3</td><td>1.62</td><td>20160918</td><td>59.6 × 1</td><td>37.8</td><td>22.2</td></tr><tr><td>HR 1189</td><td>03:48:35.48</td><td>-37:37:19.2</td><td>A1Va</td><td>12</td><td>52.64 ± 0.33</td><td>2</td><td>5.3</td><td>5.4</td><td>TucHor</td><td>45</td><td>5</td><td>2.01</td><td>20151103</td><td>59.6 × 1</td><td>37.8</td><td>60.4</td></tr><tr><td>HR 1190</td><td>03:48:35.88</td><td>-37:37:12.6</td><td>B9.5Van</td><td>12</td><td>52.70 ± 0.48</td><td>2</td><td>4.6</td><td>4.6</td><td>TucHor</td><td>45</td><td>5</td><td>2.53</td><td>20151103</td><td>59.6 × 1</td><td>36.8</td><td>38.7</td></tr><tr><td>HD 24916</td><td>03:57:28.70</td><td>-01:09:34.1</td><td>K4V</td><td>1</td><td>15.301 ± 0.011</td><td>2</td><td>7.0</td><td>5.5</td><td>...</td><td>125</td><td>9</td><td>0.80</td><td>20141214</td><td>59.6 × 1</td><td>37.8</td><td>21.8</td></tr><tr><td>HR 1249</td><td>04:02:36.74</td><td>-00:16:08.1</td><td>F7/8V</td><td>1</td><td>18.771 ± 0.041</td><td>2</td><td>4.9</td><td>4.3</td><td>AB Dor</td><td>149</td><td>5</td><td>1.28</td><td>20141217</td><td>59.6 × 1</td><td>32.8</td><td>18.9</td></tr><tr><td>39 Tau</td><td>04:05:20.26</td><td>+22:00:32.0</td><td>G5V</td><td>23</td><td>16.895 ± 0.020</td><td>2</td><td>5.3</td><td>4.5</td><td>...</td><td>300</td><td>9</td><td>1.11</td><td>20150829</td><td>59.6 × 1</td><td>37.8</td><td>12.6</td></tr><tr><td>HD 25953</td><td>04:06:41.52</td><td>+01:41:02.9</td><td>F6V</td><td>1</td><td>56.97 ± 0.15</td><td>2</td><td>7.3</td><td>6.7</td><td>AB Dor</td><td>149</td><td>5</td><td>1.25</td><td>20151221</td><td>59.6 × 1</td><td>36.8</td><td>17.6</td></tr><tr><td>γ Dor</td><td>04:16:01.59</td><td>-51:29:11.9</td><td>F1V</td><td>10</td><td>20.46 ± 0.15</td><td>11</td><td>3.8</td><td>3.5</td><td>...</td><td>894</td><td>3</td><td>1.56</td><td>20141110</td><td>59.6 × 1</td><td>37.8</td><td>28.6</td></tr><tr><td>HD 28287</td><td>04:26:38.57</td><td>-28:57:06.5</td><td>G9.5Vk:</td><td>10</td><td>35.645 ± 0.039</td><td>2</td><td>8.0</td><td>6.9</td><td>AB Dor</td><td>149</td><td>5</td><td>0.89</td><td>20141112</td><td>59.6 × 1</td><td>38.8</td><td>145.3</td></tr><tr><td>51 Eri</td><td>04:37:36.13</td><td>-02:28:24.8</td><td>F0IV</td><td>14</td><td>29.78 ± 0.12</td><td>2</td><td>4.8</td><td>4.8</td><td>β Pic</td><td>26</td><td>6</td><td>1.58</td><td>20141214</td><td>59.6 × 1</td><td>37.8</td><td>23.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20141218</td><td>59.6 × 1</td><td>37.8</td><td>23.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20150131</td><td>59.6 × 1</td><td>62.6</td><td>36.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20150831</td><td>59.6 × 1</td><td>44.7</td><td>29.9</td></tr><tr><td>V834 Tau</td><td>04:41:18.86</td><td>+20:54:05.4</td><td>K4V</td><td>24</td><td>13.2381 ± 0.0100</td><td>2</td><td>6.7</td><td>5.3</td><td>...</td><td>82</td><td>9</td><td>0.71</td><td>20160919</td><td>59.6 × 1</td><td>37.8</td><td>12.7</td></tr><tr><td>HD 30447</td><td>04:46:49.50</td><td>-26:18:09.0</td><td>F3V</td><td>22</td><td>80.54 ± 0.25</td><td>2</td><td>7.4</td><td>7.0</td><td>Columba</td><td>42</td><td>5</td><td>1.45</td><td>20160920</td><td>59.6 × 1</td><td>31.8</td><td>93.0</td></tr><tr><td>HD 31392</td><td>04:54:04.20</td><td>-35:24:16.2</td><td>G9V</td><td>10</td><td>25.765 ± 0.020</td><td>2</td><td>6.9</td><td>5.9</td><td>...</td><td>200</td><td>9</td><td>0.99</td><td>20141216</td><td>59.6 × 1</td><td>40.8</td><td>62.1</td></tr><tr><td>HR 1597</td><td>04:54:53.04</td><td>-58:32:51.5</td><td>F5V</td><td>10</td><td>30.468 ± 0.028</td><td>2</td><td>5.7</td><td>5.1</td><td>...</td><td>125</td><td>9</td><td>1.33</td><td>20160125</td><td>59.6 × 1</td><td>36.8</td><td>21.4</td></tr><tr><td>π<sup>1</sup> Ori</td><td>04:54:53.71</td><td>+10:09:04.0</td><td>A0Va IB</td><td>12</td><td>35.66 ± 0.32</td><td>11</td><td>4.4</td><td>4.5</td><td>...</td><td>424</td><td>3</td><td>1.85</td><td>20151221</td><td>59.6 × 1</td><td>23.9</td><td>9.9</td></tr><tr><td>HIP 23309</td><td>05:00:47.13</td><td>-57:15:25.5</td><td>M0Ve</td><td>4</td><td>26.901 ± 0.020</td><td>2</td><td>8.1</td><td>6.4</td><td>β Pic</td><td>26</td><td>6</td><td>0.79</td><td>20151106</td><td>59.6 × 1</td><td>37.8</td><td>23.1</td></tr><tr><td>ζ Dor</td><td>05:05:30.69</td><td>-57:28:22.7</td><td>F9VFe-0.5</td><td>10</td><td>11.646 ± 0.024</td><td>11</td><td>4.1</td><td>3.4</td><td>...</td><td>300</td><td>9</td><td>1.18</td><td>20151105</td><td>18.9 × 3</td><td>35.3</td><td>20.8</td></tr><tr><td>AS Col</td><td>05:20:38.02</td><td>-39:45:18.0</td><td>F6V</td><td>4</td><td>47.663 ± 0.073</td><td>2</td><td>6.9</td><td>6.2</td><td>Columba</td><td>42</td><td>5</td><td>1.28</td><td>20141111</td><td>59.6 × 1</td><td>38.8</td><td>43.6</td></tr><tr><td>HD 35650</td><td>05:24:30.14</td><td>-38:58:10.2</td><td>K6V</td><td>4</td><td>17.4793 ± 0.0086</td><td>2</td><td>7.9</td><td>6.1</td><td>AB Dor</td><td>149</td><td>5</td><td>0.71</td><td>20141110</td><td>59.6 × 1</td><td>39.8</td><td>24.7</td></tr><tr><td>HIP 25434</td><td>05:26:22.98</td><td>-43:22:36.4</td><td>G0</td><td>25</td><td>92.41 ± 0.26</td><td>2</td><td>8.3</td><td>7.9</td><td>Columba</td><td>42</td><td>5</td><td>1.23</td><td>20151221</td><td>59.6 × 1</td><td>36.8</td><td>45.0</td></tr><tr><td>HD 35841</td><td>05:26:36.59</td><td>-22:29:23.7</td><td>F3V</td><td>16</td><td>103.68 ± 0.30</td><td>2</td><td>8.4</td><td>7.8</td><td>Columba</td><td>42</td><td>5</td><td>1.30</td><td>20160228</td><td>59.6 × 1</td><td>46.7</td><td>46.1</td></tr><tr><td>AF Lep</td><td>05:27:04.76</td><td>-11:54:03.5</td><td>F8V(n)k:</td><td>10</td><td>26.877 ± 0.019</td><td>2</td><td>5.7</td><td>5.1</td><td>β Pic</td><td>26</td><td>6</td><td>1.15</td><td>20141214</td><td>59.6 × 1</td><td>28.8</td><td>25.8</td></tr><tr><td>HD 36435</td><td>05:27:39.35</td><td>-60:24:57.6</td><td>G5V</td><td>4</td><td>19.478 ± 0.012</td><td>2</td><td>6.3</td><td>5.3</td><td>...</td><td>125</td><td>9</td><td>0.97</td><td>20141217</td><td>59.6 × 1</td><td>39.8</td><td>22.0</td></tr><tr><td>HD 36329</td><td>05:29:24.10</td><td>-34:30:56.0</td><td>G3V</td><td>4</td><td>72.46 ± 0.15</td><td>2</td><td>7.8</td><td>7.1</td><td>Columba</td><td>42</td><td>5</td><td>1.35</td><td>20141111</td><td>59.6 × 1</td><td>35.8</td><td>40.4</td></tr><tr><td>ζ Lep</td><td>05:46:57.35</td><td>-14:49:19.0</td><td>A2IV-V(n)</td><td>10</td><td>21.608 ± 0.075</td><td>11</td><td>3.3</td><td>3.3</td><td>β Pic</td><td>26</td><td>6</td><td>2.06</td><td>20160229</td><td>29.1 × 2</td><td>37.8</td><td>37.5</td></tr><tr><td>β Pic</td><td>05:47:17.08</td><td>-51:04:00.2</td><td>A6V</td><td>10</td><td>19.440 ± 0.045</td><td>11</td><td>3.6</td><td>3.5</td><td>β Pic</td><td>26</td><td>6</td><td>1.73</td><td>20151106</td><td>59.6 × 1</td><td>36.8</td><td>28.0</td></tr><tr><td>HD 40781</td><td>05:59:13.38</td><td>-38:42:38.1</td><td>G0V</td><td>4</td><td>64.19 ± 0.71</td><td>2</td><td>7.7</td><td>7.0</td><td>...</td><td>127</td><td>9</td><td>1.21</td><td>20160227</td><td>59.6 × 1</td><td>37.8</td><td>60.4</td></tr><tr><td>HD 41278</td><td>06:00:55.34</td><td>-54:57:05.0</td><td>F5V</td><td>13</td><td>57.201 ± 0.083</td><td>2</td><td>6.9</td><td>6.4</td><td>AB Dor</td><td>149</td><td>5</td><td>1.37</td><td>20151217</td><td>59.6 × 1</td><td>36.8</td><td>27.0</td></tr><tr><td>HR 2234</td><td>06:15:29.66</td><td>-04:54:52.7</td><td>A3Vas</td><td>26</td><td>73.28 ± 0.45</td><td>2</td><td>5.8</td><td>5.9</td><td>...</td><td>440</td><td>3</td><td>1.91</td><td>20150130</td><td>59.6 × 1</td><td>40.8</td><td>26.3</td></tr><tr><td>AO Men</td><td>06:18:28.21</td><td>-72:02:41.5</td><td>K4Ve</td><td>4</td><td>39.263 ± 0.039</td><td>2</td><td>8.4</td><td>7.0</td><td>β Pic</td><td>26</td><td>6</td><td>0.89</td><td>20151221</td><td>59.6 × 1</td><td>35.8</td><td>16.6</td></tr><tr><td>V1358 Ori</td><td>06:19:08.04</td><td>-03:26:20.0</td><td>G0V</td><td>1</td><td>52.01 ± 0.15</td><td>2</td><td>7.3</td><td>6.6</td><td>Columba</td><td>42</td><td>5</td><td>1.24</td><td>20150129</td><td>59.6 × 1</td><td>29.8</td><td>18.1</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20151219</td><td>59.6 × 1</td><td>30.8</td><td>18.4</td></tr><tr><td>AB Pic</td><td>06:19:12.89</td><td>-58:03:15.8</td><td>K1V(e)</td><td>4</td><td>50.120 ± 0.073</td><td>2</td><td>8.3</td><td>7.1</td><td>Carina</td><td>45</td><td>5</td><td>1.00</td><td>20151106</td><td>59.6 × 1</td><td>35.8</td><td>19.5</td></tr><tr><td>V435 Car</td><td>06:21:50.05</td><td>-51:14:15.8</td><td>A5V</td><td>27</td><td>71.43 ± 0.17</td><td>2</td><td>6.4</td><td>6.3</td><td>...</td><td>439</td><td>3</td><td>1.64</td><td>20141217</td><td>59.6 × 1</td><td>39.8</td><td>29.9</td></tr><tr><td>HD 45270</td><td>06:22:30.96</td><td>-60:13:07.7</td><td>G1V</td><td>4</td><td>23.890 ± 0.014</td><td>2</td><td>5.9</td><td>5.2</td><td>AB Dor</td><td>149</td><td>5</td><td>1.15</td><td>20141216</td><td>59.6 × 1</td><td>39.8</td><td>22.7</td></tr></tbody></table>
<table><caption>Table 4
(Continued)</caption><thead><tr><th>Name</th><th>R.A.<br>(hh:mm:ss.ss)</th><th>Decl.<br>(dd:mm:ss.s)</th><th>SpT</th><th>Ref.</th><th>d<br>(pc)</th><th>Ref.</th><th>I<br>(mag)</th><th>H<br>(mag)</th><th>Moving Group</th><th>Age<br>(Myr)</th><th>Ref.</th><th>Mass<br>(M<sub>☉</sub>)</th><th>UT Date</th><th>t<sub>exp</sub> × n<sub>coadd</sub><br>(s)</th><th>t<sub>int</sub><br>(minutes)</th><th>ΔPA<br>(deg)</th></tr></thead><tbody><tr><td>HIP 31878</td><td>06:39:50.06</td><td>-61:28:42.2</td><td>K7V(e)</td><td>4</td><td>22.241 ± 0.014</td><td>2</td><td>8.4</td><td>6.6</td><td>AB Dor</td><td>149</td><td>5</td><td>0.69</td><td>20160128</td><td>59.6 × 1</td><td>33.8</td><td>17.4</td></tr><tr><td>HD 48524</td><td>06:42:05.52</td><td>-38:00:13.7</td><td>G0V</td><td>22</td><td>67.84 ± 0.49</td><td>2</td><td>7.8</td><td>7.1</td><td>...</td><td>125</td><td>9</td><td>1.17</td><td>20141214</td><td>59.6 × 1</td><td>33.8</td><td>74.7</td></tr><tr><td>26 Gem</td><td>06:42:24.33</td><td>+17:38:43.1</td><td>A2V</td><td>15</td><td>43.63 ± 1.28</td><td>11</td><td>5.1</td><td>5.1</td><td>Columba</td><td>42</td><td>5</td><td>1.87</td><td>20160226</td><td>59.6 × 1</td><td>37.8</td><td>14.0</td></tr><tr><td>HD 48370</td><td>06:43:01.02</td><td>-02:53:19.3</td><td>K0V(+G)</td><td>1</td><td>36.072 ± 0.068</td><td>2</td><td>7.2</td><td>6.4</td><td>...</td><td>125</td><td>9</td><td>1.06</td><td>20150201</td><td>59.6 × 1</td><td>39.8</td><td>24.2</td></tr><tr><td>HD 49855</td><td>06:43:46.32</td><td>-71:58:35.8</td><td>G6V</td><td>4</td><td>55.59 ± 0.10</td><td>2</td><td>8.2</td><td>7.4</td><td>Carina</td><td>45</td><td>5</td><td>0.93</td><td>20160320</td><td>59.6 × 1</td><td>36.8</td><td>15.3</td></tr><tr><td>HR 2562</td><td>06:50:01.01</td><td>-60:14:57.9</td><td>F5VFe+0.4</td><td>10</td><td>34.041 ± 0.048</td><td>2</td><td>5.7</td><td>5.1</td><td>...</td><td>300</td><td>9</td><td>1.43</td><td>20160125</td><td>59.6 × 1</td><td>33.8</td><td>19.4</td></tr><tr><td>HR 2558</td><td>06:51:42.45</td><td>-36:13:48.5</td><td>A3V</td><td>22</td><td>55.41 ± 0.17</td><td>2</td><td>5.6</td><td>5.5</td><td>...</td><td>687</td><td>3</td><td>1.69</td><td>20151217</td><td>59.6 × 1</td><td>35.8</td><td>74.6</td></tr><tr><td>HD 51797</td><td>06:56:23.55</td><td>-46:46:55.1</td><td>K0V(e)</td><td>4</td><td>95.65 ± 0.24</td><td>2</td><td>8.9</td><td>7.9</td><td>Columba</td><td>42</td><td>5</td><td>1.20</td><td>20160321</td><td>59.6 × 1</td><td>30.8</td><td>99.4</td></tr><tr><td>HD 53143</td><td>06:59:59.66</td><td>-61:20:10.2</td><td>G9V</td><td>4</td><td>18.3602 ± 0.0079</td><td>2</td><td>6.1</td><td>5.1</td><td>IC2391</td><td>50</td><td>17</td><td>1.00</td><td>20151219</td><td>59.6 × 1</td><td>36.8</td><td>20.3</td></tr><tr><td>HR 2720</td><td>07:12:04.08</td><td>-30:49:17.0</td><td>A8III/IV</td><td>22</td><td>47.423 ± 0.087</td><td>2</td><td>5.7</td><td>5.5</td><td>...</td><td>724</td><td>3</td><td>1.54</td><td>20160327</td><td>59.6 × 1</td><td>36.8</td><td>170.6</td></tr><tr><td>λ Gem A</td><td>07:18:05.61</td><td>+16:32:25.7</td><td>A4IV</td><td>26</td><td>30.93 ± 0.19</td><td>11</td><td>3.4</td><td>3.5</td><td>...</td><td>638</td><td>3</td><td>2.12</td><td>20160319</td><td>29.1 × 2</td><td>34.9</td><td>14.6</td></tr><tr><td>HR 2813</td><td>07:20:21.40</td><td>-52:18:41.0</td><td>F5V</td><td>10</td><td>33.15 ± 0.24</td><td>2</td><td>5.5</td><td>5.1</td><td>Carina-Near</td><td>200</td><td>7</td><td>1.43</td><td>20160320</td><td>59.6 × 1</td><td>35.8</td><td>27.0</td></tr><tr><td>HR 2814</td><td>07:20:21.86</td><td>-52:18:33.3</td><td>F9V</td><td>10</td><td>36.65 ± 0.59</td><td>2</td><td>5.9</td><td>5.3</td><td>Carina-Near</td><td>200</td><td>7</td><td>1.35</td><td>20141217</td><td>59.6 × 1</td><td>38.8</td><td>28.8</td></tr><tr><td>HD 57969</td><td>07:20:22.98</td><td>-56:17:41.3</td><td>A1V</td><td>13</td><td>72.20 ± 0.17</td><td>2</td><td>6.3</td><td>6.4</td><td>...</td><td>217</td><td>3</td><td>1.71</td><td>20151221</td><td>59.6 × 1</td><td>31.8</td><td>32.3</td></tr><tr><td>HD 58192</td><td>07:23:53.76</td><td>-17:24:48.2</td><td>F7V</td><td>16</td><td>36.594 ± 0.053</td><td>2</td><td>6.3</td><td>5.7</td><td>...</td><td>300</td><td>9</td><td>1.24</td><td>20160324</td><td>59.6 × 1</td><td>29.8</td><td>48.2</td></tr><tr><td>HD 58556</td><td>07:25:57.22</td><td>-02:14:54.5</td><td>G1V</td><td>1</td><td>33.925 ± 0.052</td><td>2</td><td>6.4</td><td>5.7</td><td>...</td><td>300</td><td>9</td><td>1.19</td><td>20160328</td><td>59.6 × 1</td><td>36.8</td><td>21.1</td></tr><tr><td>HD 59704</td><td>07:29:31.40</td><td>-38:07:21.0</td><td>F7V</td><td>22</td><td>53.370 ± 0.097</td><td>2</td><td>7.2</td><td>6.5</td><td>Carina-Near</td><td>200</td><td>7</td><td>1.26</td><td>20160326</td><td>59.6 × 1</td><td>22.9</td><td>62.5</td></tr><tr><td>HR 2882</td><td>07:30:42.57</td><td>-37:20:22.1</td><td>G3V</td><td>4</td><td>21.770 ± 0.013</td><td>2</td><td>6.0</td><td>5.3</td><td>...</td><td>200</td><td>9</td><td>1.06</td><td>20141111</td><td>59.6 × 1</td><td>38.8</td><td>76.3</td></tr><tr><td>HD 61033</td><td>07:34:28.03</td><td>-52:58:05.4</td><td>G7V(e)</td><td>4</td><td>28.93 ± 0.38</td><td>11</td><td>6.8</td><td>5.8</td><td>...</td><td>300</td><td>9</td><td>1.05</td><td>20160324</td><td>59.6 × 1</td><td>36.8</td><td>22.9</td></tr><tr><td>HD 61005</td><td>07:35:47.50</td><td>-32:12:14.7</td><td>G8Vk</td><td>10</td><td>36.486 ± 0.043</td><td>2</td><td>7.5</td><td>6.6</td><td>Argus</td><td>50</td><td>17</td><td>0.98</td><td>20150130</td><td>59.6 × 1</td><td>40.8</td><td>109.5</td></tr><tr><td>HIP 37288</td><td>07:39:23.04</td><td>+02:11:01.2</td><td>M0.0Ve</td><td>28</td><td>15.058 ± 0.012</td><td>2</td><td>8.1</td><td>6.1</td><td>...</td><td>300</td><td>9</td><td>0.67</td><td>20141218</td><td>59.6 × 1</td><td>38.8</td><td>22.6</td></tr><tr><td>HD 62850</td><td>07:42:36.10</td><td>-59:17:51.0</td><td>G3V</td><td>4</td><td>32.912 ± 0.028</td><td>2</td><td>6.6</td><td>5.9</td><td>Carina-Near</td><td>200</td><td>7</td><td>1.15</td><td>20160318</td><td>59.6 × 1</td><td>37.8</td><td>19.5</td></tr><tr><td>HD 63008</td><td>07:44:12.51</td><td>-50:27:24.2</td><td>F9VFe-0.5</td><td>10</td><td>30.084 ± 0.022</td><td>2</td><td>6.1</td><td>5.4</td><td>...</td><td>300</td><td>9</td><td>1.19</td><td>20160322</td><td>59.6 × 1</td><td>35.8</td><td>28.0</td></tr><tr><td>HD 63008 B</td><td>07:44:16.48</td><td>-50:27:59.8</td><td>G5V</td><td>10</td><td>30.080 ± 0.027</td><td>2</td><td>6.9</td><td>6.0</td><td>...</td><td>300</td><td>9</td><td>1.01</td><td>20160322</td><td>59.6 × 1</td><td>35.8</td><td>27.2</td></tr><tr><td>YZ CMi</td><td>07:44:40.17</td><td>+03:33:08.8</td><td>M4.0Ve</td><td>29</td><td>5.9874 ± 0.0021</td><td>2</td><td>8.3</td><td>6.0</td><td>...</td><td>70</td><td>9</td><td>0.29</td><td>20160225</td><td>59.6 × 1</td><td>32.8</td><td>22.5</td></tr><tr><td>HD 62883</td><td>07:45:50.88</td><td>-07:31:46.4</td><td>F5V</td><td>1</td><td>58.23 ± 0.23</td><td>2</td><td>6.8</td><td>6.2</td><td>...</td><td>30</td><td>9</td><td>1.26</td><td>20150201</td><td>59.6 × 1</td><td>40.8</td><td>43.3</td></tr><tr><td>HD 63608</td><td>07:46:17.00</td><td>-59:48:34.0</td><td>K0IV-V(k)</td><td>10</td><td>33.806 ± 0.030</td><td>2</td><td>7.5</td><td>6.5</td><td>Carina-Near</td><td>200</td><td>7</td><td>0.97</td><td>20160125</td><td>59.6 × 1</td><td>32.8</td><td>21.0</td></tr><tr><td>HR 3070</td><td>07:49:12.90</td><td>-60:17:01.0</td><td>F4V</td><td>10</td><td>33.767 ± 0.074</td><td>2</td><td>5.2</td><td>4.9</td><td>Carina-Near</td><td>200</td><td>7</td><td>1.56</td><td>20150408</td><td>59.6 × 1</td><td>41.8</td><td>22.4</td></tr><tr><td>HD 64122</td><td>07:49:29.35</td><td>-54:54:04.4</td><td>F6/7V</td><td>13</td><td>80.66 ± 0.19</td><td>2</td><td>7.9</td><td>7.3</td><td>...</td><td>30</td><td>9</td><td>1.32</td><td>20141214</td><td>59.6 × 1</td><td>32.8</td><td>25.1</td></tr><tr><td>TYC 206-419-1</td><td>08:07:09.09</td><td>+07:23:00.1</td><td>K8</td><td>15</td><td>41.469 ± 0.068</td><td>2</td><td>8.9</td><td>7.3</td><td>AB Dor</td><td>149</td><td>5</td><td>0.82</td><td>20160319</td><td>59.6 × 1</td><td>26.8</td><td>12.2</td></tr><tr><td>HD 68143</td><td>08:11:15.15</td><td>+01:16:36.4</td><td>F3/5+F6/7</td><td>1</td><td>80.15 ± 1.37</td><td>2</td><td>8.1</td><td>7.4</td><td>...</td><td>125</td><td>9</td><td>1.25</td><td>20160226</td><td>59.6 × 1</td><td>34.8</td><td>19.2</td></tr><tr><td>q Pup</td><td>08:18:33.31</td><td>-36:39:33.4</td><td>A8V</td><td>10</td><td>28.63 ± 0.15</td><td>11</td><td>4.2</td><td>4.1</td><td>...</td><td>894</td><td>3</td><td>1.71</td><td>20160321</td><td>59.6 × 1</td><td>37.8</td><td>76.6</td></tr><tr><td>HD 70703</td><td>08:21:00.48</td><td>-52:13:40.9</td><td>A0/1V</td><td>27</td><td>72.71 ± 0.14</td><td>2</td><td>6.3</td><td>6.3</td><td>Carina</td><td>45</td><td>5</td><td>1.74</td><td>20151217</td><td>59.6 × 1</td><td>34.8</td><td>30.9</td></tr><tr><td>V478 Hya</td><td>08:22:49.98</td><td>+01:51:34.0</td><td>G1/2V</td><td>1</td><td>59.28 ± 0.16</td><td>2</td><td>8.1</td><td>7.3</td><td>...</td><td>200</td><td>9</td><td>1.08</td><td>20141216</td><td>59.6 × 1</td><td>13.9</td><td>6.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20151220</td><td>59.6 × 1</td><td>9.9</td><td>8.7</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20160328</td><td>59.6 × 1</td><td>35.8</td><td>21.7</td></tr><tr><td>HR 3300</td><td>08:22:55.18</td><td>-52:07:25.6</td><td>A0V</td><td>27</td><td>73.25 ± 0.38</td><td>2</td><td>5.7</td><td>5.9</td><td>Carina</td><td>45</td><td>5</td><td>2.21</td><td>20150407</td><td>59.6 × 1</td><td>41.8</td><td>29.7</td></tr><tr><td>V592 Pup</td><td>08:25:17.72</td><td>-34:22:01.2</td><td>G1IV</td><td>4</td><td>66.98 ± 0.15</td><td>2</td><td>7.1</td><td>6.4</td><td>...</td><td>80</td><td>9</td><td>1.34</td><td>20160227</td><td>59.6 × 1</td><td>44.7</td><td>86.6</td></tr><tr><td>HR 3341</td><td>08:26:25.24</td><td>-52:48:27.0</td><td>A0V</td><td>27</td><td>69.33 ± 0.27</td><td>2</td><td>5.8</td><td>5.9</td><td>AB Dor</td><td>149</td><td>5</td><td>1.96</td><td>20151221</td><td>59.6 × 1</td><td>37.8</td><td>23.1</td></tr><tr><td>HD 72687</td><td>08:33:15.36</td><td>-29:57:23.8</td><td>G5V</td><td>4</td><td>45.436 ± 0.085</td><td>2</td><td>7.6</td><td>6.8</td><td>...</td><td>125</td><td>9</td><td>1.05</td><td>20160229</td><td>59.6 × 1</td><td>35.8</td><td>172.2</td></tr><tr><td>HR 3395</td><td>08:35:50.98</td><td>+06:37:12.8</td><td>F8V</td><td>30</td><td>26.280 ± 0.090</td><td>2</td><td>5.4</td><td>4.8</td><td>...</td><td>125</td><td>9</td><td>2.08<sup>a</sup></td><td>20150404</td><td>59.6 × 1</td><td>32.8</td><td>19.5</td></tr><tr><td>V401 Hya</td><td>08:37:50.46</td><td>-06:48:25.2</td><td>G8/K0(IV)</td><td>1</td><td>24.342 ± 0.026</td><td>2</td><td>6.1</td><td>5.3</td><td>...</td><td>300</td><td>9</td><td>1.11</td><td>20150406</td><td>59.6 × 1</td><td>36.8</td><td>29.0</td></tr><tr><td>HD 74341 A</td><td>08:40:43.54</td><td>-57:32:43.6</td><td>A4V</td><td>31</td><td>72.63 ± 0.19</td><td>2</td><td>6.4</td><td>5.9</td><td>...</td><td>320</td><td>3</td><td>1.72</td><td>20151219</td><td>59.6 × 1</td><td>37.8</td><td>24.6</td></tr></tbody></table>