RadVel is an open source Python package for modeling Keplerian orbits in radial velocity (RV) time series. RadVel provides a convenient framework to fit RVs using maximum a posteriori optimization and to compute robust confidence intervals by sampling the posterior probability density via Markov Chain Monte Carlo (MCMC). RadVel allows users to float or fix parameters, impose priors, and perform Bayesian model comparison. We have implemented realtime MCMC convergence tests to ensure adequate sampling of the posterior. RadVel can output a number of publication-quality plots and tables. Users may interface with RadVel through a convenient command-line interface or directly from Python. The code is object-oriented and thus naturally extensible. We encourage contributions from the community. Documentation is available at http://radvel.readthedocs.io.

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-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 multiplanet 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.

We used high-precision radial velocity measurements of FGKM stars to determine the occurrence of giant planets as a function of orbital separation spanning 0.03-30 au. Giant planets are more prevalent at orbital distances of 1-10 au compared to orbits interior or exterior of this range. The increase in planet occurrence at ∼1 au by a factor of ∼4 is highly statistically significant. A fall-off in giant planet occurrence at larger orbital distances is favored over models with flat or increasing occurrence. We measure -+ 14.1 1.8 2.0 giant planets per 100 stars with semimajor axes of 2-8 au and -+ 8.9 2.4 3.0 giant planets per 100 stars in the range 8-32 au, a decrease in occurrence with increasing orbital separation that is significant at the ∼2σ level. We find that the occurrence rate of sub-Jovian planets (0.1-1 Jupiter masses) is also enhanced for 1-10 au orbits. This suggests that lower-mass planets may share the formation or migration mechanisms that drive the increased prevalence near the water-ice line for their Jovian counterparts. Our measurements of cold gas giant occurrence are consistent with the latest results from direct imaging surveys and gravitational lensing surveys despite different stellar samples. We corroborate previous findings that giant planet occurrence increases with stellar mass and metallicity. Unified Astronomy Thesaurus concepts: Exoplanets (498); Exoplanet astronomy (486); Exoplanet catalogs (488); Surveys (1671); Radial velocity (1332); Exoplanet detection methods (489); Extrasolar gaseous planets (2172); Extrasolar gaseous giant planets (509)

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.