Habitability has been generally defined as the capability of an environment to support life. Ecologists have been using Habitat Suitability Models (HSMs) for more than four decades to study the habitability of Earth from local to global scales. Astrobiologists have been proposing different habitability models for some time, with little integration and consistency between them and different in function to those used by ecologists. In this white paper, we suggest a mass-energy habitability model as an example of how to adapt and expand the models used by ecologists to the astrobiology field. We propose to implement these models into a NASA Habitability Standard (NHS) to standardize the habitability objectives of planetary missions. These standards will help to compare and characterize potentially habitable environments, prioritize target selections, and study correlations between habitability and biosignatures. Habitability models are the foundation of planetary habitability science. The synergy between the methods used by ecologists and astrobiologists will help to integrate and expand our understanding of the habitability of Earth, the Solar System, and exoplanets.

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in XUV irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system, and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on a UV continuum flux measurement (centered around $\sim1900$ Angstroms) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photo-dissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 $\times$10$^{-4} L_{*}$, while the full flux including non-flaring periods is 6.1 $\times$10$^{-5} L_{*}$, when $L_{*}$ is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass-loss and tidal evolution, and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio $Q'=3Q/2k_{2}$ for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator \texttt{VPLanet} to study ocean retention for these planets using our derived system parameters.

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

Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist.

M dwarfs host most of the exoplanets in the local Milky Way. Some of these planets, ranging from sub-Earths to super-Jupiters, orbit in their stars' habitable zones (HZs), although many likely possess surface environments that preclude habitability. Moreover, exomoons around these planets could harbor life for long timescales and thus may also be targets for biosignature surveys. Here we investigate the potential habitability, stability, and detectability of exomoons around exoplanets orbiting M dwarfs. We first compile an updated list of known M-dwarf exoplanet hosts, comprising 109 stars and 205 planets. For each M dwarf, we compute and update precise luminosities with the Virtual Observatory spectral energy distribution Analyzer and Gaia DR2 parallaxes to determine inner and outer boundaries of their HZs. For each planet, we retrieve (or, when necessary, homogeneously estimate) their masses and radii, calculate the long-term dynamical stability of hypothetical moons, and identify those planets that can support habitable moons. We find that 33 exoplanet candidates are located in the HZs of their host stars and that four of them could host Moon- to Titan-mass exomoons for timescales longer than the Hubble time.

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

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

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

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

The solar tide in an ancient Venusian ocean is simulated using a dedicated numerical tidal model. Simulations with varying ocean depth and rotational periods ranging from -243 to 64 sidereal Earth days are used to calculate the tidal dissipation rates and associated tidal torque. The results show that the tidal dissipation could have varied by more than 5 orders of magnitude, from 0.001--780 Gigawatts, depending on rotational period and ocean depth. The associated tidal torque is about 2 orders of magnitude below the present day Venusian atmospheric torque, and could change the Venusian day length by up to 72 days per million years depending on rotation rate. Consequently, an ocean tide on ancient Venus could have had significant effects on the rotational history of the planet. These calculations have implications for the rotational periods of similarly close-in exoplanetary worlds and the location of the inner edge of the liquid water habitable zone.

The Dharma Planet Survey (DPS) aims to monitor about 150 nearby very bright FGKM dwarfs (within 50 pc) during 2016$-$2020 for low-mass planet detection and characterization using the TOU very high resolution optical spectrograph (R$\approx$100,000, 380-900nm). TOU was initially mounted to the 2-m Automatic Spectroscopic Telescope at Fairborn Observatory in 2013-2015 to conduct a pilot survey, then moved to the dedicated 50-inch automatic telescope on Mt. Lemmon in 2016 to launch the survey. Here we report the first planet detection from DPS, a super-Earth candidate orbiting a bright K dwarf star, HD 26965. It is the second brightest star ($V=4.4$ mag) on the sky with a super-Earth candidate. The planet candidate has a mass of 8.47$\pm0.47M_{\rm Earth}$, period of $42.38\pm0.01$ d, and eccentricity of $0.04^{+0.05}_{-0.03}$. This RV signal was independently detected by Diaz et al. (2018), but they could not confirm if the signal is from a planet or from stellar activity. The orbital period of the planet is close to the rotation period of the star (39$-$44.5 d) measured from stellar activity indicators. Our high precision photometric campaign and line bisector analysis of this star do not find any significant variations at the orbital period. Stellar RV jitters modeled from star spots and convection inhibition are also not strong enough to explain the RV signal detected. After further comparing RV data from the star's active magnetic phase and quiet magnetic phase, we conclude that the RV signal is due to planetary-reflex motion and not stellar activity.

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

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 geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

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

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

How can scientists conclude with high confidence that an exoplanet hosts life? As telescopes come on line over the next 20 years that can directly observe photons from terrestrial exoplanets, this question will dictate the activities of many scientists across many fields. The expected data will be sparse and with low signal-to-noise, which will make disentangling biosignatures from abiotic features challenging. Our Earth is not just unique in that it hosts life, it is also the only terrestrial planet with direct observations of its interior through seismic waves, and compositional evolution through field and laboratory measurements. This extensive research reveals a planet born from collisions between worlds, followed by a complicated biogeochemical evolution. Exoplanet interiors, on the other hand, can only be constrained by the following observations: 1) photometric and spectroscopic analysis of the planet's atmosphere, 2) spectroscopic and photometric analysis of the host star, and 3) companion planet properties. From these (future) data, astrobiologists must generate plausible compositional and evolutionary models that constrain a potentially habitable exoplanet's internal properties and history, provide environmental context, and rule out geochemical explanations for any putative biosignatures. The goal of this white paper is to frame the role of geophysical and geochemical processes relevant to the search for life beyond the Solar System and to identify critical, but understudied, areas of future research.

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

Potentially habitable planets can orbit close enough to their host star that the differential gravity across their diameters can fix the rotation rate at a specific frequency, a process called tidal locking. Tidally locked planets on circular orbits will rotate synchronously, but those on eccentric orbits will either librate or rotate super-synchronously. I calculate how habitable planets evolve under two commonly-used models and find, for example, that one model predicts that the Earth's rotation rate would have synchronized after 4.5 Gyr if its initial rotation period was 3 days, it had no satellites, and it always maintained the modern Earth's tidal properties. Lower mass stellar hosts will induce stronger tidal effects on potentially habitable planets, and tidal locking is possible for most planets in the habitable zones of GKM dwarf stars. For fast rotating planets, both models predict eccentricity growth and that circularization can only occur once the rotational frequency is similar to the orbital frequency. The orbits of potentially habitable planets of very late M dwarfs (<0.15 solar masses) are very likely to be circularized within 1 Gyr and hence those planets will be synchronous rotators. Proxima b is almost assuredly tidally locked, but its orbit may not have circularized yet, so the planet could be rotating super-synchronously today. The evolution of the isolated and potentially habitable Kepler planet candidates is computed and about half could be tidally locked. Finally, projected TESS planets are simulated over a wide range of assumptions, and the vast majority of all cases are found to tidally lock within 1 Gyr. These results suggest that the process of tidal locking is a major factor in the evolution of most of the potentially habitable exoplanets to be discovered in the near future. [abridged]

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

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

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

As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is $0.00 \leqslant e \leqslant 0.32$, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO$_2$ in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60$^\circ$ and 90$^\circ$). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO$_2$ is increased to 5-bar levels. In a low-CO$_2$ case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO$_2$ levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet's substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.

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

We constrain the true spin-orbit alignment of the KOI-89 system by numerically fitting the two \emph{Kepler} photometric lightcurves produced by transiting planets KOI-89.01 and KOI-89.02. The two planets have periods of 84.69 days and 207.58 days, respectively. We find that the two bodies are low-density giant planets with radii $0.45 \pm 0.03~\mathrm{R_{jup}}$ and $0.43 \pm 0.05~\mathrm{R_{jup}}$ and spin-orbit misalignments $72^{\circ} \pm 3^\circ$ and $73^{\circ+11}_{-5}$, respectively. Via dynamic stability tests we demonstrate the general trend of higher system stability with the two planets close to mutual alignment and estimate their coalignment angle to $20^\circ \pm 20^\circ$ -- i.e. the planets are misaligned with the star but may be aligned with each other. From these results, we limit KOI-89's misalignment mechanisms to star-disk-binary interactions, disk warping via planet-disk interactions, planet-planet scattering, Kozai resonance, or internal gravity waves.

Exoplanet habitability is traditionally assessed by comparing a planet's semi-major axis to the location of its host star's "habitable zone," the shell around a star for which Earth-like planets can possess liquid surface water. The Kepler space telescope has discovered numerous planet candidates near the habitable zone, and many more are expected from missions such as K2, TESS and PLATO. These candidates often require significant follow-up observations for validation, so prioritizing planets for habitability from transit data has become an important aspect of the search for life in the universe. We propose a method to compare transiting planets for their potential to support life based on transit data, stellar properties and previously reported limits on planetary emitted flux. For a planet in radiative equilibrium, the emitted flux increases with eccentricity, but decreases with albedo. As these parameters are often unconstrained, there is an "eccentricity-albedo degeneracy" for the habitability of transiting exoplanets. Our method mitigates this degeneracy, includes a penalty for large-radius planets, uses terrestrial mass-radius relationships, and, when available, constraints on eccentricity to compute a number we call the "habitability index for transiting exoplanets" that represents the relative probability that an exoplanet could support liquid surface water. We calculate it for Kepler Objects of Interest and find that planets that receive between 60-90% of the Earth's incident radiation, assuming circular orbits, are most likely to be habitable. Finally, we make predictions for the upcoming TESS and JWST missions.

The internal thermal and magnetic evolution of rocky exoplanets is critical to their habitability. We focus on the thermal-orbital evolution of Earth-mass planets around low mass M stars whose radiative habitable zone overlaps with the "tidal zone". We develop a thermal-orbital evolution model calibrated to Earth that couples tidal dissipation, with a temperature-dependent Maxwell rheology, to orbital circularization and migration. We illustrate thermal-orbital steady states where surface heat flow is balanced by tidal dissipation and cooling can be stalled for billions of years until circularization occurs. Orbital energy dissipated as tidal heat in the interior drives both inward migration and circularization, with a circularization time that is inversely proportional to the dissipation rate. We identify a peak in the internal dissipation rate as the mantle passes through a visco-elastic state at mantle temperatures near 1800 K. Planets orbiting a 0.1 solar-mass star within $0.07$ AU circularize before 10 Gyr, independent of initial eccentricity. Once circular, these planets cool monotonically and maintain dynamos similar to Earth. Planets forced into eccentric orbits can experience a super-cooling of the core and rapid core solidification, inhibiting dynamo action for planets in the habitable zone. We find that tidal heating is insignificant in the habitable zone around $0.45$ (or larger) solar mass stars because tidal dissipation is a stronger function of orbital distance than stellar mass, and the habitable zone is further from larger stars. Suppression of the planetary magnetic field exposes the atmosphere to stellar wind erosion and the surface to harmful radiation. In addition to weak magnetic fields, massive melt eruption rates and prolonged magma oceans may render eccentric planets in the habitable zone of low mass stars inhospitable for life.

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

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

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

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

In this paper, we derive the fundamental properties of 1SWASPJ011351.29+314909.7 (J0113+31), a metal-poor (-0.40 +/- 0.04 dex), eclipsing binary in an eccentric orbit (~0.3) with an orbital period of ~14.277 d. Eclipsing M dwarfs orbiting solar-type stars (EBLMs), like J0113+31, have been identified from WASP light curves and follow-up spectroscopy in the course of the transiting planet search. We present the first binary of the EBLM sample to be fully analysed, and thus, define here the methodology. The primary component with a mass of 0.945 +/- 0.045 Msun has a large radius (1.378 +/- 0.058 Rsun) indicating that the system is quite old, ~9.5 Gyr. The M-dwarf secondary mass of 0.186 +/- 0.010 Msun and radius of 0.209 +/- 0.011 Rsun are fully consistent with stellar evolutionary models. However, from the near-infrared secondary eclipse light curve, the M dwarf is found to have an effective temperature of 3922 +/- 42 K, which is ~600 K hotter than predicted by theoretical models. We discuss different scenarios to explain this temperature discrepancy. The case of J0113+31 for which we can measure mass, radius, temperature and metallicity, highlights the importance of deriving mass, radius and temperature as a function of metallicity for M dwarfs to better understand the lowest mass stars. The EBLM Project will define the relationship between mass, radius, temperature and metallicity for M dwarfs providing important empirical constraints at the bottom of the main sequence.

The diversity and quantity of moons in the Solar System suggest a manifold population of natural satellites exist around extrasolar planets. Of peculiar interest from an astrobiological perspective, the number of sizable moons in the stellar habitable zones may outnumber planets in these circumstellar regions. With technological and theoretical methods now allowing for the detection of sub-Earth-sized extrasolar planets, the first detection of an extrasolar moon appears feasible. In this review, we summarize formation channels of massive exomoons that are potentially detectable with current or near-future instruments. We discuss the orbital effects that govern exomoon evolution, we present a framework to characterize an exomoon's stellar plus planetary illumination as well as its tidal heating, and we address the techniques that have been proposed to search for exomoons. Most notably, we show that natural satellites in the range of 0.1 - 0.5 Earth mass (i) are potentially habitable, (ii) can form within the circumplanetary debris and gas disk or via capture from a binary, and (iii) are detectable with current technology.

We explore the impact of obliquity variations on planetary habitability in hypothetical systems with high mutual inclination. We show that large amplitude, high frequency obliquity oscillations on Earth-like exoplanets can suppress the ice-albedo feedback, increasing the outer edge of the habitable zone. We restrict our exploration to hypothetical systems consisting of a solar-mass star, an Earth-mass planet at 1 AU, and 1 or 2 larger planets. We verify that these systems are stable for $10^8$ years with N-body simulations, and calculate the obliquity variations induced by the orbital evolution of the Earth-mass planet and a torque from the host star. We run a simplified energy balance model on the terrestrial planet to assess surface temperature and ice coverage on the planet's surface, and we calculate differences in the outer edge of the habitable zone for planets with rapid obliquity variations. For each hypothetical system, we calculate the outer edge of habitability for two conditions: 1) the full evolution of the planetary spin and orbit, and 2) the eccentricity and obliquity fixed at their average values. We recover previous results that higher values of fixed obliquity and eccentricity expand the habitable zone, but also find that obliquity oscillations further expand habitable orbits in all cases. Terrestrial planets near the outer edge of the habitable zone may be more likely to support life in systems that induce rapid obliquity oscillations as opposed to fixed-spin planets. Such planets may be the easiest to directly characterize with space-borne telescopes.

Earth-scale planets in the classical habitable zone (HZ) are more likely to be habitable if they possess active geophysics. Without a constant internal energy source, planets cool as they age, eventually terminating tectonic activity and rendering the planet sterile to life. However, for planets orbiting low-mass stars, the presence of an outer companion could generate enough tidal heat in the HZ planet to prevent such cooling. The range of mass and orbital parameters for the companion that give adequate long-term heating of the inner HZ planet, while avoiding very early total desiccation, is probably substantial. We locate the ideal location for the outer of a pair of planets, under the assumption that the inner planet has the same incident flux as Earth, orbiting example stars: a generic late M dwarf ($T_{eff}=2670 K$) and the M9V/L0 dwarf DEN1048. Thus discoveries of Earth-scale planets in the HZ zone of old small stars should be followed by searches for outer companion planets that might be essential for current habitability.

The determination of an exoplanet as rocky is critical for the assessment of planetary habitability. Observationally, the number of small-radius, transiting planets with accompanying mass measurements is insufficient for a robust determination of the transitional mass or radius. Theoretically, models predict that rocky planets can grow large enough to become gas giants when they reach ~10 Earth-masses, but the transitional mass remains unknown. Here I show how transit data, interpreted in the context of tidal theory, can reveal the critical radius that separates rocky and gaseous exoplanets. Standard tidal models predict that rocky exoplanets' orbits are tidally circularized much more rapidly than gaseous bodies', suggesting the former will tend to be found on circular orbits at larger semi-major axes than the latter. Well-sampled transits can provide a minimum eccentricity of the orbit, allowing a measurement of this differential circularization. I show that this effect should be present in the data from the Kepler spacecraft, but is not apparent. Instead, it appears that there is no evidence of tidal circularization at any planetary radius, probably because the publicly-available data, particularly the impact parameters, are not accurate enough. I also review the bias in the transit duration toward values that are smaller than that of planets on circular orbits, stressing that the azimuthal velocity of the planet determines the transit duration. The ensemble of Kepler planet candidates may be able to determine the critical radius between rocky and gaseous exoplanets, tidal dissipation as a function of planetary radius, and discriminate between tidal models.

The Kepler space telescope has detected transits of objects as small as the Earth's Moon, and moons as small as 0.2 Earth masses can be detected in the Kepler data by transit timing and transit duration variations of their host planets. Such massive moons around giant planets in the stellar habitable zone (HZ) could serve as habitats for extraterrestrial life. We here assess the danger of exomoons to be in a runaway greenhouse (RG) state due to extensive heating from the planet. We apply pre-computed evolution tracks for giant planets to calculate the incident planetary radiation on the moon as a function of time. The total energy budget of stellar flux, illumination from the planet, and tidal heating in the satellite is compared to the critical flux for the moon to experience an RG effect. Irradiation from a 13-Jupiter-mass planet onto an Earth-sized moon at a distance of ten Jupiter radii can drive an RG state on the moon for about 200 Myr. If stellar illumination equivalent to that received by Earth from the Sun is added, the RG holds for about 500 Myr. After 1000 Myr, the planet's habitable edge has moved inward to about 6 planetary radii. Exomoons in orbits with eccentricities of 0.1 experience strong tidal heating; they must orbit a 13-Jupiter-mass host beyond 29 or 18 Jupiter radii after 100 Myr (at the inner and outer boundaries of the stellar HZ, respectively), and beyond 13 Jupiter radii (in both cases) after 1000 Myr to be habitable. If a roughly Earth-sized moon would be detected in orbit around a giant planet, and if this binary would orbit in the stellar HZ, then it will be crucial to recover the moon's orbital history. If such a moon around a 13-Jupiter-mass planet has been closer than 20 Jupiter radii to its host during the first few hundred million years at least, then it might have lost substantial amounts of its initial water reservoir and be uninhabitable today.

The detection of moons orbiting extrasolar planets ("exomoons") has now become feasible. Once they are discovered in the circumstellar habitable zone, questions about their habitability will emerge. Exomoons are likely to be tidally locked to their planet and hence experience days much shorter than their orbital period around the star and have seasons, all of which works in favor of habitability. These satellites can receive more illumination per area than their host planets, as the planet reflects stellar light and emits thermal photons. On the contrary, eclipses can significantly alter local climates on exomoons by reducing stellar illumination. In addition to radiative heating, tidal heating can be very large on exomoons, possibly even large enough for sterilization. We identify combinations of physical and orbital parameters for which radiative and tidal heating are strong enough to trigger a runaway greenhouse. By analogy with the circumstellar habitable zone, these constraints define a circumplanetary "habitable edge". We apply our model to hypothetical moons around the recently discovered exoplanet Kepler-22b and the giant planet candidate KOI211.01 and describe, for the first time, the orbits of habitable exomoons. If either planet hosted a satellite at a distance greater than 10 planetary radii, then this could indicate the presence of a habitable moon.

In a multiplanet system, tides acting on the inner planet can significantly affect the orbital evolution of the entire system. While tides usually damp eccentricities, a novel mechanism identified by Correia et al. (2012) tends to raise eccentricities as a result of the tides' effect on the inner planet's rotation. Our analytical description of this spin-driven tidal (SDT) effect shows that, while the inner planet's eccentricity undergoes pumping, the process is more completely described by an exchange of strength between the two eigenmodes of the dynamical system. Our analysis allows derivation of criteria for two-planet coplanar systems where the SDT effect can reverse tidal damping, and may preclude the effect's being significant for realistic systems. For the specific case quantified by Correia et al., the effect is strong because of the large adopted tidal time lag, which may not be appropriate for the assumed Saturn-like inner planet. On the other hand, the effective Q for any given planet in exotic circumstances is very uncertain, so the SDT effect could play a role in planetary evolution.

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.

Since low-mass stars have low luminosities, orbits at which liquid water can exist on Earth-sized planets are relatively close-in, which produces Doppler signals that are detectable using state-of-the-art Doppler spectroscopy. GJ 667C is already known to be orbited by two super-Earth candidates. We investigate whether the data supports the presence of additional companions. We obtain new Doppler measurements from HARPS extracted spectra and combined them with those obtained from the PFS and HIRES spectrographs. We used Bayesian and periodogram-based methods to re-assess the number of candidates and evaluated the confidence of each detection. Among other tests, we validated the planet candidates by analyzing correlations of each Doppler signal activity indices and investigate quasi-periodicity. Doppler measurements of GJ 667C are described better by six Keplerian-like signals: the two known candidates (b and c); three additional few-Earth mass candidates with periods of 92, 62 and 39 days (d, e and f); a cold super-Earth in a 260-day orbit (g) and tantalizing evidence of a $\sim$ 1 M$_\oplus$ object in a close-in orbit of 17 days (h). We explore whether long-term stable orbits are compatible with the data by integrating 8$\times 10^4$ solutions derived from the Bayesian samplings. The system consisting of six planets is compatible with dynamically stable configurations. As for the solar system, the most stable solutions do not contain mean-motion resonances and are described well by analytic Laplace-Lagrange solutions. The presence of a seventh planet (h) is supported by the fact that it appears squarely centered on the only island of stability left in the six-planet solution. Habitability assessments accounting for the stellar flux, as well as tidal dissipation effects, indicate that three (maybe four) planets are potentially habitable...

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of three years. We present measurements which confirm previous estimates of system parameters like the normalized semi-major axis (a/R_{*}), stellar density (\rho_{*}), impact parameter (b) and orbital inclination (i_{orb}). Our errors on system parameters like a/R_{*} and \rho_{*} have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations and we are able to rule out co-planar companions with masses \ge 0.20 \mearth\ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 sec) measured for XO-2b, and hence remain undetectable.

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 +/- 2.0 M_Jup to GSC 03546-01452 from the MARVELS radial velocity survey, which we designate as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS-6b has a mass below the hydrogen-burning limit of 0.072 M_Sun, and thus it is a high-confidence BD companion. It has a moderately long orbital period of 47.8929 +0.0063/-0.0062 days with a low eccentricty of 0.1442 +0.0078/-0.0073, and a semi-amplitude of 1644 +12/-13 m/s. Moderate resolution spectroscopy of the host star has determined the following parameters: T_eff = 5598 +/- 63, log g = 4.44 +/- 0.17, and [Fe/H] = +0.40 +/- 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of M_star = 1.11 +/- 0.11 M_Sun and R_star = 1.06 +/- 0.23 R_Sun with an age consistent with less than ~6 Gyr at a distance of 219 +/- 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot definitively rule out a transiting configuration based on our observations. There is a visual companion detected with Lucky Imaging at 7.7 arcsec from the host star, but our analysis shows that it is not bound to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for both stars and planets, making this a rare object even compared to other BDs.

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) program. Model fitting to these data allows us to improve measurements of the hot-Jupiter exoplanet WASP-2b and its orbital parameters by a factor of ~2 over prior studies; we do not find evidence for transit depth variations. We do find reduced chi^2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration cycle (e.g. continuous observations from space-based platforms). However, even in systems observed with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot-Jupiter systems. This result is more meaningful, and indicates that hot-Jupiter systems, while they are easily observable from the ground, do not appear to be currently configured in a manner favorable to the detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, and/or implementation at extreme quality observing sites to minimize atmospheric red noise.

The exquisite photometric precision of the Kepler space telescope now puts the detection of extrasolar moons at the horizon. Here, we firstly review observational and analytical techniques that have recently been proposed to find exomoons. Secondly, we discuss the prospects of characterizing potentially habitable extrasolar satellites. With moons being much more numerous than planets in the solar system and with most exoplanets found in the stellar habitable zone being gas giants, habitable moons could be as abundant as habitable planets. However, satellites orbiting planets in the habitable zones of cool stars will encounter strong tidal heating and likely appear as hot moons.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed eleven transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrESb are consistent with previous estimates to within the 2\sigma\ confidence level. We improved the errors (by 10--30%) on system parameters like the orbital inclination ($i_{\text{orb}}$), impact parameter (b) and stellar density (\rho$_{\star}$) compared to previous measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data show no evidence for transit timing variations and our timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b.

Traditionally stellar radiation has been the only heat source considered capable of determining global climate on long timescales. Here we show that terrestrial exoplanets orbiting low-mass stars may be tidally heated at high enough levels to induce a runaway greenhouse for a long enough duration for all the hydrogen to escape. Without hydrogen, the planet no longer has water and cannot support life. We call these planets "Tidal Venuses," and the phenomenon a "tidal greenhouse." Tidal effects also circularize the orbit, which decreases tidal heating. Hence, some planets may form with large eccentricity, with its accompanying large tidal heating, and lose their water, but eventually settle into nearly circular orbits (i.e. with negligible tidal heating) in the habitable zone (HZ). However, these planets are not habitable as past tidal heating desiccated them, and hence should not be ranked highly for detailed follow-up observations aimed at detecting biosignatures. Planets orbiting stars with masses <0.3 solar masses may be in danger of desiccation via tidal heating. We apply these concepts to Gl 667C c, a ~4.5 Earth-mass planet orbiting a 0.3 solar mass star at 0.12 AU. We find that it probably did not lose its water via tidal heating as orbital stability is unlikely for the high eccentricities required for the tidal greenhouse. As the inner edge of the HZ is defined by the onset of a runaway or moist greenhouse powered by radiation, our results represent a fundamental revision to the HZ for non-circular orbits. In the appendices we review a) the moist and runaway greenhouses, b) hydrogen escape, c) stellar mass-radius and mass-luminosity relations, d) terrestrial planet mass-radius relations, and e) linear tidal theories. [abridged]

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

White and brown dwarfs are astrophysical objects that are bright enough to support an insolation habitable zone (IHZ). Unlike hydrogen-burning stars, they cool and become less luminous with time, and hence their IHZ moves in with time. The inner edge of the IHZ is defined as the orbital radius at which a planet may enter a moist or runaway greenhouse, phenomena that can remove a planet's surface water forever. Thus, as the IHZ moves in, planets that enter it may no longer have any water, and are still uninhabitable. Additionally, the close proximity of the IHZ to the primary leads to concern that tidal heating may also be strong enough to trigger a runaway greenhouse, even for orbital eccentricities as small as 10^-6. Water loss occurs due to photolyzation by UV photons in the planetary stratosphere, followed by hydrogen escape. Young white dwarfs emit a large amount of these photons as their surface temperatures are over 10^4 K. The situation is less clear for brown dwarfs, as observational data do not constrain their early activity and UV emission very well. Nonetheless, both types of planets are at risk of never achieving habitable conditions, but planets orbiting white dwarfs may be less likely to sustain life than those orbiting brown dwarfs. We consider the future habitability of the planet candidates KOI 55.01 and 55.02 in these terms and find they are unlikely to become habitable.

Detections of massive extrasolar moons are shown feasible with the Kepler space telescope. Kepler's findings of about 50 exoplanets in the stellar habitable zone naturally make us wonder about the habitability of their hypothetical moons. Illumination from the planet, eclipses, tidal heating, and tidal locking distinguish remote characterization of exomoons from that of exoplanets. We show how evaluation of an exomoon's habitability is possible based on the parameters accessible by current and near-future technology.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 \sigma, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~<5 Gyr) solar-like star having a mass of 1.07 +/- 0.08 MSun and radius of 0.99 +/- 0.18 RSun. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG telescope obtained over a period of ~2 years. The best Keplerian orbital fit parameters were found to have a period of 78.994 +/- 0.012 days, an eccentricity of 0.1095 +/- 0.0023, and a semi-amplitude of 4199 +/- 11 m/s. We determine the minimum companion mass (if sin i = 1) to be 97.7 +/- 5.8 MJup. The system's companion to host star mass ratio, >0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

Parenago 1802, a member of the ~1 Myr Orion Nebula Cluster, is a double-lined, detached eclipsing binary in a 4.674 d orbit, with equal-mass components (M_2/M_1 = 0.985 \pm 0.029). Here we present extensive VIcJHKs light curves spanning ~15 yr, as well as a Keck/HIRES optical spectrum. The light curves evince a third light source that is variable with a period of 0.73 d, and is also manifested in the high-resolution spectrum, strongly indicating the presence of a third star in the system, probably a rapidly rotating classical T Tauri star. We incorporate this third light into our radial velocity and light curve modeling of the eclipsing pair, measuring accurate masses (M_1 = 0.391 \pm 0.032, M_2 = 0.385 \pm 0.032 M\odot), radii (R_1 = 1.73 \pm 0.02, R_2 = 1.62 \pm 0.02 R\odot), and temperature ratio (T_1/T_2 = 1.0924 \pm 0.0017). Thus the radii of the eclipsing stars differ by 6.9 \pm 0.8%, the temperatures differ by 9.2 \pm 0.2%, and consequently the luminosities differ by 62 \pm 3%, despite having masses equal to within 3%. This could be indicative of an age difference of ~3x10^5 yr between the two eclipsing stars, perhaps a vestige of the binary formation history. We find that the eclipsing pair is in an orbit that has not yet fully circularized, e = 0.0166 \pm 0.003. In addition, we measure the rotation rate of the eclipsing stars to be 4.629 \pm 0.006 d; they rotate slightly faster than their 4.674 d orbit. The non-zero eccentricity and super-synchronous rotation suggest that the eclipsing pair should be tidally interacting, so we calculate the tidal history of the system according to different tidal evolution theories. We find that tidal heating effects can explain the observed luminosity difference of the eclipsing pair, providing an alternative to the previously suggested age difference.

Stellar radiation has conservatively been used as the key constraint to planetary habitability. We review here the effects of tides, exerted by the host star on the planet, on the evolution of the planetary spin. Tides initially drive the rotation period and the orientation of the rotation axis into an equilibrium state but do not necessarily lead to synchronous rotation. As tides also circularize the orbit, eventually the rotation period does equal the orbital period and one hemisphere will be permanently irradiated by the star. Furthermore, the rotational axis will become perpendicular to the orbit, i.e. the planetary surface will not experience seasonal variations of the insolation. We illustrate here how tides alter the spins of planets in the traditional habitable zone. As an example, we show that, neglecting perturbations due to other companions, the Super-Earth Gl581d performs two rotations per orbit and that any primordial obliquity has been eroded.

Stellar insolation has been used as the main constraint on a planet's habitability. However, as more Earth-like planets are discovered around low-mass stars (LMSs), a re-examination of the role of tides on the habitability of exoplanets has begun. Those studies have yet to consider the misalignment between a planet's rotational axis and the orbital plane normal, i.e. the planetary obliquity. We apply two equilibrium tide theories to compute the obliquity evolution of terrestrial planets orbiting in the habitable zones around LMSs. The time for the obliquity to decrease from an Earth-like obliquity of 23.5 deg to 5 deg, the 'tilt erosion time', is compared to the traditional insolation habitable zone (IHZ) as a function of semi-major axis, eccentricity, and stellar mass. We also compute tidal heating and equilibrium rotation caused by obliquity tides. The Super-Earth Gl581d and the planet candidate Gl581g are studied as examples for tidal processes. Earth-like obliquities of terrestrial planets in the IHZ around <0.25 solar mass (M_sun) stars are eroded in <0.1 Gyr. Only terrestrial planets orbiting stars with masses >0.9 M_sun experience tilt erosion times larger than 1 Gyr throughout the IHZ. Terrestrial planets in the IHZ of stars with masses <0.25 M_sun undergo significant tidal heating due to obliquity tides. The predictions of the two tidal models diverge significantly for e>0.3. In our two-body simulations, Gl581d's obliquity is eroded to 0 and its rotation period reached its equilibrium state of half its orbital period in <0.1 Gyr. Tidal surface heating on the putative Gl581g is <150 mW/m^2 as long as its eccentricity is <0.3. Obliquity tides modify the concept of the habitable zone. Tilt erosion of terrestrial planets orbiting LMSs should be included by atmospheric modelers. Tidal heating needs to be considered by geologists.

The search for rocky exoplanets plays an important role in our quest for extra-terrestrial life. Here, we discuss the extreme physical properties possible for the first characterized rocky super-Earth, CoRoT-7b (R_pl = 1.58 \pm 0.10 R_Earth, Mpl = 6.9 \pm 1.2 M_Earth). It is extremely close to its star (a = 0.0171 AU = 4.48 R_st), with its spin and orbital rotation likely synchronized. The comparison of its location in the (Mpl, Rpl) plane with the predictions of planetary models for different compositions points to an Earth-like composition, even if the error bars of the measured quantities and the partial degeneracy of the models prevent a definitive conclusion. The proximity to its star provides an additional constraint on the model. It implies a high extreme-UV flux and particle wind, and the corresponding efficient erosion of the planetary atmosphere especially for volatile species including water. Consequently, we make the working hypothesis that the planet is rocky with no volatiles in its atmosphere, and derive the physical properties that result. As a consequence, the atmosphere is made of rocky vapours with a very low pressure (P \leq 1.5 Pa), no cloud can be sustained, and no thermalisation of the planetary is expected. The dayside is very hot (2474 \leq 71 K at the sub-stellar point) while the nightside is very cold (50 to 75 K). The sub-stellar point is as hot as the tungsten filament of an incandescent bulb, resulting in the melting and distillation of silicate rocks and the formation of a lava ocean. These possible features of CoRoT-7b could be common to many small and hot planets, including the recently discovered Kepler-10b. They define a new class of objects that we propose to name "Lava-ocean planets".

Terrestrial planets are more likely to be detected if they orbit M dwarfs due to the favorable planet/star size and mass ratios. However, M dwarf habitable zones are significantly closer to the star than the one around our Sun, which leads to different requirements for planetary habitability and its detection. We review 1) the current limits to detection, 2) the role of M dwarf spectral energy distributions on atmospheric chemistry, 3) tidal effects, stressing that tidal locking is not synonymous with synchronous rotation, 4) the role of atmospheric mass loss and propose that some habitable worlds may be the volatile-rich, evaporated cores of giant planets, and 5) the role of planetary rotation and magnetic field generation, emphasizing that slow rotation does not preclude strong magnetic fields and their shielding of the surface from stellar activity. Finally we present preliminary findings of the NASA Astrobiology Institute's workshop "Revisiting the Habitable Zone." We assess the recently-announced planet Gl 581 g and find no obvious barriers to habitability. We conclude that no known phenomenon completely precludes the habitability of terrestrial planets orbiting cool stars.

We present three transits of GJ 1214b, observed as part of the Apache Point Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE). We used APOSTLE r-band lightcurves in conjunction with previously gathered data of GJ 1214b to re-derive system parameters. By using parameters such as transit duration and ingress/egress length we are able to reduce the degeneracies between parameters in the fitted transit model, which is a preferred condition for Markov Chain Monte Carlo techniques typically used to quantify uncertainties in measured parameters. The joint analysis of this multi-wavelength dataset confirms earlier estimates of system parameters including planetary orbital period, the planet-to-star radius ratio and stellar density. We fit the photometric spectralenergy distribution of GJ 1214 to derive stellar luminosity, which we then use to derive its absolute mass and radius. From these derived stellar properties and previously published radial velocity data we were able to refine estimates of the absolute parameters for the planet GJ 1214b. Transit times derived from our study show no evidence for strong transit timing variations. Some lightcurves we present show features that we believe are due to stellar activity. During the first night we observed a rise in the out-of-eclipse flux of GJ 1214 with a characteristic fast-rise exponential decay shape commonly associated with stellar flares. On the second night we observed a minor brightening during transit, which we believe might have been caused by the planet obscuring a star-spot on the stellar disk.

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

We evaluate the orbital evolution and several plausible origins scenarios for the mutually inclined orbits of Upsilon Andromedae c and d. These two planets have orbital elements that oscillate with large amplitudes and lie close to the stability boundary. This configuration, and in particular the observed mutual inclination, demands an explanation. The planetary system may be influenced by a nearby low-mass star, Upsilon And B, which could perturb the planetary orbits, but we find it cannot modify two coplanar orbits into the observed mutual inclination of ~30 deg. However, it could incite ejections or collisions between planetary companions that subsequently raise the mutual inclination to >30 deg. Our simulated systems with large mutual inclinations tend to be further from the stability boundary than Upsilon And, but we are able to produce similar systems. We conclude that scattering is a plausible mechanism to explain the observed orbits of Upsilon And c and d, but we cannot determine whether the scattering was caused by instabilities among the planets themselves or by perturbations from Upsilon And B. We also develop a procedure to quantitatively compare numerous properties of the observed system to our numerical models. Although we only implement this procedure to Upsilon And, it may be applied to any exoplanetary system.

Tides may be crucial to the habitability of exoplanets. If such planets form around low-mass stars, then those in the circumstellar habitable zone will be close enough to their host stars to experience strong tidal forces. Tides may result in orbital decay and circularization, evolution toward zero obliquity, a fixed rotation rate (not necessarily synchronous), and substantial internal heating. Due to tidal effects, the range of habitable orbital locations may be quite different from that defined by the traditional concept of a habitable zone (HZ) based on stellar insolation, atmospheric effects, and liquid water on a planet's surface. Tidal heating may make locations within the traditional HZ too hot, while planets outside the traditional zone could be rendered quite habitable due to tides. Here we consider these effects on the exoplanet GJ 581 d.

We report the discovery of a low-mass companion orbiting the metal-rich, main sequence F star TYC 2949-00557-1 during the MARVELS (Multi-object APO Radial Velocity Exoplanet Large-area Survey) Pilot Project. The host star has an effective temperature T_eff = 6135 +/- 40 K, log(g) = 4.4 +/- 0.1 and [Fe/H] = 0.32 +/- 0.01, indicating a mass of M = 1.25 +/- 0.09 M_\odot and R = 1.15 +/- 0.15 R_\odot. The companion has an orbital period of 5.69449 +/- 0.00023 days and straddles the hydrogen burning limit with a minimum mass of 64 M_J, and may thus be an example of the rare class of brown dwarfs orbiting at distances comparable to those of "Hot Jupiters." We present relative photometry that demonstrates the host star is photometrically stable at the few millimagnitude level on time scales of hours to years, and rules out transits for a companion of radius greater than 0.8 R_J at the 95% confidence level. Tidal analysis of the system suggests that the star and companion are likely in a double synchronous state where both rotational and orbital synchronization have been achieved. This is the first low-mass companion detected with a multi-object, dispersed, fixed-delay interferometer.

CoRoT-7 b is the first confirmed rocky exoplanet, but, with an orbital semi-major axis of 0.0172 AU, its origins may be unlike any rocky planet in our solar system. In this study, we consider the roles of tidal evolution and evaporative mass loss in CoRoT-7 b's history, which together have modified the planet's mass and orbit. If CoRoT-7 b has always been a rocky body, evaporation may have driven off almost half its original mass, but the mass loss may depend sensitively on the extent of tidal decay of its orbit. As tides caused CoRoT-7 b's orbit to decay, they brought the planet closer to its host star, thereby enhancing the mass loss rate. Such a large mass loss also suggests the possibility that CoRoT-7 b began as a gas giant planet and had its original atmosphere completely evaporated. In this case, we find that CoRoT-7 b's original mass probably didn't exceed 200 Earth masses (about 2/3 of a Jupiter mass). Tides raised on the host star by the planet may have significantly reduced the orbital semi-major axis, perhaps causing the planet to migrate through mean-motion resonances with the other planet in the system, CoRoT-7 c. The coupling between tidal evolution and mass loss may be important not only for CoRoT-7 b but also for other close-in exoplanets, and future studies of mass loss and orbital evolution may provide insight into the origin and fate of close-in planets, both rocky and gaseous.

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.

(Abridged) Hubble Space Telescope (HST) Fine Guidance Sensor astrometric observations of the G4 IV star HD 38529 are combined with the results of the analysis of extensive ground-based radial velocity data to determine the mass of the outermost of two previously known companions. Our new radial velocities obtained with the Hobby-Eberly Telescope and velocities from the Carnegie-California group now span over eleven years. With these data we obtain improved RV orbital elements for both the inner companion, HD 38529 b and the outer companion, HD 38529 c. We identify a rotational period of HD 38529 (P_{rot}=31.65 +/- 0.17 d) with FGS photometry. We model the combined astrometric and RV measurements to obtain the parallax, proper motion, perturbation period, perturbation inclination, and perturbation size due to HD 38529 c. For HD 38529 c we find P = 2136.1 +/- 0.3 d, perturbation semi-major axis \alpha =1.05 +/-0.06$ mas, and inclination $i$ = 48.3 deg +/- 4 deg. Assuming a primary mass M_* = 1.48 M_{sun}, we obtain a companion mass M_c = 17.6 ^{+1.5}_{-1.2} M_{Jup}, 3-sigma above a 13 M_{Jup} deuterium burning, brown dwarf lower limit. Dynamical simulations incorporating this accurate mass for HD 38529 c indicate that a near-Saturn mass planet could exist between the two known companions. We find weak evidence of an additional low amplitude signal that can be modeled as a planetary-mass (~0.17 M$_{Jup}) companion at P~194 days. Additional observations (radial velocities and/or Gaia astrometry) are required to validate an interpretation of HD 38529 d as a planetary-mass companion. If confirmed, the resulting HD 38529 planetary system may be an example of a "Packed Planetary System".

2MASSJ05352184-0546085 (2M0535-05) is the only known eclipsing brown dwarf (BD) binary, and so may serve as an important benchmark for models of BD formation and evolution. However, theoretical predictions of the system's properties seem inconsistent with observations: i. The more massive (primary) component is observed to be cooler than the less massive (secondary) one. ii. The secondary is more luminous (by roughly 10^{24} W) than expected. We study the impact of tidal heating to the energy budget of both components. We also compare various plausible tidal models to determine a range of predicted properties. We apply two versions of two different, well-known models for tidal interaction, respectively, (i.) the 'constant-phase-lag' model and (ii.) the 'constant-time-lag' model, and incorporate the predicted tidal heating into a model of BD structure. We find that the contribution of heat from tides in 2M0535-05 alone may only be large enough to account for the discrepancies between observation and theory in an unlikely region of the parameter space. The tidal quality factor of BDs, Q_{BD}, would have to be 10^{3.5} and the secondary needs a spin-orbit misalignment greater than 50 degrees. However, tidal synchronization time scales for 2M0535-05 restrict the tidal dissipation function Q_{BD} to values greater than 10^{4.5} and rule out intense tidal heating in 2M0535-05. We provide the first constraint on Q_{BD}. Tidal heating alone is unlikely to be responsible for the surprising temperature reversal within 2M0535-05. But an evolutionary embedment of tidal effects and a coupled treatment with the structural evolution of the BDs is necessary to corroborate or refute this result.

Hill's equations are an approximation that is useful in a number of areas of astrophysics including planetary rings and planetesimal disks. We derive a symplectic method for integrating Hill's equations based on a generalized leapfrog. This method is implemented in the parallel N-body code, PKDGRAV and tested on some simple orbits. The method demonstrates a lack of secular changes in orbital elements, making it a very useful technique for integrating Hill's equations over many dynamical times. Furthermore, the method allows for efficient collision searching using linear extrapolation of particle positions.

We review how tides may impact the habitability of terrestrial-like planets. If such planets form around low-mass stars, then planets in the circumstellar habitable zone will be close enough to their host stars to experience strong tidal forces. We discuss 1) decay of semi-major axis, 2) circularization of eccentric orbits, 3) evolution toward zero obliquity, 4) fixed rotation rates (not necessarily synchronous), and 5) internal heating. We briefly describe these effects using the example of a 0.25 solar mass star with a 10 Earth-mass companion. We suggest that the concept of a habitable zone should be modified to include the effects of tides.

CoRoT-7 b, a planet about 70% larger than the Earth orbiting a Sun-like star, is the first-discovered rocky exoplanet, and hence has been dubbed a "super-Earth". Some initial studies suggested that since the planet is so close to its host star, it receives enough insolation to partially melt its surface. However, these past studies failed to take into consideration the role that tides may play in this system. Even if the planet's eccentricity has always been zero, we show that tidal decay of semi-major axis could have been large enough that the planet formed on a wider orbit which received less insolation. Moreover, CoRoT-7 b could be tidally heated at a rate that dominates its geophysics and drives extreme volcanism. In this case, CoRoT-7 b is a "super-Io" that, like Jupiter's volcanic moon, is dominated by volcanism and rapid resurfacing. Such heating could occur with an eccentricity of just 10^-5. This small value could be driven by CoRoT-7 c if its own eccentricity is larger than ~10^-4. CoRoT-7 b may be the first of a class of planetary super-Ios likely to be revealed by the CoRoT and Kepler spacecraft.

The habitable zones of main sequence stars have traditionally been defined as the range of orbits that intercept the appropriate amount of stellar flux to permit surface water on a planet. Terrestrial exoplanets discovered to orbit M stars in these zones, which are close-in due to decreased stellar luminosity, may also undergo significant tidal heating. Tidal heating may span a wide range for terrestrial exoplanets and may significantly affect conditions near the surface. For example, if heating rates on an exoplanet are near or greater than that on Io (where tides drive volcanism that resurface the planet at least every 1 Myr) and produce similar surface conditions, then the development of life seems unlikely. On the other hand, if the tidal heating rate is less than the minimum to initiate plate tectonics, then CO_2 may not be recycled through subduction, leading to a runaway greenhouse that sterilizes the planet. These two cases represent potential boundaries to habitability and are presented along with the range of the traditional habitable zone for main sequence, low-mass stars. We propose a revised habitable zone that incorporates both stellar insolation and tidal heating. We apply these criteria to GJ 581 d and find that it is in the traditional habitable zone, but its tidal heating alone may be insufficient for plate tectonics.

The distribution of the orbits of close-in exoplanets shows evidence for on-going removal and destruction by tides. Tides raised on a planet's host star cause the planet's orbit to decay, even after the orbital eccentricity has dropped to zero. Comparison of the observed orbital distribution and predictions of tidal theory show good qualitative agreement, suggesting tidal destruction of close-in exoplanets is common. The process can explain the observed cut-off in small a-values, the clustering of orbital periods near three days, and the relative youth of transiting planets. Contrary to previous considerations, a mechanism to stop the inward migration of close-in planets at their current orbits is not necessarily required. Planets nearing tidal destruction may be found with extremely small a, possibly already stripped of any gaseous envelope. The recently discovered CoRoT-Exo-7 b may be an example of such a planet and will probably be destroyed by tides within the next few Gyrs. Also, where one or more planets have already been accreted, a star may exhibit an unusual composition and/or spin rate.

We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 solar mass star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 10^5 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation.

The known extrasolar multiple-planet systems share a surprising dynamical attribute: they cluster just beyond the Hill stability boundary. Here we show that the planet-planet scattering model, which naturally explains the observed exoplanet eccentricity distribution, can reproduce the observed distribution of dynamical configurations. We calculated how each of our scattered systems would appear over an appropriate range of viewing geometries; as Hill stability is weakly dependent on the masses, the mass-inclination degeneracy does not significantly affect our results. We consider a wide range of initial planetary mass distributions and find that some are poor fits to the observed systems. In fact, many of our scattering experiments overproduce systems very close to the stability boundary. The distribution of dynamical configurations of two-planet systems actually may provide better discrimination between scattering models than the distribution of eccentricity. Our results imply that, at least in their inner regions which are weakly affected by gas or planetesimal disks, planetary systems should be "packed", with no large gaps between planets.

We study the dynamical stability of an additional, potentially habitable planet in the HD 47186 planetary system. Two planets are currently known in this system: a "hot Neptune" with a period of 4.08 days and a Saturn-mass planet with a period of 3.7 years. Here we consider the possibility that one or more undetected planets exist between the two known planets and possibly within the habitable zone in this system. Given the relatively low masses of the known planets, additional planets could have masses $\lsim 10 \mearth$, and hence be terrestrial-like and further improving potential habitability. We perform $N$-body simulations to identify the stable zone between planets $b$ and $c$ and find that much of the inner habitable zone can harbor a $10 \mearth$ planet. With the current radial-velocity threshold of $\sim 1$ m/s, an additional planet should be detectable if it lies at the inner edge of the habitable zone at 0.8 AU. We also show that the stable zone could contain two additional planets of $10 \mearth$ each if their eccentricities are lower than $\sim 0.3$.

Three planets with minimum masses less than 10 Earth masses orbit the star HD 40307, suggesting these planets may be rocky. However, with only radial velocity data, it is impossible to determine if these planets are rocky or gaseous. Here we exploit various dynamical features of the system in order to assess the physical properties of the planets. Observations allow for circular orbits, but a numerical integration shows that the eccentricities must be at least 0.0001. Also, planets b and c are so close to the star that tidal effects are significant. If planet b has tidal parameters similar to the terrestrial planets in the Solar System and a remnant eccentricity larger than 0.001, then, going back in time, the system would have been unstable within the lifetime of the star (which we estimate to be 6.1 +/- 1.6 Gyr). Moreover, if the eccentricities are that large and the inner planet is rocky, then its tidal heating may be an order of magnitude greater than extremely volcanic Io, on a per unit surface area basis. If planet b is not terrestrial, e.g. Neptune-like, these physical constraints would not apply. This analysis suggests the planets are not terrestrial-like, and are more like our giant planets. In either case, we find that the planets probably formed at larger radii and migrated early-on (via disk interactions) into their current orbits. This study demonstrates how the orbital and dynamical properties of exoplanet systems may be used to constrain the planets' physical properties.

Planet-planet scattering is the leading mechanism to explain the large eccentricities of the observed exoplanet population. However, scattering has not been considered important to the production of pairs of planets in mean motion resonances (MMRs). We present results from a large number of numerical simulations of dynamical instabilities in 3-planet systems. We show that MMRs arise naturally in about five percent of cases. The most common resonances we populate are the 2:1 and 3:1 MMRs, although a wide variety of MMRs can occur, including high-order MMRs (up to eleventh order). MMRs are generated preferentially in systems with uneven mass distributions: the smallest planet is typically ejected after a series of close encounters, leaving the remaining, more massive planets in resonance. The distribution of resonant planets is consistent with the phase-space density of resonant orbits, meaning that planets are randomly thrown into MMRs rather than being slowly pulled into them. It may be possible to distinguish between MMRs created by scattering vs. convergent migration in a gaseous disk by considering planetary mass ratios: resonant pairs of planets beyond ~1 AU with more massive outer planets are likely to have formed by scattering. In addition, scattering may be responsible for pairs of planets in high-order MMRs (3:1 and higher) that are not easily populated by migration. The frequency of MMRs from scattering is comparable to the expected survival rate of MMRs in turbulent disks. Thus, planet-planet scattering is likely to be a major contributor to the population of resonant planets.

Transiting planets are generally close enough to their host stars that tides may govern their orbital and thermal evolution of these planets. We present calculations of the tidal evolution of recently discovered transiting planets and discuss their implications. The tidal heating that accompanies this orbital evolution can be so great that it controls the planet's physical properties and may explain the large radii observed in several cases, including, for example, TrES-4. Also because a planet's transit probability depends on its orbit, it evolves due to tides. Current values depend sensitively on the physical properties of the star and planet, as well as on the system's age. As a result, tidal effects may introduce observational biases in transit surveys, which may already be evident in current observations. Transiting planets tend to be younger than non-transiting planets, an indication that tidal evolution may have destroyed many close-in planets. Also the distribution of the masses of transiting planets may constrain the orbital inclinations of non-transiting planets.

Five planets are known to orbit the star 55 Cancri. The recently-discovered planet f at 0.78 AU (Fischer et al. 2008) is located at the inner edge of a previously-identified stable zone that separates the three close-in planets from planet d at 5.9 AU. Here we map the stability of the orbital space between planets f and d using a suite of n-body integrations that include an additional, yet-to-be-discovered planet g with a radial velocity amplitude of 5 m/s (planet mass = 0.5-1.2 Saturn masses). We find a large stable zone extending from 0.9 to 3.8 AU at eccentricities below 0.4. For each system we quantify the probability of detecting planets b-f on their current orbits given perturbations from hypothetical planet g, in order to further constrain the mass and orbit of an additional planet. We find that large perturbations are associated with specific mean motion resonances (MMRs) with planets f and d. We show that two MMRs, 3f:1g (the 1:3 MMR between planets g and f) and 4g:1d cannot contain a planet g. The 2f:1g MMR is unlikely to contain a planet more massive than about 20 Earth masses. The 3g:1d and 5g:2d MMRs could contain a resonant planet but the resonant location is strongly confined. The 3f:2g, 2g:1d and 3g:2d MMRs exert a stabilizing influence and could contain a resonant planet. Furthermore, we show that the stable zone may in fact contain 2-3 additional planets, if they are ~50 Earth masses each. Finally, we show that any planets exterior to planet d must reside beyond 10 AU.

The tidal heating of hypothetical rocky (or terrestrial) extra-solar planets spans a wide range of values depending on stellar masses and initial orbits. Tidal heating may be sufficiently large (in many cases, in excess of radiogenic heating) and long-lived to drive plate tectonics, similar to the Earth's, which may enhance the planet's habitability. In other cases, excessive tidal heating may result in Io-like planets with violent volcanism, probably rendering them unsuitable for life. On water-rich planets, tidal heating may generate sub-surface oceans analogous to Europa's with similar prospects for habitability. Tidal heating may enhance the outgassing of volatiles, contributing to the formation and replenishment of a planet's atmosphere. To address these issues, we model the tidal heating and evolution of hypothetical extra-solar terrestrial planets. The results presented here constrain the orbital and physical properties required for planets to be habitable.

Tides raised on a planet by its host star's gravity can reduce a planet's orbital semi-major axis and eccentricity. This effect is only relevant for planets orbiting very close to their host stars. The habitable zones of low-mass stars are also close-in and tides can alter the orbits of planets in these locations. We calculate the tidal evolution of hypothetical terrestrial planets around low-mass stars and show that tides can evolve planets past the inner edge of the habitable zone, sometimes in less than 1 billion years. This migration requires large eccentricities (>0.5) and low-mass stars (<0.35 M_Sun). Such migration may have important implications for the evolution of the atmosphere, internal heating and the Gaia hypothesis. Similarly, a planet detected interior to the habitable zone could have been habitable in the past. We consider the past habitability of the recently-discovered, ~5 M_Earth planet, Gliese 581 c. We find that it could have been habitable for reasonable choices of orbital and physical properties as recently as 2 Gyr ago. However, when we include constraints derived from the additional companions, we see that most parameter choices that predict past habitability require the two inner planets of the system to have crossed their mutual 3:1 mean motion resonance. As this crossing would likely have resulted in resonance capture, which is not observed, we conclude that Gl 581 c was probably never habitable.

Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that would place these systems near instability. The hypothetical planet in HD 74156 was expected to have a mass about equal to that of Saturn, a semi-major axis between 0.9 and 1.4 AU, and an eccentricity less than 0.2. HD 74156 d, a planet with a mass of 1.3 Saturn masses at 1.04 AU with an eccentricity of 0.25, was recently announced. We have reanalyzed all published data on this system in order to place tighter constraints on the properties of the new planet. We find two possible orbits for this planet, one close to that already identified and another (with a slightly better fit to the data) at ~0.89 AU. We also review the current status of other planet predictions, discuss the observed single planet systems, and suggest other systems which may contain planets in between those that are already known. The confirmation of the existence of HD 74156 d suggests that planet formation is an efficient process, and planetary systems should typically contain many planets.

Extra-solar planets close to their host stars have likely undergone significant tidal evolution since the time of their formation. Tides probably dominated their orbital evolution once the dust and gas had cleared away, and as the orbits evolved there was substantial tidal heating within the planets. The tidal heating history of each planet may have contributed significantly to the thermal budget that governed the planet's physical properties, including its radius, which in many cases may be measured by observing transit events. Typically, tidal heating increases as a planet moves inward toward its star and then decreases as its orbit circularizes. Here we compute the plausible heating histories for several planets with measured radii, using the same tidal parameters for the star and planet that had been shown to reconcile the eccentricity distribution of close-in planets with other extra-solar planets. Several planets are discussed, including for example HD 209458 b, which may have undergone substantial tidal heating during the past billion years, perhaps enough to explain its large measured radius. Our models also show that GJ 876 d may have experienced tremendous heating and is probably not a solid, rocky planet. Theoretical models should include the role of tidal heating, which is large, but time-varying.

The distribution of eccentricities e of extra-solar planets with semi-major axes a > 0.2 AU is very uniform, and values for e are generally large. For a < 0.2 AU, eccentricities are much smaller (most e < 0.2), a characteristic widely attributed to damping by tides after the planets formed and the protoplanetary gas disk dissipated. We have integrated the classical coupled tidal evolution equations for e and a backward in time over the estimated age of each planet, and confirmed that the distribution of initial e values of close-in planets matches that of the general population for reasonable tidal dissipation values Q, with the best fits for stellar and planetary Q being ~ 10^5.5 and ~ 10^6.5 respectively. The current small values of a were only reached gradually due to tides over the lifetimes of the planets, i.e., the earlier gas disk migration did not bring all planets to their current orbits. As the orbits tidally evolved, there was substantial tidal heating within the planets. The past tidal heating of each planet may have contributed significantly to the thermal budget that governed the planet's physical properties, including its radius, which in many cases may be measured by observing transit events. Here we also compute the plausible heating histories for a few planets with anomalously large measured radii, including HD 209458 b. We show that they may have undergone substantial tidal heating during the past billion years, perhaps enough to explain their large radii. Theoretical models of exoplanet interiors and the corresponding radii should include the role of large and time-variable tidal heating. Our results may have important implications for planet formation models, physical models of "hot Jupiters", and the success of transit surveys.

The dynamical interactions of planetary systems may be a clue to their formation histories. Therefore, the distribution of these interactions provides important constraints on models of planet formation. We focus on each system's apsidal motion and proximity to dynamical instability. Although only ~25 multiple planet systems have been discovered to date, our analyses in these terms have revealed several important features of planetary interactions. 1) Many systems interact such that they are near the boundary between stability and instability. 2) Planets tend to form such that at least one planet's eccentricity periodically drops to near zero. 3) Mean-motion resonant pairs would be unstable if not for the resonance. 4) Scattering of approximately equal mass planets is unlikely to produce the observed distribution of apsidal behavior. 5) Resonant interactions may be identified through calculating a system's proximity to instability, regardless of knowledge of angles such as mean longitude and longitude of periastron (e.g. GJ 317 b and c are probably in a 4:1 resonance). These properties of planetary systems have been identified through calculation of two parameters that describe the interaction. The apsidal interaction can be quantified by determining how close a planet is to an apsidal separatrix (a boundary between qualitatively different types of apsidal oscillations, e.g. libration or circulation of the major axes). The proximity to instability can be measured by comparing the observed orbital elements to an analytic boundary that describes a type of stability known as Hill stability. We have set up a website dedicated to presenting the most up-to-date information on dynamical interactions: http://www.lpl.arizona.edu/~rory/research/xsp/dynamics

The distribution of eccentricities e of extra-solar planets with semi-major axes a > 0.2 AU is very uniform, and values for e are relatively large, averaging 0.3 and broadly distributed up to near 1. For a < 0.2 AU, eccentricities are much smaller (most e < 0.2), a characteristic widely attributed to damping by tides after the planets formed and the protoplanetary gas disk dissipated. Most previous estimates of the tidal damping considered the tides raised on the planets, but ignored the tides raised on the stars. Most also assumed specific values for the planets' poorly constrained tidal dissipation parameter Qp. Perhaps most important, in many studies, the strongly coupled evolution between e and a was ignored. We have now integrated the coupled tidal evolution equations for e and a over the estimated age of each planet, and confirmed that the distribution of initial e values of close-in planets matches that of the general population for reasonable Q values, with the best fits for stellar and planetary Q being ~10^5.5 and ~10^6.5, respectively. The accompanying evolution of a values shows most close-in planets had significantly larger a at the start of tidal migration. The earlier gas disk migration did not bring all planets to their current orbits. The current small values of a were only reached gradually due to tides over the lifetimes of the planets. These results may have important implications for planet formation models, atmospheric models of "hot Jupiters", and the success of transit surveys.

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets (< 5-10 Earth masses). Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form "hot Earths" and their observable predictions. Models include: 1) in situ accretion; 2) formation at larger orbital distance followed by inward "type 1" migration; 3) formation from material being "shepherded" inward by a migrating gas giant planet; 4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disk; 5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances; and 6) photo-evaporative mass loss of a close-in giant planet. Models 1-4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets (> 5 Earth masses) with very close-in perihelion distances (< 0.025 AU), and even then the net inward movement in orbital distance is at most only 0.1-0.15 AU. For planets of less than about 70 Earth masses, photo-evaporation can remove the planet's envelope and leave behind the solid core on a Gyr timescale, but only for planets inside 0.025-0.05 AU. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit-timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).

The relationship between the boundaries for Hill and Lagrange stability in orbital element space is modified in the case of resonantly interacting planets. Hill stability requires the ordering of the planets to remain constant while Lagrange stability also requires all planets to remain bound to the central star. The Hill stability boundary is defined analytically, but no equations exist to define the Lagrange boundary, so we perform numerical experiments to estimate the location of this boundary. To explore the effect of resonances, we consider orbital element space near the conditions in the HD 82943 and 55 Cnc systems. Previous studies have shown that, for non-resonant systems, the two stability boundaries are nearly coincident. However the Hill stability formula are not applicable to resonant systems, and our investigation shows how the two boundaries diverge in the presence of a mean-motion resonance, while confirming that the Hill and Lagrange boundaries are similar otherwise. In resonance the region of stability is larger than the domain defined by the analytic formula for Hill stability. We find that nearly all known resonant interactions currently lie in this extra stable region, i.e. where the orbits would be unstable according to the non-resonant Hill stability formula. This result bears on the dynamical packing of planetary systems, showing how quantifying planetary systems' dynamical interactions (such as proximity to the Hill-stability boundary) provides new constraints on planet formation models.

Planets in extrasolar systems tend to interact such that their orbits lie near a boundary between apsidal libration and circulation, a "separatrix", with one eccentricity periodically reaching near-zero. One explanation, applied to the Upsilon Andromedae system, assumed three original planets on circular orbits. One is ejected, leaving the other two with near-separatrix behavior. We test that model by integrating hundreds of hypothetical, unstable planetary systems that eject a planet. We find that the probability that the remaining planets exhibit near-separatrix motion is small (< 5% compared with nearly 50% of observed systems). Moreover, while observed librating systems are evenly divided between aligned and anti-aligned pericenter longitudes, the scattering model strongly favors alignment. Alternative scattering theories are proposed, which may provide a more satisfactory fit with observed systems.

A widely considered characteristic of extra-solar planetary systems has been a seeming tendency for major axes of adjacent orbits to librate in stable configurations. Based on a new catalog of extra-solar planets (Butler et al. 2006) and our numerical integrations, we find that such small amplitude oscillations are actually not common, but in fact quite rare; most pairs of planets' major axes a re consistenet with circulating relative to one another. However, the new result s are consistent with studies that find that two-planet systems tend to lie near a separatrix between libration and circulation. Similarly, in systems of more than two planets, many adjacent orbits lie near a separatrix that divides modes of circulation.

Two types of stability boundaries exist for any planetary system consisting of one star and two planets. Lagrange stability requires that the planets remain bound to the star, conserves the ordering of the distance from the star, and limits the variations of orbital elements like semi-major axis and eccentricity. Hill stability only requires that the ordering of the planets remain constant; the outer planet may escape to infinity. A simple formula defines a region in orbital element space that is guaranteed to be Hill stable, although Hill stable orbits may lie outside the region as well. No analytic criteria describe Lagrange stability. We compare the results of 1000 numerical simulations of planetary systems similar to 47 UMa and HD 12661 with these two types of boundaries. All cases are consistent with the analytic criterion for Hill stability. Moreover, the numerically determined Lagrange boundary lies close to the analytic boundary for Hill stability. This result suggests an analytic formulation that may describe the criterion for Lagrange stability.

Recent results have shown that many of the known extrasolar planetary systems contain regions which are stable for both Earth-mass and Saturn-mass planets. Here we simulate the formation of terrestrial planets in four planetary systems -- 55 Cancri, HD 38529, HD 37124, and HD 74156 -- under the assumption that these systems of giant planets are complete and that their orbits are well-determined. Assuming the giant planets formed and migrated quickly, then terrestrial planets may form from a second generation of planetesimals. In each case, Moon- to Mars-sized planetary embryos are placed in between the giant planets and evolved for 100 Myr. We find that planets form relatively easily in 55 Cnc, with masses up to 0.6 Earth masses and in some cases substantial water contents and orbits in the habitable zone. HD 38529 is likely to support an asteroid belt but no terrestrial planets of significant mass. No terrestrial planets form in HD 37124 and HD 74156, although in some cases 1-2 lone embryos survive for 100 Myr. If migration occurred later, depleting the planetesimal disk, then massive terrestrial planets are unlikely to form in any of these systems.

Recent results have shown that many of the known extrasolar planetary systems contain regions which are stable for massless test particles. We examine the possibility that Saturn-mass planets exist in these systems, just below the detection threshold, and attempt to predict likely orbital parameters for such unseen planets. To do this, we insert a Saturn-mass planet into the stable regions of these systems and integrate its orbit for 100 million years. We conduct 200-600 of these experiments to test parameter space in HD37124, HD38529, 55Cnc, and HD74156. In HD37124 the global maximum of the survival rate of Saturns in parameter space is at semimajor axis a = 1.03 AU, eccentricity e=0.1. In HD38529, only 5% of Saturns are unstable, and the region in which a Saturn could survive is very broad, centered on 0.5<a<0.6, e<0.2. In 55Cnc we find three maxima at (a,e) = (1.0 AU, 0.02), (2.0 AU, 0.08), and (3.0 AU, 0.17). In HD74156 we find a broad maximum with $a$ = 0.9-1.2 AU, e<=0.15. Several of these maxima are located in the habitable zones of their parent stars and are therefore of astrobiological interest. We suggest the possibility that companions may lie in these locations of parameter space, and encourage further observational investigation of these systems.

Recent work has suggested that many planetary systems lie near instability. If all systems are near instability, an additional planet must exist in stable regions of well-separated extra-solar planetary systems to push these systems to the edge of stability. We examine the known systems by placing massless test particles in between the planets and integrating for 1-10 million years. We find that some systems, HD168443 and HD74156, eject nearly all test particles within 2 million years. However we find that HD37124, HD38529, and 55Cnc have large contiguous regions in which particles survive for 10 million years. These three systems, therefore, seem the most likely candidates for additional companions. Furthermore HD74156 and HD168443 must be complete and therefore radial velocity surveys should only focus on detecting more distant companions. We also find that several systems show stable regions that only exist at nonzero eccentricities.

We present results of numerical simulations which examine the dynamical stability of known planetary systems, a star with two or more planets. First we vary the initial conditions of each system based on observational data. We then determine regions of phase space which produce stable planetary configurations. For each system we perform 1000 ~1 million year integrations. We examine upsilon And, HD83443, GJ876, HD82943, 47UMa, HD168443, and the solar system (SS). We find that the resonant systems, 2 planets in a first order mean motion resonance, (HD82943 and GJ876) have very narrow zones of stability. The interacting systems, not in first order resonance, but able to perturb each other (upsilon And, 47UMa, and SS) have broad regions of stability. The separated systems, 2 planets beyond 10:1 resonance, (we only examine HD83443 and HD168443) are fully stable. Furthermore we find that the best fits to the interacting and resonant systems place them very close to unstable regions. The boundary in phase space between stability and instability depends strongly on the eccentricities, and (if applicable) the proximity of the system to perfect resonance. In addition to million year integrations, we also examined stability on ~100 million year timescales. For each system we ran ~10 long term simulations, and find that the Keplerian fits to these systems all contain configurations which may be regular on this timescale.

Because of the high eccentricities (~0.3) of two of the possible planets about the star Upsilon Andromeda, the stability of the system requires careful study. We present results of 1000 numerical simulations which explore the orbital parameter space as constrained by the observations. The orbital parameters of each planet are chosen from a Gaussian error distribution, and the resulting configuration is integrated for 1,000,000 years. We find that 84% of these integrations are stable. Configurations in which the eccentricity of the third planet is <0.3 are always stable, but when the eccentricity is >0.45 the system is always unstable, typically producing a close encounter between the second and third planets. A similar exercise with the gas giants in our own Solar System sampled with the same error distribution was performed. Approximately 81% of these simulations were stable for 10^6 years.