We describe the current state of knowledge about Mercury's interior structure. We review the available observational constraints, including mass, size, density, gravity field, spin state, composition, and tidal response. These data enable the construction of models that represent the distribution of mass inside Mercury. In particular, we infer radial profiles of the pressure, density, and gravity in the core, mantle, and crust. We also examine Mercury's rotational dynamics and the influence of an inner core on the spin state and the determination of the moment of inertia. Finally, we discuss the wide-ranging implications of Mercury's internal structure on its thermal evolution, surface geology, capture in a unique spin-orbit resonance, and magnetic field generation.

The orbits of Pluto's four small satellites (Styx, Nix, Kerberos, and Hydra) are nearly circular and coplanar with the orbit of the large satellite Charon, with orbital periods nearly in the ratios 3:1, 4:1, 5:1, and 6:1 with Charon's orbital period. These properties suggest that the small satellites were created during the same impact event that placed Charon in orbit and had been pushed to their current positions by being locked in mean-motion resonances with Charon as Charon's orbit was expanded by tidal interactions with Pluto. Using the Pluto-Charon tidal evolution models developed by Cheng et al. (2014), we show that stable capture and transport of a test particle in multiple resonances at the same mean-motion commensurability is possible at the 5:1, 6:1, and 7:1 commensurabilities, if Pluto's zonal harmonic $J_{2P} = 0$. However, the test particle has significant orbital eccentricity at the end of the tidal evolution of Pluto-Charon in almost all cases, and there are no stable captures and transports at the 3:1 and 4:1 commensurabilities. Furthermore, a non-zero hydrostatic value of $J_{2P}$ destroys the conditions necessary for multiple resonance migration. Simulations with finite but minimal masses of Nix and Hydra also fail to yield any survivors. We conclude that the placing of the small satellites at their current orbital positions by resonant transport is extremely unlikely.

Both Pluto and its satellite Charon have rotation rates synchronous with their orbital mean motion. This is the theoretical end point of tidal evolution where transfer of angular momentum has ceased. Here we follow Pluto's tidal evolution from an initial state having the current total angular momentum of the system but with Charon in an eccentric orbit with semimajor axis $a \approx 4R_P$ (where $R_P$ is the radius of Pluto), consistent with its impact origin. Two tidal models are used, where the tidal dissipation function $Q \propto$ 1/frequency and $Q=$ constant, where details of the evolution are strongly model dependent. The inclusion of the gravitational harmonic coefficient $C_{22}$ of both bodies in the analysis allows smooth, self consistent evolution to the dual synchronous state, whereas its omission frustrates successful evolution in some cases. The zonal harmonic $J_2$ can also be included, but does not cause a significant effect on the overall evolution. The ratio of dissipation in Charon to that in Pluto controls the behavior of the orbital eccentricity, where a judicious choice leads to a nearly constant eccentricity until the final approach to dual synchronous rotation. The tidal models are complete in the sense that every nuance of tidal evolution is realized while conserving total angular momentum - including temporary capture into spin-orbit resonances as Charon's spin decreases and damped librations about the same.

The rotational evolution of Mercury's mantle and its core under conservative and dissipative torques is important for understanding the planet's spin state. Dissipation results from tides and viscous, magnetic and topographic core--mantle interactions. The dissipative core--mantle torques take the system to an equilibrium state wherein both spins are fixed in the frame precessing with the orbit, and in which the mantle and core are differentially rotating. This equilibrium exhibits a mantle spin axis that is offset from the Cassini state by larger amounts for weaker core--mantle coupling for all three dissipative core--mantle coupling mechanisms, and the spin axis of the core is separated farther from that of the mantle, leading to larger differential rotation. The relatively strong core--mantle coupling necessary to bring the mantle spin axis to its observed position close to the Cassini state is not obtained by any of the three dissipative core--mantle coupling mechanisms. For a hydrostatic ellipsoidal core--mantle boundary, pressure coupling dominates the dissipative effects on the mantle and core positions, and dissipation together with pressure coupling brings the mantle spin solidly to the Cassini state. The core spin goes to a position displaced from that of the mantle by about 3.55 arcmin nearly in the plane containing the Cassini state. With the maximum viscosity considered of $\nu\sim 15.0\,{\rm cm^2/s}$ if the coupling is by the circulation through an Ekman boundary layer or $\nu\sim 8.75\times 10^5\,{\rm cm^2/s}$ for purely viscous coupling, the core spin lags the precessing Cassini plane by 23 arcsec, whereas the mantle spin lags by only 0.055 arcsec. Larger, non hydrostatic values of the CMB ellipticity also result in the mantle spin at the Cassini state, but the core spin is moved closer to the mantle spin.

Recent discoveries of several transiting planets with clearly non-zero eccentricities and some large inclinations started changing the simple picture of close-in planets having circular and well-aligned orbits. Two major scenarios to form such planets are planet migration in a disk, and planet--planet interactions combined with tidal dissipation. The former scenario can naturally produce a circular and low-obliquity orbit, while the latter implicitly assumes an initially highly eccentric and possibly high-obliquity orbit, which are then circularized and aligned via tidal dissipation. We investigate the tidal evolution of transiting planets on eccentric orbits. We show that the current and future orbital evolution of these systems is likely dominated by tidal dissipation, and not by a more distant companion. Although most of these close-in planets experience orbital decay all the way to the Roche limit, there are two characteristic evolution paths for them, depending on the relative efficiency of tidal dissipation inside the star and the planet. We point out that the current observations may be consistent with one of them. Our results suggest that at least some of the close-in planets with non-zero orbital eccentricity may have been formed by tidally circularizing an initially eccentric orbit. We also find that even when the stellar spin-orbit misalignment is observed to be small at present, some systems could have had a highly misaligned orbit in the past. Finally, we also re-examine the recent claim by Levrard et. al., who found that all orbital and spin parameters evolve on a similar timescale to orbital decay.

This is a scientific strategy for the detection and characterization of extrasolar planets; that is, planets orbiting other stars. As such, it maps out over a 15-year horizon the techniques and capabilities required to detect and measure the properties of planets as small as Earth around stars as large as our own Sun. It shows how the technology pieces and their development fit together to achieve the primary goal of the strategy: if planets like Earth exist around stars within some tens of light years of our own Solar System, those planets will be found and their basic properties characterized. Essential to this strategy is not only the search for and examination of individual planets, but also a knowledge of the arrangement, or architecture, of planetary systems around as large a number of stars as possible; this is the second goal of the strategy. The final goal of the strategy is the study of disks around stars, important both to understand the implications of the variety of exoplanet systems for planet formation, and to determine how many nearby stars have environments around them clean enough of debris that planets may be sought and, if found, characterized.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

Tsiganis et al. (2005) have proposed that the current orbital architecture of the outer solar system could have been established if it was initially compact and Jupiter and Saturn crossed the 2:1 orbital resonance by divergent migration. The crossing led to close encounters among the giant planets, but the orbital eccentricities and inclinations were damped to their current values by interactions with planetesimals. Brunini (2006) has presented widely publicized numerical results showing that the close encounters led to the current obliquities of the giant planets. We present a simple analytic argument which shows that the change in the spin direction of a planet relative to an inertial frame during an encounter between the planets is very small and that the change in the obliquity (which is measured from the orbit normal) is due to the change in the orbital inclination. Since the inclinations are damped by planetesimal interactions on timescales much shorter than the timescales on which the spins precess due to the torques from the Sun, especially for Uranus and Neptune, the obliquities should return to small values if they are small before the encounters. We have performed simulations using the symplectic integrator SyMBA, modified to include spin evolution due to the torques from the Sun and mutual planetary interactions. Our numerical results are consistent with the analytic argument for no significant remnant obliquities.

(Abridged) The orbits of the recently discovered satellites of Pluto, S/2005 P2 and S/2005 P1, are significantly non-Keplerian, even if P2 and P1 have negligible masses, because the mass ratio of Charon-Pluto is ~0.1. We present an analytic theory with P2 and P1 treated as test particles. This analytic theory shows that the azimuthal periods of P2 and P1 are shorter than the Keplerian orbital periods and that the periapse and ascending node of each of the satellites precess at nearly equal rates in opposite directions. The deviation from Kepler's third law is already detected in the unperturbed Keplerian fit of Buie and coworkers. We also present direct numerical orbit integrations with different assumed masses for P2 and P1 within the ranges allowed by the albedo uncertainties. If the albedos are as high as that of Charon, the masses of P2 and P1 are sufficiently low that their orbits are well described by the analytic theory. There is at present no evidence that P2 has any significant epicyclic eccentricity. However, the orbit of P1 has a significant epicyclic eccentricity, and its prograde periapse precession with a period of 5300 days should be easily detectable. If the albedos are as low as that of comets, the large inferred masses induce significant variations in the epicyclic eccentricities and/or periapse longitudes on the 400-500-day timescales, due to the proximity of P2 and P1 to the 3:2 mean-motion commensurability. In fact, for the maximum inferred masses, P2 and P1 may be in the 3:2 mean-motion resonance, with the resonance variable involving the periapse longitude of P1 librating. Observations that sample the orbits of P2 and P1 well on the 400-500-day timescales should provide strong constraints on the masses of P2 and P1 in the near future.

In determining Mercury's core structure from its rotational properties, the value of the normalized moment of inertia, $C/MR^2$, from the location of Cassini 1 is crucial. If Mercury's spin axis occupies Cassini state 1, its position defines the location of the state. The spin might be displaced from the Cassini state if the spin is unable to follow the changes in the state position induced by the variations in the orbital parameters and the geometry of the solar system. The spin axis is expected to follow the Cassini state for orbit variations with time scales long compared to the 1000 year precession period of the spin about the Cassini state because the solid angle swept out by the spin axis as it precesses is an adiabatic invariant. Short period variations in the orbital elements of small amplitude should cause displacements that are commensurate with the amplitudes of the short period terms. By following simultaneously the spin position and the Cassini state position during long time scale orbital variations over past 3 million years (Quinn {\it et al.}, 1991) and short time scale variations from JPL Ephemeris DE 408 (Standish, 2005) we show that the spin axis will remain within one arcsec of the Cassini state after it is brought there by dissipative torques. We thus expect Mercury's spin to occupy Cassini state 1 well within the uncertainties for both radar and spacecraft measurements, with correspondingly tight constraints on $C/MR^2$.

An analysis based on the direct torque equations including tidal dissipation and a viscous core-mantle coupling is used to determine the damping time scales of O(10^5) years for free precession of the spin about the Cassini state and free libration in longitude for Mercury. The core-mantle coupling dominates the damping over the tides by one to two orders of magnitude for the plausible parameters chosen. The short damping times compared with the age of the solar system means we must find recent or on-going excitation mechanisms if such free motions are found by the current radar experiments or the future measurement by the MESSENGER and BepiColombo spacecraft that will orbit Mercury. We also show that the average precession rate is increased by about 30% over that obtained from the traditional precession constant because of a spin-orbit resonance induced contribution by the C_{22} term in the expansion of the gravitational field. The C_{22} contribution also causes the path of the spin during the precession to be slightly elliptical with a variation in the precession rate that is a maximum when the obliquity is a minimum. An observable free precession will compromise the determination of obliquity of the Cassini state and hence of C/MR^2 for Mercury, but a detected free libration will not compromise the determination of the forced libration amplitude and thus the verification of a liquid core

The two planets about the star GJ 876 appear to have undergone extensive migration from their point of origin in the protoplanetary disk -- both because of their close proximity to the star (30 and 60 day orbital periods) and because of their occupying three stable orbital resonances at the 2:1 mean-motion commensurability. The resonances were most likely established by converging differential migration of the planets leading to capture into the resonances. A problem with this scenario is that continued migration of the system while it is trapped in the resonances leads to orbital eccentricities that rapidly exceed the observational upper limits of e_1 = 0.31 and e_2 = 0.05. As seen in forced 3-body simulations, lower eccentricities would persist during migration only for an applied eccentricity damping. Here we explore the evolution of the GJ 876 system using two-dimensional hydrodynamical simulations that include viscous heating and radiative effects. We find that a hydrodynamic evolution within the resonance, where only the outer planet interacts with the disk, always rapidly leads to large values of eccentricities that exceed those observed. Only if mass is removed from the disk on a time scale of the order of the migration time scale (before there has been extensive migration after capture), as might occur for photoevaporation in the late phases of planet formation, can we end up with eccentricities that are consistent with the observations.

The Microlensing Planet Finder (MPF) is a proposed Discovery mission that will complete the first census of extrasolar planets with sensitivity to planets like those in our own solar system. MPF will employ a 1.1m aperture telescope, which images a 1.3 sq. deg. field-of-view in the near-IR, in order to detect extrasolar planets with the gravitational microlensing effect. MPF's sensitivity extends down to planets of 0.1 Earth masses, and MPF can detect Earth-like planets at all separations from 0.7AU to infinity. MPF's extrasolar planet census will provide critical information needed to understand the formation and frequency of extra solar planetary systems similar to our own.

(Abridged) We investigate the dynamical evolution of coplanar hierarchical two-planet systems where the ratio of the orbital semimajor axes alpha=a_1/a_2 is small. The orbital parameters obtained from a multiple Kepler fit to the radial velocity variations of a star are best interpreted as Jacobi coordinates and Jacobi coordinates should be used in any analyses of hierarchical planetary systems. An approximate theory that can be applied to coplanar hierarchical two-planet systems with a wide range of masses m_j and orbital eccentricities e_j is the octupole-level secular perturbation theory (OSPT). The OSPT shows that if the ratio of the maximum orbital angular momenta, lambda \approx (m_1/m_2) alpha^{1/2}, for given a_j is approximately equal to a critical value lambda_{crit}, then libration of the difference in the longitudes of periapse, w_1-w_2, about either 0 or 180 deg. is almost certain, with possibly large amplitude variations of both e_j. We establish that the OSPT is highly accurate for systems with alpha<0.1 and reasonably accurate even for systems with alpha as large as 1/3, provided that alpha is not too close to a significant mean-motion commensurability or above the stability boundary. The HD 168443 system is not in a secular resonance and its w_1-w_2 circulates. The HD 12661 system is the first extrasolar planetary system found to have w_1-w_2 librating about 180 deg. The libration of w_1-w_2 and the large-amplitude variations of both e_j in the HD 12661 system are consistent with the analytic results on systems with lambda \approx lambda_{crit}. The HD 12661 system with the best- fit orbital parameters and sin i = 1 is affected by the close proximity to the 11:2 commensurability, but small changes in the outer orbital period can result in configurations that are not affected by mean-motion commensurabilities.

Understanding the origin of the orbital resonances of the Galilean satellites of Jupiter will constrain the longevity of the extensive volcanism on Io, may explain a liquid ocean on Europa, and may guide studies of the dissipative properties of stars and Jupiter-like planets. The differential migration of the newly formed Galilean satellites due to interactions with a circumjovian disk can lead to the primordial formation of the Laplace relation n_1 - 3 n_2 + 2 n_3 = 0, where the n_i are the mean orbital angular velocities of Io, Europa, and Ganymede, respectively. This contrasts with the formation of the resonances by differential expansion of the orbits from tidal torques from Jupiter.

The Galactic Exoplanet Survey Telescope (GEST) will observe a 2 square degree field in the Galactic bulge to search for extra-solar planets using a gravitational lensing technique. This gravitational lensing technique is the only method employing currently available technology that can detect Earth-mass planets at high signal-to-noise, and can measure the frequency of terrestrial planets as a function of Galactic position. GEST's sensitivity extends down to the mass of Mars, and it can detect hundreds of terrestrial planets with semi-major axes ranging from 0.7 AU to infinity. GEST will be the first truly comprehensive survey of the Galaxy for planets like those in our own Solar System.

The 2:1 orbital resonances of the GJ 876 system can be easily established by the differential planet migration due to planet-nebula interaction. Significant eccentricity damping is required to produce the observed orbital eccentricities. The geometry of the GJ 876 resonance configuration differs from that of the Io-Europa pair, and this difference is due to the magnitudes of the eccentricities involved. We show that a large variation in the configuration of 2:1 and 3:1 resonances and, in particular, asymmetric librations can be expected among future discoveries.

(Abridged) A dynamical fit has placed the two planets about the star GJ 876 in coplanar orbits deep in 3 resonances at the 2:1 mean-motion commensurability with small libration amplitudes. The libration of both lowest order mean-motion resonance variables, theta_1 and theta_2, and the secular resonance variable, theta_3, about 0 deg. differs from the familiar geometry of the Io-Europa pair, where theta_2 and theta_3 librate about 180 deg. By considering a condition for stable simultaneous librations of theta_1 and theta_2, we show that the GJ 876 geometry results because of the large orbital eccentricities e_i, whereas the very small e_i in the Io-Europa system lead to the latter's geometry. Surprisingly, the GJ 876 resonance configuration remains stable for e_1 up to 0.86 and for amplitude of libration of theta_1 approaching 45 deg. with the current e_i. We find that inward migration of the outer planet of the GJ 876 system results in certain capture into the observed resonances if initially e_1 <0.06 and e_2<0.03 and the migration rate |(da_2/dt)/a_2| < 0.03(a_2/AU)^{-3/2} yr^{-1}. The bound on the migration rate is easily satisfied by migration due to planet-nebula interaction. If there is no eccentricity damping, eccentricity growth is rapid with continued migration within the resonance, with e_i exceeding the observed values after a further reduction in the semi-major axes a_i of only 7%. With eccentricity damping (de_i/dt)/e_i = -K|(da_i/dt)/a_i|, the e_i reach equilibrium values that remain constant for arbitrarily long migration within the resonances. The equilibrium e_i are close to the observed e_i for K=100 (K=10) if there is migration and damping of the outer planet only (of both planets). It is as yet unclear that planet-nebula interaction can produce the large value of K required to obtain the observed eccentricities.

The probability of detecting a planetary companion of a lensing star during a microlensing event toward the Galactic center, averaged over all relevant event and galactic parameters, when the planet-star mass ratio $q=0.001$ has a maximum exceeding 10% at an orbit semimajor axis $a$ near 1.5 AU for a uniform distribution of impact parameters. The maximum probability is raised to more than 20% for a distribution of source-lens impact parameters that is determined by the efficiency of event detection. The averaging procedures are carefully defined, and they determinine the dependence of the detection probabilities on several properties of the Galaxy. The probabilities scale approximately as $\sqrt{q}$. A planet is assumed detectable if the perturbation of the single lens light curve exceeds $2/(S/N)$ for at least 20 consecutive photometric points sometime during the event. Two meter telescopes with 60 second integrations in I-band with high time resolution photometry throughout the duration of an ongoing event are assumed. The probabilities are derived as a function of $a$, where they remain significant for $0.6<a<10$ AU. Dependence of the detection probabilities on the lens mass function, luminosity function of the source stars as modified by extinction, distribution of source-lens impact parameters, and the line of sight to the source are also determined, and the probabilities are averaged over the distribution of the projected planet position, the lens mass function, the distribution of impact parameters, the lens and source distances as weighted by their distributions along the line of sight and over the $I$-band apparent luminosity function of the sources. The extraction of the probabilility as a function of $a$ for a particular $q$ from empirical data is indicated.

The 2-Micron All Sky Survey (2MASS) (Skrutskie et al. 1997) and the DEep Near Infrared Survey of the southern sky (DENIS) (Epchtein et al. 1997) have revealed a heretofore unknown population of free brown dwarfs that has extended the local mass function down to as small as 0.01M_sun (Reid et al. 1999). If this local proportion of brown dwarfs extends throughout the Galaxy---in particular in the Galactic bulge---one expects an increase in the predicted fraction of short time scale microlensing events in directions toward the Galactic bulge. Zhao et al.(1996) have indicated that a mass function with 30-60% of the lens mass in brown dwarfs is not consistent with empirical microlensing data. Here we show that even the much lower mass fraction (~ 10%) of brown dwarfs inferred from the new discoveries appears inconsistent with the data. The added brown dwarfs do indeed increase the expected number of short time scale events, but they appear to drive the peak in the time scale frequency distribution to time scales smaller than that observed, and do not otherwise match the observed distribution. A reasonably good match to the empirical data (Alcock et al. 1996) is obtained by increasing the fraction of stars in the range 0.08<m<0.7M_sun considerably above that deduced from several star counts. However, all inferences from microlensing about the appropriate stellar mass function must be qualified by the meagerness of the microlensing data and the uncertainties in the Galactic model.

The dependence of microlensing time scale frequency distributions and optical depth toward the galactic center on galactic model parameters is explored in detail for a distribution of stars consisting of the Zhao (1996) bar and nucleus and the Bahcall and Soneira (1980) double exponential disk. The high sensitivity of these two microlensing measures to the circular velocity model, velocity dispersions, bulge mass, direction of the line of sight, bar axis orientation, star spatial distribution and the stellar mass function means no single galaxy property can be constrained very well without constraining most of the others. The model time scale frequency distributions are compared with that determined empirically by the MACHO group throughout. Although the MACHO empirical data are matched quite well with a nominal velocity model and with a mass function only of hydrogen burning stars that varies as m^{-2.2} to m^{-2.5} in the M star region, uncertainties in galactic structure, kinematics and content together with the paucity of published microlensing data preclude any claim of the model representing the real world. A variation of the mass function \sim m^{-1} in the M star region obtained from recent star counts, both local and in the galactic bulge, fails to yield a sufficient number of short time scale events compared to the MACHO data. The high sensitivity of the microlensing measures to the direction of the line of sight may mean that sufficient microlensing data is already in hand to constrain the bar distribution of stars. The procedure developed here for determining the time scale frequency distribution is particularly convenient for rapidly incorporating model changes as data from all sources continues to accumulate.

The statistical distribution of the masses of planets about stars between the Sun and the center of the galaxy is constrained to within a factor of three by an intensive search for planets during microlensing events. Projected separations in terms of the lens Einstein ring radius yield a rough estimate of the distribution of planetary semimajor axes with planetary mass. The search consists of following ongoing stellar microlensing events involving sources in the center of the galaxy lensed by intervening stars with high time resolution, 1% photometry in two colors in an attempt to catch any short time scale planetary perturbations of the otherwise smooth light curve. It is assumed that 3000 events are followed over an 8 year period, but with half of the lenses, those that are members of binary systems, devoid of planets. The remaining 1500 lenses have solar-system-like distributions of 4 or 5 planets. The expectations from the microlensing search are extremely assumption dependent with 56, 138, and 81 planets being detected for three sets of assumptions involving how the planetary masses and separations vary with lens mass. The events can be covered from 54% to 62% of the time on average by high time resolution photometry from a system of three or four dedicated two meter telescopes distributed in longitude, so 38% to 46% of the detectable small mass planets (very short perturbations of the light curve) will be missed. But perturbations comparable to a day in length means all of the detectable Jupiters and Saturns will in fact be detected as well as a large fraction of the Uranuses. Although meaningful statistics on planetary masses and separations can be inferred from such an intensive search, they, like the inferred data set, will be dependent on the assumed nature of the systems.