We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

We present self-consistent three-dimensional climate simulations of possible habitable states for the newly discovered Habitable Zone Earth-sized planet, TOI-700 d. We explore a variety of atmospheric compositions, pressures, and rotation states for both ocean-covered and completely desiccated planets in order to assess the planet's potential for habitability. For all 20 of our simulated cases, we use our climate model outputs to synthesize transmission spectra, combined-light spectra, and integrated broadband phase curves. These climatologically-informed observables will help the community assess the technological capabilities necessary for future characterization of this planet - as well as similar transiting planets discovered in the future - and will provide a guide for distinguishing possible climate states if one day we do obtain sensitive spectral observations of a habitable planet around a M-star. We find that TOI-700 d is a strong candidate for a habitable world and can potentially maintain temperate surface conditions under a wide variety of atmospheric compositions. Unfortunately, the spectral feature depths from the resulting transmission spectra and the peak flux and variations from our synthesized phase curves for TOI-700 d do not exceed 10 ppm. This will likely prohibit the James Webb Space Telescope (JWST) from characterizing its atmosphere; however, this motivates the community to invest in future instrumentation that perhaps can one day reveal the true nature of TOI-700 d, and to continue to search for similar planets around less distant stars.

We report on the validation of two planets orbiting the nearby (36pc) M2 dwarf TOI-1266 observed by the TESS mission. The inner planet is sub-Neptune-sized ($R=2.46 \pm 0.08 R_\oplus$) with an orbital period of 10.9 days. The outer planet has a radius of $1.67_{-0.11}^{+0.09} R_\oplus$ and resides in the exoplanet Radius Valley---the transition region between rocky and gaseous planets. With an orbital period of 18.8 days, the outer planet receives an insolation flux of 2.4 times that of Earth, similar to the insolation of Venus. Using precision near-infrared radial velocities with the Habitable-zone Planet Finder Spectrograph, we place upper mass limits of $15.9 M_\oplus$ and $6.4 M_\oplus$ at 95% confidence for the inner and outer planet, respectively. A more precise mass constraint of planet c, achievable with current RV instruments given the host star brightness (V=12.9, J=9.7), will yield further insights into the dominant processes sculpting the exoplanet Radius Valley.

We present $Spitzer$ 4.5$\mu$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s).

Explaining the evidence for surface liquid water on early Mars has been a challenge for climate modelers, as the sun was ~30% less luminous during the late-Noachian. We propose that the additional greenhouse forcing of CO2-H2 collision-induced absorption is capable of bringing the surface temperature above freezing and can put early Mars into a limit-cycling regime. Limit cycles occur when insolation is low and CO2 outgassing rates are unable to balance with the rapid drawdown of CO2 during warm weathering periods. Planets in this regime will alternate between global glaciation and transient warm climate phases. This mechanism is capable of explaining the geomorphological evidence for transient warm periods in the martian record. Previous work has shown that collision-induced absorption of CO2-H2 was capable of deglaciating early Mars, but only with high H2 outgassing rates (greater than ~600 Tmol/yr) and at high surface pressures (between 3 to 4 bars). We used new theoretically derived collision-induced absorption coefficients for CO2-H2 to reevaluate the climate limit cycling hypothesis for early Mars. Using the new and stronger absorption coefficients in our 1-dimensional radiative convective model as well as our energy balance model, we find that limit cycling can occur with an H2 outgassing rate as low as ~300 Tmol/yr at surface pressures below 3 bars. Our results agree more closely with paleoparameters for early martian surface pressure and hydrogen abundance.

Discoveries of terrestrial, Earth-sized exoplanets that lie within the habitable zone (HZ) of their host stars continue to occur at increasing rates. Transit spectroscopy can potentially enable the detection of molecular signatures from such worlds, providing an indication of the presence of an atmosphere and its chemical composition, including gases potentially indicative of a biosphere. Such planets around nearby M-dwarf stars - such as TRAPPIST-1 - provide relatively good signal, high signal/noise, and frequent transits for follow-up spectroscopy. However, even with these advantages, transit spectroscopy of terrestrial planets in the HZ of nearby M-stars will still be a challenge. Herein, we examine the potential for future space observatories to conduct such observations, using a Global Climate Model (GCM), a photochemical model, and a radiative transfer suite to simulate modern-Earth-like atmospheric boundary conditions on TRAPPIST-1e. The detectability of biosignatures on such an atmosphere via transmission spectroscopy is modeled for JWST, LUVOIR, HabEx, and Origins. We show that for any of these observatories, only CO2 would be detectable at the 3 sigma level in transmission spectroscopy, when clouds are included in our simulations. This is because the impacts of clouds on scale height strongly limits the detectability of molecules in the atmosphere. Synergies between space- and ground-based spectroscopy may be essential in order to overcome these difficulties.

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.

Robust atmospheric and radiative transfer modeling will be required to properly interpret reflected light and thermal emission spectra of terrestrial exoplanets. This will help break observational degeneracies between the numerous atmospheric, planetary, and stellar factors that drive planetary climate. Here we simulate the climates of Earth-like worlds around the Sun with increasingly slow rotation periods, from Earth-like to fully Sun-synchronous, using the ROCKE-3D general circulation model. We then provide these results as input to the Spectral Planet Model (SPM), which employs the SMART radiative transfer model to simulate the spectra of a planet as it would be observed from a future space-based telescope. We find that the primary observable effects of slowing planetary rotation rate are the altered cloud distributions, altitudes, and opacities which subsequently drive many changes to the spectra by altering the absorption band depths of biologically-relevant gas species (e.g., H2O, O2, and O3). We also identify a potentially diagnostic feature of synchronously rotating worlds in mid-infrared H2O absorption/emission lines.

Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well.

The search for spectroscopic biosignatures with the next-generation of space telescopes could provide observational constraints on the abundance of exoplanets with signs of life. An extension of this spectroscopic characterization of exoplanets is the search for observational evidence of technology, known as technosignatures. Current mission concepts that would observe biosignatures from ultraviolet to near-infrared wavelengths could place upper limits on the fraction of planets in the galaxy that host life, although such missions tend to have relatively limited capabilities of constraining the prevalence of technosignatures at mid-infrared wavelengths. Yet search-ing for technosignatures alongside biosignatures would provide important knowledge about the future of our civilization. If planets with technosignatures are abundant, then we can increase our confidence that the hardest step in planetary evolution--the Great Filter--is probably in our past. But if we find that life is commonplace while technosignatures are absent, then this would in-crease the likelihood that the Great Filter awaits to challenge us in the future.

The search for water-rich Earth-sized exoplanets around low-mass stars is rapidly gaining attention because they represent the best opportunity to characterize habitable planets in the near future. Understanding the atmospheres of these planets and determining the optimal strategy for characterizing them through transmission spectroscopy with our upcoming instrumentation is essential in order to constrain their environments. For this study, we present simulated transmission spectra of tidally locked Earth-sized ocean-covered planets around late-M to mid-K stellar spectral types, utilizing GCM modeling results previously published by Kopparapu et al. (2017) as inputs for our radiative transfer calculations performed using NASA's Planetary Spectrum Generator (psg.gsfc.nasa.gov; Villanueva et al. (2018)). We identify trends in the depth of H$_2$O spectral features as a function of planet surface temperature and rotation rate. These trends allow us to calculate the exposure times necessary to detect water vapor in the atmospheres of aquaplanets through transmission spectroscopy with the upcoming James Webb Space Telescope (JWST) as well as several future flagship space telescope concepts under consideration (LUVOIR and OST) for a target list constructed from the TESS Input Catalog (TIC). Our calculations reveal that transmission spectra for water-rich Earth-sized planets around low-mass stars will be dominated by clouds, with spectral features < 20 ppm, and only a small subset of TIC stars would allow for the characterization of an ocean planet in the Habitable Zone. We thus present a careful prioritization of targets that are most amenable to follow-up characterizations with next-generation instrumentation, in order to assist the community in efficiently utilizing precious telescope time.

The collision-induced fundamental vibration-rotation band at 6.4 um is the most significant absorption feature from O2 in the infrared (Timofeyev and Tonkov, 1978; Rinslandet al., 1982, 1989), yet it has not been previously incorporated into exoplanet spectral analyses for several reasons. Either CIAs were not included or incomplete/obsolete CIA databases were used. Also, the current version of HITRAN does not include CIAs at 6.4 um with other collision partners (O2-X). We include O2-X CIA features in our transmission spectroscopy simulations by parameterizing the 6.4 um O2-N2 CIA based on Rinsland et al.(1989) and the O2-CO2 CIA based on Baranov et al. (2004). Here we report that the O2-X CIA may be the most detectable O2 feature for transit observations. For a potentialTRAPPIST-1e analogue system within 5 pc of the Sun, it could be the only O2 detectable signature with JWST (using MIRI LRS) for a modern Earth-like cloudy atmosphere with biological quantities of O2. Also, we show that the 6.4 um O2-X CIA would be prominent for O2-rich desiccated atmospheres (Luger and Barnes, 2015) and could be detectable with JWST in just a few transits. For systems beyond 5 pc, this feature could therefore be a powerful discriminator of uninhabited planets with non-biological "false positive" O2 in their atmospheres - as they would only be detectable at those higher O2 pressures.

The TRAPPIST-1 system, consisting of an ultra-cool host star having seven known Earth-size planets will be a prime target for atmospheric characterization with JWST. However, the detectability of atmospheric molecular species may be severely impacted by the presence of clouds and/or hazes. In this work, we perform 3-D General Circulation Model (GCM) simulations with the LMD Generic model supplemented by 1-D photochemistry simulations at the terminator with the Atmos model to simulate several possible atmospheres for TRAPPIST-1e, 1f and 1g: 1) modern Earth, 2) Archean Earth, and 3) CO2-rich atmospheres. JWST synthetic transit spectra were computed using the GSFC Planetary Spectrum Generator (PSG). We find that TRAPPIST-1e, 1f and 1g atmospheres, with clouds and/or hazes, could be detected using JWST's NIRSpec prism from the CO2 absorption line at 4.3 um in less than 15 transits at 3 sigma or less than 35 transits at 5 sigma. However, our analysis suggests that other gases would require hundreds (or thousands) of transits to be detectable. We also find that H2O, mostly confined in the lower atmosphere, is very challenging to detect for these planets or similar systems if the planets' atmospheres are not in a moist greenhouse state. This result demonstrates that the use of GCMs, self-consistently taking into account the effect of clouds and sub-saturation, is crucial to evaluate the detectability of atmospheric molecules of interest as well as for interpreting future detections in a more global (and thus robust and relevant) approach.

Planets that revolve around a binary pair of stars are known as circumbinary planets. The orbital motion of the stars around their center of mass causes a periodic variation in the total instellation incident upon a circumbinary planet. This study uses both an analytic and numerical energy balance model to calculate the extent to which this effect can drive changes in surface temperature on circumbinary terrestrial planets. We show that the amplitude of the temperature variation is largely constrained by the effective heat capacity, which corresponds to the ocean-to-land ratio on the planet. Planets with large ocean fractions should experience only modest warming and cooling of only a few degrees, which suggests that habitability cannot be precluded for such circumbinary planets. Planets with large land fractions that experience extreme periodic forcing could be prone to changes in temperature of tens of degrees or more, which could drive more extreme climate changes that inhibit continuously habitable conditions.

Habitability is a measure of an environment's potential to support life, and a habitable exoplanet supports liquid water on its surface. However, a planet's success in maintaining liquid water on its surface is the end result of a complex set of interactions between planetary, stellar, planetary system and even Galactic characteristics and processes, operating over the planet's lifetime. In this chapter, we describe how we can now determine which exoplanets are most likely to be terrestrial, and the research needed to help define the habitable zone under different assumptions and planetary conditions. We then move beyond the habitable zone concept to explore a new framework that looks at far more characteristics and processes, and provide a comprehensive survey of their impacts on a planet's ability to acquire and maintain habitability over time. We are now entering an exciting era of terrestiral exoplanet atmospheric characterization, where initial observations to characterize planetary composition and constrain atmospheres is already underway, with more powerful observing capabilities planned for the near and far future. Understanding the processes that affect the habitability of a planet will guide us in discovering habitable, and potentially inhabited, planets.

We use a one-dimensional (1-D) cloud-free climate model to estimate habitable zone (HZ) boundaries for terrestrial planets of masses 0.1 M$_{E}$ and 5 M$_{E}$ around circumbinary stars of various spectral type combinations. Specifically, we consider binary systems with host spectral types F-F, F-G, F-K, F-M, G-G, G-K, G-M, K-K, K-M and M-M. Scaling the background N2 atmospheric pressure with the radius of the planet, we find that the inner edge of the HZ moves inwards towards the star for 5ME compared to 0.1ME planets for all spectral types. This is because the water-vapor column depth is smaller for larger planets and higher temperatures are needed before water vapor completely dominates the outgoing longwave radiation. The outer edge of the HZ changes little due to competing effects of the albedo and greenhouse effect. While these results are broadly consistent with the trend of single star HZ results for different mass planets, there are significant differences between single star and binary star systems for the inner edge of the HZ. Interesting combinations of stellar pairs from our 1-D model results can be used to explore for in-depth climate studies with 3-D climate models. We identify a common HZ stellar flux domain for all circumbinary spectral types

This White Paper presents a mission concept called MIRECLE - the Mid-InfraRed Exoplanet CLimate Explorer. With a moderately sized aperture of 2 meters, broad wavelength coverage (4 - 25 um), and next generation instruments, MIRECLE will be capable of efficiently characterizing a statistically significant sample of terrestrial planets, many of which will be in their host stars's habitable zones. Spectroscopic characterization of terrestrial atmospheres will provide constraints for the distribution of planets with tenuous vs. substantial atmospheres, on the inner and outer edges of the habitable zone, and climate models to assess the potential for habitability. For the few brightest targets, the detection of specific combinations of molecules would provide evidence of biosignatures. For all other targets, this comprehensive survey would filter out the airless, desiccated, or lifeless worlds, thus providing a subset of potentially habitable worlds ready for in-depth atmospheric characterization using a larger aperture telescope.

Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location of water vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.

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.

Several of NASA missions (TESS, JWST, WFIRST, etc.) and mission concepts (LUVOIR, HabEx, and OST) emphasize the exploration and characterization of exoplanets, and the study of the interstellar medium. We anticipate that a much broader set of chemical environments exists on exoplanets, necessitating data from a correspondingly broader set of chemical reactions. Similarly, the conditions that exist in astrophysical environments are very different from those traditionally probed in laboratory chemical kinetics studies. These are areas where quantum mechanical theory, applied to important reactions via well-validated chemical kinetics models, can fill a critical knowledge gap. Quantum chemical calculations are also introduced to study interior of planets, photochemical escape, and many critical chemical pathways (e.g. prebiotic environments, contaminations, etc.) After years of development of the relevant quantum chemical theories and significant advances in computational power, quantum chemical simulations have currently matured enough to describe real systems with an accuracy that competes with experiments. These approaches, therefore, have become the best possible alternative in many circumstances where performing experiments is too difficult, too expensive, or too dangerous, or simply not possible. In this white paper, several existing quantum chemical studies supporting exoplanetary science, planetary astronomy, and astrophysics are described, and the potential positive impacts of improved models associated with scientific goals of missions are addressed. In the end, a few recommendations from the scientific community to strengthen related research efforts at NASA are provided.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

As we begin to discover rocky planets in the habitable zone of nearby stars with missions like TESS and CHEOPS, we will need quick advancements on instrumentation and observational techniques that will enable us to answer key science questions, such as What are the atmospheric characteristics of habitable zone rocky planets? How common are Earth-like biosignatures in rocky planets?} How similar or dissimilar are those planets to Earth? Over the next decade we expect to have discovered several Earth-analog candidates, but we will not have the tools to study the atmospheres of all of them in detail. Ground-based ELTs can identify biosignatures in the spectra of these candidate exo-Earths and understand how the planets' atmospheres compare to the Earth at different epochs. Transit spectroscopy, high-resolution spectroscopy, and reflected-light direct imaging on ELTs can identify multiple biosignatures for habitable zone, rocky planets around M stars at optical to near-infrared wavelengths. Thermal infrared direct imaging can detect habitable zone, rocky planets around AFGK stars, identifying ozone and motivating reflected-light follow-up observations with NASA missions like HabEx/LUVOIR. Therefore, we recommend that the Astro2020 Decadal Survey Committee support: (1) the search for Earth-like biosignatures on rocky planets around nearby stars as a key science case; (2) the construction over the next decade of ground-based Extremely Large Telecopes (ELTs), which will provide the large aperture and spatial resolution necessary to start revealing the atmospheres of Earth-analogues around nearby stars; (3) the development of instrumentation that optimizes the detection of biosignatures; and (4) the generation of accurate line lists for potential biosignature gases, which are needed as model templates to detect those molecules.

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.

While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote sensing of planetary spectroscopic and temporal signals, along with careful fitting of this data to advanced models of planets and their atmospheres. Terrestrial planets are small and significantly more challenging to observe compared to their larger gaseous brethren; however observatories coming on-line in the coming decade will begin to allow their characterization. Still, it is not enough to invest only in observational endeavors. Comprehensive modeling of planetary atmospheres is required in order to fully understand what it is that our grand telescopes see in the night-sky. In our quest to characterize habitable, and possibly inhabited worlds, 3D general circulation models (GCMs) should be used to evaluate potential climate states and their associated temporal and spatial dependent observable signals. 3D models allow for coupled, self-consistent, multi-dimensional simulations, which can realistically simulate the climates of terrestrial extrasolar planets. A complete theoretical understanding of terrestrial exoplanetary atmospheres, gained through comprehensive 3D modeling, is critical for interpreting spectra of exoplanets taken from current and planned instruments, and is critical for designing future missions that aim to measure spectra of potentially habitable exoplanets as one of their key science goals. We recommend continued institutional support for 3D GCM modeling teams that focus on planetary and exoplanetary applications.

Transit spectroscopy of terrestrial planets around nearby M dwarfs is a primary goal of space missions in coming decades. 3-D climate modeling has shown that slow-synchronous rotating terrestrial planets may develop thick clouds at the substellar point, increasing the albedo. For M dwarfs with Teff > 3000 K, such planets at the inner habitable zone (IHZ) have been shown to retain moist greenhouse conditions, with enhanced stratospheric water vapor (fH2O > 1E-3) and low Earth-like surface temperatures. However, M dwarfs also possess strong UV activity, which may effectively photolyze stratospheric H2O. Prior modeling efforts have not included the impact of high stellar UV activity on the H2O. Here, we employ a 1-D photochemical model with varied stellar UV, to assess whether H2O destruction driven by high stellar UV would affect its detectability in transmission spectroscopy. Temperature and water vapor profiles are taken from published 3-D climate model simulations for an IHZ Earth-sized planet around a 3300 K M dwarf with an N2-H2O atmosphere; they serve as self-consistent input profiles for the 1-D model. We explore additional chemical complexity within the 1-D model by introducing other species into the atmosphere. We find that as long as the atmosphere is well-mixed up to 1 mbar, UV activity appears to not impact detectability of H2O in the transmission spectrum. The strongest H2O features occur in the JWST MIRI instrument wavelength range and are comparable to the estimated systematic noise floor of ~50 ppm.

A planet's atmospheric constituents (e.g., O$_2$, O$_3$, H$_2$O$_v$, CO$_2$, CH$_4$, N$_2$O) can provide clues to its surface habitability, and may offer biosignature targets for remote life detection efforts. The plethora of rocky exoplanets found by recent transit surveys (e.g., the Kepler mission) indicates that potentially habitable systems orbiting K- and M-dwarf stars may have very different orbital and atmospheric characteristics than Earth. To assess the physical distribution and observational prospects of various biosignatures and habitability indicators, it is important to understand how they may change under different astrophysical and geophysical configurations, and to simulate these changes with models that include feedbacks between different subsystems of a planet's climate. Here we use a three-dimensional (3D) Chemistry-Climate model (CCM) to study the effects of changes in stellar spectral energy distribution (SED), stellar activity, and planetary rotation on Earth-analogs and tidally-locked planets. Our simulations show that, apart from shifts in stellar SEDs and UV radiation, changes in illumination geometry and rotation-induced circulation can influence the global distribution of atmospheric biosignatures. We find that the stratospheric day-to-night side mixing ratio differences on tidally-locked planets remain low ($<20\%$) across the majority of the canonical biosignatures. Interestingly however, secondary photosynthetic biosignatures (e.g., C$_2$H$_6$S) show much greater (${\sim}67\%$) day-to-night side differences, and point to regimes in which tidal-locking could have observationally distinguishable effects on phase curve, transit, and secondary eclipse measurements. Overall, this work highlights the potential and promise for 3D CCMs to study the atmospheric properties and habitability of terrestrial worlds.

1D climate models are less sophisticated than 3D global circulation models (GCMs), however their computational time is much less expensive, allowing a large number of runs in a short period of time to explore a wide parameter space. Exploring parameter space is particularly important for predicting the observable properties of exoplanets, for which few parameters are known with certainty. Therefore, 1D climate models are still very useful tools for planetary studies. In most of these 1D models, clouds are not physically represented in the atmosphere, despite having a well-known, significant impact on a planetary radiative budget. This impact is simulated by artificially raising surface albedo, in order to reproduce the observed-averaged surface temperature (i.e. 288 K for modern Earth) and a radiative balance at the top of the atmosphere. This non-physical representation of clouds, causes atmospheric longwave and shortwaves fluxes to not match observational data. Additionally, this technique represents a parameter that is highly-tuned to modern Earth's climate, and may not be appropriate for planets that deviate from modern Earth's climate conditions. In this paper, we present an update to the climate model within the Atmos 1D atmospheric modeling package with a physical representation of clouds. We show that this physical representation of clouds in the atmosphere allows both longwave and shortwave fluxes to match observational data. This improvement will allow us to study the energy fluxes for a variety of cloudy rocky planets, and increase our confidence in future simulations of temperature profile and net energy balance.

The Kepler mission found hundreds of planet candidates within the habitable zones (HZ) of their host star, including over 70 candidates with radii larger than 3 Earth radii ($R_\oplus$) within the optimistic habitable zone (OHZ) (Kane et al. 2016). These giant planets are potential hosts to large terrestrial satellites (or exomoons) which would also exist in the HZ. We calculate the occurrence rates of giant planets ($R_p =$~3.0--25~$R_\oplus$) in the OHZ and find a frequency of $(6.5 \pm 1.9)\%$ for G stars, $(11.5 \pm 3.1)\%$ for K stars, and $(6 \pm 6)\%$ for M stars. We compare this with previously estimated occurrence rates of terrestrial planets in the HZ of G, K and M stars and find that if each giant planet has one large terrestrial moon then these moons are less likely to exist in the HZ than terrestrial planets. However, if each giant planet holds more than one moon, then the occurrence rates of moons in the HZ would be comparable to that of terrestrial planets, and could potentially exceed them. We estimate the mass of each planet candidate using the mass-radius relationship developed by Chen & Kipping (2016). We calculate the Hill radius of each planet to determine the area of influence of the planet in which any attached moon may reside, then calculate the estimated angular separation of the moon and planet for future imaging missions. Finally, we estimate the radial velocity semi-amplitudes of each planet for use in follow up observations.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

The discovery of a truly habitable exoplanet would be one of the most important events in the history of science. However, the nature and distribution of habitable environments on exoplanets is currently unconstrained. The exoplanet revolution teaches us to expect surprises. Thus, versatile, capable observatories, and multiple observation techniques are needed to study the full diversity of habitable environments. Here, we summarize the challenges and opportunities of observing planets orbiting M dwarf vs. FGK dwarfs, which are best targeted with different methods.

Future NASA concept missions that are currently under study, like Habitable Exoplanet Imaging Mission (HabEx) & Large Ultra-Violet Optical Infra Red (LUVOIR) Surveyor, would discover a large diversity of exoplanets. We propose here a classification scheme that distinguishes exoplanets into different categories based on their size and incident stellar flux, for the purpose of providing the expected number of exoplanets observed (yield) with direct imaging missions. The boundaries of this classification can be computed using the known chemical behavior of gases and condensates at different pressures and temperatures in a planetary atmosphere. In this study, we initially focus on condensation curves for sphalerite ZnS, H2O, CO2 and CH4. The order in which these species condense in a planetary atmosphere define the boundaries between different classes of planets. Broadly, the planets are divided into rocky (0.5 - 1.0RE), super-Earths (1.0- 1.75RE), sub-Neptunes (1.75-3.5RE), sub-Jovians (3.5 - 6.0RE) and Jovians (6-14.3RE) based on their planet sizes, and 'hot', 'warm' and 'cold' based on the incident stellar flux. We then calculate planet occurrence rates within these boundaries for different kinds of exoplanets, \eta_{planet}, using the community co-ordinated results of NASA's Exoplanet Program Analysis Group's Science Analysis Group-13 (SAG-13). These occurrence rate estimates are in turn used to estimate the expected exoplanet yields for direct imaging missions of different telescope diameter.

We investigate the atmospheric dynamics of terrestrial planets in synchronous rotation within the habitable zone of low-mass stars using the Community Atmosphere Model (CAM). The surface temperature contrast between day and night hemispheres decreases with an increase in incident stellar flux, which is opposite the trend seen on gas giants. We define three dynamical regimes in terms of the equatorial Rossby deformation radius and the Rhines length. The slow rotation regime has a mean zonal circulation that spans from day to night side, with both the Rossby deformation radius and the Rhines length exceeding planetary radius, which occurs for planets around stars with effective temperatures of 3300 K to 4500 K (rotation period > 20 days). Rapid rotators have a mean zonal circulation that partially spans a hemisphere and with banded cloud formation beneath the substellar point, with the Rossby deformation radius is less than planetary radius, which occurs for planets orbiting stars with effective temperatures of less than 3000 K (rotation period < 5 days). In between is the Rhines rotation regime, which retains a thermally-direct circulation from day to night side but also features midlatitude turbulence-driven zonal jets. Rhines rotators occur for planets around stars in the range of 3000 K to 3300 K (rotation period ~ 5 to 20 days), where the Rhines length is greater than planetary radius but the Rossby deformation radius is less than planetary radius. The dynamical state can be observationally inferred from comparing the morphology of the thermal emission phase curves of synchronously rotating planets.

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.

Terrestrial planets in the habitable zones (HZs) of low-mass stars and cool dwarfs have received significant scrutiny recently because their shorter orbital periods increase their chances of detection and characterization compared to planets around G-dwarfs. As these planets are likely tidal-locked, improved 3D numerical simulations of such planetary atmospheres are needed to guide target selection. Here we use a 3-D climate system model, updated with new water-vapor absorption coefficients derived from the HITRAN 2012 database, to study ocean covered planets at the inner edge of the HZ around late-M to mid-K stars ($2600$ K <= Teff <= 4500K). Our results indicate that these updated water-vapor coefficients result in significant warming compared to previous studies, so the inner HZ around M-dwarfs is not as close as suggested by earlier work. Assuming synchronously rotating planets with background 1 bar N2 atmospheres, we find that planets at the inner HZ of stars with Teff > 3000K undergo the classical "moist-greenhouse" (H2O mixing ratio > 10-3 in the stratosphere) at significantly lower surface temperature (~ 280K) in our 3-D model compared with 1-D climate models (~ 340K). This implies that some planets around low mass stars can simultaneously undergo water-loss and remain habitable. However, for star with Teff <= 3000K, planets at the inner HZ may directly transition to a runaway state, while bypassing the moist greenhouse water-loss entirely. We analyze transmission spectra of planets in a moist green-house regime, and find that there are several prominent H2O features within JWST instruments range. Thus, relying only upon standard Earth-analog spectra with 24-hour rotation period around M-dwarfs for habitability studies will miss the strong H2O features that one would expect to see on synchronously rotating planets around M-dwarf stars, with JWST.

This chapter draws upon astronomical observations and modeling to constrain the prevalence of communicative civilizations in the galaxy. We discuss the dependence of the Drake equation parameters on the spectral type of the host star and the time since the galaxy formed, which allow us to examine trajectories for the emergence of communicative civilizations over the history of the galaxy. We suggest that the maximum lifetime of communicative civilizations depends on the spectral type of the host star, which implies that F- and G-dwarf stars are the best places to search for signs of technological intelligence today.

Conventional definitions of habitability require abundant liquid surface water to exist continuously over geologic timescales. Water in each of its thermodynamic phases interacts with solar and thermal radiation and is the cause for strong climatic feedbacks. Thus, assessments of the habitable zone require models to include a complete treatment of the hydrological cycle over geologic time. Here, we use the Community Atmosphere Model from the National Center for Atmospheric Research to study the evolution of climate for an Earth-like planet at constant CO2, under a wide range of stellar fluxes from F-, G-, and K-dwarf main sequence stars. Around each star we find four stable climate states defined by mutually exclusive global mean surface temperatures (Ts); snowball (Ts < 235 K), waterbelt (235 K < Ts < 250 K), temperate (275 K < Ts < 315 K), and moist greenhouse (Ts > 330 K). Each is separated by abrupt climatic transitions. Waterbelt, temperate, and cooler moist greenhouse climates can maintain open-ocean against both sea-ice albedo and hydrogen escape processes respectively, and thus constitute habitable worlds. We consider the warmest possible habitable planet as having Ts ~ 355 K, at which point diffusion limited water-loss could remove an Earth ocean in ~1 Gyr. Without long timescale regulation of non-condensable greenhouse species at Earth-like temperatures and pressures, such as CO2, habitability can be maintained for an upper limit of ~2.2, ~2.4 and ~4.7 Gyr around F-, G- and K-dwarf stars respectively due to main sequence brightening.

For decades, scientists have tried to explain the evidence for fluvial activity on early Mars, but a consensus has yet to emerge regarding the mechanism for producing it. One hypothesis suggests early Mars was warmed by a thick greenhouse atmosphere. Another suggests that early Mars was generally cold but was warmed occasionally by impacts or by episodes of enhanced volcanism. These latter hypotheses struggle to produce the amounts of rainfall needed to form the martian valleys, but are consistent with inferred low rates of weathering compared to Earth. Here, we provide a geophysical mechanism that could have induced cycles of glaciation and deglaciation on early Mars. Our model produces dramatic climate cycles with extended periods of glaciation punctuated by warm periods lasting up to 10 Myr, much longer than those generated in other episodic warming models. The cycles occur because stellar insolation was low, and because CO2 outgassing is not able to keep pace with CO2 consumption by silicate weathering followed by deposition of carbonates. While CO2 by itself is not able to deglaciate early Mars in our model, we assume that the greenhouse effect is enhanced by substantial amounts of H2 outgassed from Mars' reduced crust and mantle. Our hypothesis can be tested by future Mars exploration that better establishes the time scale for valley formation.

The NASA Kepler mission has discovered thousands of new planetary candidates, many of which have been confirmed through follow-up observations. A primary goal of the mission is to determine the occurrance rate of terrestrial-size planets within the Habitable Zone (HZ) of their host stars. Here we provide a list of HZ exoplanet candidates from the Kepler Data Release 24 Q1-Q17 data vetting process. This work was undertaken as part of the Kepler Habitable Zone Working Group. We use a variety of criteria regarding HZ boundaries and planetary sizes to produce complete lists of HZ candidates, including a catalog of 104 candidates within the optimistic HZ and 20 candidates with radii less than two Earth radii within the conservative HZ. We cross-match our HZ candidates with the Data Release 25 stellar properties and confirmed planet properties to provide robust stellar parameters and candidate dispositions. We also include false positive probabilities recently calculated by Morton et al. (2016) for each of the candidates within our catalogs to aid in their validation. Finally, we performed dynamical analysis simulations for multi-planet systems that contain candidates with radii less than two Earth radii as a step toward validation of those systems.

The population of known extrasolar planets includes giant and terrestrial planets that closely orbit their host star. Such planets experience significant tidal distortions that can force the planet into synchronous rotation. The combined effects of tidal deformation and centripetal acceleration induces significant asphericity in the shape of these planets, compared to the mild oblateness of Earth, with maximum gravitational acceleration at the poles. Here we show that this latitudinal variation in gravitational acceleration is relevant for modeling the climate of oblate planets including Jovian planets within the solar system, closely-orbiting hot Jupiters, and planets within the habitable zone of white dwarfs. We compare first- and third-order approximations for gravitational acceleration on an oblate spheroid and calculate the geostrophic wind that would result from this asphericity on a range of solar system planets and exoplanets. Third-order variations in gravitational acceleration are negligible for Earth but become significant for Jupiter, Saturn, and Jovian exoplanets. This latitudinal variation in gravitational acceleration can be measured remotely, and the formalism presented here can be implemented for use in general circulation climate modeling studies of exoplanet atmospheres.

The liquid water habitable zone (HZ) describes the orbital distance at which a terrestrial planet can maintain above-freezing conditions through regulation by the carbonate-silicate cycle. Recent calculations have suggested that planets in the outer regions of the habitable zone cannot maintain stable, warm climates, but rather should oscillate between long, globally glaciated states and shorter periods of climatic warmth. Such conditions, similar to 'Snowball Earth' episodes experienced on Earth, would be inimical to the development of complex land life, including intelligent life. Here, we build upon previous studies with an updated an energy balance climate model to calculate this 'limit cycle' region of the habitable zone where such cycling would occur. We argue that an abiotic Earth would have a greater CO$_2$ partial pressure than today because plants and other biota help to enhance the storage of CO$_2$ in soil. When we tune our abiotic model accordingly, we find that limit cycles can occur but that previous calculations have overestimated their importance. For G stars like the Sun, limit cycles occur only for planets with CO$_2$ outgassing rates less than that on modern Earth. For K and M star planets, limit cycles should not occur; however, M-star planets may be inhospitable to life for other reasons. Planets orbiting late G-type and early K-type stars retain the greatest potential for maintaining warm, stable conditions. Our results suggest that host star type, planetary volcanic activity, and seafloor weathering are all important factors in determining whether planets will be prone to limit cycling.

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

Terrestrial planets at the inner edge of the habitable zone of late-K and M-dwarf stars are expected to be in synchronous rotation, as a consequence of strong tidal interactions with their host stars. Previous global climate model (GCM) studies have shown that, for slowly-rotating planets, strong convection at the substellar point can create optically thick water clouds, increasing the planetary albedo, and thus stabilizing the climate against a thermal runaway. However these studies did not use self-consistent orbital/rotational periods for synchronously rotating planets placed at different distances from the host star. Here we provide new estimates of the inner edge of the habitable zone for synchronously rotating terrestrial planets around late-K and M-dwarf stars using a 3-D Earth-analog GCM with self-consistent relationships between stellar metallicity, stellar effective temperature, and the planetary orbital/rotational period. We find that both atmospheric dynamics and the efficacy of the substellar cloud deck are sensitive to the precise rotation rate of the planet. Around mid-to-late M-dwarf stars with low metallicity, planetary rotation rates at the inner edge of the HZ become faster, and the inner edge of the habitable zone is farther away from the host stars than in previous GCM studies. For an Earth-sized planet, the dynamical regime of the substellar clouds begins to transition as the rotation rate approaches ~10 days. These faster rotation rates produce stronger zonal winds that encircle the planet and smear the substellar clouds around it, lowering the planetary albedo, and causing the onset of the water-vapor greenhouse climatic instability to occur at up to ~25% lower incident stellar fluxes than found in previous GCM studies. For mid-to-late M-dwarf stars with high metallicity and for mid-K to early-M stars, we agree with previous studies.

A radiative-convective climate model is used to calculate stratospheric temperatures and water vapor concentrations for ozone-free atmospheres warmer than that of modern Earth. Cold, dry stratospheres are predicted at low surface temperatures, in agreement with recent 3-D calculations. However, at surface temperatures above 350 K, the stratosphere warms and water vapor becomes a major upper atmospheric constituent, allowing water to be lost by photodissociation and hydrogen escape. Hence, a 'moist greenhouse' explanation for loss of water from Venus, or some exoplanet receiving a comparable amount of stellar radiation, remains a viable hypothesis. Temperatures in the upper parts of such atmospheres are well below those estimated for a gray atmosphere, and this factor should be taken into account when performing 'inverse' climate calculations to determine habitable zone boundaries using 1-D models.

Earth-like planets within the liquid water habitable zone of M type stars may evolve into synchronous rotators. On these planets, the sub-stellar hemisphere experiences perpetual daylight while the opposing anti-stellar hemisphere experiences perpetual darkness. Because the night-side hemisphere has no direct source of energy, the air over this side of the planet is prone to freeze out and deposit on the surface, which could result in atmospheric collapse. However, general circulation models (GCMs) have shown that atmospheric dynamics can counteract this problem and provide sufficient energy transport to the anti-stellar side. Here we use an idealized GCM to consider the impact of geothermal heating on the habitability of synchronously rotating planets. Geothermal heating may be expected due to tidal interactions with the host star, and the effects of geothermal heating provide additional habitable surface area and may help to induce melting of ice on the anti-stellar hemisphere. We also explore the persistence of atmospheric asymmetries between the northern and southern hemispheres, and we find that the direction of the meridional circulation (for rapidly rotating planets) or the direction of zonal wind (for slowly rotating planets) reverses on either side of the sub-stellar point. We show that the zonal circulation approaches a theoretical state similar to a Walker circulation only for slowly rotating planets, while rapidly rotating planets show a zonal circulation with the opposite direction. We find that a cross-polar circulation is present in all cases and provides an additional mechanism of mass and energy transport from the sub-stellar to anti-stellar point. Characterization of the atmospheres of synchronously rotating planets should include consideration of hemispheric differences in meridional circulation and examination of transport due to cross-polar flow.

The field of exoplanetary science has seen a dramatic improvement in sensitivity to terrestrial planets over recent years. Such discoveries have been a key feature of results from the {\it Kepler} mission which utilizes the transit method to determine the size of the planet. These discoveries have resulted in a corresponding interest in the topic of the Habitable Zone (HZ) and the search for potential Earth analogs. Within the Solar System, there is a clear dichotomy between Venus and Earth in terms of atmospheric evolution, likely the result of the large difference ($\sim$ factor of two) in incident flux from the Sun. Since Venus is 95\% of the Earth's radius in size, it is impossible to distinguish between these two planets based only on size. In this paper we discuss planetary insolation in the context of atmospheric erosion and runaway greenhouse limits for planets similar to Venus. We define a ``Venus Zone'' (VZ) in which the planet is more likely to be a Venus analog rather than an Earth analog. We identify 43 potential Venus analogs with an occurrence rate ($\eta_{\venus}$) of $0.32^{+0.05}_{-0.07}$ and $0.45^{+0.06}_{-0.09}$ for M dwarfs and GK dwarfs respectively.

The presence of valleys on ancient terrains of Mars suggest that liquid water flowed on the martian surface 3.8 billion years ago or before. The above-freezing temperatures required to explain valley formation could have been transient, in response to frequent large meteorite impacts on early Mars, or they could have been caused by long-lived greenhouse warming. Climate models that consider only the greenhouse gases carbon dioxide and water vapor have been unable to recreate warm surface conditions, given the lower solar luminosity at that time. Here we use a one-dimensional climate model to demonstrate that an atmosphere containing 1.3-4 bar of CO2 and water vapor, along with 5 to 20 percent H2, could have raised the mean surface temperature of early Mars above the freezing point of water. Vigorous volcanic outgassing from a highly reduced early martian mantle is expected to provide sufficient atmospheric H2 and CO2, the latter from the photochemical oxidation of outgassed CH4 and CO, to form a CO2-H2 greenhouse. Such a dense early martian atmosphere is consistent with independent estimates of surface pressure based on cratering data.

The ongoing discoveries of extrasolar planets are unveiling a wide range of terrestrial mass (size) planets around their host stars. In this letter, we present estimates of habitable zones (HZs) around stars with stellar effective temperatures in the range 2600 K - 7200 K, for planetary masses between 0.1 ME and 5 ME. Assuming H2O (inner HZ) and CO2 (outer HZ) dominated atmospheres, and scaling the background N2 atmospheric pressure with the radius of the planet, our results indicate that larger planets have wider HZs than do smaller ones. Specifically, with the assumption that smaller planets will have less dense atmospheres, the inner edge of the HZ (runaway greenhouse limit) moves outward (~10% lower than Earth flux) for low mass planets due to larger greenhouse effect arising from the increased H2O column depth. For larger planets, the H2O column depth is smaller, and higher temperatures are needed before water vapor completely dominates the outgoing longwave radiation. Hence the inner edge moves inward (7% higher than Earth's flux). The outer HZ changes little due to the competing effects of the greenhouse effect and an increase in albedo. New, 3-D climate model results from other groups are also summarized, and we argue that further, independent studies are needed to verify their predictions. Combined with our previous work, the results presented here provide refined estimates of HZs around main-sequence stars and provide a step towards a more comprehensive analysis of HZs.

The habitable zone (HZ) around a star is typically defined as the region where a rocky planet can maintain liquid water on its surface. That definition is appropriate, because this allows for the possibility that carbon-based, photosynthetic life exists on the planet in sufficient abundance to modify the planet's atmosphere in a way that might be remotely detected. Exactly what conditions are needed, however, to maintain liquid water remains a topic for debate. Historically, modelers have restricted themselves to water-rich planets with CO2 and H2O as the only important greenhouse gases. More recently, some researchers have suggested broadening the definition to include arid, 'Dune' planets on the inner edge and planets with captured H2 atmospheres on the outer edge, thereby greatly increasing the HZ width. Such planets could exist, but we demonstrate that an inner edge limit of 0.59 AU or less is physically unrealistic. We further argue that conservative HZ definitions should be used for designing future space-based telescopes, but that optimistic definitions may be useful in interpreting the data from such missions. In terms of effective solar flux, Seff, the recently recalculated HZ boundaries are: recent Venus-1.78, runaway greenhouse-1.04, moist greenhouse-1.01, maximum greenhouse-0.35, early Mars-0.32. Based on a combination of different HZ definitions, the frequency of potentially Earth-like planets around late-K and M stars observed by Kepler is in the range of 0.4-0.5.

This paper has been withdrawn by the author in favor of a private review.

If mutual gravitational scattering among exoplanets occurs, then it may produce unique orbital properties. For example, two-planet systems that lie near the boundary between circulation and libration of their periapses could result if planet-planet scattering ejected a former third planet quickly, leaving one planet on an eccentric orbit and the other on a circular orbit. We first improve upon previous work that examined the apsidal behavior of known multiplanet systems by doubling the sample size and including observational uncertainties. This analysis recovers previous results that demonstrated that many systems lay on the apsidal boundary between libration and circulation. We then performed over 12,000 three-dimensional N-body simulations of hypothetical three-body systems that are unstable, but stabilize to two-body systems after an ejection. Using these synthetic two-planet systems, we test the planet-planet scattering hypothesis by comparing their apsidal behavior, over a range of viewing angles, to that of the observed systems and find that they are statistically consistent regardless of the multiplicity of the observed systems. Finally, we combine our results with previous studies to show that, from the sampled cases, the most likely planetary mass function prior to planet-planet scattering follows a power law with index -1.1. We find that this pre-scattering mass function predicts a mutual inclination frequency distribution that follows an exponential function with an index between -0.06 and -0.1.

Because of their large numbers, low mass stars may be the most abundant planet hosts in our Galaxy. Furthermore, terrestrial planets in the habitable zones (HZs) around M-dwarfs can potentially be characterized in the near future and hence may be the first such planets to be studied. Recently Dressing & Charbonneau(2013) used Kepler data and calculated the frequency of terrestrial planets in the HZ of cool stars to be 0.15^{+0.13}_{-0.06} per star for Earth-size planets (0.5-1.4 R_{Earth}). However, this estimate was derived using the Kasting et al.(1993) HZ limits, which were not valid for stars with effective temperatures lower than 3700 K. Here we update their result using new HZ limits from Kopparapu et al.(2013) for stars with effective temperatures between 2600 K and 7200 K, which includes the cool M stars in the Kepler target list. The new habitable zone boundaries increase the number of planet candidates in the habitable zone. Assuming Earth-size planets as 0.5 - 1.4 R_{Earth}, when we reanalyze their results, we obtain a terrestrial planet frequency of 0.48^{+0.12}_{-0.24} and 0.53^{+0.08}_{-0.17} planets per M-dwarf star for conservative and optimistic limits of the HZ boundaries, respectively. Assuming Earth-size planets as 0.5 - 2 R_{Earth}, the frequency increases to 0.51^{+0.10}_{-0.20} per star for the conservative estimate and to 0.61^{+0.07}_{-0.15} per star for the optimistic estimate. Within uncertainties, our optimistic estimates are in agreement with a similar optimistic estimate from the radial velocity survey of M-dwarfs (0.41^{+0.54}_{-0.13}, Bonfils et al.(2011)). So, the potential for finding Earth-like planets around M stars may be higher than previously reported.

Identifying terrestrial planets in the habitable zones (HZs) of other stars is one of the primary goals of ongoing radial velocity and transit exoplanet surveys and proposed future space missions. Most current estimates of the boundaries of the HZ are based on 1-D, cloud-free, climate model calculations by Kasting et al.(1993). The inner edge of the HZ in Kasting et al.(1993) model was determined by loss of water, and the outer edge was determined by the maximum greenhouse provided by a CO2 atmosphere. A conservative estimate for the width of the HZ from this model in our Solar system is 0.95-1.67 AU. Here, an updated 1-D radiative-convective, cloud-free climate model is used to obtain new estimates for HZ widths around F, G, K and M stars. New H2O and CO2 absorption coefficients, derived from the HITRAN 2008 and HITEMP 2010 line-by-line databases, are important improvements to the climate model. According to the new model, the water loss (inner HZ) and maximum greenhouse (outer HZ) limits for our Solar System are at 0.99 AU and 1.70 AU, respectively, suggesting that the present Earth lies near the inner edge. Additional calculations are performed for stars with effective temperatures between 2600 K and 7200 K, and the results are presented in parametric form, making them easy to apply to actual stars. The new model indicates that, near the inner edge of the HZ, there is no clear distinction between runaway greenhouse and water loss limits for stars with T_{eff} ~< 5000 K which has implications for ongoing planet searches around K and M stars. To assess the potential habitability of extrasolar terrestrial planets, we propose using stellar flux incident on a planet rather than equilibrium temperature. Our model does not include the radiative effects of clouds; thus, the actual HZ boundaries may extend further in both directions than the estimates just given.

Within the next decade, ground based gravitational wave detectors are in principle capable of determining the compact object merger rate per unit volume of the local universe to better than 20% with more than 30 detections. We argue that the stellar models are sensitive to heterogeneities (in age and metallicity at least) in such a way that the predicted merger rates are subject to an additional 30-50% systematic errors unless these heterogeneities are taken into account. Without adding new electromagnetic constraints on massive binary evolution or relying on more information from each merger (e.g., binary masses and spins), as few as the $\simeq 5$ merger detections could exhaust the information available in a naive comparison to merger rate predictions. As a concrete example, we use a nearby-universe catalog to demonstrate that no one tracer of stellar content can constrain merger rates without introducing a systematic error of order $O(30%)$ at 90% confidence. More generally, we argue that theoretical binary evolution can depend sufficiently sensitively on star-forming conditions -- even assuming no uncertainty in binary evolution model -- that the \emph{distribution} of star forming conditions must be incorporated to reduce the systematic error in merger rate predictions below roughly 40%. (Abridged)

We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date.

Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes.

Binaries that contain a hot subdwarf (sdB) star and a main sequence companion may have interacted in the past. This binary population has historically helped determine our understanding of binary stellar evolution. We have computed a grid of binary population synthesis models using different assumptions about the minimum core mass for helium ignition, the envelope binding energy, the common envelope ejection efficiency, the amount of mass and angular momentum lost during stable mass transfer, and the criteria for stable mass transfer on the red giant branch and in the Hertzsprung gap. These parameters separately and together can significantly change the entire predicted population of sdBs. Nonetheless, several different parameter sets can reproduce the observed subpopulation of sdB + white dwarf and sdB + M dwarf binaries, which has been used to constrain these parameters in previous studies. The period distribution of sdB + early F dwarf binaries offers a better test of different mass transfer scenarios for stars that fill their Roche lobes on the red giant branch.

Extraterrestrial technology may exist in the Solar System without our knowledge. This is because the vastness of space, combined with our limited searches to date, implies that any remote unpiloted exploratory probes of extraterrestrial origin would likely remain unnoticed. Here we develop a probabilistic approach to quantify our certainty (or uncertainty) of the existence of such technology in the Solar System. We discuss some possible strategies for improving this uncertainty that include analysis of moon- and Mars-orbiting satellite data as well as continued exploration of the Solar System.

The hot Jupiter WASP-12b is a heavily irradiated exoplanet in a short period orbit around a G0-star with twice the metallicity of the Sun. A recent thermochemical equilibrium analysis based on Spitzer and ground-based infrared observations suggests that the presence of $\ch4$ in its atmosphere and the lack of $\h2o$ features can only be explained if the carbon-to-oxygen ratio in the planet's atmosphere is much greater than the solar ratio ($\ctoo = 0.54$). Here, we use a 1-D photochemical model to study the effect of disequilibrium chemistry on the observed abundances of $\h2o, \com, \co2$ and $\ch4$ in the WASP-12b atmosphere. We consider two cases: one with solar $\ctoo$ and another with $\ctoo = 1.08$. The solar case predicts that $\h2o$ and $\com$ are more abundant than $\co2$ and $\ch4$, as expected, whereas the high $\ctoo$ model shows that $\com$, C$_{2}$H$_{2}$ and HCN are more abundant. This indicates that the extra carbon from the high $\ctoo$ model is in hydrocarbon species. $\h2o$ photolysis is the dominant disequilibrium mechanism that alters the chemistry at higher altitudes in the solar $\ctoo$ case, whereas photodissociation of C$_{2}$H$_{2}$ and HCN is significant in the super-solar case. Furthermore, our analysis indicates that $\c2h2$ is the major absorber in the atmosphere of WASP-12b and the absorption features detected near 1.6 and 8 micron may be arising from C$_{2}$H$_{2}$ rather than $\ch4$. The Hubble Space Telescope's WFC3 can resolve this discrepancy, as $\c2h2$ has absorption between $1.51 - 1.54$ microns, while $\ch4$ does not.

The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data from the LIGO S5 science run: one optimized for pointlike sources and one for arbitrary extended sources. The latter result is the first of its kind. Finding no evidence to support the detection of GWs, we present 90% confidence level (CL) upper-limit maps of GW strain power with typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2 Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on pointlike sources constitute a factor of 30 improvement over the previous best limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain from interesting targets including Sco X-1, SN1987A and the Galactic Center as low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These limits are the most constraining to date and constitute a factor of 5 improvement over the previous best limits.

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures and excitation of non-radial modes which would emit gravitational waves (GWs). We present the results of a search for GW bursts from six galactic magnetars that is sensitive to neutron star f-modes, thought to be the most efficient GW emitting oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1 kpc from Earth, an order of magnitude closer than magnetars targeted in previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an estimated isotropic energy >10^{44} erg which is comparable to the giant flares. We find no evidence of GWs associated with a sample of 1279 electromagnetic triggers from six magnetars occurring between November 2006 and June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest model-dependent GW emission energy upper limits for band- and time-limited white noise bursts in the detector sensitive band, and for f-mode ringdowns (at 1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively, where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516. These limits on GW emission from f-modes are an order of magnitude lower than any previous, and approach the range of electromagnetic energies seen in SGR giant flares for the first time.

We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis and value of the wave polarization angle of, respectively, $1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95% degree of belief. These limits are below the indirect {\it spin-down limit} of $3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred from observed decrease in Vela's spin frequency, and correspond to a limit on the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but still well below the spin-down limit, are obtained assuming the star's spin axis inclination and the wave polarization angles are unknown.

The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer Gravitational-wave observatory (LIGO) were operational and taking data as part of the fifth LIGO science run (S5). We present the first direct search for the gravitational wave emission associated with oscillations of the fundamental quadrupole mode excited by a pulsar timing glitch. No gravitational wave detection candidate was found. We place Bayesian 90% confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain amplitude of gravitational wave ring-down signals, depending on which spherical harmonic mode is excited. The corresponding range of energy upper limits is 5.0e44 to 1.3e45 erg.

We present a search for periodic gravitational waves from the neutron star in the supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day interval taken from the fifth science run of the Laser Interferometer Gravitational-Wave Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers a wide range of first and second frequency derivatives appropriate for the age of the remnant and for different spin-down mechanisms. No gravitational wave signal was detected. Within the range of search frequencies, we set 95% confidence upper limits of 0.7--1.2e-24 on the intrinsic gravitational wave strain, 0.4--4e-4 on the equatorial ellipticity of the neutron star, and 0.005--0.14 on the amplitude of r-mode oscillations of the neutron star. These direct upper limits beat indirect limits derived from energy conservation and enter the range of theoretical predictions involving crystalline exotic matter or runaway r-modes. This is the first gravitational wave search to present upper limits on r-modes.

We study the dynamical stability of planetary systems consisting of one hypothetical terrestrial mass planet ($1 $ or $10 \mearth$) and one massive planet ($10 \mearth - 10 \mjup$). We consider masses and orbits that cover the range of observed planetary system architectures (including non-zero initial eccentricities), determine the stability limit through N-body simulations, and compare it to the analytic Hill stability boundary. We show that for given masses and orbits of a two planet system, a single parameter, which can be calculated analytically, describes the Lagrange stability boundary (no ejections or exchanges) but which diverges significantly from the Hill stability boundary. However, we do find that the actual boundary is fractal, and therefore we also identify a second parameter which demarcates the transition from stable to unstable evolution. We show the portions of the habitable zones of $\rho$ CrB, HD 164922, GJ 674, and HD 7924 which can support a terrestrial planet. These analyses clarify the stability boundaries in exoplanetary systems and demonstrate that, for most exoplanetary systems, numerical simulations of the stability of potentially habitable planets are only necessary over a narrow region of parameter space. Finally we also identify and provide a catalog of known systems which can host terrestrial planets in their habitable zones.

We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs.

We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.