The Impact of Bayesian Hyperpriors on the Population-level Eccentricity Distribution of Imaged Planets

Orbital eccentricities directly trace the formation mechanisms and dynamical histories of substellar companions. Here, we study the effect of hyperpriors on the population-level eccentricity distributions inferred for the sample of directly imaged substellar companions (brown dwarfs and cold Jupiters) from hierarchical Bayesian modeling (HBM). We find that the choice of hyperprior can have a significant impact on the population-level eccentricity distribution inferred for imaged companions, an effect that becomes more important as the sample size and orbital coverage decrease to values that mirror the existing sample. We reanalyze the current observational sample of imaged giant planets in the 5-100 au range from Bowler et al. and find that the underlying eccentricity distribution implied by the imaged planet sample is broadly consistent with the eccentricity distribution for close-in exoplanets detected using radial velocities. Furthermore, our analysis supports the conclusion from that study that long-period giant planets and brown dwarf eccentricity distributions differ by showing that it is robust to the choice of hyperprior. We release our HBM and forward-modeling code in an open-source Python package, ePop!, and make it freely available to the community.