The Yarkovsky effect is a thermal process acting upon the orbits of small celestial bodies, which can cause these orbits to slowly expand or contract with time. The effect is subtle (da/dt ~ 10^-4 au/My for a 1 km diameter object) and is thus generally difficult to measure. We analyzed both optical and radar astrometry for 600 near-Earth asteroids (NEAs) for the purpose of detecting and quantifying the Yarkovsky effect. We present 247 NEAs with measured drift rates, which is the largest published set of Yarkovsky detections. This large sample size provides an opportunity to examine the Yarkovsky effect in a statistical manner. In particular, we describe two independent population-based tests that verify the measurement of Yarkovsky orbital drift. First, we provide observational confirmation for the Yarkovsky effect's theoretical size dependence of 1/D, where D is diameter. Second, we find that the observed ratio of negative to positive drift rates in our sample is 2.34, which, accounting for bias and sampling uncertainty, implies an actual ratio of $2.7^{+0.3}_{-0.7}$. This ratio has a vanishingly small probability of occurring due to chance or statistical noise. The observed ratio of retrograde to prograde rotators is two times lower than the ratio expected from numerical predictions from NEA population studies and traditional assumptions about the sense of rotation of NEAs originating from various main belt escape routes. We also examine the efficiency with which solar energy is converted into orbital energy and find a median efficiency in our sample of 12%. We interpret this efficiency in terms of NEA spin and thermal properties.

Advancing the scientific frontier in the search for life in the universe requires support of searches for both biosignatures and technosignatures. A modest budgetary increment can expand the search for life in the universe from primitive to complex life and from the solar neighborhood to the entire Galaxy.

As part of our ongoing search for technosignatures, we collected over three terabytes of data in May 2017 with the L-band receiver (1.15-1.73 GHz) of the 100 m diameter Green Bank Telescope. These observations focused primarily on planetary systems in the Kepler field, but also included scans of the recently discovered TRAPPIST-1 and LHS 1140 systems. We present the results of our search for narrowband signals in this data set with techniques that are generally similar to those described by Margot et al. (2018). Our improved data processing pipeline classified over $98\%$ of the $\sim$ 6 million detected signals as anthropogenic Radio Frequency Interference (RFI). Of the remaining candidates, 30 were detected outside of densely populated frequency regions attributable to RFI. These candidates were carefully examined and determined to be of terrestrial origin. We discuss the problems associated with the common practice of ignoring frequency space around candidate detections in radio technosignature detection pipelines. These problems include inaccurate estimates of figures of merit and unreliable upper limits on the prevalence of technosignatures. We present an algorithm that mitigates these problems and improves the efficiency of the search. Specifically, our new algorithm increases the number of candidate detections by a factor of more than four compared to Margot et al. (2018).

Planetary radar observations offer the potential for probing the properties of characteristics of solid bodies throughout the inner solar system and at least as far as the orbit of Saturn. In addition to the direct scientific value, precise orbital determinations can be obtained from planetary radar observations, which are in turn valuable for mission planning or spacecraft navigation and planetary defense. The next-generation Very Large Array would not have to be equipped with a transmitter to be an important asset in the world's planetary radar infrastructure. Bistatic radar, in which one antenna transmits (e.g., Arecibo or Goldstone) and another receives, are used commonly today, with the Green Bank Telescope (GBT) serving as a receiver. The improved sensitivity of the ngVLA relative to the GBT would improve the signal-to-noise ratios on many targets and increase the accessible volume specifically for asteroids. Goldstone-ngVLA bistatic observations would have the potential of rivaling the sensitivity of Arecibo, but with much wider sky access.

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 primary gravity science objective of NASA's Clipper mission to Europa is to confirm the presence or absence of a global subsurface ocean beneath Europa's Icy crust. Gravity field measurements obtained with a radio science investigation can reveal much about Europa's interior structure. Here, we conduct extensive simulations of the radio science measurements with the anticipated spacecraft trajectory and attitude (17F12v2) and assets on the spacecraft and the ground, including antenna orientations and beam patterns, transmitter characteristics, and receiver noise figures. In addition to two-way Doppler measurements, we also include radar altimeter crossover range measurements. We concentrate on +/-2 hour intervals centered on the closest approach of each of the 46 flybys. Our covariance analyses reveal the precision with which the tidal Love number k2, second-degree gravity coefficients C20 and C22, and higher-order gravity coefficients can be determined. The results depend on the Deep Space Network (DSN) assets that are deployed to track the spacecraft. We find that some DSN allocations are sufficient to conclusively confirm the presence or absence of a global ocean. Given adequate crossover range performance, it is also possible to evaluate whether the ice shell is hydrostatic.

Analysis of Kepler mission data suggests that the Milky Way includes billions of Earth-like planets in the habitable zone of their host star. Current technology enables the detection of technosignatures emitted from a large fraction of the Galaxy. We describe a search for technosignatures that is sensitive to Arecibo-class transmitters located within ~420 ly of Earth and transmitters that are 1000 times more effective than Arecibo within ~13 000 ly of Earth. Our observations focused on 14 planetary systems in the Kepler field and used the L-band receiver (1.15-1.73 GHz) of the 100 m diameter Green Bank Telescope. Each source was observed for a total integration time of 5 minutes. We obtained power spectra at a frequency resolution of 3 Hz and examined narrowband signals with Doppler drift rates between +/-9 Hz/s. We flagged any detection with a signal-to-noise ratio in excess of 10 as a candidate signal and identified approximately 850 000 candidates. Most (99%) of these candidate signals were automatically classified as human-generated radio-frequency interference (RFI). A large fraction (>99%) of the remaining candidate signals were also flagged as anthropogenic RFI because they have frequencies that overlap those used by global navigation satellite systems, satellite downlinks, or other interferers detected in heavily polluted regions of the spectrum. All 19 remaining candidate signals were scrutinized and none were attributable to an extraterrestrial source.

The population of large asteroids is thought to be primordial and they are the most direct witnesses of the early history of our Solar System. Those satellites allow study of the mass, and hence density and internal structure. We study here the properties of the triple asteroid (107) Camilla from lightcurves, stellar occultations, optical spectroscopy, and high-contrast and high-angular-resolution images and spectro-images. Using 80 positions over 15 years, we determine the orbit of its larger satellite to be circular, equatorial, and prograde, with RMS residuals of 7.8 mas. From 11 positions in three epochs only, in 2015 and 2016, we determine a preliminary orbit for the second satellite. We find the orbit to be somewhat eccentric and slightly inclined to the primary's equatorial plane, reminiscent of the inner satellites of other asteroid triple systems. Comparison of the near-infrared spectrum of the larger satellite reveals no significant difference with Camilla. Hence, these properties argue for a formation of the satellites by excavation from impact and re-accumulation of ejecta. We determine the spin and 3-D shape of Camilla. The model fits well each data set. We determine Camilla to be larger than reported from modeling of mid-infrared photometry, with a spherical-volume-equivalent diameter of 254 $\pm$ 36 km (3 $\sigma$ uncertainty), in agreement with recent results from shape modeling (Hanus2017+). Combining the mass of (1.12 $\pm$ 0.01) $\times$ 10$^{19}$ kg determined from the dynamics of the satellites and the volume from the 3-D shape model, we determine a density of 1,280 $\pm$ 130 SI. From this density, and considering Camilla's spectral similarities with (24) Themis and (65) Cybele (for which water ice coating on surface grains was reported), we infer a silicate-to-ice mass ratio of 1-6, with a 10-30% macroporosity.

We evaluated the prospects of quantifying the parameterized post-Newtonian parameter beta and solar quadrupole moment J2 with observations of near-Earth asteroids with large orbital precession rates (9 to 27 arcsec century$^{-1}$). We considered existing optical and radar astrometry, as well as radar astrometry that can realistically be obtained with the Arecibo planetary radar in the next five years. Our sensitivity calculations relied on a traditional covariance analysis and Monte Carlo simulations. We found that independent estimates of beta and J2 can be obtained with precisions of $6\times10^{-4}$ and $3\times10^{-8}$, respectively. Because we assumed rather conservative observational uncertainties, as is the usual practice when reporting radar astrometry, it is likely that the actual precision will be closer to $2\times10^{-4}$ and $10^{-8}$, respectively. A purely dynamical determination of solar oblateness with asteroid radar astronomy may therefore rival the helioseismology determination.

Near-Earth asteroid (NEA) 1566 Icarus ($a=1.08$ au, $e=0.83$, $i=22.8^{\circ}$) made a close approach to Earth in June 2015 at 22 lunar distances (LD). Its detection during the 1968 approach (16 LD) was the first in the history of asteroid radar astronomy. A subsequent approach in 1996 (40 LD) did not yield radar images. We describe analyses of our 2015 radar observations of Icarus obtained at the Arecibo Observatory and the DSS-14 antenna at Goldstone. These data show that the asteroid is a moderately flattened spheroid with an equivalent diameter of 1.44 km with 18% uncertainties, resolving long-standing questions about the asteroid size. We also solve for Icarus' spin axis orientation ($\lambda=270^{\circ}\pm10^{\circ}, \beta=-81^{\circ}\pm10^{\circ}$), which is not consistent with the estimates based on the 1968 lightcurve observations. Icarus has a strongly specular scattering behavior, among the highest ever measured in asteroid radar observations, and a radar albedo of $\sim$2\%, among the lowest ever measured in asteroid radar observations. The low cross-section suggests a high-porosity surface, presumably related to Icarus' cratering, spin, and thermal histories. Finally, we present the first use of our orbit determination software for the generation of observational ephemerides, and we demonstrate its ability to determine subtle perturbations on NEA orbits by measuring Icarus' orbit-averaged drift in semi-major axis ($(-4.62\pm0.48) \times 10^{-4}$ au/My, or $\sim$60 m per revolution). Our Yarkovsky rate measurement resolves a discrepancy between two published rates that did not include the 2015 radar astrometry.

Hi'iaka is the larger outer satellite of the dwarf planet Haumea. Using relative photometry from the Hubble Space Telescope and Magellan and a phase dispersion minimization analysis, we have identified the rotation period of Hi'iaka to be ~9.8 hrs (double-peaked). This is ~120 times faster than its orbital period, creating new questions about the formation of this system and possible tidal evolution. The rapid rotation suggests that Hi'iaka could have a significant obliquity and spin precession that could be visible in light curves within a few years. We then turn to an investigation of what we learn about the (presently unclear) formation of the Haumea system and family based on this unexpectedly rapid rotation rate. We explore the importance of the initial semi-major axis and rotation period in tidal evolution theory and find they strongly influence the time required to despin to synchronous rotation, relevant to understanding a wide variety of satellite and binary systems. We find that despinning tides do not necessarily lead to synchronous spin periods for Hi'iaka, even if it formed near the Roche limit. Therefore the short rotation period of Hi'iaka does not rule out significant tidal evolution. Hi'iaka's spin period is also consistent with formation near its current location and spin up due to Haumea-centric impactors.

We analyze radio tracking data obtained during 1311 orbits of the MESSENGER spacecraft in the period March 2011 to April 2014. A least-squares minimization of the residuals between observed and computed values of two-way range and Doppler allows us to solve for a model describing Mercury's gravity, tidal response, and spin state. We use a spherical harmonic representation of the gravity field to degree and order 40 and report error bars corresponding to 10 times the formal uncertainties of the fit. Our estimate of the product of Mercury's mass and the gravitational constant, $GM = (22031.87404 \pm 9 \times 10^{-4})$ km$^{3}$s$^{-2}$, is in excellent agreement with published results. Our solution for the geophysically important second-degree coefficients ($\bar{C}_{2,0} = -2.25100 \times 10^{-5} \pm 1.3 \times 10^{-9}$, $\bar{C}_{2,2} = 1.24973 \times 10^{-5} \pm 1.2 \times 10^{-9}$) confirms previous estimates to better than 0.4\% and, therefore, inferences about Mercury's moment of inertia and interior structure. Our estimate of the tidal Love number $k_2 = 0.464 \pm 0.023$ indicates that Mercury's mantle may be hotter and weaker than previously thought. Our spin state solution suggests that gravity-based estimates of Mercury's spin axis orientation are marginally consistent with previous measurements of the orientation of the crust.

We evaluated the planetary radar capabilities at Arecibo, the Goldstone 70-m DSS-14 and 34-m DSS-13 antennas, the 70-m DSS-43 antenna at Canberra, the Green Bank Telescope, and the Parkes Radio Telescope in terms of their relative sensitivities and the number of known near-Earth asteroids (NEAs) detectable per year in monostatic and bistatic configurations. In the 2015 calendar year, monostatic observations with Arecibo and DSS-14 were capable of detecting 253 and 131 NEAs respectively, with signal-to-noise ratios (SNRs) greater than 30/track. Combined, the two observatories were capable of detecting 276 NEAs. Of these, Arecibo detected 77 and Goldstone detected 32, or 30% and 24% the numbers that were possible. The two observatories detected an additional 18 and 7 NEAs respectively, with SNRs of less than 30/track. This indicates that a substantial number of potential targets are not being observed. The bistatic configuration with DSS-14 transmitting and the Green Bank Telescope receiving was capable of detecting about 195 NEAs, or ~50% more than with monostatic observations at DSS-14. Most of the detectable asteroids were targets of opportunity that were discovered less than 15 days before the end of their observing windows. About 50% of the detectable asteroids have absolute magnitudes > 25, which corresponds diameters < ~30 m.

A simple metric can be used to determine whether a planet or exoplanet can clear its orbital zone during a characteristic time scale, such as the lifetime of the host star on the main sequence. This criterion requires only estimates of star mass, planet mass, and orbital period, making it possible to immediately classify 99% of all known exoplanets. All 8 planets and all classifiable exoplanets satisfy the criterion. This metric may be useful in generalizing and simplifying the definition of a planet.

In the past decade, the number of known binary near-Earth asteroids has more than quadrupled and the number of known large main belt asteroids with satellites has doubled. Half a dozen triple asteroids have been discovered, and the previously unrecognized populations of asteroid pairs and small main belt binaries have been identified. The current observational evidence confirms that small (<20 km) binaries form by rotational fission and establishes that the YORP effect powers the spin-up process. A unifying paradigm based on rotational fission and post-fission dynamics can explain the formation of small binaries, triples, and pairs. Large (>20 km) binaries with small satellites are most likely created during large collisions.

We describe our implementation of a global-parameter optimizer and Square Root Information Filter (SRIF) into the asteroid-modelling software SHAPE. We compare the performance of our new optimizer with that of the existing sequential optimizer when operating on various forms of simulated data and actual asteroid radar data. In all cases, the new implementation performs substantially better than its predecessor: it converges faster, produces shape models that are more accurate, and solves for spin axis orientations more reliably. We discuss potential future changes to improve SHAPE's fitting speed and accuracy.

Potentially hazardous asteroid (185851) 2000 DP107 was the first binary near-Earth asteroid to be imaged. Radar observations in 2000 provided images at 75 m resolution that revealed the shape, orbit, and spin-up formation mechanism of the binary. The asteroid made a more favorable flyby of the Earth in 2008, yielding images at 30 m resolution. We used these data to obtain shape models for the two components and to improve the estimates of the mutual orbit, component masses, and spin periods. The primary has a sidereal spin period of 2.7745 +/- 0.0007 h and is roughly spheroidal with an equivalent diameter of 863 m +/- 5 %. It has a mass of 4.656 +/- 0.43 x 10^11 kg and a density of 1381 +/- 244 kg m^{-3}. It exhibits an equatorial ridge similar to the (66391) 1999 KW4 primary, however the equatorial ridge in this case is not as regular and has a ~300 m diameter concavity on one side. The secondary has a sidereal spin period of 1.77 +/- 0.02 days commensurate with the orbital period. The secondary is slightly elongated and has overall dimensions of 377 x 314 x 268 m (6 % uncertainties). Its mass is 0.178 +/- 0.021 x 10^{11} kg and its density is 1047 +/- 230 kg m^{-3}. The mutual orbit has a semi-major axis of 2.659 +/- 0.08 km, an eccentricity of 0.019 +/- 0.01, and a period of 1.7556 +/- 0.0015 days. The normalized total angular momentum of this system exceeds the amount required for the expected spin-up formation mechanism. An increase of angular momentum from non-gravitational forces after binary formation is a possible explanation. The two components have similar radar reflectivity, suggesting a similar composition consistent with formation by spin-up. The secondary appears to exhibit a larger circular polarization ratio than the primary, suggesting a rougher surface or subsurface at radar wavelength scales.

We develop a fourth-order numerical integrator to simulate the coupled spin and orbital motions of two rigid bodies having arbitrary mass distributions under the influence of their mutual gravitational potential. We simulate the dynamics of components in well-characterized binary and triple near-Earth asteroid systems and use surface of section plots to map the possible spin configurations of the satellites. For asynchronous satellites, the analysis reveals large regions of phase space where the spin state of the satellite is chaotic. For synchronous satellites, we show that libration amplitudes can reach detectable values even for moderately elongated shapes. The presence of chaotic regions in the phase space has important consequences for the evolution of binary asteroids. It may substantially increase spin synchronization timescales, explain the observed fraction of asynchronous binaries, delay BYORP-type evolution, and extend the lifetime of binaries. The variations in spin rate due to large librations also affect the analysis and interpretation of lightcurve and radar observations.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

We derive the locations of the fully synchronous end states of tidal evolution for binary asteroid systems having one spherical component and one oblate- or prolate-spheroid component. Departures from a spherical shape, at levels observed among binary asteroids, can result in the lack of a stable tidal end state for particular combinations of the system mass fraction and angular momentum, in which case the binary must collapse to contact. We illustrate our analytical results with near-Earth asteroids (8567) 1996 HW1, (66391) 1999 KW4, and 69230 Hermes.

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.

We observed near-Earth asteroid (162421) 2000 ET70 using the Arecibo and Goldstone radar systems over a period of 12 days during its close approach to the Earth in February 2012. We obtained continuous wave spectra and range-Doppler images with range resolutions as fine as 15 m. Inversion of the radar images yields a detailed shape model with an effective spatial resolution of 100 m. The asteroid has overall dimensions of 2.6 km X 2.2 km X 2.1 km (5% uncertainties) and a surface rich with kilometer-scale ridges and concavities. This size, combined with absolute magnitude measurements, implies an extremely low albedo (~2%). It is a principal axis rotator and spins in a retrograde manner with a sidereal spin period of 8.96 +/- 0.01 hours. In terms of gravitational slopes evaluated at scales of 100 m, the surface seems mostly relaxed with over 99% of the surface having slopes less than 30 degrees, but there are some outcrops at the north pole that may have steeper slopes. Our precise measurements of the range and velocity of the asteroid, combined with optical astrometry, enables reliable trajectory predictions for this potentially hazardous asteroid in the interval 460-2813.

We used a sample of Kepler candidate planets with orbital periods less than 200 days and radii between 1.5 and 30 Earth radii to determine the typical dynamical spacing of neighboring planets. To derive the intrinsic (i.e., free of observational bias) dynamical spacing of neighboring planets, we generated populations of planetary systems following various dynamical spacing distributions, subjected them to synthetic observations by the Kepler spacecraft, and compared the properties of observed planets in our simulations with actual Kepler detections. We found that, on average, neighboring planets are spaced 21.7 mutual Hill radii apart with a standard deviation of 9.5. This dynamical spacing distribution is consistent with that of adjacent planets in the Solar System. To test the packed planetary systems hypothesis, the idea that all planetary systems are dynamically packed or filled to capacity, we determined the fraction of systems that are dynamically packed by performing long-term (10^8 years) numerical simulations. In each simulation, we integrated a system with planets spaced according to our best-fit dynamical spacing distribution but containing an additional planet on an intermediate orbit. The fraction of simulations exhibiting signs of instability provides an approximate lower bound on the fraction of systems that are dynamically packed; we found that >31%, >35%, and >45% of 2-planet, 3-planet, and 4-planet systems are dynamically packed, respectively. Such sizeable fractions suggest that many planetary systems are indeed filled to capacity. This feature of planetary systems is another profound constraint that formation and evolution models must satisfy.

We investigated the underlying architecture of planetary systems by deriving the distribution of planet multiplicity (number of planets) and the distribution of orbital inclinations based on the sample of planet candidates discovered by the Kepler mission. The scope of our study included solar-like stars and planets with orbital periods less than 200 days and with radii between 1.5 and 30 Earth radii, and was based on Kepler planet candidates detected during Quarters 1 through 6. We created models of planetary systems with different distributions of planet multiplicity and inclinations, simulated observations of these systems by Kepler, and compared the properties of the transits of detectable objects to actual Kepler planet detections. Specifically, we compared with both the Kepler sample's transit numbers and normalized transit duration ratios in order to determine each model's goodness-of-fit. We did not include any constraints from radial velocity surveys. Based on our best-fit models, 75-80% of planetary systems have 1 or 2 planets with orbital periods less than 200 days. In addition, over 85% of planets have orbital inclinations less than 3 degrees (relative to a common reference plane). This high degree of coplanarity is comparable to that seen in our Solar System. These results have implications for planet formation and evolution theories. Low inclinations are consistent with planets forming in a protoplanetary disk, followed by evolution without significant and lasting perturbations from other bodies capable of increasing inclinations.

Sylvia is a triple asteroid system located in the main belt. We report new adaptive optics observations of this system that extend the baseline of existing astrometric observations to a decade. We present the first fully dynamical 3-body model for this system by fitting to all available astrometric measurements. This model simultaneously fits for individual masses, orbits, and primary oblateness. We find that Sylvia is composed of a dominant central mass surrounded by two satellites orbiting at 706.5 +/- 2.5 km and 1357 +/- 4.0 km, i.e., about 5 and nearly 10 primary radii. We derive individual masses of 1.484 -0.014/+0.016 x 10^19 kg for the primary (corresponding to a density of 1.29 +/- 0.39 g cm^-3), 7.33 -2.3/+4.7 x 10^14 kg for the inner satellite, and 9.32 -8.3/+20.7 x 10^14 kg for the outer satellite. The oblateness of the primary induces substantial precession and the J_2 value can be constrained to the range of 0.0985-0.1. The orbits of the satellites are relatively circular with eccentricities less than 0.04. The spin axis of the primary body and the orbital poles of both satellites are all aligned within about two degrees of each other, indicating a nearly coplanar configuration and suggestive of satellite formation in or near the equatorial plane of the primary. We also investigate the past orbital evolution of the system by simulating the effects of a recent passage through 3:1 mean-motion eccentricity-type resonances. In some scenarios this allow us to place constraints on interior structure and past eccentricities.

We have identified and quantified semi-major axis drifts in Near-Earth Asteroids (NEAs) by performing orbital fits to optical and radar astrometry of all numbered NEAs. We focus on a subset of 54 NEAs that exhibit some of the most reliable and strongest drift rates. Our selection criteria include a Yarkovsky sensitivity metric that quantifies the detectability of semi-major axis drift in any given data set, a signal-to-noise metric, and orbital coverage requirements. In 42 cases, the observed drifts (~10^-3 AU/Myr) agree well with numerical estimates of Yarkovsky drifts. This agreement suggests that the Yarkovsky effect is the dominant non-gravitational process affecting these orbits, and allows us to derive constraints on asteroid physical properties. In 12 cases, the drifts exceed nominal Yarkovsky predictions, which could be due to inaccuracies in our knowledge of physical properties, faulty astrometry, or modeling errors. If these high rates cannot be ruled out by further observations or improvements in modeling, they would be indicative of the presence of an additional non-gravitational force, such as that resulting from a loss of mass of order a kilogram per second. We define the Yarkovsky efficiency f_Y as the ratio of the change in orbital energy to incident solar radiation energy, and we find that typical Yarkovsky efficiencies are ~10^-5.

We investigate whether any multi-planet systems among Kepler candidates (2011 February release) can harbor additional terrestrial-mass planets or smaller bodies. We apply the "packed planetary systems" hypothesis that suggests all planetary systems are filled to capacity, and use a Hill stability criterion to identify eight 2-planet systems with significant gaps between the innermost and outermost planets. For each of these systems, we perform long-term numerical integrations of 10^7 years to investigate the stability of 4000-8000 test particles injected into the gaps. We map out stability regions in orbital parameter space, and therefore quantify the ranges of semi-major axes and eccentricities of stable particles. Strong mean-motion resonances can add additional regions of stability in otherwise unstable parameter space. We derive simple expressions for the extent of the stability regions, which is related to quantities such as the dynamical spacing Delta, the separation between two planets in units of their mutual Hill radii. Our results suggest that planets with separation Delta < 10 are unlikely to host extensive stability regions, and that about 95 out of a total of 115 two-planet systems in the Kepler sample may have sizeable stability regions. We predict that Kepler candidate systems including KOI 433, KOI 72/Kepler-10, KOI 555, KOI 1596, KOI 904, KOI 223, KOI 1590, and KOI 139 can harbor additional planets or low-mass bodies between the inner and outer detected planets. These predicted planets may be detected by future observations.

We investigate the Kozai mechanism in the context of near-Earth binaries and the Sun. The Kozai effect can lead to changes in eccentricity and inclination of the binary orbit, but it can be weakened or completely suppressed by other sources of pericenter precession, such as the oblateness of the primary body. Through numerical integrations including primary oblateness and 3 bodies (the two binary components and the Sun), we show that Kozai cycles cannot occur for the closely-separated near-Earth binaries in our sample. We demonstrate that this is due to pericenter precession around the oblate primary, even for very small oblateness values. Since the majority of observed near-Earth binaries are not well-separated, we predict that Kozai cycles do not play an important role in the orbital evolution of most near-Earth binaries. For a hypothetical wide binary modeled after 1998 ST27, the separation is large at 16 primary radii and so the orbital effects of primary oblateness are lessened. For this wide binary, we illustrate the possible excursions in eccentricity and inclination due to Kozai cycles as well as depict stable orientations for the binary's orbital plane. Unstable orientations lead to collisions between binary components, and we suggest that the Kozai effect acting in wide binaries may be a route to the formation of near-Earth contact binaries.

Many asteroids that make close encounters with terrestrial planets are in a binary configuration. Here we calculate the relevant encounter timescales and investigate the effects of encounters on a binary's mutual orbit. We use a combination of analytical and numerical approaches with a wide range of initial conditions. Our test cases include generic binaries with close, moderate, and wide separations, as well as seven well-characterized near-Earth binaries. We find that close approaches (<10 Earth radii) occur for almost all binaries on 1-10 million year timescales. At such distances, our results suggest substantial modifications to a binary's semi-major axis, eccentricity, and inclination, which we quantify. Encounters within 30 Earth radii typically occur on sub-million year timescales and significantly affect the wider binaries. Important processes in the lives of near-Earth binaries, such as tidal and radiative evolution, can be altered or stopped by planetary encounters.

In the near-Earth asteroid population, binary and triple systems have been discovered with mutual orbits that have significant eccentricities as well as large semi-major axes. All known systems with eccentric orbits and all widely-separated primary-satellite pairs have rapidly rotating satellites. Here we study processes that can elucidate the origin of these spin-orbital properties. Binary formation models based on rotational fissioning can reproduce asynchronous satellites on orbits with high eccentricities and a wide range of separations, but do not match observed properties. We explore whether any evolutionary mechanisms can link the spin and orbital parameters expected from post-fission dynamics to those observed today. We investigate four processes: tidal torques, radiative perturbations (BYORP), close planetary encounters, and Kozai oscillations. We find that a combination of post-fission dynamics and tidal evolution can explain nearly all the spin-orbit properties in a sample of nine well-characterized near-Earth binaries and triples. The other mechanisms may act but are not required to explain the observed data. Lastly, we describe evolutionary pathways between observed spin-orbital states including synchronous and circular, asynchronous and circular, and asynchronous and eccentric configurations.

The locations of the fully despun, double synchronous end states of tidal evolution are derived for spherical components. With the exception of nearly equal-mass binaries, binary asteroid systems are in the midst of lengthy tidal evolutions, far from their fully synchronous tidal end states. Calculations of material strength indicate that binaries in the main belt with 100-km-scale primary components are consistent with being made of monolithic or fractured rock as expected for binaries likely formed from sub-catastrophic impacts in the early solar system. To tidally evolve in their dynamical lifetime, near-Earth binaries with km-scale primaries or smaller must be much weaker mechanically than their main-belt counterparts even if formed in the main belt prior to injection into the near-Earth region. Small main-belt binaries with primary components less than 10 km in diameter, depending on their ages, could either be as strong as large main-belt binaries or as weak as near-Earth binaries because the inherent uncertainty in the age of a binary system can affect the calculation of material strength by orders of magnitude. Several other issues are considered, though these typically affect the calculation of material strength by no more than a factor of two. We also find indirect evidence within all three groups of binary asteroids that the inter-component separation may evolve via another mechanism(s) with the binary YORP effect being a likely candidate.

Three-body model fits to Arecibo and Goldstone radar data reveal the nature of two near-Earth asteroid triples. Triple-asteroid system 2001 SN263 is characterized by a primary of ~10^13 kg, an inner satellite ~1% as massive orbiting at ~3 primary radii in ~0.7 days, and an outer satellite ~2.5% as massive orbiting at ~13 primary radii in ~6.2 days. 1994 CC is a smaller system with a primary of mass ~2.6 \times 10^11 kg and two satellites ~2% and ~1% as massive orbiting at distances of ~5.5 and ~19.5 primary radii. Their orbital periods are ~1.2 and ~8.4 days. Examination of resonant arguments shows that the satellites are not currently in a mean-motion resonance. Precession of the apses and nodes are detected in both systems (2001 SN263 inner body: d{\varpi}/dt ~1.1 deg/day, 1994 CC inner body: d{\varpi}/dt ~ -0.2 deg/day), which is in agreement with analytical predictions of the secular evolution due to mutually interacting orbits and primary oblateness. Nonzero mutual inclinations between the orbital planes of the satellites provide the best fits to the data in both systems (2001 SN263: ~14 degrees, 1994 CC: ~16 degrees). Our best-fit orbits are consistent with nearly circular motion, except for 1994 CC's outer satellite which has an eccentric orbit of e ~ 0.19. We examine several processes that can generate the observed eccentricity and inclinations, including the Kozai and evection resonances, past mean-motion resonance crossings, and close encounters with terrestrial planets. In particular, we find that close planetary encounters can easily excite the eccentricities and mutual inclinations of the satellites' orbits to the currently observed values.

We observed the near-Earth asteroid 2008 EV5 with the Arecibo and Goldstone planetary radars and the Very Long Baseline Array during December 2008. EV5 rotates retrograde and its overall shape is a 400 /pm 50 m oblate spheroid. The most prominent surface feature is a ridge parallel to the asteroid's equator that is broken by a concavity 150 m in diameter. Otherwise the asteroid's surface is notably smooth on decameter scales. EV5's radar and optical albedos are consistent with either rocky or stony-iron composition. The equatorial ridge is similar to structure seen on the rubble-pile near-Earth asteroid (66391) 1999 KW4 and is consistent with YORP spin-up reconfiguring the asteroid in the past. We interpret the concavity as an impact crater. Shaking during the impact and later regolith redistribution may have erased smaller features, explaining the general lack of decameter-scale surface structure.

We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations less than five times the radius of the larger primary component, the tidal potential due to the presence of a smaller secondary component is expanded in terms of Legendre polynomials to arbitrary order rather than truncated at leading order as is typically done in studies of well-separated system like the Earth and Moon. The equations of tidal evolution including tidal torques, the changes in spin rates of the components, and the change in semimajor axis (orbital separation) are then derived for binary asteroid systems with circular and equatorial mutual orbits. Accounting for higher-order terms in the tidal potential serves to speed up the tidal evolution of the system leading to underestimates in the time rates of change of the spin rates, semimajor axis, and mean motion in the mutual orbit if such corrections are ignored. Special attention is given to the effect of close orbits on the calculation of material properties of the components, in terms of the rigidity and tidal dissipation function, based on the tidal evolution of the system. It is found that accurate determinations of the physical parameters of the system, e.g., densities, sizes, and current separation, are typically more important than accounting for higher-order terms in the potential when calculating material properties. In the scope of the long-term tidal evolution of the semimajor axis and the component spin rates, correcting for close orbits is a small effect, but for an instantaneous rate of change in spin rate, semimajor axis, or mean motion, the close-orbit correction can be on the order of tens of percent.

Observations with the adaptive optics system on the Very Large Telescope reveal that outer main belt asteroid (702) Alauda has a small satellite with primary to secondary diameter ratio of $\sim$56. The secondary revolves around the primary in 4.9143 $\pm$ 0.007 days at a distance of 1227 $\pm$ 24 km, yielding a total system mass of (6.057 $\pm$ 0.36) $\times$ 10$^{18}$ kg. Combined with an IRAS size measurement, our data yield a bulk density for this B-type asteroid of 1570 $\pm$ 500 kg~m$^{-3}$.

Detecting heat from minor planets in the outer solar system is challenging, yet it is the most efficient means for constraining the albedos and sizes of Kuiper Belt Objects (KBOs) and their progeny, the Centaur objects. These physical parameters are critical, e.g., for interpreting spectroscopic data, deriving densities from the masses of binary systems, and predicting occultation tracks. Here we summarize Spitzer Space Telescope observations of 47 KBOs and Centaurs at wavelengths near 24 and 70 microns. We interpret the measurements using a variation of the Standard Thermal Model (STM) to derive the physical properties (albedo and diameter) of the targets. We also summarize the results of other efforts to measure the albedos and sizes of KBOs and Centaurs. The three or four largest KBOs appear to constitute a distinct class in terms of their albedos. From our Spitzer results, we find that the geometric albedo of KBOs and Centaurs is correlated with perihelion distance (darker objects having smaller perihelia), and that the albedos of KBOs (but not Centaurs) are correlated with size (larger KBOs having higher albedos). We also find hints that albedo may be correlated with with visible color (for Centaurs). Interestingly, if the color correlation is real, redder Centaurs appear to have higher albedos. Finally, we briefly discuss the prospects for future thermal observations of these primitive outer solar system objects.

Binaries have played a crucial role many times in the history of modern astronomy and are doing so again in the rapidly evolving exploration of the Kuiper Belt. The large fraction of transneptunian objects that are binary or multiple, 48 such systems are now known, has been an unanticipated windfall. Separations and relative magnitudes measured in discovery images give important information on the statistical properties of the binary population that can be related to competing models of binary formation. Orbits, derived for 13 systems, provide a determination of the system mass. Masses can be used to derive densities and albedos when an independent size measurement is available. Angular momenta and relative sizes of the majority of binaries are consistent with formation by dynamical capture. The small satellites of the largest transneptunian objects, in contrast, are more likely formed from collisions. Correlations of the fraction of binaries with different dynamical populations or with other physical variables have the potential to constrain models of the origin and evolution of the transneptunian population as a whole. Other means of studying binaries have only begun to be exploited, including lightcurve, color, and spectral data. Because of the several channels for obtaining unique physical information, it is already clear that binaries will emerge as one of the most useful tools for unraveling the many complexities of transneptunian space.

We measured the system-integrated thermal emission of the binary Kuiper Belt Object 1999 TC36 at wavelengths near 24 and 70 microns using the Spitzer space telescope. We fit these data and the visual magnitude using both the Standard Thermal Model and thermophysical models. We find that the effective diameter of the binary is 405 km, with a range of 350 -- 470 km, and the effective visible geometric albedo for the system is 0.079 with a range of 0.055 -- 0.11. The binary orbit, magnitude contrast between the components, and system mass have been determined from HST data (Margot et al., 2004; 2005a; 2005b). Our effective diameter, combined with that system mass, indicate an average density for the objects of 0.5 g/cm3, with a range 0.3 -- 0.8 g/cm3. This density is low compared to that of materials expected to be abundant in solid bodies in the trans-Neptunian region, requiring 50 -- 75% of the interior of 1999 TC36 be taken up by void space. This conclusion is not greatly affected if 1999 TC36 is ``differentiated'' (in the sense of having either a rocky or just a non-porous core). If the primary is itself a binary, the average density of that (hypothetical) triple system would be in the range 0.4 -- 1.1 g/cm3, with a porosity in the range 15 -- 70%.

Arecibo S-band (lambda=13cm) radar observations of Comet C/2001 A2 (LINEAR) on 2001 July 7-9 showed a strong echo from large coma grains. This echo was significantly depolarized. This is the first firm detection of depolarization in a grain-coma radar echo and indicates that the largest grains are at least lambda / 2 or 2 cm in radius. The grains are moving at tens of m/s with respect to the nucleus. The non-detection of the nucleus places an upper limit of 3 km on its diameter. The broad, asymmetric echo power spectrum suggests a fan of grains that have a steep (differential number ~ a^-4) size distribution at cm-scales, though the observed fragmentation of this comet complicates that picture.