An Oort cloud origin for the high-inclination, high-perihelion Centaurs

We analyse the origin of three Centaurs with perihelia in the range 15 AU to 30 AU, inclinations above 70 deg and semi-major axes shorter than 100 AU. Based on long-term numerical simulations we conclude that these objects most likely originate from the Oort cloud rather than the Kuiper Belt or Scattered Disc. We estimate that there are currently between 1 and 200 of these high-inclination, high-perihelion Centaurs with absolute magnitude H<8.Comment: Accepted for publication in MNRA

Dynamical evolution of escaped plutinos, another source of Centaurs

It was shown in previous works the existence of weakly chaotic orbits in the plutino population that diffuse very slowly. These orbits correspond to long-term plutino escapers and then represent the plutinos that are escaping from the resonance at present. In this paper we perform numerical simulations in order to explore the dynamical evolution of plutinos recently escaped from the resonance. The numerical simulations were divided in two parts. In the first one we evolved 20,000 test particles in the resonance in order to detect and select the long-term escapers. In the second one, we numerically integrate the selected escaped plutinos in order to study their dynamical post escaped behavior. Our main results include the characterization of the routes of escape of plutinos and their evolution in the Centaur zone. We obtained a present rate of escape of plutinos between 1 and 10 every 10 years. The escaped plutinos have a mean lifetime in the Centaur zone of 108 Myr and their contribution to the Centaur population would be a fraction of less than 6 % of the total Centaur population. In this way, escaped plutinos would be a secondary source of Centaurs.Comment: Accepted for publication in A&

The formation of Uranus and Neptune among Jupiter and Saturn

The outer giant planets, Uranus and Neptune, pose a challenge to theories of planet formation. They exist in a region of the Solar System where long dynamical timescales and a low primordial density of material would have conspired to make the formation of such large bodies ($\sim$ 15 and 17 times as massive as the Earth, respectively) very difficult. Previously, we proposed a model which addresses this problem: Instead of forming in the trans-Saturnian region, Uranus and Neptune underwent most of their growth among proto-Jupiter and -Saturn, were scattered outward when Jupiter acquired its massive gas envelope, and subsequently evolved toward their present orbits. We present the results of additional numerical simulations, which further demonstrate that the model readily produces analogues to our Solar System for a wide range of initial conditions. We also find that this mechanism may partly account for the high orbital inclinations observed in the Kuiper belt.Comment: Submitted to AJ; 38 pages, 16 figure

Symplectic integrators with adaptive time steps

In recent decades, there have been many attempts to construct symplectic integrators with variable time steps, with rather disappointing results. In this paper we identify the causes for this lack of performance, and find that they fall into two categories. In the first, the time step is considered a function of time alone, \Delta=\Delta(t). In this case, backwards error analysis shows that while the algorithms remain symplectic, parametric instabilities arise because of resonance between oscillations of \Delta(t) and the orbital motion. In the second category the time step is a function of phase space variables \Delta=\Delta(q,p). In this case, the system of equations to be solved is analyzed by introducing a new time variable \tau with dt=\Delta(q,p) d\tau. The transformed equations are no longer in Hamiltonian form, and thus are not guaranteed to be stable even when integrated using a method which is symplectic for constant \Delta. We analyze two methods for integrating the transformed equations which do, however, preserve the structure of the original equations. The first is an extended phase space method, which has been successfully used in previous studies of adaptive time step symplectic integrators. The second, novel, method is based on a non-canonical mixed-variable generating function. Numerical trials for both of these methods show good results, without parametric instabilities or spurious growth or damping. It is then shown how to adapt the time step to an error estimate found by backward error analysis, in order to optimize the time-stepping scheme. Numerical results are obtained using this formulation and compared with other time-stepping schemes for the extended phase space symplectic method.Comment: 23 pages, 9 figures, submitted to Plasma Phys. Control. Fusio

Implications of the Small Spin Changes Measured for Large Jupiter-Family Comet Nuclei

Rotational spin-up due to outgassing of comet nuclei has been identified as a possible mechanism for considerable mass-loss and splitting. We report a search for spin changes for three large Jupiter-family comets (JFCs): 14P/Wolf, 143P/Kowal-Mrkos, and 162P/Siding Spring. None of the three comets has detectable period changes, and we set conservative upper limits of 4.2 (14P), 6.6 (143P) and 25 (162P) minutes per orbit. Comparing these results with all eight other JFCs with measured rotational changes, we deduce that none of the observed large JFCs experiences significant spin changes. This suggests that large comet nuclei are less likely to undergo rotationally-driven splitting, and therefore more likely to survive more perihelion passages than smaller nuclei. We find supporting evidence for this hypothesis in the cumulative size distributions of JFCs and dormant comets, as well as in recent numerical studies of cometary orbital dynamics. We added 143P to the sample of 13 other JFCs with known albedos and phase-function slopes. This sample shows a possible correlation of increasing phase-function slopes for larger geometric albedos. Partly based on findings from recent space missions to JFCs, we hypothesise that this correlation corresponds to an evolutionary trend for JFCs. We propose that newly activated JFCs have larger albedos and steeper phase functions, which gradually decrease due to sublimation-driven erosion. If confirmed, this could be used to analyse surface erosion from ground and to distinguish between dormant comets and asteroids

Testing the comet nature of main belt comets. The spectra of 133P/Elst-Pizarro and 176P/LINEAR

We present the visible spectrum of MBCs 133P/Elst-Pizarro and 176P/LINEAR, as well as three Themis family asteroids: (62) Erato, (379), Huenna and (383) Janina, obtained in 2007 using three telescopes at "El Roque de los Muchachos"' Observatory, in La Palma, Spain, and the 8m Kueyen (UT2) VLT telescope at Cerro Paranal, Chile. The spectra of 133P and 176P resemble best those of B-type asteroid and are very similar to those of Themis family members and are significantly different from the spectrum of comet 162P/Siding-Spring and most of the observed cometary nuclei. CN gas emission is not detected in the spectrum of 133P. We determine an upper limit for the CN production rate Q(CN) = $= 2.8 \times 10^{21}$ mol/s, three orders of magnitude lower than the Q(CN) of Jupiter family comets observed at similar heliocentric distances. The spectra of 133P/Elst-Pizarro and 176P/LINEAR confirm that they are likely members of the Themis family of asteroids, fragments that probably retained volatiles, and unlikely have a cometary origin in the trans-neptunian belt or the Oort cloud.Comment: Paper sumbmited to A&A. 7 pages and 6 figure

Water in the Near IR spectrum of Comet 8P/Tuttle

High resolution spectra of Comet 8P/Tuttle were obtained in the frequency range 3440.6-3462.6 cm-1 on 3 January 2008 UT using CGS4 with echelle grating on UKIRT. In addition to recording strong solar pumped fluorescent (SPF) lines of H2O, the long integration time (152 miutes on target) enabled eight weaker H2O features to be assigned, most of which had not previously been identified in cometary spectra. These transitions, which are from higher energy upper states, are similar in character to the so-called 'SH' lines recorded in the post Deep Impact spectrum of comet Tempel 1 (Barber et al., 2007). We have identified certain characteristics that these lines have in common, and which in addition to helping to define this new class of cometary line, give some clues to the physical processes involved in their production. Finally, we derive an H2O rotational temperature of 62+/- K and a water production rate of (1.4+/-0.3)E28 molecules/s.Comment: Paper has been accepted for publication by MNRAS (11/06/09

The orbit, mass, size, albedo, and density of (65489) Ceto/Phorcys: A tidally-evolved binary Centaur

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