We present preliminary results of terrestrial planet formation using on the one hand classical numerical integration of hundreds of small bodies on CPUs and on the other hand -- for comparison reasons -- the results of our GPU code with thousands of
small bodies which then merge to larger ones. To be able to determine the outcome of collision events we use our smooth particle hydrodynamics (SPH) code which tracks how water is lost during such events.
As part of a national scientific network Pathways to Habitability the formation of planets and the delivery of water onto these planets is a key question as water is essential for the development of life. In the first part of the paper we summarize t
he state of the art of planet formation - which is still under debate in the astronomical community - before we show our results on this topic. The outcome of our numerical simulations depends a lot on the choice of the initial distribution of planetesimals and planetary embryos after gas disappeared in the protoplanetary disk. We also take into account that some of these planetesimals of sizes in the order of the mass of the Moon already contained water; the quantity depends on the distance from the Sun - close-by bodies are dry, but starting from a distance of about 2 AU they can contain substantial amounts of water. We assume that the gas giants and terrestrial planets are already formed when we check the collisions of the small bodies containing water (in the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the crust and on the surface of Earth. In the second part we discuss in more detail how the formation of larger bodies after a collision may happen as the outcome depends on parameters like collision velocity, impact angle, and the materials involved. We present results obtained by SPH (Smooth Particle Hydrodynamics) simulations. We briefly describe this method and show different scenarios with respect to the formed bodies, possible fragmentation and the water content before and after the collision. In an appendix we discuss detection methods for extrasolar planets (close to 2000 such objects have been discovered so far).
We investigate the outcome of collisions of Ceres-sized planetesimals composed of a rocky core and a shell of water ice. These collisions are not only relevant for explaining the formation of planetary embryos in early planetary systems, but also pro
vide insight into the formation of asteroid families and possible water transport via colliding small bodies. Earlier studies show characteristic collision velocities exceeding the bodies mutual escape velocity which - along with the distribution of the impact angles - cover the collision outcome regimes partial accretion, erosion, and hit-and-run leading to different expected fragmentation scenarios. Existing collision simulations use bodies composed of strengthless material; we study the distribution of fragments and their water contents considering the full elasto-plastic continuum mechanics equations also including brittle failure and fragmentation.
We investigate the distribution of encounter velocities and impact angles describing collisions in the habitable zone of the early planetary system. Here we present a catalogue of collision characteristics for a particular mass ratio of the colliding
bodies and seven different planetesimal masses ranging from a tenth of Ceres mass to 10 times the mass of the Moon. We show that there are virtually no collisions with impact speeds lower than the surface escape velocity and a similar velocity-impact angle distribution for different planetesimal masses if velocities are normalized using the escape velocity. An additional perturbing Jupiter-like object distorts the collision velocity and impact picture in the sense that grazing impacts at higher velocities are promoted if the perturbers orbit is close to the habitable zone whereas a more distant perturber has more the effect of a mere widening of the velocity dispersion.
In a previous paper, we have presented a global view of the stability of Neptune Trojan (NT hereafter) on inclined orbit. We discuss in this paper the dependence of stability of NT orbits on the eccentricity. High-resolution dynamical maps are constr
ucted using the results of extensive numerical integrations of orbits initialized on the fine grids of initial semimajor axis (a0) versus eccentricity (e0). The extensions of regions of stable orbits on the (a0, e0) plane at different inclinations are shown. The maximum eccentricities of stable orbits in three most stable regions at low (0, 12deg.), medium (22,36deg.) and high (51, 59deg.) inclination, are found to be 0.10, 0.12 and 0.04, respectively. The fine structures in the dynamical maps are described. Via the frequency analysis method, the mechanisms that portray the dynamical maps are revealed. The secondary resonances, concerning the frequency of the librating resonant angle and the frequency of the quasi 2:1 mean motion resonance between Neptune and Uranus, are found deeply involved in the motion of NTs. Secular resonances are detected and they also contribute significantly to the triggering of chaos in the motion. Particularly, the effects of the secular resonance v8, v18 are clarified. We also investigate the orbital stabilities of six observed NTs by checking the orbits of hundreds clones of them generated within the observing error bars. We conclude that four of them, except 2001 QR322 and 2005 TO74, are deeply inside the stable region. The 2001 QR322 is in the close vicinity of the most significant secondary resonance. The 2005 TO74 locates close to the boundary separating stable orbits from unstable ones, and it may be influenced by a secular resonance.
The stability of Trojan type orbits around Neptune is studied. As the first part of our investigation, we present in this paper a global view of the stability of Trojans on inclined orbits. Using the frequency analysis method based on the FFT techniq
ue, we construct high resolution dynamical maps on the plane of initial semimajor axis $a_0$ versus inclination $i_0$. These maps show three most stable regions, with $i_0$ in the range of $(0^circ,12^circ), (22^circ,36^circ)$ and $(51^circ,59^circ)$ respectively, where the Trojans are most probably expected to be found. The similarity between the maps for the leading and trailing triangular Lagrange points $L_4$ and $L_5$ confirms the dynamical symmetry between these two points. By computing the power spectrum and the proper frequencies of the Trojan motion, we figure out the mechanisms that trigger chaos in the motion. The Kozai resonance found at high inclination varies the eccentricity and inclination of orbits, while the $ u_8$ secular resonance around $i_0sim44^circ$ pumps up the eccentricity. Both mechanisms lead to eccentric orbits and encounters with Uranus that introduce strong perturbation and drive the objects away from the Trojan like orbits. This explains the clearance of Trojan at high inclination ($>60^circ$) and an unstable gap around $44^circ$ on the dynamical map. An empirical theory is derived from the numerical results, with which the main secular resonances are located on the initial plane of $(a_0,i_0)$. The fine structures in the dynamical maps can be explained by these secular resonances.
The dynamics of artificial asteroids on the Trojan-like orbits around Neptune is investigated in this paper. We describe the dependence of the orbital stability on the initial semimajor axis a and inclination i by constructing a dynamical map on the
(a,i)-plane. Rich details are revealed in the dynamical map, especially a unstable gap at i=45 deg is determined and the mechanism triggering chaos in this region is figured out. Our investigation can be used to guide the observations.
