N-Body Simulations of Late Stage Planetary Formation with a Simple Fragmentation Model

We present results of two-dimensional gravitationalN-body simulations of the late stage of planetary formation. This stage is characterized by the direct accretion of hundreds of lunar-sized planetesimals into planetary bodies. Our simulation code is based on the Hermite Individual Timestep integration algorithm, and gravitational interactions among all bodies are included throughout the simulations. We compare our simulation with earlier works that do not include all interactions, and we find very good agreement. A previously published collisional fragmentation model is included in our simulation to study the effects of the production of fragments on the subsequent evolution of the larger planetary bodies. It is found that for realistic two-body collisions that, according to this model, both bodies will suffer fragmentation, and that the outcome of the collision will be a relatively large core containing most of the mass and a few small fragments. We present the results of simulations that include this simple fragmentation model. They indicate that the presence of small fragments have only a small effect on the growth or orbital evolution of the large planet-sized bodies.     

Numerical Modeling of the Largest Terrestrial Meteorite Craters

Multi-ring impact basins have been found on the surfaces of almost all planetary bodies in the Solar system with solid crusts. The details of their formation mechanism are still unclear. We present results of our numerical modeling of the formation of the largest known terrestrial impact craters. The geological and geophysical data on these structures accumulated over many decades are used to place constraints on the parameters of available numerical models with a dual purpose: (i) to choose parameters in available mechanical models for the crustal response of planetary bodies to a large impact and (ii) to use numerical modeling to refine the possible range of original diameters and the morphology of partially eroded terrestrial craters. We present numerical modeling results for the Vredefort, Sudbury, Chicxulub, and Popigai impact craters and compare these results with available geological and geophysical information.     

Thermal evolution of trans‐Neptunian objects,icy satellites,and minor icy planets in the early solar system

Numerical simulations are performed to understand the early thermal evolution and planetary scale differentiation of icy bodies with the radii in the range of 100–2500 km. These icy bodies include trans‐Neptunian objects, minor icy planets (e.g., Ceres, Pluto); the icy satellites of Jupiter, Saturn, Uranus, and Neptune; and probably the icy‐rocky cores of these planets. The decay energy of the radionuclides, ²⁶Al, ⁶⁰Fe, ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th, along with the impact‐induced heating during the accretion of icy bodies were taken into account to thermally evolve these planetary bodies. The simulations were performed for a wide range of initial ice and rock (dust) mass fractions of the icy bodies. Three distinct accretion scenarios were used. The sinking of the rock mass fraction in primitive water oceans produced by the substantial melting of ice could lead to planetary scale differentiation with the formation of a rocky core that is surrounded by a water ocean and an icy crust within the initial tens of millions of years of the solar system in case the planetary bodies accreted prior to the substantial decay of ²⁶Al. However, over the course of billions of years, the heat produced due to ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th could have raised the temperature of the interiors of the icy bodies to the melting point of iron and silicates, thereby leading to the formation of an iron core. Our simulations indicate the presence of an iron core even at the center of icy bodies with radii ≥500 km for different ice mass fractions.     

Meteoritical and dynamical constraints on the growth mechanisms and formation times of asteroids and Jupiter

Thermal models and radiometric ages for meteorites show that the peak temperatures inside their parent bodies were closely linked to their accretion times. Most iron meteorites come from bodies that accreted <0.5 Myr after CAIs formed and were melted by ²⁶Al and ⁶⁰Fe, probably inside 2 AU. Rare carbon-rich differentiated meteorites like ureilites probably also come from bodies that formed <1 Myr after CAIs, but in the outer part of the asteroid belt. Chondrite groups accreted intermittently from diverse batches of chondrules and other materials over a 4 Myr period starting 1 Myr after CAI formation when planetary embryos may already have formed at ∼1 AU. Meteorite evidence precludes accretion of late-forming chondrites on the surface of early-formed bodies; instead chondritic and non-chondritic meteorites probably formed in separate planetesimals. Maximum metamorphic temperatures in chondrite groups are correlated with mean chondrule age, as expected if ²⁶Al and ⁶⁰Fe were the predominant heat sources. Because late-forming bodies could not accrete close to large, early-formed bodies, planetesimal formation may have spread across the nebula from regions where the differentiated bodies formed. Dynamical models suggest that the asteroids could not have accreted in the main belt if Jupiter formed before the asteroids. Therefore Jupiter probably reached its current mass >3-5 Myr after CAIs formed. This precludes formation of Jupiter via a gravitational instability <1 Myr after the solar nebula formed, and strongly favors core accretion. Jupiter probably formed too late to make chondrules by generating shocks directly, or indirectly by scattering Ceres-sized bodies across the belt. Nevertheless, shocks formed by gravitational instabilities or Ceres-sized bodies scattered by planetary embryos may have produced some chondrules. The minimum lifetime for the solar nebula of 3-5 Myr inferred from the total spread of CAI and chondrule ages may exceed the median lifetime of 3 Myr for protoplanetary disks, but is well within the 1-10 Myr observed range. Shorter formation times for extrasolar planets may help to explain their unusual orbits compared to those of solar giant planets.     

On the mechanism of reactions in rarefied gas: Processes in the solar system

The present study develops a previously suggested physico-chemical theory, according to which the regular structure of planetary and satellite systems is explained on the basis of notions regarding periodic condensation of gaseous matter in space and time during formation of the centre body. The author examines specific chemical transformation models that may have been determinative factors in the formation of primary rings of planetary systems. A simplified model is based on the concept that condensed iron compounds, which formed as a result of reactions between the protocloud (protodisk) matter and iron carbonyl hydrides diffusing from the neighbourhood of the centre body, were essentially the primary ring embryos. A concept has been suggested according to which optically active substances (stereoisomers) can form in protoplanetary clouds in the magnetic and gravitational fields of centre bodies.     

Star-like activity from a very young 'isolated planet'

The discovery of isolated bodies of planetary mass has challenged the paradigm that planets form only as companions to stars. To determine whether 'isolated planets', brown dwarfs and stars can have a common origin, we have made deep submillimetre observations of part of the ρ Oph B star formation region. Spectroscopy of the 9-Jupiter-mass core Oph B-11 has revealed carbon monoxide line wings such as those of a protostar. Moreover, the estimated mass of outflowing gas lies on the force versus core-mass relation for protostars and protobrown dwarfs. This is evidence for a common process that can form any object between planetary and stellar masses in a molecular cloud. In a submillimetre continuum map, six compact cores in ρ Oph B were found to have masses presently below the deuterium-burning limit, extending the core mass function down to  0.01 M_⊙  with the approximate form  d N /d M ∝ M ^−3/2  . If these lowest-mass cores are not transient and can collapse under gravity, then isolated planets should be very common in ρ Oph in the future, as is the case in the Orion star formation region. In fact, the isolated planetary objects that may form from these cores would outnumber the massive planets that have been found as companions to stars.     

Accumulation processes in the primitive solar nebula

Particle accumulation processes are discussed for a variety of physical environments, ranging from the collapse phase of an interstellar cloud to the different parts of the models of the primitive solar nebula constructed by Cameron and Pine. Because of turbulence in the collapsing interstellar gas, it is concluded that interstellar grains accumulate into bodies with radii of a few tens of centimeters before the outer parts of the solar nebula are formed. These bodies can descend quite rapidly through the gas toward midplane of the nebula, and accumulation to planetary size can occur in a few thousand years. Substantial modifications of these processes take place in the outer convection zone of the solar nebula, but again it is concluded that bodies in that zone can grow to planetary size in a few thousand years.From the discussion of the interstellar collapse phase it is concluded that the angular momentum of the primitive solar nebula was predominantly of random turbulent origin, and that it is plausible that the primitive solar nebula should have possessed satellite nebulae in highly elliptical orbits. It is proposed that the comets were formed in these satellite nebulae.A number of other detailed conclusions are drawn from the analysis. It is shown to be plausible that an iron-rich planet should be formed in the inner part of the outer nebular convection zone. Discussions are given of the processes of planetary gas accretion, the formation of satellites, the T Tauri solar wind, and the dissipation of excess condensed material after the nebular gases have been removed by the T Tauri solar wind. It is shown that the present radial distances of the planets (but not Bode's Law) should be predicted reasonably well by a solar nebula model intermediate between the uniform and linear cases of Cameron and Pine.     

Numerical modeling of protocore destabilization during planetary accretion: Methodology and results

We developed and tested an efficient 2D numerical methodology for modeling gravitational redistribution processes in a quasi spherical planetary body based on a simple Cartesian grid. This methodology allows one to implement large viscosity contrasts and to handle properly a free surface and self-gravitation. With this novel method we investigated in a simplified way the evolution of gravitationally unstable global three-layer structures in the interiors of large metal-silicate planetary bodies like those suggested by previous models of cold accretion [Sasaki, S., Nakazawa, K., 1986. J. Geophys. Res. 91, 9231-9238; Karato, S., Murthy, V.R., 1997. Phys. Earth Planet Interios 100, 61-79; Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819]: an innermost solid protocore (either undifferentiated or partly differentiated), an intermediate metal-rich layer (either continuous or disrupted), and an outermost silicate-rich layer. Long-wavelength (degree-one) instability of this three-layer structure may strongly contribute to core formation dynamics by triggering planetary-scale gravitational redistribution processes. We studied possible geometrical modes of the resulting planetary reshaping using scaled 2D numerical experiments for self-gravitating planetary bodies with Mercury-, Mars- and Earth-size. In our simplified model the viscosity of each material remains constant during the experiment and rheological effects of gravitational energy dissipation are not taken into account. However, in contrast to a previously conducted numerical study [Honda, R., Mizutani, H., Yamamoto, T., 1993. J. Geophys. Res. 98, 2075-2089] we explored a freely deformable planetary surface and a broad range of viscosity ratios between the metallic layer and the protocore (0.001-1000) as well as between the silicate layer and the protocore (0.001-1000). An important new prediction from our study is that realistic modes of planetary reshaping characterized by a high viscosity protocore and low viscosity molten silicate and metal [Senshu, H., Kuramoto, K., Matsui, T., 2002. J. Geophys. Res. 107 (E12), 5118. 10.1029/2001JE001819] may result in the transient exposure of the protocore to the planetary surface and a strongly (up to 8% of the planetary diameter) aspherical deviation of the planetary shape during the early stages of core formation. Exposure of the protocore might happen in the early stages of iron core formation. This process may conceivably convert a large amount of potential energy into temperature increase and a transient strongly non-uniform depth of the magma ocean around the protoplanet. Our simplified model also predicts that the time for metallic core formation out of the metal-rich layer depends mainly on the dynamics of the deformation of the solid strong protocore. In nature this dynamics will be strongly dependent on the effective viscosity of the protocore, which should generally have non-Newtonian pressure-, temperature-, and stress-dependent rheology with strong thermomechanical feedbacks from gravitational energy dissipation.     

Simulation of the Magnetic Anomaly Associated with a Complex Crater Using the Example of the Bosumtwi Crater

Solar System Research - The impact crater formation on the surface of the Earth and other planetary bodies is accompanied by the action of shock waves on rocks and their displacement into a new...     

Stochastic coalescence model for terrestrial planetary accretion

A nonequilibrium stochastic coalescence model for terrestrial planetary accretion is developed by using an approximation to the Safronov-Golovin solution for the scalar transport equation with linear kernel. According to this model, formation of comparatively massive objects occurs quite rapidly during the early stages of accretive evolution in a given terrestrial planetesimal population, while during late growth stages, an increasingly substantial fraction of total population mass becomes incorporated into progressively fewer, relatively very large bodies. The model also implies that the (conservative) growth rate of the population's largest member varies directly as its mass, and further suggests that this growth rate may not decline significantly until very nearly final planetary mass is attained.     

天体运动的轨道偏心率研究概述---从恒星系统到行星系统

偏心率是描述天体运动轨道的重要参数之一, 能够为揭示天体的动力学演化提供重要线索, 进而帮助理解天体形成与演化的过程及背后的物理机制. 随着天文观测技术的不断发展, 人们对于天体运动轨道的研究已经走出太阳系, 包含的系统也从大质量端的恒星系统延伸到了低质量端的行星系统. 聚焦天体轨道偏心率研究, 回顾了目前在恒星系统(包括主序恒星、褐矮星以及致密星)和行星系统(包括太阳系外巨行星以及``超级地球''、``亚海王星''等小质量系外行星)方面取得的进展, 总结了不同尺度结构下偏心率研究的一些共同之处和待解决的问题. 并结合当下和未来的相关天文观测设备和项目, 对未来天体轨道偏心率方面的研究工作进行了展望.     

Orbital resonances in the solar nebula: Implications for planetary accretion

The influence of gas drag and gravitational perturbations by a planetary embryo on the orbit of a planetesimal in the solar nebula was examined. Non-Keplerian rotation of the gas causes secular decay of the orbit. If the planetesimal's orbit is exterior to the perturber's, resonant perturbations oppose this drag and can cause it to be trapped in a stable orbit at a commensurability of order j/(j + 1), where j is an integer. Numerical and analytical demonstrations show that resonant trapping occurs for wide ranges of perturbing mass, planetesimal size, and j. Induced eccentricities are large, causing overlap of orbits for bodies in different resonances with j > 2. Collisions between planetesimals in different resonances, or between resonant and nonresonant bodies, result in their disruption. Fragments smaller than a critical size can pass through resonances under the influence of drag and be accreted by the embryo. This effect speeds accretion and tends to prevent dynamical isolation of planetary embryos, making gas-rich scenarios for planetary formation more plausible.     

Clues on the importance of comets in the origin and evolution of the atmospheres of Titan and Earth

Earth and Titan are two planetary bodies formed far from each other. Nevertheless the chemical composition of their atmospheres exhibits common indications of being produced by the accretion, plus ulterior in-situ processing of cometary materials. This is remarkable because while the Earth formed in the inner part of the disk, presumably from the accretion of rocky planetesimals depleted in oxygen and exhibiting a chemical similitude with enstatite chondrites, Titan formed within Saturn's sub-nebula from oxygen- and volatile-rich bodies, called cometesimals. From a cosmochemical and astrobiological perspective, the study of the H, C, N, and O isotopes on Earth and Titan could be the key to decipher the processes occurred in the early stages of formation of both planetary bodies. The main goal of this paper is to quantify the presumable ways of chemical evolution of both planetary bodies, in particular the abundance of CO and N₂ in their early atmospheres. In order to do that the primeval atmospheres and evolution of Titan and Earth have been analyzed from a thermodynamic point of view. The most relevant chemical reactions involving these species and presumably important at their early stages are discussed. Then, we have interpreted the results of this study in light of the results obtained by the Cassini–Huygens mission on these species and their isotopes. Given that H, C, N, and O were preferentially depleted from inner disk materials that formed our planet, the observed similitude of their isotopic fractionation, and subsequent close evolution of Earth's and Titan's atmospheres points towards a cometary origin of Earth atmosphere. Consequently, our scenario also supports the key role of late veneers (comets and water-rich carbonaceous asteroids) enriching the volatile content of the Earth at the time of the late heavy bombardment of terrestrial planets.     

Origin of the atmospheres of the terrestrial planets

The monotonic decrease in the atmospheric abundance of ³⁶Ar per gram of planet in the sequence, Venus, Earth, and Mars has been assumed to reflect some conditions in the primitive solar nebula at the time of formation of the planetary atmospheres, having to do either with the composition of the nebula itself or the composition of the trapped gases in small solid bodies in the nebula. Behind such hypotheses lies the assumption that planetary atmospheres steadily gain components. However, not only can gases enter atmospheres; they may also be lost from atmospheres both by adsorption into the planetary interior and by loss into space as a result of collisions with minor and major planetesimals. In this paper a necessarily qualitative discussion is given of the problem of collisions with minor planetesimals, a process called atmospheric cratering or atmospheric erosion, and a discussion is given of atmospheric loss accompanying collision of a planet with a major planetesimal, such as may have produced the Earth's Moon.     

Accretion of the asteroids: Implications for their thermal evolution

Thermal models of asteroids generally assume that they accreted either instantaneously or over an extended interval with a prescribed growth rate. It is conventionally assumed that the onset of accretion of chondrite parent bodies was delayed until a substantial fraction of the initial ²⁶Al had decayed. However, this interval is not consistent with the early melting, and differentiation of parent bodies of iron meteorites. Formation time scales are tested by dynamical simulations of accretion from small primary planetesimals. Gravitational accretion yields rapid runaway growth of large planetary embryos until most smaller bodies are depleted. In a given simulation, all asteroid‐sized bodies have comparable growth times, regardless of size. For plausible parameters, growth times are shorter than the lifetime of ²⁶Al, consistent with thermal models that assume instantaneous accretion. Rapid growth after planetesimal formation is consistent with differentiation of parent bodies of iron meteorites, but not with the assumed delay in formation of chondritic bodies. After the initial growth stage, there is an interval of slower evolution until the belt is stirred and the embryos are dynamically removed. During this interval, a fraction of asteroid‐sized bodies experience large accretional impacts, allowing bodies of the same final size to have very different histories of radius versus time. Accretion from small primary planetesimals leaves some fraction of material in bodies small enough to preserve CAIs while avoiding heating by ²⁶Al. Unheated material can be a significant fraction of the mass that remains after large embryos are removed from the Main Belt.     

Possible patterns in the distribution of planetary formation regions

Eris, an object larger than Pluto, is known to reside in the transneptunian region further away than Pluto. One can wonder whether its semimajor orbital axis fits in a generalized Titius–Bode law, in the same way as Pluto does. We performed a new least-squares fit to a generalized Titius–Bode law including Eris and found that not only does Eris fit in the trend, but also the correlation coefficient improves. In addition, there is a remarkable symmetry of the location of the planetary formation regions with respect to Jupiter when the natural logarithm of the heliocentric distance is used as the metric. The issue of whether the observed patterns have some physical meaning or are due to mere chance is addressed using a Monte Carlo approach identical to that of Lynch. Although the probability of chance occurrence is highly dependent on the way in which the random configurations of synthetic planetary systems are selected, we find that in all reasonable scenarios of random planetary systems the probability of chance occurrence of the observed patterns is small (below 1 per cent in most cases). If the trend were used as a prediction tool, one might expect another planet or dwarf planet or a swarm of bodies with semimajor orbital axis of 120 ± 20 au. Simple calculations show that the protoplanetary nebula most likely had enough mass to allow the accretion of at least a dwarf planet at that distance. We also found that if the surface density of the nebula decayed with heliocentric distance ( r ) as a power of −2, the regular spacing in ln  r in the Solar system could be a natural consequence of the existence of a threshold mass for planetary formation.     

High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

The origin of the solar system: Implications for transneptunian planets and the nature of the long-period comets

If the solar system origin is considered within the framework of the author's hypothesis on the binary stars formation as a result of rotational-exchange break-up of the rotating protostar, then difficulties involved in the usual nebular hypotheses are automatically removed (unclear aspects of the possibility of formation of the gas disc proper, the problems of the angular momentum including slow rotation of the Sun and coplanarity of the planetary orbits, of differences in planetary masses and composition, the need, for the disc remnants to be swept out, the long time of planetary formation as compared with the possible lifetime of a turbulized disc etc.).The major stages of division and evolution of the Jupiter-Sun system are described. Similarities between the massive rotating proto-Jupiter (PJ) and the classical protoplanetary discs are pointed out. The process of planetoid condensation inside PJ is discussed. The most probable site of the condensation is the region of the first Lagrangian point. The planetoids condensed were lost by PJ as a result of its fast mass decrease. A gas dynamic consideration of the motion of planetoids in PJ yields 1000–3000 yr as a time scale for the PJ's mass loss. The number of the moonlike bodies lost (the remaining Galilean satellites fixing their lower mass limit) could reach 10⁴.Evolution of such interacting bodies results in the formation beyond Neptune of a cloud (up to 10³) of moonlike (and more massive) planets.The excess concentration of the long-period comets aphelia in this area implies their genetic relation to the planets. A concept of a joint planeto-cometary cloud is introduced. A concrete hydrodynamic mechanism of ice ejection from planets into space, viz. the formation of cumulative (Monroe) jets, is pointed out.A program of further investigations is outlined and recommendations given for an experimental check on the implications of the new cosmogonic concepts.     

The influence of shapes of the bodies on the dynamics of planetary and star systems

The dynamics of planetary and star systems including perturbing forces due to the flattering and quadrupole distortion of the bodies is studied. The analytical model describing the perturbations which influence on the orbital motion of extrasolar planetary systems is presented. The calculations of the secular evolution of the mean orbital elements have shown that the effects related to the shape of the body are more important than the ones due to the quadrupole distortion.     

The first 800 million years: Environmental models for early earth

Space exploration has changed our views not only on the properties of celestial bodies and the interplanetary medium but also our perspective on the formation of such bodies, including the Earth. However, on few points do the new hetegonic insights approach certainty. More appropriately, they can be said to raise physically meaningful questions and to decrease the degrees of freedom allowable in any given scenario, which necessarily extends into the early history of the planetary crust, ocean and atmosphere. An attempt is made here to outline the basic schemes envisaged in modern studies of the origin of the Earth and the ramification of consequences for its early evolution that currently appear allowable and supported by observation.