This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

The study of planet occurrence as a function of stellar mass is important for a better understanding of planet formation. Estimating stellar mass, especially in the red giant regime, is difficult. In particular, stellar masses of a sample of evolved planet-hosting stars based on spectroscopy and grid-based modelling have been put to question over the past decade with claims they were overestimated. Although efforts have been made in the past to reconcile this dispute using asteroseismology, results were inconclusive. In an attempt to resolve this controversy, we study four more evolved planet-hosting stars in this paper using asteroseismology, and we revisit previous results to make an informed study of the whole ensemble in a self-consistent way. For the four new stars, we measure their masses by locating their characteristic oscillation frequency, $\mathrm{\nu}_{\mathrm{max}}$, from their radial velocity time series observed by SONG. For two stars, we are also able to measure the large frequency separation, $\mathrm{\Delta\nu}$, helped by extended SONG single-site and dual-site observations and new TESS observations. We establish the robustness of the $\mathrm{\nu}_{\mathrm{max}}$-only-based results by determining the stellar mass from $\mathrm{\Delta\nu}$, and from both $\mathrm{\Delta\nu}$ and $\mathrm{\nu}_{\mathrm{max}}$. We then compare the seismic masses of the full ensemble of 16 stars with the spectroscopic masses from three different literature sources. We find an offset between the seismic and spectroscopic mass scales that is mass-dependent, suggesting that the previously claimed overestimation of spectroscopic masses only affects stars more massive than about 1.6 M$_\mathrm{\odot}$.

We investigate the properties of 1306 red giant stars with high photospheric abundances of lithium observed by the GALAH, K2-HERMES and TESS-HERMES surveys, and discuss them in the context of proposed mechanisms for lithium enrichment in giant stars. We confirm that Li-rich giants are rare, making up only 1.1 per cent of our giant star sample. We use GALAH+ DR3 stellar parameters and a Bayesian isochrone analysis to divide the sample into first-ascent red giant branch and red clump stars, and confirm these classifications using asteroseismic data from K2. We find that red clump stars are more than three times as likely to be lithium-rich as red giant branch stars, and the occurrence rate of lithium richness with metallicity is quite different between the two populations. The probability for a rapidly rotating giant to be lithium-rich is distinctly higher as compared to a non rapidly rotating giant, independent of its evolutionary state. There is also a clear correspondence on the red clump between rapid rotation and extreme levels of lithium enrichment. There are almost no stars on the secondary red clump with lithium abundances above the primordial level, indicating a mass dependence in the red clump lithium enrichment mechanism(s). The complex distribution of lithium-rich giants across evolutionary phase, metallicity, rotation rate and mass implies multiple independent mechanisms for producing lithium enrichment.

We report the Transiting Exoplanet Survey Satellite discovery of three small planets transiting one of the nearest and brightest M dwarf hosts to date, TOI-270 (TIC 259377017; K-mag 8.3; 22.5 parsec). The M3V-type star is transited by the super-Earth-sized TOI-270 b (1.247+0.089-0.083 R_earth) and the sub-Neptune-sized exoplanets TOI-270 c (2.42+-0.13 R_earth) and TOI-270 d (2.13+-0.12 R_earth). The planets orbit close to a mean-motion resonant chain, with periods (3.36, 5.66, and 11.38 days) near ratios of small integers (5:3 and 2:1). TOI-270 is a prime target for future studies since: 1) its near-resonance allows detecting transit timing variations for precise mass measurements and dynamical studies; 2) its brightness enables independent radial velocity mass measurements; 3) the outer planets are ideal for atmospheric characterisation via transmission spectroscopy; and 4) the quiet star enables future searches for habitable zone planets. Altogether, very few systems with small, temperate exoplanets are as suitable for such complementary and detailed characterisation as TOI-270.

All-sky photometric time-series missions have allowed for the monitoring of thousands of young ($t_{\rm age} < 800$Myr) to understand the evolution of stellar activity. Here we developed a convolutional neural network (CNN), $\texttt{stella}$, specifically trained to find flares in $\textit{Transiting Exoplanet Survey Satellite}$ ($\textit{TESS}$) short-cadence data. We applied the network to 3200 young stars to evaluate flare rates as a function of age and spectral type. The CNN takes a few seconds to identify flares on a single light curve. We also measured rotation periods for 1500 of our targets and find that flares of all amplitudes are present across all spot phases, suggesting high spot coverage across the entire surface. Additionally, flare rates and amplitudes decrease for stars $t_{\rm age} > 50$Myr across all temperatures $T_{\rm eff} \geq 4000$K, while stars from $2300 \leq T_{\rm eff} < 4000$K show no evolution across 800 Myr. Stars of $T_{\rm eff} \leq 4000$K also show higher flare rates and amplitudes across all ages. We investigate the effects of high flare rates on photoevaporative atmospheric mass loss for young planets. In the presence of flares, planets lose 4-7% more atmosphere over the first 1 Gyr. $\texttt{stella}$ is an open-source Python tool-kit hosted on GitHub and PyPI.

Magnetic activity of the Sun and other stars causes their brightness to vary. We investigate how typical the Sun's variability is compared to other solar-like stars, i.e. those with near-solar effective temperatures and rotation periods. By combining four years of photometric observations from the Kepler space telescope with astrometric data from the Gaia spacecraft, we measure photometric variabilities of 369 solar-like stars. Most of the solar-like stars with well-determined rotation periods show higher variability than the Sun and are therefore considerably more active. These stars appear nearly identical to the Sun, except for their higher variability. Their existence raises the question of whether the Sun can also experience epochs of such high variability.

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

The abundance of short-period planetary systems with high orbital obliquities relative to the spin of their host stars is often taken as evidence that scattering processes play important roles in the formation and evolution of these systems. More recent studies have suggested that wide binary companions can tilt protoplanetary disks, inducing a high stellar obliquity that form through smooth processes like disk migration. DS Tuc Ab, a transiting planet with an 8.138 day period in the 40 Myr Tucana-Horologium association, likely orbits in the same plane as its now-dissipated protoplanetary disk, enabling us to test these theories of disk physics. Here, we report on Rossiter-McLaughlin observations of one transit of DS Tuc Ab with the Planet Finder Spectrograph on the Magellan Clay Telescope at Las Campanas Observatory. We confirm the previously detected planet by modeling the planet transit and stellar activity signals simultaneously. We test multiple models to describe the stellar activity-induced radial velocity variations over the night of the transit, finding the obliquity to be low: $\lambda = 12 \pm 13$ degrees, suggesting that this planet likely formed through smooth disk processes and its protoplanetary disk was not significantly torqued by DS Tuc B. The specific stellar activity model chosen affects the results at the $\approx 5$ degree level. This is the youngest planet to be observed using this technique; we provide a discussion on best practices to accurately measure the observed signal of similar young planets.

In recent years, dedicated extreme-precision radial velocity (EPRV) spectrographs have produced vast quantities of high-resolution, high-signal-to-noise time-series spectra for bright stars. These data contain valuable information for the dual purposes of planet detection via the measured RVs and stellar characterization via the co-added spectra. However, considerable data analysis challenges exist in extracting these data products from the observed spectra at the highest possible precision, including the issue of poorly-characterized telluric absorption features and the common use of an assumed stellar spectral template. In both of these examples, precision-limiting reliance on external information can be sidestepped using the data directly. Here we propose a data-driven method to simultaneously extract precise RVs and infer the underlying stellar and telluric spectra using a linear model (in the log of flux). The model employs a convex objective and convex regularization to keep the optimization of the spectral components fast. We implement this method in wobble, an open-source python package which uses TensorFlow in one of its first non-neural-network applications to astronomical data. In this work, we demonstrate the performance of wobble on archival HARPS spectra. We recover the canonical exoplanet 51 Pegasi b, detect the secular RV evolution of the M dwarf Barnard's Star, and retrieve the Rossiter-McLaughlin effect for the Hot Jupiter HD 189733b. The method additionally produces extremely high-S/N composite stellar spectra and detailed time-variable telluric spectra, which we also present here.

K2-146 is a mid-M dwarf ($M_\star = 0.331 \pm 0.009 M_\odot$; $R_\star = 0.330 \pm 0.010 R_\odot$), observed in Campaigns 5, 16, and 18 of the K2 mission. In Campaign 5 data, a single planet was discovered with an orbital period of $2.6$~days and large transit timing variations due to an unknown perturber. Here we analyze data from Campaigns 16 and 18, detecting the transits of a second planet, c, with an orbital period of $4.0$~days, librating in a 3:2 resonance with planet b. Large, anti-correlated timing variations of both planets exist due to their resonant perturbations. The planets have a mutual inclination of $2.40^\circ\pm0.25^\circ$, which torqued planet c more closely into our line-of-sight. Planet c was grazing in Campaign 5 and thus missed in previous searches; in Campaigns 16 and 18 it is fully transiting, and its transit depth is three times larger. We improve the stellar properties using data from Gaia DR2, and using dynamical fits find that both planets are sub-Neptunes: their masses are $5.77\pm0.18$ and $7.50\pm0.23 M_{\oplus}$ and their radii are $2.04\pm0.06$ and $2.19\pm0.07$ R$_\oplus$, respectively. These mass constraints set the precision record for small exoplanets (a few gas giants have comparable relative precision). These planets lie in the photoevaporation valley when viewed in Radius-Period space, but due to the low-luminosity M-dwarf host star, they lie among the atmosphere-bearing planets when viewed in Radius-Irradiation space. This, along with their densities being 60%-80% that of Earth, suggests that they may both have retained a substantial gaseous envelope.

During its two year prime mission the Transiting Exoplanet Survey Satellite (TESS) will perform a time-series photometric survey covering over 80% of the sky. This survey comprises observations of 26 24 x 96 degree sectors that are each monitored continuously for approximately 27 days. The main goal of TESS is to find transiting planets around 200,000 pre-selected stars for which fixed aperture photometry is recorded every two minutes. However, TESS is also recording and delivering Full-Frame Images (FFIs) of each detector at a 30 minute cadence. We have created an open-source tool, eleanor, to produce light curves for objects in the TESS FFIs. Here, we describe the methods used in eleanor to produce light curves that are optimized for planet searches. The tool performs background subtraction, aperture and PSF photometry, decorrelation of instrument systematics, and cotrending using principal component analysis. We recover known transiting exoplanets in the FFIs to validate the pipeline and perform a limited search for new planet candidates in Sector 1. Our tests indicate that eleanor produces light curves with significantly less scatter than other tools that have been used in the literature. Cadence-stacked images, and raw and detrended eleanor light curves for each analyzed star will be hosted on MAST, with planet candidates on ExoFOP-TESS as Community TESS Objects of Interest (CTOIs). This work confirms the promise that the TESS FFIs will enable the detection of thousands of new exoplanets and a broad range of time domain astrophysics.

Exoplanet surveys of evolved stars have provided increasing evidence that the formation of giant planets depends not only on stellar metallicity ([Fe/H]), but also the mass ($M_\star$). However, measuring accurate masses for subgiants and giants is far more challenging than it is for their main-sequence counterparts, which has led to recent concerns regarding the veracity of the correlation between stellar mass and planet occurrence. In order to address these concerns we use HIRES spectra to perform a spectroscopic analysis on an sample of 245 subgiants and derive new atmospheric and physical parameters. We also calculate the space velocities of this sample in a homogeneous manner for the first time. When reddening corrections are considered in the calculations of stellar masses and a -0.12 M$_{\odot}$ offset is applied to the results, the masses of the subgiants are consistent with their space velocity distributions, contrary to claims in the literature. Similarly, our measurements of their rotational velocities provide additional confirmation that the masses of subgiants with $M_\star \geq 1.6$ M$_{\odot}$ (the "Retired A Stars") have not been overestimated in previous analyses. Using these new results for our sample of evolved stars, together with an updated sample of FGKM dwarfs, we confirm that giant planet occurrence increases with both stellar mass and metallicity up to 2.0 M$_{\odot}$. We show that the probability of formation of a giant planet is approximately a one-to-one function of the total amount of metals in the protoplanetary disk $M_\star 10^{[Fe/H]}$. This correlation provides additional support for the core accretion mechanism of planet formation.

We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system.

We examine the ability of the Transiting Exoplanet Survey Satellite (TESS) to detect and improve our understanding of planetary systems in the Kepler field. By modeling the expected transits of all confirmed and candidate planets detected by Kepler as expected to be observed by TESS, we find that TESS has a greater than 50% chance of detecting 260 of these planets at the 3 sigma level in one sector of observations and an additional 120 planets in two sectors. Most of these are large planets in short orbits around their host stars, although a small number of rocky planets are expected to be recovered. Most of these systems have only one known transiting planet; in only ~5 percent of known multiply-transiting systems do we anticipate more than one planet to be recovered. When these planets are recovered, we expect TESS to be a powerful tool to characterizing transit timing variations. Using Kepler-88 (KOI-142) as an example, we show that TESS will improve measurements of planet-star mass ratios and orbital parameters, and significantly reduce the transit timing uncertainty in future years. Since TESS will be most sensitive to hot Jupiters, we research whether TESS will be able to detect tidal orbital decay in these systems. We find two confirmed planetary systems (Kepler-2 b and Kepler-13 b) and five candidate systems that will be good candidates to detect tidal decay.

A Neptune-sized exomoon candidate was recently announced by Teachey & Kipping, orbiting a 287 day gas giant in the Kepler-1625 system. However, the system is poorly characterized and needs more observations to be confirmed, with the next potential transit in 2019 May. In this Letter, we aid observational follow up by analyzing the transit signature of exomoons. We derive a simple analytic equation for the transit probability and use it to demonstrate how exomoons may frequently avoid transit if their orbit is larger than the stellar radius and sufficiently misaligned. The nominal orbit for the moon in Kepler-1625 has both of these characteristics, and we calculate that it may only transit roughly 40% of the time. This means that approximately six non-transits would be required to rule out the moon's existence at 95% confidence. When an exomoon's impact parameter is displaced off the star, the planet's impact parameter is displaced the other way, so larger planet transit durations are typically positively correlated with missed exomoon transits. On the other hand, strong correlations do not exist between missed exomoon transits and transit timing variations of the planet. We also show that nodal precession does not change an exomoon's transit probability and that it can break a prograde-retrograde degeneracy.

The Wide Field Infrared Survey Telescope (WFIRST) will monitor $\sim 2$ deg$^2$ toward the Galactic bulge in a wide ($\sim 1-2~\mu$m) W149 filter at 15-minute cadence with exposure times of $\sim$50s for 6 seasons of 72 days each, for a total $\sim$41,000 exposures taken over $\sim$432 days, spread over the 5-year prime mission. This will be one of the deepest exposures of the sky ever taken, reaching a photon-noise photometric precision of 0.01 mag per exposure and collecting a total of $\sim 10^9$ photons over the course of the survey for a W149$_{\rm AB}\sim 21$ star. Of order $4 \times 10^7$ stars will be monitored with W149$_{\rm AB}$<21, and 10$^8$ stars with W145$_{\rm AB}$<23. The WFIRST microlensing survey will detect $\sim$54,000 microlensing events, of which roughly 1% ($\sim$500) will be due to isolated black holes, and $\sim$3% ($\sim$1600) will be due to isolated neutron stars. It will be sensitive to (effectively) isolated compact objects with masses as low as the mass of Pluto, thereby enabling a measurement of the compact object mass function over 10 orders of magnitude. Assuming photon-noise limited precision, it will detect $\sim 10^5$ transiting planets with sizes as small as $\sim 2~R_\oplus$, perform asteroseismology of $\sim 10^6$ giant stars, measure the proper motions to $\sim 0.3\%$ and parallaxes to $\sim 10\%$ for the $\sim 6 \times 10^6$ disk and bulge stars in the survey area, and directly detect $\sim 5 \times 10^3$ Trans-Neptunian objects (TNOs) with diameters down to $\sim 10$ km, as well as detect $\sim 10^3$ occulations of stars by TNOs during the survey. All of this science will completely serendipitous, i.e., it will not require modifications of the WFIRST optimal microlensing survey design. Allowing for some minor deviation from the optimal design, such as monitoring the Galactic center, would enable an even broader range of transformational science.

By combining ground-based spectrographic observations of variability in the chromospheric emission from Sun-like stars with the variability seen in their eigenmode frequencies, it is possible to relate the changes observed at the surfaces of these stars to the changes taking place in the interior. By further comparing this variability to changes in the relative flux from the stars, one can obtain an expression for how these activity indicators relate to the energy output from the stars. Such studies become very pertinent when the variability can be related to stellar cycles as they can then be used to improve our understanding of the solar cycle and its effect on the energy output from the Sun. Here we present observations of chromospheric emission in 20 Sun-like stars obtained over the course of the nominal 4-year Kepler mission. Even though 4 years is too short to detect stellar equivalents of the 11-year solar cycle, observations from the Kepler mission can still be used to analyse the variability of the different activity indicators thereby obtaining information of the physical mechanism generating the variability. The analysis reveals no strong correlation between the different activity indicators, except in very few cases. We suggest that this is due to the sparse sampling of our ground-based observations on the one hand and that we are likely not tracing cyclic variability on the other hand. We also discuss how to improve the situation.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

We observed eclipses of the transiting brown dwarf CWW 89Ab at 3.6um and 4.5um using Spitzer/IRAC. The CWW 89 binary system is a member of the $3.0\pm0.25$ Gyr-old open cluster Ruprecht 147, and is composed of a Sun-like primary and an early M-dwarf secondary separated by a projected distance of 25 AU. CWW 89Ab has a radius of $0.937\pm0.042$ RJ and a mass of $36.5\pm0.1$ MJ, and is on a 5.3 day orbit about CWW 89A with a non-zero eccentricity of $e=0.19$ (Curtis et al. 2016). We strongly detect the eclipses of CWW 89Ab in both Spitzer channels as $\delta_{3.6}=1147\pm213$ ppm and $\delta_{4.5}=1097\pm225$ ppm after correcting for the dilution from CWW 89B. After accounting for the irradiation that CWW 89Ab receives from its host star, these measurements imply that the brown dwarf has an internal luminosity of $\log(\mathrm{L_{bol}/\mathrm{L}_\odot})=-4.19\pm0.14$. This is 16 times, or $9.3\,\sigma$, higher than model predictions given the known mass, radius, and age of CWW 89Ab. As we discuss, this over-luminosity is not explainable by an inaccurate age determination, additional stellar heating, nor tidal heating. Instead, we suggest that the anomalous luminosity of CWW 89Ab is caused by a dayside temperature inversion -- though a significant error in the evolutionary models is also a possibility. Importantly, a temperature inversion would require a super-stellar C/O ratio in CWW 89Ab's atmosphere. If this is indeed the case, it implies that CWW 89Ab is a 36.5 MJ object that formed via core accretion processes. Finally, we use our measurement of CWW 89Ab's orbital eccentricity, improved via these observations, to constrain the tidal quality factors of the brown dwarf and the host star CWW 89A to be $Q_{BD}>10^{4.15}$ and $Q_*>10^{9}$, respectively.

We present moderate resolution near-infrared spectra in $H, J$ and $K$ band of M dwarf hosts to candidate transiting exoplanets discovered by NASA's K2 mission. We employ known empirical relationships between spectral features and physical stellar properties to measure the effective temperature, radius, metallicity, and luminosity of our sample. Out of an initial sample of 56 late-type stars in K2, we identify 35 objects as M dwarfs. For that sub-sample, we derive temperatures ranging from $2,870$ to $4,187$ K, radii of $0.09-0.83$ $R_{\odot}$, luminosities of $-2.67<log L/L_{\odot}<-0.67$ and [Fe/H] metallicities between $-0.49$ and $0.83$ dex. We then employ the stellar properties derived from spectra, in tandem with the K2 lightcurves, to characterize their planets. We report 33 exoplanet candidates with orbital periods ranging from 0.19 to 21.16 days, and median radii and equilibrium temperatures of 2.3 $R_{\oplus}$ and 986 K, respectively. Using planet mass-radius relationships from the literature, we identify 7 exoplanets as potentially rocky, although we conclude that probably none reside in the habitable zone of their parent stars.

Data from the TESS mission will be used to discover hundreds of small planets amenable to radial velocity (RV) followup. Often, RV observations are obtained until a particular fractional precision on the inferred mass is achieved. I show that when the decision to stop collecting RV observations depends on the mass inferred from already-collected data, this will bias mass measurements upward, particularly when the fixed length of RV observing campaigns is considered. I suggest that observing teams should determine their observing strategy for each star before any data are collected, and all stopping criteria should not directly depend on the inferred mass. Observing teams should explicitly publish both their criteria for observing targets and deciding to end their observations, as well as their mass non-detections to avoid introducing biases into the masses---and thus inferences on the densities, compositions, and atmospheres---of transiting planets.

The enigmatic star KIC 8462852, informally known as "Boyajian's Star", has exhibited unexplained variability from both short timescale (days) dimming events, and years-long fading in the Kepler mission. No single physical mechanism has successfully explained these observations to date. Here we investigate the ultraviolet variability of KIC 8462852 on a range of timescales using data from the GALEX mission that occurred contemporaneously with the Kepler mission. The wide wavelength baseline between the Kepler and GALEX data provides a unique constraint on the nature of the variability. Using 1600 seconds of photon-counting data from four GALEX visits spread over 70 days in 2011, we find no coherent NUV variability in the system on 10-100 second or months timescales. Comparing the integrated flux from these 2011 visits to the 2012 NUV flux published in the GALEX-CAUSE Kepler survey, we find a 3% decrease in brightness for KIC 8462852. We find this level of variability is significant, but not necessarily unusual for stars of similar spectral type in the GALEX data. This decrease coincides with the secular optical fading reported by Montet & Simon (2016). We find the multi-wavelength variability is somewhat inconsistent with typical interstellar dust absorption, but instead favors a R$_V$ = 5.0 $\pm$ 0.9 reddening law potentially from circumstellar dust.

We present ~800 days of photometric monitoring of Boyajian's Star (KIC 8462852) from the All-Sky Automated Survey for Supernovae (ASAS-SN) and ~4000 days of monitoring from the All Sky Automated Survey (ASAS). We show that from 2015 to the present the brightness of Boyajian's Star has steadily decreased at a rate of 6.3 +/- 1.4 mmag yr^-1, such that the star is now 1.5% fainter than it was in February 2015. Moreover, the longer time baseline afforded by ASAS suggests that Boyajian's Star has also undergone two brightening episodes in the past 11 years, rather than only exhibiting a monotonic decline. We analyze a sample of ~1000 comparison stars of similar brightness located in the same ASAS-SN field and demonstrate that the recent fading is significant at >99.4% confidence. The 2015-2017 dimming rate is consistent with that measured with Kepler data for the time period from 2009 to 2013. This long-term variability is difficult to explain with any of the physical models for the star's behavior proposed to date.

Photometry from the Kepler mission is optimized to detect small, short duration signals like planet transits at the expense of long-term trends. This long-term variability can be recovered in photometry from the Full Frame Images (FFIs), a set of calibration data collected approximately monthly during the Kepler mission. Here, we present f3, an open-source package to perform photometry on the Kepler FFIs in order to detect changes in the brightness of stars in the Kepler field of view over long time baselines. We apply this package to a sample of 4,000 Sun-like stars with measured rotation periods. We find $\approx 10\%$ of these targets have long-term variability in their observed flux. For the majority of targets we identify the luminosity variations are either correlated or anticorrelated with the short-term variability due to starspots on the stellar surface. We find a transition between anticorrelated (starspot-dominated) variability and correlated (facula-dominated) variability between rotation periods of 15 and 25 days, suggesting the transition between the two modes is complete for stars at the age of the Sun. We also identify a sample of stars with apparently complete cycles, as well as a collection of short-period binaries with extreme photometric variation over the Kepler mission.

Observations of Sun-like stars over the last half-century have improved our understanding of how magnetic dynamos, like that responsible for the 11-year solar cycle, change with rotation, mass and age. Here we show for the first time how metallicity can affect a stellar dynamo. Using the most complete set of observations of a stellar cycle ever obtained for a Sun-like star, we show how the solar analog HD 173701 exhibits solar-like differential rotation and a 7.4-year activity cycle. While the duration of the cycle is comparable to that generated by the solar dynamo, the amplitude of the brightness variability is substantially stronger. The only significant difference between HD 173701 and the Sun is its metallicity, which is twice the solar value. Therefore, this provides a unique opportunity to study the effect of the higher metallicity on the dynamo acting in this star and to obtain a comprehensive understanding of the physical mechanisms responsible for the observed photometric variability. The observations can be explained by the higher metallicity of the star, which is predicted to foster a deeper outer convection zone and a higher facular contrast, resulting in stronger variability.

Kepler photometry of the hot Neptune host star HAT-P-11 suggests that its spot latitude distribution is comparable to the Sun's near solar maximum. We search for evidence of an activity cycle in the CaII H & K chromospheric emission $S$-index with archival Keck/HIRES spectra and observations from the echelle spectrograph on the ARC 3.5 m Telescope at APO. The chromospheric emission of HAT-P-11 is consistent with a $\gtrsim 10$ year activity cycle, which plateaued near maximum during the Kepler mission. In the cycle that we observed, the star seemed to spend more time near active maximum than minimum. We compare the $\log R^\prime_{HK}$ normalized chromospheric emission index of HAT-P-11 with other stars. HAT-P-11 has unusually strong chromospheric emission compared to planet-hosting stars of similar effective temperature and rotation period, perhaps due to tides raised by its planet.

We image 104 newly identified low-mass (mostly M-dwarf) pre-main sequence members of nearby young moving groups with Magellan Adaptive Optics (MagAO) and identify 27 binaries with instantaneous projected separation as small as 40 mas. 15 were previously unknown. The total number of multiple systems in this sample including spectroscopic and visual binaries from the literature is 36, giving a raw multiplicity rate of at least $35^{+5}_{-4}\%$ for this population. In the separation range of roughly 1 - 300 AU in which infrared AO imaging is most sensitive, the raw multiplicity rate is at least $24^{+5}_{-4}\%$ for binaries resolved by the MagAO infrared camera (Clio). The M-star sub-sample of 87 stars yields a raw multiplicity of at least $30^{+5}_{-4}\%$ over all separations, $21^{+5}_{-4}\%$ for secondary companions resolved by Clio from 1 to 300 AU ($23^{+5}_{-4}\%$ for all known binaries in this separation range). A combined analysis with binaries discovered by the Search for Associations Containing Young stars shows that multiplicity fraction as a function of mass and age over the range of 0.2 to 1.2 $M_\odot$ and 10 - 200 Myr appears to be linearly flat in both parameters and across YMGs. This suggests that multiplicity rates are largely set by 100 Myr without appreciable evolution thereafter. After bias corrections are applied, the multiplicity fraction of low-mass YMG members ($< 0.6 M_\odot$) is in excess of the field.

Measurements of radial velocity variations from the spectroscopic monitoring of stars and their companions are essential for a broad swath of astrophysics, providing access to the fundamental physical properties that dictate all phases of stellar evolution and facilitating the quantitative study of planetary systems. The conversion of those measurements into both constraints on the orbital architecture and individual component spectra can be a serious challenge, however, especially for extreme flux ratio systems and observations with relatively low sensitivity. Gaussian processes define sampling distributions of flexible, continuous functions that are well-motivated for modeling stellar spectra, enabling proficient search for companion lines in time-series spectra. We introduce a new technique for spectral disentangling, where the posterior distributions of the orbital parameters and intrinsic, rest-frame stellar spectra are explored simultaneously without needing to invoke cross-correlation templates. To demonstrate its potential, this technique is deployed on red-optical time-series spectra of the mid-M dwarf binary LP661-13. We report orbital parameters with improved precision compared to traditional radial velocity analysis and successfully reconstruct the primary and secondary spectra. We discuss potential applications for other stellar and exoplanet radial velocity techniques and extensions to time-variable spectra. The code used in this analysis is freely available as an open source Python package.

The WFIRST microlensing mission will measure precise light curves and relative parallaxes for millions of stars, giving it the potential to characterize short-period transiting planets all along the line of sight and into the galactic bulge. These light curves will enable the detection of more than 100,000 transiting planets whose host stars have measured distances. Although most of these planets cannot be followed up, several thousand hot Jupiters can be confirmed directly by detection of their secondary eclipses in the WFIRST data. Additionally, some systems of small planets may be confirmed by detecting transit timing variations over the duration of the WFIRST microlensing survey. Finally, many more planets may be validated by ruling out potential false positives. The combination of WFIRST transits and microlensing will provide a complete picture of planetary system architectures, from the very shortest periods to unbound planets, as a function of galactocentric distance.

KIC 8462852 is a superficially ordinary main sequence F star for which Kepler detected an unusual series of brief dimming events. We obtain accurate relative photometry of KIC 8462852 from the Kepler full frame images, finding that the brightness of KIC 8462852 monotonically decreased over the four years it was observed by Kepler. Over the first ~1000 days, KIC 8462852 faded approximately linearly at a rate of 0.341 +/- 0.041%/yr, for a total decline of 0.9%. KIC 8462852 then dimmed much more rapidly in the next ~200 days, with its flux dropping by more than 2%. For the final ~200 days of Kepler photometry the magnitude remained approximately constant, although the data are also consistent with the decline rate measured for the first 2.7 yr. Of a sample of 193 nearby comparison stars and 355 stars with similar stellar parameters, none exhibit the rapid decline by 2% or the cumulative fading by 3% of KIC 8462852. Moreover, of these comparison stars, only one changes brightness as quickly as the 0.341%/yr measured for KIC 8462852 during the first three years of the Kepler mission. We examine whether the rapid decline could be caused by a cloud of transiting circumstellar material, finding while such a cloud could evade detection in sub-mm observations, the transit ingress and duration cannot be explained by a simple cloud model. Moreover, this model cannot account for the observed longer-term dimming. No known or proposed stellar phenomena can fully explain all aspects of the observed light curve.

While the vast majority of multiple-planet systems have their orbital angular momentum axes aligned with the spin axis of their host star, Kepler-56 is an exception: its two transiting planets are coplanar yet misaligned by at least 40 degrees with respect to their host star. Additional follow-up observations of Kepler-56 suggest the presence of a massive, non-transiting companion that may help explain this misalignment. We model the transit data along with Keck/HIRES and HARPS-N radial velocity data to update the masses of the two transiting planets and infer the physical properties of the third, non-transiting planet. We employ a Markov Chain Monte Carlo sampler to calculate the best-fitting orbital parameters and their uncertainties for each planet. We find the outer planet has a period of 1002 $\pm$ 5 days and minimum mass of 5.61 $\pm$ 0.38 Jupiter masses. We also place a 95% upper limit of 0.80 m/s/yr on long-term trends caused by additional, more distant companions.

The mass and age of substellar objects are degenerate parameters leaving the evolutionary state of brown dwarfs ambiguous without additional information. Theoretical models are normally used to help distinguish between old, massive brown dwarfs and young, low mass brown dwarfs but these models have yet to be properly calibrated. We have carried out an infrared high-contrast imaging program with the goal of detecting substellar objects as companions to nearby stars to help break degeneracies in inferred physical properties such as mass, age, and composition. Rather than using imaging observations alone, our targets are pre-selected based on the existence of dynamical accelerations informed from years of stellar radial velocity (RV) measurements. In this paper, we present the discovery of a rare benchmark brown dwarf orbiting the nearby ($d=18.69\pm0.19$ pc), solar-type (G9V) star HD 4747 ([Fe/H]=$-0.22\pm0.04$) with a projected separation of only $\rho=11.3\pm0.2$ AU ($\theta \approx$ 0.6"). Precise Doppler measurements taken over 18 years reveal the companion's orbit and allow us to place strong constraints on its mass using dynamics ($m \sin(i) = 55.3\pm1.9M_J$). Relative photometry ($\Delta K_s=9.05\pm0.14$, $M_{K_s}=13.00\pm0.14$, $K_s - L' = 1.34\pm0.46$) indicates that HD 4747 B is most-likely a late-type L-dwarf and, if near the L/T transition, an intriguing source for studying cloud physics, variability, and polarization. We estimate a model-dependent mass of $m=72^{+3}_{-13}M_J$ for an age of $3.3^{+2.3}_{-1.9}$ Gyr based on gyrochronology. Combining astrometric measurements with RV data, we calculate the companion dynamical mass ($m=60.2\pm3.3M_J$) and orbit ($e=0.740\pm0.002$) directly. As a new mass, age, and metallicity benchmark, HD 4747 B will serve as a laboratory for precision astrophysics to test theoretical models that describe the emergent radiation of brown dwarfs.

There is currently no evidence for life on any known exoplanet. Here, we propose a form of "galactic anthropology" to detect not only the existence of life on transiting exoplanets, but also the existence of environmentalism movements. By observing the planet's atmosphere over long time baselines, the destruction and recovery of a hole in an exoplanet's ozone layer may be observable. While not readily detectable for any one system with JWST, by binning together observations of hundreds of systems we can finally determine the occurrence rate of environmental movements on Earthlike planets in the galaxy, a number we term eta-Green-Earth.

There are no field brown dwarf analogs with measured masses, radii, and luminosities, precluding our ability to connect the population of transiting brown dwarfs with measurable masses and radii and field brown dwarfs with measurable luminosities and atmospheric properties. LHS 6343 C, a weakly-irradiated brown dwarf transiting one member of an M+M binary in the Kepler field, provides the first opportunity to probe the atmosphere of a non-inflated brown dwarf with a measured mass and radius. Here, we analyze four Spitzer observations of secondary eclipses of LHS 6343 C behind LHS 6343 A. Jointly fitting the eclipses with a Gaussian process noise model of the instrumental systematics, we measure eclipse depths of 1.06 \pm 0.21 ppt at 3.6 microns and 2.09 \pm 0.08 ppt at 4.5 microns, corresponding to brightness temperatures of 1026 \pm 57 K and 1249 \pm 36 K, respectively. We then apply brown dwarf evolutionary models to infer a bolometric luminosity log(L_star / L_sun) = -5.16 \pm 0.04. Given the known physical properties of the brown dwarf and the two M dwarfs in the LHS 6343 C system, these depths are consistent with models of a 1100 K T dwarf at an age of 5 Gyr and empirical observations of field T5-6 dwarfs with temperatures of 1070 \pm 130 K. We investigate the possibility that the orbit of LHS 6343 C has been altered by the Kozai-Lidov mechanism and propose additional astrometric or Rossiter-McLaughlin measurements of the system to probe the dynamical history of the system.

We measure carbon radio recombination line (RRL) emission at 5.3 GHz toward four HII regions with the Green Bank Telescope (GBT) to determine the magnetic field strength in the photodissociation region (PDR) that surrounds the ionized gas. Roshi (2007) suggests that the non-thermal line widths of carbon RRLs from PDRs are predominantly due to magneto-hydrodynamic (MHD) waves, thus allowing the magnetic field strength to be derived. We model the PDR with a simple geometry and perform the non-LTE radiative transfer of the carbon RRL emission to solve for the PDR physical properties. Using the PDR mass density from these models and the carbon RRL non-thermal line width we estimate total magnetic field strengths of B ~ 100-300 micro Gauss in W3 and NGC6334A. Our results for W49 and NGC6334D are less well constrained with total magnetic field strengths between B ~ 200-1000 micro Gauss. HI and OH Zeeman measurements of the line-of-sight magnetic field strength (B_los), taken from the literature, are between a factor of ~0.5-1 of the lower bound of our carbon RRL magnetic field strength estimates. Since |B_los| <= B, our results are consistent with the magnetic origin of the non-thermal component of carbon RRL widths.

We combine new high resolution imaging and spectroscopy from Keck/NIRC2, Discovery Channel Telescope/DSSI, and Keck/HIRES with published astrometry and radial velocities to measure individual masses and orbital elements of the GJ 3305 AB system, a young (~20 Myr) M+M binary (unresolved spectral type M0) member of the beta Pictoris moving group comoving with the imaged exoplanet host 51 Eri. We measure a total system mass of 1.11 \pm 0.04 M_sun, a period of 29.03 \pm 0.50$ yr, a semimajor axis of 9.78 \pm 0.14 AU, and an eccentricity of 0.19 \pm 0.02. The primary component has a dynamical mass of 0.67 \pm 0.05 M_sun and the secondary has a mass of 0.44 \pm 0.05 M_sun. The recently updated BHAC15 models are consistent with the masses of both stars to within 1.5 sigma. Given the observed masses the models predict an age of the GJ 3305 AB system of 37 \pm 9 Myr. Based on the the observed system architecture and our dynamical mass measurement, it is unlikely that the orbit of 51 Eri b has been significantly altered by the Kozai-Lidov mechanism.

The T Tauri star PTFO 8-8695 exhibits periodic fading events that have been interpreted as the transits of a giant planet on a precessing orbit. Here we present three tests of the planet hypothesis. First, we sought evidence for the secular changes in light-curve morphology that are predicted to be a consequence of orbital precession. We observed 28 fading events spread over several years, and did not see the expected changes. Instead we found that the fading events are not strictly periodic. Second, we attempted to detect the planet's radiation, based on infrared observations spanning the predicted times of occultations. We ruled out a signal of the expected amplitude. Third, we attempted to detect the Rossiter-McLaughlin effect by performing high-resolution spectroscopy throughout a fading event. No effect was seen at the expected level, ruling out most (but not all) possible orientations for the hypothetical planetary orbit. Our spectroscopy also revealed strong, time-variable, high-velocity H{\alpha} and Ca H & K emission features. All these observations cast doubt on the planetary hypothesis, and suggest instead that the fading events represent starspots, eclipses by circumstellar dust, or occultations of an accretion hotspot.

The extended Kepler mission, K2, is now providing photometry of new fields every three months in a search for transiting planets. In a recent study, Foreman-Mackey and collaborators presented a list of 36 planet candidates orbiting 31 stars in K2 Campaign 1. In this contribution, we present stellar and planetary properties for all systems. We combine ground-based seeing-limited survey data and adaptive optics imaging with an automated transit analysis scheme to validate 21 candidates as planets, 17 for the first time, and identify 6 candidates as likely false positives. Of particular interest is K2-18 (EPIC 201912552), a bright (K=8.9) M2.8 dwarf hosting a 2.23 \pm 0.25 R_Earth planet with T_eq = 272 \pm 15 K and an orbital period of 33 days. We also present two new open-source software packages which enable this analysis. The first, isochrones, is a flexible tool for fitting theoretical stellar models to observational data to determine stellar properties using a nested sampling scheme to capture the multimodal nature of the posterior distributions of the physical parameters of stars that may plausibly be evolved. The second is vespa, a new general-purpose procedure to calculate false positive probabilities and statistically validate transiting exoplanets.

Photometry of stars from the K2 extension of NASA's Kepler mission is afflicted by systematic effects caused by small (few-pixel) drifts in the telescope pointing and other spacecraft issues. We present a method for searching K2 light curves for evidence of exoplanets by simultaneously fitting for these systematics and the transit signals of interest. This method is more computationally expensive than standard search algorithms but we demonstrate that it can be efficiently implemented and used to discover transit signals. We apply this method to the full Campaign 1 dataset and report a list of 36 planet candidates transiting 31 stars, along with an analysis of the pipeline performance and detection efficiency based on artificial signal injections and recoveries. For all planet candidates, we present posterior distributions on the properties of each system based strictly on the transit observables.

Kepler-296 is a binary star system with two M-dwarf components separated by 0.2 arcsec. Five transiting planets have been confirmed to be associated with the Kepler-296 system; given the evidence to date, however, the planets could in principle orbit either star. This ambiguity has made it difficult to constrain both the orbital and physical properties of the planets. Using both statistical and analytical arguments, this paper shows that all five planets are highly likely to orbit the primary star in this system. We performed a Markov-Chain Monte Carlo simulation using a five transiting planet model, leaving the stellar density and dilution with uniform priors. Using importance sampling, we compared the model probabilities under the priors of the planets orbiting either the brighter or the fainter component of the binary. A model where the planets orbit the brighter component, Kepler-296A, is strongly preferred by the data. Combined with our assertion that all five planets orbit the same star, the two outer planets in the system, Kepler-296 Ae and Kepler-296 Af, have radii of 1.53 +/- 0.26 and 1.80 +/- 0.31 R_earth, respectively, and receive incident stellar fluxes of 1.40 +/- 0.23 and 0.62 +/- 0.10 times the incident flux the Earth receives from the Sun. This level of irradiation places both planets within or close to the circumstellar habitable zone of their parent star.

We present optical and near-infrared adaptive optics (AO) imaging and spectroscopy of 13 ultracool (>M6) companions to late-type stars (K7-M4.5), most of which have recently been identified as candidate members of nearby young moving groups (YMGs; 8-120 Myr) in the literature. The inferred masses of the companions (~10-100 Mjup) are highly sensitive to the ages of the primary stars so we critically examine the kinematic and spectroscopic properties of each system to distinguish bona fide YMG members from old field interlopers. 2MASS J02155892-0929121 C is a new M7 substellar companion (40-60 Mjup) with clear spectroscopic signs of low gravity and hence youth. The primary, possibly a member of the ~40 Myr Tuc-Hor moving group, is visually resolved into three components, making it a young low-mass quadruple system in a compact (<100 AU) configuration. In addition, Li 1 $\lambda$6708 absorption in the intermediate-gravity M7.5 companion 2MASS J15594729+4403595 B provides unambiguous evidence that it is young (<200 Myr) and resides below the hydrogen burning limit. Three new close-separation (<1") companions (2MASS J06475229-2523304 B, PYC J11519+0731 B, and GJ 4378 Ab) orbit stars previously reported as candidate YMG members, but instead are likely old (>1 Gyr) tidally-locked spectroscopic binaries without convincing kinematic associations with any known moving group. The high rate of false positives in the form of old active stars with YMG-like kinematics underscores the importance of radial velocity and parallax measurements to validate candidate young stars identified via proper motion and activity selection alone. Finally, we spectroscopically confirm the cool temperature and substellar nature of HD 23514 B, a recently discovered M8 benchmark brown dwarf orbiting the dustiest-known member of the Pleiades. [Abridged]

The coolest dwarf stars targeted by the Kepler Mission constitute a relatively small but scientifically valuable subset of the Kepler target stars, and provide a high-fidelity and nearby sample of transiting planetary systems. Using archival Kepler data spanning the entire primary mission we perform a uniform analysis to extract, confirm and characterize the transit signals discovered by the Kepler pipeline toward M-type dwarf stars. We recover all but two of the signals reported in a recent listing from the Exoplanet Archive resulting in 165 planet candidates associated with a sample of 106 low-mass stars. We fitted the observed light curves to transit models using Markov Chain Monte Carlo and we have made the posterior samples publicly available to facilitate further studies. We fitted empirical transit times to individual transit signals with significantly non-linear ephemerides for accurate recovery of transit parameters and measuring precise transit timing variations. We also provide the physical parameters for the stellar sample, including new measurements of stellar rotation, allowing the conversion of transit parameters into planet radii and orbital parameters.

We present an updated analysis of LHS 6343, a triple system in the Kepler field which consists of a brown dwarf transiting one member of a widely-separated M+M binary system. By analyzing the full Kepler dataset and 34 Keck/HIRES radial velocity observations, we measure both the observed transit depth and Doppler semiamplitude to 0.5% precision. With Robo-AO and Palomar/PHARO adaptive optics imaging as well as TripleSpec spectroscopy, we measure a model-dependent mass for LHS 6343 C of 62.1 +/- 1.2 M_Jup and a radius of 0.783 +/- 0.011 R_Jup. We detect the secondary eclipse in the Kepler data at 3.5 sigma, measuring e cos omega = 0.0228 +/- 0.0008. We also derive a method to measure the mass and radius of a star and transiting companion directly, without any reliance on stellar models. The mass and radius of both objects depend only on the orbital period, stellar density, reduced semimajor axis, Doppler semiamplitude, eccentricity, and inclination, as well as the knowledge that the primary star falls on the main sequence. With this method, we calculate a model-independent mass and radius for LHS 6343 C to a precision of 3% and 2%, respectively.

We report the first planet discovery from the two-wheeled Kepler (K2) mission: HIP 116454 b. The host star HIP 116454 is a bright (V = 10.1, K = 8.0) K1-dwarf with high proper motion, and a parallax-based distance of 55.2 +/- 5.4 pc. Based on high-resolution optical spectroscopy, we find that the host star is metal-poor with [Fe/H] = -.16 +/- .18, and has a radius R = 0.716 +/- .0024 R_sun and mass M = .775 +/- .027 Msun. The star was observed by the Kepler spacecraft during its Two-Wheeled Concept Engineering Test in February 2014. During the 9 days of observations, K2 observed a single transit event. Using a new K2 photometric analysis technique we are able to correct small telescope drifts and recover the observed transit at high confidence, corresponding to a planetary radius of Rp = 2.53 +/- 0.18 Rearth. Radial velocity observations with the HARPS-N spectrograph reveal a 11.82 +/- 1.33 Mearth planet in a 9.1 day orbit, consistent with the transit depth, duration, and ephemeris. Follow-up photometric measurements from the MOST satellite confirm the transit observed in the K2 photometry and provide a refined ephemeris, making HIP 116454 b amenable for future follow-up observations of this latest addition to the growing population of transiting super-Earths around nearby, bright stars.

We present high spatial resolution images that demonstrate the hot Jupiters WASP-12b and HAT-P-8b orbit the primary star of hierarchical triple star systems. In each case, two distant companions with colors and brightness consistent with M dwarfs co-orbit the planet host as well as one another. Our adaptive optics images spatially resolve the secondary around WASP-12, previously identified by Bergfors et al. 2011 and Crossfield et al. 2012, into two distinct sources separated by 84.3+/-0.6 mas (21 +/- 3 AU). We find that the secondary to HAT-P-8, also identified by Bergfors et al. 2011, is in fact composed of two stars separated by 65.3+/-0.5 mas (15+/-1 AU). Our follow-up observations demonstrate physical association through common proper-motion. HAT-P-8 C has a particularly low mass, which we estimate to be 0.18+/-0.02Msun using photometry. Due to their hierarchy, WASP-12 BC and HAT-P-8 BC will enable the first dynamical mass determination for hot Jupiter stellar companions. These previously well-studied planet hosts now represent higher-order multi-star systems with potentially complex dynamics, underscoring the importance of diffraction-limited imaging and providing additional context for understanding the migrant population of transiting hot Jupiters.

In this paper we search for distant massive companions to known transiting gas giant planets that may have influenced the dynamical evolution of these systems. We present new radial velocity observations for a sample of 51 planets obtained using the Keck HIRES instrument, and find statistically significant accelerations in fifteen systems. Six of these systems have no previously reported accelerations in the published literature: HAT-P-10, HAT-P-22, HAT-P-29, HAT-P-32, WASP-10, and XO-2. We combine our radial velocity fits with Keck NIRC2 adaptive optics (AO) imaging data to place constraints on the allowed masses and orbital periods of the companions responsible for the detected accelerations. The estimated masses of the companions range between 1-500 M_Jup, with orbital semi-major axes typically between 1-75 AU. A significant majority of the companions detected by our survey are constrained to have minimum masses comparable to or larger than those of the transiting planets in these systems, making them candidates for influencing the orbital evolution of the inner gas giant. We estimate a total occurrence rate of 51 +/- 10% for companions with masses between 1-13 M_Jup and orbital semi-major axes between 1-20 AU in our sample. We find no statistically significant difference between the frequency of companions to transiting planets with misaligned or eccentric orbits and those with well-aligned, circular orbits. We combine our expanded sample of radial velocity measurements with constraints from transit and secondary eclipse observations to provide improved measurements of the physical and orbital characteristics of all of the planets included in our survey.

Doppler-based planet surveys have discovered numerous giant planets but are incomplete beyond several AU. At larger star-planet separations, direct planet detection through high-contrast imaging has proven successful, but this technique is sensitive only to young planets and characterization relies upon theoretical evolution models. Here we demonstrate that radial velocity measurements and high-contrast imaging can be combined to overcome these issues. The presence of widely separated companions can be deduced by identifying an acceleration (long-term trend) in the radial velocity of a star. By obtaining high spatial resolution follow-up imaging observations, we rule out scenarios in which such accelerations are caused by stellar binary companions with high statistical confidence. We report results from an analysis of Doppler measurements of a sample of 111 M-dwarf stars with a median of 29 radial velocity observations over a median time baseline of 11.8 yr. By targeting stars that exhibit a radial velocity acceleration ("trend") with adaptive optics imaging, we determine that 6.5% +/- 3.0% of M-dwarf stars host one or more massive companions with 1 < m/M_Jupiter < 13 and 0 < a < 20 AU. These results are lower than analyses of the planet occurrence rate around higher-mass stars. We find the giant planet occurrence rate is described by a double power law in stellar mass M and metallicity F = [Fe/H] such that f(M,F) = 0.039(+0.056,-0.028) M^(0.8(+1.1,-0.9)) 10^((3.8 +/- 1.2)F). Our results are consistent with gravitational microlensing measurements of the planet occurrence rate; this study represents the first model-independent comparison with microlensing observations.

Stars hosting hot Jupiters are often observed to have high obliquities, whereas stars with multiple co-planar planets have been seen to have low obliquities. This has been interpreted as evidence that hot-Jupiter formation is linked to dynamical disruption, as opposed to planet migration through a protoplanetary disk. We used asteroseismology to measure a large obliquity for Kepler-56, a red giant star hosting two transiting co-planar planets. These observations show that spin-orbit misalignments are not confined to hot-Jupiter systems. Misalignments in a broader class of systems had been predicted as a consequence of torques from wide-orbiting companions, and indeed radial-velocity measurements revealed a third companion in a wide orbit in the Kepler-56 system.

Kepler's immense photometric precision to date was maintained through satellite stability and precise pointing. In this white paper, we argue that image modeling--fitting the Kepler-downlinked raw pixel data--can vastly improve the precision of Kepler in pointing-degraded two-wheel mode. We argue that a non-trivial modeling effort may permit continuance of photometry at 10-ppm-level precision. We demonstrate some baby steps towards precise models in both data-driven (flexible) and physics-driven (interpretably parameterized) modes. We demonstrate that the expected drift or jitter in positions in the two-weel era will help with constraining calibration parameters. In particular, we show that we can infer the device flat-field at higher than pixel resolution; that is, we can infer pixel-to-pixel variations in intra-pixel sensitivity. These results are relevant to almost any scientific goal for the repurposed mission; image modeling ought to be a part of any two-wheel repurpose for the satellite. We make other recommendations for Kepler operations, but fundamentally advocate that the project stick with its core mission of finding and characterizing Earth analogs. [abridged]

In today's mailing, Hogg et al. propose image modeling techniques to maintain 10-ppm-level precision photometry in Kepler data with only two working reaction wheels. While these results are relevant to many scientific goals for the repurposed mission, all modeling efforts so far have used a toy model of the Kepler telescope. Because the two-wheel performance of Kepler remains to be determined, we advocate for the consideration of an alternate strategy for a >1 year program that maximizes the science return from the "low-torque" fields across the ecliptic plane. Assuming we can reach the precision of the original Kepler mission, we expect to detect 800 new planet candidates in the first year of such a mission. Our proposed strategy has benefits for transit timing variation and transit duration variation studies, especially when considered in concert with the future TESS mission. We also expect to help address the first key science goal of Kepler: the frequency of planets in the habitable zone as a function of spectral type.

The Kepler space telescope has opened new vistas in exoplanet discovery space by revealing populations of Earth-sized planets that provide a new context for understanding planet formation. Approximately 70% of all stars in the Galaxy belong to the diminutive M dwarf class, several thousand of which lie within Kepler's field of view, and a large number of these targets show planet transit signals. Kepler-32 is a typical star in the Kepler M dwarf sample that presents us with a rare opportunity: five planets transit this star giving us an expansive view of its architecture. All five planets of this compact system orbit their host star within a distance one third the size of Mercury's orbit with the innermost planet positioned a mere 4.3 stellar radii from the stellar photosphere. New observations limit possible false positive scenarios allowing us to validate the entire Kepler-32 system making it the richest known system of transiting planets around an M dwarf. Based on considerations of the stellar dust sublimation radius, the near period commensurability of three adjacent planets, a minimum mass protoplanetary nebula, and the volatile content inferred for the planets, we propose that the Kepler-32 planets formed at larger orbital radii and migrated inward to their present locations. The volatile content inferred for the Kepler-32 planets and order of magnitude estimates for the disk migration rates suggest these planets may have formed beyond the snow line and migrated in the presence of a gaseous disk. The Kepler-32 planets are representative of the full ensemble of planet candidates orbiting the Kepler M dwarfs for which we calculate an occurrence rate of 1.0 +/- 0.1 planet per star. The formation of the Kepler-32 planets therefore offers a plausible blueprint for the formation of one of the largest known populations of planets in our Galaxy.

Precise exoplanet characterization requires precise classification of exoplanet host stars. The masses of host stars are commonly estimated by comparing their spectra to those predicted by stellar evolution models. However, spectroscopically determined properties are difficult to measure accurately for stars that are substantially different from the Sun, such as M-dwarfs and evolved stars. Here, we propose a new method to dynamically measure the masses of transiting planets near mean-motion resonances and their host stars by combining observations of transit timing variations with radial velocity measurements. We derive expressions to analytically determine the mass of each member of the system and demonstrate the technique on the Kepler-18 system. We compare these analytic results to numerical simulations and find the two are consistent. We identify eight systems for which our technique could be applied if follow-up radial velocity measurements are collected. We conclude this analysis would be optimal for systems discovered by next generation missions similar to TESS or PLATO, which will target bright stars that are amenable to efficient RV follow-up.

We present initial results from a new high-contrast imaging program dedicated to stars that exhibit long-term Doppler radial velocity accelerations (or "trends"). The goal of the TRENDS (TaRgetting bENchmark-objects with Doppler Spectroscopy and) imaging survey is to directly detect and study the companions responsible for accelerating their host star. In this first paper of the series, we report the discovery of low-mass stellar companions orbiting HD 53665, HD 68017, and HD 71881 using NIRC2 adaptive optics (AO) observations at Keck. Follow-up imaging demonstrates association through common proper-motion. These co-moving companions have red colors with estimated spectral-types of K7--M0, M5, and M3--M4 respectively. We determine a firm lower-limit to their mass from Doppler and astrometric measurements. In the near future, it will be possible to construct three-dimensional orbits and calculate the dynamical mass of HD 68017 B and possibly HD 71881 B. We already detect astrometric orbital motion of HD 68017 B, which has a projected separation of 13.0 AU. Each companion is amenable to AO-assisted direct spectroscopy. Further, each companion orbits a solar-type star, making it possible to infer metallicity and age from the primary. Such benchmark objects are essential for testing theoretical models of cool dwarf atmospheres.