Terrestrial planet and asteroid formation in the presence of giant planets. I. Relative velocities of planetesimals subject to Jupiter and Saturn perturbations.

We investigate the orbital evolution of 10(13)- to 10(25) -g planetesimals near 1 AU and in the asteroid belt (near 2.6 AU) prior to the stage of evolution when the mutual perturbations between the planetesimals become important. We include nebular gas drag and the effects of Jupiter and Saturn at their present masses and in their present orbits. Gas drag introduces a size-dependent phasing of the secular perturbations, which leads to a pronounced dip in encounter velocities (Venc) between bodies of similar mass. Plantesimals of identical mass have Venc approximately 1 and approximately 10 m s-1 (near 1 and 2.6 AU, respectively) while bodies differing by approximately 10 in mass have Venc approximately 10 and approximately 100 m s-1 (near 1 and 2.6 AU, respectively). Under these conditions, growth, rather than erosion, will occur only by collisions of bodies of nearly the same mass. There will be essentially no gravitational focusing between bodies less than 10(22) to 10(25) g, allowing growth of planetary embryos in the terrestrial planet region to proceed in a slower nonrunaway fashion. The environment in the asteroid belt will be even more forbidding and it is uncertain whether even the severely depleted present asteroid belt could form under these conditions. The perturbations of Jupiter and Saturn are quite sensitive to their semi-major axes and decrease when the planets' heliocentric distances are increased to allow for protoplanet migration. It is possible, though not clearly demonstrated, that this could produce a depleted asteroid belt but permit formation of a system of terrestrial planet embryos on a approximately 10(6)-year timescale, initially by nonrunaway growth and transitioning to runaway growth after approximately 10(5) years. The calculations reported here are valid under the condition that the relative velocities of the bodies are determined only by Jupiter and Saturn perturbations and by gas drag, with no mutual perturbations between planetesimals. If, while subject to these conditions, the bodies become large enough for their mutual perturbations to influence their velocity and size evolution significantly, the problem becomes much more complex. This problem is under investigation.     

Theoretical profile of Saturn's rings

Steady-state solutions for the optical thickness of Saturn's rings are studied in terms of Hämeen-Anttila's (1983) theory of bimodal gravitating systems. The elastic properties of particles determine the behaviour of the rarefied mode (gaps), while the dense mode (ringlets) depends on the size and the internal density of the particles. In the outer parts of the rings the dense mode is unstable against the growth of gravitational perturbations. Inside the Roche distance this produces only very narrow ring-shaped configurations with helical orbits around them, and the system is not destroyed. The outer boundary of the rings corresponds to the distance beyond which the gravitational instability transforms the dense mode into strictly local condensations (moons). The inner boundary of the ring system is caused by the absence of dense mode near Saturn. The rarefied mode is stable in a larger region.     

Stability of the solar system: Evidence from the asteroids

An analysis of the distribution of the orbital periods of the asteroids has shown that there is a preference for these periods to be near-commensurate with that of Mars. We suggest that this preference is associated with a formation process and implies that the orbital period of Mars has not changed greatly since the time of asteroid formation. We deduce from this that the solar system is highly stable and long-period gravitational perturbations have probably had little influence on the gross evolution of the solar system.Paper presented to the NATO Advanced Study Institute on Lunar Studies, Patras, Greece, September 1971.     

A former asteroidal planet as the origin of comets

The idea of a missing planet between Mars and Jupiter has been with us since the formulation of the Titius-Bode law. The discovery of the asteroid belt in that location led to speculation about a planetary breakup event. Both ideas remained conjectures until Ovenden's finding in 1972, from which it could be derived that the mass of the missing planet was about 90 Earth masses and that its breakup was astronomically recent. Apparently much of that mass was blown out of the solar system during the disruption of the planet. Because of the action of planetary perturbations, only two types of orbits of surviving fragments could remain at present-asteroid orbits and once-around very-long-period elliptical orbits. Objects in the latter type of orbit are known to exist-the very-long-period comets. A large number of these are on elliptical trajectories with periods of revolution of 5 million years; yet they are known to have made no more than one revolution in an orbit passing close to the Sun. By direct calculation it is possible to predict the distribution of the orbital elements of objects moving on long-period ellipses which might have originated in a breakup event in the asteroid belt 5 million years ago. The comet orbits have the predicted distribution in every case where a measure is possible. Some of the distribution anomalies, such as a bias in the directions of perihelion passage, are statistically strong and would be difficult to explain in any other uncontrived way. In addition, a relative deficiency of orbits with perihelia less than 1 AU indicates that the comets must have had small perihelion distances since their origin, rather than that they have been perturbed into small perihelion orbits from a distant “cloud” of comets by means of stellar encounters. The comet orbital data lead to the conclusion that all comets originated in a breakup event in the asteroid belt (5.5±0.6) × 10⁶ years ago. Asteroid and meteoritic evidence can now be interpreted in a way which not only is supportive but also provides fresh insights into understanding their physical, chemical, and dynamical properties. Particularily noteworthy are the young cosmic-ray exposure ages of meteorites, evidence of a previous high-temperature/pressure environment and of chemical differentiation of the parent body, and compositional similarities among comets, asteroids, and meteorites. Certain “explosion signatures” in asteroid orbital element distributions are likewise indicative. Tektites may also have originated in the same event; but if so, there are important implications regarding the absolute accuracy of certain geological dating methods. Little is known about possible planetary breakup mechanisms of the requisite type, though some speculations are offered. In any case, the asteroid belt is an existing fact; and the arguments presented here that a large planet did disintegrate 5 million years ago must be judged on their merits, even in the absence of a suitable theory of planetary explosions.     

The Poynting-Robertson effect in the Newtonian potential with a Yukawa correction

We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.     

A preliminary analysis of the orbit of the Mars Trojan asteroid (5261) Eureka

Observations and results of orbit determination of the first known Mars Trojan asteroid (5261) Eureka are presented. We have numerically calculated the evolution of the orbital elements, and have analyzed the behavior of the motion during the next 2 Myr. Strong perturbations by planets other than Mars seem to stabilize the eccentricity of the asteroid by stirring the high order resonances present in the elliptic restricted problem. As a result, the orbit appears stable at least on megayear timescales. The difference of the mean longitudes of Mars and Eureka and the semimajor axis of the asteroid form a pair of variables that essentially behave in an adiabatic manner, while the evolution of the other orbital elements is largely determined by the perturbations due to other planets.     

Dynamical Effects of Mars on Asteroidal Dust Particles

Asteroidal dust particles resulting from family-forming events migrate from their source locations in the asteroid belt inwards towards the Sun under the effect of Poynting-Robertson (PR) drag. Understanding the distribution of these dust particle orbits in the inner solar system is of great importance to determining the asteroidal contribution to the zodiacal cloud, the accretion rate by the Earth, and the threat that these particles pose to spacecraft and satellites in near-Earth space. In order to correctly describe this distribution of orbits in the inner solar system, we must track the dynamical perturbations that the dust particle orbits experience as they migrate inwards. In a seminal paper Öpik (1951) determines that very few of the μm-cm sized dust particles suffer a collision with the planet face as they decay inwards past Mars. Here we re-analyze this problem, considering additionally the likelihood that the dust particle orbits pass through the Hill sphere of Mars (to various depths) and experience potentially significant perturbations to their orbits. We find that a considerable fraction of dust particle orbits will enter the Hill sphere of Mars. Furthermore, we find that there is a bias with inclination, particle size, and eccentricity of the particle orbits that enter the Martian Hill sphere. In particular the bias with inclination may create a bias towards higher-inclination sources in the proportions of asteroid family particles that reach near-Earth space.     

Artificial Martian frozen orbits and Sun-Synchronous orbits using continuous low-thrust control

This paper analyses three types of artificial orbits around Mars pushed by continuous low-thrust control: artificial frozen orbits, artificial Sun-Synchronous orbits and artificial Sun-Synchronous frozen orbits. These artificial orbits have similar characteristics to natural frozen orbits and Sun-Synchronous orbits, and their orbital parameters can be selected arbitrarily by using continuous low-thrust control. One control strategy to achieve the artificial frozen orbit is using both the transverse and radial continuous low-thrust control, and another to achieve the artificial Sun-Synchronous orbit is using the normal continuous low-thrust control. These continuous low-thrust control strategies consider J ₂, J ₃, and J ₄ perturbations of Mars. It is proved that both control strategies can minimize characteristic velocity. Relevant formulas are derived, and numerical results are presented. Given the same initial orbital parameters, the control acceleration and characteristic velocity taking into account J ₂, J ₃, and J ₄ perturbations are similar to those taking into account J ₂ perturbations for both Mars and the Earth. The control thrust of the orbit around Mars is smaller than that around the Earth. The magnitude of the control acceleration of ASFOM-4 (named as Artificial Sun-Synchronous Frozen Orbit Method 4) is the lowest among these strategies and the characteristic velocity within one orbital period is only 0.5219 m/s for the artificial Sun-Synchronous frozen orbit around Mars. It is evident that the relationship among the control thrusts and the primary orbital parameters of Martian artificial orbits is always similar to that of the Earth. Simulation shows that the control scheme extends the orbital parameters’ selection range of three types of orbits around Mars, compared with the natural frozen orbit and Sun-Synchronous orbit.     

Slow Diffusion and Effective Stability of Dust Particles Orbiting Asteroids

This paper presents quantitative estimates of orbital lifetimes for small dust particles in the vicinity of a spherical asteroid orbiting the Sun on a circular orbit. In particular, the modeling accounts for gravitational interactions with the asteroid and the Sun and radiation pressure from the Sun. The analysis focuses on slow drift of particles in a state-space neighborhood of weakly unstable periodic orbits corresponding to the most tightly bound orbits about the asteroid. Through a Hamiltonian normal-form computation, Nekhoroshev-type estimates are derived that yield regions in state space of effective stability, that is, regions in which the grains would remain over significant periods of time in the absence of other depletion and replenishment mechanisms.     

Determining asteroid masses from perturbations on Mars

The orbit of Mars is perturbed more than 5 m, a value compatible with the accuracy of the Viking lander ranging data, by about three dozen asteroids. In addition to larger asteroids throughout the belt, significant perturbations of long period are generated by smaller objects near commensurabilities with Mars. The largest periodic terms induced by 1 Ceres and 2 Pallas have amplitudes of 0.8 and 0.2 km, respectively, both with 10-year periods. Due to a near commensurability, 4 Vesta causes a 5-km, 52-year term. While the Viking ranges will yield significant mass determinations for the largest three asteroids, and some of the smaller bodies should be detectable, it will be difficult to seperate the smaller bodies with useful accuracies. Accurate discrimination must await range data from future missions to Mars or other bodies in the neighborhood of the asteroid belt. The Viking ranges can also yield improved masses for the outer planets (except Pluto), an application which is being exploited by groups analyzing these data. Uncertainties in the asteroid masses limit the ultimate accuracy of the Viking determinations of both the long time scale motion of the system the inner four planets with respect to an inertial frame and the rate of change of the gravitational constant.     

The Use of Lie Series in the Construction of a Perturbation Theory and Some Recent Results in the Theory of the Motion of Hyperion

This paper begins with a brief review of a form of the Lie series transformation, and then reports some new results in the study, using Lie series methods, of the orbit of Saturn's satellite Hyperion. In particular, improved expressions are given for the long-period perturbations of the orbital elements which describe the motion in the orbit plane, and also first results for expressions for the short-period perturbations in the apse longitude, derived from the Lie series generating function. This revised version was published online in July 2006 with corrections to the Cover Date.     

Precession of the orbital nodes of Jupiter and Saturn triggered by the mutual perturbation: A model of two rings

The problem of the precession of the orbital planes of Jupiter and Saturn under the influence of mutual gravitational perturbations was formulated and solved using a simple dynamical model. Using the Gauss method, the planetary orbits are modeled by material circular rings, intersecting along the diameter at a small angle α. The planet masses, semimajor axes and inclination angles of orbits correspond to the rings. What is new is that each ring has an angular momentum equal to the orbital angular momentum of the planet. Contrary to popular belief, it was proved that the orbital resonance 5: 2 does not preclude the use of the ring model. Moreover, the period of averaging of the disturbing force (T ≈ 1332 yr) proves to be appreciably greater than a conventionally used period (≈900 yr). The mutual potential energy of rings and the torque of gravitational forces between the rings were calculated. We compiled and solved the system of differential equations for the spatial motion of rings. It was established that a perturbing torque causes the precession and simultaneous rotation of the orbital planes of Jupiter and Saturn. Moreover, the opposite orbit nodes on the Laplace plane coincide and perform a secular movement in retrograde direction with the same velocity of 25.6″/yr and the period T _J = T _S ≈ 50687 yr. These results are close to those obtained in the general theory (25.93″/yr), which confirms the adequacy of the developed model. It was found that the vectors of the angular velocity of orbital rings move counterclockwise over circular cones and describe circles on the celestial sphere with radii β₁ ≈ 0.8403504° (Saturn) and β₂ ≈ 0.3409296° (Jupiter) around the point which is located at an angular distance of 1.647607° from the ecliptic pole.     

Model orbits of asteroid 3040 Kozai

The orbital evolutions of the asteroid 3040 Kozai and model asteroids with similar orbits have been investigated. Their osculating orbits for an epoch 1991 December 10 were numerically integrated forward within the interval of 20,000 years, using a dynamical model of the solar system consisting of all inner planets, Jupiter, and Saturn.The orbit of the asteroid Kozai is stable. Its motion is affected only by long-period perturbations of planets. With change of the argument of perihelion of the asteroid Kozai, the evolution of the model asteroid orbits changes essentially, too. The model orbits with the argument of perihelion changed by the order of 10% show that asteroids with such orbital parameters may approach the Earth orbit, while asteroids with larger changes may even cross it, at least after 10,000 years. Long-term orbital evolution of asteroids with these orbital parameters is very sensitive on their angular elements.     

Computer simulation of the formation of asteroid belt structure

The system of two gravitational centers with variable separation between components one of which (the primary) loses its mass onto another (the secondary) is investigated under condition of total mass and angular momentum conservation. When the primary/secondary mass ratio becomes about that of Jupiter/Sun the small bodies ejected with the gaseous matter through the inner Lagrange point from the Roche lobe of the primary form a ring similar to the asteroid belt of the solar system. The formation of ring structure is calculated by numerical integration of Newtonian equations of N-body problem in orbital plane of the gravitational centers. The results are compared with the planar subsystem of the asteroid belt. The presence of the main gaps in the distribution of their mean motions at 2/1, 3/1, 5/2 and some other commensurabilities with the primary mean motion is found. More fine details of the belt structure are obtained, e.g. the gap asymmetry and a qualitative agreement with the eccentricity distribution. Within the scope of the same model the external part of the ring is investigated all the pairwise interactions being included. The clustering of bodies near 3/2 commensurability isolated from the main belt by the wide gap centered at 5/3 commensurability is obtained. It is supposed that the ring structure and the interplanetary spacing law for the terrestrial planets are due to the same mechanism.     

On asteroid distribution

The existing explanations for the asteroid distribution in the main belt (between the orbits of Mars and Jupiter) are based on numerical integration of resonance orbits in models with more than two degrees of freedom. We suggest an approach based on the investigation of the families of periodic solutions of the planar circular restricted three-body problem, i.e., a model with two degrees of freedom. This work shows that (a) the distribution of asteroids near the (p + 1)/p resonances and position of the outer boundary of the main asteroid belt can be explained within the planar circular restricted three-body problem and (b) this problem does not explain the asteroid distribution near other resonances.     

Approximate analytic method for high-apogee twelve-hour orbits of artificial Earth’s satellites

We propose an approach to the study of the evolution of high-apogee twelve-hour orbits of artificial Earth’s satellites. We describe parameters of the motion model used for the artificial Earth’s satellite such that the principal gravitational perturbations of the Moon and Sun, nonsphericity of the Earth, and perturbations from the light pressure force are approximately taken into account. To solve the system of averaged equations describing the evolution of the orbit parameters of an artificial satellite, we use both numeric and analytic methods. To select initial parameters of the twelve-hour orbit, we assume that the path of the satellite along the surface of the Earth is stable. Results obtained by the analytic method and by the numerical integration of the evolving system are compared. For intervals of several years, we obtain estimates of oscillation periods and amplitudes for orbital elements. To verify the results and estimate the precision of the method, we use the numerical integration of rigorous (not averaged) equations of motion of the artificial satellite: they take into account forces acting on the satellite substantially more completely and precisely. The described method can be applied not only to the investigation of orbit evolutions of artificial satellites of the Earth; it can be applied to the investigation of the orbit evolution for other planets of the Solar system provided that the corresponding research problem will arise in the future and the considered special class of resonance orbits of satellites will be used for that purpose.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets - Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model”. We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size.     

Orbital evolution of the Gefion and Adeona asteroid families: close encounters with massive asteroids and the Yarkovsky effect

Asteroid families are groupings of minor planets identified by clustering in their proper orbital elements; these objects have spectral signatures consistent with an origin in the break-up of a common parent body. From the current values of proper semimajor axes a of family members one might hope to estimate the ejection velocities with which the fragments left the putative break-up event (assuming that the pieces were ejected isotropically). However, the ejection velocities so inferred are consistently higher than N-body and hydro-code simulations, as well as laboratory experiments, suggest. To explain this discrepancy between today’s orbital distribution of asteroid family members and their supposed launch velocities, we study whether asteroid family members might have been ejected from the collision at low speeds and then slowly drifted to their current positions, via one or more dynamical processes. Studies show that the proper a of asteroid family members can be altered by two mechanisms: (i) close encounters with massive asteroids, and (ii) the Yarkovsky non-gravitational effect. Because the Yarkovsky effect for kilometer-sized bodies decreases with asteroid diameter D, it is unlikely to have appreciably moved large asteroids (say those with D > 15 km) over the typical family age (1-2 Gyr).For this reason, we numerically studied the mobility of family members produced by close encounters with main-belt, non-family asteroids that were thought massive enough to significantly change their orbits over long timescales. Our goal was to learn the degree to which perturbations might modify the proper a values of all family members, including those too large to be influenced by the Yarkovsky effect. Our initial simulations demonstrated immediately that very few asteroids were massive enough to significantly alter relative orbits among family members. Thus, to maximize gravitational perturbations in our 500-Myr integrations, we investigated the effect of close encounters on two families, Gefion and Adeona, that have high encounter probabilities with 1 Ceres, by far the largest asteroid in the main belt. Our results show that members of these families spreads in a of less than 5% since their formation. Thus gravitational interactions cannot account for the large inferred escape velocities.The effect of close encounters with massive asteroids is, however, not entirely negligible. For about 10% of the simulated bodies, close encounters increased the “inferred” ejection velocities from sub-100 m/s to values greater than 100 m/s, beyond what hydro-code and N-body simulations suggest are the maximum possible initial ejection velocity for members of Adeona and Gefion with D > 15 km. Thus this mechanism of mobility may be responsible for the unusually high inferred ejection speeds of a few of the largest members of these two families.To understand the orbital evolution of the entire family, including smaller members, we also performed simulations to account for the drift of smaller asteroids caused by the Yarkovsky effect. Our two sets of simulations suggest that the two families we investigated are relatively young compared to larger families like Koronis and Themis, which have estimated ages of about 2 Byr. The Adeona and Gefion families seems to be no more than 600 and 850 Myr old, respectively.     

Orbital history of the Martian satellites with inferences on their origin

Recent Viking results indicate the Martian satellites are composed of carbonaceous chondritic material, suggesting that Phobos and Deimos were once asteroids captured by Mars. On the other hand, the low eccentricities and inclinations of their orbits on the equator of Mars argue against that hypothesis. This paper presents detailed calculations of the tidal evolution of Phobos and Deimos, considering dissipation in both Mars and its satellites simultaneously and using a new method applicable for any value of the eccentricity. In particular, including precession of the satellites' orbits indicates that they have always remained close to their Laplacian plane, so that the orbital planes of Phobos and Deimos switched from near the Martian orbital plane to the Martian equator once the perturbations due to the planetary oblateness dominated the solar perturbations, as they do presently. The results show that Deimos has been little affected by tides, but several billion (10⁹) years ago, Phobos was in a highly eccentric orbit lying near the common plane of the solar system. This outcome is obtained for very reasonable values of dissipation inside Mars and inside Phobos. Implications for the origin of the Martian satellites are discussed.     

Problems of origin of asteroids and comets

Abstract— Various hypotheses of the origin of asteroids and comets are briefly discussed. Interaction of planetesimals in the asteroid zone (AZ) with the gas, their perturbations by proto-Jupiter, and sweeping them out by more massive Jupiter zone bodies when they penetrated the AZ are considered. If the gas was turbulent, it could prevent a settling of dust particles to the equatorial plane of the disk and formation of dust condensations due to gravitational instability. Then particles grew by sticking upon collision. Gas moved radially due to turbulent viscosity and its dissipation. Small particles moved more-or-less together with the gas. As a result of gas drag, larger particles and bodies moved relative to the gas in the direction of increasing gas pressure. Gas would remove much of the solid material from the AZ if most bodies larger than a few km disintegrated by collisions into fragments smaller than a few tens of meters. Most of these fragments would then move into the Martian zone, and the small mass of Mars would have no explanation. Resonant perturbations of asteroids by Jupiter are discussed. In the model of a small mass disk they could scan through the asteroid belt due to changes in Jupiter's distance from the Sun that occurred when this planet accreted the gas and ejected the bodies from the solar system. Such a scanning considerably accelerated the removal of asteroids from the AZ. Massive Jupiter zone bodies with large orbital eccentricities that crossed the AZ were probably efficient at sweeping out bodies. Larger bodies increased the random velocities of the remaining asteroids at close encounters to the present values ～ 5 km/s. Restrictions on the runaway growth of giant planets, on the relative velocities of bodies and the disk surface density that follow from the consideration of the origin of the asteroid belt and the cometary cloud are considered.