We report the masses, sizes, and orbital properties of 86 planets orbiting 55 stars observed by NASA’s K2 Mission with follow-up Doppler measurements by the HIRES spectrometer at the W. M. Keck Observatory and the Automated Planet Finder at Lick Observatory. Eighty-one of the planets were discovered from their transits in the K2 photometry, while five were found based on subsequent Doppler measurements of transiting planet-host stars. The sizes of the transiting planets range from Earth-size to larger than Jupiter (1–3 R _⊕ is typical), while the orbital periods range from less than a day to a few months. For 32 of the planets, the Doppler signal was detected with significance greater than 5 σ (51 were detected with >3 σ significance). An important characteristic of this catalog is the use of uniform analysis procedures to determine stellar and planetary properties. This includes the transit search and fitting procedures applied to the K2 photometry, the Doppler fitting techniques applied to the radial velocities (RVs), and the spectral modeling to determine bulk stellar parameters. Such a uniform treatment will make the catalog useful for statistical studies of the masses, densities, and system architectures of exoplanetary systems. This work also serves as a data release for all previously unpublished RVs and associated stellar activity indicators obtained by our team for these systems, along with derived stellar and planet parameters.

With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet's low surface gravity and the star's brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4 yr Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of the gas envelope of WASP-107b. We reveal that WASP-107b's mass is only 1.8 Neptune masses (M-b = 30.5 1.7 M-circle plus). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of >85% unless it is substantially inflated. The corresponding core mass of M-circle plus at 3 sigma is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at greater than or similar to 1 au onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day a(b) = 0.0566 0.0017 au. Beyond WASP-107b, we also detect a second, more massive planet ((M-c sin i = 0.36 +/- 0.04M(J)) on a wide eccentric orbit (e(c) = 0.28 +/- 0.07) that may have influenced the orbital migration and spin-orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system's formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion.

We explore the impact of outer stellar companions on the occurrence rate of giant planets detected with radial velocities. We searched for stellar and planetary companions to a volume-limited sample of solar-type stars within 25 pc. Using adaptive optics imaging observations from the Lick 3 m and Palomar 200″ Telescopes, we characterized the multiplicity of our sample stars, down to the bottom of the main sequence. With these data, we confirm field star multiplicity statistics from previous surveys. We additionally combined three decades of radial velocity (RV) data from the California Planet Search with newly collected RV data from Keck/HIRES and the Automated Planet Finder/Levy Spectrometer to search for planetary companions in these same systems. Using an updated catalog of both stellar and planetary companions, as well as detailed injection/recovery tests to determine our sensitivity and completeness, we measured the occurrence rate of planets among the single- and multiple-star systems. We found that planets with masses in the range of 0.1–10 M J and with semimajor axes of 0.1–10 au have an occurrence rate of planets per star when they orbit single stars and an occurrence rate of 0.12 ± 0.04 planets per star when they orbit a star in a binary system. Breaking the sample down by the binary separation, we found that only one planet-hosting binary system had a binary separation <100 au, and none had a separation <50 au. These numbers yielded planet occurrence rates of planets per star for binaries with separation a B > 100 au and planets per star for binaries with separation a B < 100 au. The similarity in the planet occurrence rate around single stars and wide primaries implies that wide binary systems should actually host more planets than single-star systems, since they have more potential host stars. We estimated a system-wide planet occurrence rate of 0.3 planets per wide binary system for binaries with separations a B > 100 au. Finally, we found evidence that giant planets in binary systems have a different semimajor-axis distribution than their counterparts in single-star systems. The planets in the single-star sample had a significantly higher occurrence rate outside of 1 au than inside 1 au by nearly 4σ, in line with expectations that giant planets are most common near the snow line. However, the planets in the wide binary systems did not follow this distribution, but rather had equivalent occurrence rates interior and exterior to 1 au. This may point to binary-mediated planet migration acting on our sample, even in binaries wider than 100 au.

We present the discovery of Kepler-88 d (Pd = 1403 +/- 14 days, M sin i(d) = 965 +/- 44 M-circle plus = 3.04 +/- 0.13 M-J, e(d) = 0.41 +/- 0.03) based on six years of radial velocity (RV) follow-up from the W. M. Keck Observatory High Resolution Echelle Spectrometer spectrograph. Kepler-88 has two previously identified planets. Kepler-88 b (KOI-142.01) transits in the NASA Kepler photometry and has very large transit timing variations (TTVs). Nesvorny et al. performed a dynamical analysis of the TTVs to uniquely identify the orbital period and mass of the perturbing planet (Kepler-88 c), which was later was confirmed with RVs from the Observatoire de Haute-Provence (OHP). To fully explore the architecture of this system, we performed photodynamical modeling on the Kepler photometry combined with the RVs from Keck and OHP and stellar parameters from spectroscopy and Gaia. Planet d is not detectable in the photometry, and long-baseline RVs are needed to ascertain its presence. A photodynamical model simultaneously optimized to fit the RVs and Kepler photometry yields the most precise planet masses and orbital properties yet for b and c: P-b = 10.91647 +/- 0.00014 days, M-b = 9.5 +/- 1.2 M-circle plus, P-c = 22.2649 +/- 0.0007 days, and M-c = 214.0 +/- 5.3 M-circle plus. The photodynamical solution also finds that planets b and c have low eccentricites and low mutual inclination, are apsidally anti-aligned, and have conjunctions on the same hemisphere of the star. Continued RV follow-up of systems with small planets will improve our understanding of the link between inner planetary system architectures and giant planets.

K2-291 (EPIC 247418783) is a solar-type star with a radius of R_star = 0.899 $\pm$ 0.034 R_sun and mass of M_star=0.934 $\pm$ 0.038 M_sun. From K2 C13 data, we found one super-Earth planet (R_p = 1.589+0.095-0.072 R_Earth) transiting this star on a short period orbit (P = 2.225177 +6.6e-5 -6.8e-5 days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using HIRES on Keck I and HARPS-N on Telescopio Nazionale Galileo, we constrained the mass of EPIC 247418783b to M_p = 6.49 $\pm$ 1.16 M_Earth. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process in order to more accurately model the effect of stellar noise on our data; the addition of the Gaussian process also improved the precision of this mass measurement. With a bulk density of 8.84+2.50-2.03 g cm-3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photo-evaporation during the first billion years of the star's lifetime.

We present the most sensitive direct imaging and radial velocity (RV) exploration of ϵ Eridani to date. ϵ Eridani is an adolescent planetary system, reminiscent of the early solar system. It is surrounded by a prominent and complex debris disk that is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 yr ago, but has met with scrutiny due to possible confusion with stellar noise. We confirm the planet with a new compilation and analysis of precise RV data spanning 30 yr, and combine it with upper limits from our direct imaging search, the most sensitive ever performed. The deep images were taken in the Ms band (4.7 μm) with the vortex coronagraph recently installed in W.M. Keck Observatory’s infrared camera NIRC2, which opens a sensitive window for planet searches around nearby adolescent systems. The RV data and direct imaging upper limit maps were combined in an innovative joint Bayesian analysis, providing new constraints on the mass and orbital parameters of the elusive planet. ϵ Eridani b has a mass of MJup and is orbiting ϵ Eridani at about 3.48 ± 0.02 au with a period of 7.37 ± 0.07 yr. The eccentricity of ϵ Eridani b’s orbit is , an order of magnitude smaller than early estimates and consistent with a circular orbit. We discuss our findings from the standpoint of planet–disk interactions and prospects for future detection and characterization with the James Webb Space Telescope.

K2-291 is a solar-type star with a radius of R* = 0.899 ± 0.034 R⊙ and mass of M* = 0.934 ± 0.038 M⊙. From the K2 C13 data, we found one super-Earth planet (Rp = R⊕) transiting this star on a short period orbit (P = days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using High Resolution Echelle Spectrometer on Keck I and High Accuracy Radial velocity Planet Searcher in the northern hemisphere on Telescopio Nazionale Galileo, we constrained the mass of K2-291 b to Mp = 6.49 ± 1.16 M⊕. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process (GP) in order to more accurately model the effect of stellar noise on our data; the addition of the GP also improved the precision of this mass measurement. With a bulk density of ρ = 8.84 g cm−3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photoevaporation during the first billion years of the star’s lifetime.

We report improved masses, radii, and densities for four planets in two bright M-dwarf systems, K2-3 and GJ3470, derived from a combination of new radial velocity and transit observations. Supplementing K2 photometry with follow-up Spitzer transit observations refined the transit ephemerides of K2-3 b, c, and d by over a factor of 10. We analyze ground-based photometry from the Evryscope and Fairborn Observatory to determine the characteristic stellar activity timescales for our Gaussian Process fit, including the stellar rotation period and activity region decay timescale. The stellar rotation signals for both stars are evident in the radial velocity data and is included in our fit using a Gaussian process trained on the photometry. We find the masses of K2-3 b, K2-3 c, and GJ3470 b to be 6.48(-0.93)(+0.99), 2.14(-1.04)(+1.08), and 12.58(-1.28)(+1.31) M-circle plus, respectively. K2-3 d was not significantly detected and has a 3 sigma upper limit of 2.80 M-circle plus. These two systems are training cases for future TESS systems; due to the low planet densities (rho < 3.7 g cm(-3)) and bright host stars (K < 9 mag), they are among the best candidates for transmission spectroscopy in order to characterize the atmospheric compositions of small planets.

We report on the discovery of a white dwarf companion to the nearby late G dwarf star, HD 159062. The companion is detected in 14 years of precise radial velocity (RV) data, and in high-resolution imaging observations. RVs of HD 159062 from 2003 to 2018 reveal an acceleration of -13.3 +/- 0.12 m s(-1 )yr(-1), indicating that it hosts a companion with a long-period orbit. Subsequent imaging observations with the ShaneAO system on the Lick Observatory 3 m Shane telescope, the PHARO AO system on the Palomar Observatory 5 m telescope, and the NIRC2 AO system at the Keck II 10 m telescope reveal a faint companion 2 ''.7 from the primary star. We performed relative photometry, finding Delta J = 10.09 +/- 0.38 mag, Delta K-S = 10.06 +/- 0.22 mag, and Delta L' = 9.67 +/- 0.08 mag for the companion from these observations. Analysis of the radial velocities, astrometry, and photometry reveals that the combined data set can only be reconciled for the scenario where HD 159062 B is a white dwarf. A full Bayesian analysis of the RV and imaging data to obtain the cooling age, mass, and orbital parameters of the white dwarf indicates that the companion is an old M-B = 0.65(-0.04)(+0.12)M(circle dot) white dwarf with an orbital period of P = 250(-7)(6)(+130) yr, and a cooling age of tau = 8.2(-0.5)(+0.3) Gyr.