Envelope instability in giant planet formation

We compute the growth of isolated gaseous giant planets for several values of the density of the protoplanetary disk, several distances from the central star and two values for the (fixed) radii of accreted planetesimals. Calculations were performed in the frame of the core instability mechanism and the solids accretion rate adopted is that corresponding to the oligarchic growth regime. We find that for massive disks and/or for protoplanets far from the star and/or for large planetesimals, the planetary growth occurs smoothly. However, notably, there are some cases for which we find an envelope instability in which the planet exchanges gas with the surrounding protoplanetary nebula. The timescale of this instability shows that it is associated with the process of planetesimals accretion. The presence of this instability makes it more difficult the formation of gaseous giant planets.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

Making other earths: dynamical simulations of terrestrial planet formation and water delivery

We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×10²⁴ g), and vary in mass between 0.23M_⊕ and 3.85M_⊕. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M_⊕ protoplanets which can serve as the accretion seeds for 2+M_⊕ planets with large water contents.     

Accretion rates of planetesimals by protoplanets embedded in nebular gas

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals’ eccentricity, inclination, planet’s radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals’ kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with 1-10 Earth masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.     

Oligarchic growth of giant planets

Runaway growth ends when the largest protoplanets dominate the dynamics of the planetesimal disk; the subsequent self-limiting accretion mode is referred to as “oligarchic growth.” Here, we begin by expanding on the existing analytic model of the oligarchic growth regime. From this, we derive global estimates of the planet formation rate throughout a protoplanetary disk. We find that a relatively high-mass protoplanetary disk (∼10 × minimum-mass) is required to produce giant planet core-sized bodies (∼10 M_⊕) within the lifetime of the nebular gas (?10 million years). However, an implausibly massive disk is needed to produce even an Earth mass at the orbit of Uranus by 10 Myrs. Subsequent accretion without the dissipational effect of gas is even slower and less efficient. In the limit of noninteracting planetesimals, a reasonable-mass disk is unable to produce bodies the size of the Solar System’s two outer giant planets at their current locations on any timescale; if collisional damping of planetesimal random velocities is sufficiently effective, though, it may be possible for a Uranus/Neptune to form in situ in less than the age of the Solar System. We perform numerical simulations of oligarchic growth with gas and find that protoplanet growth rates agree reasonably well with the analytic model as long as protoplanet masses are well below their estimated final masses. However, accretion stalls earlier than predicted, so that the largest final protoplanet masses are smaller than those given by the model. Thus the oligarchic growth model, in the form developed here, appears to provide an upper limit for the efficiency of giant planet formation.     

A semi-analytic model for oligarchic growth

A new semi-analytic model for the oligarchic growth phase of planetary accretion is developed. The model explicitly calculates damping and excitation of planetesimal eccentricities e and inclinations i due to gas drag and perturbations from embryos. The effects of planetesimal fragmentation, enhanced embryo capture cross sections due to atmospheres, inward planetesimal drift, and embryo-embryo collisions are also incorporated. In the early stages of oligarchic growth, embryos grow rapidly as e and i fall below their equilibrium values. The formation of planetesimal collision fragments also speeds up embryo growth as fragments have low-e, low-i orbits, thereby optimizing gravitational focussing. At later times, the presence of thick atmospheres captured from the nebula aids embryo growth by increasing their capture cross sections. Planetesimal drift due to gas drag can lead to substantial inward transport of solid material. However, inward drift is greatly reduced when embryo atmospheres are present, as the drift timescale is no longer short compared to the accretion timescale. Embryo-embryo collisions increase embryo growth rates by 50% compared to the case where growth is solely due to accretion of planetesimals. Formation of 0.1-Earth-mass protoplanets at 1 AU and 10-Earth-mass cores at 5 AU requires roughly 0.1 and 1 million years respectively, in a nebula where the local solid surface density is 7 g cm⁻² at each of these locations.     

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

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.     

The effects of ablation on the cross section of planetary envelopes at capturing planetesimals

We explore the cross section of giant planet envelopes at capturing planetesimals of different sizes. For this purpose we employ two sets of realistic planetary envelope models (computed assuming for the protoplanetary nebula masses of 10 and 5 times the mass of the minimum mass solar nebula), account for drag and ablation effects and study the trajectories along which planetesimals move. The core accretion of these models has been computed in the oligarchic growth regime [Fortier, A., Benvenuto, O.G., Brunini, A., 2007. Astron. Astrophys. 473, 311-322], which has also been considered for the velocities of the incoming planetesimals. This regime predicts velocities larger that those used in previous studies of this problem. As the rate of ablation is dependent on the third power of velocity, ablation is more important in the oligarchic growth regime. We compute energy and mass deposition, fractional ablated masses and the total cross section of planets for a wide range of values of the critical parameter of ablation. In computing the total cross section of the planet we have included the contributions due to mass deposited by planetesimals moving along unbound orbits. Our results indicate that, for the case of small planetary cores and low velocities for the incoming planetesimals, ablation has a negligible impact on the capture cross section in agreement with the results presented in Inaba and Ikoma [Inaba, S., Ikoma, M., 2003. Astron. Astrophys. 410, 711-723]. However for the case of larger cores and high velocities of the incoming planetesimals as predicted by the oligarchic growth regime, we find that ablation is important in determining the planetary cross section, being several times larger than the value corresponding ignoring ablation. This is so regardless of the size of the incoming planetesimals.     

Metallicity of the massive protoplanets around HR 8799 If formed by gravitational instability

The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their ‘pre-collapse’ stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be and 5.67×10³ years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions.     

Formation of early water oceans on rocky planets

Terrestrial planets, with silicate mantles and metallic cores, are likely to obtain water and carbon compounds during accretion. Here I examine the conditions that allow early formation of a surface water ocean (simultaneous with cooling to clement surface conditions), and the timeline of degassing the planetary interior into the atmosphere. The greatest fraction of a planet’s initial volatile budget is degassed into the atmosphere during the end of magma ocean solidification, leaving only a small fraction of the original volatiles to be released into the atmosphere through later volcanism. Rocky planets that accrete with water in their bulk mantle have two mechanisms for producing an early water ocean: First, if they accrete with at least 1 to 3 mass% of water in their bulk composition, liquid water may be extruded onto the planetary surface at the end of magma ocean solidification. Second, at initial water contents as low as 0.01 mass% or lower, during solidification a massive supercritical fluid and steam atmosphere is produced that collapses into a water ocean upon cooling. The low water contents required for this process indicate that rocky super-Earth exoplanets may be expected to commonly produce water oceans within tens to hundreds of millions of years of their last major accretionary impact, through collapse of their atmosphere.     

Planetesimals to protoplanets – II. Effect of debris on terrestrial planet formation

In this paper, we extend our numerical method for simulating terrestrial planet formation to include dynamical friction from the unresolved debris component. In the previous work, we implemented a rubble pile planetesimal collision model into direct N -body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10 per cent or more of the total disc mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead, they grow concurrently, similar to oligarchic growth. The dynamical friction from the unresolved debris damps the eccentricities of the planetesimals, reducing the mean impact speeds and causing all collisions to result in merging with no mass loss. As a result, there is no debris production. The mass in debris slowly decreases with time. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions (for which the initial mass in unresolved debris is at least equal to the mass in resolved planetesimals).     

Formation of the first oxidized iron in the solar system

For fayalite formation times of several thousand years, and systems enriched in water by a factor of ten relative to solar composition, 1 μm radius olivine grains could reach 2 mole% fayalite and 0.1 μm grains 5 mole% by nebular condensation, well short of the values appropriate for precursors of most chondrules and the values found in the matrices of unequilibrated ordinary chondrites. Even 10 μm olivine crystals could reach 30 mole% fayalite above 1100 K in solar gas if condensation of metallic nickel‐iron were delayed sufficiently by supersaturation. Consideration of the surface tensions of several phases with equilibrium condensation temperatures above that of metallic iron shows that, even if they were supersaturated, they would still nucleate homogeneously above the equilibrium condensation temperature of metallic iron. This phenomenon would have provided nuclei for heterogeneous nucleation of metallic nickel‐iron, thus preventing the latter from supersaturating significantly and preventing olivine from becoming fayalitic. Unless a way is found to make nebular regions far more oxidizing than in existing models, it is unlikely that chondrule precursors or the matrix olivine grains of unequilibrated ordinary chondrites obtained their fayalite contents by condensation processes. Perhaps stabilization of FeO occurred after condensation of water ice and accretion of icy planetesimals, during heating of the planetesimals and/or in hot, dense, water‐rich vapor plumes generated by impacts on them. This would imply that FeO is a relatively young feature of nebular materials.     

Core instability models of giant planet accretion – II. Forming planetary systems

We develop a simple model for computing planetary formation based on the core instability model for the gas accretion and the oligarchic growth regime for the accretion of the solid core. In this model several planets can form simultaneously in the disc, a fact that has important implications especially for the changes in the dynamic of the planetesimals and the growth of the cores since we consider the collision between them as a source of potential growth. The type I and type II migration of the embryos and the migration of the planetesimals due to the interaction with the disc of gas are also taken into account. With this model we consider different initial conditions to generate a variety of planetary systems and analyse them statistically. We explore the effects of using different type I migration rates on the final number of planets formed per planetary system such as on the distribution of masses and semimajor axis of extrasolar planets, where we also analyse the implications of considering different gas accretion rates. A particularly interesting result is the generation of a larger population of habitable planets when the gas accretion rate and type I migration are slower.     

Core formation in giant gaseous protoplanets

Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 MSaturn to 10 MJupiter. Giant protoplanets with masses of 5 MJupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ∼0.25 M_⊕ and ∼0.5 M_⊕, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M_⊕. Saturn's core mass is found to be ∼8 M_⊕.     

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

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.     

On oligarchic growth of planets in protoplanetary disks

In this paper we present a new semianalytical model of oligarchic growth of planets considering a distribution of planetesimal sizes, fragmentation of planetesimals in mutual collisions, sublimation of ices through the snow line, random velocities out of equilibrium and merging of planetary embryos. We show that the presence of several planetary embryos growing simultaneously at different locations in the protoplanetary disk affects the whole accretion history, specially for the innermost planets. The results presented here clearly indicate the relevance of considering a distribution of planetesimal sizes. Fragmentation occurring during planetesimal-planetesimal collisions represent only a marginal effect in shaping the surface density of solid material in the protoplanetary disc.     

A hybrid scenario for gas giant planet formation in rings

The core-accretion mechanism for gas giant formation may be too slow to create all observed gas giant planets during reasonable gas disk lifetimes, but it has yet to be firmly established that the disk instability model can produce permanent bound gaseous protoplanets under realistic conditions. Based on our recent simulations of gravitational instabilities in disks around young stars, we suggest that, even if instabilities due to disk self-gravity do not produce gaseous protoplanets directly, they may create persistent dense rings that are conducive to accelerated growth of gas giants through core accretion. The rings occur at and near the boundary between stable and unstable regions of the disk and appear to be produced by resonances with discrete spiral modes on the unstable side.     

Co-accretion of the Earth-Moon system after the giant impact: reduction of the total angular momentum by lunar impact ejecta

Though the Moon is considered to have been formed by the so-called giant impact, the mass of the Earth immediately after the impact is still controversial. If the Moon was formed during the Earth's accretion, a subsequent accretion of residual heliocentric planetesimals onto the protoearth and the protomoon must have occurred. In this co-accretion stage, a significant amount of lunar-impact-ejecta would be ejected to circumterrestrial orbits, since the mean impact velocity of the planetesimals with the protomoon is much larger than the escape velocity of the protomoon. Orbital calculations of test particles ejected from the protomoon, whose semimajor axis is smaller than that of the present Moon, reveal that most of the particles escaping from the protomoon also escape from the Hill sphere of the protoearth and reduce the planetocentric angular momentum of the primordial Earth-Moon system. Using the results of the ejecta simulations, we investigate the evolution of the mass ratio and the total angular momentum (Earth's spin angular momentum + Moon's orbital angular momentum) of the Earth-Moon system during the co-accretion. We find that the mass of the protomoon is almost constant or rather decreases and the total angular momentum decreases significantly, if the random velocity of planetesimals is as large as the escape velocity of the protoearth. On the other hand, if the random velocity is the half of the escape velocity of the protoearth, the mass ratio is kept to be almost as large as the present value and the decrease of the total angular momentum is not so significant. Comparing with the results of giant impact simulations, we find that the mass of the protoearth immediately after the Moon-forming impact was 0.7-0.8 times the present value if the impactor-to-target mass ratio was 3:7, whereas the giant impact occurred almost in the end of the Earth's accretion if the impactor-to-target mass ratio was 1:9.     

The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

Implications of an anti-glitch in AXP/SGR

The recently observed anti-glitch of AXP 1E 2259+586 may be explained as the consequence of sudden accretion of retrograde matter. AXP/SGR are explained as single neutron stars accompanied by fallback matter from their natal supernovæ, including rocky or metallic planetesimals. The upper bounds on a glitch or anti-glitch in the giant outburst of SGR 1806-20 pose a problem for any model.