
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 publicationquality plots and tables. Users may interface
with RadVel through a convenient commandline interface or directly from
Python. The code is objectoriented and thus naturally extensible. We encourage
contributions from the community. Documentation is available at
http://radvel.readthedocs.io.

We present evidence that the recently discovered, directlyimaged planet HD
131399 Ab is a background star with nonzero proper motion. From new JHK1L'
photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE,
and Keck/NIRC2, and a reanalysis of the discovery data obtained with
VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The
broader wavelength coverage and higher data quality allow us to reinvestigate
its status. Its nearinfrared spectral energy distribution excludes spectral
types later than L0 and is consistent with a K or M dwarf, which are the most
likely candidates for a background object in this direction at the apparent
magnitude observed. If it were a physically associated object, the projected
velocity of HD 131399 Ab would exceed escape velocity given the mass and
distance to HD 131399 A. We show that HD 131399 Ab is also not following the
expected track for a stationary background star at infinite distance. Solving
for the proper motion and parallax required to explain the relative motion of
HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to
predicted background objects drawn from a galactic model, we find this proper
motion to be high, but consistent with the top 4% fastestmoving background
stars. From our analysis we conclude that HD 131399 Ab is a background K or M
dwarf.

Our Keck/NIRC2 imaging survey searches for stellar companions around 144
systems with radial velocity (RV) detected giant planets to determine whether
stellar binaries influence the planets' orbital parameters. This survey, the
largest of its kind to date, finds eight confirmed binary systems and three
confirmed triple systems. These include three new multistellar systems (HD
30856, HD 86081, and HD 207832) and three multistellar systems with newly
confirmed common proper motion (HD 43691, HD 116029, and HD 164509). We combine
these systems with seven RV planethosting multistellar systems from the
literature in order to test for differences in the properties of planets with
semimajor axes ranging between 0.15 au in single vs multistellar systems. We
find no evidence that the presence or absence of stellar companions alters the
distribution of planet properties in these systems. Although the observed
stellar companions might influence the orbits of more distant planetary
companions in these systems, our RV observations currently provide only weak
constraints on the masses and orbital properties of planets beyond 5 au. In
order to aid future efforts to characterize long period RV companions in these
systems, we publish our contrast curves for all 144 targets. Using four years
of astrometry for six hierarchical triple star systems hosting giant planets,
we fit the orbits of the stellar companions in order to characterize the
orbital architecture in these systems. We find that the orbital plane of the
secondary and tertiary companions are inconsistent with an edgeon orbit in
four out of six cases.

We describe a Bayesian rejection sampling algorithm designed to efficiently
compute posterior distributions of orbital elements for data covering short
fractions of longperiod 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 MetropolisHastings 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 longperiod 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.

We present new observations of the lowmass companion to HD 984 taken with
the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey
campaign. Images of HD 984 B were obtained in the J (1.121.3 micron) and H
(1.501.80 micron) bands. Combined with archival epochs from 2012 and 2014, we
fit the first orbit to the companion to find an 18 AU (70 year) orbit with a
68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a
68% confidence interval between 0.05 and 0.47, and an inclination of 119
degrees with a 68% confidence interval between 114 degrees and 125 degrees. To
address considerable spectral covariance in both spectra, we present a method
of splitting the spectra into low and high frequencies to analyze the spectral
structure at different spatial frequencies with the proper spectral noise
correlation. Using the split spectra, we compare to known spectral types using
field brown dwarf and lowmass star spectra and find a best fit match of a
field gravity M6.5+/1.5 spectral type with a corresponding temperature of
2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun)
= 2.88+/0.07 dex, using DUSTY models. Mass estimates, again from DUSTY
models, find an agedependent mass of 34+/1 to 95+/4 M_Jup. These results are
consistent with previous measurements of the object.