THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendon\c{c}a et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendon\c{c}a et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR.

A planet's climate can be strongly affected by its orbital eccentricity and obliquity. Here we use a 1-dimensional energy balance model modified to include a simple runaway greenhouse (RGH) parameterization to explore the effects of these two parameters on the climate of Earth-like aqua planets - completely ocean-covered planets - orbiting F-, G-, K-, and M-dwarf stars. We find that the range of instellations for which planets exhibit habitable surface conditions throughout an orbit decreases with increasing eccentricity. However, the appearance of temporarily habitable conditions during an orbit creates an eccentric habitable zone (EHZ) that is sensitive to orbital eccentricity and obliquity, planetary latitude, and host star spectral type. We find that the fraction of a planet's orbit over which it exhibits habitable surface conditions is larger on eccentric planets orbiting M-dwarf stars, due to the lower broadband planetary albedos of these planets. Planets with larger obliquities have smaller EHZs, but exhibit warmer climates if they do not enter a snowball state during their orbits. We also find no transient runaway greenhouse state on planets at all eccentricities. Rather, planets spend their entire orbits either in a RGH or not. For G-dwarf planets receiving 100% of the modern solar constant and with eccentricities above 0.55, an entire Earth ocean inventory can be lost in 3.6 Gyr. M-dwarf planets, due to their larger incident XUV flux, can become desiccated in only 690 Myr with eccentricities above 0.38. This work has important implications for eccentric planets that may exhibit surface habitability despite technically departing from the traditional habitable zone as they orbit their host stars.

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

Exoplanet hunting efforts have revealed the prevalence of exotic worlds with diverse properties, including Earth-sized bodies, which has fueled our endeavor to search for life beyond the Solar System. Accumulating experiences in astrophysical, chemical, and climatological characterization of uninhabitable planets are paving the way to characterization of potentially habitable planets. In this paper, we review our possibilities and limitations in characterizing temperate terrestrial planets with future observational capabilities through 2030s and beyond, as a basis of a broad range of discussions on how to advance "astrobiology" with exoplanets. We discuss the observability of not only the proposed biosignature candidates themselves, but also of more general planetary properties that provide circumstantial evidence, since the evaluation of any biosignature candidate relies on their context. Characterization of temperate Earth-size planets in the coming years will focus on those around nearby late-type stars. JWST and later 30 meter-class ground-based telescopes will empower their chemical investigations. Spectroscopic studies of potentially habitable planets around solar-type stars will likely require a designated spacecraft mission for direct imaging, leveraging technologies that are already being developed and tested as part of the WFIRST mission. Successful initial characterization of a few nearby targets will be an important touchstone toward a more detailed scrutiny and a larger survey that are envisioned beyond 2030. The broad outlook this paper presents may help develop new observational techniques to detect relevant features as well as frameworks to diagnose planets based on the observables.

The chemical composition of an exoplanet is a key ingredient in constraining its formation history. Iron is the most abundant transition metal, but has never been directly detected in an exoplanet due to its highly refractory nature. KELT-9b (HD 195689b) is the archetype of the class of ultra-hot Jupiters that straddle the transition between stars and gas-giant exoplanets and serve as distinctive laboratories for studying atmospheric chemistry, because of its high equilibrium temperature of 4050 +/- 180 K. These properties imply that its atmosphere is a tightly constrained chemical system that is expected to be nearly in chemical equilibrium and cloud-free. It was previously predicted that the spectral lines of iron will be detectable in the visible range of wavelengths. At these high temperatures, iron and several other transition metals are not sequestered in molecules or cloud particles and exist solely in their atomic forms. Here, we report the direct detection of atomic neutral and singly-ionized iron (Fe and Fe+), and singly-ionized titanium (Ti+) in KELT-9b via the cross-correlation technique applied to high-resolution spectra obtained during the primary transit of the exoplanet.

The atmospheric temperatures of the ultra-hot Jupiter KELT-9b straddle the transition between gas giants and stars, and therefore between two traditionally distinct regimes of atmospheric chemistry. Previous theoretical studies assume the atmosphere of KELT-9b to be in chemical equilibrium. Despite the high ultraviolet flux from KELT-9, we show using photochemical kinetics calculations that the observable atmosphere of KELT-9b is predicted to be close to chemical equilibrium, which greatly simplifies any theoretical interpretation of its spectra. It also makes the atmosphere of KELT-9b, which is expected to be cloudfree, a tightly constrained chemical system that lends itself to a clean set of theoretical predictions. Due to the lower pressures probed in transmission (compared to emission) spectroscopy, we predict the abundance of water to vary by several orders of magnitude across the atmospheric limb depending on temperature, which makes water a sensitive thermometer. Carbon monoxide is predicted to be the dominant molecule under a wide range of scenarios, rendering it a robust diagnostic of the metallicity when analyzed in tandem with water. All of the other usual suspects (acetylene, ammonia, carbon dioxide, hydrogen cyanide, methane) are predicted to be subdominant at solar metallicity, while atomic oxygen, iron and magnesium are predicted to have relative abundances as high as 1 part in 10,000. Neutral atomic iron is predicted to be seen through a forest of optical and near-infrared lines, which makes KELT-9b suitable for high-resolution ground-based spectroscopy with HARPS-N or CARMENES. We summarize future observational prospects of characterizing the atmosphere of KELT-9b.

Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball states within the habitable zone. The large amplitude changes in obliquity and eccentricity cause the ice edge, the lowest latitude extent of the ice caps, to become unstable and grow to the equator. We apply an analytical theory of the ice edge latitude to show that obliquity is the primary driver of the instability. The thermal inertia of the ice sheets and the spectral energy distribution of the G dwarf star increase the sensitivity of the model to triggering runaway glaciation. Finally, we apply a machine learning algorithm to demonstrate how this technique can be used to extend the power of climate models. This work illustrates the importance of orbital evolution for habitability in dynamically rich planetary systems. We emphasize that as potentially habitable planets are discovered around G dwarfs, we need to consider orbital dynamics.

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.

For the first time in human history, we will soon be able to apply the scientific method to the question "Are We Alone?" The rapid advance of exoplanet discovery, planetary systems science, and telescope technology will soon allow scientists to search for life beyond our Solar System through direct observation of extrasolar planets. This endeavor will occur alongside searches for habitable environments and signs of life within our Solar System. While the searches are thematically related and will inform each other, they will require separate observational techniques. The search for life on exoplanets holds potential through the great diversity of worlds to be explored beyond our Solar System. However, there are also unique challenges related to the relatively limited data this search will obtain on any individual world. This white paper reviews the scientific community's ability to use data from future telescopes to search for life on exoplanets. This material summarizes products from the Exoplanet Biosignatures Workshop Without Walls (EBWWW). The EBWWW was constituted by a series of online and in person activities, with participation from the international exoplanet and astrobiology communities, to assess state of the science and future research needs for the remote detection of life on planets outside our Solar System.

The obliquity of the Earth, which controls our seasons, varies by only ~2.5 degrees over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth's ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semi-analytic model that allows us to map broad regions of parameter space. We find that in general, the largest drivers of obliquity variations are secular spin-orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.

Here we review how environmental context can be used to interpret whether O2 is a biosignature in extrasolar planetary observations. This paper builds on the overview of current biosignature research discussed in Schwieterman et al. (2017), and provides an in-depth, interdisciplinary example of biosignature identification and observation that serves as a basis for the development of the general framework for biosignature assessment described in Catling et al., (2017). O2 is a potentially strong biosignature that was originally thought to be an unambiguous indicator for life at high-abundance. We describe the coevolution of life with the early Earth's environment, and how the interplay of sources and sinks in the planetary environment may have resulted in suppression of O2 release into the atmosphere for several billion years, a false negative for biologically generated O2. False positives may also be possible, with recent research showing potential mechanisms in exoplanet environments that may generate relatively high abundances of atmospheric O2 without a biosphere being present. These studies suggest that planetary characteristics that may enhance false negatives should be considered when selecting targets for biosignature searches. Similarly our ability to interpret O2 observed in an exoplanetary atmosphere is also crucially dependent on environmental context to rule out false positive mechanisms. We describe future photometric, spectroscopic and time-dependent observations of O2 and the planetary environment that could increase our confidence that any observed O2 is a biosignature, and help discriminate it from potential false positives. By observing and understanding O2 in its planetary context we can increase our confidence in the remote detection of life, and provide a model for biosignature development for other proposed biosignatures.

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 N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.

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.

We analyzed Kepler short-cadence M dwarf observations. Spectra from the ARC 3.5m telescope identify magnetically active (H$\alpha$ in emission) stars. The active stars are of mid-M spectral type, have numerous flares, and well-defined rotational modulation due to starspots. The inactive stars are of early-M type, exhibit less starspot signature, and have fewer flares. A Kepler to U-band energy scaling allows comparison of the Kepler flare frequency distributions with previous ground-based data. M dwarfs span a large range of flare frequency and energy, blurring the distinction between active and inactive stars designated solely by the presence of H$\alpha$. We analyzed classical and complex (multiple peak) flares on GJ 1243, finding strong correlations between flare energy, amplitude, duration and decay time, with only a weak dependence on rise time. Complex flares last longer and have higher energy at the same amplitude, and higher energy flares are more likely to be complex. A power law fits the energy distribution for flares with log $E_{K_p} >$ 31 ergs, but the predicted number of low energy flares far exceeds the number observed, at energies where flares are still easily detectable, indicating that the power law distribution may flatten at low energy. There is no correlation of flare occurrence or energy with starspot phase; the flare waiting time distribution is consistent with flares occurring randomly in time; and the energies of consecutive flares are uncorrelated. These observations support a scenario where many independent active regions on the stellar surface are contributing to the observed flare rate.