Nebular gas drag and planetary accretion with eccentric high-mass planets

Aims.We investigate the dynamics of pebbles immersed in a gas disk interacting with a planet on an eccentric orbit. The model has a prescribed gap in the disk around the location of the planetary orbit, as is expected for a giant planet with a mass in the range of 0.1-1 Jupiter masses. The pebbles with sizes in the range of 1 cm to 3 m are placed in a ring outside of the giant planet orbit at distances between 10 and 30 planetary Hill radii. The process of the accumulation of pebbles closer to the gap edge, its possible implication for the planetary accretion, and the importance of the mass and the eccentricity of the planet in this process are the motivations behind the present contribution. Methods. We used the Bulirsch-Stoer numerical algorithm, which is computationally consistent for close approaches, to integrate the Newtonian equations of the planar (2D), elliptical restricted three-body problem. The angular velocity of the gas disk was determined by the appropriate balance between the gravity, centrifugal, and pressure forces, such that it is sub-Keplerian in regions with a negative radial pressure gradient and super-Keplerian where the radial pressure gradient is positive. Results. The results show that there are no trappings in the 1:1 resonance around the L 4 and L5 Lagrangian points for very low planetary eccentricities (e2 < 0.07). The trappings in exterior resonances, in the majority of cases, are because the angular velocity of the disk is super-Keplerian in the gap disk outside of the planetary orbit and because the inward drift is stopped. Furthermore, the semi-major axis location of such trappings depends on the gas pressure profile of the gap (depth) and is a = 1.2 for a planet of 1 MJ. A planet on an eccentric orbit interacts with the pebble layer formed by these resonances. Collisions occur and become important for planetary eccentricity near the present value of Jupiter (e 2 = 0.05). The maximum rate of the collisions onto a planet of 0.1 MJ occurs when the pebble size is 37.5 cm ≤ s < 75 cm; for a planet with the mass of Jupiter, it is15 cm ≤ s < 30 cm. The accretion stops when the pebble size is less than 2 cm and the gas drag dominates the motion. © 2013 ESO.

A frequency approach to identifying asteroid families II. Families interacting with nonlinear secular resonances and low-order mean-motion resonances

Aims. In an earlier paper we introduced a new method for determining asteroid families where families were identified in the proper frequency domain (n, g, g + s) ( where n is the mean-motion, and g and s are the secular frequencies of the longitude of pericenter and nodes, respectively), rather than in the proper element domain (a, e, sin(i)) (semi-major axis, eccentricity, and inclination). Here we improve our techniques for reliably identifying members of families that interact with nonlinear secular resonances of argument other than g or g + s and for asteroids near or in mean-motion resonant configurations. Methods. We introduce several new distance metrics in the frequency space optimal for determining the diffusion in secular resonances of argument 2g - s, 3g - s, g - s, s, and 2s. We also regularize the dependence of the g frequency as a function of the n frequency (Vesta family) or of the eccentricity e (Hansa family). Results. Our new approaches allow us to recognize as family members objects that were lost with previous methods, while keeping the advantages of the Carruba & Michtchenko (2007, A& A, 475, 1145) approach. More important, an analysis in the frequency domain permits a deeper understanding of the dynamical evolution of asteroid families not always obtainable with an analysis in the proper element domain.

Physical and dynamical characterization of (5201) Ferraz-Mello, a possible extinct Jupiter family comet

Context. The subject of asteroids in cometary orbits (ACOs) has been of growing interest lately. These objects have the orbital characteristics typical of comets, but are asteroidal in appearance, i.e., show no signs of a coma at any part of their orbits. At least a fraction of these objects are thought to be comets that have either exhausted all their volatile content or developed a refractory crust that prevents sublimation. In particular, the asteroid ( 5201) Ferraz-Mello has, since its discovery, been suspected to be an extinct Jupiter family comet due to the peculiar nature of its orbit. Aims. The aim of this work is to put constraints on the possible origin of ( 5201) Ferraz-Mello by means of spectroscopic characterization and a study of the dynamics of this asteroid. Methods. We used the SOAR Optical Imager (SOI) to obtain observations of ( 5201) Ferraz-Mello using four SDSS filters. These observations were compared to asteroids listed in the Sloan Moving objects catalog and also to photometry of cometary nuclei, Centaurs, and TNOs. The orbital evolution of ( 5201) Ferraz-Mello and of a sample of asteroids and comets that are close to that object in the a - e plane were simulated using a pure N-body code for 4 000 years forward and 4 000 years backward in time. Results. The reflectance spectrum obtained from its colors in the SDSS system is unusual, with a steep spectral gradient that is comparable to TNOs and Centaurs, but with an increase in the reflectance in the g band that is not common in those populations. A similar behavior is seen in cometary nuclei that were observed in the presence of a faint dust coma. The dynamical results confirm the very chaotic evolution found previously and its dynamical similarity to the chaotic evolution of some comets. The asteroid is situated in the very stochastic layer at the border of the 2/1 resonance, and it has a very short Lyapunov time ( 30 - 40) years. Together, the spectral characteristcs and the dynamical evolution suggest that ( 5201) Ferraz-Mello is a dormant or extinct comet.

Dynamic picture of the inner asteroid belt: implications for the density, size and taxonomic distributions of real objects

In this paper, we construct a dynamic portrait of the inner asteroidal belt. We use information about the distribution of test particles, which were initially placed on a perfectly rectangular grid of initial conditions, after 4.2 Myr of gravitational interactions with the Sun and five planets, from Mars to Neptune. Using the spectral analysis method introduced by Michtchenko et al., the asteroidal behaviour is illustrated in detail on the dynamical, averaged and frequency maps. On the averaged and frequency maps, we superpose information on the proper elements and proper frequencies of real objects, extracted from the data base, AstDyS, constructed by Milani and Knezevic. A comparison of the maps with the distribution of real objects allows us to detect possible dynamical mechanisms acting in the domain under study; these mechanisms are related to mean-motion and secular resonances. We note that the two- and three-body mean-motion resonances and the secular resonances (strong linear and weaker non-linear) have an important role in the diffusive transportation of the objects. Their long-lasting action, overlaid with the Yarkovsky effect, may explain many observed features of the density, size and taxonomic distributions of the asteroids.

Secular dynamics and family identification among highly inclined asteroids in the Euphrosyne region

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The stable archipelago in the region of the Pallas and Hansa dynamical families

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

On the stability of the satellites of asteroid 87 Sylvia

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Exploring S-type orbits in the Pluto-Charon binary system

This work generates, through a sample of numerical simulations of the restricted three-body problem, diagrams of semimajor axis and eccentricity which defines stable and unstable zones for particles in S-type orbits around Pluto and Charon. Since we consider initial conditions with 0 <= e <= 0.99, we found several new stable regions. We also identified the nature of each one of these newly found stable regions. They are all associated to families of periodic orbits derived from the planar circular restricted three-body problem. We have shown that a possible eccentricity of the Pluto-Charon system slightly reduces, but does not destroy, any of the stable regions.

Dynamical erosion of asteroid groups in the region of the Phocaea family

In a previous paper, the current state of knowledge of the region containing the Phocaea dynamical family was revised. Here, the dynamical evolution and possible origin of the Phocaea dynamical family and asteroid groups in the region are investigated. First, I study the case of asteroids at high eccentricity (e > 0.31). I find that these objects are unstable because of encounters with Mars on time-scales of up to 270 Myr. The minimum time needed by members of the Phocaea classical family to reach the orbital locations of these objects, 370 Myr, can be used to set a lower limit on the age of the Phocaea family.Next, attention is focused on the chaotic layer previously identified near the nu(6) secular resonance border. Using analytical and numerical tools, I find that the presence of the nu(6) secular resonance forces asteroids with vertical bar g-g(6)vertical bar < 2.55 arcsec yr(-1) to reach eccentricities high enough to allow them to experience deep, close encounters with Mars. Results of the analytical model of Yoshikawa and of my numerical simulations fully explain the low-inclination chaotic region found by Carruba.Finally, I investigate the long-term stability of the minor families and clumps identified in the previous paper, with particular emphasis on a clump only identifiable in the domain of proper frequencies (n, g, g - s) around (6246) Komurotoru. I find that while the clumps identified in the space of proper elements quickly disperse when the Yarkovsky effect is considered, the family around (19536) is still observable for time-scales of more than 50 Myr. The (6246) clump, characterized by its interaction with the nu(5) + nu(16) and 2 nu(6) - nu(16) secular resonances, is robust on time-scales of 50 Myr. I confirm that this group may be the first clump ever detected in the frequency domain that can be associated with a real collisional event.

Dynamical erosion of asteroid groups in the region of the Pallas family

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

On the Emmenthal distribution of highly inclined asteroids

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Chaotic diffusion caused by close encounters with several massive asteroids The (4) Vesta case

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

An analysis of the region of the Phocaea dynamical family

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

On the V-type asteroids outside the Vesta family - II. Is (21238) 1995 WV7 a fragment of the long-lost basaltic crust of (15) Eunomia?

Context. The V-type asteroids are associated with basaltic composition. Apart from ( 1459) Magnya, an asteroid that is clearly dynamically and mineralogically unconnected to the Vesta family, all currently known V-type asteroids are either members of the Vesta family, or are hypothesized to be former members of the dynamical family that migrated to their current orbital positions. The recent identification of ( 21238) 1995 WV7 as a V-type asteroid introduces the possibility that a second basaltic asteroid not connected with the Vesta family exists. This asteroid is on the opposite side of the 3: 1 mean motion resonance with respect to Vesta, and it would be very unlikely that a member of the Vesta family of its size (D > 5km) migrating via either the Yarkovsky effect or repeated close encounters with Vesta survived the passage through such a resonance.Aims. In this work we investigate the possibility that ( 21238) 1995 WV7 originated as a fragment of the parent body of the Eunomia family and then migrated via the interplay of the Yarkovsky effect and some powerful nonlinear secular resonances, such as the (s - s(6)) - ( g(5) - g(6)). If (15) Eunomia is, as claimed, a differentiated object whose originally pyroxene-enriched crust layer was lost in a collision that either created the Eunomia family or preceded its formation, can (21238) be a fragment of its long-lost basaltic crust that migrated to the current position?Methods. We mapped the phase space around (21238) and determined which of the nonlinear secular resonances that we identified are stronger and more capable of having caused the current difference in proper i between (21238) and members of the Eunomia family. We simulated the Yarkovsky effect by using the SWIFT-RMVSY integrator.Results. Our results suggest that it is possible to migrate from the Eunomia dynamical family to the current orbital location of ( 21238) via the interplay of the Yarkovsky effect and the (s - s6) - (g5 - g6) nonlinear secular resonance, on time-scales of at least 2.6 Gyr.Conclusions. (15) Eunomia might be the third currently known parent body for V-type asteroids.