We investigate the influence of the turbulence forcing on the mass distributions of gravitationally unstable cores by postprocessing data from simulations of non-selfgravitating isothermal supersonic turbulence with varying resolution. In one set of simulations solenoidal forcing is applied, while the second set uses purely compressive forcing to excite turbulent motions. From the resulting density field, we compute the mass distribution of gravitationally unstable cores by means of a clump-finding algorithm. Using the time-averaged probability density functions of the mass density, semi-analytic mass distributions are calculated from analytical theories. We apply stability criteria that are based on the Bonnor-Ebert mass resulting from the thermal pressure and from the sum of thermal and turbulent pressure. Although there are uncertainties in the application of the clump-finding algorithm, we find systematic differences in the mass distributions obtained from solenoidal and compressive forcing. Compressive forcing produces a shallower slope in the high-mass power-law regime compared to solenoidal forcing. The mass distributions also depend on the Jeans length resulting from the choice of the mass in the computational box, which is freely scalable for non-selfgravitating isothermal turbulence. Provided that all cores are numerically resolved and most cores are small compared to the length scale of the forcing, the normalised core mass distributions are found to be close to the semi-analytic models. Especially for the high-mass tails, the Hennebelle-Chabrier theory implies that the additional support due to turbulent pressure is important.

We performed numerical simulations of supersonic isothermal turbulence driven by mostly compressive large-scale forcing, using both a static grid and adaptive mesh refinement with an effective resolution N=768^3. After a transient phase dominated by shocks, turbulence evolves into a steady state with an RMS Mach number about 2.5, in which cloud-like structures of over-dense gas are surrounded by highly rarefied gas. The index of the turbulence energy spectrum function beta = 2.0 in the shock-dominated phase. As the flow approaches statistical equilibrium, the spectrum flattens, with beta = 1.9. For the scaling exponent of the root mean square velocity fluctuation, we obtain gamma = 0.43 from the velocity structure functions of second order. These results are well within the range of observed scaling properties for the velocity dispersion in molecular clouds. Calculating structure functions of order p=1,...,5, we find for all scaling exponents significant deviations from the Kolmogorov-Burgers model proposed by Boldyrev. Our results are very well described by a general log-Poisson model with a higher degree of intermittency, which implies an influence of the forcing on the scaling properties. Contrary to previous numerical results for isothermal turbulence, we obtain a skewed probability density function of the mass density fluctuations that is not consistent with log-normal statistics and entails a substantially higher fraction of mass in the density peaks than implied by the Padoan-Nordlund relation between the variance of the density fluctuations and the Mach number. In conclusion, it seems necessary to account for the production mechanism of turbulence in the ISM.

We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971 +/- 0.006 days and a semimajor axis of 1628 +/- 29 km, implying a system mass of (9.49 +/- 0.52) x 10^17 kg. The eccentricity of the orbit is 0.526 +/- 0.015. Combined with a radiometric size determined from Spitzer Space Telescope data and the assumption that Typhon and Echidna both have the same albedo, we estimate that their radii are 76 +14/-16 and 42 +8/-9 km, respectively. These numbers give an average bulk density of only 0.44 +0.44/-0.17 g cm^-3, consistent with very low bulk densities recently reported for two other small transneptunian binaries.

We have searched 101 Classical transneptunian objects for companions with the Hubble Space Telescope. Of these, at least 21 are binary. The heliocentric inclinations of the objects we observed range from 0.6-34 degrees. We find a very strong anticorrelation of binaries with inclination. Of the 58 targets that have inclinations of less than 5.5 degrees, 17 are binary, a binary fraction of 29 +7/-6%. All 17 are similar-brightness systems. On the contrary, only 4 of the 42 objects with inclinations greater than 5.5 degrees have satellites and only 1 of these is a similar-brightness binary. This striking dichotomy appears to agree with other indications that the low eccentricity, non-resonant Classical transneptunian objects include two overlapping populations with significantly different physical properties and dynamical histories.

Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX128) reveal it to be a close binary system. The mutual orbit has a period of 9.554 +/- 0.011 days and a semimajor axis of 1840 +/- 48 km. These values enable computation of a system mass of (5.41 +/- 0.42) 10^18 kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 microns are combined with visible photometry to constrain the system's effective radius (109 +10/-11 km) and geometric albedo (0.084 +0.021/-0.014). We estimate the average bulk density to be 1.37 +0.66/-0.32 g cm^-3, consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM165, and (47171) 1999 TC36, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity <= 0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the solar system.

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.

The Deep Ecliptic Survey is a project whose goal is to survey a large area of the near-ecliptic region to a faint limiting magnitude (R ~ 24) in search of objects in the outer solar system. We are collecting a large homogeneous data sample from the Kitt Peak Mayall 4-m and Cerro Tololo Blanco 4-m telescopes with the Mosaic prime-focus CCD cameras. Our goal is to collect a sample of 500 objects with good orbits to further our understanding of the dynamical structure of the outer solar system. This survey has been in progress since 1998 and is responsible for 272 designated discoveries as of March 2003. We summarize our techniques, highlight recent results, and describe publically available resources.

As part of our ongoing Deep Ecliptic Survey (DES) of the Kuiper belt, we report on the occupation of the 1:1 (Trojan), 4:3, 3:2, 7:4, 2:1, and 5:2 Neptunian mean-motion resonances (MMRs). The occupation of the 1:1 and 5:2 MMRs is not easily understood within the standard model of resonance sweeping by a migratory Neptune over an initially dynamically cold belt. Our dynamically hot, 5:2 resonant objects can librate with modest amplitudes of 90 deg within the resonance for at least 1 Gyr. Their trajectories cannot be explained by close encounters with Neptune alone, given the latter's current orbit. The hot orbits of such 5:2 resonant KBOs, unlike hot orbits of previously known resonant KBOs, may imply that these objects were pre-heated to large inclination and large eccentricity prior to resonance capture by a migratory Neptune. Our first discovered Neptunian Trojan, 2001QR322, may not owe its existence to Neptune's migration at all. The trajectory of 2001QR322 is remarkably stable; the object can undergo tadpole-type libration about Neptune's leading Lagrange (L4) point for at least 1 Gyr with a libration amplitude of 24 deg. Trojan capture probably occurred while Neptune accreted the bulk of its mass. For an assumed albedo of 12--4%, our Trojan is 130--230 km in diameter. Model-dependent estimates place the total number of Neptune Trojans resembling 2001QR322 at 20--60. Their existence might rule out violent orbital histories for Neptune.