Population-level Eccentricity Distributions of Imaged Exoplanets and Brown Dwarf Companions: Dynamical Evidence for Distinct Formation Channels

The orbital eccentricities of directly imaged exoplanets and brown dwarf companions provide clues about their formation and dynamical histories. We combine new high-contrast imaging observations of substellar companions obtained primarily with Keck/NIRC2 together with astrometry from the literature to test for differences in the population-level eccentricity distributions of 27 long-period giant planets and brown dwarf companions between 5 and 100 au using hierarchical Bayesian modeling. Orbit fits are performed in a uniform manner for companions with short orbital arcs; this typically results in broad constraints for individual eccentricity distributions, but together as an ensemble, these systems provide valuable insight into their collective underlying orbital patterns. The shape of the eccentricity distribution function for our full sample of substellar companions is approximately flat from e = 0–1. When subdivided by companion mass and mass ratio, the underlying distributions for giant planets and brown dwarfs show significant differences. Low mass ratio companions preferentially have low eccentricities, similar to the orbital properties of warm Jupiters found with radial velocities and transits. We interpret this as evidence for in situ formation on largely undisturbed orbits within massive extended disks. Brown dwarf companions exhibit a broad peak at e ≈ 0.6–0.9 with evidence for a dependence on orbital period. This closely resembles the orbital properties and period-eccentricity trends of wide (1–200 au) stellar binaries, suggesting that brown dwarfs in this separation range predominantly form in a similar fashion. We also report evidence that the "eccentricity dichotomy" observed at small separations extends to planets on wide orbits: the mean eccentricity for the multi-planet system HR 8799 is lower than for systems with single planets. In the future, larger samples and continued astrometric orbit monitoring will help establish whether these eccentricity distributions correlate with other parameters such as stellar host mass, multiplicity, and age

Population-level Eccentricity Distributions of Imaged Exoplanets and Brown Dwarf Companions: Dynamical Evidence for Distinct Formation Channels

The orbital eccentricities of directly imaged exoplanets and brown dwarf companions provide clues about their formation and dynamical histories. We combine new high-contrast imaging observations of substellar companions obtained primarily with Keck/NIRC2 together with astrometry from the literature to test for differences in the population-level eccentricity distributions of 27 long-period giant planets and brown dwarf companions between 5 and 100 au using hierarchical Bayesian modeling. Orbit fits are performed in a uniform manner for companions with short orbital arcs; this typically results in broad constraints for individual eccentricity distributions, but together as an ensemble, these systems provide valuable insight into their collective underlying orbital patterns. The shape of the eccentricity distribution function for our full sample of substellar companions is approximately flat from e = 0–1. When subdivided by companion mass and mass ratio, the underlying distributions for giant planets and brown dwarfs show significant differences. Low mass ratio companions preferentially have low eccentricities, similar to the orbital properties of warm Jupiters found with radial velocities and transits. We interpret this as evidence for in situ formation on largely undisturbed orbits within massive extended disks. Brown dwarf companions exhibit a broad peak at e ≈ 0.6–0.9 with evidence for a dependence on orbital period. This closely resembles the orbital properties and period-eccentricity trends of wide (1–200 au) stellar binaries, suggesting that brown dwarfs in this separation range predominantly form in a similar fashion. We also report evidence that the "eccentricity dichotomy" observed at small separations extends to planets on wide orbits: the mean eccentricity for the multi-planet system HR 8799 is lower than for systems with single planets. In the future, larger samples and continued astrometric orbit monitoring will help establish whether these eccentricity distributions correlate with other parameters such as stellar host mass, multiplicity, and age

LACEwING: A New Moving Group Analysis Code

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars

Orbits for the Impatient: A Bayesian Rejection Sampling Method for Quickly Fitting the Orbits of Long-Period Exoplanets

We describe a Bayesian rejection sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of MCMC in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to ten long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.Comment: 32 pages, 28 figures, Accepted to A

Searching for Scatterers: High-Contrast Imaging of Young Stars Hosting Wide-Separation Planetary-Mass Companions

We have conducted an angular differential imaging survey with NIRC2 at Keck in search of close-in substellar companions to a sample of seven systems with confirmed planetary-mass companions (PMCs) on wide orbits (>50 au). These wide-separation PMCs pose significant challenges to all three possible formation mechanisms: core accretion plus scattering, disk instability, and turbulent fragmentation. We explore the possibility that these companions formed closer in and were scattered out to their present-day locations by searching for other massive bodies at smaller separations. The typical sensitivity for this survey is ΔK ~ 12.5 at 1". We identify eight candidate companions, whose masses would reach as low as one Jupiter mass if gravitationally bound. From our multi-epoch astrometry we determine that seven of these are conclusively background objects, while the eighth near DH Tau is ambiguous and requires additional monitoring. We rule out the presence of >7 M_(Jup) bodies in these systems down to 15–50 au that could be responsible for scattering. This result combined with the totality of evidence suggests that dynamical scattering is unlikely to have produced this population of PMCs. We detect orbital motion from the companions ROXs 42B b and ROXs 12 b, and from this determine 95% upper limits on the companions' eccentricities of 0.58 and 0.83 respectively. Finally, we find that the 95% upper limit on the occurrence rate of additional planets with masses between 5 and 15 M_(Jup) outside of 40 au in systems with PMCs is 54%

LACEwING: A New Moving Group Analysis Code

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars

Searching for Scatterers: High-Contrast Imaging of Young Stars Hosting Wide-Separation Planetary-Mass Companions

We have conducted an angular differential imaging survey with NIRC2 at Keck in search of close-in substellar companions to a sample of seven systems with confirmed planetary-mass companions (PMCs) on wide orbits (>50 au). These wide-separation PMCs pose significant challenges to all three possible formation mechanisms: core accretion plus scattering, disk instability, and turbulent fragmentation. We explore the possibility that these companions formed closer in and were scattered out to their present-day locations by searching for other massive bodies at smaller separations. The typical sensitivity for this survey is ΔK ~ 12.5 at 1". We identify eight candidate companions, whose masses would reach as low as one Jupiter mass if gravitationally bound. From our multi-epoch astrometry we determine that seven of these are conclusively background objects, while the eighth near DH Tau is ambiguous and requires additional monitoring. We rule out the presence of >7 M_(Jup) bodies in these systems down to 15–50 au that could be responsible for scattering. This result combined with the totality of evidence suggests that dynamical scattering is unlikely to have produced this population of PMCs. We detect orbital motion from the companions ROXs 42B b and ROXs 12 b, and from this determine 95% upper limits on the companions' eccentricities of 0.58 and 0.83 respectively. Finally, we find that the 95% upper limit on the occurrence rate of additional planets with masses between 5 and 15 M_(Jup) outside of 40 au in systems with PMCs is 54%

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

We present astrometric monitoring of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager between 2013 and 2016. A small but significant position angle change is detected at constant separation; the orbital motion is confirmed with literature measurements. Efficient Monte Carlo techniques place preliminary constraints on the orbital parameters of HD 95086 b. With 68% confidence, a semimajor axis of 61.7^{+20.7}_{-8.4} au and an inclination of 153.0^{+9.7}_{-13.5} deg are favored, with eccentricity less than 0.21. Under the assumption of a co-planar planet-disk system, the periastron of HD 95086 b is beyond 51 au with 68% confidence. Therefore HD 95086 b cannot carve the entire gap inferred from the measured infrared excess in the SED of HD 95086. We use our sensitivity to additional planets to discuss specific scenarios presented in the literature to explain the geometry of the debris belts. We suggest that either two planets on moderately eccentric orbits or three to four planets with inhomogeneous masses and orbital properties are possible. The sensitivity to additional planetary companions within the observations presented in this study can be used to help further constrain future dynamical simulations of the planet-disk system.Comment: Accepted for publication in ApJ