Massive identification of asteroids in three-body resonances

An essential role in the asteroidal dynamics is played by the mean motion resonances. Two-body planet–asteroid resonances are widely known, due to the Kirkwood gaps. Besides, so-called three-body mean motion resonances exist, in which an asteroid and two planets participate. Identification of asteroids in three-body (namely, Jupiter–Saturn–asteroid) resonances was initially accomplished by Nesvorný and Morbidelli (Nesvorný D., Morbidelli, A. [1998]. Astron. J. 116, 3029–3037), who, by means of visual analysis of the time behaviour of resonant arguments, found 255 asteroids to reside in such resonances. We develop specialized algorithms and software for massive automatic identification of asteroids in the three-body, as well as two-body, resonances of arbitrary order, by means of automatic analysis of the time behaviour of resonant arguments. In the computation of orbits, all essential perturbations are taken into account. We integrate the asteroidal orbits on the time interval of 100,000 yr and identify main-belt asteroids in the three-body Jupiter–Saturn–asteroid resonances up to the 6th order inclusive, and in the two-body Jupiter–asteroid resonances up to the 9th order inclusive, in the set of ～250,000 objects from the “Asteroids – Dynamic Site” (AstDyS) database. The percentages of resonant objects, including extrapolations for higher-order resonances, are determined. In particular, the observed fraction of pure-resonant asteroids (those exhibiting resonant libration on the whole interval of integration) in the three-body resonances up to the 6th order inclusive is ≈0.9% of the whole set; and, using a higher-order extrapolation, the actual total fraction of pure-resonant asteroids in the three-body resonances of all orders is estimated as ≈1.1% of the whole set.     

The high-eccentricity libration theory revisted

Using the local asymmetric expansion of the disturbing function for the planar elliptic restricted three-body problem up to degree 1 ine ₁ we develop a small amplitude libration theory. We review the laws that characterize the mean-motion resonances of asteroids with Jupiter and we obtain other new laws. Special attention is paid to the second forced mode whose equations are reformulated and new consequences of this component are discussed. An analytical expression for the trajectories in the phase space is obtained. The predictions are compared with numerical experiments which confirm the results.     

Flattening,pole, and albedo features of 4 vesta from photometric data

A number of large asteroids show irregular lightcurves of relatively small amplitude and/or ambiguous rotational periods. These observations and the fact that their strong gravitational binding probably results in quasi-equilibrium shapes lead to model these bodies as axisymmetric, biaxial ellipsoids covered by albedo markings. We developed a general numerical algorithm for obtaining simulated lightcurves of “spotted” asteroids and varied the most critical geometrical and physical parameters (albedo contrast, size, and position of the spots; polar coordinates, and shape of the asteroid). We then analyzed the case of 4 Vesta by assuming an axisymmetric ellipsoidal shape with a large brighter region on one hemisphere, in agreement with the results of photometric and polarimetric observations. Fitting the numerical simulations to the available data, we obtained the flattening of the ellipsoid (0.79 ± 0.03), the albedo contrast and geometry of the brighter region, and the orientation of the polar axis. If the derived flattenning corresponds to the equilibrium shape of a nearly homogeneous body, a density of 2.4 ± 0.3 g cm⁻³ can be inferred. These results show satisfactory agreement with values by different techniques. We plan to apply the same method both to other large asteroids and to smaller, irregularly shaped ones; in the latter case, this will allow us to test the uncertainties in current pole determination methods.     

Cratering rate on the jovian system: the contribution from Hilda asteroids

We study the dynamical evolution of the Hilda group of asteroids trough numerical methods, performing also a collisional pseudo-evolution of the present population, in order to calculate the rate of evaporation and its contribution to the cratering history of the Galilean satellites. If the present population of small asteroids in the Hilda's region follows the same size distribution observed at larger radii, we find that this family is the main contributor to the production of small craters (i.e., crater with diameters d∼4 km) on the Galilean system, overcoming the production by Jupiter Family Comets and by Trojan asteroids. The results of this investigation encourage further observational campaigns, in order to determine the size distribution function of small Hilda asteroids.     

On Stable Chaos in the Asteroid Belt

The twenty most chaotic objects found among first hundred of numbered asteroids are studied. Lyapunov time calculated with and without inner planets indicates that for eleven of those asteroids the strongest chaotic effect results from the resonances with Mars. The filtered semimajor axis displays an abrupt variation only when a close approach to Mars takes place. The study of the behaviour of the critical argument for candidate resonances can reveal which is responsible for the semimajor axis variation. We have determined these resonances for the asteroids in question. For the asteroids chaotic even without the inner planets we have determined the most important resonances with Jupiter, or three-body resonances. This revised version was published online in July 2006 with corrections to the Cover Date.     

Probing the Nekhoroshev Stability of Asteroids

We apply the spectral formulation of the Nekhoroshev theorem to investigate the long-term stability of real main belt asteroids. We find numerical indication that some asteroids are in the so-called Nekhoroshev stability regime, that is they are on chaotic orbits but their motion is stable over very long times. We have analyzed the motion of bodies in different regions of the belt, to assess the sensitivity of our method. We found that it allows us to clearly discriminate between different dynamical regimes, such as the one described by the Nekhoroshev stability, the one well described by the KAM theory, and the unstable chaotic regime in which diffusion in phase space can be detected over time spans much shorter than the age of the solar system.     

The resonant structure of Jupiter's Trojan asteroids – I. Long-term stability and diffusion

We study the global dynamics of the jovian Trojan asteroids by means of the frequency map analysis. We find and classify the main resonant structures that serve as skeleton of the phase space near the Lagrangian points. These resonances organize and control the long-term dynamics of the Trojans. Besides the secondary and secular resonances, that have already been found in other asteroid sets in mean motion resonance (e.g. main belt, Kuiper belt), we identify a new type of resonance that involves secular frequencies and the frequency of the great inequality, but not the libration frequency. Moreover, this new family of resonances plays an important role in the slow transport mechanism that drives Trojans from the inner stable region to eventual ejections. Finally, we relate this global view of the dynamics with the observed Trojans, identify the asteroids that are close to these resonances and study their long-term behaviour.     

The Fragmentation of Asteroids – A Synopsis of Analytical Methods (Mitteilungen der Universitäts-Sternwarte Jena Nr. 119)

Basing on the author's work a review of the possibilities as well as the limits of treating the problem of the collisional history of the asteroids by analytical methods is given. Using empirical data on rock fragmentation and general principles like symmetry and mass conservation the distribution function of the fragments arising from a single collision is analytically formulated. The size distribution of asteroids adjusting when crushing collisions have taken place a sufficiently long time can be obtained as the solution of an integrodifferential equation with partial derivatives (equation of fragmentation). Quasi-stationary solutions of the equation of fragmentation are discussed for particular cases. The problem of the steady state is reduced to the solution of a transcendental equation. The results obtained show that analytical methods already offer a good theoretical understanding of the observed size distribution of the asteroids. They should be, therefore, a useful basis of carrying out numerical calculations.     

The collisional and dynamical evolution of the main-belt and NEA size distributions

The size distribution of main belt of asteroids is determined primarily by collisional processes. Large asteroids break up and form smaller asteroids in a collisional cascade, with the outcome controlled by the strength-size relationship for asteroids. In addition to collisional processes, the non-collisional removal of asteroids from the main belt (and their insertion into the near-Earth asteroid (NEA) population) is critical, and involves several effects: strong resonances increase the orbital eccentricity of asteroids and cause them to enter the inner planet region; chaotic diffusion by numerous weak resonances causes a slow leak of asteroids into the Mars- and Earth-crossing populations; and the Yarkovsky effect, a radiation force on asteroids, is the primary process that drives asteroids into these resonant escape routes. Yarkovsky drift is size-dependent and can modify the main-belt size distribution. The NEA size distribution is primarily determined by its source, the main-belt population, and by the size-dependent processes that deliver bodies from the main belt. All of these effects are simulated in a numerical collisional evolution model that incorporates removal by non-collisional processes. We test our model against a wide range of observational constraints, such as the observed main-belt and NEA size distributions, the number of asteroid families, the preserved basaltic crust of Vesta and its large south-pole impact basin, the cosmic ray exposure ages of meteorites, and the cratering records on asteroids. We find a strength-size relationship for main-belt asteroids and non-collisional removal rates from the main belt such that our model fits these constraints as best as possible within the parameter space we explore. Our results are consistent with other independent estimates of strength and removal rates.     

Frequency modified fourier transform and its application to asteroids

Recently a method has been suggested to analyze the chaotic behaviour of a conservative dynamical system by numerical analysis of the fundamental frequencies. Frequencies and amplitudes are determined step by step. As the frequencies are not generally orthogonal, a Gramm-Schmidt orthogonalization is made and for each new frequency the old amplitudes of previously determined frequencies are corrected. For a chaotic trajectory variations of the frequencies and amplitudes determined over different time periods are expected. The change of frequencies in such a calculation is a measure of the chaoticity of the trajectory. While amplitudes are corrected, the frequencies (once determined) are constant. We suggest here simple linear corrections of frequencies for the effect of other close frequencies. The improvement of frequency determination is demonstrated on a model case. This method is applied to the first fifty numbered asteroids.     

Asteroid lightcurve phase shift from rough‐surface shadowing

We have simulated asteroid lightcurves for simple shape models using a realistic surface scattering law. The scattering law includes a shadowing function computed with numerical ray‐tracing. We computed lightcurves in a variety of illumination geometries for both the traditional Lommel–Seeliger law and our seminumerical law. We observe a shift in the rotational phase of the lightcurves, which depends on the parameters of the scattering law as well as the illumination geometry and the direction of the spin axis of the asteroid. This phase shift is always zero at opposition, and can be as large as 10° for illumination geometries typical for Main Belt asteroids. The phase shift has implications on the accuracy of other results which are based on asteroid lightcurve analysis, such as spin‐state or shape determination.     

Asteroid observations at low phase angles. I. 50 Virginia, 91 Aegina and 102 Miriam

We present observations of magnitude-phase dependences of three low-albedo asteroids down to phase angles of 0.1–0.2°. Data were obtained during 40 nights from 1994 to 1995 within the joint observational program at ESO and Kharkiv Astronomical Observatories with the aim to reach as low phase angles as possible. All three low-albedo asteroids may display a small nonlinear increase in magnitude-phase dependence at subdegree phase angles. The phase curves of 50 Virginia and 102 Miriam are poorly approximated by the HG function. Rotation periods of the asteroids were also determined: 14.310±0.010 hours for 50 Virginia, 6.030±0.001 h for 91 Aegina and 15.789± 0.003 h for 102 Miriam.     

On the 3/1 planar and circular resonant problem

The simplest model of a resonant problem of second order is the planar and circular case. Simplification like this is very old and for 3/1 resonance, several authors have studied this problem with different purposes. In this work, we test this model for the available asteroids, by applying Hori's perturbation method. Explicit solutions of the intermediate orbit are obtained. In the plane of two constants of the problem, all types of motion are described. By testing the model, it is shown that, in general, one can confirm results of numerical integrations indicating libration for a few number of asteroids and circulation for most of them. However, agreement in numerical values for amplitude and period of librations seems to be not possible mainly if Jupiter's eccentricity is neglected. On the other hand, even though there might be some physical reasons determining that only asteroids with high eccentricity may librate, it is shown that, from mathematical point of view, libration may occur even in the case of small eccentricities provided that some relations are satisfied.     

The reducing transformation and Apocentric Librators

We propose a canonical transformation reducing the averaged planar planetary problem near resonance to a one degree of freedom problem when the perturbation is truncated at the first order in the eccentricities.This reducing transformation leads to a very simple explanation of the puzzling behaviour of the Apocentric Librators, a class of asteroids identified by Franklinet al. (1975).An exploration of the phase space of the average problem with the use of the mapping technique shows that the alternation of two libration mechanism is a common feature for initial conditions near, but not inside, the deep resonance region.     

Atmospheric shock waves after impacts of cosmic bodies up to 1000 m in diameter

The results of numerical modeling of meteoroids' interaction with Earth's atmosphere are presented. We model the entry in two dimensions and then interpolate the results into a 3‐D model to calculate interaction of shock waves with the surface. Maximum shock pressures, wind speeds, and areas subjected to substantial overpressure are calculated for oblique impacts of asteroids and comets. We show that vertical impacts produce a smaller damage zone on the surface than oblique ones. Damage caused by shock waves covers an order of magnitude larger area than any other hazardous effects. The function of energy release in the atmosphere, which is traditionally used in meteoritics, has a limited application if cosmic bodies are larger than tens of meters in diameter: at each time moment energy is smoothed along a substantial length of the trajectory; both emitted radiation (routinely used for calibration of semi‐analytical models) and shock wave amplitude are complex functions of temperature–density distributions in atmosphere.     

The main belt as a source of near-Earth asteroids

We investigate the flux of main-belt asteroid fragments into resonant orbits converting them into near-Earth asteroids (NEAs), and the variability of this flux due to chance interasteroidal collisions. A numerical model is used, based on collisional physics consistent with the results of laboratory impact experiments. The assumed main-belt asteroid size distribution is derived from that of known asteroids extrapolated down to sizes of ≈ 40 cm, modified in such a way to yield a quasi-stationary fragment production rate over times ≈ 100 Myr. The results show that the asteroid belt can supply a few hundred km-sized NEAs per year, well enough to sustain the current population of such bodies. On the other hand, if our collisional physics is correct, the number of existing 10-km objects implies that these objects either have very long-lived orbits, or must come from a different source (i.e., comets). Our model predicts that the fragments supplied from the asteroid belt have initially a power-law size distribution somewhat steeper than the observed one, suggesting preferential removal of small objects. The component of the NEA population with dynamical lifetimes shorter than or of the order of 1 Myr can vary by a factor reaching up to a few tens, due to single large-scale collisions in the main belt; these fluctuations are enhanced for smaller bodies and faster evolutionary time scales. As a consequence, the Earth's cratering rate can also change by about an order of magnitude over the 0.1 to 1 Myr time scales. Despite these sporadic spikes, when averaged over times of 10 Myr or longer the fluctuations are unlikely to exceed a factor two.     

Planet crossing asteroids and parallel computing: Project spaceguard

The orbits of the asteroids crossing the orbit of the Earth and other planets are chaotic and cannot be computed in a deterministic way for a time span long enough to study the probability of collisions. It is possible to study the statistical behaviour of a large number of such orbits over a long span of time, provided enough computing resources and intelligent post processing software are available. The former, problem can be handled by exploiting the enormous power of parallel computing systems. The orbit of the asteroids can be studied as a restricted (N+M)-body problem which is suitable for the use of parallel processing, by using one processor to compute the orbits of the planets and the others to compute the orbits the asteroids. This scheme has been implemented on LCAP-2, an array of IBM and FPS processors with shared memory designed by E. Clementi (IBM). The parallelisation efficiency has been over 80%, and the overall speed over 90 MegaFLOPS; the orbits of all the asteroids with perihelia lower than the aphelion of Mars (410 objects) have been computed for 200,000, years (Project SPACEGUARD). The most difficult step of the project is the post processing of the very large amount of output data and to gather qualitative information on the behaviour of so many orbits without resorting to the traditional technique, i.e. human examination of output in graphical form. Within Project SPACEGUARD we have developed a qualitative classification of the orbits of the planet crossers. To develop an entirely automated classification of the qualitative orbital behaviour-including crossing behaviour, resonances (mean motion and secular), and protection mechanisms avoiding collisions-is a challenge to be met.     

Collisional evolution of the early asteroid belt

We present numerical results obtained by a simulation of the collisional process between asteroids and scattered comets from the Uranus–Neptune zone. This mechanism allows the use of single exponent incremental size distributions for the initial belt reaching a final distribution that matches the observed population very well. Since the cometary bombardment was extremely efficient removing mass from the primordial asteroid belt in a very short time, we always obtained belts with total masses less than 0.001 M_⊕ after ≈ 2×10⁷ yrs. This result allows processes with an important initial mass preserving Vestas basaltic crust.     

Frequency analysis of the stability of asteroids in the framework of the restricted three-body problem

The stability of some asteroids, in the framework of the restricted three-body problem, has been recently proved in (Celletti and Chierchia, 2003), by developing an isoenergetic KAM theorem. More precisely, having fixed a level of energy related to the motion of the asteroid, the stability can be obtained by showing the existence of nearby trapping invariant tori existing at the same energy level. The analytical results are compatible with the astronomical observations, since the theorem is valid for the realistic mass-ratio of the primaries. The model adopted in (Celletti and Chierchia, 2003), is the planar, circular, restricted three-body model, in which only the most significant contributions of the Fourier development of the perturbation are retained. In this paper we investigate numerically the stability of the same asteroids considered in (Celletti and Chierchia, 2003), (namely, Iris, Victoria and Renzia). In particular, we implement the nowadays standard method of frequency-map analysis and we compare our investigation with the analytical results on the planar, circular model with the truncated perturbing function. By means of frequency analysis, we study the behaviour of the bounding tori and henceforth we infer stability properties on the dynamics of the asteroids. In order to test the validity of the truncated Hamiltonian, we consider also the complete expression of the perturbing function on which we perform again frequency analysis. We investigate also more realistic models, taking into account the eccentricity of the trajectory of Jupiter (planar-elliptic problem) or the relative inclination of the orbits (circular-inclined model). We did not find a relevant discrepancy among the different models.