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.

This paper reports the first OH 18-cm line observation of the first detected interstellar object 1I/2017 U1 (`Oumuamua) using the Green Bank Telescope. We have observed the OH lines at 1665.402 MHz, 1667.359, and 1720.53 MHz frequencies with a spectral resolution of 357 Hz (approximately 0.06 km-s^{-1}). At the time of the observation, `Oumuamua was at topocentric distance and velocity of 1.07 au and 63.4 km-s^{-1}, respectively, or at heliocentric distance and velocity of 1.8 au and 39 km-s^{-1}, respectively. Based on a detailed data reduction and an analogy-based inversion, our final results confirm the asteroidal origin of `Oumuamua (as discussed in Meech et al., 2017) with an upper bound of OH production of Q[OH] < 0.17 x 10^{28} s^{-1}.

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.

As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by representing it as a cohesionless self-gravitating granular aggregate subject to rotational acceleration. To achieve this goal, a soft-sphere discrete element model (SSDEM) capable of simulating granular systems in quasi-static states is implemented and a quasi-static spin-up procedure is carried out. We devise three critical spin limits for the simulated aggregates to indicate their critical states triggered by reshaping and surface shedding, internal structural deformation, and shear failure, respectively. The failure condition and mode, and shear strength of an aggregate can all be inferred from the three critical spin limits. The effects of arrangement and size distribution of constituent particles, bulk density, spin-up path, and interparticle friction are numerically explored. The results show that the shear strength of a spinning self-gravitating aggregate depends strongly on both its internal configuration and material parameters, while its failure mode and mechanism are mainly affected by its internal configuration. Additionally, this study provides some constraints on the possible physical properties of the Didymos primary based on observational data and proposes a plausible formation mechanism for this binary system. With a bulk density consistent with observational uncertainty and close to the maximum density allowed for the asteroid, the Didymos primary in certain configurations can remain geo-statically stable without including cohesion.

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.

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.

An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to feedback essential informations to AIDA's ongoing phase-A study. We present a detailed dynamic model for the simulation of an ejecta cloud from a binary asteroid that synthesizes all relevant forces based on a previous analysis of the mechanical environment. We apply our method to gain insight into the expected response of Didymos to the AIDA impact, including the subsequent evolution of debris and dust. The crater scaling relations from laboratory experiments are employed to approximate the distributions of ejecta mass and launching speed. The size composition of fragments is modeled with a power law fitted from observations of real asteroid surface. A full-scale demonstration is simulated using parameters specified by the mission. We report the results of the simulation, which include the computed spread of the ejecta cloud and the recorded history of ejecta accretion and escape. The violent period of the ejecta evolution is found to be short, and is followed by a stage where the remaining ejecta is gradually cleared. Solar radiation pressure proves to be efficient in cleaning dust-size ejecta, and the simulation results after two weeks shows that large debris on polar orbits (perpendicular to the binary orbital plane) has a survival advantage over smaller ejecta and ejecta that keep to low latitudes.

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.

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 propose a new equation of state for the Dark Energy component of the Universe. It is modeled on the equation of state $p=w(\rho-\rho_{*})$ which can describe a liquid, for example water. We show that its energy density naturally decomposes into a component that behaves as a cosmological constant and one whose energy density scales as $a^{-3(1+w)}$, and fit the parameters specifying the equation of state to the new SNIa data, as well as WMAP and 2dF data. We find that reasonable values of the parameters can be found that give our model the same $\chi^2$ as that of $\Lambda$CDM. A remarkable feature of the model is that we can do all this with $w>0$.