Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals     

J.A. Fernández  W.-H. Ip 《Icarus》1984,58(1):109-120

The final stage of the accretion of Uranus and Neptune is numerically investigated. The four Jovian planets are considered with Jupiter and Saturn assumed to have reached their present sizes, whereas Uranus and Neptune are taken with initial masses 0.2 of their present ones. Allowance is made for the orbital variation of the Jovian planets due to the exchange of angular momentum with interacting bodies (“planetesimals”). Two possible effects that may have contributed to the accretion of Uranus and Neptune are incorporated in our model: (1) an enlarged cross section for accretion of incoming planetesimals due to the presence of extended gaseous envelopes and/or circumplanetary swarms of bodies; and (2) intermediate protoplanets in mid-range orbits between the Jovian planets. Significant radial displacements are found for Uranus and Neptune during their accretion and scattering of planetesimals. The orbital angular momentum budgets of Neptune, Uranus, and Saturn turn out to be positive; i.e., they on average gain orbital angular momentum in their interactions with planetesimals and hence they are displaced outwardly. Instead, Jupiter as the main ejector of bodies loses orbital angular momentum so it moves sunward. The gravitational stirring of planetesimals caused by the introduction of intermediate protoplanets has the effect that additional solid matter is injected into the accretion zones of Uranus and Neptune. For moderate enlargements of the radius of the accretion cross section (2–4 times), the accretion time scale of Uranus and Neptune are found to be a few 10⁸ years and the initial amount of solid material required to form them of a few times their present masses. Given the crucial role played by the size of the accretion cross section, questions as to when Uranus and Neptune acquired their gaseous envelopes, when the envelopes collapsed onto the solid cores, and how massive they were are essential in defining the efficiency and time scale of accretion of the two outer Jovian planets.  相似文献   

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets     

Ryuji Morishima  Joachim Stadel 《Icarus》2010,207(2):517-535

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.  相似文献   

Asteroid (21) Lutetia as a remnant of Earth’s precursor planetesimals     

P. Vernazza  P. Lamy  O. Groussin  T. Hiroi  L. Jorda  P.L. King  M.R.M. Izawa  F. Marchis  M. Birlan  R. Brunetto 《Icarus》2011,216(2):650-659

Isotopic and chemical compositions of meteorites, coupled with dynamical simulations, suggest that the main belt of asteroids between Mars and Jupiter contains objects formed in situ as well as a population of interlopers. These interlopers are predicted to include the building blocks of the terrestrial planets as well as objects that formed beyond Neptune ( [Bottke et al., 2006] , [Levison et al., 2009] and [Walsh et al., 2011] ). Here we report that the main belt asteroid (21) Lutetia – encountered by the Rosetta spacecraft in July 2010 – has spectral (from 0.3 to 25 μm) and physical (albedo, density) properties quantitatively similar to the class of meteorites known as enstatite chondrites. The chemical and isotopic compositions of these chondrites indicate that they were an important component of the formation of Earth and other terrestrial planets. This meteoritic association implies that Lutetia is a member of a small population of planetesimals that formed in the terrestrial planet region and that has been scattered in the main belt by emerging protoplanets (Bottke et al. 2006) and/or by the migration of Jupiter (Walsh et al. 2011) early in its history. Lutetia, along with a few other main-belt asteroids, may contains part of the long-sought precursor material (or closely related materials) from which the terrestrial planets accreted.  相似文献   

Orbital migrations in planetesimal discs: N-body simulations and the resonant dynamical friction     

R. G. Cionco †  A. Brunini ‡ 《Monthly notices of the Royal Astronomical Society》2002,334(1):77-86

We have performed N -body numerical simulations of the exchange of angular momentum between a massive planet and a 3D Keplerian disc of planetesimals. Our interest is directed at the study of the classical analytical expressions of the lineal theory of density waves, as representative of the dynamical friction in discs 'dominated by the planet' and the orbital migration of the planets with regard to this effect. By means of a numerical integration of the equations of motion, we have carried out a set of numerical experiments with a large number of particles  ( N ≥10 000)  , and planets with the mass of Jupiter, Saturn and one core mass of the giant planets in the Solar system  ( M _c=10 M_⊕)  . The torque, measured in a phase in which a 'steady forcing' is clearly measurable, yields inward migration in a minimum-mass solar disc  (Σ∼10 g cm^-2  ), with a characteristic drift time of ∼ a few 10⁶ yr. The planets predate the disc, but the orbital decay rate is not sufficient to allow accretion in a time-scale relevant to the formation of giant planets. We found reductions of the measured torque on the planet, with respect to the linear theory, by a factor of 0.38 for M _c, 0.04 for Saturn and 0.01 for Jupiter, due to the increase in the perturbation on the disc. The behaviour of planets whose mass is larger than M _c is similar to the one of type II migrators in gaseous discs. Our results suggest that, in a minimum mass, solar planetesimals disc, type I migrations occur for masses smaller than M _c, whereas for this mass value it could be a transition zone between the two types of migration.  相似文献   

Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn     

Junko Kominami  Shigeru Ida 《Icarus》2004,167(2):231-243

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.  相似文献   

Building the terrestrial planets: Constrained accretion in the inner Solar System   总被引：1，自引：0，他引：1  

Sean N. Raymond  David P. O’Brien  Nathan A. Kaib 《Icarus》2009,203(2):644-19170

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.  相似文献   

Formation of terrestrial planets in disks with different surface density profiles     

Nader Haghighipour  Othon C. Winter 《Celestial Mechanics and Dynamical Astronomy》2016,124(3):235-268

We present the results of an extensive study of the final stage of terrestrial planet formation in disks with different surface density profiles and for different orbital configurations of Jupiter and Saturn. We carried out simulations in the context of the classical model with disk surface densities proportional to \({r^{-0.5}}, {r^{-1}}\) and \({r^{-1.5}}\), and also using partially depleted, non-uniform disks as in the recent model of Mars formation by Izidoro et al. (Astrophys J 782:31, 2014). The purpose of our study is to determine how the final assembly of planets and their physical properties are affected by the total mass of the disk and its radial profile. Because as a result of the interactions of giant planets with the protoplanetary disk, secular resonances will also play important roles in the orbital assembly and properties of the final terrestrial planets, we will study the effect of these resonances as well. In that respect, we divide this study into two parts. When using a partially depleted disk (Part 1), we are particularly interested in examining the effect of secular resonances on the formation of Mars and orbital stability of terrestrial planets. When using the disk in the classical model (Part 2), our goal is to determine trends that may exist between the disk surface density profile and the final properties of terrestrial planets. In the context of the depleted disk model, results of our study show that in general, the \(\nu _5\) resonance does not have a significant effect on the dynamics of planetesimals and planetary embryos, and the final orbits of terrestrial planets. However, \(\nu _6\) and \(\nu _{16}\) resonances play important roles in clearing their affecting areas. While these resonances do not alter the orbits of Mars and other terrestrial planets, they strongly deplete the region of the asteroid belt ensuring that no additional mass will be scattered into the accretion zone of Mars so that it can maintain its mass and orbital stability. In the context of the classical model, the effects of these resonances are stronger in disks with less steep surface density profiles. Our results indicate that when considering the classical model (Part 2), the final planetary systems do not seem to show a trend between the disk surface density profile and the mean number of the final planets, their masses, time of formation, and distances to the central star. Some small correlations were observed where, for instance, in disks with steeper surface density profiles, the final planets were drier, or their water contents decreased when Saturn was added to the simulations. However, in general, the final orbital and physical properties of terrestrial planets seem to vary from one system to another and depend on the mass of the disk, the spatial distribution of protoplanetary bodies (i.e., disk surface density profile), and the initial orbital configuration of giant planets. We present results of our simulations and discuss their implications for the formation of Mars and other terrestrial planets, as well as the physical properties of these objects such as their masses and water contents.  相似文献   

Reassessing the origin of Triton     

E. Nogueira  R. Gomes 《Icarus》2011,214(1):113-130

The origin of Neptune’s large, circular but retrograde satellite Triton has remained largely unexplained. There is an apparent consensus that its origin lies in it being captured, but until recently no successful capture mechanism has been found. Agnor and Hamilton (Agnor, C.B., Hamilton, D.P. [2006]. Nature 441, 192-194) demonstrated that the disruption of a trans-neptunian binary object which had Triton as a member, and which underwent a very close encounter with Neptune, was an effective mechanism to capture Triton while its former partner continued on a hyperbolic orbit. The subsequent evolution of Triton’s post-capture orbit to its current one could have proceeded through gravitational tides (Correia, A.C.M. [2009]. Astrophys. J. 704, L1-L4), during which time Triton was most likely semi-molten (McKinnon, W.B. [1984]. Nature 311, 355-358). However, to date, no study has been performed that considered both the capture and the subsequent tidal evolution. Thus it is attempted here with the use of numerical simulations. The study by Agnor and Hamilton (Agnor, C.B., Hamilton, D.P. [2006]. Nature 441, 192-194) is repeated in the framework of the Nice model (Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461) to determine the post-capture orbit of Triton. After capture Triton is then subjected to tidal evolution using the model of Mignard (Mignard, F. [1979]. Moon Planets 20, 301-315; Mignard, F. [1980]. Moon Planets 23, 185-201). The perturbations from the Sun and the figure of Neptune are included. The perturbations from the Sun acting on Triton just after its capture cause it to spend a long time in its high-eccentricity phase, usually of the order of 10 Myr, while the typical time to circularise to its current orbit is some 200 Myr, consistent with earlier studies. The current orbit of Triton is consistent with an origin through binary capture and tidal evolution, even though the model prefers Triton to be closer to Neptune than it is today. The probability of capturing Triton in this manner is approximately 0.7%. Since the capture of Triton was at most a 50% event - since only Neptune has one, but Uranus does not - we deduce that in the primordial trans-neptunian disc there were some 100 binaries with at least one Triton-sized member. Morbidelli et al. (Morbidelli, A., Levison, H.F., Bottke, W.F., Dones, L., Nesvorný, D. [2009]. Icarus 202, 310-315) concludes there were some 1000 Triton-sized bodies in the trans-neptunian proto-planetary disc, so the primordial binary fraction with at least one Triton-sized member is 10%. This value is consistent with theoretical predictions, but at the low end. If Triton was captured at the same time as Neptune’s irregular satellites, the far majority of these, including Nereid, would be lost. This suggests either that Triton was captured on an orbit with a small semi-major axisa ? 50R_N (a rare event), or that it was captured before the dynamical instability of the Nice model, or that some other mechanism was at play. The issue of keeping the irregular satellites remains unresolved.  相似文献