The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.

We model the long-term XUV luminosity of TRAPPIST-1 to constrain the evolving high-energy radiation environment experienced by its planetary system. Using Markov Chain Monte Carlo (MCMC), we derive probabilistic constraints for TRAPPIST-1's stellar and XUV evolution that account for observational uncertainties, degeneracies between model parameters, and empirical data of low-mass stars. We constrain TRAPPIST-1's mass to $m_{\star} = 0.089 \pm{0.001}$ M$_{\odot}$ and find that its early XUV luminosity likely saturated at $\log_{10}(L_{XUV}/L_{bol}) = -3.03^{+0.23}_{-0.12}$. From the posterior distribution, we infer that there is a ${\sim}40\%$ chance that TRAPPIST-1 is still in the saturated phase today, suggesting that TRAPPIST-1 has maintained high activity and $L_{XUV}/L_{bol} \approx 10^{-3}$ for several Gyrs. TRAPPIST-1's planetary system therefore likely experienced a persistent and extreme XUV flux environment, potentially driving significant atmospheric erosion and volatile loss. The inner planets likely received XUV fluxes ${\sim}10^3 - 10^4\times$ that of the modern Earth during TRAPPIST-1's ${\sim}1$ Gyr-long pre-main sequence phase. Deriving these constraints via MCMC is computationally non-trivial, so scaling our methods to constrain the XUV evolution of a larger number of M dwarfs that harbor terrestrial exoplanets would incur significant computational expenses. We demonstrate that approxposterior, an open source Python machine learning package for approximate Bayesian inference using Gaussian processes, accurately and efficiently replicates our analysis using $980\times$ less computational time and $1330\times$ fewer simulations than MCMC sampling using emcee. We find that approxposterior derives constraints with mean errors on the best fit values and $1\sigma$ uncertainties of $0.61\%$ and $5.5\%$, respectively, relative to emcee.

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.

We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formulae are crafted to be numerically stable over the expected range of usage. We additionally present analytic formulae for the partial derivatives of instantaneous flux with respect to the radius ratio, impact parameter, and limb darkening coefficients. These expressions are rapid to evaluate, and compare quite favorably in speed and accuracy to existing transit light curve codes. We also use these expressions to numerically compute the first partial derivatives of exposure-time averaged transit light curves with respect to all model parameters. An additional application is modeling eclipsing binary or eclipsing multiple star systems in cases where the stars may be treated as spherically symmetric. We provide code which implements these formulae in C++, Python, IDL, and Julia, with tests and examples of usage.

Exoplanets orbiting pre-main sequence stars are laboratories for studying planet evolution processes, including atmospheric loss, orbital migration, and radiative cooling. V1298 Tau, a young solar analog with an age of 23 $\pm$ 4 Myr, is one such laboratory. The star is already known to host a Jupiter-sized planet on a 24 day orbit. Here, we report the discovery of three additional planets --- all between the size of Neptune and Saturn --- based on our analysis of K2 Campaign 4 photometry. Planets c and d have sizes of 5.6 and 6.4 $R_\oplus$, respectively and with orbital periods of 8.25 and 12.40 days reside 0.25% outside of the nominal 3:2 mean-motion resonance. Planet e is 8.7 $R_\oplus$ in size but only transited once in the K2 time series and thus has a period longer than 36 days, but likely shorter than 223 days. The V1298 Tau system may be a precursor to the compact multiplanet systems found to be common by the Kepler mission. However, the large planet sizes stand in sharp contrast to the vast majority of Kepler multis which have planets smaller than 3 $R_\oplus$. Simple dynamical arguments suggest total masses of $<$28 $M_\oplus$ and $<$120 $M_\oplus$ for the c-d and d-b planet pairs, respectively. The implied low masses suggest that the planets may still be radiatively cooling and contracting, and perhaps losing atmosphere. The V1298 Tau system offers rich prospects for further follow-up including atmospheric characterization by transmission or eclipse spectroscopy, dynamical characterization through transit-timing variations, and measurements of planet masses and obliquities by radial velocities.

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.

We implement a search for exoplanets in campaigns zero through eight (C0-8) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C0-8 search on $1.5 \times 10^5$ target stars with magnitudes in the range of Kp = 9-15. We detect 818 transiting planet candidates, of which 374 were undiscovered by prior searches, with {64,15,5,2,1} in {2,3,4,5,6}-planet multi-planet candidate systems, respectively. Of the new planets detected, 100 orbit M dwarfs, including one which is potentially rocky and in the habitable zone. 154 of our candidates reciprocally transit with our Solar System: they are geometrically aligned to see at least one Solar System planet transit. We find candidates which display transit timing variations and dozens of candidates on both period extremes with single transits or ultra-short periods. We point to evidence that our candidates display similar patterns in frequency and size-period relation as confirmed planets, such as tentative evidence for the radius gap. Confirmation of these planet candidates with follow-up studies will increase the number of K2 planets by up to 50%, and characterization of their host stars will improve statistical studies of planet properties. Our sample includes many planets orbiting bright stars amenable for radial velocity follow-up and future characterization with JWST. We also list the 579 eclipsing binary systems detected as part of this search.

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.

We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-ubiquitous all-sky catalog photometry and Gaia parallaxes. EXOFASTv2 can model the star by itself, too. Multi-planet systems are modeled self-consistently with the same underlying stellar mass that defines their semi-major axes through Kepler's law and the planetary period. Transit timing, duration, and depth variations can be modeled with a simple command line option. We explain our methodology and rationale as well as provide an improved version of the core transit model that is both 25\% faster and more accurate. We highlight several potential pitfalls in exoplanet modeling, including the handling of eccentricity in transit-only fits, that the standard exoplanet convention for $\omega$ uses a left-handed coordinate system, contrary to most modern textbooks, how to avoid an important degeneracy when allowing negative companion masses, and a widely unappreciated, potential 10-minute ambiguity in the reported transit times. EXOFASTv2 is available at https://github.com/jdeast/EXOFASTv2 . The code is written in IDL, and includes an executable that can be run freely and legally without an IDL license or any knowledge of the language. Extensive documentation and tutorials are included in the distribution for a variety of example fits. Advanced amateurs and undergrads have successfully performed sophisticated global fits of complex planetary systems with EXOFASTv2. It is therefore a powerful tool for education and outreach as well as the broader professional community.

We suggest to use the Hubble Space Telescople (HST) follow-up observations of the TESS targets for detecting possible plasma tori along the orbits of exoplanets orbiting M dwarfs. The source of the torus could be planetary volcanic activity due to tidal or electromagnetic induction heating. Fast losses to space for planets orbiting these active stars can lead to the lost material forming a torus along the planetary orbit, similar to the Io plasma torus. We show that such torus would be potentially detectable by the HST in the UV.

We examine how tides, stellar evolution, and magnetic braking shape the rotation period (P$_{rot}$) evolution of low-mass stellar binaries up to orbital periods (P$_{orb}$) of 100 d across a wide range tidal dissipation parameters using two common equilibrium tidal models. We find that many binaries with P$_{orb} < 20$ d tidally lock, and most with $P_{orb} < 4$ d tidally lock into synchronous rotation on circularized orbits. At short P$_{orb}$, tidal torques produce a population of fast rotators that single-star only models of magnetic braking fail to produce. In many cases, we show that the competition between magnetic braking and tides produces a population of subsynchronous rotators that persists for Gyrs, even in short P$_{orb}$ binaries, qualitatively reproducing the subsynchronous eclipsing binaries (EBs) discovered in the Kepler field by Lurie et al. (2017). Both equilibrium tidal models predict that binaries can tidally-interact out to P$_{orb} \approx 80$ d, while the Constant Phase Lag tidal model predicts that binaries can tidally lock out to P$_{orb} \approx 100$ d. Tidal torques often force the P$_{rot}$ evolution of stellar binaries to depart from the long-term magnetic braking-driven spin down experienced by single stars, revealing that P$_{rot}$ is not be a valid proxy for age in all cases, i.e. gyrochronology can underpredict ages by up to $300\%$ unless one accounts for binarity. We suggest that accurate determinations of orbital eccentricties and P$_{rot}$ can be used to discriminate between which equilibrium tidal models best describes tidal interactions in low-mass binary stars.

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.

Observations of the Kepler-1625 system with the Kepler and Hubble Space Telescopes have suggested the presence of a candidate exomoon, Kepler-1625b I, a Neptune-radius satellite orbiting a long-period Jovian planet. Here we present a new analysis of the Hubble observations, using an independent data reduction pipeline. We find that the transit light curve is well fit with a planet-only model, with a best-fit $\chi^2_\nu$ equal to 1.01. The addition of a moon does not significantly improve the fit quality. We compare our results directly with the original light curve from Teachey & Kipping (2018), and find that we obtain a better fit to the data using a model with fewer free parameters (no moon). We discuss possible sources for the discrepancy in our results, and conclude that the lunar transit signal found by Teachey & Kipping (2018) was likely an artifact of the data reduction. This finding highlights the need to develop independent pipelines to confirm results that push the limits of measurement precision.

To date, a handful of exoplanets have been photometrically mapped using phase-modulated reflection or emission from their surfaces, but the small amplitudes of such signals have limited previous maps almost exclusively to coarse dipolar features on hot giant planets. In this work, we uncover a signal using recently released data from the Transiting Exoplanet Survey Satellite (TESS), which we show corresponds to time-variable reflection from a terrestrial planet with a rotation period of 0.9972696 days. Using a spherical harmonic-based reflection model developed as an extension of the STARRY package, we are able to reconstruct the surface features of this rocky world. We recover a time-variable albedo map of the planet including persistent regions which we interpret as oceans and cloud banks indicative of continental features. We argue that this planet represents the most promising detection of a habitable world to date, although the potential intelligence of any life on it is yet to be determined.

We present scope (Simulated CCD Observations for Photometric Experimentation), a Python package to create a forward model of telescope detectors and simulate stellar targets with motion relative to the CCD. The primary application of this package is the simulation of the Kepler Space Telescope detector to predict and characterize increased instrumental noise in the spacecraft's final campaigns of observation. As the fuel powering the spacecraft's stabilizing thrusters ran out and thruster fires began to sputter and fail, stellar Point Spread Functions (PSFs) experienced more extreme and less predictable motion relative to regions of varied sensitivity on the spacecraft detector, generating more noise in transiting exoplanet light curves. Using our simulations, we demonstrate that current de-trending techniques effectively capture and remove systematics caused by sensitivity variation for spacecraft motion as high as about ten times that typically experienced by K2. The scope package is open-source and has been generalized to allow custom detector and stellar target parameters. Future applications include simulating observations made by the Transiting Exoplanet Survey Satellite (TESS) and ground based observations with synthetic atmospheric interference as testbeds for noise-removal techniques.

Rotational mapping and glint are two proposed methods to directly detect liquid water on the surface of habitable exoplanets. However, false positives for both methods may prevent the unambiguous detection of exoplanet oceans. We use simulations of Earth as an exoplanet to introduce a combination of multiwavelength, multiphase, time-series direct-imaging observations and accompanying analyses that may improve the robustness of exoplanet ocean detection by spatially mapping the ocean glint signal. As the planet rotates, the glint spot appears to "blink" as Lambertian scattering continents interrupt the specular reflection from the ocean. This manifests itself as a strong source of periodic variability in crescent-phase reflected light curves. We invert these light curves to constrain the longitudinal slice maps and apparent albedo of two surfaces at both quadrature and crescent phase. At crescent phase, the retrieved apparent albedo of ocean-bearing longitudinal slices is increased by a factor of 5, compared to the albedo at quadrature phase, due to the contribution from glint. The land-bearing slices exhibit no significant change in apparent albedo with phase. The presence of forward-scattering clouds in our simulated observation increases the overall reflectivity toward crescent, but clouds do not correlate with any specific surfaces, thereby allowing for the phase-dependent glint effect to be interpreted as distinct from cloud scattering. Retrieving the same longitudinal map at quadrature and crescent phases may be used to tie changes in the apparent albedo with phase back to specific geographic surfaces, although this requires ideal geometries. We estimate that crescent-phase time-dependent glint detections are feasible for between 1-10 habitable zone exoplanets orbiting the nearest G, K, and M dwarfs using a space-based, high-contrast, direct-imaging telescope with a diameter >6 m.

We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet-planet/planet-moon occultation light curves, as well as thermal (rotational) phase curves. In this paper we also introduce STARRY, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in STARRY is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. STARRY also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map.

The TRAPPIST-1 planetary system provides an unprecedented opportunity to study terrestrial exoplanet evolution with the James Webb Space Telescope (JWST) and ground-based observatories. Since M dwarf planets likely experience extreme volatile loss, the TRAPPIST-1 planets may have highly-evolved, possibly uninhabitable atmospheres. We used a versatile, 1D terrestrial-planet climate model with line-by-line radiative transfer and mixing length convection (VPL Climate) coupled to a terrestrial photochemistry model to simulate environmental states for the TRAPPIST-1 planets. We present equilibrium climates with self-consistent atmospheric compositions, and observational discriminants of post-runaway, desiccated, 10-100 bar O2- and CO2-dominated atmospheres, including interior outgassing, as well as for water-rich compositions. Our simulations show a range of surface temperatures, most of which are not habitable, although an aqua-planet TRAPPIST-1 e could maintain a temperate surface given Earth-like geological outgassing and CO2. We find that a desiccated TRAPPIST-1 h may produce habitable surface temperatures beyond the maximum greenhouse distance. Potential observational discriminants for these atmospheres in transmission and emission spectra are influenced by photochemical processes and aerosol formation, and include collision-induced oxygen absorption (O2-O2), and O3, CO, SO2, H2O, and CH4 absorption features, with transit signals of up to 200 ppm. Our simulated transmission spectra are consistent with K2, HST, and Spitzer observations of the TRAPPIST-1 planets. For several terrestrial atmospheric compositions, we find that TRAPPIST-1 b is unlikely to produce aerosols. These results can inform JWST observation planning and data interpretation for the TRAPPIST-1 system and other M dwarf terrestrial planets.

We outline a mechanism that explains the observed lack of circumbinary planets (CBPs) via coupled stellar-tidal evolution of isolated binary stars. Tidal forces between low-mass, short-period binary stars on the pre-main sequence slow the stellar rotations, transferring rotational angular momentum to the orbit as the stars approach the tidally locked state. This transfer increases the binary orbital period, expanding the region of dynamical instability around the binary, and destabilizing CBPs that tend to preferentially orbit just beyond the initial dynamical stability limit. After the stars tidally lock, we find that angular momentum loss due to magnetic braking can significantly shrink the binary orbit, and hence the region of dynamical stability, over time impacting where surviving CBPs are observed relative to the boundary. We perform simulations over a wide range of parameter space and find that the expansion of the instability region occurs for most plausible initial conditions and that in some cases, the stability semi-major axis doubles from its initial value. We examine the dynamical and observable consequences of a CBP falling within the dynamical instability limit by running N-body simulations of circumbinary planetary systems and find that typically, at least one planet is ejected from the system. We apply our theory to the shortest period {\it Kepler} binary that possesses a CBP, Kepler-47, and find that its existence is consistent with our model. Under conservative assumptions, we find that coupled stellar-tidal evolution of pre-main sequence binary stars removes at least one close-in CBP in $87\%$ of multi-planet circumbinary systems.

We analyze the evolution of the potentially habitable planet Proxima Centauri b to identify environmental factors that affect its long-term habitability. We consider physical processes acting on size scales ranging from the galactic to the stellar system to the planet's core. We find that there is a significant probability that Proxima Centauri has had encounters with its companion stars, Alpha Centauri A and B, that are close enough to destabilize an extended planetary system. If the system has an additional planet, as suggested by the discovery data, then it may perturb planet b's eccentricity and inclination, possibly driving those parameters to non-zero values, even in the presence of strong tidal damping. We also model the internal evolution of the planet, evaluating the roles of different radiogenic abundances and tidal heating and find that magnetic field generation is likely for billions of years. We find that if planet b formed in situ, then it experienced 169 +/- 13 million years in a runaway greenhouse as the star contracted during its formation. This early phase could remove up to 5 times as much water as in the modern Earth's oceans, possibly producing a large abiotic oxygen atmosphere. On the other hand, if Proxima Centauri b formed with a substantial hydrogen atmosphere (0.01 - 1% of the planet's mass), then this envelope could have shielded the water long enough for it to be retained before being blown off itself. After modeling this wide range of processes we conclude that water retention during the host star's pre-main sequence phase is the biggest obstacle for Proxima b's habitability. These results are all obtained with a new software package called VPLANET.

We explore the occurrence and detectability of planet-planet occultations (PPOs) in exoplanet systems. These are events during which a planet occults the disk of another planet in the same system, imparting a small photometric signal as its thermal or reflected light is blocked. We focus on the planets in TRAPPIST-1, whose orbital planes we show are aligned to within 0.3 degrees at 90% confidence. We present a photodynamical model for predicting and computing PPOs in TRAPPIST-1 and other systems for various assumptions of the planets' atmospheric states. When marginalizing over the uncertainties on all orbital parameters, we find that the rate of PPOs in TRAPPIST-1 is about 1.4 per day. We investigate the prospects for detection of these events with the James Webb Space Telescope, finding that ~10-20 occultations per year of b and c should be above the noise level at 12-15 microns. Joint modeling of several of these PPOs could lead to a robust detection. Alternatively, observations with the proposed Origins Space Telescope should be able to detect individual PPOs at high signal-to-noise. We show how PPOs can be used to break transit timing variation degeneracies, imposing strong constraints on the eccentricities and masses of the planets, as well as to constrain the longitudes of nodes and thus the complete three-dimensional structure of the system. We further show how modeling of these events can be used to reveal a planet's day/night temperature contrast and construct crude surface maps. We make our photodynamical code available on github.

The target of many astronomical studies is the recovery of tiny astrophysical signals living in a sea of uninteresting (but usually dominant) noise. In many contexts (i.e., stellar time-series, or high-contrast imaging, or stellar spectroscopy), there are structured components in this noise caused by systematic effects in the astronomical source, the atmosphere, the telescope, or the detector. More often than not, evaluation of the true physical model for these nuisances is computationally intractable and dependent on too many (unknown) parameters to allow rigorous probabilistic inference. Sometimes, housekeeping data---and often the science data themselves---can be used as predictors of the systematic noise. Linear combinations of simple functions of these predictors are often used as computationally tractable models that can capture the nuisances. These models can be used to fit and subtract systematics prior to investigation of the signals of interest, or they can be used in a simultaneous fit of the systematics and the signals. In this Note, we show that if a Gaussian prior is placed on the weights of the linear components, the weights can be marginalized out with an operation in pure linear algebra, which can (often) be made fast. We illustrate this model by demonstrating the applicability of a linear model for the non-linear systematics in K2 time-series data, where the dominant noise source for many stars is spacecraft motion and variability.

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf.

We present an update to the EVEREST K2 pipeline that addresses various limitations in the previous version and improves the photometric precision of the de-trended light curves. We develop a fast regularization scheme for third order pixel level decorrelation (PLD) and adapt the algorithm to include the PLD vectors of neighboring stars to enhance the predictive power of the model and minimize overfitting, particularly for faint stars. We also modify PLD to work for saturated stars and improve its performance on extremely variable stars. On average, EVEREST 2.0 light curves have 10-20% higher photometric precision than those in the previous version, yielding the highest precision light curves at all Kp magnitudes of any publicly available K2 catalog. For most K2 campaigns, we recover the original Kepler precision to at least Kp = 14, and to at least Kp = 15 for campaigns 1, 5, and 6. We also de-trend all short cadence targets observed by K2, obtaining even higher photometric precision for these stars. All light curves for campaigns 0-8 are available online in the EVEREST catalog, which will be continuously updated with future campaigns. EVEREST 2.0 is open source and is coded in a general framework that can be applied to other photometric surveys, including Kepler and the upcoming TESS mission.

We examine the feasibility of detecting auroral emission from the potentially habitable exoplanet Proxima Centauri b. Detection of aurorae would yield an independent confirmation of the planet's existence, constrain the presence and composition of its atmosphere, and determine the planet's eccentricity and inclination, thereby breaking the mass-inclination degeneracy. If Proxima Centauri b is a terrestrial world with an Earth-like atmosphere and magnetic field, we estimate the power at the 5577\AA\ OI auroral line is on the order of 0.1 TW under steady-state stellar wind, or ${\sim} 100 {\times}$ stronger than that on Earth. This corresponds to a planet-star contrast ratio of $10^{-6}-10^{-7}$ in a narrow band about the 5577\AA\ line, although higher contrast ($10^{-4}-10^{-5}$) may be possible during periods of strong magnetospheric disturbance (auroral power $1-10$ TW). We searched the Proxima Centauri b HARPS data for the 5577\AA\ line and for other prominent oxygen and nitrogen lines, but find no signal, indicating that the OI auroral line contrast must be lower than $2\times 10^{-2}$ (with power $\lesssim$ 3,000 TW), consistent with our predictions. We find that observations of 0.1 TW auroral emission lines are likely infeasible with current and planned telescopes. However, future observations with a space-based coronagraphic telescope or a ground-based extremely large telescope (ELT) with a coronagraph could push sensitivity down to terawatt oxygen aurorae (contrast $7\times 10^{-6}$) with exposure times of ${\sim} 1$ day. If a coronagraph design contrast of $10^{-7}$ can be achieved with negligible instrumental noise, a future concept ELT could observe steady-state auroral emission in a few nights.

Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here we use 1D coupled climate-photochemical models to generate self-consistent atmospheres for evolutionary scenarios predicted in our companion paper (Barnes et al., 2016). These include high-O2, high-CO2, and more Earth-like atmospheres, with either oxidizing or reducing compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b, but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere, and JWST observations longward of 7 microns could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST, but may be within the reach of larger aperture telescopes. Direct imaging spectra may detect O4, which is diagnostic of massive water loss and O2 retention, rather than a photosynthesis. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 {\mu}m would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4 and either photosynthetically produced O2 or a hydrocarbon haze layer.

We present EVEREST, an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation (PLD) to remove systematics introduced by the spacecraft's pointing error and a Gaussian process (GP) to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0-7, yielding light curves with precision comparable to that of the original Kepler mission for stars brighter than $K_p \approx 13$, and within a factor of two of the Kepler precision for fainter targets. We perform cross-validation and transit injection and recovery tests to validate the pipeline, and compare our light curves to the other de-trended light curves available for download at the MAST High Level Science Products archive. We find that EVEREST achieves the highest average precision of any of these pipelines for unsaturated K2 stars. The improved precision of these light curves will aid in exoplanet detection and characterization, investigations of stellar variability, asteroseismology, and other photometric studies. The EVEREST pipeline can also easily be applied to future surveys, such as the TESS mission, to correct for instrumental systematics and enable the detection of low signal-to-noise transiting exoplanets. The EVEREST light curves and the source code used to generate them are freely available online.

O2 and O3 have been long considered the most robust individual biosignature gases in a planetary atmosphere, yet multiple mechanisms that may produce them in the absence of life have been described. However, these abiotic planetary mechanisms modify the environment in potentially identifiable ways. Here we briefly discuss two of the most detectable spectral discriminants for abiotic O2/O3: CO and O4. We produce the first explicit self-consistent simulations of these spectral discriminants as they may be seen by JWST. If JWST-NIRISS and/or NIRSpec observe CO (2.35, 4.6 um) in conjunction with CO2 (1.6, 2.0, 4.3 um) in the transmission spectrum of a terrestrial planet it could indicate robust CO2 photolysis and suggest that a future detection of O2 or O3 might not be biogenic. Strong O4 bands seen in transmission at 1.06 and 1.27 um could be diagnostic of a post-runaway O2-dominated atmosphere from massive H-escape. We find that for these false positive scenarios, CO at 2.35 um, CO2 at 2.0 and 4.3 um, and O4 at 1.27 um are all stronger features in transmission than O2/O3 and could be detected with SNRs $\gtrsim$ 3 for an Earth-size planet orbiting a nearby M dwarf star with as few as 10 transits, assuming photon-limited noise. O4 bands could also be sought in UV/VIS/NIR reflected light (at 0.345, 0.36, 0.38, 0.445, 0.475, 0.53, 0.57, 0.63, 1.06, and 1.27 um) by a next generation direct-imaging telescope such as LUVOIR/HDST or HabEx and would indicate an oxygen atmosphere too massive to be biologically produced.

We show that photoevaporation of small gaseous exoplanets ("mini-Neptunes") in the habitable zones of M dwarfs can remove several Earth masses of hydrogen and helium from these planets and transform them into potentially habitable worlds. We couple X-ray/extreme ultraviolet (XUV)-driven escape, thermal evolution, tidal evolution and orbital migration to explore the types of systems that may harbor such "habitable evaporated cores" (HECs). We find that HECs are most likely to form from planets with $\sim 1 M_\oplus$ solid cores with up to about 50% H/He by mass, though whether or not a given mini-Neptune forms a HEC is highly dependent on the early XUV evolution of the host star. As terrestrial planet formation around M dwarfs by accumulation of local material is likely to form planets that are small and dry, evaporation of small migrating mini-Neptunes could be one of the dominant formation mechanisms for volatile-rich Earths around these stars.

We show that terrestrial planets in the habitable zones of M dwarfs older than $\sim$ 1 Gyr could have been in runaway greenhouses for several hundred Myr following their formation due to the star's extended pre-main sequence phase, provided they form with abundant surface water. Such prolonged runaway greenhouses can lead to planetary evolution divergent from that of Earth. During this early runaway phase, photolysis of water vapor and hydrogen/oxygen escape to space can lead to the loss of several Earth oceans of water from planets throughout the habitable zone, regardless of whether the escape is energy-limited or diffusion-limited. We find that the amount of water lost scales with the planet mass, since the diffusion-limited hydrogen escape flux is proportional to the planet surface gravity. In addition to undergoing potential desiccation, planets with inefficient oxygen sinks at the surface may build up hundreds to thousands of bars of abiotically produced O$_2$, resulting in potential false positives for life. The amount of O$_2$ that builds up also scales with the planet mass; we find that O$_2$ builds up at a constant rate that is controlled by diffusion: $\sim$ 5 bars/Myr on Earth-mass planets and up to $\sim$ 25 bars/Myr on super-Earths. As a result, some recently discovered super-Earths in the habitable zone such as GJ 667Cc could have built up as many as 2000 bars of O$_2$ due to the loss of up to 10 Earth oceans of water. The fate of a given planet strongly depends on the extreme ultraviolet flux, the duration of the runaway regime, the initial water content, and the rate at which oxygen is absorbed by the surface. In general, we find that the initial phase of high luminosity may compromise the habitability of many terrestrial planets orbiting low mass stars.

The Upsilon Andromedae system is the first exoplanetary system to have the relative inclination of two planets' orbital planes directly measured, and therefore offers our first window into the 3-dimensional configurations of planetary systems. We present, for the first time, full 3-dimensional, dynamically stable configurations for the 3 planets of the system consistent with all observational constraints. While the outer 2 planets, c and d, are inclined by about 30 degrees, the inner planet's orbital plane has not been detected. We use N-body simulations to search for stable 3-planet configurations that are consistent with the combined radial velocity and astrometric solution. We find that only 10 trials out of 1000 are robustly stable on 100 Myr timescales, or about 8 billion orbits of planet b. Planet b's orbit must lie near the invariable plane of planets c and d, but can be either prograde or retrograde. These solutions predict b's mass is in the range 2 - 9 $M_{Jup}$ and has an inclination angle from the sky plane of less than 25 degrees. Combined with brightness variations in the combined star/planet light curve ("phase curve"), our results imply that planet b's radius is about 1.8 $R_{Jup}$, relatively large for a planet of its age. However, the eccentricity of b in several of our stable solutions reaches values greater than 0.1, generating upwards of $10^{19}$ watts in the interior of the planet via tidal dissipation, possibly inflating the radius to an amount consistent with phase curve observations.