The Nysa-Polana complex is a group of low-inclination asteroid families in the inner main belt, bounded in semimajor axis by the Mars-crossing region and the Jupiter 3:1 mean-motion resonance. This group is important as the most likely source region for the target of the OSIRIS-REx mission, (101955) Bennu; however, family membership in the region is complicated by the presence of several dynamically overlapping families with a range of surface properties. The large S-type structure in the region appears to be associated with the parent body (135) Hertha, and displays an ($e_\text{P},a_\text{P}$) correlation consistent with a collision event near true anomaly of ~180 degrees with ejecta velocity $v_\text{ej} \sim 285$ m/s. The ejecta distribution from a collision with these orbital properties is predicted to have a maximum semimajor axis dispersion of $\delta a_{ej} = 0.005 \pm 0.008$ AU, which constitutes only a small fraction (7\%) of the observed semimajor axis dispersion, the rest of which is attributed to the Yarkovsky effect. The age of the family is inferred from the Yarkovsky dispersion to be $300^{+60}_{-50}$ My. Objects in a smaller cluster that overlaps the large Hertha family in proper orbital element space have reflectance properties more consistent with the X-type (135) Hertha than the surrounding S-type family. These objects form a distinct Yarkovsky "V" signature in ($a_\text{P},H$) space, consistent with a more recent collision, apparently also dynamically connected to (135) Hertha. The Nysa-Polana complex also contains a low-albedo family associated with (142) Polana (called "New Polana" by Walsh et al. 2013), and two other low-albedo families associated with (495) Eulalia. The second Eulalia family may be a high-$a_\text{P}$, low-$e_\text{P}$, low-$i_\text{P}$ component of the first Eulalia family-forming collision, possibly explained by an anisotropic ejection field.

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

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

The steady-state population of bodies resulting from a collisional cascade depends on how material strength varies with size. We find a simple expression for the power-law index of the population, given a power law that describes how material strength varies with size. This result is extended to the case relevant for the asteroid belt and Kuiper belt, in which the material strength is described by 2 separate power laws--one for small bodies and one for larger bodies. We find that the power-law index of the small body population is unaffected by the strength law for the large bodies, and vice versa. Simple analytical expressions describe a wave that is superimposed on the large body population because of the transition between the two power laws describing the strength. These analytical results yield excellent agreement with a numerical simulation of collisional evolution. These results will help to interpret observations of the asteroids and KBOs, and constrain the strength properties of those objects.

The Flora family resides in the densely populated inner main belt, bounded in semimajor axis by the $\nu_6$ secular resonance and the Jupiter 3:1 mean motion resonance. The presence of several large families that overlap dynamically with the Floras (e.g., the Vesta, Baptistina, and Nysa-Polana families), and the removal of a significant fraction of Floras via the nearby $\nu_6$ resonance complicates the Flora family's distinction in both proper orbital elements and reflectance properties. Here we use orbital information from the Asteroids Dynamic Site, color information from the Sloan Digital Sky Survey, and albedo information from the Wide-field Infrared Survey Explorer to obtain the median orbital and reflectance properties of the Floras by sampling the core of the family in multidimensional phase space. We find the median Flora SDSS colors to be $a^*$ = 0.126 $\pm$ 0.007 and $i-z = -0.037 \pm 0.007$; the median Flora albedo is $p_V$ = 0.291 $\pm$ 0.012. These properties allow us to define ranges for the Flora family in orbital and reflectance properties, as required for a detailed dynamical study. We use the young Karin family, for which we have an age determined via direct backward integration of members' orbits, to calibrate the Yarkovsky drift rates for the Flora family without having to estimate the Floras' material properties. The size-dependent dispersion of the Flora members in semimajor axis (the "V" plot) then yields an age for the family of $910^{+160}_{-120}$ My, with the uncertainty dominated by the uncertainty in the material properties of the family members (e.g., density and surface thermal properties). We discuss the effects on our age estimate of two independent processes that both introduce obliquity variations among the family members on short (My) timescales: 1) the capture of Flora members in spin-orbit resonance, and 2) YORP-driven obliquity variation.

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.

Classical secular theory can be a powerful tool to describe the qualitative character of multi-planet systems and offer insight into their histories. The eigenmodes of the secular behavior, rather than current orbital elements, can help identify tidal effects, early planet-planet scattering, and dynamical coupling among the planets, for systems in which mean-motion resonances do not play a role. Although tidal damping can result in aligned major axes after all but one eigenmode have damped away, such alignment may simply be fortuitous. An example of this is 55 Cancri (orbital solution of Fischer et al., 2008) where multiple eigenmodes remain undamped. Various solutions for 55 Cancri are compared, showing differing dynamical groupings, with implications for the coupling of eccentricities and for the partitioning of damping among the planets. Solutions for orbits that include expectations of past tidal evolution with observational data, must take into account which eigenmodes should be damped, rather than expecting particular eccentricities to be near zero. Classical secular theory is only accurate for low eccentricity values, but comparison with other results suggests that it can yield useful qualitative descriptions of behavior even for moderately large eccentricity values, and may have advantages for revealing underlying physical processes and, as large numbers of new systems are discovered, for triage to identify where more comprehensive dynamical studies should have priority.

Although warm jupiters are generally too far from their stars for tides to be important, the presence of an inner planetary companion to a warm jupiter can result in tidal evolution of the system. Insight into the process and its effects comes form classical secular theory of planetary perturbations. The lifetime of the inner planet may be shorter than the age of the system, because the warm jupiter maintains its eccentricity and hence promotes tidal migration into the star. Thus a warm jupiter observed to be alone in its system might have previously cleared away any interior planets. Before its demise, even if an inner planet is of terrestrial scale, it may promote damping of the warm jupiter's eccentricity. Thus any inferences of the initial orbit of an observed warm jupiter must include the possibility of a greater initial eccentricity than would be estimated by assuming it had always been alone. Tidal evolution involving multiple planets also enhances the internal heating of the planets, which readily exceeds that of stellar radiation for the inner planet, and may be great enough to affect the internal structure of warm jupiters. Secular theory gives insight into the tidal processes, providing, among other things, a way to constrain eccentricities of transiting planets based on estimates of the tidal parameter Q.

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.

Tidal damping of one of the orbits in a planetary system can lead to aligned major-axes (the so-called "fixed-point" condition), but currently aligned major axes do not necessarily imply such a history. An example is the nominal orbital solution for the 61 Virginis system where two orbits librate about alignment, but evaluation of the eigenmodes of the secular theory shows it could not be the result of tidal damping but rather of initial conditions. Nevertheless, the amplitudes of the eigenmodes suggest that this system may have undergone some degree of tidal damping.

In a multi-planet system, a gradual change in one planet's semi-major axis will affect the eccentricities of all the planets, as angular momentum is distributed via secular interactions. If tidal dissipation in the planet is the cause of the change in semi-major axis, it also damps that planet's eccentricity, which in turn also contributes to the evolution of all the eccentricities. Formulae quantifying the combined effects on the whole system due to semi-major axis changes, as well as eccentricity damping, are derived here for a two-planet system. The CoRoT 7 system is considered as an example.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.